title: Well-behaved Online Load Balancing Against Strategic Jobs

abstract: In the online load balancing problem on related machines, we have a set of jobs (with different sizes) arriving online, and we need to assign each job to a machine immediately upon its arrival, so as to minimize the makespan, i.e., the maximum completion time. In classic mechanism design problems, we assume that the jobs are controlled by selfish agents, with the sizes being their private information. Each job (agent) aims at minimizing its own cost, which is its completion time plus the payment charged by the mechanism. Truthful mechanisms guaranteeing that every job minimizes its cost by reporting its true size have been well-studied [Aspnes et al. JACM 1997, Feldman et al. EC 2017]. In this paper, we study truthful online load balancing mechanisms that are well-behaved [Epstein et al., MOR 2016]. Well-behavior is important as it guarantees fairness between machines, and implies truthfulness in some cases when machines are controlled by selfish agents. Unfortunately, existing truthful online load balancing mechanisms are not well-behaved. We first show that to guarantee producing a well-behaved schedule, any online algorithm (even non-truthful) has a competitive ratio at least $\Omega(\sqrt{m})$, where m is the number of machines. Then we propose a mechanism that guarantees truthfulness of the online jobs, and produces a schedule that is almost well-behaved. We show that our algorithm has a competitive ratio of $O(\log m)$. Moreover, for the case when the sizes of online jobs are bounded, the competitive ratio of our algorithm improves to $O(1)$. Interestingly, we show several cases for which our mechanism is actually truthful against selfish machines.


title: Bounds on expected propagation time of probabilistic zero forcing

abstract: Probabilistic zero forcing is a coloring game played on a graph where the goal is to color every vertex blue starting with an initial blue vertex set. As long as the graph is connected, if at least one vertex is blue then eventually all of the vertices will be colored blue. The most studied parameter in probabilistic zero forcing is the expected propagation time starting from a given vertex of $G.$ In this paper we improve on upper bounds for the expected propagation time by Geneson and Hogben and Chan et al. in terms of a graph's order and radius. In particular, for a connected graph $G$ of order $n$ and radius $r,$ we prove the bound $\text{ept}(G) = O(r\log(n/r)).$ We also show using Doob's Optional Stopping Theorem and a combinatorial object known as a cornerstone that $\text{ept}(G) \le n/2 + O(\log n).$ Finally, we derive an explicit lower bound $\text{ept}(G)\ge \log_2 \log_2 n.$


title: Properties of the Null Distance and Spacetime Convergence

abstract: The null distance for Lorentzian manifolds was recently introduced by Sormani and Vega. Under mild assumptions on the time function of the spacetime, the null distance gives rise to an intrinsic, conformally invariant metric that induces the manifold topology. We show when warped products of low regularity and globally hyperbolic spacetimes endowed with the null distance are (local) integral current spaces. This metric and integral current structure sets the stage for investigating convergence analogous to Riemannian geometry. Our main theorem is a general convergence result for warped product spacetimes relating uniform, Gromov--Hausdorff and Sormani--Wenger intrinsic flat convergence of the corresponding null distances. In addition, we show that non-uniform convergence of warping functions in general leads to distinct limiting behavior, such as limits that disagree.


title: A self-consistent dynamical system with multiple absolutely continuous invariant measures

abstract: In this paper we study a class of \emph{self-consistent dynamical systems}, self-consistent in the sense that the discrete time dynamics is different in each step depending on current statistics. The general framework admits popular examples such as coupled map systems. Motivated by an example of M. Blank, we concentrate on a special case where the dynamics in each step is a $\beta$-map with some $\beta \geq 2$. Included in the definition of $\beta$ is a parameter $\varepsilon > 0$ controlling the strength of self-consistency. We show such a self-consistent system which has a unique absolutely continuous invariant measure (acim) for $\varepsilon=0$, but at least two for any $\varepsilon > 0$. With a slight modification, we transform this system into one which produces a phase transition-like behavior: it has a unique acim for $0< \varepsilon < \varepsilon^*$, and multiple for sufficiently large values of $\varepsilon$. We discuss the stability of the invariant measures by the help of computer simulations employing the numerical representation of the self-consistent transfer operator.


title: VACL: Variance-Aware Cross-Layer Regularization for Pruning Deep Residual Networks

abstract: Improving weight sparsity is a common strategy for producing light-weight deep neural networks. However, pruning models with residual learning is more challenging. In this paper, we introduce Variance-Aware Cross-Layer (VACL), a novel approach to address this problem. VACL consists of two parts, a Cross-Layer grouping and a Variance Aware regularization. In Cross-Layer grouping the $i^{th}$ filters of layers connected by skip-connections are grouped into one regularization group. Then, the Variance-Aware regularization term takes into account both the first and second-order statistics of the connected layers to constrain the variance within a group. Our approach can effectively improve the structural sparsity of residual models. For CIFAR10, the proposed method reduces a ResNet model by up to 79.5% with no accuracy drop and reduces a ResNeXt model by up to 82% with less than 1% accuracy drop. For ImageNet, it yields a pruned ratio of up to 63.3% with less than 1% top-5 accuracy drop. Our experimental results show that the proposed approach significantly outperforms other state-of-the-art methods in terms of overall model size and accuracy.


title: Data Science in Biomedicine

abstract: We highlight the role of Data Science in Biomedicine. Our manuscript goes from the general to the particular, presenting a global definition of Data Science and showing the trend for this discipline together with the terms of cloud computing and big data. In addition, since Data Science is mostly related to areas like economy or business, we describe its importance in biomedicine. Biomedical Data Science (BDS) presents the challenge of dealing with data coming from a range of biological and medical research, focusing on methodologies to advance the biomedical science discoveries, in an interdisciplinary context.


title: Equivariant Morse theory on Vietoris-Rips complexes & universal spaces for proper actions

abstract: We formalize an equivariant version of Bestvina-Brady discrete Morse theory, and apply it to Vietoris-Rips complexes in order to exhibit finite universal spaces for proper actions for all asymptotically CAT(0) groups.


title: Shadowing the rotating annulus. Part I: Measuring candidate trajectory shadowing times

abstract: An intuitively necessary requirement of models used to provide forecasts of a system's future is the existence of shadowing trajectories that are consistent with past observations of the system: given a system-model pair, do model trajectories exist that stay reasonably close to a sequence of observations of the system? Techniques for finding such trajectories are well-understood in low-dimensional systems, but there is significant interest in their application to high-dimensional weather and climate models. We build on work by Smith et al. [2010, Phys. Lett. A, 374, 2618-2623] and develop a method for measuring the time that individual "candidate" trajectories of high-dimensional models shadow observations, using a model of the thermally-driven rotating annulus in the perfect model scenario. Models of the annulus are intermediate in complexity between low-dimensional systems and global atmospheric models. We demonstrate our method by measuring shadowing times against artificially-generated observations for candidate trajectories beginning a fixed distance from truth in one of the annulus' chaotic flow regimes. The distribution of candidate shadowing times we calculated using our method corresponds closely to (1) the range of times over which the trajectories visually diverge from the observations and (2) the divergence time using a simple metric based on the distance between model trajectory and observations. An empirical relationship between the expected candidate shadowing times and the initial distance from truth confirms that the method behaves reasonably as parameters are varied.

categories: physics.data-an math.DS nlin.CD physics.ao-ph physics.flu-dyn


title: A LOFAR-IRAS cross-match study: the far-infrared radio correlation and the 150-MHz luminosity as a star-formation rate

abstract: Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift z~0:05) and improve the use of the rest-frame 150-MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction. Methods. We cross-match the 60-um selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150-MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers. Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFR_Halpha based on the Halpha line luminosity, SFR_60 based on the rest-frame 60-um luminosity and SFR_IR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (Lsun) = 1.35(0.06) x log SFR_Halpha (Msun/yr) + 3.20(0.06), log L150 (Lsun) = 1.31(0.05) x log SFR_60 (Msun/yr) + 3.14(0.06), and log L150 (Lsun) = 1.37(0.05) x log SFR_IR (Msun/yr) + 3.09(0.05), which show excellent agreement with each other.


title: Gravitational-wave-driven tidal secular instability in neutron star binaries

abstract: We report the existence of a gravitational-wave-driven secular instability in neutron star binaries, acting on the equilibrium tide. The instability is similar to the classic Chandrasekhar-Friedman-Schutz (CFS) instability of normal modes and is active when the spin of the primary star exceeds the orbital frequency of the companion. Modeling the neutron star as a Newtonian n=1 polytrope, we calculate the instability time scale, which can be as low as a few seconds at small orbital separations but still larger than the inspiral time scale. The implications for orbital and spin evolution are also briefly explored, where it is found that the instability slows down the inspiral and decreases the stellar spin.


title: Graph-based data clustering via multiscale community detection

abstract: We present a graph-theoretical approach to data clustering, which combines the creation of a graph from the data with Markov Stability, a multiscale community detection framework. We show how the multiscale capabilities of the method allow the estimation of the number of clusters, as well as alleviating the sensitivity to the parameters in graph construction. We use both synthetic and benchmark real datasets to compare and evaluate several graph construction methods and clustering algorithms, and show that multiscale graph-based clustering achieves improved performance compared to popular clustering methods without the need to set externally the number of clusters.


title: Pluggable Social Artificial Intelligence for Enabling Human-Agent Teaming

abstract: As intelligent systems are increasingly capable of performing their tasks without the need for continuous human input, direction, or supervision, new human-machine interaction concepts are needed. A promising approach to this end is human-agent teaming, which envisions a novel interaction form where humans and machines behave as equal team partners. This paper presents an overview of the current state of the art in human-agent teaming, including the analysis of human-agent teams on five dimensions; a framework describing important teaming functionalities; a technical architecture, called SAIL, supporting social human-agent teaming through the modular implementation of the human-agent teaming functionalities; a technical implementation of the architecture; and a proof-of-concept prototype created with the framework and architecture. We conclude this paper with a reflection on where we stand and a glance into the future showing the way forward.


