title: Distributed Stochastic Optimization With Unbounded Subgradients Over Randomly Time-Varying Networks

abstract: Motivated by distributed statistical learning over uncertain communication networks, we study distributed stochastic optimization by networked nodes to cooperatively minimize a sum of convex cost functions. The network is modeled by a sequence of time-varying random digraphs with each node representing a local optimizer and each edge representing a communication link. We consider the distributed subgradient optimization algorithm with noisy measurements of local cost functions' subgradients, additive and multiplicative noises among information exchanging between each pair of nodes. By stochastic Lyapunov method, convex analysis, algebraic graph theory and martingale convergence theory, we prove that if the local subgradient functions grow linearly and the sequence of digraphs is conditionally balanced and uniformly conditionally jointly connected, then proper algorithm step sizes can be designed so that all nodes' states converge to the global optimal solution almost surely.


title: Dp-minimal profinite groups and valuations on the integers

abstract: We study dp-minimal infinite profinite groups that are equipped with a uniformly definable fundamental system of open subgroups. We show that these groups have an open subgroup $A$ such that either $A$ is a direct product of countably many copies of $\mathbb{F}_p$ for some prime $p$, or $A$ is of the form $A \cong \prod_p \mathbb{Z}_p^{\alpha_p} \times A_p$ where $\alpha_p < \omega$ and $A_p$ is a finite abelian $p$-group for each prime $p$. Moreover, we show that if $A$ is of this form, then there is a fundamental system of open subgroups such that the expansion of $A$ by this family of subgroups is dp-minimal. Our main ingredient is a quantifier elimination result for a class of valued abelian groups. We also apply it to $(\mathbb{Z},+)$ and we show that if we expand $(\mathbb{Z},+)$ by any chain of subgroups $(B_i)_{i<\omega}$, we obtain a dp-minimal structure. This structure is distal if and only if the size of the quotients $B_i/B_{i+1}$ is bounded.


title: Existence of EFX for Two Additive Valuations

abstract: Fair division of indivisible items is a well-studied topic in Economics and Computer Science.The objective is to allocate items to agents in a fair manner, where each agent has a valuation for each subset of items. Envy-freeness is one of the most widely studied notions of fairness. Since complete envy-free allocations do not always exist when items are indivisible, several relaxations have been considered. Among them, possibly the most compelling one is envy-freeness up to any item (EFX), where no agent envies another agent after the removal of any single item from the other agent's bundle. However, despite significant efforts by many researchers for several years, it is known that a complete EFX allocation always exists only in limited cases. In this paper, we show that a complete EFX allocation always exists when each agent is of one of two given types, where agents of the same type have identical additive valuations. This is the first such existence result for non-identical valuations when there are any number of agents and items and no limit on the number of distinct values an agent can have for individual items. We give a constructive proof, in which we iteratively obtain a Pareto dominating (partial) EFX allocation from an existing partial EFX allocation.


title: Emergent particles and gauge fields in quantum matter

abstract: I give a pedagogical introduction to some of the many particles and gauge fields that can emerge in correlated matter. The standard model of materials is built on Landau's foundational principles: adiabatic continuity and spontaneous symmetry breaking. These ideas lead to quasiparticles that inherit their quantum numbers from fundamental particles, Nambu-Goldstone bosons, the Anderson-Higgs mechanism, and topological defects in order parameters. I then describe the modern discovery of physics beyond the standard model. Here, quantum correlations (entanglement) and topology play key roles in defining the properties of matter. This can lead to fractionalised quasiparticles that carry only a fraction of the quantum numbers that define fundamental particles. These particles can have exotic properties: for example Majorana fermions are their own antiparticles, anyons have exchange statistics that are neither bosonic nor fermionic, and magnetic monopoles do not occur in the vacuum. Gauge fields emerge naturally in the description of highly correlated matter and can lead to gauge bosons. Relationships to the standard model of particle physics are discussed.

categories: cond-mat.str-el cond-mat.quant-gas cond-mat.supr-con hep-ph quant-ph


title: A Metabolite Specific 3D Stack-of-Spiral bSSFP Sequence for Improved Lactate Imaging in Hyperpolarized [1-$^{13}$C]Pyruvate Studies on a 3T Clinical Scanner

abstract: Purpose: The balanced steady-state free precession sequence has been previously explored to improve the efficient use of non-recoverable hyperpolarized $^{13}$C magnetization, but suffers from poor spectral selectivity and long acquisition time. The purpose of this study was to develop a novel metabolite-specific 3D bSSFP ("MS-3DSSFP") sequence with stack-of-spiral readouts for improved lactate imaging in hyperpolarized [1-$^{13}$C]pyruvate studies on a clinical 3T scanner. Methods: Simulations were performed to evaluate the spectral response of the MS-3DSSFP sequence. Thermal $^{13}$C phantom experiments were performed to validate the MS-3DSSFP sequence. In vivo hyperpolarized [1-$^{13}$C]pyruvate studies were performed to compare the MS-3DSSFP sequence with metabolite specific gradient echo ("MS-GRE") sequences for lactate imaging. Results: Simulations, phantom and in vivo studies demonstrate that the MS-3DSSFP sequence achieved spectrally selective excitation on lactate while minimally perturbing other metabolites. Compared with MS-GRE sequences, the MS-3DSSFP sequence showed approximately a 2.5-fold SNR improvement for lactate imaging in rat kidneys, prostate tumors in a mouse model and human kidneys. Conclusions: Improved lactate imaging using the MS-3DSSFP sequence in hyperpolarized [1-$^{13}$C]pyruvate studies was demonstrated in animals and humans. The MS-3DSSFP sequence could be applied for other clinical applications such as in the brain or adapted for imaging other metabolites such as pyruvate and bicarbonate.


title: Preferential evaporation in atom probe tomography: an analytical approach

abstract: Atom probe tomography (APT) analysis conditions play a major role in the composition measurement accuracy. Preferential evaporation, which significantly biases apparent composition, more than other well-known phenomena in APT, is strongly connected to those analysis conditions. One way to optimize them, in order to have the most accurate measurement, is therefore to be able to predict and then to estimate their influence on the apparent composition. An analytical model is proposed to quantify the preferential evaporation. This model is applied to three different alloys: NiCu, FeCrNi and FeCu. The model explains not only the analysis temperature dependence, as in already existing model, but also the dependence to the pulse fraction and the pulse frequency. Moreover, the model can also provide energetic constant directly linked to energy barrier, required to field evaporate atom from the sample surface. 2


title: Garside groups and geometry

abstract: This article in memory of Patrick Dehornoy (1952-2019) is an invitation to Garside theory for mainstream geometric group theorists interested in mapping class groups, curve complexes, and the geometry of Artin-Tits groups.


title: Probing Fundamental Symmetries of Deformed Nuclei in Symmetric Top Molecules

abstract: Precision measurements of Schiff moments in heavy, deformed nuclei are sensitive probes of beyond Standard Model $T,P$-violation in the hadronic sector. While the most sensitive limits on Schiff moments to date are set with diamagnetic atoms, polar polyatomic molecules can offer higher sensitivities with unique experimental advantages. In particular, symmetric top molecular ions possess $K$-doublets of opposite parity with especially small splittings, leading to full polarization at low fields, internal co-magnetometer states useful for rejection of systematic effects, and the ability to perform sensitive searches for $T,P$-violation using a small number of trapped ions containing heavy exotic nuclei. We consider the symmetric top cation $^{225}$RaOCH$_3^+$ as a prototypical and candidate platform for performing sensitive nuclear Schiff measurements and characterize in detail its internal structure using relativistic ab initio methods. The combination of enhancements from a deformed nucleus, large polarizability, and unique molecular structure make this molecule a promising platform to search for fundamental symmetry violation even with a single trapped ion.


title: Assessing the Quality-of-Experience of Adaptive Bitrate Video Streaming

