title: Analytic I-Love-C relations for realistic neutron stars

abstract: Recent observations of neutron stars (NS) with radio, X-rays and gravitational waves have begun to constrain the equation of state (EoS) for nuclear matter beyond the nuclear saturation density. To one's surprise, there exist approximate universal relations connecting certain bulk properties of NSs that are insensitive to the underlying EoS and having important applications on probing fundamental physics including nuclear and gravitational physics. To date, analytic works on universal relations for realistic NSs are lacking, which may lead to a better understanding of the origin of the universality. Here, we focus on the universal relations between the compactness(C), moment of inertia (I), and tidal deformability(related to the Love number), and derive analytic, approximate I-Love-C relations. To achieve this, we construct slowly-rotating/ tidally-deformed NS solutions analytically starting from an extended Tolman VII model that accurately describes non-rotating realistic NSs, allowing us to extract the moment of inertia and the tidal deformability on top of the compactness. We solve the field equations analytically by expanding them about the Newtonian limit and keeping up to 6th order in the stellar compactness. Based on these analytic solutions, we can mathematically demonstrate the O(10%) EoS variation in the I-C and Love-C relations and the O(1%) variation in the I-Love relation that have previously been found numerically. Our new analytic relations agree more accurately with numerical results for realistic NSs(especially the I-C and Love-C ones) than the analytic relations for constant density stars derived in previous work. Based on these analytic findings, we attribute a possible origin of the universality for the I-C and Love-C relations to the fact that the energy density of realistic NSs can be approximated as a quadratic function, as is the case for the Tolman VII solution.


title: New incompressible symmetric tensor categories in positive characteristic

abstract: We propose a method of constructing abelian envelopes of symmetric rigid monoidal Karoubian categories over an algebraically closed field $\bf k$. If ${\rm char}({\bf k})=p>0$, we use this method to construct generalizations ${\rm Ver}_{p^n}$, ${\rm Ver}_{p^n}^+$ of the incompressible abelian symmetric tensor categories defined in arXiv:1807.05549 for $p=2$ and by Gelfand-Kazhdan and Georgiev-Mathieu for $n=1$. Namely, ${\rm Ver}_{p^n}$ is the abelian envelope of the quotient of the category of tilting modules for $SL_2(\bf k)$ by the $n$-th Steinberg module, and ${\rm Ver}_{p^n}^+$ is its subcategory generated by $PGL_2(\bf k)$-modules. We show that ${\rm Ver}_{p^n}$ are reductions to characteristic $p$ of Verlinde braided tensor categories in characteristic zero, which explains the notation. We study the structure of these categories in detail, and in particular show that they categorify the real cyclotomic rings $\mathbb{Z}[2\cos(2\pi/p^n)]$, and that ${\rm Ver}_{p^n}$ embeds into ${\rm Ver}_{p^{n+1}}$. We conjecture that every symmetric tensor category of moderate growth over $\bf k$ admits a fiber functor to the union ${\rm Ver}_{p^\infty}$ of the nested sequence ${\rm Ver}_{p}\subset {\rm Ver}_{p^2}\subset\cdots$. This would provide an analog of Deligne's theorem in characteristic zero and a generalization of the result of arXiv:1503.01492, which shows that this conjecture holds for fusion categories, and then moreover the fiber functor lands in ${\rm Ver}_p$.


title: Systematic Analysis of Distributed Optimization Algorithms over Jointly-Connected Networks

abstract: We consider the distributed optimization problem, where a group of agents work together to optimize a common objective by communicating with neighboring agents and performing local computations. For a given algorithm, we use tools from robust control to systematically analyze the performance in the case where the communication network is time-varying. In particular, we assume only that the network is jointly connected over a finite time horizon (commonly referred to as B-connectivity), which does not require connectivity at each time instant. When applied to the distributed algorithm DIGing, our bounds are orders of magnitude tighter than those available in the literature.


title: Holography of geodesic flows, harmonizing metrics, and billiards' dynamics

abstract: Let $(M, g)$ be a Riemannian manifold with boundary, where $g$ is a non-trapping metric. Let $SM$ be the space of the spherical tangent to $M$ bundle, and $v^g$ the geodesic vector field on $SM$. We study the scattering maps $C_{v^g}: \partial^+_1SM \to \partial^-_1SM$, generated by the $v^g$-flow, and the dynamics of the billiard maps $B_{v^g, \tau}: \partial^+_1SM \to \partial^+_1SM$, where $\tau$ denotes an involution, mimicking the elastic reflection from the the boundary $\partial M$. We getting a variety of holography theorems that tackle the inverse scattering problems for $C_{v^g}$ and theorems that describe the dynamics of $B_{v^g, \tau}$. Our main tools are a Lyapunov function $F: SM \to \mathbb R$ for $v^g$ and a special harmonizing Riemannian metrics $g^\bullet$ on $SM$, a metric in which $dF$ is harmonic. For such metrics $g^\bullet$, we get a family of isoperimetric inequalities of the type $vol_{g^\bullet}(SM) \leq vol_{g^\bullet |}(\partial(SM))$ and formulas for the average volume of the minimal hypesufaces $\{F^{-1}(c)\}_{c \in F(SM)}$. We investigate the interplay between the harmonizing metrics $g^\bullet$ and the classical Sasaki metric $gg$ on $SM$. Assuming ergodicity of $B_{v^g, \tau}$, we also get Santal\'{o}-Chernov type formulas for the average length of free geodesic segments in $M$ and for the average variation of the Lyapunov function $F$ along the $v^g$-trajectories.


title: A closed-form solution for optimal mean-reverting trading strategies

abstract: When prices reflect all available information, they oscillate around an equilibrium level. This oscillation is the result of the temporary market impact caused by waves of buyers and sellers. This price behavior can be approximated through an Ornstein-Uhlenbeck (O-U) process. Market makers provide liquidity in an attempt to monetize this oscillation. They enter a long position when a security is priced below its estimated equilibrium level, and they enter a short position when a security is priced above its estimated equilibrium level. They hold that position until one of three outcomes occur: (1) they achieve the targeted profit; (2) they experience a maximum tolerated loss; (3) the position is held beyond a maximum tolerated horizon. All market makers are confronted with the problem of defining profit-taking and stop-out levels. More generally, all execution traders acting on behalf of a client must determine at what levels an order must be fulfilled. Those optimal levels can be determined by maximizing the trader's Sharpe ratio in the context of O-U processes via Monte Carlo experiments. This paper develops an analytical framework and derives those optimal levels by using the method of heat potentials.


title: A mildly relativistic outflow from the energetic, fast-rising blue optical transient CSS161010 in a dwarf galaxy

abstract: We present X-ray and radio observations of the Fast Blue Optical Transient (FBOT) CRTS-CSS161010 J045834-081803 (CSS161010 hereafter) at t=69-531 days. CSS161010 shows luminous X-ray ($L_x\sim5\times 10^{39}\,\rm{erg\,s^{-1}}$) and radio ($L_{\nu}\sim10^{29}\,\rm{erg\,s^{-1}Hz^{-1}}$) emission. The radio emission peaked at ~100 days post transient explosion and rapidly decayed. We interpret these observations in the context of synchrotron emission from an expanding blastwave. CSS161010 launched a mildly relativistic outflow with velocity $\Gamma\beta c\ge0.55c$ at ~100 days. This is faster than the non-relativistic AT2018cow ($\Gamma\beta c\sim0.1c$) and closer to ZTF18abvkwla ($\Gamma\beta c\ge0.3c$ at 63 days). The inferred initial kinetic energy of CSS161010 ($E_k\gtrsim10^{51}$ erg) is comparable to that of long Gamma Ray Bursts (GRBs), but the ejecta mass that is coupled to the mildly relativistic outflow is significantly larger ($\sim0.01-0.1\,\rm{M_{\odot}}$). This is consistent with the lack of observed gamma-rays. The luminous X-rays were produced by a different emission component to the synchrotron radio emission. CSS161010 is located at ~150 Mpc in a dwarf galaxy with stellar mass $M_{*}\sim10^{7}\,\rm{M_{\odot}}$ and specific star formation rate $sSFR\sim 0.3\,\rm{Gyr^{-1}}$. This mass is among the lowest inferred for host-galaxies of explosive transients from massive stars. Our observations of CSS161010 are consistent with an engine-driven aspherical explosion from a rare evolutionary path of a H-rich stellar progenitor, but we cannot rule out a stellar tidal disruption event on a centrally-located intermediate mass black hole. Regardless of the physical mechanism, CSS161010 establishes the existence of a new class of rare (rate $<0.4\%$ of the local core-collapse supernova rate) H-rich transients that can launch mildly relativistic outflows.


title: Challenges of Bridging the Gap between Mass People and Welfare Organizations in Bangladesh

abstract: Computing for the development of marginalized communities is a big deal of challenges for researchers. Different social organizations are working to develop the conditions of a specialized marginalized community namely Street Children, one of the most underprivileged communities in Bangladesh. However, lack of proper engagement among different social welfare organizations, donors, and the mass community limits the goal of the development of street children. Developing a virtual organization hub can eliminate communication gap as well as the information gap by involving people of all communities. However, some human imposed stigmas may often limit the rate of success of potential virtual computing solutions intended for organizations working with the marginalized communities, which we also face in our case. After a partial successful deployment, the design itself needs to be self comprehensive and trustworthy in order to overcome the stigmas that demand a reasonable amount of time. Moreover, after a wide scalable deployment, it is yet to be investigated whether the design of our computational solution can attain the goal for the facilitation of the organizations so that those organizations can become more effective for the development of street children than before.