title: Context-aware Deep Model for Entity Recommendation in Search Engine at Alibaba

abstract: Entity recommendation, providing search users with an improved experience via assisting them in finding related entities for a given query, has become an indispensable feature of today's search engines. Existing studies typically only consider the queries with explicit entities. They usually fail to handle complex queries that without entities, such as "what food is good for cold weather", because their models could not infer the underlying meaning of the input text. In this work, we believe that contexts convey valuable evidence that could facilitate the semantic modeling of queries, and take them into consideration for entity recommendation. In order to better model the semantics of queries and entities, we learn the representation of queries and entities jointly with attentive deep neural networks. We evaluate our approach using large-scale, real-world search logs from a widely used commercial Chinese search engine. Our system has been deployed in ShenMa Search Engine and you can fetch it in UC Browser of Alibaba. Results from online A/B test suggest that the impression efficiency of click-through rate increased by 5.1% and page view increased by 5.5%.


title: Energy Efficient Power Allocation for Device-to-Device Communications Underlaid Cellular Networks Using Stochastic Geometry

abstract: In this paper, we study an energy efficiency maximization problem in uplink for D2D communications underlaid with cellular networks on multiple bands. Utilizing stochastic geometry, we derive closed-form expressions for the average sum rate, successful transmission probability, and energy efficiency of cellular and D2D users. Then, we formulate an optimization problem to jointly maximize the energy efficiency of D2D and cellular users and obtain optimum transmission power of both D2D and cellular users. In the optimization problem, we guarantee the QoS of users by taking into account the success transmission probability on each link. To solve the problem, first we convert the problem into canonical convex form. Afterwards, we solve the problem in two phases, energy efficiency maximization of devices and energy efficiency maximization of cellular users. In the first phase, we maximize the energy efficiency of D2D users and feed the solution to the second phase where we maximize the energy efficiency of cellular users. Simulation results reveal that significant energy efficiency can be attained e.g., 10% energy efficiency improvement compared to fix transmission power in high density scenario.


title: Natural Adversarial Sentence Generation with Gradient-based Perturbation

abstract: This work proposes a novel algorithm to generate natural language adversarial input for text classification models, in order to investigate the robustness of these models. It involves applying gradient-based perturbation on the sentence embeddings that are used as the features for the classifier, and learning a decoder for generation. We employ this method to a sentiment analysis model and verify its effectiveness in inducing incorrect predictions by the model. We also conduct quantitative and qualitative analysis on these examples and demonstrate that our approach can generate more natural adversaries. In addition, it can be used to successfully perform black-box attacks, which involves attacking other existing models whose parameters are not known. On a public sentiment analysis API, the proposed method introduces a 20% relative decrease in average accuracy and 74% relative increase in absolute error.


title: Assessing Fashion Recommendations: A Multifaceted Offline Evaluation Approach

abstract: Fashion is a unique domain for developing recommender systems (RS). Personalization is critical to fashion users. As a result, highly accurate recommendations are not sufficient unless they are also specific to users. Moreover, fashion data is characterized by a large majority of new users, so a recommendation strategy that performs well only for users with prior interaction history is a poor fit to the fashion problem. Critical to addressing these issues in fashion recommendation is an evaluation strategy that: 1) includes multiple metrics that are relevant to fashion, and 2) is performed within segments of users with different interaction histories. Here, we present our multifaceted offline strategy for evaluating fashion RS. Using our proposed evaluation methodology, we compare the performance of three different algorithms, a most popular (MP) items strategy, a collaborative filtering (CF) strategy, and a content-based (CB) strategy. We demonstrate that only by considering the performance of these algorithms across multiple metrics and user segments can we determine the extent to which each algorithm is likely to fulfill fashion users' needs.


title: Equity2Vec: End-to-end Deep Learning Framework for Cross-sectional Asset Pricing

abstract: Pricing assets has attracted significant attention from the financial technology community. We observe that the existing solutions overlook the cross-sectional effects and not fully leveraged the heterogeneous data sets, leading to sub-optimal performance. To this end, we propose an end-to-end deep learning framework to price the assets. Our framework possesses two main properties: 1) We propose Equity2Vec, a graph-based component that effectively captures both long-term and evolving cross-sectional interactions. 2) The framework simultaneously leverages all the available heterogeneous alpha sources including technical indicators, financial news signals, and cross-sectional signals. Experimental results on datasets from the real-world stock market show that our approach outperforms the existing state-of-the-art approaches. Furthermore, market trading simulations demonstrate that our framework monetizes the signals effectively.


title: New X-ray observations toward PSR J1826-1256

abstract: We present the results of XMM-Newton observations toward the pulsar PSR J1826$-$1256, which lies at 5$^{\prime}$.4 from the centroid of the TeV source HESS J1826$-$130. These data show an elongated nebula with the pulsar located in the northern border with the emission in the direction of the peak of the very high energy source. The spectral study of the X-ray emission shows a clear softening of the photon index with increasing distance from the pulsar. We discuss the connection among the pulsar, its nebula, and the HESS source.


title: Countering Language Drift via Visual Grounding

abstract: Emergent multi-agent communication protocols are very different from natural language and not easily interpretable by humans. We find that agents that were initially pretrained to produce natural language can also experience detrimental language drift: when a non-linguistic reward is used in a goal-based task, e.g. some scalar success metric, the communication protocol may easily and radically diverge from natural language. We recast translation as a multi-agent communication game and examine auxiliary training constraints for their effectiveness in mitigating language drift. We show that a combination of syntactic (language model likelihood) and semantic (visual grounding) constraints gives the best communication performance, allowing pre-trained agents to retain English syntax while learning to accurately convey the intended meaning.


title: Radial pulsations of stars on the stage of the final helium flash

abstract: Stellar evolution calculations to the stage of the cooling white dwarf were done for population~I stars with masses on the main sequence $1M_\odot\le M_0\le 1.5M_\odot$. The final helium flash LTP is shown to occur in post--AGB stars with initial masses $1.3M_\odot\le M_0\le 1.32M_\odot$ for the overshooting parameter $f=0.016$. In the case of more effective overshooting ($f=0.018$) the final helium flash occurs at initial masses $1.28M_\odot\le M_0\le 1.3M_\odot$. Fivefold variations of the parameter responsible for the mass loss rate during the post--AGB stage do not affect occurrence of the final helium flash but lead to perceptible changes of the evolutionary time. Selected models of two evolutionary sequences with initial mass $M_0 = 1.3M_\odot$ computed with overshooting parameters $f=0.016$ and $f=0.018$ were used as initial conditions in solution of the equations of hydrodynamics describing radial oscillations of stars on the stage of the final helium flash at effective temperatures less than $10^4$ K. The maximum pulsation period $\Pi=117$ day determined for the evolutionary sequence $M_0=1.3M_\odot$, $f=0.016$ is in a good agreement with observational estimates of the period of FG Sge. The mass, the radius and the effective temperature of the star are $M=0.565M_\odot$, $R=126R_\odot$ and $T_\mathrm{eff}=4445$ K, respectively. At the same time the average period change rate of FG Sge from 1960 to 1990 is nearly three time larger than its theoretical estimate.


title: A Study of Deep Learning for Network Traffic Data Forecasting

abstract: We present a study of deep learning applied to the domain of network traffic data forecasting. This is a very important ingredient for network traffic engineering, e.g., intelligent routing, which can optimize network performance, especially in large networks. In a nutshell, we wish to predict, in advance, the bit rate for a transmission, based on low-dimensional connection metadata ("flows") that is available whenever a communication is initiated. Our study has several genuinely new points: First, it is performed on a large dataset (~50 million flows), which requires a new training scheme that operates on successive blocks of data since the whole dataset is too large for in-memory processing. Additionally, we are the first to propose and perform a more fine-grained prediction that distinguishes between low, medium and high bit rates instead of just "mice" and "elephant" flows. Lastly, we apply state-of-the-art visualization and clustering techniques to flow data and show that visualizations are insightful despite the heterogeneous and non-metric nature of the data. We developed a processing pipeline to handle the highly non-trivial acquisition process and allow for proper data preprocessing to be able to apply DNNs to network traffic data. We conduct DNN hyper-parameter optimization as well as feature selection experiments, which clearly show that fine-grained network traffic forecasting is feasible, and that domain-dependent data enrichment and augmentation strategies can improve results. An outlook about the fundamental challenges presented by network traffic analysis (high data throughput, unbalanced and dynamic classes, changing statistics, outlier detection) concludes the article.


title: Relativistic Lattice Boltzmann Methods: Theory and Applications

abstract: We present a systematic account of recent developments of the relativistic Lattice Boltzmann method (RLBM) for dissipative hydrodynamics. We describe in full detail a unified, compact and dimension-independent procedure to design relativistic LB schemes capable of bridging the gap between the ultra-relativistic regime, $k_{\rm B} T \gg mc^2$, and the non-relativistic one, $k_{\rm B} T \ll mc^2$. We further develop a systematic derivation of the transport coefficients as a function of the kinetic relaxation time in $d=1,2,3$ spatial dimensions. The latter step allows to establish a quantitative bridge between the parameters of the kinetic model and the macroscopic transport coefficients. This leads to accurate calibrations of simulation parameters and is also relevant at the theoretical level, as it provides neat numerical evidence of the correctness of the Chapman-Enskog procedure. We present an extended set of validation tests, in which simulation results based on the RLBMs are compared with existing analytic or semi-analytic results in the mildly-relativistic ($k_{\rm B} T \sim mc^2$) regime for the case of shock propagations in quark-gluon plasmas and laminar electronic flows in ultra-clean graphene samples. It is hoped and expected that the material collected in this paper may allow the interested readers to reproduce the present results and generate new applications of the RLBM scheme.