abstract: The diversity of video delivery pipeline poses a grand challenge to the evaluation of adaptive bitrate (ABR) streaming algorithms and objective quality-of-experience (QoE) models. Here we introduce so-far the largest subject-rated database of its kind, namely WaterlooSQoE-IV, consisting of 1350 adaptive streaming videos created from diverse source contents, video encoders, network traces, ABR algorithms, and viewing devices. We collect human opinions for each video with a series of carefully designed subjective experiments. Subsequent data analysis and testing/comparison of ABR algorithms and QoE models using the database lead to a series of novel observations and interesting findings, in terms of the effectiveness of subjective experiment methodologies, the interactions between user experience and source content, viewing device and encoder type, the heterogeneities in the bias and preference of user experiences, the behaviors of ABR algorithms, and the performance of objective QoE models. Most importantly, our results suggest that a better objective QoE model, or a better understanding of human perceptual experience and behaviour, is the most dominating factor in improving the performance of ABR algorithms, as opposed to advanced optimization frameworks, machine learning strategies or bandwidth predictors, where a majority of ABR research has been focused on in the past decade. On the other hand, our performance evaluation of 11 QoE models shows only a moderate correlation between state-of-the-art QoE models and subjective ratings, implying rooms for improvement in both QoE modeling and ABR algorithms. The database is made publicly available at: \url{https://ece.uwaterloo.ca/~zduanmu/waterloosqoe4/}.


title: Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

abstract: In this paper, we develop an adaptive finite element method for the nonlinear steady-state Poisson-Nernst-Planck equations, where the spatial adaptivity for geometrical singularities and boundary layer effects are mainly considered. As a key contribution, the steady-state Poisson-Nernst-Planck equations are studied systematically and rigorous analysis for a residual-based a posteriori error estimate of the nonlinear system is presented. With the help of Schauder fixed point theorem, we show the solution existence and uniqueness of the linearized system derived by taking $G-$derivatives of the nonlinear system, followed by the proof of the relationship between the error of solution and the a posteriori error estimator $\eta$. Numerical experiments are given to validate the efficiency of the a posteriori error estimator and demonstrate the expected rate of convergence. In the further tests, adaptive mesh refinements for geometrical singularities and boundary layer effects are successfully observed.


title: Varying Annotations in the Steps of the Visual Analysis

abstract: Annotations in Visual Analytics (VA) have become a common means to support the analysis by integrating additional information into the VA system. That additional information often depends on the current process step in the visual analysis. For example, the data preprocessing step has data structuring operations while the data exploration step focuses on user interaction and input. Describing suitable annotations to meet the goals of the different steps is challenging. To tackle this issue, we identify individual annotations for each step and outline their gathering and design properties for the visual analysis of heterogeneous clinical data. We integrate our annotation design into a visual analysis tool to show its applicability to data from the ophthalmic domain. In interviews and application sessions with experts we asses its usefulness for the analysis of patients with different medications.


title: Not one but many Tradeoffs: Privacy Vs. Utility in Differentially Private Machine Learning

abstract: Data holders are increasingly seeking to protect their user's privacy, whilst still maximizing their ability to produce machine models with high quality predictions. In this work, we empirically evaluate various implementations of differential privacy (DP), and measure their ability to fend off real-world privacy attacks, in addition to measuring their core goal of providing accurate classifications. We establish an evaluation framework to ensure each of these implementations are fairly evaluated. Our selection of DP implementations add DP noise at different positions within the framework, either at the point of data collection/release, during updates while training of the model, or after training by perturbing learned model parameters. We evaluate each implementation across a range of privacy budgets, and datasets, each implementation providing the same mathematical privacy guarantees. By measuring the models' resistance to real world attacks of membership and attribute inference, and their classification accuracy. we determine which implementations provide the most desirable tradeoff between privacy and utility. We found that the number of classes of a given dataset is unlikely to influence where the privacy and utility tradeoff occurs. Additionally, in the scenario that high privacy constraints are required, perturbing input training data does not trade off as much utility, as compared to noise added later in the ML process.


title: BGC: Multi-Agent Group Belief with Graph Clustering

abstract: Recent advances have witnessed that value decomposed-based multi-agent reinforcement learning methods make an efficient performance in coordination tasks. Most current methods assume that agents can make communication to assist decisions, which is impractical in some situations. In this paper, we propose a semi-communication method to enable agents can exchange information without communication. Specifically, we introduce a group concept to help agents learning a belief which is a type of consensus. With this consensus, adjacent agents tend to accomplish similar sub-tasks to achieve cooperation. We design a novel agent structure named Belief in Graph Clustering(BGC), composed of an agent characteristic module, a belief module, and a fusion module. To represent each agent characteristic, we use an MLP-based characteristic module to generate agent unique features. Inspired by the neighborhood cognitive consistency, we propose a group-based module to divide adjacent agents into a small group and minimize in-group agents' beliefs to accomplish similar sub-tasks. Finally, we use a hyper-network to merge these features and produce agent actions. To overcome the agent consistent problem brought by GAT, a split loss is introduced to distinguish different agents. Results reveal that the proposed method achieves a significant improvement in the SMAC benchmark. Because of the group concept, our approach maintains excellent performance with an increase in the number of agents.


title: Canted Antiferromagnetism in the Quasi-1D Iron Chalcogenide BaFe$_{2}$Se$_{4}$

abstract: We report the synthesis and physical properties studies of quais-1D iron chalcogenide $\rm BaFe_2Se_4$ which shares the $\rm FeSe_4$ tetrahedra building motif commonly seen in the iron chalcogenide superconductors. A high-quality polycrystalline sample was achieved by solid-state reaction method and characterized by X-ray diffraction, electrical resistivity, magnetic susceptibility and neutron diffraction measurements. $\rm BaFe_2Se_4$ is a narrow gap semiconductor that magnetically orders at $\sim$ 310 K. Both neutron powder diffraction results and isothermal M-H loops suggest a canted antiferromagnetic structure where Fe sublattice are antiferromagnetically ordered along the c-axis quasi-1D chain direction, resulting in a net ferromagnetic moment in the perpendicular direction along the a-axis with tilted angle of 18.7$^\circ$ towards the b-axis.


title: An Experimental Study of Deep Neural Network Models for Vietnamese Multiple-Choice Reading Comprehension

abstract: Machine reading comprehension (MRC) is a challenging task in natural language processing that makes computers understanding natural language texts and answer questions based on those texts. There are many techniques for solving this problems, and word representation is a very important technique that impact most to the accuracy of machine reading comprehension problem in the popular languages like English and Chinese. However, few studies on MRC have been conducted in low-resource languages such as Vietnamese. In this paper, we conduct several experiments on neural network-based model to understand the impact of word representation to the Vietnamese multiple-choice machine reading comprehension. Our experiments include using the Co-match model on six different Vietnamese word embeddings and the BERT model for multiple-choice reading comprehension. On the ViMMRC corpus, the accuracy of BERT model is 61.28% on test set.


title: Faster Heuristics for Graph Burning

abstract: Graph burning is a process of information spreading through the network by an agent in discrete steps. The problem is to find an optimal sequence of nodes which have to be given information so that the network is covered in least number of steps. Graph burning problem is NP-Hard for which two approximation algorithms and a few heuristics have been proposed in the literature. In this work, we propose three heuristics, namely, Backbone Based Greedy Heuristic (BBGH), Improved Cutting Corners Heuristic (ICCH) and Component Based Recursive Heuristic (CBRH). These are mainly based on Eigenvector centrality measure. BBGH finds a backbone of the network and picks vertex to be burned greedily from the vertices of the backbone. ICCH is a shortest path based heuristic and picks vertex to burn greedily from best central nodes. The burning number problem on disconnected graphs is harder than on the connected graphs. For example, burning number problem is easy on a path where as it is NP-Hard on disjoint paths. In practice, large networks are generally disconnected and moreover even if the input graph is connected, during the burning process the graph among the unburned vertices may be disconnected. For disconnected graphs, ordering of the components is crucial. Our CBRH works well on disconnected graphs as it prioritizes the components. All the heuristics have been implemented and tested on several bench-mark networks including large networks of size more than $50$K nodes. The experimentation also includes comparison to the approximation algorithms. The advantages of our algorithms are that they are much simpler to implement and also several orders faster than the heuristics proposed in the literature.


title: Expressing Diverse Human Driving Behavior with Probabilistic Rewards and Online Inference

abstract: In human-robot interaction (HRI) systems, such as autonomous vehicles, understanding and representing human behavior are important. Human behavior is naturally rich and diverse. Cost/reward learning, as an efficient way to learn and represent human behavior, has been successfully applied in many domains. Most of traditional inverse reinforcement learning (IRL) algorithms, however, cannot adequately capture the diversity of human behavior since they assume that all behavior in a given dataset is generated by a single cost function.In this paper, we propose a probabilistic IRL framework that directly learns a distribution of cost functions in continuous domain. Evaluations on both synthetic data and real human driving data are conducted. Both the quantitative and subjective results show that our proposed framework can better express diverse human driving behaviors, as well as extracting different driving styles that match what human participants interpret in our user study.


title: Polarization amplification by spin-doping in nanomagnetic/graphene hybrid systems

abstract: The generation of non-equilibrium electron spin polarization, spin transport, and spin detection are fundamental in many quantum devices. We demonstrate that a lattice of magnetic nanodots enhances the electron spin polarization in monolayer graphene via carrier exchange. We probed the spin polarization through a resistively-detected variant of electron spin resonance (ESR) and observed resonance amplification mediated by the presence of the nanodots. Each nanodot locally injects a surplus of spin-polarized carriers into the graphene, and the ensemble of all "spin hot spots" generates a non-equilibrium electron spin polarization in the graphene layer at macroscopic lengths. This occurs whenever the interdot distance is comparable or smaller than the spin diffusion length.


title: An RNA-Based Theory of Natural Universal Computation

abstract: Life is confronted with computation problems in a variety of domains including animal behavior, single-cell behavior, and embryonic development. Yet we currently do not know of a naturally existing biological system that is capable of universal computation, i.e., Turing-equivalent in scope. Generic finite-dimensional dynamical systems (which encompass most models of neural networks, intracellular signaling cascades, and gene regulatory networks) fall short of universal computation, but are assumed to be capable of explaining cognition and development. I present a class of models that bridge two concepts from distant fields: combinatory logic (or, equivalently, lambda calculus) and RNA molecular biology. A set of basic RNA editing rules can make it possible to compute any computable function with identical algorithmic complexity to that of Turing machines. The models do not assume extraordinarily complex molecular machinery or any processes that radically differ from what we already know to occur in cells. Distinct independent enzymes can mediate each of the rules and RNA molecules solve the problem of parenthesis matching through their secondary structure. In the most plausible of these models all of the editing rules can be implemented with merely cleavage and ligation operations at fixed positions relative to predefined motifs. This demonstrates that universal computation is well within the reach of molecular biology. It is therefore reasonable to assume that life has evolved - or possibly began with - a universal computer that yet remains to be discovered. The variety of seemingly unrelated computational problems across many scales can potentially be solved using the same RNA-based computation system. Experimental validation of this theory may immensely impact our understanding of memory, cognition, development, disease, evolution, and the early stages of life.


title: A Generalized Framework for Domain Adaptation of PLDA in Speaker Recognition

abstract: This paper proposes a generalized framework for domain adaptation of Probabilistic Linear Discriminant Analysis (PLDA) in speaker recognition. It not only includes several existing supervised and unsupervised domain adaptation methods but also makes possible more flexible usage of available data in different domains. In particular, we introduce here the two new techniques described below. (1) Correlation-alignment-based interpolation and (2) covariance regularization. The proposed correlation-alignment-based interpolation method decreases minCprimary up to 30.5% as compared with that from an out-of-domain PLDA model before adaptation, and minCprimary is also 5.5% lower than with a conventional linear interpolation method with optimal interpolation weights. Further, the proposed regularization technique ensures robustness in interpolations w.r.t. varying interpolation weights, which in practice is essential.


title: The high helium abundance and charge states of the interplanetary CME and its material source on the Sun

abstract: Identifying the source of the material within coronal mass ejections (CMEs) and understanding CME onset mechanisms are fundamental issues in solar and space physics. Parameters relating to plasma composition, such as charge states and He abundance (\ahe), may be different for plasmas originating from differing processes or regions on the Sun. Thus, it is crucial to examine the relationship between in-situ measurements of CME composition and activity on the Sun. We study the CME that erupted on 2014 September 10, in association with an X1.6 flare, by analyzing AIA imaging and IRIS spectroscopic observations and its in-situ signatures detected by Wind and ACE. We find that during the slow expansion and intensity increase of the sigmoid, plasma temperatures of 9 MK, and higher, first appear at the footpoints of the sigmoid, associated with chromospheric brightening. Then the high-temperature region extends along the sigmoid. IRIS observations confirm that this extension is caused by transportation of hot plasma upflow. Our results show that chromospheric material can be heated to 9 MK, and above, by chromospheric evaporation at the sigmoid footpoints before flare onset. The heated chromospheric material can transport into the sigmoidal structure and supply mass to the CME. The aforementioned CME mass supply scenario provides a reasonable explanation for the detection of high charge states and elevated \ahe\ in the associated ICME. The observations also demonstrate that the quasi-steady evolution in the precursor phase is dominated by magnetic reconnection between the rising flux rope and the overlying magnetic field structure.


title: Grasping Detection Network with Uncertainty Estimation for Confidence-Driven Semi-Supervised Domain Adaptation

abstract: Data-efficient domain adaptation with only a few labelled data is desired for many robotic applications, e.g., in grasping detection, the inference skill learned from a grasping dataset is not universal enough to directly apply on various other daily/industrial applications. This paper presents an approach enabling the easy domain adaptation through a novel grasping detection network with confidence-driven semi-supervised learning, where these two components deeply interact with each other. The proposed grasping detection network specially provides a prediction uncertainty estimation mechanism by leveraging on Feature Pyramid Network (FPN), and the mean-teacher semi-supervised learning utilizes such uncertainty information to emphasizing the consistency loss only for those unlabelled data with high confidence, which we referred it as the confidence-driven mean teacher. This approach largely prevents the student model to learn the incorrect/harmful information from the consistency loss, which speeds up the learning progress and improves the model accuracy. Our results show that the proposed network can achieve high success rate on the Cornell grasping dataset, and for domain adaptation with very limited data, the confidence-driven mean teacher outperforms the original mean teacher and direct training by more than 10% in evaluation loss especially for avoiding the overfitting and model diverging.


title: Ensemble learning reveals dissimilarity between rare-earth transition metal binary alloys with respect to the Curie temperature

abstract: We propose a data-driven method to extract dissimilarity between materials, with respect to a given target physical property. The technique is based on an ensemble method with Kernel ridge regression as the predicting model; multiple random subset sampling of the materials is done to generate prediction models and the corresponding contributions of the reference training materials in detail. The distribution of the predicted values for each material can be approximated by a Gaussian mixture model. The reference training materials contributed to the prediction model that accurately predicts the physical property value of a specific material, are considered to be similar to that material, or vice versa. Evaluations using synthesized data demonstrate that the proposed method can effectively measure the dissimilarity between data instances. An application of the analysis method on the data of Curie temperature (TC) of binary 3d transition metal 4f rare earth binary alloys also reveals meaningful results on the relations between the materials. The proposed method can be considered as a potential tool for obtaining a deeper understanding of the structure of data, with respect to a target property, in particular.


title: Exclusive $B^0 \to \pi^- \ell^+ \nu_\ell$ Decays with Hadronic Full Event Interpretation Tagging in 34.6 fb$^{-1}$ of Belle II Data

abstract: We present the results of the re-discovery of the decay $B^0 \to \pi^- \ell^+ \nu_\ell$ in 34.6 fb$^{-1}$ of Belle II data using hadronic $B$-tagging via the Full Event Interpretation algorithm. We observe 21 signal events on a background of 155 in a fit to the distribution of the square of the missing mass, $M_{\mathrm{miss}}^2$, with a significance of 5.69$\sigma$, and determine a total branching fraction of (1.58 $\pm$ 0.43$_{\mathrm{stat}}$ $\pm$ 0.07$_{\mathrm{sys}}$) $\times 10^{-4}$.


title: A realization result for systems of sets of lengths

abstract: Let $\mathcal L^*$ be a family of finite subsets of $\mathbb N_0$ having the following properties. (a). $\{0\}, \{1\} \in \mathcal L^*$ and all other sets of $\mathcal L^*$ lie in $\mathbb N_{\ge 2}$. (b). If $L_1, L_2 \in \mathcal L^*$, then the sumset $L_1 + L_2 \in \mathcal L^*$. We show that there is a Dedekind domain $D$ whose system of sets of lengths equals $\mathcal L^*$.


title: VAIM: Visual Analytics for Influence Maximization

abstract: In social networks, individuals' decisions are strongly influenced by recommendations from their friends and acquaintances. The influence maximization (IM) problem asks to select a seed set of users that maximizes the influence spread, i.e., the expected number of users influenced through a stochastic diffusion process triggered by the seeds. In this paper, we present VAIM, a visual analytics system that supports users in analyzing the information diffusion process determined by different IM algorithms. By using VAIM one can: (i) simulate the information spread for a given seed set on a large network, (ii) analyze and compare the effectiveness of different seed sets, and (iii) modify the seed sets to improve the corresponding influence spread.


title: A Simple and Fast Algorithm for Computing the $N$-th Term of a Linearly Recurrent Sequence

abstract: We present a simple and fast algorithm for computing the $N$-th term of a given linearly recurrent sequence. Our new algorithm uses $O(\mathsf{M}(d) \log N)$ arithmetic operations, where $d$ is the order of the recurrence, and $\mathsf{M}(d)$ denotes the number of arithmetic operations for computing the product of two polynomials of degree $d$. The state-of-the-art algorithm, due to Charles Fiduccia (1985), has the same arithmetic complexity up to a constant factor. Our algorithm is simpler, faster and obtained by a totally different method. We also discuss several algorithmic applications, notably to polynomial modular exponentiation, powering of matrices and high-order lifting.


title: DronePose: Photorealistic UAV-Assistant Dataset Synthesis for 3D Pose Estimation via a Smooth Silhouette Loss

abstract: In this work we consider UAVs as cooperative agents supporting human users in their operations. In this context, the 3D localisation of the UAV assistant is an important task that can facilitate the exchange of spatial information between the user and the UAV. To address this in a data-driven manner, we design a data synthesis pipeline to create a realistic multimodal dataset that includes both the exocentric user view, and the egocentric UAV view. We then exploit the joint availability of photorealistic and synthesized inputs to train a single-shot monocular pose estimation model. During training we leverage differentiable rendering to supplement a state-of-the-art direct regression objective with a novel smooth silhouette loss. Our results demonstrate its qualitative and quantitative performance gains over traditional silhouette objectives. Our data and code are available at https://vcl3d.github.io/DronePose


title: Unified Rules of Renewable Weighted Sums for Various Online Updating Estimations

abstract: This paper establishes unified frameworks of renewable weighted sums (RWS) for various online updating estimations in the models with streaming data sets. The newly defined RWS lays the foundation of online updating likelihood, online updating loss function, online updating estimating equation and so on. The idea of RWS is intuitive and heuristic, and the algorithm is computationally simple. This paper chooses nonparametric model as an exemplary setting. The RWS applies to various types of nonparametric estimators, which include but are not limited to nonparametric likelihood, quasi-likelihood and least squares. Furthermore, the method and the theory can be extended into the models with both parameter and nonparametric function. The estimation consistency and asymptotic normality of the proposed renewable estimator are established, and the oracle property is obtained. Moreover, these properties are always satisfied, without any constraint on the number of data batches, which means that the new method is adaptive to the situation where streaming data sets arrive perpetually. The behavior of the method is further illustrated by various numerical examples from simulation experiments and real data analysis.


title: Efficient and Accurate Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem for Crystalline Systems

abstract: Optical properties of materials related to light absorption and scattering are explained by the excitation of electrons. The Bethe-Salpeter equation is the state-of-the-art approach to describe these processes from first principles (ab initio), i.e. without the need for empirical data in the model. To harness the predictive power of the equation, it is mapped to an eigenvalue problem via an appropriate discretization scheme. The eigenpairs of the resulting large, dense, structured matrix can be used to compute dielectric properties of the considered crystalline or molecular system. The matrix always shows a $2\times 2$ block structure. Additionally, certain definiteness properties typically hold. One form can be acquired for crystalline systems, another one is more general and can for example be used to study molecules. In this work, we present new theoretical results characterizing the structure of the two forms in the language of non-standard scalar products. These results enable us to develop a new perspective on the state-of-the-art solution approach for crystalline systems. This new viewpoint is used to develop two new methods for solving the eigenvalue problem. One requires less computational effort while providing the same degree of accuracy. The other one improves the expected accuracy, compared to methods currently in use, with a comparable performance. Both methods are well suited for high performance environments and only rely on basic numerical linear algebra building blocks.


title: Simultaneous Detection and Tracking with Motion Modelling for Multiple Object Tracking

abstract: Deep learning-based Multiple Object Tracking (MOT) currently relies on off-the-shelf detectors for tracking-by-detection.This results in deep models that are detector biased and evaluations that are detector influenced. To resolve this issue, we introduce Deep Motion Modeling Network (DMM-Net) that can estimate multiple objects' motion parameters to perform joint detection and association in an end-to-end manner. DMM-Net models object features over multiple frames and simultaneously infers object classes, visibility, and their motion parameters. These outputs are readily used to update the tracklets for efficient MOT. DMM-Net achieves PR-MOTA score of 12.80 @ 120+ fps for the popular UA-DETRAC challenge, which is better performance and orders of magnitude faster. We also contribute a synthetic large-scale public dataset Omni-MOT for vehicle tracking that provides precise ground-truth annotations to eliminate the detector influence in MOT evaluation. This 14M+ frames dataset is extendable with our public script (Code at Dataset <https://github.com/shijieS/OmniMOTDataset>, Dataset Recorder <https://github.com/shijieS/OMOTDRecorder>, Omni-MOT Source <https://github.com/shijieS/DMMN>). We demonstrate the suitability of Omni-MOT for deep learning with DMMNet and also make the source code of our network public.


title: Plane Spanning Trees in Edge-Colored Simple Drawings of $K_n$

abstract: K\'{a}rolyi, Pach, and T\'{o}th proved that every 2-edge-colored straight-line drawing of the complete graph contains a monochromatic plane spanning tree. It is open if this statement generalizes to other classes of drawings, specifically, to simple drawings of the complete graph. These are drawings where edges are represented by Jordan arcs, any two of which intersect at most once. We present two partial results towards such a generalization. First, we show that the statement holds for cylindrical simple drawings. (In a cylindrical drawing, all vertices are placed on two concentric circles and no edge crosses either circle.) Second, we introduce a relaxation of the problem in which the graph is $k$-edge-colored, and the target structure must be hypochromatic, that is, avoid (at least) one color class. In this setting, we show that every $\lceil (n+5)/6\rceil$-edge-colored monotone simple drawing of $K_n$ contains a hypochromatic plane spanning tree. (In a monotone drawing, every edge is represented as an $x$-monotone curve.)


title: On the Use of Quasiorders in Formal Language Theory

abstract: In this thesis we use quasiorders on words to offer a new perspective on two well-studied problems from Formal Language Theory: deciding language inclusion and manipulating the finite automata representations of regular languages. First, we present a generic quasiorder-based framework that, when instantiated with different quasiorders, yields different algorithms (some of them new) for deciding language inclusion. We then instantiate this framework to devise an efficient algorithm for searching with regular expressions on grammar-compressed text. Finally, we define a framework of quasiorder-based automata constructions to offer a new perspective on residual automata.


title: Basis divisors and balanced metrics

abstract: Using log canonical thresholds and basis divisors Fujita--Odaka introduced purely algebro-geometric invariants $\delta_m$ whose limit in $m$ is now known to characterize uniform K-stability on a Fano variety. As shown by Blum-Jonsson this carries over to a general polarization, and together with work of Berman, Boucksom, and Jonsson, it is now known that the limit of these $\delta_m$-invariants characterizes uniform Ding stability. A basic question since Fujita-Odaka's work has been to find an analytic interpretation of these invariants. We show that each $\delta_m$ is the coercivity threshold of a quantized Ding functional on the $m$-th Bergman space and thus characterizes the existence of balanced metrics. This approach has a number of applications. The most basic one is that it provides an alternative way to compute these invariants, which is new even for $\mathbb{P}^n$. Second, it allows us to introduce algebraically defined invariants that characterize the existence of K\"ahler-Ricci solitons (and the more general $g$-solitons of Berman-Witt Nystr\"om), as well as coupled versions thereof. Third, it leads to approximation results involving balanced metrics in the presence of automorphisms that extend some results of Donaldson.


title: Optimal Load Balancing in Bipartite Graphs

abstract: Applications in cloud platforms motivate the study of efficient load balancing under job-server constraints and server heterogeneity. In this paper, we study load balancing on a bipartite graph where left nodes correspond to job types and right nodes correspond to servers, with each edge indicating that a job type can be served by a server. Thus edges represent locality constraints, i.e., each job can only be served at servers which contained certain data and/or machine learning (ML) models. Servers in this system can have heterogeneous service rates. In this setting, we investigate the performance of two policies named Join-the-Fastest-of-the-Shortest-Queue (JFSQ) and Join-the-Fastest-of-the-Idle-Queue (JFIQ), which are simple variants of Join-the-Shortest-Queue and Join-the-Idle-Queue, where ties are broken in favor of the fastest servers. Under a "well-connected" graph condition, we show that JFSQ and JFIQ are asymptotically optimal in the mean response time when the number of servers goes to infinity. In addition to asymptotic optimality, we also obtain upper bounds on the mean response time for finite-size systems. We further show that the well-connectedness condition can be satisfied by a random bipartite graph construction with relatively sparse connectivity.


title: The MOSDEF Survey: Neon as a Probe of ISM Physical Conditions at High Redshift

abstract: We present results on the properties of neon emission in $z\sim2$ star-forming galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Doubly-ionized neon ([NeIII]3869) is detected at $\geq3\sigma$ in 61 galaxies, representing $\sim$25% of the MOSDEF sample with H$\alpha$, H$\beta$, and [OIII]$5007$ detections at similar redshifts. We consider the neon emission-line properties of both individual galaxies with [NeIII]3869 detections and composite $z\sim2$ spectra binned by stellar mass. With no requirement of [NeIII]3869 detection, the latter provide a more representative picture of neon emission-line properties in the MOSDEF sample. The [NeIII]3869/[OII]3727 ratio (Ne3O2) is anti-correlated with stellar mass in $z\sim2$ galaxies, as expected based on the mass-metallicity relation. It is also positively correlated with the [OIII]$5007$/[OII]$3727$ ratio (O32), but $z\sim2$ line ratios are offset towards higher Ne3O2 at fixed O32, compared with both local star-forming galaxies and individual H~II regions. Despite the offset towards higher Ne3O2 at fixed O32 at $z\sim2$, biases in inferred Ne3O2-based metallicity are small. Accordingly, Ne3O2 may serve as an important metallicity indicator deep into the reionization epoch. Analyzing additional rest-optical line ratios including [NeIII]$3869$/[OIII]$5007$ (Ne3O3) and [OIII]$5007$/H$\beta$ (O3H$\beta$), we conclude that the nebular emission-line ratios of $z\sim2$ star-forming galaxies suggest a harder ionizing spectrum (lower stellar metallicity, i.e., Fe/H) at fixed gas-phase oxygen abundance, compared to systems at $z\sim0$. These new results based on neon lend support to the physical picture painted by oxygen, nitrogen, hydrogen, and sulfur emission, of an ionized ISM in high-redshift star-forming galaxies irradiated by chemically young, $\alpha$-enhanced massive stars.


title: Period and Signal Reconstruction from the Curve of Sample-Sequences

abstract: A finite sequence of equidistant samples (a sample train) of a periodic signal can be identified with a point in a multi-dimensional space. Such a point depends on the sampled signal, the sampling period, and the starting time of the sequence. If the starting time varies, then the corresponding point moves along a closed curve. We prove that such a curve, i.e., the set of all sample trains of a given length, determines the period of the sampled signal, provided that the sampling period is known. This is true even if the trains are short, and if the samples comprising trains are taken at a sub-Nyquist rate. The presented result is proved with a help of the theory of rotation numbers developed by Poincar\'e. We also prove that the curve of sample trains determines the sampled signal up to a time shift, provided that the ratio of the sampling period to the period of the signal is irrational. Eventually, we give an example which shows that the assumption on incommensurability of the periods cannot be dropped.


title: Spectrum and pseudospectrum for quadratic polynomials in Ginibre matrices

abstract: For a fixed quadratic polynomial $\mathfrak{p}$ in $n$ non-commuting variables, and $n$ independent $N\times N$ complex Ginibre matrices $X_1^N,\dots, X_n^N$, we establish the convergence of the empirical spectral distribution of $P^N =\mathfrak{p}(X_1^N,\dots, X_n^N)$ to the Brown measure of $\mathfrak{p}$ evaluated at $n$ freely independent circular elements $c_1,\dots, c_n$ in a non-commutative probability space. The main step of the proof is to obtain quantitative control on the pseudospectrum of $P^N$. Via the well-known linearization trick this hinges on anti-concentration properties for certain matrix-valued random walks, which we find can fail for structural reasons of a different nature from the arithmetic obstructions that were illuminated in works on the Littlewood--Offord problem for discrete scalar random walks.


title: Validating inversions for toroidal flows using normal-mode coupling

abstract: Normal-mode coupling is a helioseismic technique that uses measurements of mode eigenfunctions to infer the interior structure of the Sun. This technique has led to insights into the evolution and structure of toroidal flows in the solar interior. Here, we validate an inversion algorithm for normal-mode coupling by generating synthetic seismic measurements associated with input flows and comparing the input and inverted velocities. We study four different cases of input toroidal flows and compute synthetics that take into account the partial visibility of the Sun. We invert the synthetics using Subtractive Optimally Localized Averages (SOLA) and also try to mitigate the systematics of mode leakage. We demonstrate that, ultimately, inversions are only as good as the model we assume for the correlation between flow velocities.


title: EGO-Planner: An ESDF-free Gradient-based Local Planner for Quadrotors

abstract: Gradient-based planners are widely used for quadrotor local planning, in which a Euclidean Signed Distance Field (ESDF) is crucial for evaluating gradient magnitude and direction. Nevertheless, computing such a field has much redundancy since the trajectory optimization procedure only covers a very limited subspace of the ESDF updating range. In this paper, an ESDF-free gradient-based planning framework is proposed, which significantly reduces computation time. The main improvement is that the collision term in the penalty function is formulated by comparing the colliding trajectory with a collision-free guiding path. The resulting obstacle information will be stored only if the trajectory hits new obstacles, making the planner only extract necessary obstacle information. Then, we lengthen the time allocation if dynamical feasibility is violated. An anisotropic curve fitting algorithm is introduced to adjust higher-order derivatives of the trajectory while maintaining the original shape. Benchmark comparisons and real-world experiments verify its robustness and high-performance. The source code is released as ROS packages.


title: Elimination of fast variables in stochastic nonlinear kinetics

abstract: A reduced chemical scheme involving a small number of variables is often sufficient to account for the deterministic evolution of the concentrations of the main species contributing to a reaction. However its predictions are questionable in small systems used for example in fluorescence correlation spectroscopy (FCS) or in explosive systems involving strong nonlinearities such as autocatalytic steps. We make precise dynamical criteria defining the validity domain of the quasi-steady-state approximation and the elimination of a fast concentration in deterministic dynamics. Designing two different three-variable models converging toward the same two-variable model, we show that the variances and covariance of the fluctuations of the slow variables are not correctly predicted by the two-variable model, even in the limit of a large system size. The more striking weaknesses of the reduced scheme are figured out in mesoscaled systems containing a small number of molecules. The results of two stochastic approaches are compared and the shortcomings of the Langevin equations with respect to the master equation are pointed out. We conclude that the description of the fluctuations and their coupling with nonlinearities of deterministic dynamics escape reduced chemical schemes.


title: Uncertainty Estimation in Medical Image Denoising with Bayesian Deep Image Prior

abstract: Uncertainty quantification in inverse medical imaging tasks with deep learning has received little attention. However, deep models trained on large data sets tend to hallucinate and create artifacts in the reconstructed output that are not anatomically present. We use a randomly initialized convolutional network as parameterization of the reconstructed image and perform gradient descent to match the observation, which is known as deep image prior. In this case, the reconstruction does not suffer from hallucinations as no prior training is performed. We extend this to a Bayesian approach with Monte Carlo dropout to quantify both aleatoric and epistemic uncertainty. The presented method is evaluated on the task of denoising different medical imaging modalities. The experimental results show that our approach yields well-calibrated uncertainty. That is, the predictive uncertainty correlates with the predictive error. This allows for reliable uncertainty estimates and can tackle the problem of hallucinations and artifacts in inverse medical imaging tasks.


title: Training Matters: Unlocking Potentials of Deeper Graph Convolutional Neural Networks

abstract: The performance limit of Graph Convolutional Networks (GCNs) and the fact that we cannot stack more of them to increase the performance, which we usually do for other deep learning paradigms, are pervasively thought to be caused by the limitations of the GCN layers, including insufficient expressive power, etc. However, if so, for a fixed architecture, it would be unlikely to lower the training difficulty and to improve performance by changing only the training procedure, which we show in this paper not only possible but possible in several ways. This paper first identify the training difficulty of GCNs from the perspective of graph signal energy loss. More specifically, we find that the loss of energy in the backward pass during training nullifies the learning of the layers closer to the input. Then, we propose several methodologies to mitigate the training problem by slightly modifying the GCN operator, from the energy perspective. After empirical validation, we confirm that these changes of operator lead to significant decrease in the training difficulties and notable performance boost, without changing the composition of parameters. With these, we conclude that the root cause of the problem is more likely the training difficulty than the others.


title: Delta-T noise in the Kondo regime

abstract: We study the delta-T noise in the Kondo regime, which implies the charge current noise under the temperature bias for the SU(2) Kondo quantum dot. We propose an experimentally measurable quantity to quantify the low-temperature properties in the delta-T noise: $S_{\ell}=S(T_{\mathrm{L}},T_{\mathrm{R}}) - (1/2)[S(T_{\mathrm{L}},T_{\mathrm{L}}) + S(T_{\mathrm{R}},T_{\mathrm{R}})]$, which yields the shot noise expression in the noninteracting limit. We calculate this quantity for the SU(2) Kondo quantum dot in the particle-hole symmetric case. We found that the $S_{\ell}$ exhibits qualitatively the same behavior in both the electrochemical potential biased case and the temperature biased case. The quantitative difference appears as a difference of the coefficients of the noises, which reflects the difference of the Fermi distribution function: electrochemical potential biased or temperature biased.


title: Image quality assessment for closed-loop computer-assisted lung ultrasound

abstract: We describe a novel, two-stage computer assistance system for lung anomaly detection using ultrasound imaging in the intensive care setting to improve operator performance and patient stratification during coronavirus pandemics. The proposed system consists of two deep-learning-based models: a quality assessment module that automates predictions of image quality, and a diagnosis assistance module that determines the likelihood-oh-anomaly in ultrasound images of sufficient quality. Our two-stage strategy uses a novelty detection algorithm to address the lack of control cases available for training the quality assessment classifier. The diagnosis assistance module can then be trained with data that are deemed of sufficient quality, guaranteed by the closed-loop feedback mechanism from the quality assessment module. Using more than 25000 ultrasound images from 37 COVID-19-positive patients scanned at two hospitals, plus 12 control cases, this study demonstrates the feasibility of using the proposed machine learning approach. We report an accuracy of 86% when classifying between sufficient and insufficient quality images by the quality assessment module. For data of sufficient quality - as determined by the quality assessment module - the mean classification accuracy, sensitivity, and specificity in detecting COVID-19-positive cases were 0.95, 0.91, and 0.97, respectively, across five holdout test data sets unseen during the training of any networks within the proposed system. Overall, the integration of the two modules yields accurate, fast, and practical acquisition guidance and diagnostic assistance for patients with suspected respiratory conditions at point-of-care.


title: Distributions of Bubble Lifetimes and Bubble Lengths in DNA

abstract: We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these distributions. For both models we find that base pair-dependent (GC vs AT) thresholds for considering complementary nucleotides to be separated are able to reproduce the observed dependence of the melting temperature on the GC content of the DNA sequence. Using these thresholds for base pair openings, we obtain bubble lifetime distributions for bubbles of lengths up to ten base pairs as the GC content of the sequences is varied, which are accurately fitted with stretched exponential functions. We find that for both models the average bubble lifetime decreases with increasing either the bubble length or the GC content. In addition, the obtained bubble length distributions are also fitted by appropriate stretched exponential functions and our results show that short bubbles have similar likelihoods for any GC content, but longer ones are substantially more likely to occur in AT-rich sequences. We also show that the ePBD model permits more, longer-lived, bubbles than the PBD system.

categories: physics.bio-ph cond-mat.stat-mech physics.comp-ph q-bio.BM


title: Encounters involving planetary systems in birth environments: the significant role of binaries

abstract: Most stars form in a clustered environment. Both single and binary stars will sometimes encounter planetary systems in such crowded environments. Encounter rates for binaries may be larger than for single stars, even for binary fractions as low as 10-20 per cent. In this work, we investigate scatterings between a Sun-Jupiter pair and both binary and single stars as in young clusters. We first perform a set of simulations of encounters involving wide ranges of binaries and single stars, finding that wider binaries have larger cross sections for the planet's ejection. Secondly, we consider such scatterings in a realistic population, drawing parameters for the binaries and single stars from the observed population. The scattering outcomes are diverse, including ejection, capture/exchange and collision. The binaries are more effective than single stars by a factor of several or more in causing the planet's ejection and collision. Hence, in a cluster, as long as the binary fraction is larger than about 10 per cent, the binaries will dominate the scatterings in terms of these two outcomes. For an open cluster of a stellar density 50 pc$^{-3}$, a lifetime 100 Myr and a binary fraction 0.5, we estimate that of the order of 1 per cent of the Jupiters are ejected, 0.1 per cent collide with a star, 0.1 per cent change ownership and 10 per cent of the Sun-Jupiter pairs acquire a stellar companion during scatterings. These companions are typically 1000s of au distant and in half of the cases (so 5 per cent of all Sun-Jupiter pairs), they can excite the planet's orbit through Kozai--Lidov mechanism before stripped by later encounters. Our result suggests that the Solar System may have once had a companion in its birth cluster.


title: Plenty of big projections imply big pieces of Lipschitz graphs

abstract: I prove that a closed $n$-regular set $E \subset \mathbb{R}^{d}$ with plenty of big projections has big pieces of Lipschitz graphs. This answers a question of David and Semmes.


title: Complete the Missing Half: Augmenting Aggregation Filtering with Diversification for Graph Convolutional Networks

abstract: The core operation of current Graph Neural Networks (GNNs) is the aggregation enabled by the graph Laplacian or message passing, which filters the neighborhood node information. Though effective for various tasks, in this paper, we show that they are potentially a problematic factor underlying all GNN methods for learning on certain datasets, as they force the node representations similar, making the nodes gradually lose their identity and become indistinguishable. Hence, we augment the aggregation operations with their dual, i.e. diversification operators that make the node more distinct and preserve the identity. Such augmentation replaces the aggregation with a two-channel filtering process that, in theory, is beneficial for enriching the node representations. In practice, the proposed two-channel filters can be easily patched on existing GNN methods with diverse training strategies, including spectral and spatial (message passing) methods. In the experiments, we observe desired characteristics of the models and significant performance boost upon the baselines on 9 node classification tasks.


title: Real-Time Time-Dependent Density Functional Theory within FHI-aims

abstract: Real-Time Time-Dependent Density Functional Theory (TDDFT) has become an attractive tool to model quantum dynamics on a first-principles Density Functional Theory level. In recent years, several developments and applications in this field were achieved and hopefully lead to new insights. We present here our versatile and efficient Real-Time TDDFT implementation into the all-electron numerical basis-set DFT code package FHI-aims. This article is meant as a short overview on how we performed this task and what can be done with our implementation. We further shed light on the connection of the basis set size to the accuracy of absorption spectrum simulation results.


title: Spectral analysis for a multi-dimensional split-step quantum walk with a defect

abstract: This paper studies the spectrum of a multi-dimensional split-step quantum walk with a defect that cannot be analysed in the previous papers. To this end, we have developed a new technique which allow us to use a spectral mapping theorem for the one-defect model. We also derive the time-averaged limit measure for one-dimensional case as an application of the spectral analysis.


title: Yet Another Intermediate-Level Attack

abstract: The transferability of adversarial examples across deep neural network (DNN) models is the crux of a spectrum of black-box attacks. In this paper, we propose a novel method to enhance the black-box transferability of baseline adversarial examples. By establishing a linear mapping of the intermediate-level discrepancies (between a set of adversarial inputs and their benign counterparts) for predicting the evoked adversarial loss, we aim to take full advantage of the optimization procedure of multi-step baseline attacks. We conducted extensive experiments to verify the effectiveness of our method on CIFAR-100 and ImageNet. Experimental results demonstrate that it outperforms previous state-of-the-arts considerably. Our code is at https://github.com/qizhangli/ila-plus-plus.


title: Using Atmospheric Impact Data to Model Meteoroid Close Encounters

abstract: Based on telescopic observations of Jupiter-family comets (JFCs), there is predicted to be a paucity of objects at sub-kilometre sizes. However, several bright fireballs and some meteorites have been tenuously linked to the JFC population, showing metre-scale objects do exist in this region. In 2017, the Desert Fireball Network (DFN) observed a grazing fireball that redirected a meteoroid from an Apollo-type orbit to a JFC-like orbit. Using orbital data collected by the DFN, in this study, we have generated an artificial dataset of close terrestrial encounters that come within $1.5$ lunar distances (LD) of the Earth in the size-range of $0.01-100$kg. This range of objects is typically too small for telescopic surveys to detect, so using atmospheric impact flux data from fireball observations is currently one of the only ways to characterise these close encounters. Based on this model, we predict that within the considered size-range $2.5\times 10^{8}$ objects ($0.1\%$ of the total flux) from asteroidal orbits ($T_{J}>3$) are annually sent onto JFC-like orbits ($2<T_{J}<3$), with a steady-state population of about $8\times 10^{13}$ objects. Close encounters with the Earth provide another way to transfer material to the JFC region. Additionally, using our model, we found that approximately $1.96\times 10^{7}$ objects are sent onto Aten-type orbits and $\sim10^{4}$ objects are ejected from the Solar System annually via a close encounter with the Earth.


title: The Curse of Shared Knowledge: Recursive Belief Reasoning in a Coordination Game with Imperfect Information

abstract: Common knowledge is a necessary condition for safe group coordination. When common knowledge can not be obtained, humans routinely use their ability to attribute beliefs and intentions in order to infer what is known. But such shared knowledge attributions are limited in depth and therefore prone to coordination failures, because any finite-order knowledge attribution allows for an even higher order attribution that may change what is known by whom. In three separate experiments we investigate to which degree human participants (N=802) are able to recognize the difference between common knowledge and nth-order shared knowledge. We use a new two-person coordination game with imperfect information that is able to cast the recursive game structure and higher-order uncertainties into a simple, everyday-like setting. Our results show that participants have a very hard time accepting the fact that common knowledge is not reducible to shared knowledge. Instead, participants try to coordinate even at the shallowest depths of shared knowledge and in spite of huge payoff penalties.


title: Isotopically distinct terrestrial planets via local accretion

abstract: Combining isotopic constraints from meteorite data with dynamical models of planet formation proves to be advantageous in identifying the best model for terrestrial planet formation. Prior studies have shown that the probability of reproducing the distinct isotopic compositions of the Earth and Mars for both classical and Grand Tack models is very low. In the framework of the Grand Tack model, for Mars to be isotopically different from the Earth, it had to form under very specific conditions. Here, we subjected a fairly new and unexplored model--the depleted disc model--to the test. It presupposes that the region in the inner protoplanetary disc from Mars' orbit and beyond is depleted in mass such that Mars is left with insufficient material to grow to a larger size. Our aim is to test the whether the distinct isotopic compositions of the Earth and Mars are a natural outcome of this model. We found that the terrestrial planets accrete material mostly locally and have feeding zones that are sufficiently distinct. The Earth and Mars, and by extension, Venus, can have distinct isotopic compositions if there is an isotopic gradient in the terrestrial planet region of the protoplanetary disc. Our results suggest that the material in the inner Solar System most likely did not undergo substantial mixing that homogenised the potential isotopic gradient, in contrast to the Grand Tack model where the feeding zones of the terrestrial planets are nearly identical due to the mixing of material by Jupiter's migration.


title: Decay of streaks and rolls in plane Couette-Poiseuille flow

abstract: We report the results of an experimental investigation into the decay of turbulence in plane Couette-Poiseuille flow using 'quench' experiments where the flow laminarises after a sudden reduction in Reynolds number $Re$. Specifically, we study the velocity field in the streamwise-spanwise plane. We show that the spanwise velocity containing rolls, decays faster than the streamwise velocity, which displays elongated regions of higher or lower velocity called streaks. At final Reynolds numbers above 425, the decay of streaks displays two stages: first a slow decay when rolls are present and secondly a more rapid decay of streaks alone. The difference in behaviour results from the regeneration of streaks by rolls, called the lift-up effect. We define the turbulent fraction as the portion of the flow containing turbulence and this is estimated by thresholding the spanwise velocity component. It decreases linearly with time in the whole range of final $Re$. The corresponding decay slope increases linearly with final $Re$. The extrapolated value at which this decay slope vanishes is $Re_{a_z}\approx 656\pm10$, close to $Re_g\approx 670$ at which turbulence is self-sustained. The decay of the energy computed from the spanwise velocity component is found to be exponential. The corresponding decay rate increases linearly with $Re$, with an extrapolated vanishing value at $Re_{A_z}\approx 688\pm10$. This value is also close to the value at which the turbulence is self-sustained, showing that valuable information on the transition can be obtained over a wide range of $Re$.


title: On the Kodaira dimension of the moduli space of curves of genus 16

abstract: We show that the moduli space of curves of genus 16 is NOT of general type.


title: Relaxation of the Boussinesq system and applications to the Rayleigh-Taylor instability

abstract: We consider the evolution of two incompressible fluids with homogeneous densities $\rho_-<\rho_+$ subject to gravity described by the inviscid Boussinesq equations and provide the explicit relaxation of the associated differential inclusion. The existence of a subsolution to the relaxation allows one to conclude the existence of turbulently mixing solutions to the original Boussinesq system. As a specific application we investigate subsolutions emanating from the classical Rayleigh-Taylor initial configuration where the two fluids are separated by a horizontal interface with the heavier fluid being on top of the lighter. It turns out that among all self-similar subsolutions the criterion of maximal initial energy dissipation selects a linear density profile and a quadratic growth of the mixing zone. The subsolution selected this way can be extended in an admissible way to exist for all times. We provide two possible extensions with different long-time limits. The first one corresponds to a total mixture of the two fluids, the second corresponds to a full separation with the lighter fluid on top of the heavier. There is no motion in either of the limit states.


title: Checkworthiness in Automatic Claim Detection Models: Definitions and Analysis of Datasets

abstract: Public, professional and academic interest in automated fact-checking has drastically increased over the past decade, with many aiming to automate one of the first steps in a fact-check procedure: the selection of so-called checkworthy claims. However, there is little agreement on the definition and characteristics of checkworthiness among fact-checkers, which is consequently reflected in the datasets used for training and testing checkworthy claim detection models. After elaborate analysis of checkworthy claim selection procedures in fact-check organisations and analysis of state-of-the-art claim detection datasets, checkworthiness is defined as the concept of having a spatiotemporal and context-dependent worth and need to have the correctness of the objectivity it conveys verified. This is irrespective of the claim's perceived veracity judgement by an individual based on prior knowledge and beliefs. Concerning the characteristics of current datasets, it is argued that the data is not only highly imbalanced and noisy, but also too limited in scope and language. Furthermore, we believe that the subjective concept of checkworthiness might not be a suitable filter for claim detection.


title: Magnetic Two-Dimensional Chromium Trihalides: A Theoretical Perspective

abstract: The discovery of ferromagnetic order in monolayer 2D crystals has opened a new venue in the field of two dimensional (2D) materials. 2D magnets are not only interesting on their own, but their integration in van der Waals heterostructures allows for the observation of new and exotic effects in the ultrathin limit. The family of Chromium trihalides, CrI$_3$, CrBr$_3$ and CrCl$_3$, is, so far, the most studied among magnetic 2D crystals. In this mini-review, we provide a perspective of the state of the art of the theoretical understanding of magnetic 2D trihalides, most of which will also be relevant for other 2D magnets, such as vanadium trihalides. We discuss both the well-established facts, such as the origin of the magnetic moment and magnetic anisotropy and address as well open issues such as the nature of the anisotropic spin couplings and the magnitude of the magnon gap. Recent theoretical predictions on Moir\' e magnets and magnetic skyrmions are also discussed. Finally, we give some prospects about the future interest of these materials and possible device applications.


title: Fast recognition of some parametric graph families

abstract: We identify all $[1, \lambda, 8]$-cycle regular $I$-graphs and all $[1, \lambda, 8]$-cycle regular double generalized Petersen graphs. As a consequence we describe linear recognition algorithms for these graph families. Using structural properties of folded cubes we devise a $o(N \log N)$ recognition algorithm for them. We also study their $[1,\lambda,4]$, $[1,\lambda,6]$ and $[2, \lambda, 6]$-cycle regularity and settle the value of parameter $\lambda$.


title: Simple Analysis of Johnson-Lindenstrauss Transform under Neuroscience Constraints

abstract: The paper re-analyzes a version of the celebrated Johnson-Lindenstrauss Lemma, in which matrices are subjected to constraints that naturally emerge from neuroscience applications: a) sparsity and b) sign-consistency. This particular variant was studied first by Allen-Zhu, Gelashvili, Micali, Shavit and more recently by Jagadeesan (RANDOM'19). The contribution of this work is a novel proof, which in contrast to previous works a) uses the modern probability toolkit, particularly basics of sub-gaussian and sub-gamma estimates b) is self-contained, with no dependencies on subtle third-party results c) offers explicit constants. At the heart of our proof is a novel variant of Hanson-Wright Lemma (on concentration of quadratic forms). Of independent interest are also auxiliary facts on sub-gaussian random variables.


title: On a class of random sets of positive integers

abstract: In this note, we study a class of random subsets of positive integers induced by Bernoulli random variables. We obtain sufficient conditions such that the random set is almost surely lacunary, does not have bounded gaps and contains infinitely many arithmetic progressions, respectively.


title: Model-based Automated Testing of Mobile Applications: An Industrial Case Study

abstract: Automatic testing of mobile applications has been a well-researched area in recent years. However, testing in industry is still a very manual practice, as research results have not been fully transferred and adopted. Considering mobile applications, manual testing has the additional burden of adequate testing posed by a large number of available devices and different configurations, as well as the maintenance and setup of such devices. In this paper, we propose and evaluate the use of a model-based test generation approach, where generated tests are executed on a set of cloud-hosted real mobile devices. By using a model-based approach we generate dynamic, less brittle, and implementation simple test cases. The test execution on multiple real devices with different configurations increase the confidence in the implementation of the system under test. Our evaluation shows that the used approach produces a high coverage of the parts of the application related to user interactions. Nevertheless, the inclusion of external services in test generation is required in order to additionally increase the coverage of the complete application. Furthermore, we present the lessons learned while transferring and implementing this approach in an industrial context and applying it to the real product.


title: Global existence and spatial analyticity for a nonlocal flux with fractional diffusion

abstract: In this paper, we study a one dimensional nonlinear equation with diffusion $-\nu(-\partial_{xx})^{\frac{\alpha}{2}}$ for $0\leq \alpha\leq 2$ and $\nu>0$. We use a viscous-splitting algorithm to obtain global nonnegative weak solutions in space $L^1(\mathbb{R})\cap H^{1/2}(\mathbb{R})$ when $0\leq\alpha\leq 2$. For subcritical $1<\alpha\leq 2$ and critical case $\alpha=1$, we obtain global existence and uniqueness of nonnegative spatial analytic solutions. We use a fractional bootstrap method to improve the regularity of mild solutions in Bessel potential spaces for subcritical case $1<\alpha\leq 2$. Then, we show that the solutions are spatial analytic and can be extended globally. For the critical case $\alpha=1$, if the initial data $\rho_0$ satisfies $-\nu<\inf\rho_0<0$, we use the characteristics methods for complex Burgers equation to obtain a unique spatial analytic solution to our target equation in some bounded time interval. If $\rho_0\geq0$, the solution exists globally and converges to steady state.


title: The SuperFGD Prototype Charged Particle Beam Tests

abstract: A novel scintillator detector, the SuperFGD, has been selected as the main neutrino target for an upgrade of the T2K experiment ND280 near detector. The detector design will allow nearly 4{\pi} coverage for neutrino interactions at the near detector and will provide lower energy thresholds, significantly reducing systematic errors for the experiment. The SuperFGD is made of optically-isolated scintillator cubes of size 10x10x10 mm^3, providing the required spatial and energy resolution to reduce systematic uncertainties for future T2K runs. The SuperFGD for T2K will have close to two million cubes in a 1920x560x1840 mm^3 volume. A prototype made of 24x8x48 cubes was tested at a charged particle beamline at the CERN PS facility. The SuperFGD Prototype was instrumented with readout electronics similar to the future implementation for T2K. Results on electronics and detector response are reported in this paper, along with a discussion of the 3D reconstruction capabilities of this type of detector. Several physics analyses with the prototype data are also discussed, including a study of stopping protons.


title: Unsupervised Learning Facial Parameter Regressor for Action Unit Intensity Estimation via Differentiable Renderer

abstract: Facial action unit (AU) intensity is an index to describe all visually discernible facial movements. Most existing methods learn intensity estimator with limited AU data, while they lack generalization ability out of the dataset. In this paper, we present a framework to predict the facial parameters (including identity parameters and AU parameters) based on a bone-driven face model (BDFM) under different views. The proposed framework consists of a feature extractor, a generator, and a facial parameter regressor. The regressor can fit the physical meaning parameters of the BDFM from a single face image with the help of the generator, which maps the facial parameters to the game-face images as a differentiable renderer. Besides, identity loss, loopback loss, and adversarial loss can improve the regressive results. Quantitative evaluations are performed on two public databases BP4D and DISFA, which demonstrates that the proposed method can achieve comparable or better performance than the state-of-the-art methods. What's more, the qualitative results also demonstrate the validity of our method in the wild.


title: On measuring divergence for magnetic field modeling

abstract: A physical magnetic field has a divergence of zero. Numerical error in constructing a model field and computing the divergence, however, introduces a finite divergence into these calculations. A popular metric for measuring divergence is the average fractional flux $\langle |f_{i}| \rangle$. We show that $\langle |f_{i}| \rangle$ scales with the size of the computational mesh, and may be a poor measure of divergence because it becomes arbitrarily small for increasing mesh resolution, without the divergence actually decreasing. We define a modified version of this metric that does not scale with mesh size. We apply the new metric to the results of DeRosa et al. (2015), who measured $\langle |f_{i}| \rangle$ for a series of Nonlinear Force-Free Field (NLFFF) models of the coronal magnetic field based on solar boundary data binned at different spatial resolutions. We compute a number of divergence metrics for the DeRosa et al. (2015) data and analyze the effect of spatial resolution on these metrics using a non-parametric method. We find that some of the trends reported by DeRosa et al. (2015) are due to the intrinsic scaling of $\langle |f_{i}| \rangle$. We also find that different metrics give different results for the same data set and therefore there is value in measuring divergence via several metrics.


title: Bigrassmannian permutations and Verma modules

abstract: We show that bigrassmannian permutations determine the socle of the cokernel of an inclusion of Verma modules in type $A$. All such socular constituents turn out to be indexed by Weyl group elements from the penultimate two-sided cell. Combinatorially, the socular constituents in the cokernel of the inclusion of a Verma module indexed by $w\in S_n$ into the dominant Verma module are shown to be determined by the essential set of $w$ and their degrees in the graded picture are shown to be computable in terms of the associated rank function. As an application, we compute the first extension from a simple module to a Verma module.


title: Using Multi-Resolution Feature Maps with Convolutional Neural Networks for Anti-Spoofing in ASV

abstract: This paper presents a simple but effective method that uses multi-resolution feature maps with convolutional neural networks (CNNs) for anti-spoofing in automatic speaker verification (ASV). The central idea is to alleviate the problem that the feature maps commonly used in anti-spoofing networks are insufficient for building discriminative representations of audio segments, as they are often extracted by a single-length sliding window. Resulting trade-offs between time and frequency resolutions restrict the information in single spectrograms. The proposed method improves both frequency resolution and time resolution by stacking multiple spectrograms that are extracted using different window lengths. These are fed into a convolutional neural network in the form of multiple channels, making it possible to extract more information from input signals while only marginally increasing computational costs. The efficiency of the proposed method has been conformed on the ASVspoof 2019 database. We show that the use of the proposed multiresolution inputs consistently outperforms that of score fusion across different CNN architectures. Moreover, computational cost remains small.


title: The Magnetic Topology and Eruption Mechanism of a Multiple-ribbon Flare

abstract: Multiple-ribbon flares are usually complex in their magnetic topologies and eruption mechanisms. In this paper, we investigate an X2.1 flare (SOL2015-03-11T16:22) that occurred in active region 12297 near the center of the solar disk by both potential and nonlinear force-free field models extrapolated with the data observed by the Helioseismic and Magnetic Imager (HMI) on board Solar Dynamics Observatory (SDO). We calculate the three-dimensional squashing degree distribution. The results reveal that there are two flux ropes in this active region, covered by a large scale hyperbolic flux tube (HFT), which is the intersection of quasi-separatrix layers with a null point embedded in it. When the background magnetic field diminishes due to the separation of the northwest dipole and the flux cancellation, the central flux rope rises up forming the two brightest central ribbons. It then squeezes the upper lying HFT structure to generate further brightenings. This very energetic flare with a complex shape is accompanied by a coronal mass ejection (CME). We adopt the simplified line-tied force-balance equation of the current ring model and assign the observed value of the decay index to the equation to simulate the acceleration profile of the CME in the early stage. It is found that the path with an inclination of $45^\circ$ from radial best fits the profile of the actual acceleration.


title: Negative correlations can play a positive role in disordered quantum walks

abstract: We investigate the emerging properties of quantum walks with temporal disorder engineered from a binary Markov chain with tailored correlation, $C$, and disorder strength, $r$. We show that when the disorder is weak -- $r \ll 1$ -- the introduction of negative correlation leads to a counter-intuitive higher production of spin-lattice entanglement entropy, $S_e$, than the setting with positive correlation, that is $S_e(-|C|)>S_e(|C|)$. These results show that negatively correlated disorder plays a more important role in quantum entanglement than it has been assumed in the literature.


title: Tunneling splittings of vibrationally excited states using general instanton paths

abstract: A multidimensional semiclassical method for calculating tunneling splittings in vibrationally excited states of molecules using Cartesian coordinates is developed. It is an extension of the theory by Mil'nikov and Nakamura [$\textit{ J. Chem. Phys.}$ $\textbf{122}$, $124311$ $(2005)$] to asymmetric paths that are necessary for calculating tunneling splitting patterns in multi-well systems, such as water clusters. Additionally, new terms are introduced in the description of the semiclassical wavefunction that drastically improve the splitting estimates for certain systems. The method is based on the instanton theory and builds the semiclassical wavefunction of the vibrationally excited states from the ground-state instanton wavefunction along the minimum action path and its harmonic neighborhood. The splittings of excited states are thus obtained at a negligible added numerical effort. The cost is concentrated, as for the ground-state splittings, in the instanton path optimization and the hessian evaluation along the path. The method can thus be applied without modification to many mid-sized molecules in full dimensionality and in combination with on-the-fly evaluation of electronic potentials. The tests were performed on several model potentials and on the water dimer.


title: Hemodynamic Bigger Hydrostatic Pressure Instead of Lower Shear Stress Aggravates Atherosclerosis

abstract: The traditional view is that low shear stress of blood flow accelerates atherosclerosis. How does a low shear stress damage the vessel or aggravate atherosclerosis? It is a puzzling question. Furthermore, shear stress of veins is similar to that of Relatively Vulnerable Zones (RVZ) of arteries, while intimal thickening of veins is not obvious as the arteries. There are many places in arterioles with lower shear stress as the decrease of blood velocity, while the thickening of arterioles is not as obvious as that of the medium-sized arteries. These issues make it unconvincing for us to explain atherogenesis with lower shear stress. According to Bernoulli fluid equation (P+1/2\r{ho}v^2+\r{ho}gh=constant), in a closed pipe full of flowing blood, at the same height, the hydrostatic pressure is inversely proportional to the second power of blood velocity. At any point of per unit mass of fluid micro cluster here, the reduction of 1/2\r{ho}v^2 would be converted into P. Therefore, hydrostatic pressure(P) is inversely proportional to v^2 in a very short distance or in the same transection of the artery. When blood micro cluster flows over a very short distance or the same transection of the artery, previous studies did not consider the conversion between 1/2\r{ho}v^2 and P. Therefore, low shear stress aggravates atherosclerosis is an appearance, and the essence is that these areas with smaller blood velocity have much bigger hydrostatic pressure, which aggravates atherosclerosis. This perfectly explains the predisposing sites of atherosclerosis.


title: Octupole correlations in light actinides from the interacting boson model based on the Gogny energy density functional

abstract: The quadrupole-octupole coupling and the related spectroscopic properties have been studied for the even-even light actinides $^{218-238}$Ra and $^{220-240}$Th. The Hartree-Fock-Bogoliubov approximation, based on the Gogny-D1M energy density functional, has been employed as a microscopic input, i.e., to obtain (axially symmetric) mean-field potential energy surfaces as functions of the quadrupole and octupole deformation parameters. The mean-field potential energy surfaces have been mapped onto the corresponding bosonic potential energy surfaces using the expectation value of the $sdf$ Interacting Boson Model (IBM) Hamiltonian in the boson condensate state. The strength parameters of the $sdf$-IBM Hamiltonian have been determined via this mapping procedure. The diagonalization of the mapped IBM Hamiltonian provides energies for positive- and negative-parity states as well as wave functions which are employed to obtain transitional strengths. The results of the calculations compare well with available data from Coulomb excitation experiments and point towards a pronounced octupole collectivity around $^{224}$Ra and $^{226}$Th.


title: Deep learning-based transformation of the H&E stain into special stains

abstract: Pathology is practiced by visual inspection of histochemically stained slides. Most commonly, the hematoxylin and eosin (H&E) stain is used in the diagnostic workflow and it is the gold standard for cancer diagnosis. However, in many cases, especially for non-neoplastic diseases, additional "special stains" are used to provide different levels of contrast and color to tissue components and allow pathologists to get a clearer diagnostic picture. In this study, we demonstrate the utility of supervised learning-based computational stain transformation from H&E to different special stains (Masson's Trichrome, periodic acid-Schiff and Jones silver stain) using tissue sections from kidney needle core biopsies. Based on evaluation by three renal pathologists, followed by adjudication by a fourth renal pathologist, we show that the generation of virtual special stains from existing H&E images improves the diagnosis in several non-neoplastic kidney diseases sampled from 58 unique subjects. A second study performed by three pathologists found that the quality of the special stains generated by the stain transformation network was statistically equivalent to those generated through standard histochemical staining. As the transformation of H&E images into special stains can be achieved within 1 min or less per patient core specimen slide, this stain-to-stain transformation framework can improve the quality of the preliminary diagnosis when additional special stains are needed, along with significant savings in time and cost, reducing the burden on healthcare system and patients.


title: The two particle-hole pairs contribution to the dynamic correlation functions of quantum integrable models

abstract: We consider the problem of computing dynamic correlation functions of quantum integrable models employing the thermodynamic form-factor approach. Specifically, we focus on correlations of local operators that conserve the number of particles and consider the 2-particle-hole contribution to their two-point functions. With the method developed being generally applicable to any finite energy and entropy state, our primary focus is on the thermal states. To exemplify this approach, we chose the Lieb-Liniger model and study the leading contribution from 2-particle-hole excitations at small momenta to the dynamic density-density correlation function. We also consider analogous contributions to two-point functions of higher local conserved densities and currents present in integrable theories.


title: Thermal fluctuations in superconducting phases with chiral $ d + id $ and $ s $ symmetry on a triangular lattice

abstract: The behavior of thermal fluctuations of a superconducting order parameter with extended $ s $ and chiral $ d + id $ symmetry is investigated. The study is carried out on a triangular lattice within the framework of the quasi-two-dimensional single-band model with attraction between electrons at neighboring sites. The method of consistent consideration of the order parameter fluctuations and the charge carrier scattering by fluctuations of coupled electron pairs, based on the theory of functional integration is used. The distribution functions of the phase fluctuation probabilities depending on temperature and charge carrier concentration are obtained. The temperature dependences of the amplitudes of the averaged superconducting order parameter are calculated. A phase diagram of superconducting states is constructed for the entire range of variation in the charge carrier concentration $ 0 <n <2 $. Near the boundaries of this range, topologically trivial superconducting states with extended $ s $ symmetry are realized, while a superconducting state with topologically nontrivial chiral $ d + id $ symmetry is realized between them. The calculated anomalous self-energies are compared with the experimental ones obtained using machine learning techniques.


title: Solitons in a box-shaped wavefield with noise: perturbation theory and statistics

abstract: We investigate the fundamental problem of the nonlinear wavefield scattering data corrections in response to a perturbation of initial condition using inverse scattering transform theory. We present a complete theoretical linear perturbation framework to evaluate first-order corrections of the full set of the scattering data within the integrable one-dimensional focusing nonlinear Schr\"odinger (NLSE) equation. The general scattering data portrait reveals nonlinear coherent structures - solitons - playing the key role in the wavefield evolution. Applying the developed theory to a classic box-shaped wavefield we solve the derived equations analytically for a single Fourier mode acting as a perturbation to the initial condition, thus, leading to the sensitivity closed-form expressions for basic soliton characteristics, i.e. the amplitude, velocity, phase and its position. With the appropriate statistical averaging we model the soliton noise-induced effects resulting in compact relations for standard deviations of soliton parameters. Relying on a concept of a virtual soliton eigenvalue we derive the probability of a soliton emergence or the opposite due to noise and illustrate these theoretical predictions with direct numerical simulations of the NLSE evolution. The presented framework can be generalised to other integrable systems and wavefield patterns.


title: Bilirubin lowering effect and safety of a prototype low cost blue light emitting diode (LED) phototherapy device in the treatment of indirect hyperbilirubinemia among healthy term infants in a tertiary government hospital: a pilot study

abstract: Objective: This pilot study was done to evaluate the capability of a prototype low cost blue light emitting diode (LED) phototherapy device in lowering bilirubin levels among healthy term infants diagnosed with indirect hyperbilirubinemia. Methods: Experimental study on term infants diagnosed with indirect hyperbilirubinemia in Ospital ng Makati from May 2016 to November 2016 who underwent phototherapy using the low cost blue LED phototherapy prototype. Results: After 24 hours of phototherapy under the prototype LED phototherapy unit, 16% of the total patients completed treatment as they were already classified in the low risk zone, and another 36% of patients completed treatment after 48 hours. The total bilirubin significantly decreased from baseline bilirubin levels after 24 hours by 16.5% (p = 0.0001). The mean percentage of change of bilirubin reduced after 48 hours of 29.9% was also significant. The proportion of subjects in the high risk zone during baseline to 24th hour went down significantly from 80% to 28% (p = 0.0003), while comparing baseline to 48th hour, the percentage of high risk zone went down from 80% to 9.5% (p = 0.0001). No subjects were reported to have rebound hyperbilirubinemia after discontinuation of phototherapy treatment under the LED prototype. No patient experienced any complication while on phototherapy treatment. Conclusion: The prototype low cost blue light emitting diode (LED) phototherapy was able to lower total serum bilirubin among healthy term infants with indirect hyperbilirubinemia and was safe to use.


title: Mechanics and Thermodynamics of a New Minimal Model of the Atmosphere

abstract: The understanding of the fundamental properties of the climate system has long benefitted from the use of simple numerical models able to parsimoniously represent the essential ingredients of its processes. Here we introduce a new model for the atmosphere that is constructed by supplementing the now-classic Lorenz '96 one-dimensional lattice model with temperature-like variables. The model features an energy cycle that allows for conversion between the kinetic and potential forms and for introducing a notion of efficiency. The model's evolution is controlled by two contributions - a quasi-symplectic and a gradient one, which resemble (yet not conforming to) a metriplectic structure. After investigating the linear stability of the symmetric fixed point, we perform a systematic parametric investigation that allows us to define regions in the parameters space where at steady state stationary, quasi-periodic, and chaotic motions are realised, and study how the terms responsible for defining the energy budget of the system depend on the external forcing injecting energy in the kinetic and in the potential energy reservoirs. Finally, we find preliminary evidence that the model features extensive chaos. We also introduce a more complex version of the model that is able to accommodate for multiscale dynamics and that features an energy cycle that more closely mimics the one of the Earth's atmosphere.

categories: nlin.CD cond-mat.stat-mech physics.ao-ph physics.flu-dyn


title: Multidimensional RFI Framework for Characterising Radio Astronomy Observatories

abstract: Radio Frequency Interference (RFI) has historically plagued radio astronomy, worsening with the rapid spread of electronics and increasing telescope sensitivity. We present a multi-dimensional probabilistic framework for characterising the RFI environment around a radio astronomy site that uses automatically flagged data from the array itself. We illustrate the framework using about 1500 hours of commissioning data from the MeerKAT radio telescope; producing a 6-dimensional array that yields both average RFI occupancy as well as confidence intervals around the mean as a function of key variables (frequency, direction, baseline, time). Our results provide the first detailed view of the MeerKAT RFI environment at high sensitivity as a function of direction, frequency, time of day and baseline. They allow us to track the historical evolution of the RFI and to quantify fluctuations which can be used for alerting on new RFI. As expected we find the major RFI contributors for MeerKAT site are from Global Positioning System (GPS) satellites, flight Distance Measurement Equipment (DME) and the Global System for Mobile (GSM) Communications. Beyond characterising RFI environments our approach allows observers access to the prior probability of RFI in any combination of tracked variables, allowing for more efficient observation planning and data excision.