title: Simulating Solar Flare Irradiance with Multithreaded Models of Flare Arcades

abstract: Understanding how energy is released in flares is one of the central problems of solar and stellar astrophysics. Observations of high temperature flare plasma hold many potential clues as to the nature of this energy release. It is clear, however, that flares are not composed of a few impulsively heated loops, but are the result of heating on many small-scale threads that are energized over time, making it difficult to compare observations and numerical simulations in detail. Several previous studies have shown that it is possible to reproduce some aspects of the observed emission by considering the flare as a sequence of independently heated loops, but these studies generally focus on small-scale features while ignoring the global features of the flare. In this paper, we develop a multithreaded model that encompasses the time-varying geometry and heating rate for a series of successively-heated loops comprising an arcade. To validate, we compare with spectral observations of five flares made with the MinXSS CubeSat as well as light curves measured with GOES/XRS and SDO/AIA. We show that this model can successfully reproduce the light curves and quasi-periodic pulsations in GOES/XRS, the soft X-ray spectra seen with MinXSS, and the light curves in various AIA passbands. The AIA light curves are most consistent with long duration heating, but elemental abundances cannot be constrained with the model. Finally, we show how this model can be used to extrapolate to spectra of extreme events that can predict irradiance across a wide wavelength range including unobserved wavelengths.


title: Peeking into occluded joints: A novel framework for crowd pose estimation

abstract: Although occlusion widely exists in nature and remains a fundamental challenge for pose estimation, existing heatmap-based approaches suffer serious degradation on occlusions. Their intrinsic problem is that they directly localize the joints based on visual information; however, the invisible joints are lack of that. In contrast to localization, our framework estimates the invisible joints from an inference perspective by proposing an Image-Guided Progressive GCN module which provides a comprehensive understanding of both image context and pose structure. Moreover, existing benchmarks contain limited occlusions for evaluation. Therefore, we thoroughly pursue this problem and propose a novel OPEC-Net framework together with a new Occluded Pose (OCPose) dataset with 9k annotated images. Extensive quantitative and qualitative evaluations on benchmarks demonstrate that OPEC-Net achieves significant improvements over recent leading works. Notably, our OCPose is the most complex occlusion dataset with respect to average IoU between adjacent instances. Source code and OCPose will be publicly available.


title: A Comparison of Discrete and Polyhedral Uncertainty Sets for Robust Network Design

abstract: We consider a network design and expansion problem, where we need to make a capacity investment now, such that uncertain future demand can be satisfied as closely as possible. To use a robust optimization approach, we need to construct an uncertainty set that contains all scenarios that we believe to be possible. In this paper we discuss how to construct two common types of uncertainty sets, which are discrete and polyhedral uncertainty, using real-world data. We employ clustering to generate a discrete uncertainty set, and place hyperplanes sequentially to generate a polyhedral uncertainty set. We then compare the performance of the resulting robust solutions for these two types of models on real-world data. Our results indicate that polyhedral models, while being popular in the recent network design literature, are less effective than discrete models both in terms of computational burden and solution quality with respect to the performance measure considered.


title: How do academic topics shift across altmetric sources? A case study of the research area of Big Data

abstract: Taking the research area of Big Data as a case study, we propose an approach for exploring how academic topics shift through the interactions among audiences across different altmetric sources. Data used is obtained from Web of Science (WoS) and Altmetric.com, with a focus on Blog, News, Policy, Wikipedia, and Twitter. Author keywords from publications and terms from online events are extracted as the main topics of the publications and the online discussion of their audiences at Altmetric. Different measures are applied to determine the (dis)similarities between the topics put forward by the publication authors and those by the online audiences. Results show that overall there are substantial differences between the two sets of topics around Big Data scientific research. The main exception is Twitter, where high-frequency hashtags in tweets have a stronger concordance with the author keywords in publications. Among the online communities, Blogs and News show a strong similarity in the terms commonly used, while Policy documents and Wikipedia articles exhibit the strongest dissimilarity in considering and interpreting Big Data related research. Specifically, the audiences not only focus on more easy-to-understand academic topics related to social or general issues, but also extend them to a broader range of topics in their online discussions. This study lays the foundations for further investigations about the role of online audiences in the transformation of academic topics across altmetric sources, and the degree of concern and reception of scholarly contents by online communities.


title: Defect dynamics in growing bacterial colonies

abstract: Colonies of rod-shaped bacteria constitute a system of colloidal active matter with nematic properties. As a single initial bacterium multiplies through repeated divisions, the resulting colony quickly loses long-range orientational order, but retains locally ordered domains. At the boundaries of these domains, topological defects emerge, which move around randomly as the colony grows. In both experiments and simulations, we find that these defects are created at a rate that corresponds to the exponential growth of the colony, resulting in a stable defect density. Both this defect density and the colony's correlation length are regulated by the aspect ratio of the rod-shaped particles. Moreover, we find that the defect dynamics are well described by a Gamma distribution, which is due to repeated divisions and subsequent re-orientations of the bacteria.


title: Accurate laser frequency locking to optical frequency combs under low-signal-to-noise-ratio conditions

abstract: We demonstrate a method for accurately locking the frequency of a continuous-wave laser to an optical frequency comb in conditions where the signal-to-noise ratio is low, too low to accommodate other methods. Our method is typically orders of magnitude more accurate than conventional wavemeters and can considerably extend the usable wavelength range of a given optical frequency comb. We illustrate our method by applying it to the frequency control of a dipole lattice trap for an optical lattice clock, a representative case where our method provides significantly better accuracy than other methods.


title: Efficiently computing logical noise in quantum error correcting codes

abstract: Quantum error correction protocols have been developed to offset the high sensitivity to noise inherent in quantum systems. However, much is still unknown about the behaviour of a quantum error-correcting code under general noise, including noisy measurements. This lack of knowledge is largely due to the computational cost of simulating quantum systems large enough to perform nontrivial encodings. In this paper, we develop general methods for incorporating noisy measurement operations into simulations of quantum error-correcting codes and show that measurement errors on readout qubits manifest as a renormalization on the effective logical noise. We also derive general methods for reducing the computational complexity of calculating the exact effective logical noise by many orders of magnitude. This reduction is achieved by determining when different recovery operations produce equivalent logical noise. These methods could also be used to better approximate soft decoding schemes for concatenated codes or to reduce the size of a lookup table to speed up the error correction step in implementations of quantum error-correcting codes. We give examples of such reductions for the three-qubit, five-qubit, Steane, concatenated, and toric codes.


title: Runaway electron synchrotron radiation in a vertically translated plasma

abstract: Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron diagnostic Soft. This experimental validation lends confidence to the theory underlying the synthetic diagnostics which are used for benchmarking theoretical models of and probing runaway dynamics. We present a comparison of synchrotron measurements in TCV with predictions of kinetic theory for runaway dynamics in uniform magnetic fields. We find that to explain the detected synchrotron emission, significant non-collisional pitch angle scattering as well as radial transport of REs would be needed. Such effects could be caused by the presence of magnetic perturbations, which should be further investigated in future TCV experiments.


title: Fault Attacks on Secure Embedded Software: Threats, Design and Evaluation

abstract: Embedded software is developed under the assumption that hardware execution is always correct. Fault attacks break and exploit that assumption. Through the careful introduction of targeted faults, an adversary modifies the control-flow or data-flow integrity of software. The modified program execution is then analyzed and used as a source of information leakage, or as a mechanism for privilege escalation. Due to the increasing complexity of modern embedded systems, and due to the difficulty of guaranteeing correct hardware execution even under a weak adversary, fault attacks are a growing threat. For example, the assumption that an adversary has to be close to the physical execution of software, in order to inject an exploitable fault into hardware, has repeatedly been shown to be incorrect. This article is a review on hardware-based fault attacks on software, with emphasis on the context of embedded systems. We present a detailed discussion of the anatomy of a fault attack, and we make a review of fault attack evaluation techniques. The paper emphasizes the perspective from the attacker, rather than the perspective of countermeasure development. However, we emphasize that improvements to countermeasures often build on insight into the attacks.


title: Independent components of human brain morphology

abstract: Quantification of brain morphology has become an important cornerstone in understanding brain structure. Measures of cortical morphology such as thickness and surface area are frequently used to compare groups of subjects or characterise longitudinal changes. However, such measures are often treated as independent from each other. A recently described scaling law, derived from a statistical physics model of cortical folding, demonstrates that there is a tight covariance between three commonly used cortical morphology measures: cortical thickness, total surface area, and exposed surface area. We show that assuming the independence of cortical morphology measures can hide features and potentially lead to misinterpretations. Using the scaling law, we account for the covariance between cortical morphology measures and derive novel independent measures of cortical morphology. By applying these new measures, we show that new information can be gained; in our example we show that distinct morphological alterations underlie healthy ageing compared to temporal lobe epilepsy, even on the coarse level of a whole hemisphere. We thus provide a conceptual framework for characterising cortical morphology in a statistically valid and interpretable manner, based on theoretical reasoning about the shape of the cortex.


title: Influence of Surface Hydrophilicity and Hydration on the Rotational Relaxation of Supercooled Water on Graphene Oxide Surfaces