categories: hep-lat gr-qc hep-ph hep-th physics.comp-ph physics.flu-dyn


title: ArduCode: Predictive Framework for Automation Engineering

abstract: Automation engineering is the task of integrating, via software, various sensors, actuators, and controls for automating a real-world process. Today, automation engineering is supported by a suite of software tools including integrated development environments (IDE), hardware configurators, compilers, and runtimes. These tools focus on the automation code itself, but leave the automation engineer unassisted in their decision making. This can lead to increased time for software development because of imperfections in decision making leading to multiple iterations between software and hardware. To address this, this paper defines multiple challenges often faced in automation engineering and propose solutions using machine learning to assist engineers tackle such challenges. We show that machine learning can be leveraged to assist the automation engineer in classifying automation, finding similar code snippets, and reasoning about the hardware selection of sensors and actuators. We validate our architecture on two real datasets consisting of 2,927 Arduino projects, and 683 Programmable Logic Controller (PLC) projects. Our results show that paragraph embedding techniques can be utilized to classify automation using code snippets with precision close to human annotation, giving an F1-score of 72%. Further, we show that such embedding techniques can help us find similar code snippets with high accuracy. Finally, we use autoencoder models for hardware recommendation and achieve a p@3 of 0.79 and p@5 of 0.95.


title: PySPH: a Python-based framework for smoothed particle hydrodynamics

abstract: PySPH is an open-source, Python-based, framework for particle methods in general and Smoothed Particle Hydrodynamics (SPH) in particular. PySPH allows a user to define a complete SPH simulation using pure Python. High-performance code is generated from this high-level Python code and executed on either multiple cores, or on GPUs, seamlessly. It also supports distributed execution using MPI. PySPH supports a wide variety of SPH schemes and formulations. These include, incompressible and compressible fluid flow, elastic dynamics, rigid body dynamics, shallow water equations, and other problems. PySPH supports a variety of boundary conditions including mirror, periodic, solid wall, and inlet/outlet boundary conditions. The package is written to facilitate reuse and reproducibility. This paper discusses the overall design of PySPH and demonstrates many of its features. Several example results are shown to demonstrate the range of features that PySPH provides.


title: A probabilistic proof of the spherical excess formula

abstract: This note offers a probabilistic proof of Girard's angle excess formula for the area of a spherical triangle, based on the observation that an unbounded 3-dimensional convex cone, with single vertex at the origin, has only three kinds of 2-dimensional orthogonal projections: a 2-dimensional convex cone with one vertex, a 2-dimensional half-plane (this is an outcome of probability zero), and a 2-dimensional plane.


title: Anisotropic Optical Shock Waves in Isotropic Media with Giant Nonlocal Nonlinearity

abstract: Dispersive shock waves in thermal optical media belong to the third-order nonlinear phenomena, whose intrinsic irreversibility is described by time asymmetric quantum mechanics. Recent studies demonstrated that nonlocal wave breaking evolves in an exponentially decaying dynamics ruled by the reversed harmonic oscillator, namely, the simplest irreversible quantum system in the rigged Hilbert spaces. The generalization of this theory to more complex scenarios is still an open question. In this work, we use a thermal third-order medium with an unprecedented giant Kerr coefficient, the M-Cresol/Nylon mixed solution, to access an extremely-nonlinear highly-nonlocal regime and realize anisotropic shock waves. We prove that a superposition of the Gamow vectors in an ad hoc rigged Hilbert space describes the nonlinear beam propagation beyond the shock point. Specifically, the resulting rigged Hilbert space is a tensorial product between the reversed and the standard harmonic oscillators spaces. The anisotropy turns out from the interaction of trapping and antitrapping potentials in perpendicular directions. Our work opens the way to a complete description of novel intriguing shock phenomena, and those mediated by extreme nonlinearities.


title: Studying the Solar system dynamics using pulsar timing arrays and the LINIMOSS dynamical model

abstract: Pulsar timing arrays (PTAs) can be used to study the Solar-system ephemeris (SSE), the errors of which can lead to correlated timing residuals and significantly contribute to the PTA noise budget. Most Solar-system studies with PTAs assume the dominance of the term from the shift of the Solar-system barycentre (SSB). However, it is unclear to which extent this approximation can be valid, since the perturbations on the planetary orbits may become important as data precision keeps increasing. To better understand the effects of SSE uncertainties on pulsar timing, we develop the LINIMOSS dynamical model of the Solar system, based on the SSE of Guangyu Li. Using the same input parameters as DE435, the calculated planetary positions by LINIMOSS are compatible with DE435 at centimetre level over a 20-year timespan, which is sufficiently precise for pulsar-timing applications. We utilize LINIMOSS to investigate the effects of SSE errors on pulsar timing in a fully dynamical way, by perturbing one SSE parameter per trial and examining the induced timing residuals. For the outer planets, the timing residuals are dominated by the SSB shift, as assumed in previous work. For the inner planets, the variations in the orbit of the Earth are more prominent, making previously adopted assumptions insufficient. The power spectra of the timing residuals have complex structures, which may introduce false signals in the search of gravitational waves. We also study how to infer the SSE parameters using PTAs, and calculate the accuracy of parameter estimation.


title: A remark on embedding of a cylinder on a real commutative Banach algebra

abstract: Let $A$ be a real commutative Banach algebra with unity. Let $a_0\in A\setminus\{0\}$. Let $\mathbb Z a_0:=\{na_0\}_{n\in \mathbb Z}$. Then, $\mathbb Z a_0$ is a discrete subgroup of $A$. For any $n\in \mathbb Z$, the Frechet derivative of the mapping $$x \, \in \, A \ \ \ \mapsto \ \ \ x+na_0 \, \in \, A$$ is the identity map on $A$ and, especially, an $A$-linear transformation on $A$. So, the quotient group $A/(\mathbb Z a_0)$ is a $1$-dimensional $A$-manifold and the covering projection $$x \, \in \, A \ \ \ \mapsto \ \ \ x+\mathbb Z a_0 \, \in \, A/(\mathbb Z a_0)$$ is an $A$-map. We call $A/(\mathbb Z a_0)$ the $1$-dimensional $A$-cylinder by $a_0$. Let $T$ be a compact Hausdorff space. Suppose that there exist $t_1\in T$ and $t_2\in T$ such that $t_1\not=t_2$ holds. Then, the set $C(T;\mathbb R)$ of all real-valued continuous functions on $T$ is a real commutative Banach algebra with unity and $\mathbb R \, \subsetneq \, C(T;\mathbb R)$ holds. In this paper, we show that there exists $a_0 \, \in \, C(T;\mathbb R)\setminus \mathbb R$ such that for any $k\, \in \, \mathbb N$, the $1$-dimensional $C(T;\mathbb R)$-cylinder $(C(T;\mathbb R))/(\mathbb Z a_0)$ by $a_0$ cannot be embedded in the finite direct product space $(C(T;\mathbb R))^k$ as a $C(T;\mathbb R)$-submanifold.


title: Unrolling Ternary Neural Networks

abstract: The computational complexity of neural networks for large scale or real-time applications necessitates hardware acceleration. Most approaches assume that the network architecture and parameters are unknown at design time, permitting usage in a large number of applications. This paper demonstrates, for the case where the neural network architecture and ternary weight values are known a priori, that extremely high throughput implementations of neural network inference can be made by customising the datapath and routing to remove unnecessary computations and data movement. This approach is ideally suited to FPGA implementations as a specialized implementation of a trained network improves efficiency while still retaining generality with the reconfigurability of an FPGA. A VGG style network with ternary weights and fixed point activations is implemented for the CIFAR10 dataset on Amazon's AWS F1 instance. This paper demonstrates how to remove 90% of the operations in convolutional layers by exploiting sparsity and compile-time optimizations. The implementation in hardware achieves 90.9 +/- 0.1% accuracy and 122 k frames per second, with a latency of only 29 us, which is the fastest CNN inference implementation reported so far on an FPGA.


title: Emergence of global synchronization in directed excitatory networks of type I neurons

abstract: The collective behaviour of neural networks depends on the cellular and synaptic properties of the neurons. The phase-response curve (PRC) is an experimentally obtainable measure of cellular properties that quantifies the shift in the next spike time of a neuron as a function of the phase at which stimulus is delivered to that neuron. The neuronal PRCs can be classified as having either purely positive values (type I) or distinct positive and negative regions (type II). Networks of type 1 PRCs tend not to synchronize via mutual excitatory synaptic connections. We study the synchronization properties of identical type I and type II neurons, assuming unidirectional synapses. Performing the linear stability analysis and the numerical simulation of the extended Kuramoto model, we show that feedforward loop motifs favour synchronization of type I excitatory and inhibitory neurons, while feedback loop motifs destroy their synchronization tendency. Moreover, large directed networks, either without feedback motifs or with many of them, have been constructed from the same undirected backbones, and a high synchronization level is observed for directed acyclic graphs with type I neurons. It has been shown that, the synchronizability of type I neurons depends on both the directionality of the network connectivity and the topology of its undirected backbone. The abundance of feedforward motifs enhances the synchronizability of the directed acyclic graphs.


title: Untangling Optical Emissions of the Jet and Accretion Disk in the Flat-Spectrum Radio Quasar 3C 273 with Reverberation Mapping Data