abstract: Hydration or interfacial water present in biomolecules and inorganic solids have been shown to exhibit a dynamical transition upon supercooling. However, an understanding of the extent of the underlying surface hydrophilicity as well as the local distribution of hydrophilic/hydrophobic patches on the dynamical transition is unexplored. Here, we use molecular dynamics simulations with a TIP4P/2005 water model to study translational and rotational relaxation dynamics of interfacial water on graphene surfaces. The purpose of this study is to investigate the influence of both surface chemistry as well as the extent of hydration on the rotational transitions of interfacial water on graphene oxide (GO) surfaces in the deeply supercooled region. We have considered three graphene-based surfaces; a GO surface with equal proportions of oxidized and pristine graphene regions in a striped topology, a fully oxidized surface and a pristine graphene surface. The dipole relaxation time of interfacial water shows a strong-to-strong, strong, and a fragile-to-strong transition on these surfaces, respectively, in the temperature range of 210-298 K. In contrast, bulk water shows a fragile-to-strong transition upon supercooling. In all these cases at high hydration, interfacial water co-exists with a thick water film with bulk-like properties. To investigate the influence of bulk water on dynamical transitions, we simulated a low hydration regime where only bound water (surface water) is present on the GO surfaces and found that the rotational relaxation of surface water on both the GO and fully oxidized surfaces show a single Arrhenius behavior. Our results indicate that not only does the local extent of surface hydrophilicity play a role in determining the energy landscape explored by the water molecules upon supercooling, but the presence of bulk water also modulates the dynamic transition.


title: Global asymptotic dynamics of Cosmological Einsteinian Cubic Gravity

abstract: In this paper we investigate the cosmological dynamics of an up to cubic curvature correction to General Relativity (GR) known as Cosmological Einsteinian Cubic Gravity (CECG), whose vacuum spectrum consists of the graviton exclusively and its cosmology is well-posed as an initial value problem. We are able to uncover the global asymptotic structure of the phase space of this theory. It is revealed that an inflationary matter-dominated bigbang is the global past attractor which means that inflation is the starting point of any physically meaningful cosmic history. Given that higher order curvature corrections to GR are assumed to influence the cosmological dynamics at early times -- high energies/large curvature limit -- the late-time inflation can not be a consequence of the up to cubic order curvature modifications. We confirm this assumption by showing that late-time acceleration of the expansion in the CECG model is possible only if add a cosmological constant term.


title: Multivariate Matrix Mittag--Leffler distributions

abstract: We extend the construction principle of multivariate phase-type distributions to establish an analytically tractable class of heavy-tailed multivariate random variables whose marginal distributions are of Mittag-Leffler type with arbitrary index of regular variation. The construction can essentially be seen as allowing a scalar parameter to become matrix-valued. The class of distributions is shown to be dense among all multivariate positive random variables and hence provides a versatile candidate for the modelling of heavy-tailed, but tail-independent, risks in various fields of application.


title: Covid-19: Data analysis of the Lombardy region and the provinces of Bergamo and Brescia

abstract: The data analysis on deaths in the Lombardy Region and of both the provinces of Bergamo and Brescia shows a twofold aspect on the trend of the epidemic: - all the data show a bias linked to the event of March 10th (day for which the Lombardy region data is partial) and the subsequent change in the way in which positive cases and deaths are calculated; - following the containment measures of the Prime Minister's Decree of March 11th, the spread of the epidemic, although still exponential in nature, has a reduced multiplication coefficient. Our analysis concludes that the situation is not yet compatible with a plateau trend and allows us to predict the trend in the number of deaths in the Lombardy region. We therefore conclude that the containment measures put in place by the government on March 11th will allow a reduction in deaths from around 8000 to just over 6500 for March 27th.


title: On the Heisenberg limit for detecting vacuum birefringence

abstract: Quantum electrodynamics predicts the vacuum to behave as a non-linear medium, including effects such as birefringence. However, for experimentally available field strengths, this vacuum polarizability is extremely small and thus very hard to measure. In analogy to the Heisenberg limit in quantum metrology, we study the minimum requirements for such a detection in a given strong field (the pump field). Using a laser pulse as the probe field, we find that its energy must exceed a certain threshold depending on the interaction time. However, a detection at that threshold, i.e., the Heisenberg limit, requires highly non-linear measurement schemes - while for ordinary linear-optics schemes, the required energy (Poisson or shot noise limit) is much larger. Finally, we discuss several currently considered experimental scenarios from this point of view.


title: Neural Game Engine: Accurate learning of generalizable forward models from pixels