abstract: 3C 273 is an intensively monitored flat-spectrum radio quasar with both a beamed jet and blue bump together with broad emission lines. The coexistence of the comparably prominent jet and accretion disk leads to complicated variability properties. Recent reverberation mapping monitoring for 3C 273 revealed that the optical continuum shows a distinct long-term trend that does not have a corresponding echo in the Hbeta fluxes. We compile multi-wavelength monitoring data from the Swift archive and other ground-based programs and clearly find two components of emissions at optical wavelength. One component stems from the accretion disk itself and the other component can be ascribed to the jet contribution, which also naturally accounts for the non-echoed trend in reverberation mapping data. We develop an approach to decouple the optical emissions from the jet and accretion disk in 3C 273 with the aid of multi-wavelength monitoring data. By assuming the disk emission has a negligible polarization in consideration of the low inclination of the jet, the results show that the jet contributes a fraction of ~10% at the minimum and up to ~40% at the maximum to the total optical emissions. This is the first time to provide a physical interpretation to the "detrending" manipulation conventionally adopted in reverberation mapping analysis. Our work also illustrates the importance of appropriately analyzing variability properties in cases of coexisting jets and accretion disks.


title: Incorporating long-range physics in atomic-scale machine learning

abstract: The most successful and popular machine learning models of atomic-scale properties derive their transferability from a locality ansatz. The properties of a large molecule or a bulk material are written as a sum over contributions that depend on the configurations within finite atom-centered environments. The obvious downside of this approach is that it cannot capture non-local, non-additive effects such as those arising due to long-range electrostatics or quantum interference. We propose a solution to this problem by introducing non-local representations of the system that are remapped as feature vectors that are defined locally and are equivariant in O(3). We consider in particular one form that has the same asymptotic behavior as the electrostatic potential. We demonstrate that this framework can capture non-local, long-range physics by building a model for the electrostatic energy of randomly distributed point-charges, for the unrelaxed binding curves of charged organic molecular dimers, and for the electronic dielectric response of liquid water. By combining a representation of the system that is sensitive to long-range correlations with the transferability of an atom-centered additive model, this method outperforms current state-of-the-art machine-learning schemes, and provides a conceptual framework to incorporate non-local physics into atomistic machine learning.


title: Group-Testing on Hypergraphs with Variable-Cost Tests: A Power Systems Case Study

abstract: Assessing risk of cascading failure in an electrical grid requires identifying many small "defective" subsets of the N elements in a power system, such that the simultaneous failure of all elements in a defective set triggers a large cascading blackout. Most supersets of defective sets will also initiate cascades. This property was leveraged by Random Chemistry (RC), a state-of-the-art algorithm for efficiently finding minimal defective sets. While not previously acknowledged as such, RC is a group-testing algorithm, in that it tests pools of candidate solutions. RC then reduces the size of defective pools until a minimal defective set is found. Using a synthetic model of the Western United States power grid, we show that run times are minimized with a much smaller initial pool size than previously recommended, nearly doubling the speed of RC. We compare RC to a proposed new alternative group-testing algorithm (SIGHT). Both algorithms have their advantages; while SIGHT is shown to require fewer total tests than RC, RC requires fewer tests of defective sets. Thus, the relative speed of each algorithm depends on the relative cost of defective vs. non defective tests, a factor not previously considered in the group-testing literature. In the power systems application, tests of defective sets are much more computationally costly, so RC is faster. However, with only a single control parameter (compared to many in RC), SIGHT is easier to optimize, and in applications where the cost ratio of defective:non-defective tests is low enough, SIGHT will be faster than RC.


title: Real Numbers are the Hidden Variables of Classical Mechanics

abstract: Do scientific theories limit human knowledge? In other words, are there physical variables hidden by essence forever? We argue for negative answers and illustrate our point on chaotic classical dynamical systems. We emphasize parallels with quantum theory and conclude that the common real numbers are, de facto, the hidden variables of classical physics. Consequently, real numbers should not be considered as "physically real" and classical mechanics, like quantum physics, is indeterministic.


title: Elements of resolvent methods in fluid mechanics: notes for an introductory short course

abstract: This is a collection of notes for part of a short course on modal methods in fluid mechanics held at DAMTP, University of Cambridge, in the summer of 2019. These notes introduce the reader to resolvent analysis as it is currently used in fluid mechanics. They are aimed at the average beginning PhD student, and will serve as a base from which to explore the literature.


title: Ring Rydberg Composites

abstract: The properties and behaviour of a Ring Rydberg Composite are explicated. This system consists of a ring of ground state atoms centered on a Rydberg atom, whose electron elastically scatters off the ground state atoms. We transform the electronic Hamiltonian into a discrete tight-binding representation in which the on-site energies and long-range hopping between sites are controlled and mediated by the Rydberg electron. From this new representation, which to a large extent enables an analytic treatment, we derive scaling laws and analytic expressions for the wave functions and eigenspectrum. The interface between ring and Rydberg geometries leads to a range of rich properties which can be tuned as a function of ring size, number of scatterers, and principal quantum number.


title: Transport of probe particles in polymer network: effects of probe size, network rigidity and probe-polymer interaction

abstract: Fundamental understanding of the effect of microscopic parameters on the dynamics of probe particles in different complex environments has wide implications. Examples include diffusion of proteins in the biological hydrogels, porous media, polymer matrix, etc. Here, we use extensive molecular dynamics simulations to investigate the dynamics of the probe particle in a polymer network on a diamond lattice which provides substantial crowding to mimic the cellular environment. Our simulations show that the dynamics of the probe increasingly becomes restricted, non-Gaussian and subdiffusive on increasing the network rigidity, binding affinity and the probe size. In addition, the velocity autocorrelation functions show negative dips owing to the viscoelasticity and caging due to surrounding network. These observations go with general experimental findings. Surprisingly for a probe particle of size comparable to the mesh size, unrestricted motion engulfing large length scales has been observed. This happens with the more flexible polymer network, which is easily pushed by the bigger probe. Our study gives a general qualitative picture of transport of probes in gel like medium, as encountered in different contexts.


title: Virtual organelle self-coding for fluorescence imaging via adversarial learning

abstract: Fluorescence microscopy plays a vital role in understanding the subcellular structures of living cells. However, it requires considerable effort in sample preparation related to chemical fixation, staining, cost, and time. To reduce those factors, we present a virtual fluorescence staining method based on deep neural networks (VirFluoNet) to transform fluorescence images of molecular labels into other molecular fluorescence labels in the same field-of-view. To achieve this goal, we develop and train a conditional generative adversarial network (cGAN) to perform digital fluorescence imaging demonstrated on human osteosarcoma U2OS cell fluorescence images captured under Cell Painting staining protocol. A detailed comparative analysis is also conducted on the performance of the cGAN network between predicting fluorescence channels based on phase contrast or based on another fluorescence channel using human breast cancer MDA-MB-231 cell line as a test case. In addition, we implement a deep learning model to perform autofocusing on another human U2OS fluorescence dataset as a preprocessing step to defocus an out-focus channel in U2OS dataset. A quantitative index of image prediction error is introduced based on signal pixel-wise spatial and intensity differences with ground truth to evaluate the performance of prediction to high-complex and throughput fluorescence. This index provides a rational way to perform image segmentation on error signals and to understand the likelihood of mis-interpreting biology from the predicted image. In total, these findings contribute to the utility of deep learning image regression for fluorescence microscopy datasets of biological cells, balanced against savings of cost, time, and experimental effort. Furthermore, the approach introduced here holds promise for modeling the internal relationships between organelles and biomolecules within living cells.


title: Theoretical and numerical considerations on the polar (exponential) form of octonions and elements of higher-order Cayley-Dickson algebras

abstract: The article is devoted to the issue of the polar form of octonions. This is a~continuation of the works initiated by Hahn and Snopek in their articles from 2011. The results presented in the article show errors made in previous considerations and suggest the possibility of their improvement. The presented numerical method gives promising results, which as a result of further work can give analytical formulas for angles in the polar representation.


title: Computed stereo lensless X-ray imaging

abstract: The ability to gain insights into the 3D properties of artificial or biological systems is often critical. However, 3D structures are difficult to retrieve at low dose and with extremely fast processing, as most techniques are based on acquiring and computing hundreds of 2D angular projections. This is even more challenging with ultrashort X-rays which allow realizing nanometre scale studies and ultrafast time resolved 2D movies. Here we show that computer stereo vision concepts can be transposed to X-rays. We demonstrate nanoscale three-dimensional reconstruction from a single ultrafast acquisition. Two diffraction patterns are recorded simultaneously on a single CCD camera and after phase retrieval two stereo images are reconstructed. A 3D representation of the sample is then computed from quantitative disparity maps with about 130x130x380nm3 voxel resolution in a snapshot of 20 femtoseconds. We extend our demonstration to phase contrast X-ray stereo imaging and reveal hidden 3D features of a sample. Computed phase stereo imaging will find scientific applications at X-ray free electron lasers, synchrotrons and laser-based sources, but also in fast industrial and medical 3D diagnostics.


title: Is culture related to strong science? An empirical investigation

abstract: National culture is among those societal factors which could influence research and innovation activities. In this study, we investigated the associations of two national culture models with citation impact of nations (measured by the proportion of papers belonging to the 10% and 1% most cited papers in the corresponding fields, PPtop 10% and PPtop 1%). Bivariate statistical analyses showed that of six Hofstede's national culture dimensions (HNCD), uncertainty avoidance and power distance had a statistically significant negative associa-tion, while individualism and indulgence had a statistically significant positive associationwith both citation impact indicators. The study also revealed that of two Inglehart-Welzel cultural values (IWCV), the value survival versus self-expression is statistically significantly related to citation impact indicators. We additionally calculated multiple regression analyses controlling for the possible effects of confounding factors including national self-citations, international co-authorships, invest-ments in research and development, international migrant stock, number of researchers ofeach nation, language, and productivity. The results revealed that the statistically significant associations of HNCD with citation impact indicators disappeared. But the statistically significant relationship between survivals versus self-expression values and citation impact indicators remained stable even after controlling for the confounding variables. Thus, the freedom of expression and trust in society might contribute to better scholarly communication systems, higher level of international collaborations, and further quality research.


title: The Shear Viscosity to Entropy Density Ratio of Hagedorn States

abstract: The fireball concept of Rolf Hagedorn, developed in the 1960's, is an alternative description of hadronic matter. Using a recently derived mass spectrum, we use the transport model GiBUU to calculate the shear viscosity of a gas of such Hagedorn states, applying the Green-Kubo method to Monte-Carlo calculations. Since the entropy density is rising ad infinitum near $T_H$, this leads to a very low shear viscosity to entropy density ratio near $T_H$. Further, by comparing our results with analytic expressions, we find a nice extrapolation behavior, indicating that a gas of Hagedorn states comes close or even below the boundary $1/4\pi$ from AdS-CFT.


title: Li Evolution and the Open Cluster NGC 6819: A Correlation Between Li Depletion and Spindown in Dwarfs More Massive Than the F-Dwarf Li-Dip

abstract: Spectroscopy of 333 NGC 6819 stars and Gaia astrometry are used to map Li evolution from the giant branch tip to 0.5 mag below the Li dip. Isochrone comparison with [Fe/H] = -0.04, based upon neural network spectroscopic analysis, produces an age of 2.25 (2.4) Gyr for E(B-V) = 0.16 (0.14) and (m-M) = 12.40 (12.29). Despite originating outside the Li dip, only 10% of single subgiants/giants have measurable Li. Above the Li dip, the limiting A(Li) for single stars is 3.2 +/- 0.1 but the lower range is comparable to that found within the dip. The F-dwarf Li dip profile agrees with the Hyades/Praesepe, evolved forward. The Li level among stars populating the plateau fainter than the Li dip is A(Li) = 2.83 +/- 0.16; the dispersion is larger than expected from spectroscopic error alone. Comparison of Li and V_ROT distributions among turnoff stars in NGC 7789, NGC 2506, NGC 3680, and NGC 6819 indicates that rotational spindown from the main sequence is critical in defining the boundaries of the Li dip. For higher mass dwarfs, spindown is likewise correlated with Li depletion, creating a second dip, but at higher mass and on a longer time scale. The Li distribution among evolved stars of NGC 6819 is more representative of the older M67, where subgiant and giant stars emerge from within the Li dip, than the younger NGC 7789, where a broad range in V_ROT among the turnoff stars likely produces a range in mass among the giants.


title: Kondo effect in a Aharonov-Casher interferometer

abstract: We consider a model describing a spin field-effect transistor based on a quantum nanowire with a tunable spin-orbit interaction embedded between two ferromagnetic leads with anticollinear magnetization. We investigate a regime of a strong interplay between resonance Kondo scattering and interference associated with the Aharonov-Casher effect. Using the Keldysh technique at weak coupling regime we calculate perturbatively the charge current. It is predicted that the effects of the spin-orbit interaction result in a non-vanishing current for any spin polarization of the leads including the case of fully polarized anti-collinear contacts. We analyze the influence of the Aharonov-Casher phase and degree of spin polarization in the leads onto a Kondo temperature.


title: Skin cancer detection based on deep learning and entropy to detect outlier samples

abstract: We describe our methods that achieved the 3rd and 4th places in tasks 1 and 2, respectively, at ISIC challenge 2019. The goal of this challenge is to provide the diagnostic for skin cancer using images and meta-data. There are nine classes in the dataset, nonetheless, one of them is an outlier and is not present on it. To tackle the challenge, we apply an ensemble of classifiers, which has 13 convolutional neural networks (CNN), we develop two approaches to handle the outlier class and we propose a straightforward method to use the meta-data along with the images. Throughout this report, we detail each methodology and parameters to make it easy to replicate our work. The results obtained are in accordance with the previous challenges and the approaches to detect the outlier class and to address the meta-data seem to be work properly.


title: Random Walk Equivalence to the Compressible Baker Map and the Kaplan-Yorke Approximation to Its Information Dimension

abstract: Simple time-reversible systems can generate {\it irreversible} flows satisfying the Second Law of Thermodynamics. Maps, and equivalent random walks, can also do this. We study a pair of time-reversible Baker Maps, $N2$ and $N3$, which generate dissipative{\it fractal} phase-space structures. Steadily decreasing phase-space volumes correspond to the dissipation associated with entropy production. Like three smooth reversible dissipative one-body phase-space flows developed in the 1980s and 1990s our maps generate fractal distributions, but in two dimensions rather than three, simplifying visualization and analyses. The continuity equation, which quantifies phase-volume loss, motivates study of the fractals' reduced ``information dimensions'', which were approximated by Kaplan and Yorke in terms of two-dimensional maps' two Lyapunov exponents. The maps studied here generate fractal (fractional dimensional) distributions in their phase spaces. By mapping uniformly dense grids of points, fractal dimensions can be determined by ``area-wise'' mappings. Beginning with a uniform grid area-wise mapping of the $N2$ Baker Map provides an information dimension of 1.78969. Alternatively, as many as a trillion iterations, starting from an arbitrary point, gives a smaller ``point-wise'' dimensionality, $1.741_5$. Neither of these precisely determined estimates matches the Kaplan-Yorke conjecture value, 1.7337. In the course of studying these three different approaches to information dimension we developed random walk equivalents to both mappings, which greatly simplifies analyses. We found that for the older $ N2$ Baker map the three approaches all disagree with one another! We later discovered that for the newer $N3$ Baker mapping the three approaches to information dimension, area-wise, point-wise and Kaplan-Yorke, agree.


title: On the theoretical prospects of multiport devices for photon-number-resolving detections

abstract: Ideal photon-number-resolving detectors form a class of important optical components in quantum optics and quantum information theory. In this article, we theoretically investigate the potential of multiport devices having reconstruction performances approaching that of the Fock-state measurement. By recognizing that all multiport devices are minimally complete, we first provide a general analytical framework to describe the tomographic accuracy (or quality) of these devices. Next, we show that a perfect multiport device with an infinite number of output ports functions as either the Fock-state measurement when photon losses are absent or binomial mixtures of Fock-state measurements when photon losses are present, and derive their respective expressions for the tomographic transfer function. This function is the scaled asymptotic mean squared-error of the reconstructed photon-number distributions uniformly averaged over all distributions in the probability simplex. We then supply more general analytical formulas for the transfer function for finite numbers of output ports in both the absence and presence of photon losses. The effects of photon losses on the photon-number resolving power of both infinite- and finite-size multiport devices are also investigated.


title: Real-space imaging of confined magnetic skyrmion tubes

abstract: Magnetic skyrmions are topologically nontrivial particles with a potential application as information elements in future spintronic device architectures. While they are commonly portrayed as two dimensional objects, in reality magnetic skyrmions are thought to exist as elongated, tube-like objects extending through the thickness of the sample. Study of this skyrmion tube (SkT) state is highly relevant for investigating skyrmion metastability and for implementation in recently proposed magnonic computing. However, direct experimental imaging of skyrmion tubes has yet to be reported. Here, we demonstrate the first real-space observation of skyrmion tubes in a lamella of FeGe using resonant magnetic x-ray imaging and comparative micromagnetic simulations, confirming their extended structure. The formation of these structures at the edge of the sample highlights the importance of confinement and edge effects in the stabilisation of the SkT state, opening the door to further investigations into this unexplored dimension of the skyrmion spin texture.


title: Local Large Deviation Principle, Large Deviation Principle and Information theory for the Signal -to- Interference -Plus- Noise Ratio Graph Models