abstract: Access to a fast and easily copied forward model of a game is essential for model-based reinforcement learning and for algorithms such as Monte Carlo tree search, and is also beneficial as a source of unlimited experience data for model-free algorithms. Learning forward models is an interesting and important challenge in order to address problems where a model is not available. Building upon previous work on the Neural GPU, this paper introduces the Neural Game Engine, as a way to learn models directly from pixels. The learned models are able to generalise to different size game levels to the ones they were trained on without loss of accuracy. Results on 10 deterministic General Video Game AI games demonstrate competitive performance, with many of the games models being learned perfectly both in terms of pixel predictions and reward predictions. The pre-trained models are available through the OpenAI Gym interface and are available publicly for future research here: \url{https://github.com/Bam4d/Neural-Game-Engine}


title: An entrainment-based model for annular wakes, with applications to airborne wind energy

abstract: Several novel wind energy systems produce wakes with annular cross-sections, which are qualitatively different from the wakes with circular cross-sections commonly generated by conventional horizontal-axis wind turbines and by compact obstacles. Since wind farms use arrays of hundreds of turbines, good analytical wake models are essential for efficient wind farm planning. Several models already exist for circular wakes; however, none have yet been proposed for annular wakes, making it impossible to estimate their array performance. We use the entrainment hypothesis to develop a reduced-order model for the shape and flow velocity of an annular wake from a generic annular obstacle. Our model consists of a set of three ordinary differential equations, which we solve numerically. In addition, by assuming that the annular wake does not drift radially, we further reduce the problem to a model comprising only two differential equations, which we solve analytically. Both of our models are in good agreement with previously published large eddy simulation results.


title: On the preconditioning of three-by-three block saddle point problems

abstract: We establish a new iterative method for solving a class of large and sparse linear systems of equations with three-by-three block coefficient matrices having saddle point structure. Convergence properties of the proposed method are studied in details and its induced preconditioner is examined for accelerating the convergence speed of generalized minimal residual (GMRES) method. More precisely, we analyze the eigenvalue distribution of the preconditioned matrix. Numerical experiments are reported to demonstrate the effectiveness of the proposed preconditioner.


title: Neural Networks and Polynomial Regression. Demystifying the Overparametrization Phenomena

abstract: In the context of neural network models, overparametrization refers to the phenomena whereby these models appear to generalize well on the unseen data, even though the number of parameters significantly exceeds the sample sizes, and the model perfectly fits the in-training data. A conventional explanation of this phenomena is based on self-regularization properties of algorithms used to train the data. In this paper we prove a series of results which provide a somewhat diverging explanation. Adopting a teacher/student model where the teacher network is used to generate the predictions and student network is trained on the observed labeled data, and then tested on out-of-sample data, we show that any student network interpolating the data generated by a teacher network generalizes well, provided that the sample size is at least an explicit quantity controlled by data dimension and approximation guarantee alone, regardless of the number of internal nodes of either teacher or student network. Our claim is based on approximating both teacher and student networks by polynomial (tensor) regression models with degree depending on the desired accuracy and network depth only. Such a parametrization notably does not depend on the number of internal nodes. Thus a message implied by our results is that parametrizing wide neural networks by the number of hidden nodes is misleading, and a more fitting measure of parametrization complexity is the number of regression coefficients associated with tensorized data. In particular, this somewhat reconciles the generalization ability of neural networks with more classical statistical notions of data complexity and generalization bounds. Our empirical results on MNIST and Fashion-MNIST datasets indeed confirm that tensorized regression achieves a good out-of-sample performance, even when the degree of the tensor is at most two.


title: Transfer Learning in Automated Gamma Spectral Identification

abstract: The models and weights of prior trained Convolutional Neural Networks (CNN) created to perform automated isotopic classification of time-sequenced gamma-ray spectra, were utilized to provide source domain knowledge as training on new domains of potential interest. The previous results were achieved solely using modeled spectral data. In this work we attempt to transfer the knowledge gained to the new, if similar, domain of solely measured data. The ability to train on modeled data and predict on measured data will be crucial in any successful data-driven approach to this problem space.


title: Modelling Network Interference with Multi-valued Treatments: the Causal Effect of Immigration Policy on Crime Rates

abstract: Policy evaluation studies, which intend to assess the effect of an intervention, face some statistical challenges: in real-world settings treatments are not randomly assigned and the analysis might be further complicated by the presence of interference between units. Researchers have started to develop novel methods that allow to manage spillover mechanisms in observational studies; recent works focus primarily on binary treatments. However, many policy evaluation studies deal with more complex interventions. For instance, in political science, evaluating the impact of policies implemented by administrative entities often implies a multivariate approach, as a policy towards a specific issue operates at many different levels and can be defined along a number of dimensions. In this work, we extend the statistical framework about causal inference under network interference in observational studies, allowing for a multi-valued individual treatment and an interference structure shaped by a weighted network. The estimation strategy is based on a joint multiple generalized propensity score and allows one to estimate direct effects, controlling for both individual and network covariates. We follow the proposed methodology to analyze the impact of the national immigration policy on the crime rate. We define a multi-valued characterization of political attitudes towards migrants and we assume that the extent to which each country can be influenced by another country is modeled by an appropriate indicator, summarizing their cultural and geographical proximity. Results suggest that implementing a highly restrictive immigration policy leads to an increase of the crime rate and the estimated effects is larger if we take into account interference from other countries.


title: Hessian metric via transport information geometry

abstract: We propose to study the Hessian metric of a functional on the space of probability measures endowed with the Wasserstein $2$-metric. We name it transport Hessian metric, which contains and extends the classical Wasserstein-$2$ metric. We formulate several dynamical systems associated with transport Hessian metrics. Several connections between transport Hessian metrics and mathematical physics equations are discovered. E.g., the transport Hessian gradient flow, including Newton's flow, formulates a mean-field kernel Stein variational gradient flow; The transport Hessian Hamiltonian flow of Boltzmann-Shannon entropy forms the Shallow water equation; The transport Hessian gradient flow of Fisher information leads to the heat equation. Several examples and closed-form solutions for transport Hessian distances are presented.


title: Gravitational Waves and Possible Fast Radio Bursts from Axion Clumps

abstract: The axion objects such as axion mini-clusters and axion clouds around spinning black holes induce parametric resonances of electromagnetic waves through the axion-photon interaction. In particular, it has been known that the resonances from the axion with the mass around $10^{-6}$eV may explain the observed fast radio bursts (FRBs). Here we argue that similar bursts of high frequency gravitational waves, which we call the fast gravitational wave bursts (FGBs), are generated from axion clumps with the presence of gravitational Chern-Simons (CS) coupling. The typical frequency is half of the axion mass, which in general can range from kHz to GHz. We also discuss the secondary gravitational wave production associated with FRB, as well as the possible host objects of the axion clouds, such as primordial black holes with typical masses around $10^{-5}M_{\odot}$. Future detections of FGBs together with the observed FRBs are expected to provide more evidence for the axion.


title: A Simultaneous Inference Procedure to Identify Subgroups from RCTs with Survival Outcomes: Application to Analysis of AMD Progression Studies

abstract: With the uptake of targeted therapies, instead of the "one-fits-all" approach, modern randomized clinical trials (RCTs) often aim to develop treatments that target a subgroup of patients. Motivated by analyzing the Age-Related Eye Disease Study (AREDS) data, a large RCT to study the efficacy of nutritional supplements in delaying the progression of an eye disease, age-related macular degeneration (AMD), we develop a simultaneous inference procedure to identify and infer subgroups with differential treatment efficacy in RCTs with survival outcome. Specifically, we formulate the multiple testing problem through contrasts and construct their simultaneous confidence intervals, which control both within- and across- marker multiplicity appropriately. Realistic simulations are conducted using real genotype data to evaluate the method performance under various scenarios. The method is then applied to AREDS to assess the efficacy of antioxidants and zinc combination in delaying AMD progression. Multiple gene regions including ESRRB-VASH1 on chromosome 14 have been identified with subgroups showing differential efficacy. We further validate our findings in an independent subsequent RCT, AREDS2, by discovering consistent differential treatment responses in the targeted and non-targeted subgroups been identified from AREDS. This simultaneous inference approach provides a step forward to confidently identify and infer subgroups in modern drug development.


title: Generic symmetry-forced infinitesimal rigidity: translations and rotations

abstract: We characterize the combinatorial types of symmetric frameworks in the plane that are minimally generically symmetry-forced infinitesimally rigid when the symmetry group consists of rotations and translations. Along the way, we use tropical geometry to show how a construction of Edmonds that associates a matroid to a submodular function can be used to give a description of the algebraic matroid of a Hadamard product of two linear spaces in terms of the matroids of each linear space. This leads to new, short, proofs of Laman's theorem, and a theorem of Jord{\'a}n, Kaszanitzky, and Tanigawa, and Malestein and Theran characterizing the minimally generically symmetry-forced rigid graphs in the plane when the symmetry group contains only rotations.


title: Understanding the fate of corona virus transmission using a simple model

abstract: We propose a simple model for understanding the kinetics of corona virus transmission. Our model assume spreading of corona virus can happen from one to another only, if someone without enough protection comes close contact to a person carrying the corona virus. Therefore this virus spreads on a large scale within a short time through chains of such events. Using our model we provide an estimation of the number of people affected by this virus within reasonable duration of time. We choose values of different parameters of our model by non-linear least square fit of the real time data and we predict fate of this corona virus transmission using our model.


title: Crowdsourced Smartphone Sensing for Localization in Metro Trains

abstract: Traditional fingerprint based localization techniques mainly rely on infrastructure support such as RFID, Wi-Fi or GPS. They operate by war-driving the entire space which is both time-consuming and labor-intensive. In this paper, we present MLoc, a novel infrastructure-free localization system to locate mobile users in a metro line. It does not rely on any Wi-Fi infrastructure, and does not need to war-drive the metro line. Leveraging crowdsourcing, we collect accelerometer,magnetometer and barometer readings on smartphones, and analyze these sensor data to extract patterns. Through advanced data manipulating techniques, we build the pattern map for the entire metro line, which can then be used for localization. We conduct field studies to demonstrate the accuracy, scalability, and robustness of M-Loc. The results of our field studies in 3 metro lines with 55 stations show that M-Loc achieves an accuracy of 93% when travelling 3 stations, 98% when travelling 5 stations.


title: A simple iterative map forecast of the COVID-19 pandemic

abstract: We develop a simple 3-dimensional iterative map model to forecast the global spread of the coronavirus disease. Our model contains at most two fitting parameters, which we determine from the data supplied by the world health organisation for the total number of cases and new cases each day. We find that our model provides a surprisingly good fit to the currently-available data, which exhibits a cross-over from exponential to power-law growth, as lock-down measures begin to take effect. Before these measures, our model predicts exponential growth from day 30 to 69, starting from the date on which the world health organisation provided the first `Situation report' (21 January 2020 $-$ day 1). Based on this initial data the disease may be expected to infect approximately 23% of the global population, i.e. about 1.76 billion people, taking approximately 83 million lives. Under this scenario, the global number of new cases is predicted to peak on day 133 (about the middle of May 2020), with an estimated 60 million new cases per day. If current lock-down measures can be maintained, our model predicts power law growth from day 69 onward. Such growth is comparatively slow and would have to continue for several decades before a sufficient number of people (at least 23% of the global population) have developed immunity to the disease through being infected. Lock-down measures appear to be very effective in postponing the unimaginably large peak in the daily number of new cases that would occur in the absence of any interventions. However, should these measure be relaxed, the spread of the disease will most likely revert back to its original exponential growth pattern. As such, the duration and severity of the lock-down measures should be carefully timed against their potentially devastating impact on the world economy.


title: Formation of the first stars and black holes

abstract: We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.


title: A new paradigm for accelerating clinical data science at Stanford Medicine

abstract: Stanford Medicine is building a new data platform for our academic research community to do better clinical data science. Hospitals have a large amount of patient data and researchers have demonstrated the ability to reuse that data and AI approaches to derive novel insights, support patient care, and improve care quality. However, the traditional data warehouse and Honest Broker approaches that are in current use, are not scalable. We are establishing a new secure Big Data platform that aims to reduce time to access and analyze data. In this platform, data is anonymized to preserve patient data privacy and made available preparatory to Institutional Review Board (IRB) submission. Furthermore, the data is standardized such that analysis done at Stanford can be replicated elsewhere using the same analytical code and clinical concepts. Finally, the analytics data warehouse integrates with a secure data science computational facility to support large scale data analytics. The ecosystem is designed to bring the modern data science community to highly sensitive clinical data in a secure and collaborative big data analytics environment with a goal to enable bigger, better and faster science.


title: Real forms of complex Lie superalgebras and supergroups

abstract: We investigate the notion of real form of complex Lie superalgebras and supergroups, both in the standard and graded version. Our functorial approach allows most naturally to go from the superalgebra to the supergroup and retrieve the real forms as fixed points, as in the ordinary setting. We also introduce a more general notion of compact real form for Lie superalgebras and supergroups, and we prove some existence results for Lie superalgebras that are simple contragredient and their associated connected simply connected supergroups.


title: ProGraML: Graph-based Deep Learning for Program Optimization and Analysis

abstract: The increasing complexity of computing systems places a tremendous burden on optimizing compilers, requiring ever more accurate and aggressive optimizations. Machine learning offers significant benefits for constructing optimization heuristics but there remains a gap between what state-of-the-art methods achieve and the performance of an optimal heuristic. Closing this gap requires improvements in two key areas: a representation that accurately captures the semantics of programs, and a model architecture with sufficient expressiveness to reason about this representation. We introduce ProGraML - Program Graphs for Machine Learning - a novel graph-based program representation using a low level, language agnostic, and portable format; and machine learning models capable of performing complex downstream tasks over these graphs. The ProGraML representation is a directed attributed multigraph that captures control, data, and call relations, and summarizes instruction and operand types and ordering. Message Passing Neural Networks propagate information through this structured representation, enabling whole-program or per-vertex classification tasks. ProGraML provides a general-purpose program representation that equips learnable models to perform the types of program analysis that are fundamental to optimization. To this end, we evaluate the performance of our approach first on a suite of traditional compiler analysis tasks: control flow reachability, dominator trees, data dependencies, variable liveness, and common subexpression detection. On a benchmark dataset of 250k LLVM-IR files covering six source programming languages, ProGraML achieves an average 94.0 F1 score, significantly outperforming the state-of-the-art approaches. We then apply our approach to two high-level tasks - heterogeneous device mapping and program classification - setting new state-of-the-art performance in both.


title: Higher order singular value decomposition and the reduced density matrices of three qubits

abstract: In this paper, we demonstrate that higher order singular value decomposition (HOSVD) can be used to identify special states in three qubits by local unitary (LU) operations. Since the matrix unfoldings of three qubits are related to their reduced density matrices, HOSVD simultaneously diagonalizes the one-body reduced density matrices of three qubits. From the all-orthogonality conditions of HOSVD, we computed the special states of three qubits. Furthermore, we showed that it is possible to construct a polytope that encapsulates all the special states of three qubits by LU operations with HOSVD.


title: The Internet of Things as a Deep Neural Network

abstract: An important task in the Internet of Things (IoT) is field monitoring, where multiple IoT nodes take measurements and communicate them to the base station or the cloud for processing, inference, and analysis. This communication becomes costly when the measurements are high-dimensional (e.g., videos or time-series data). The IoT networks with limited bandwidth and low power devices may not be able to support such frequent transmissions with high data rates. To ensure communication efficiency, this article proposes to model the measurement compression at IoT nodes and the inference at the base station or cloud as a deep neural network (DNN). We propose a new framework where the data to be transmitted from nodes are the intermediate outputs of a layer of the DNN. We show how to learn the model parameters of the DNN and study the trade-off between the communication rate and the inference accuracy. The experimental results show that we can save approximately 96% transmissions with only a degradation of 2.5% in inference accuracy. Our findings have the potentiality to enable many new IoT data analysis applications generating large amount of measurements.


title: The topological period-index conjecture

abstract: We prove the topological analogue of the period-index conjecture in each dimension away from a small set of primes.


title: Data-driven models and computational tools for neurolinguistics: a language technology perspective

abstract: In this paper, our focus is the connection and influence of language technologies on the research in neurolinguistics. We present a review of brain imaging-based neurolinguistic studies with a focus on the natural language representations, such as word embeddings and pre-trained language models. Mutual enrichment of neurolinguistics and language technologies leads to development of brain-aware natural language representations. The importance of this research area is emphasized by medical applications.


title: Physical realization of complex dynamical pattern formation in magnetic active feedback rings

abstract: We report the clean experimental realization of cubic-quintic complex Ginzburg-Landau physics in a single driven, damped system. Four numerically predicted categories of complex dynamical behavior and pattern formation are identified for bright and dark solitary waves propagating around an active magnetic thin film-based feedback ring: (1) periodic breathing; (2) complex recurrence; (3) spontaneous spatial shifting; and (4) intermittency. These nontransient, long lifetime behaviors are observed in microwave spin wave envelopes circulating within a dispersive, nonlinear yttrium iron garnet waveguide operating in a ring geometry where the net losses are directly compensated for via linear amplification on each round trip O(100 ns). The behaviors exhibit periods ranging from tens to thousands of round trip times O($\mu$ s) and are stable for 1000s of periods O(ms). We present 10 observations of these dynamical behaviors which span the experimentally accessible ranges of attractive cubic nonlinearity, dispersion, and external field strength that support the self-generation of backward volume spin waves in a four-wave-mixing dominant regime. Three-wave splitting is not explicitly forbidden and is treated as an additional source of nonlinear losses. Each of these long lifetime behaviors of bright solitary waves was previously numerically predicted to be observable. The dynamical pattern formation of dark solitary waves in attractive nonlinearity, however, is entirely novel and is reported for both the periodic breather and complex recurrence behaviors. All behaviors are identified in the group velocity co-moving frame.


title: Source Depth of Three-minute Umbral Oscillations

abstract: We infer the depth of the internal sources giving rise to three-minute umbral oscillations. Recent observations of ripple-like velocity patterns of umbral oscillations supported the notion that there exist internal sources exciting the umbral oscillations. We adopt the hypothesis that the fast magnetohydrodynamic (MHD) waves generated at a source below the photospheric layer propagate along different paths, reach the surface at different times, and excited slow MHD waves by mode conversion. These slow MHD waves are observed as the ripples that apparently propagate horizontally. The propagation distance of the ripple given as a function of time is strongly related to the depth of the source. Using the spectral data of the Fe I 5435A line taken by the Fast Imaging Solar Spectrograph of the Goode Solar Telescope at Big Bear Solar Observatory, we identified five ripples and determined the propagation distance as a function of time in each ripple. From the model fitting to these data, we obtained the depth between 1000 km and 2000 km. Our result will serve as an observational constraint to understanding the detailed processes of magnetoconvection and wave generation in sunspots.


title: Automated Detection of Cribriform Growth Patterns in Prostate Histology Images

abstract: Cribriform growth patterns in prostate carcinoma are associated with poor prognosis. We aimed to introduce a deep learning method to detect such patterns automatically. To do so, convolutional neural network was trained to detect cribriform growth patterns on 128 prostate needle biopsies. Ensemble learning taking into account other tumor growth patterns during training was used to cope with heterogeneous and limited tumor tissue occurrences. ROC and FROC analyses were applied to assess network performance regarding detection of biopsies harboring cribriform growth pattern. The ROC analysis yielded a mean area under the curve up to 0.81. FROC analysis demonstrated a sensitivity of 0.9 for regions larger than 0.0150 mm2 with on average 7.5 false positives. To benchmark method performance for intra-observer annotation variability, false positive and negative detections were re-evaluated by the pathologists. Pathologists considered 9% of the false positive regions as cribriform, and 11% as possibly cribriform; 44% of the false negative regions were not annotated as cribriform. As a final experiment, the network was also applied on a dataset of 60 biopsy regions annotated by 23 pathologists. With the cut-off reaching highest sensitivity, all images annotated as cribriform by at least 7/23 of the pathologists, were all detected as cribriform by the network and 9/60 of the images were detected as cribriform whereas no pathologist labelled them as such. In conclusion, the proposed deep learning method has high sensitivity for detecting cribriform growth patterns at the expense of a limited number of false positives. It can detect cribriform regions that are labelled as such by at least a minority of pathologists. Therefore, it could assist clinical decision making by suggesting suspicious regions.


title: A Multiple Power Law Distribution for Initial Mass Functions

abstract: We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framework to demonstrate that the interplay between a mass-dependent and a time-dependent step-like dropout rate from accretion leads to IMFs that exhibit multiple power laws for an exponential mass growth. While the mass-dependent accretion and its dropout is intrinsic to each star, the time-dependent dropout might be tied to a specific history such as the rapid consumption of nebular material by nearby stars or the sweeping away of some material by shock waves. The time-dependent dropout folded into the mass-dependent process of star formation is shown to have a significant influence on the IMFs.


title: Potassium Isotope Compositions of Carbonaceous and Ordinary Chondrites: Implications on the Origin of Volatile Depletion in the Early Solar System

abstract: Solar system materials are variably depleted in moderately volatile elements (MVEs) relative to the proto-solar composition. To address the origin of this MVE depletion, we conducted a systematic study of high-precision K isotopic composition on 16 carbonaceous chondrites (CCs) of types CM1-2, CO3, CV3, CR2, CK4-5 and CH3 and 28 ordinary chondrites (OCs) covering petrological types 3 to 6 and chemical groups H, L, and LL. We observed significant overall K isotope (delta41K) variations (-1.54 to 0.70 permil). The K isotope compositions of CCs are largely higher than the Bulk Silicate Earth (BSE) value, whereas OCs show typically lower values than BSE. Neither CCs nor OCs show resolvable correlations between K isotopes and chemical groups, petrological types, shock levels, exposure ages, fall or find occurrence, or terrestrial weathering. The lack of a clear trend between K isotopes and K content indicates that the K isotope fractionations were decoupled from the relative elemental K depletions. The range of K isotope variations in the CCs is consistent with a four-component (chondrule, refractory inclusion, matrix and water) mixing model that is able to explain the bulk elemental and isotopic compositions of the main CC groups, but requires a fractionation in K isotopic compositions in chondrules. We propose that the major control of the isotopic compositions of group averages is condensation or vaporization in nebular environments that is preserved in the compositional variation of chondrules. Parent-body processes (aqueous alteration, thermal metamorphism, and metasomatism) can mobilize K and affect the K isotopes in individual samples. In the case of the OCs, the full range of K isotopic variations can only be explained by the combined effects of the size and relative abundances of chondrules, parent-body aqueous and thermal alteration.


title: Forensic Analysis of Residual Information in Adobe PDF Files

abstract: In recent years, as electronic files include personal records and business activities, these files can be used as important evidences in a digital forensic investigation process. In general, the data that can be verified using its own application programs is largely used in the investigation of document files. However, in the case of the PDF file that has been largely used at the present time, certain data, which include the data before some modifications, exist in electronic document files unintentionally. Because such residual information may present the writing process of a file, it can be usefully used in a forensic viewpoint. This paper introduces why the residual information is stored inside the PDF file and explains a way to extract the information. In addition, we demonstrate the attributes of PDF files can be used to hide data.


title: An advanced hybrid deep adversarial autoencoder for parameterized nonlinear fluid flow modelling

abstract: Considering the high computation cost produced in conventional computation fluid dynamic simulations, machine learning methods have been introduced to flow dynamic simulations in recent years. However, most of studies focus mainly on existing fluid fields learning, the prediction of spatio-temporal nonlinear fluid flows in varying parameterized space has been neglected. In this work, we propose a hybrid deep adversarial autoencoder (DAA) to integrate generative adversarial network (GAN) and variational autoencoder (VAE) for predicting parameterized nonlinear fluid flows in spatial and temporal space. High-dimensional inputs are compressed into the low-representation representations by nonlinear functions in a convolutional encoder. In this way, the predictive fluid flows reconstructed in a convolutional decoder contain the dynamic flow physics of high nonlinearity and chaotic nature. In addition, the low-representation representations are applied into the adversarial network for model training and parameter optimization, which enables a fast computation process. The capability of the hybrid DAA is demonstrated by varying inputs on a water collapse example. Numerical results show that this hybrid DAA has successfully captured the spatio-temporal flow features with CPU speed-up of three orders of magnitude. Promising results suggests that the hybrid DAA can play a critical role in efficiently and accurately predicting complex flows in future.


title: spsurv: An R package for semi-parametric survival analysis

abstract: Software development innovations and advances in computing have enabled more complex and less costly computations in medical research (survival analysis), engineering studies (reliability analysis), and social sciences event analysis (historical analysis). As a result, many semi-parametric modeling efforts emerged when it comes to time-to-event data analysis. In this context, this work presents a flexible Bernstein polynomial (BP) based framework for survival data modeling. This innovative approach is applied to existing families of models such as proportional hazards (PH), proportional odds (PO), and accelerated failure time (AFT) models to estimate unknown baseline functions. Along with this contribution, this work also presents new automated routines in R, taking advantage of algorithms available in Stan. The proposed computation routines are tested and explored through simulation studies based on artificial datasets. The tools implemented to fit the proposed statistical models are combined and organized in an R package. Also, the BP based proportional hazards (BPPH), proportional odds (BPPO), and accelerated failure time (BPAFT) models are illustrated in real applications related to cancer trial data using maximum likelihood (ML) estimation and Markov chain Monte Carlo (MCMC) methods.


title: The influence of initial perturbation power spectra on the growth of a turbulent mixing layer induced by Richtmyer-Meshkov instability

abstract: This paper investigates the influence of different broadband perturbations on the evolution of a Richtmyer--Meshkov turbulent mixing layer initiated by a Mach 1.84 shock traversing a perturbed interface separating gases with a density ratio of 3:1. Both the bandwidth of modes in the interface perturbation, as well as their relative amplitudes, are varied in a series of carefully designed numerical simulations at grid resolutions up to $3.2\times10^9$ cells. Three different perturbations are considered, characterised by a power spectrum of the form $P(k)\propto k^m$ where $m=-1$, $-2$ and $-3$. The growth of the mixing layer is shown to strongly depend on the initial conditions, with the growth rate exponent $\theta$ found to be $0.5$, $0.63$ and $0.75$ for each value of $m$ at the highest grid resolution. The asymptotic values of the molecular mixing fraction $\Theta$ are also shown to vary significantly with $m$; at the latest time considered $\Theta$ is $0.56$, $0.39$ and $0.20$ respectively. Turbulent kinetic energy (TKE) is also analysed in both the temporal and spectral domains. The temporal decay rate of TKE is found not to match the predicted value of $n=2-3\theta$, which is shown to be due to a time-varying {normalised dissipation rate $C_\epsilon$}. In spectral space, the data follow the theoretical scaling of $k^{(m+2)/2}$ at low wavenumbers and tend towards $k^{-3/2}$ and $k^{-5/3}$ scalings at high wavenumbers for the spectra of transverse and normal velocity components respectively. The results represent a significant extension of previous work on the Richtmyer--Meshkov instability evolving from broadband initial perturbations and provide useful benchmarks for future research.


title: Regret and Belief Complexity Trade-off in Gaussian Process Bandits via Information Thresholding

abstract: Bayesian optimization is a framework for global search via maximum a posteriori updates rather than simulated annealing, and has gained prominence for decision-making under uncertainty. In this work, we cast Bayesian optimization as a multi-armed bandit problem, where the payoff function is sampled from a Gaussian process (GP). Further, we focus on action selections via upper confidence bound (UCB) or expected improvement (EI) due to their prevalent use in practice. Prior works using GPs for bandits cannot allow the iteration horizon $T$ to be large, as the complexity of computing the posterior parameters scales cubically with the number of past observations. To circumvent this computational burden, we propose a simple statistical test: only incorporate an action into the GP posterior when its conditional entropy exceeds an $\epsilon$ threshold. Doing so permits us to precisely characterize the trade-off between regret bounds of GP bandit algorithms and complexity of the posterior distributions depending on the compression parameter $\epsilon$ for both discrete and continuous action sets. To best of our knowledge, this is the first result which allows us to obtain sublinear regret bounds while still maintaining sublinear growth rate of the complexity of the posterior which is linear in the existing literature. Moreover, a provably finite bound on the complexity could be achieved but the algorithm would result in $\epsilon$-regret which means $\textbf{Reg}_T/T \rightarrow \mathcal{O}(\epsilon)$ as $T\rightarrow \infty$. Experimentally, we observe state of the art accuracy and complexity trade-offs for GP bandit algorithms applied to global optimization, suggesting the merits of compressed GPs in bandit settings.


title: G-Net: A Deep Learning Approach to G-computation for Counterfactual Outcome Prediction Under Dynamic Treatment Regimes

abstract: Counterfactual prediction is a fundamental task in decision-making. G-computation is a method for estimating expected counterfactual outcomes under dynamic time-varying treatment strategies. Existing G-computation implementations have mostly employed classical regression models with limited capacity to capture complex temporal and nonlinear dependence structures. This paper introduces G-Net, a novel sequential deep learning framework for G-computation that can handle complex time series data while imposing minimal modeling assumptions and provide estimates of individual or population-level time varying treatment effects. We evaluate alternative G-Net implementations using realistically complex temporal simulated data obtained from CVSim, a mechanistic model of the cardiovascular system.


title: Axion oscillations in binary systems: angle-action surgery

abstract: Scalar, tensor waves induce oscillatory perturbations in Keplerian systems that can be probed with measurements of pulsar timing residuals. In this paper, we consider the imprint of coherent oscillations produced by ultralight axion dark matter on the Roemer time delay. We use the angle-action formalism to calculate the time evolution of the observed signal and its dependence on the orbital parameters and the axion phase. We derive exact analytical expressions for arbitrary binary pulsar mass ratio and eccentricity, alleviating the need for long numerical integrations. We emphasize the similarity of the expected signal-to-noise ratio with the response of a harmonic oscillator to an external oscillatory driving. We validate our theoretical predictions with numerical simulations. Our results furnish a useful benchmark for numerical codes and analysis procedures and, hopefully, will motivate the search for such imprints in real data.

categories: gr-qc astro-ph.CO astro-ph.GA astro-ph.HE hep-ph


title: RoboMem: Giving Long Term Memory to Robots

abstract: Robots have the potential to improve health monitoring outcomes for the elderly by providing doctors, and caregivers with information about the person's behavior, health activities and their surrounding environment. Over the years, less work has been done to enable robots to preserve information for longer periods of time, on the order of months and years of data, and use this contextual information to answer queries. Time complexity to process this massive sensor data in a timely fashion, inability to anticipate the future queries in advance and imprecision involved in the results have been the main impediments in making progress in this area. We make a contribution by introducing RoboMem, a query answering system for health-care assistance of elderly over long term; continuous data feeds that intends to overcome the challenges of giving long term memory to robots. The design for our framework preprocesses the sensor data and stores this preprocessed data into the database. This data is updated in the database by going through successive refinements, improving its accuracy for responding to queries. If data in the database is not enough to answer a query, a small set of relevant frames (also obtained from the database) will be reprocessed to obtain the answer. [Our initial prototype of RoboMem stores 3.5MB of data in the database as compared to 535.8MB of actual video frames and with minimal data in the database it is able to fetch information fundamental to respond to queries in 0.0002 seconds on average].


title: Elle: Inferring Isolation Anomalies from Experimental Observations

abstract: Users who care about their data store it in databases, which (at least in principle) guarantee some form of transactional isolation. However, experience shows [Kleppmann 2019, Kingsbury and Patella 2019a] that many databases do not provide the isolation guarantees they claim. With the recent proliferation of new distributed databases, demand has grown for checkers that can, by generating client workloads and injecting faults, produce anomalies that witness a violation of a stated guarantee. An ideal checker would be sound (no false positives), efficient (polynomial in history length and concurrency), effective (finding violations in real databases), general (analyzing many patterns of transactions), and informative (justifying the presence of an anomaly with understandable counterexamples). Sadly, we are aware of no checkers that satisfy these goals. We present Elle: a novel checker which infers an Adya-style dependency graph between client-observed transactions. It does so by carefully selecting database objects and operations when generating histories, so as to ensure that the results of database reads reveal information about their version history. Elle can detect every anomaly in Adya et al's formalism [Adya et al. 2000] (except for predicates), discriminate between them, and provide concise explanations of each. This paper makes the following contributions: we present Elle, demonstrate its soundness, measure its efficiency against the current state of the art, and give evidence of its effectiveness via a case study of four real databases.


title: ELECTRA: Pre-training Text Encoders as Discriminators Rather Than Generators

abstract: Masked language modeling (MLM) pre-training methods such as BERT corrupt the input by replacing some tokens with [MASK] and then train a model to reconstruct the original tokens. While they produce good results when transferred to downstream NLP tasks, they generally require large amounts of compute to be effective. As an alternative, we propose a more sample-efficient pre-training task called replaced token detection. Instead of masking the input, our approach corrupts it by replacing some tokens with plausible alternatives sampled from a small generator network. Then, instead of training a model that predicts the original identities of the corrupted tokens, we train a discriminative model that predicts whether each token in the corrupted input was replaced by a generator sample or not. Thorough experiments demonstrate this new pre-training task is more efficient than MLM because the task is defined over all input tokens rather than just the small subset that was masked out. As a result, the contextual representations learned by our approach substantially outperform the ones learned by BERT given the same model size, data, and compute. The gains are particularly strong for small models; for example, we train a model on one GPU for 4 days that outperforms GPT (trained using 30x more compute) on the GLUE natural language understanding benchmark. Our approach also works well at scale, where it performs comparably to RoBERTa and XLNet while using less than 1/4 of their compute and outperforms them when using the same amount of compute.


title: On the dynamics of a seventh-order generalized H\'enon-Heiles potential

abstract: This paper deals with the derivation and analysis of a seventh-order generalization of the H\'enon-Heiles potential. The new potential has axial and reflection symmetries, and finite escape energy with three channels of escape. Based on SALI indicator and exits basins, the dynamic behavior of the seventh-order system is investigated qualitatively in cases of bounded and unbounded movement. Moreover, a quantitative analysis is carried out through the percentage of chaotic orbits and the basin entropy, respectively. After classifying large sets of initial conditions of orbits for several values of the energy constant in both regimes, we observe that when the energy moves away from the critical value, the chaoticity of the system decreases and the basin structure becomes simpler with sharper and well defined bounds. Our results suggest that when the seventh-order contributions of the potential are taken into account, the system becomes less ergodic in comparison with the classical version of the H\'enon-Heiles system.


title: ScrabbleGAN: Semi-Supervised Varying Length Handwritten Text Generation

abstract: Optical character recognition (OCR) systems performance have improved significantly in the deep learning era. This is especially true for handwritten text recognition (HTR), where each author has a unique style, unlike printed text, where the variation is smaller by design. That said, deep learning based HTR is limited, as in every other task, by the number of training examples. Gathering data is a challenging and costly task, and even more so, the labeling task that follows, of which we focus here. One possible approach to reduce the burden of data annotation is semi-supervised learning. Semi supervised methods use, in addition to labeled data, some unlabeled samples to improve performance, compared to fully supervised ones. Consequently, such methods may adapt to unseen images during test time. We present ScrabbleGAN, a semi-supervised approach to synthesize handwritten text images that are versatile both in style and lexicon. ScrabbleGAN relies on a novel generative model which can generate images of words with an arbitrary length. We show how to operate our approach in a semi-supervised manner, enjoying the aforementioned benefits such as performance boost over state of the art supervised HTR. Furthermore, our generator can manipulate the resulting text style. This allows us to change, for instance, whether the text is cursive, or how thin is the pen stroke.


title: Perspective picture from Visual Sphere: a new approach to image rasterization

abstract: In this paper alternative method for real-time 3D model rasterization is given. Surfaces are drawn in perspective-map space which acts as a virtual camera lens. It can render single-pass 360{\deg} angle of view (AOV) image of unlimited shape, view-directions count and unrestrained projection geometry (e.g. direct lens distortion, projection mapping, curvilinear perspective), natively aliasing-free. In conjunction to perspective vector map, visual-sphere perspective model is proposed. A model capable of combining pictures from sources previously incompatible, like fish-eye camera and wide-angle lens picture. More so, method is proposed for measurement and simulation of a real optical system variable no-parallax point (NPP). This study also explores philosophical and historical aspects of picture perception and presents a guide for perspective design.


title: Asymptotic theory of quantum channel estimation

abstract: The quantum Fisher information (QFI), as a function of quantum states, measures the amount of information that a quantum state carries about an unknown parameter. The (entanglement-assisted) QFI of a quantum channel is defined to be the maximum QFI of the output state assuming an entangled input state over a single probe and an ancilla. In quantum metrology, people are interested in calculating the QFI of $N$ identical copies of a quantum channel when $N \rightarrow \infty$, which is called the asymptotic QFI. Over the years, researchers found various types of upper bounds of the asymptotic QFI, but they were proven achievable only in several specific situations. It was known that the asymptotic QFI of an arbitrary quantum channel grows either linearly or quadratically with $N$. Here we show that a simple criterion can determine whether the scaling is linear or quadratic. In both cases, the asymptotic QFI and a quantum error correction protocol to achieve it are computable via a semidefinite program. When the scaling is quadratic, the Heisenberg limit, a feature of noiseless quantum channels, is recovered. When the scaling is linear, we show the asymptotic QFI is still in general larger than $N$ times the single-channel QFI and furthermore, sequential estimation strategies provide no advantage over parallel ones.


title: Attract More Miners to Join in Blochchain Construction for Internet of Things

abstract: The world-changing blockchain technique provides a novel method to establish a secure, trusted and decentralized system for solving the security and personal privacy problems in Industrial Internet of Things (IIoT) applications. The mining process in blockchain requires miners to solve a proof-of-work puzzle, which requires high computational power. However, the lightweight IIoT devices cannot directly participate in the mining process due to the limitation of power and computational resources. The edge computing service makes it possible for IIoT applications to build a blockchain network, in which IIoT devices purchase computational resources from edge servers and thus can offload their computational tasks. The amount of computational resource purchased by IIoT devices depends on how many profits they can get in the mining process, and will directly affect the security of the blockchain network. In this paper, we investigate the incentive mechanism for the blockchain platform to attract IIoT devices to purchase more computational power from edge servers to participate in the mining process, thereby building a more secure blockchain network. We model the interaction between the blockchain platform and IIoT devices as a two-stage Stackelberg game, where the blockchain platform act as the leader, and IIoT devices act as followers. We analyze the existence and uniqueness of the Stackelberg equilibrium, and propose an efficient algorithm to compute the Stackelberg equilibrium point. Furthermore, we evaluate the performance of our algorithm through extensive simulations, and analyze the strategies of blockchain platform and IIoT devices under different situations.


title: Direct Numerical Simulations of turbulent flows using high-order Asynchrony-Tolerant schemes: accuracy and performance

abstract: Direct numerical simulations (DNS) are an indispensable tool for understanding the fundamental physics of turbulent flows. Because of their steep increase in computational cost with Reynolds number ($R_{\lambda}$), well-resolved DNS are realizable only on massively parallel supercomputers, even at moderate $R_{\lambda}$. However, at extreme scales, the communications and synchronizations between processing elements (PEs) involved in current approaches become exceedingly expensive and are expected to be a major bottleneck to scalability. In order to overcome this challenge, we developed algorithms using the so-called Asynchrony-Tolerant (AT) schemes that relax communication and synchronization constraints at a mathematical level, to perform DNS of decaying and solenoidally forced compressible turbulence. Asynchrony is introduced using two approaches, one that avoids synchronizations and the other that avoids communications. These result in periodic and random delays, respectively, at PE boundaries. We show that both asynchronous algorithms accurately resolve the large-scale and small-scale motions of turbulence, including instantaneous and intermittent fields. We also show that in asynchronous simulations the communication time is a relatively smaller fraction of the total computation time, especially at large processor count, compared to standard synchronous simulations. As a consequence, we observe improved parallel scalability up to $262144$ processors for both asynchronous algorithms.


title: The dual properties of chiral and isospin asymmetric dense quark matter formed of two-color quarks

abstract: In this paper the phase structure of dense baryon matter composed of $u$ and $d$ quarks with two colors has been investigated in the presence of baryon $\mu_B$, isospin $\mu_I$ and chiral isospin $\mu_{I5}$ chemical potentials in the framework of Nambu--Jona-Lasinio model. In the chiral limit, it has been shown that the duality between phases with spontaneous chiral symmetry breaking and condensation of charged pions, found in the three color case, remains valid in the two color case. In addition, it has been shown that there are two more dualities in the phase diagram in two color case, namely (as in the case with $\mu_{I5}=0$), at $\mu_{I5}\ne 0$ the general $(\mu,\mu_I,\mu_{I5})$-phase portrait of the model has dual symmetry between the phase with condensation of charged pions and the phase with diquark condensation. This duality stays exact even in the physical point, $m_0\ne 0$. And at $m_0=0$ the $(\mu,\mu_I,\mu_{I5})$-phase portrait becomes even more symmetrical, since dual symmetry between phases with spontaneous chiral symmetry breaking and diquark condensation appears. It is shown that due to the dualities the phase diagram is extremely symmetric and has interlacing structure. One can show that the phase portrait of two-color NJL model can be obtained just by duality properties from the results of investigations of three-color NJL model (it was noticed only after the numerical calculations have been performed). Three-color case shares only one duality of the two color one, and one can only see a facet of this enormously symmetric picture in the case of three colors. Using dualities only, it is possible to show that there are no mixed phases (phases with two non-zero condensates). This prediction of dualities is of great use, because for sure it can be shown by the direct calculations but it would be enormously more complicated and time-consuming numerically.


title: Resilient Distributed Diffusion in Networks with Adversaries

abstract: In this paper, we study resilient distributed diffusion for multi-task estimation in the presence of adversaries where networked agents must estimate distinct but correlated states of interest by processing streaming data. We show that in general diffusion strategies are not resilient to malicious agents that do not adhere to the diffusion-based information processing rules. In particular, by exploiting the adaptive weights used for diffusing information, we develop time-dependent attack models that drive normal agents to converge to states selected by the attacker. We show that an attacker that has complete knowledge of the system can always drive its targeted agents to its desired estimates. Moreover, an attacker that does not have complete knowledge of the system including streaming data of targeted agents or the parameters they use in diffusion algorithms, can still be successful in deploying an attack by approximating the needed information. The attack models can be used for both stationary and non-stationary state estimation.In addition, we present and analyze a resilient distributed diffusion algorithm that is resilient to any data falsification attack in which the number of compromised agents in the local neighborhood of a normal agent is bounded. The proposed algorithm guarantees that all normal agents converge to their true target states if appropriate parameters are selected. We also analyze trade-off between the resilience of distributed diffusion and its performance in terms of steady-state mean-square-deviation (MSD) from the correct estimates. Finally, we evaluate the proposed attack models and resilient distributed diffusion algorithm using stationary and non-stationary multi-target localization.


title: Improving Yor\`ub\'a Diacritic Restoration

abstract: Yor\`ub\'a is a widely spoken West African language with a writing system rich in orthographic and tonal diacritics. They provide morphological information, are crucial for lexical disambiguation, pronunciation and are vital for any computational Speech or Natural Language Processing tasks. However diacritic marks are commonly excluded from electronic texts due to limited device and application support as well as general education on proper usage. We report on recent efforts at dataset cultivation. By aggregating and improving disparate texts from the web and various personal libraries, we were able to significantly grow our clean Yor\`ub\'a dataset from a majority Bibilical text corpora with three sources to millions of tokens from over a dozen sources. We evaluate updated diacritic restoration models on a new, general purpose, public-domain Yor\`ub\'a evaluation dataset of modern journalistic news text, selected to be multi-purpose and reflecting contemporary usage. All pre-trained models, datasets and source-code have been released as an open-source project to advance efforts on Yor\`ub\'a language technology.


title: A K-Nearest Neighbor Heuristic for Real-Time DC Optimal Transmission Switching

abstract: While transmission switching is known to reduce power generation costs, the difficulty of solving even DC optimal transmission switching (DCOTS) has prevented optimal transmission switching from becoming commonplace in real-time power systems operation. In this paper, we present a k-nearest neighbors (KNN) heuristic for DCOTS which relies on the insight that, for routine operations on a fixed network, the DCOTS solutions for similar load profiles and generation cost profiles will likely turn off similar sets of lines. We take a data-driven approach and assume that we have DCOTS solutions for many historical instances, which is realistic given that the problem is solved every 5 minutes in practice. Given a new instance, we find a set of "close" instances from the past and return the best of their solutions for the new instance. We present a case study on 7 test networks ranging in size from 118 to 3,375 buses. We compare the proposed heuristic to DCOTS heuristics from the literature, commercial solver heuristics, and a simple greedy local search algorithm. In most cases, we find better quality solutions in less computational time. In addition, the computational time is within the limits imposed by real-time operations, even on larger networks. Last, we present an empirical study of our training data to understand why the heuristic works well.


title: Broad Area Search and Detection of Surface-to-Air Missile Sites Using Spatial Fusion of Component Object Detections from Deep Neural Networks

abstract: Here we demonstrate how Deep Neural Network (DNN) detections of multiple constitutive or component objects that are part of a larger, more complex, and encompassing feature can be spatially fused to improve the search, detection, and retrieval (ranking) of the larger complex feature. First, scores computed from a spatial clustering algorithm are normalized to a reference space so that they are independent of image resolution and DNN input chip size. Then, multi-scale DNN detections from various component objects are fused to improve the detection and retrieval of DNN detections of a larger complex feature. We demonstrate the utility of this approach for broad area search and detection of Surface-to-Air Missile (SAM) sites that have a very low occurrence rate (only 16 sites) over a ~90,000 km^2 study area in SE China. The results demonstrate that spatial fusion of multi-scale component-object DNN detections can reduce the detection error rate of SAM Sites by $>$85% while still maintaining a 100% recall. The novel spatial fusion approach demonstrated here can be easily extended to a wide variety of other challenging object search and detection problems in large-scale remote sensing image datasets.


title: Geometrical Aspects of An Abelian (2,0) Action

abstract: We explore various geometrical aspects of an action for six-dimensional chiral 2-forms based on the formalism of 1903.12196. We elucidate the coupling to general backgrounds and construct the full supersymmetric completion to an abelian (2,0) superconformal lagrangian including matter. We investigate the non-standard diffeomorphism properties of the fields and their relation to the hamiltonian formulation. We also test the action by considering compactifications on a circle, K3 and a Riemann surface. The results are consistent with expectations for an action describing the low-energy physics of an M5-brane in M-theory.


title: Free idempotent generated semigroups: The word problem and structure via gain graphs

abstract: Building on the previous extensive study of Yang, Gould and the present author, we provide a more precise insight into the group-theoretical ramifications of the word problem for free idempotent generated semigroups over finite biordered sets. We prove that such word problems are in fact equivalent to the problem of computing intersections of cosets of certain subgroups of direct products of maximal subgroups of the free idempotent generated semigroup in question, thus providing decidability of those word problems under group-theoretical assumptions related to the Howson property and the coset intersection property. We also provide a basic sketch of the global semigroup-theoretical structure of an arbitrary free idempotent generated semigroup, including the characterisation of Green's relations and the key parameters of non-regular $\mathscr{D}$-classes. In particular, we prove that all Sch\"utzenberger groups of $\mathsf{IG}(\mathcal{E})$ for a finite biordered set $\mathcal{E}$ must be among the divisors of the maximal subgroups of $\mathsf{IG}(\mathcal{E})$.


title: Universal conductivity dependence of pure water polyelectrolyte solutions

abstract: In order to understand how live matter functions one needs to understand the interaction between polyelectrolytes. We discover a general dependence of polyelectrolyte conductivity valid in at least nine decades of polyelectrolyte concentration spanning dilute and semidilute pure water solutions. Furthermore, we showed that current state of the art theories can not explain polyelectrolyte conductivity and suggest the path in transport theories which needs to be taken in order to explain polyelectrolyte conductivity.


title: On the complexity of Broadcast Domination and Multipacking in digraphs

abstract: We study the complexity of the two dual covering and packing distance-based problems Broadcast Domination and Multipacking in digraphs. A dominating broadcast of a digraph $D$ is a function $f:V(D)\to\mathbb{N}$ such that for each vertex $v$ of $D$, there exists a vertex $t$ with $f(t)>0$ having a directed path to $v$ of length at most $f(t)$. The cost of $f$ is the sum of $f(v)$ over all vertices $v$. A multipacking is a set $S$ of vertices of $D$ such that for each vertex $v$ of $D$ and for every integer $d$, there are at most $d$ vertices from $S$ within directed distance at most $d$ from $v$. The maximum size of a multipacking of $D$ is a lower bound to the minimum cost of a dominating broadcast of $D$. Let Broadcast Domination denote the problem of deciding whether a given digraph $D$ has a dominating broadcast of cost at most $k$, and Multipacking the problem of deciding whether $D$ has a multipacking of size at least $k$. It is known that Broadcast Domination is polynomial-time solvable for the class of all undirected graphs (that is, symmetric digraphs), while polynomial-time algorithms for Multipacking are known only for a few classes of undirected graphs. We prove that Broadcast Domination and Multipacking are both NP-complete for digraphs, even for planar layered acyclic digraphs of small maximum degree. Moreover, when parameterized by the solution cost/solution size, we show that the problems are W-hard. We also show that Broadcast Domination is FPT on acyclic digraphs, and that it does not admit a polynomial kernel for such inputs, unless the polynomial hierarchy collapses to its third level. In addition, we show that both problems are FPT when parameterized by the solution cost/solution size together with the maximum out-degree, and as well, by the vertex cover number. Finally, we give for both problems polynomial-time algorithms for some subclasses of acyclic digraphs.


title: Demo: Closed-Loop Control over Wireless -- Remotely Balancing an Inverted Pendulum on Wheels

abstract: Achieving closed-loop control over wireless is crucial in realizing the vision of Industry 4.0 and beyond. This demonstration shows the viability of closed-loop control over wireless through a high-performance wireless solution. The closed-loop control problem involves remote balancing of a two-wheeled robot that represents an inverted pendulum on wheels.


title: Characterizing the Roles of Contributors in Open-source Scientific Software Projects

abstract: The development of scientific software is, more than ever, critical to the practice of science, and this is accompanied by a trend towards more open and collaborative efforts. Unfortunately, there has been little investigation into who is driving the evolution of such scientific software or how the collaboration happens. In this paper, we address this problem. We present an extensive analysis of seven open-source scientific software projects in order to develop an empirically-informed model of the development process. This analysis was complemented by a survey of 72 scientific software developers. In the majority of the projects, we found senior research staff (e.g. professors) to be responsible for half or more of commits (an average commit share of 72%) and heavily involved in architectural concerns (seniors were more likely to interact with files related to the build system, project meta-data, and developer documentation). Juniors (e.g.graduate students) also contribute substantially -- in one studied project, juniors made almost 100% of its commits. Still, graduate students had the longest contribution periods among juniors (with 1.72 years of commit activity compared to 0.98 years for postdocs and 4 months for undergraduates). Moreover, we also found that third-party contributors are scarce, contributing for just one day for the project. The results from this study aim to help scientists to better understand their own projects, communities, and the contributors' behavior, while paving the road for future software engineering research


title: Small-Scale Dynamical Coherence Accompanied with Galaxy Conformity

abstract: The discovery of the coherence between galaxy rotation and neighbor motion in 1-Mpc scales has been reported recently. Following up the discovery, we investigate whether the neighbors in such dynamical coherence also present galaxy conformity, using the Calar Alto Legacy Integral Field Area Survey (CALIFA) data and the NASA-Sloan Atlas (NSA) catalog. We measure the correlation coefficient of g - r colors between the CALIFA galaxies and their neighbors, as a quantitative indicator of galaxy conformity. The neighbors are divided into coherently moving and anti-coherently moving ones, the correlation coefficients from which are compared with each other, in various bins of relative luminosity and projected distance. In most cases, the CALIFA galaxies and their neighbors show positive correlation coefficients in g - r color, even for the anti-coherent neighbors. However, we find statistically significant (2.6-sigma) difference between coherent and anti-coherent neighbors, when the neighbor galaxies are bright (Delta Mr <= -1.0) and close (D <= 400 kpc). That is, when they are bright and close to the CALIFA galaxies, the coherently moving neighbors show stronger conformity with the CALIFA galaxies than the anti-coherently moving neighbors. This result supports that the small-scale dynamical coherence may originate from galaxy interactions as galaxy conformity is supposed to do, which agrees with the conclusion of the previous study.