abstract: Given devices space $D$, an intensity measure $\lambda m\in(0,\infty)$, a transition kernel $Q$ from the space $D$ to positive real numbers $(0,\infty,$ a path-loss function (which depends on the Euclidean distance between the devices and a positive constant $\alpha$), we define a Marked Poisson Point process (MPPP). For a given MPPP and technical constants $\tau_{\lambda},\gamma_{\lambda}:(0,\,\infty)\to (0,\infty),$ we define a Marked Signal-to- Interference and Noise Ratio (SINR) graph, and associate with it two empirical measures; the \emph{empirical marked measure} and the \emph{empirical connectivity measure}. For a class of marked SINR graphs, we prove a joint \emph{ large deviation principle}(LDP) for these empirical measures, with speed $\lambda$ in the $\tau$-topology. From the joint large deviation principle for the empirical marked measure and the empirical connectivity measure, we obtain an Asymptotic Equipartition Property(AEP) for network structured data modelled as a marked SINR graph. Specifically, we show that for large dense marked SINR graph one require approximately about $\lambda^{2}H(Q\times Q)/\log 2$ bits to transmit the information contained in the network with high probability, where $ H(Q\times Q)$ is a properly defined entropy for the exponential transition kernel with parameter $c$. Further, we prove a \emph {local large deviation principle} (LLDP) for the class of marked SINR graphs on $D,$ where $\lambda[\tau_{\lambda}(a)\gamma_{\lambda}(a)+\lambda\tau_{\lambda}(b)\gamma_{\lambda}(b)]\to \beta(a,b),$ $ a,b\in (0,\infty)$, with speed $\lambda$ from a \emph{ spectral potential} point. From the LLDP we derive a conditional LDP for the marked SINR graphs.


title: Non-linear diffusion of cosmic rays escaping from supernova remnants: Cold partially neutral atomic and molecular phases

abstract: We aim to elucidate cosmic ray (CR) propagation in the weakly ionized environments of supernova remnants (SNRs) basing our analysis on the cosmic ray cloud (CRC) model developed by \citet{2013ApJ...768...73M} and \citet{2016MNRAS.461.3552N}. We solved two transport equations simultaneously: one for the CR pressure and one for the Alfv\'en wave energy density where CRs are initially confined in the SNR shock. Cosmic rays trigger a streaming instability and produce slab-type resonant Alfv\'en modes. The self-generated turbulence is damped by ion-neutral collisions and by noncorrelated interaction with Alfv\'en modes generated at large scales. We show that CRs leaking in cold dense phases such as those found in cold neutral medium (CNM) and diffuse molecular medium (DiM) can still be confined over distances of a few tens of parsecs from the CRC center for a few thousand years. At 10 TeV, CR diffusion can be suppressed by two or three orders of magnitude. This effect results from a reduced ion-neutral collision damping in the decoupled regime. We calculate the grammage of CRs in these environments. We find that in both single and multi-phase setups at 10 GeV, CNM and DiM media can produce grammage in the range 10-20 $\rm{g/cm^2}$ in the CNM and DiM phases. At 10 TeV, because of nonlinear propagation the grammage increases to values in the range 0.5-20 $\rm{g/cm^2}$ in these two phases. We also present preliminary calculations in inhomogeneous interstellar medium combining two or three different phases where we obtain the same trends.


title: Optimizing magnetoresistive sensor signal-to-noise via pinning field tuning

abstract: The presence of magnetic noise in magnetoresistive-based magnetic sensors degrades their detection limit at low frequencies. In this paper, different ways of stabilizing the magnetic sensing layer to suppress magnetic noise are investigated by applying a pinning field, either by an external field, internally in the stack or by shape anisotropy. We show that these three methods are equivalent, could be combined and that there is a competition between noise suppression and sensitivity reduction, which results in an optimum total pinning field, for which the detection limit of the sensor is improved up to a factor of ten.


title: Byzantine-Resilient Stochastic Gradient Descent for Distributed Learning: A Lipschitz-Inspired Coordinate-wise Median Approach

abstract: In this work, we consider the resilience of distributed algorithms based on stochastic gradient descent (SGD) in distributed learning with potentially Byzantine attackers, who could send arbitrary information to the parameter server to disrupt the training process. Toward this end, we propose a new Lipschitz-inspired coordinate-wise median approach (LICM-SGD) to mitigate Byzantine attacks. We show that our LICM-SGD algorithm can resist up to half of the workers being Byzantine attackers, while still converging almost surely to a stationary region in non-convex settings. Also, our LICM-SGD method does not require any information about the number of attackers and the Lipschitz constant, which makes it attractive for practical implementations. Moreover, our LICM-SGD method enjoys the optimal $O(md)$ computational time-complexity in the sense that the time-complexity is the same as that of the standard SGD under no attacks. We conduct extensive experiments to show that our LICM-SGD algorithm consistently outperforms existing methods in training multi-class logistic regression and convolutional neural networks with MNIST and CIFAR-10 datasets. In our experiments, LICM-SGD also achieves a much faster running time thanks to its low computational time-complexity.


title: On list 3-dynamic coloring of near-triangulations

abstract: An $r$-dynamic $k$-coloring of a graph $G$ is a proper $k$-coloring such that for any vertex $v$, there are at least $\min\{r, deg_G(v) \}$ distinct colors in $N_G(v)$. The $r$-dynamic chromatic number $\chi_r^d(G)$ of a graph $G$ is the least $k$ such that there exists an $r$-dynamic $k$-coloring of $G$. The list $r$-dynamic chromatic number of a graph $G$ is denoted by $ch_r^d(G)$. Loeb et al. $[11]$ showed that $ch_3^d(G)\leq 10$ for every planar graph $G$, and there is a planar graph $G$ with $\chi_3^d(G)= 7$. In this paper, we study a special class of planar graphs which have better upper bounds of $ch_3^d(G)$. We prove that $ch_3^d(G) \leq 6$ if $G$ is a planar graph which is near-triangulation, where a near-triangulation is a planar graph whose bounded faces are all 3-cycles.


title: The International Pulsar Timing Array: Second data release

abstract: In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.


title: The galaxy stellar mass function from CCSNe with improved photo-z techniques

abstract: In Sedgwick et al. (2019) we introduced and utilised a method to combat surface brightness and mass biases in galaxy sample selection, using core-collapse supernovae (CCSNe) as pointers towards their host galaxies, in order to: (i) search for low-surface brightness galaxies (LSBGs); (ii) assess the contributions of galaxies at a given mass to the star-formation-rate density (SFRD); and (iii) infer from this, using estimates of specific star-formation (SF) rate, the form of the SF-galaxy stellar mass function (GSMF). A CCSN-selection of SF-galaxies allows a probe of the form of the SFRD and GSMF deep into the dwarf galaxy mass regime. In the present work, we give improved constraints on our estimates of the SFRD and star-forming GSMF, in light of improved photometric redshift estimates required for estimates of galaxy stellar mass. The results are consistent with a power-law increase to SF-galaxy number density down to our low stellar mass limit of $\sim 10^{6.2}$ M$_{\odot}$. No deviation from the high-mass version of the surface brightness - mass relation is found in the dwarf mass regime. These findings imply no truncation to galaxy formation processes at least down to $\sim 10^{6.2}$ M$_{\odot}$.


title: ODEA: Orbital Dynamics in a complex Evolving Architecture -- Application to the planetary system HD 106906

abstract: Mixed-variable symplectic integrators are widely used in orbital dynamics. However, they have been developed for Solar system-type architectures, and can not handle evolving hierarchy, in particular in systems with two or more stellar components. Such configuration may have occurred in the history of HD 106906, a tight pair of F-type stars surrounded by a debris disk and a planetary-mass companion on a wide orbit. We present the new algorithm ODEA, based on the symplectic algorithm SWIFT HJS, that can model any system (binary,...) with unstable architecture. We study the peculiar system HD 106906 as a testcase for the code. We define and compute a criterion based on acceleration ratios to indicate when the initial hierarchy is not relevant anymore. A new hierarchy is then computed. The code is applied to study the two fly-bys that occurred on system HD 106906, recently evidenced by De Rosa & Kalas (2019), to determine if they could account for the wide orbit of the planet. Thousands of simulations have been performed to account for the uncertainty on the perturbers coordinates and velocities. The algorithm is able to handle any change of hierarchy, temporary or not. We used it to fully model HD 106906 encounters. The simulations confirm that the fly-bys could have stabilized the planet orbit, and show that it can account for the planet probable misalignment with respect to the disk plane as well as the disk morphology. However, that requires a small distance at closest approach (< 0.05 pc), and this configuration is not guaranteed. ODEA is a very good choice for the study of non-Solar type architecture. It can now adapt to an evolving hierarchy, and is thus suitable to study capture of planets and dust. Further observations of the perturbers, in particular their radial velocity, are required to conclude on the effects of the fly-by on system HD 106906.


title: Prediction of the strength and timing of sunspot cycle 25 reveal decadal-scale space environmental conditions

abstract: The Sun's activity cycle governs the radiation, particle and magnetic flux in the heliosphere creating hazardous space weather. Decadal-scale variations define space climate and force the Earth's atmosphere. However, predicting the solar cycle is challenging. Current understanding indicates a short window for prediction best achieved at previous cycle minima. Utilizing magnetic field evolution models for the Sun's surface and interior we perform the first century-scale, data-driven simulations of solar activity and present a scheme for extending the prediction window to a decade. Our ensemble forecast indicates cycle 25 would be similar or slightly stronger than the current cycle and peak around 2024. Sunspot cycle 25 may thus reverse the substantial weakening trend in solar activity which has led to speculation of an imminent Maunder-like grand minimum and cooling global climate. Our simulations demonstrate fluctuation in the tilt angle distribution of sunspots is the dominant mechanism responsible for solar cycle variability.


title: DeepPrivacy: A Generative Adversarial Network for Face Anonymization

abstract: We propose a novel architecture which is able to automatically anonymize faces in images while retaining the original data distribution. We ensure total anonymization of all faces in an image by generating images exclusively on privacy-safe information. Our model is based on a conditional generative adversarial network, generating images considering the original pose and image background. The conditional information enables us to generate highly realistic faces with a seamless transition between the generated face and the existing background. Furthermore, we introduce a diverse dataset of human faces, including unconventional poses, occluded faces, and a vast variability in backgrounds. Finally, we present experimental results reflecting the capability of our model to anonymize images while preserving the data distribution, making the data suitable for further training of deep learning models. As far as we know, no other solution has been proposed that guarantees the anonymization of faces while generating realistic images.


title: Efficient Interleaved Batch Matrix Solvers for CUDA

abstract: In this paper we present a new methodology for data accesses when solving batches of Tridiagonal and Pentadiagonal matrices that all share the same LHS matrix. By only storing one copy of this matrix there is a significant reduction in storage overheads and the authors show that there is also a performance increase in terms of compute time. These two results combined lead to an overall more efficient implementation over the current state of the art algorithms cuThomasBatch and cuPentBatch, allowing for a greater number of systems to be solved on a single GPU.


title: Stirring up an embedded star cluster with a moving gas filament

abstract: We perform simulations to test the effects of a moving gas filament on a young star cluster (i.e. the "Slingshot" Model). We model Orion Nebula Cluster-like clusters as Plummer spheres and the Integral Shaped Filament gas as a cylindrical potential. We observe that in a static filament, an initially spherical cluster evolves naturally into an elongated distribution of stars. For sinusoidal moving filaments, we observe different remnants, and classify them into 4 categories.%: 3 different objects and one transition object. "Healthy" clusters, where almost all the stars stay inside the filament and the cluster; "destroyed" clusters are the opposite case, with almost no particles in the filament or near the centre of density of the clusters; "ejected" clusters, where a large fraction of stars are close to the centre of density of the stars , but almost none of them in the filament; and "transition" clusters, where roughly the same number of particles is ejected from the cluster and from the filament. An {{Orion Nebula Cluster-like}} cluster might stay inside the filament or be ejected, but it will not be destroyed.


title: The Keck Planet Imager and Characterizer: Demonstrating advanced exoplanet characterization techniques for future extremely large telescopes

abstract: The Keck Planet Imager and Characterizer (KPIC) is an upgrade to the Keck II adaptive optics system enabling high contrast imaging and high-resolution spectroscopic characterization of giant exoplanets in the mid-infrared (2-5 microns). The KPIC instrument will be developed in phases. Phase I entails the installation of an infrared pyramid wavefront sensor (PyWFS) based on a fast, low-noise SAPHIRA IR-APD array. The ultra-sensitive infrared PyWFS will enable high contrast studies of infant exoplanets around cool, red, and/or obscured targets in star forming regions. In addition, the light downstream of the PyWFS will be coupled into an array of single-mode fibers with the aid of an active fiber injection unit (FIU). In turn, these fibers route light to Keck's high-resolution infrared spectrograph NIRSPEC, so that high dispersion coronagraphy (HDC) can be implemented for the first time. HDC optimally pairs high contrast imaging and high-resolution spectroscopy allowing detailed characterization of exoplanet atmospheres, including molecular composition, spin measurements, and Doppler imaging. Here we provide an overview of the instrument, its science scope, and report on recent results from on-sky commissioning of Phase I. The instrument design and techniques developed will be key for more advanced instrument concepts needed for the extremely large telescopes of the future.


title: Integrating cross-modality hallucinated MRI with CT to aid mediastinal lung tumor segmentation

abstract: Lung tumors, especially those located close to or surrounded by soft tissues like the mediastinum, are difficult to segment due to the low soft tissue contrast on computed tomography images. Magnetic resonance images contain superior soft-tissue contrast information that can be leveraged if both modalities were available for training. Therefore, we developed a cross-modality educed learning approach where MR information that is educed from CT is used to hallucinate MRI and improve CT segmentation. Our approach, called cross-modality educed deep learning segmentation (CMEDL) combines CT and pseudo MR produced from CT by aligning their features to obtain segmentation on CT. Features computed in the last two layers of parallelly trained CT and MR segmentation networks are aligned. We implemented this approach on U-net and dense fully convolutional networks (dense-FCN). Our networks were trained on unrelated cohorts from open-source the Cancer Imaging Archive CT images (N=377), an internal archive T2-weighted MR (N=81), and evaluated using separate validation (N=304) and testing (N=333) CT-delineated tumors. Our approach using both networks were significantly more accurate (U-net $P <0.001$; denseFCN $P <0.001$) than CT-only networks and achieved an accuracy (Dice similarity coefficient) of 0.71$\pm$0.15 (U-net), 0.74$\pm$0.12 (denseFCN) on validation and 0.72$\pm$0.14 (U-net), 0.73$\pm$0.12 (denseFCN) on the testing sets. Our novel approach demonstrated that educing cross-modality information through learned priors enhances CT segmentation performance


title: Moduli Problem of Hitchin Pairs over Deligne-Mumford Stack

abstract: We define the moduli problem of Hitchin pairs over Deligne-Mumford Stack and prove this moduli problem is represented by a separated and locally finitely presented algebraic space, which is considered as the moduli space of Hitchin pairs over Deligne-Mumford stack.


title: A Loosely Self-stabilizing Protocol for Randomized Congestion Control with Logarithmic Memory

abstract: We consider congestion control in peer-to-peer distributed systems. The problem can be reduced to the following scenario: Consider a set $V$ of $n$ peers (called clients in this paper) that want to send messages to a fixed common peer (called server in this paper). We assume that each client $v \in V$ sends a message with probability $p(v) \in [0,1)$ and the server has a capacity of $\sigma \in \mathbb{N}$, i.e., it can recieve at most $\sigma$ messages per round and excess messages are dropped. The server can modify these probabilities when clients send messages. Ideally, we wish to converge to a state with $\sum p(v) = \sigma$ and $p(v) = p(w)$ for all $v,w \in V$. We propose a loosely self-stabilizing protocol with a slightly relaxed legimate state. Our protocol lets the system converge from any initial state to a state where $\sum p(v) \in \left[\sigma \pm \epsilon\right]$ and $|p(v)-p(w)| \in O(\frac{1}{n})$. This property is then maintained for $\Omega(n^{\mathfrak{c}})$ rounds in expectation. In particular, the initial client probabilities and server variables are not necessarily well-defined, i.e., they may have arbitrary values. Our protocol uses only $O(W + \log n)$ bits of memory where $W$ is length of node identifers, making it very lightweight. Finally we state a lower bound on the convergence time an see that our protocol performs asymptotically optimal (up to some polylogarithmic factor).


title: SN2018kzr: a rapidly declining transient from the destruction of a white dwarf

abstract: We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks.


title: Crossover from fast to slow dynamics in quantum Ising chains with long range interactions

abstract: Quantum many body systems with long range interactions are known to display many fascinating phenomena experimentally observable in trapped ions, Rydberg atoms and polar molecules. Among these are dynamical phase transitions which occur after an abrupt quench in spin chains with interactions decaying as $r^{-\alpha}$ and whose critical dynamics depend crucially on the power $\alpha$: for systems with $\alpha<1$ the transition is sharp while for $\alpha>1$ it fans out in a chaotic crossover region. In this paper we explore the fate of critical dynamics in Ising chains with long-range interactions when the transverse field is ramped up with a finite speed. While for abrupt quenches we observe a chaotic region that widens as $\alpha$ is increased, the width of the crossover region diminishes as the time of the ramp increases, suggesting that chaos will disappear altogether and be replaced by a sharp transition in the adiabatic limit.


title: What do Deep Networks Like to Read?

abstract: Recent research towards understanding neural networks probes models in a top-down manner, but is only able to identify model tendencies that are known a priori. We propose Susceptibility Identification through Fine-Tuning (SIFT), a novel abstractive method that uncovers a model's preferences without imposing any prior. By fine-tuning an autoencoder with the gradients from a fixed classifier, we are able to extract propensities that characterize different kinds of classifiers in a bottom-up manner. We further leverage the SIFT architecture to rephrase sentences in order to predict the opposing class of the ground truth label, uncovering potential artifacts encoded in the fixed classification model. We evaluate our method on three diverse tasks with four different models. We contrast the propensities of the models as well as reproduce artifacts reported in the literature.


title: PREMA: A Predictive Multi-task Scheduling Algorithm For Preemptible Neural Processing Units

abstract: To amortize cost, cloud vendors providing DNN acceleration as a service to end-users employ consolidation and virtualization to share the underlying resources among multiple DNN service requests. This paper makes a case for a "preemptible" neural processing unit (NPU) and a "predictive" multi-task scheduler to meet the latency demands of high-priority inference while maintaining high throughput. We evaluate both the mechanisms that enable NPUs to be preemptible and the policies that utilize them to meet scheduling objectives. We show that preemptive NPU multi-tasking can achieve an average 7.8x, 1.4x, and 4.8x improvement in latency, throughput, and SLA satisfaction, respectively.


title: Probing Flavor Non-Universal Theories Through Higgs Physics at the LHC and Future Colliders

abstract: We explored the possibility that Higgs coupling to new physics violates flavor universality. In particular, we parameterize such models with dimension-six effective operators which modify the coupling between the first generation quarks, Higgs boson, and $Z$ boson. Through the use of boosted Higgsstrahlung events at both the HL-LHC and potential future hadron colliders, as well as existing ATLAS data for background estimates, we projected constraints on the scale of new physics as a function of the Wilson coefficient. The high energy $Zh$ process will provide unique information about these class of operators, and the sensitivity is competitive with the LEP electroweak precision measurements. We include different scenarios of the overall systematic uncertainties and the PDF uncertainties when presenting the projected sensitivities. We also discuss the constraints from FCNCs to these flavor-violating models and the complementarity of the exotic Higgs decay to the $Zh$ process.


title: MBWU: Benefit Quantification for Data Access Function Offloading

abstract: The storage industry is considering new kinds of storage devices that support data access function offloading, i.e. the ability to perform data access functions on the storage device itself as opposed to performing it on a separate compute system to which the storage device is connected. But what is the benefit of offloading to a storage device that is controlled by an embedded platform, very different from a host platform? To quantify the benefit, we need a measurement methodology that enables apple-to-apple comparisons between different platforms. We propose a Media-based Work Unit (MBWU, pronounced "MibeeWu"), and an MBWU-based measurement methodology to standardize the platform efficiency evaluation so as to quantify the benefit of offloading. To demonstrate the merit of this methodology, we implemented a prototype to automate quantifying the benefit of offloading the key-value data access function.


title: TMA: Tera-MACs/W Neural Hardware Inference Accelerator with a Multiplier-less Massive Parallel Processor

abstract: Computationally intensive Inference tasks of Deep neural networks have enforced revolution of new accelerator architecture to reduce power consumption as well as latency. The key figure of merit in hardware inference accelerators is the number of multiply-and-accumulation operations per watt (MACs/W), where, the state-of-the-arts MACs/W remains several hundreds Giga-MACs/W. We propose a Tera-MACS/W neural hardware inference Accelerator (TMA) with 8-bit activations and scalable integer weights less than 1-byte. The architectures main feature is configurable neural processing element for matrix-vector operations. The proposed neural processing element has Multiplier-less Massive Parallel Processor to work without any multiplications, which makes it attractive for energy efficient high-performance neural network applications. We benchmark our systems latency, power, and performance using Alexnet trained on ImageNet. Finally, we compared our accelerators throughput and power consumption to the prior works. The proposed accelerator outperforms the state of the art in terms of energy and area achieving 2.3 TMACS/W@1.0 V, 65 nm CMOS technology.


title: Strong converse inequalities for the weighted multivariate Bernstein-Durrmeyer operator on the simplex via multipliers

abstract: It is demonstrated that multiplier methods naturally yield better constants in strong converse inequalities for the Bernstein-Durrmeyer operator. The absolute constants obtained in some of the inequalities are independent of the weight and the dimension. The estimates are stated in terms of the $K$-functional that is naturally associated to the operator.


title: Maximum likelihood degree of the two-dimensional linear Gaussian covariance model

abstract: In algebraic statistics, the maximum likelihood degree of a statistical model is the number of complex critical points of its log-likelihood function. A priori knowledge of this number is useful for applying techniques of numerical algebraic geometry to the maximum likelihood estimation problem. We compute the maximum likelihood degree of a generic two-dimensional subspace of the space of $n\times n$ Gaussian covariance matrices. We use the intersection theory of plane curves to show that this number is $2n-3$.


title: NoCs in Heterogeneous 3D SoCs: Co-Design of Routing Strategies and Microarchitectures

abstract: Heterogeneous 3D System-on-Chips (3D SoCs) are the most promising design paradigm to combine sensing and computing within a single chip. A special characteristic of communication networks in heterogeneous 3D SoCs is the varying latency and throughput in each layer. As shown in this work, this variance drastically degrades the network performance. We contribute a co-design of routing algorithms and router microarchitecture that allows to overcome these performance limitations. We analyze the challenges of heterogeneity: Technology-aware models are proposed for communication and thereby identify layers in which packets are transmitted slower. The communication models are precise for latency and throughput under zero load. The technology model has an area error and a timing error of less than 7.4% for various commercial technologies from 90 to 28nm. Second, we demonstrate how to overcome limitations of heterogeneity by proposing two novel routing algorithms called Z+(XY)Z- and ZXYZ that enhance latency by up to 6.5x compared to conventional dimension order routing. Furthermore, we propose a high vertical-throughput router microarchitecture that is adjusted to the routing algorithms and that fully overcomes the limitations of slower layers. We achieve an increased throughput of 2 to 4x compared to a conventional router. Thereby, the dynamic power of routers is reduced by up to 41.1% and we achieve improved flit latency of up to 2.26x at small total router area costs between 2.1% and 10.4% for realistic technologies and application scenarios.


title: Anisotropic magnetic excitations of a frustrated bilinear-biquadratic spin model -- Implications for spin waves of detwinned iron pnictides

abstract: Elucidating the nature of spin excitations is important to understanding the mechanism of superconductivity in the iron pnictides. Motivated by recent inelastic neutron scattering measurements in the nearly 100% detwinned BaFe$_{2}$As$_{2}$, we study the spin dynamics of an $S=1$ frustrated bilinear-biquadratic Heisenberg model in the antiferromagnetic phase with wavevector $(\pi,0)$. The biquadratic interactions are treated in a dynamical way using a flavor-wave theory in an $SU(3)$ representation. Besides the magnon (dipolar) excitations, the biquadratic interactions give rise to quadrupolar excitations at high energies. We find that the quadrupolar wave significantly influences, in an energy dependent way, the anisotropy between the spin excitation spectra along the $(\pi,0)$ and $(0,\pi)$ directions in the wave vector space. Our theoretical results capture the essential behavior of the spin dynamics measured in the antiferromagnetic phase of the detwinned BaFe$_2$As$_2$. More generally, our results underscore the importance of electron correlation effects for the microscopic physics of the iron pnictides.


title: Human Languages in Source Code: Auto-Translation for Localized Instruction

abstract: Computer science education has promised open access around the world, but access is largely determined by what human language you speak. As younger students learn computer science it is less appropriate to assume that they should learn English beforehand. To that end we present CodeInternational, the first tool to translate code between human languages. To develop a theory of non-English code, and inform our translation decisions, we conduct a study of public code repositories on GitHub. The study is to the best of our knowledge the first on human-language in code and covers 2.9 million Java repositories. To demonstrate CodeInternational's educational utility, we build an interactive version of the popular English-language Karel reader and translate it into 100 spoken languages. Our translations have already been used in classrooms around the world, and represent a first step in an important open CS-education problem.


title: Quasi-steady uptake and bacterial community assembly in a mathematical model of soil-phosphorus mobility

abstract: We mathematically model the uptake of phosphorus by a soil community consisting of a plant and two bacterial groups: copiotrophs and oligotrophs. Four equilibrium states emerge, one for each of the species monopolising the resource and dominating the community and one with coexistence of all species. We show that the dynamics are controlled by the ratio of chemical adsorption to bacterial death permitting either oscillatory states or quasi-steady uptake. We show how a steady state can emerge which has soil and plant nutrient content unresponsive to increased fertilization. However, the additional fertilization supports the copiotrophs leading to community reassembly. Our results demonstrate the importance of time-series measurements in nutrient uptake experiments.


title: Star catalog position and proper motion corrections in asteroid astrometry II: The Gaia era

abstract: Astrometric positions of moving objects in the Solar System have been measured using a variety of star catalogs in the past. Previous work has shown that systematic errors in star catalogs can affect the accuracy of astrometric observations. That, in turn, can influence the resulting orbit fits for minor planets. In order to quantify these systematic errors, we compare the positions and proper motion of stellar sources in the most utilized star catalogs to the second release of the Gaia star catalog. The accuracy of Gaia astrometry allows us to unambiguously identify local biases and derive a scheme that can be used to correct past astrometric observations of solar system objects. Here we provide a substantially improved debiasing scheme for 26 astrometric catalogs that were extensively used in minor planet astrometry. Revised corrections near the galactic center eliminate artifacts that could be traced back to reference catalogs used in previous debiasing schemes. Median differences in stellar positions between catalogs now tend to be on the order of several tens of milliarcseconds (mas) but can be as large as 175 mas. Median stellar proper motion corrections scatter around 0.3 mas/yr and range from 1 to 4 mas/yr for star catalogs with and without proper motion, respectively. The tables in this work are meant to be applied to existing optical observations. They are not intended to correct new astrometric measurments as those should make use of the Gaia astrometric catalog. Since previous debiasing schemes already reduced systematics in past observations to a large extent, corrections beyond the current work may not be needed in the foreseeable future.


title: Learning Hierarchically Structured Concepts

abstract: We study the question of how concepts that have structure get represented in the brain. Specifically, we introduce a model for hierarchically structured concepts and we show how a biologically plausible neural network can recognize these concepts, and how it can learn them in the first place. Our main goal is to introduce a general framework for these tasks and prove formally how both (recognition and learning) can be achieved. We show that both tasks can be accomplished even in presence of noise. For learning, we analyze Oja's rule formally, a well-known biologically-plausible rule for adjusting the weights of synapses. We complement the learning results with lower bounds asserting that, in order to recognize concepts of a certain hierarchical depth, neural networks must have a corresponding number of layers.


title: Digital Beings as an option to study gut flora evolution and adaptation

abstract: In this work, we introduce a computational model for the study of the host-bacteria interaction and the influence of the intestinal microbiota on the behavior and feeding pattern of an individual. The model is based on digital entities, which we've called Digital Beings (DBs), modeled using dynamic systems and genetic algorithms. We have successfully tested the use of the DBs by reproducing observation in previously made studies using rats and humans. Among these studies, we highlight those on how the bacteria in an individual's stomach could influence their eating behavior and how a controlled and continuous diet can affect the longevity of a certain population. Our results point that the Digital Beings can be used as a tool for supporting the devising of experiments and corroborating with theoretical hypotheses, reducing the number of in vivo tests.


title: Tight-binding bond parameters for dimers across the periodic table from density-functional theory

abstract: We obtain parameters for non-orthogonal and orthogonal TB models from two-atomic molecules for all combinations of elements of period 1 to 6 and group 3 to 18 of the periodic table. The TB bond parameters for 1711 homoatomic and heteroatomic dimers show clear chemical trends. In particular, using our parameters we compare to the rectangular d-band model, the reduced sp TB model as well as canonical TB models for sp- and d-valent systems which have long been used to gain qualitative insight into the interatomic bond. The transferability of our dimer-based TB bond parameters to bulk systems is discussed exemplarily for the bulk ground-state structures of Mo and Si. Our dimer-based TB bond parameters provide a well-defined and promising starting point for developing refined TB parameterizations and for making the insight of TB available for guiding materials design across the periodic table.


title: In-gap states and strain-tuned band convergence in layered structure trivalent iridate K0.75Na0.25IrO2

abstract: Iridium oxides (iridates) provide a good platform to study the delicate interplay between spin-orbit coupling (SOC) interactions, electron correlation effects, Hund's coupling and lattice degree of freedom. However, overwhelming investigations primarily focus on tetravalent (Ir4+, 5d5) and pentavalent (Ir5+, 5d4) iridates, far less attention has been paid to iridates with other valence states. Here, we pay our attention to a less-explored trivalent (Ir3+, 5d6) iridates, K0.75Na0.25IrO2, crystalizing in a triangular lattice with edge-sharing IrO6 octahedra and alkali metal ions intercalated [IrO2]- layers. We theoretically determine the preferred occupied positions of the alkali metal ions from energetic viewpoints and reproduce the experimentally observed semiconducting behavior and nonmagnetic (NM) properties. The SOC interactions play a critical role in the band dispersion, resulting in NM Jeff = 0 states. More intriguingly, our electronic structure not only uncovers the presence of in-gap states and explains the abnormal low activation energy in K0.75Na0.25IrO2, but also predicts the band edge can be effectively modulated by mechanical strain. Especially, the in-gap states feature with enhanced band-convergence characteristics by 6% compressive strain, which will greatly enhance the electrical conductivity of K0.75Na0.25IrO2. Present work sheds new lights on the unconventional electronic structures of the trivalent iridates, indicating its promising application as nanoelectronic and thermoelectric material.


title: Flooding and Diameter in General Weighted Random Graphs

abstract: We study in this paper, the first passage percolation on a random graph model, the configuration model. We first introduce, the notions of weighted diameter, which is the maximum of the weighted lengths of all optimal paths between any two vertices in the graph, and the flooding time, which represents the time (weighted length) needed to reach all the vertices in the graph starting from a uniformly chosen vertex. Our result consists of describing the asymptotic behavior of the diameter and the flooding time, as the number of vertices n tends to infinity, in the case where the weight distribution G has an exponential tail behavior, and proving that this category of distributions is the largest possible for which the asymptotic behavior holds.