title: Nodeless Superconductivity and its Evolution with Pressure in the Layered Dirac Semimetal 2M-WS_2

abstract: The interaction between superconductivity and band topology can lead to various unconventional superconducting (SC) states, and represents a new frontier in condensed matter physics research. Recently, the transition metal dichalcogenide (TMD) system 2M-WS2 has been identified as a Dirac semimetal exhibiting both superconductivity with the highest Tc = 8.5 K among all the TMD materials and topological surface states with a single Dirac cone. Here we report on muon spin rotation (muSR) and density functional theory studies of microscopic SC properties and the electronic structure in 2M-WS2 at ambient and under hydrostatic pressures (p_max = 1.9 GPa). The SC order parameter in 2M-WS2 is determined to have single-gap s-wave symmetry. We further show a strong negative pressure effect on Tc and on the SC gap. This may be partly caused by the pressure induced reduction of the size of the electron pocket around the Gamma-point, at which a band inversion appears up to the highest applied pressure. We also find that the superfluid density ns is very weakly affected by pressure. The absence of a strong pressure effect on the superfluid density and the absence of a correlation between ns and Tc in 2M-WS2, in contrast to the other SC TMDs Td-MoTe2 and 2H-NbSe2, is explained in terms of its location in the optimal (ambient pressure) and above the optimal (under pressure) superconducting regions of the phase diagram and its large distance to the other possible competing or cooperating orders. These results hint towards a complex nature of the superconductivity in TMDs, despite the observed s-wave order parameter.


title: Cyclically covering subspaces in $\mathbb{F}_2^n$

abstract: A subspace of $\mathbb{F}_2^n$ is called cyclically covering if every vector in $\mathbb{F}_2^n$ has a cyclic shift which is inside the subspace. Let $h_2(n)$ denote the largest possible codimension of a cyclically covering subspace of $\mathbb{F}_2^n$. We show that $h_2(p)= 2$ for every prime $p$ such that 2 is a primitive root modulo $p$, which, assuming Artin's conjecture, answers a question of Peter Cameron from 1991. We also prove various bounds on $h_2(ab)$ depending on $h_2(a)$ and $h_2(b)$ and extend some of our results to a more general set-up proposed by Cameron, Ellis and Raynaud.


title: On Using Retrained and Incremental Machine Learning for Modeling Performance of Adaptable Software: An Empirical Comparison

abstract: Given the ever-increasing complexity of adaptable software systems and their commonly hidden internal information (e.g., software runs in the public cloud), machine learning based performance modeling has gained momentum for evaluating, understanding and predicting software performance, which facilitates better informed self-adaptations. As performance data accumulates during the run of the software, updating the performance models becomes necessary. To this end, there are two conventional modeling methods: the retrained modeling that always discard the old model and retrain a new one using all available data; or the incremental modeling that retains the existing model and tunes it using one newly arrival data sample. Generally, literature on machine learning based performance modeling for adaptable software chooses either of those methods according to a general belief, but they provide insufficient evidences or references to justify their choice. This paper is the first to report on a comprehensive empirical study that examines both modeling methods under distinct domains of adaptable software, 5 performance indicators, 8 learning algorithms and settings, covering a total of 1,360 different conditions. Our findings challenge the general belief, which is shown to be only partially correct, and reveal some of the important, statistically significant factors that are often overlooked in existing work, providing evidence-based insights on the choice.


title: The Role of Interfacial Inherent Structures in the Fast Crystal Growth from Molten Salts and Metals

abstract: Molecular dynamics simulations of the temperature dependent crystal growth rates of the salts, NaCl and ZnS, from their melts are reported, along with those of a number of pure metals. The growth rate of NaCl and the FCC-forming metals show little evidence of activated control, while that of ZnS and Fe, a BCC forming metal, exhibit activation barriers similar to those observed for diffusion in the melt. Unlike ZnS and Fe, the interfacial inherent structures of NaCl and Cu and Ag are found to be crystalline. We calculate the median displacement between the interfacial liquid and crystalline states and show that this distance is smaller than the cage length, demonstrating that crystal growth in the fast crystallizers can occur via local vibrations and so largely avoid the activated kinetics associated with the larger displacements associated with particle transport.


title: Probing Unseen Planet Populations with Resolved Debris Disk Structures

abstract: Thousands of exoplanets have been found with many widely different from the ones in our own system. Despite the success, systems with planets in wide orbits analogous to those of Jupiter and Saturn, in the critical first several hundred million years of evolution, are virtually unexplored. Where are the low-mass planets that are hidden from our exoplanet detection techniques? Is our Solar System's planetary architecture unique? High-fidelity debris disk images offer an effective method to answer these questions. We can use them to study the formation and evolution of low-mass planets from youth to the age of the Solar System, providing snapshots of the complex processes and valuable insights into the formation and migration history of giant planets at wide orbits. This white paper focuses on resolving debris structures in thermal emission that is applicable to a large unbiased sample. We summarize the properties of the known debris disks and assess the feasibility of resolving them within our current and future infrared and millimeter facilities by adopting uniform criteria. JWST and the 9-m Origins Space Telescope are the most promising missions in the coming decades to resolve almost half of the known disks at high fidelity. Resolved debris structures at multiple wavelengths and at all stages of evolution would reveal the properties of unseen planet populations, enabling a unique demographic study of overall planet formation and evolution.


title: The dynamical evolution of molecular clouds near the Galactic Centre -- III. Tidally--induced star formation in protocluster clouds

abstract: As part of a series of papers aimed at understanding the evolution of the Milky Way's Central Molecular Zone (CMZ), we present hydrodynamical simulations of turbulent molecular clouds orbiting in an accurate model of the gravitational potential extant there. We consider two sets of model clouds differing in the energy content of their velocity fields. In the first, self--virialised set, the turbulent kinetic energies are chosen to be close in magnitude to the clouds' self--gravitational potential energies. Comparison with isolated clouds evolving without an external potential shows that the self--virialised clouds are unable to withstand the compressive tidal field of the CMZ and rapidly collapse, forming stars much faster and reaching gas exhaustion after a small fraction of a Galactocentric orbit. In the second, tidally--virialised, set of simulations, the clouds' turbulent kinetic energies are in equilibrium with the external tidal field. These models are better supported against the field and the stronger turbulence suppresses star formation. Our results strongly support the inference that anomalously low star formation rates in the CMZ are due primarily to high velocity dispersions in the molecular gas. The clouds follow open, eccentric orbits oscillating in all three spatial coordinates. We examine the consequences of the orbital dynamics, particularly pericentre passage, by performing companion simulations of clouds on circular orbits. The increased tidal forces at pericentre produce transient accelerations in star formation rates of at most a factor of 2.7. Our results demonstrate that modelling star formation in galactic centres requires the inclusion of tidal forces.


title: Faster Random $k$-CNF Satisfiability

abstract: We describe an algorithm to solve the problem of Boolean CNF-Satisfiability when the input formula is chosen randomly. We build upon the algorithms of Sch{\"{o}}ning 1999 and Dantsin et al.~in 2002. The Sch{\"{o}}ning algorithm works by trying many possible random assignments, and for each one searching systematically in the neighborhood of that assignment for a satisfying solution. Previous algorithms for this problem run in time $O(2^{n (1- \Omega(1)/k)})$. Our improvement is simple: we count how many clauses are satisfied by each randomly sampled assignment, and only search in the neighborhoods of assignments with abnormally many satisfied clauses. We show that assignments like these are significantly more likely to be near a satisfying assignment. This improvement saves a factor of $2^{n \Omega(\lg^2 k)/k}$, resulting in an overall runtime of $O(2^{n (1- \Omega(\lg^2 k)/k)})$ for random $k$-SAT.


title: Lecture notes on quantitative unique continuation for solutions of second order elliptic equations

abstract: In these lectures we present some useful techniques to study quantitative properties of solutions of elliptic PDEs. Our aim is to outline a proof of a recent result on propagation of smallness. The ideas are also useful in the study of the zero sets of eigenfunctions of Laplace-Beltrami operator and we discuss the connection. Some basic facts about second order elliptic PDEs in divergent form are collected in the Appendix at the end of the notes.


title: An Optimal Graph-Search Method for Secure State Estimation

abstract: The growing complexity of modern Cyber-Physical Systems (CPS) and the frequent communication between their components make them vulnerable to malicious attacks. As a result, secure state estimation is a critical requirement for the control of these systems. Many existing secure state estimation methods suffer from combinatorial complexity which grows with the number of states and sensors in the system. This complexity can be mitigated using optimization-based methods that relax the original state estimation problem, although at the cost of optimality as these methods often identify attack-free sensors as attacked. In this paper, we propose a new optimal graph-search algorithm to correctly identify malicious attacks and to securely estimate the states even in large-scale CPS modeled as linear time-invariant systems. The graph consists of layers, each one containing two nodes capturing a truth assignment of any given sensor, and directed edges connecting adjacent layers only. Then, our algorithm searches the layers of this graph incrementally, favoring directions at higher layers with more attack-free assignments, while actively managing a repository of nodes to be expanded at later iterations. The proposed search bias and the ability to revisit nodes in the repository and self-correct, allow our graph-search algorithm to reach the optimal assignment faster and tackle larger problems. We show that our algorithm is complete and optimal provided that process and measurement noises do not dominate the attack signal. Moreover, we provide numerical simulations that demonstrate the ability of our algorithm to correctly identify attacked sensors and securely reconstruct the state. Our simulations show that our method outperforms existing algorithms both in terms of optimality and execution time.


title: Data-driven Decision Making with Probabilistic Guarantees (Part 1): A Schematic Overview of Chance-constrained Optimization

abstract: Uncertainties from deepening penetration of renewable energy resources have posed critical challenges to the secure and reliable operations of future electric grids. Among various approaches for decision making in uncertain environments, this paper focuses on chance-constrained optimization, which provides explicit probabilistic guarantees on the feasibility of optimal solutions. Although quite a few methods have been proposed to solve chance-constrained optimization problems, there is a lack of comprehensive review and comparative analysis of the proposed methods. Part I of this two-part paper reviews three categories of existing methods to chance-constrained optimization: (1) scenario approach; (2) sample average approximation; and (3) robust optimization based methods. Data-driven methods, which are not constrained by any particular distributions of the underlying uncertainties, are of particular interest. Part II of this two-part paper provides a literature review on the applications of chance-constrained optimization in power systems. Part II also provides a critical comparison of existing methods based on numerical simulations, which are conducted on standard power system test cases.


title: Large-scale quantum-dynamics with matrix product states

abstract: Dynamical electronic- and vibrational-structure theories have received a growing interest in the last years due to their ability to simulate spectra recorded with ultrafast experimental techniques. The exact time evolution of a molecular system can, in principle, be obtained from the time-dependent version of full configuration interaction. Such an approach is, however, limited to few-atom systems due to the exponential increase of its cost with the system dimension. In the present work, we overcome this unfavorable scaling by employing the time-dependent density matrix renormalization group (TD-DMRG) which parametrizes the time-dependent wavefunction as a matrix product state. The time-dependent Schroedinger equation is then integrated with a sweep-based algorithm, as in standard time-independent DMRG. Unlike other TD-DMRG approaches, the one presented here leads to a set of coupled equations that can be integrated exactly. The resulting theory enables us to study real- and imaginary-time evolutions of Hamiltonians comprising more than 20 degrees of freedom that are challenging for current state-of-the-art quantum dynamics algorithms. We apply our algorithm to the simulation of quantum dynamics of models of increasing complexity, ranging from simple excitonic Hamiltonians to more complex ab-initio vibronic ones.


title: Attitude- and Cruise Control of a VTOL Tiltwing UAV

abstract: This paper presents the mathematical modeling, controller design, and flight-testing of an over-actuated Vertical Take-off and Landing (VTOL) tiltwing Unmanned Aerial Vehicle (UAV). Based on simplified aerodynamics and first-principles, a dynamical model of the UAV is developed which captures key aerodynamic effects including propeller slipstream on the wing and post-stall characteristics of the airfoils. The model-based steady-state flight envelope and the corresponding trim-actuation is analyzed and the overactuation of the UAV solved by optimizing for, e.g., power-optimal trims. The developed control system is composed of two controllers: First, a low-level attitude controller based on dynamic inversion and a daisy-chaining approach to handle allocation of redundant actuators. Secondly, a higher-level cruise controller to track a desired vertical velocity. It is based on a linearization of the system and look-up tables to determine the strong and nonlinear variation of the trims throughout the flight-envelope. We demonstrate the performance of the control-system for all flight phases (hover, transition, cruise) in extensive flight-tests.


title: Surface expression of a wall fountain: application to subglacial discharge plumes

abstract: We use laboratory experiments and theoretical modeling to investigate the surface expression of a subglacial discharge plume, as occurs at many fjords around Greenland. The experiments consider a fountain that is released vertically into a homogeneous fluid, adjacent either to a vertical or a sloping wall, that then spreads horizontally at the free surface before sinking back to the bottom. We present a model that separates the fountain into two separate regions: a vertical fountain and a horizontal, negatively buoyant jet. The model is compared to laboratory experiments that are conducted over a range of volume fluxes, density differences, and ambient fluid depths. It is shown that the non-dimensionalised length, width, and aspect ratio of the surface expression are dependent on the Froude number, calculated at the start of the negatively buoyant jet. The model is applied to observations of the surface expression from a Greenland subglacial discharge plume. In the case where the discharge plume reaches the surface with negative buoyancy the model can be used to estimate the discharge properties at the base of the glacier.


title: Neural Grammatical Error Correction with Finite State Transducers

abstract: Grammatical error correction (GEC) is one of the areas in natural language processing in which purely neural models have not yet superseded more traditional symbolic models. Hybrid systems combining phrase-based statistical machine translation (SMT) and neural sequence models are currently among the most effective approaches to GEC. However, both SMT and neural sequence-to-sequence models require large amounts of annotated data. Language model based GEC (LM-GEC) is a promising alternative which does not rely on annotated training data. We show how to improve LM-GEC by applying modelling techniques based on finite state transducers. We report further gains by rescoring with neural language models. We show that our methods developed for LM-GEC can also be used with SMT systems if annotated training data is available. Our best system outperforms the best published result on the CoNLL-2014 test set, and achieves far better relative improvements over the SMT baselines than previous hybrid systems.


title: The role of physiological complexity changes in resting-state EEG in clinical effectiveness of rTMS and tDCS in treatments of resistant depression

abstract: The present literature about possible mechanisms behind the effectivity of noninvasive electromagnetic stimulation in major depressive disorder (MDD) is not very rich. Despite extensive research in applications for clinical practice, the exact effects are yet not clear. We are comparing our previous results about the complexity changes induced by repetitive Transcranial Magnetic Stimulation (rTMS), and transcranial direct current stimulation (tDCS) which are known to modulate neural dynamics. Also, we are reviewing different biomarkers of complexity changes connected to depression, and how they change with the stimulation. TDCS is low-intensity TES, known to have polarity specific effects (neuromodulatory effects), and rTMS is inducing an electric field in the tissue circumstantially via Faraday's law. Both nonlinear modalities of electromagnetic stimulation may affect the levels of physiological complexity in the brain. We also compare the changes of complexity in electroencephalogram (EEG) and electrocardiogram (ECG), as potential future predictors of therapy outcome.


title: Extreme Debris Disk Variability -- Exploring the Diverse Outcomes of Large Asteroid Impacts During the Era of Terrestrial Planet Formation

abstract: The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-size bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system argue for two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.


title: A numerical method for an inverse source problem for parabolic equations and its application to a coefficient inverse problem

abstract: Two main aims of this paper are to develop a numerical method to solve an inverse source problem for parabolic equations and apply it to solve a nonlinear coefficient inverse problem. The inverse source problem in this paper is the problem to reconstruct a source term from external observations. Our method to solve this inverse source problem consists of two stages. We first establish an equation of the derivative of the solution to the parabolic equation with respect to the time variable. Then, in the second stage, we solve this equation by the quasi-reversibility method. The inverse source problem considered in this paper is the linearization of a nonlinear coefficient inverse problem. Hence, iteratively solving the inverse source problem provides the numerical solution to that coefficient inverse problem. Numerical results for the inverse source problem under consideration and the corresponding nonlinear coefficient inverse problem are presented.


title: Rapidly-exploring Random Trees-based Test Generation for Autonomous Vehicles

abstract: Autonomous vehicles are in an intensive research and development stage, and the organizations developing these systems are targeting to deploy them on public roads in a very near future. One of the expectations from fully-automated vehicles is never to cause an accident. However, an automated vehicle may not be able to avoid all collisions, e.g., the collisions caused by other road occupants. Hence, it is important for the system designers to understand the boundary case scenarios where an autonomous vehicle can no longer avoid a collision. In this paper, an automated test generation approach that utilizes Rapidly-exploring Random Trees is presented. A comparison of the proposed approach with an optimization-guided falsification approach from the literature is provided. Furthermore, a cost function that guides the test generation toward almost-avoidable collisions or near-misses is proposed.


title: Diversifying Reply Suggestions using a Matching-Conditional Variational Autoencoder

abstract: We consider the problem of diversifying automated reply suggestions for a commercial instant-messaging (IM) system (Skype). Our conversation model is a standard matching based information retrieval architecture, which consists of two parallel encoders to project messages and replies into a common feature representation. During inference, we select replies from a fixed response set using nearest neighbors in the feature space. To diversify responses, we formulate the model as a generative latent variable model with Conditional Variational Auto-Encoder (M-CVAE). We propose a constrained-sampling approach to make the variational inference in M-CVAE efficient for our production system. In offline experiments, M-CVAE consistently increased diversity by ~30-40% without significant impact on relevance. This translated to a 5% gain in click-rate in our online production system.


title: More on the extremal number of subdivisions

abstract: Given a graph $H$, the extremal number $\mathrm{ex}(n,H)$ is the largest number of edges in an $H$-free graph on $n$ vertices. We make progress on a number of conjectures about the extremal number of bipartite graphs. First, writing $K'_{s,t}$ for the subdivision of the bipartite graph $K_{s,t}$, we show that $\mathrm{ex}(n, K'_{s,t}) = O(n^{3/2 - \frac{1}{2s}})$. This proves a conjecture of Kang, Kim and Liu and is tight up to the implied constant for $t$ sufficiently large in terms of $s$. Second, for any integers $s, k \geq 1$, we show that $\mathrm{ex}(n, L) = \Theta(n^{1 + \frac{s}{sk+1}})$ for a particular graph $L$ depending on $s$ and $k$, answering another question of Kang, Kim and Liu. This result touches upon an old conjecture of Erd\H{o}s and Simonovits, which asserts that every rational number $r \in (1,2)$ is realisable in the sense that $\mathrm{ex}(n,H) = \Theta(n^r)$ for some appropriate graph $H$, giving infinitely many new realisable exponents and implying that $1 + 1/k$ is a limit point of realisable exponents for all $k \geq 1$. Writing $H^k$ for the $k$-subdivision of a graph $H$, this result also implies that for any bipartite graph $H$ and any $k$, there exists $\delta > 0$ such that $\mathrm{ex}(n,H^{k-1}) = O(n^{1 + 1/k - \delta})$, partially resolving a question of Conlon and Lee. Third, extending a recent result of Conlon and Lee, we show that any bipartite graph $H$ with maximum degree $r$ on one side which does not contain $C_4$ as a subgraph satisfies $\mathrm{ex}(n, H) = o(n^{2 - 1/r})$.


title: Dark Matter Through the Quark Vector Current Portal

abstract: We consider models of light dark matter coupled to quarks through a vector current interaction. For low energies, these models must be treated through the effective couplings to mesons, which are implemented here through the chiral Lagrangian. We find the rates of dark matter annihilation and decay to the light mesons, and find the expected photon spectrum from the decay of the hadrons. We compare to current and future observations, and show that there is a significant discovery reach for these models.


title: Archives of CrAO spectral observations. Catalogues of objects and images

abstract: The work described in this article is a continuation of the previously initiated research on archival spectral observations carried out in the Crimea. It covers a time interval of about 90 years and contains information about spectroscopy using various facilities: from the wide-angle astrographs with an objective prism to the main CrAO telescope - ZTSh. A brief history of telescopes and their equipment are presented. The article is illustrated with the possibilities of network access to the catalogues of observations taken with various instruments in the interactive Aladin Sky Atlas with the redirection to the original spectrograms. For this aim, the linear coordinates of the scanned negatives were converted into a scale which corresponds to the wavelengths. The possibilities of taking into account the spectral sensitivity of the recorded images by the absolute energy distribution are shown. A feature of this work is the connection of digitized original observations and the results of their independent processing with the data published for objects in the "Izvestiya of the Crimean Astrophysical Observatory" journals.


title: A $SU(5)\times Z_2$ kink solution and its local stability

abstract: A non-abelian kink inducing asymptotically the breaking pattern $SU(5)\times Z_2\rightarrow SU(4)\times U(1)/Z_4$ is obtained. We consider a fourth order Higgs potential in a $1+1$ theory where the scalar field is in the adjoint representation of $SU(5)$. The perturbative stability of the kink is also evaluated. A Schr\"odinger-like equation for the excitations along each $SU(5)$ generator is determined, and in none of the cases negative eigenvalues compromising the stability of solution are found. In particular, several bounded scalar states are found, being one of them the translational zero mode of the flat space $SU(5)\times Z_2$ kink.


title: Federated Learning Of Out-Of-Vocabulary Words

abstract: We demonstrate that a character-level recurrent neural network is able to learn out-of-vocabulary (OOV) words under federated learning settings, for the purpose of expanding the vocabulary of a virtual keyboard for smartphones without exporting sensitive text to servers. High-frequency words can be sampled from the trained generative model by drawing from the joint posterior directly. We study the feasibility of the approach in two settings: (1) using simulated federated learning on a publicly available non-IID per-user dataset from a popular social networking website, (2) using federated learning on data hosted on user mobile devices. The model achieves good recall and precision compared to ground-truth OOV words in setting (1). With (2) we demonstrate the practicality of this approach by showing that we can learn meaningful OOV words with good character-level prediction accuracy and cross entropy loss.


title: On Fundamental Bounds of Failure Identifiability by Boolean Network Tomography

abstract: Boolean network tomography is a powerful tool to infer the state (working/failed) of individual nodes from path-level measurements obtained by egde-nodes. We consider the problem of optimizing the capability of identifying network failures through the design of monitoring schemes. Finding an optimal solution is NP-hard and a large body of work has been devoted to heuristic approaches providing lower bounds. Unlike previous works, we provide upper bounds on the maximum number of identifiable nodes, given the number of monitoring paths and different constraints on the network topology, the routing scheme, and the maximum path length. These upper bounds represent a fundamental limit on identifiability of failures via Boolean network tomography. Our analysis provides insights on how to design topologies and related monitoring schemes to achieve the maximum identifiability under various network settings. Through analysis and experiments we demonstrate the tightness of the bounds and efficacy of the design insights for engineered as well as real networks.


title: A Tutorial on Sim-ATAV: Simulation-based Adversarial Testing Framework for Autonomous Vehicles

abstract: Testing autonomous vehicles in simulation environments is crucial. Sim-ATAV is an open-source framework developed for experimenting with different test generation techniques in simulation environments for research purposes. This document provides a tutorial on Sim-ATAV with a running example.


title: Implications on the first observation of charm CPV at LHCb

abstract: Very recently, the LHCb Collaboration observed the $CP$ violation (CPV) in the charm sector for the first time, with $\Delta A_{CP}^{\rm dir}\equiv A_{CP}(D^0\to K^+K^-)-A_{CP}(D^0\to \pi^+\pi^-)=(-1.54\pm0.29)\times10^{-3}$. This result is consistent with our prediction of $\Delta A_{CP}^{\rm SM}=(-0.57\sim -1.87)\times 10^{-3}$ obtained in the factorization-assisted topological-amplitude (FAT) approach in [PRD86,036012(2012)]. It implies that the current understanding of the penguin dynamics in charm decays in the Standard Model is reasonable. Motivated by the success of the FAT approach, we further suggest to measure the $D^+\to K^+K^-\pi^+$ decay, which is the next potential mode to reveal the CPV of the same order as $10^{-3}$.


title: Asymmetric quantum shot noise in magnon transport

abstract: We study a frequency-dependent noise-to-current ratio for asymmetric, symmetric, and noncommutative current-noise in a ferromagnetic insulating junction, and extract quantum-mechanical properties of magnon transport at low temperatures. We demonstrate that the noncommutative noise, vanished in the dc-limit (i.e., a classical regime), increases monotonically as a function of frequency, and show that the noncommutative noise associated directly with quantum fluctuations of magnon currents breaks through the classical upper limit determined by the symmetric noise and realizes asymmetric quantum shot noise. Finally, we show that our theoretical predictions are within experimental reach with current device and measurement scheme while exploiting low temperatures. Our work provides a platform toward experimental access to quantum fluctuations of magnon currents.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.quant-gas cond-mat.stat-mech


title: Nonlinear eigenvalue problems for generalized Painlev\'e equations

abstract: Eigenvalue problems for linear differential equations, such as time-independent Schr\"odinger equations, can be generalized to eigenvalue problems for nonlinear differential equations. In the nonlinear context a separatrix plays the role of an eigenfunction and the initial conditions that give rise to the separatrix play the role of eigenvalues. Previously studied examples of nonlinear differential equations that possess discrete eigenvalue spectra are the first-order equation $y'(x)=\cos[\pi xy(x)]$ and the first, second, and fourth Painlev\'e transcendents. It is shown here that the differential equations for the first and second Painlev\'e transcendents can be generalized to large classes of nonlinear differential equations, all of which have discrete eigenvalue spectra. The large-eigenvalue behavior is studied in detail, both analytically and numerically, and remarkable new features, such as hyperfine splitting of eigenvalues, are described quantitatively.


title: INFER: INtermediate representations for FuturE pRediction

abstract: In urban driving scenarios, forecasting future trajectories of surrounding vehicles is of paramount importance. While several approaches for the problem have been proposed, the best-performing ones tend to require extremely detailed input representations (eg. image sequences). But, such methods do not generalize to datasets they have not been trained on. We propose intermediate representations that are particularly well-suited for future prediction. As opposed to using texture (color) information, we rely on semantics and train an autoregressive model to accurately predict future trajectories of traffic participants (vehicles) (see fig. above). We demonstrate that using semantics provides a significant boost over techniques that operate over raw pixel intensities/disparities. Uncharacteristic of state-of-the-art approaches, our representations and models generalize to completely different datasets, collected across several cities, and also across countries where people drive on opposite sides of the road (left-handed vs right-handed driving). Additionally, we demonstrate an application of our approach in multi-object tracking (data association). To foster further research in transferrable representations and ensure reproducibility, we release all our code and data.


title: Spacetime instability due to quantum gravity

abstract: We show that quantum gravity yields exponentially growing gravitational waves. Without a mechanism to stop these modes from growing, the universe would go through a gravitational collapse. For Minkowski background, we propose a solution by choosing an integration contour in Fourier space that does not enclose the problematic modes, thus preventing them from showing up in the effective theory. It turns out that this is only possible when the modes are removed altogether. For an expanding universe, we argue that the runaway modes can be managed accordingly to the dynamics of the Hubble constant, leading to important implications for astrophysics.


title: Study of Activity-Aware Multiple Feedback Successive Interference Cancellation for Massive Machine-Type Communications

abstract: In this work, we propose an activity-aware low-complexity multiple feedback successive interference cancellation (AA-MF-SIC) strategy for massive machine-type communications. The computational complexity of the proposed AA-MF-SIC is as low as the conventional SIC algorithm with very low additional complexity added. The simulation results show that the algorithm significantly outperforms the conventional SIC schemes and other proposals.


title: Detecting and Receiving Phase Modulated Signals with a Rydberg Atom-Based Mixer

abstract: Recently, we introduced a Rydberg-atom based mixer capable of detecting and measuring the phase of a radio-frequency field through the electromagnetically induced transparency (EIT) and Autler-Townes (AT) effect. The ability to measure phase with this mixer allows for an atom-based receiver to detect digital modulated communication signals. In this paper, we demonstrate detection and reception of digital modulated signals based on various phase-shift keying approaches. We demonstrate Rydberg atom-based digital reception of binary phase-shift keying (BPSK), quadrature phase-shift keying (QPSK), and quadrature amplitude (QAM) modulated signals over a 19.626~GHz carrier to transmit and receive a bit stream in cesium vapor. We present measured values of Error Vector Magnitude (EVM, a common communication metric used to assess how accurate a symbol or bit stream is received) as a function of symbol rate for BPSK, QPSK, 16QAM, 32QAM, and 64QAM modulation schemes. These results allow us to discuss the bandwidth of a Rydberg-atom based receiver system.


title: An Alarm System For Segmentation Algorithm Based On Shape Model

abstract: It is usually hard for a learning system to predict correctly on rare events that never occur in the training data, and there is no exception for segmentation algorithms. Meanwhile, manual inspection of each case to locate the failures becomes infeasible due to the trend of large data scale and limited human resource. Therefore, we build an alarm system that will set off alerts when the segmentation result is possibly unsatisfactory, assuming no corresponding ground truth mask is provided. One plausible solution is to project the segmentation results into a low dimensional feature space; then learn classifiers/regressors to predict their qualities. Motivated by this, in this paper, we learn a feature space using the shape information which is a strong prior shared among different datasets and robust to the appearance variation of input data.The shape feature is captured using a Variational Auto-Encoder (VAE) network that trained with only the ground truth masks. During testing, the segmentation results with bad shapes shall not fit the shape prior well, resulting in large loss values. Thus, the VAE is able to evaluate the quality of segmentation result on unseen data, without using ground truth. Finally, we learn a regressor in the one-dimensional feature space to predict the qualities of segmentation results. Our alarm system is evaluated on several recent state-of-art segmentation algorithms for 3D medical segmentation tasks. Compared with other standard quality assessment methods, our system consistently provides more reliable prediction on the qualities of segmentation results.


title: Increasing Iterate Averaging for Solving Saddle-Point Problems

abstract: Many problems in machine learning and game theory can be formulated as saddle-point problems, for which various first-order methods have been developed and proven efficient in practice. Under the general convex-concave assumption, most first-order methods only guarantee an ergodic convergence rate, that is, the uniform averages of the iterates converge at a $O(1/T)$ rate in terms of the saddle-point residual. However, numerically, the iterates themselves can often converge much faster than the uniform averages. This observation motivates increasing averaging schemes that put more weight on later iterates, in contrast to the usual uniform averaging. We show that such increasing averaging schemes, applied to various first-order methods, are able to preserve the $O(1/T)$ convergence rate with no additional assumptions or computational overhead. Extensive numerical experiments on zero-sum game solving, market equilibrium computation and image denoising demonstrate the effectiveness of the proposed schemes. In particular, the increasing averages consistently outperform the uniform averages in all test problems by orders of magnitude. When solving matrix and extensive-form games, increasing averages consistently outperform the last iterates as well. For matrix games, a first-order method equipped with increasing averaging outperforms the highly competitive CFR$^+$ algorithm.


title: On the containment problem for fat points ideals

abstract: In this note we show that Harbourne's conjecture is true for symbolic powers of ideals of points, we check that the stable version of this conjecture is valid for ideals of very general points (resp. generic points) in $\mathbb P_{\mathbb K}^N$ (resp. $\mathbb P_{\mathbb K(\underline{z})}^N$). We also show that this conjecture and the Harbourne-Huneke conjecture are true for a class of ideals $I$ defining fat points obtained from line arrangements in $\mathbb P_{\mathbb K}^2$.


title: Plasmon Resonances and Tachyon Ghost Modes in Highly Conducting Sheets

abstract: Plasmon-polariton modes in two-dimensional electron gases have a dual field-matter nature that endows them with unusual properties when electrical conductivity exceeds a certain threshold set by the speed of light. In this regime plasmons display an interesting relation with tachyons, the hypothetical faster-than-light particles. While not directly observable, tachyons directly impact properties of plasmon modes. Namely, in the ``tachyon'' regime, plasmon resonances remain sharp even when the carrier collision rate $\gamma$ exceeds plasmon resonance frequency. Resonances feature a recurrent behavior as $\gamma$ increases, first broadening and then narrowing and acquiring asymmetric non-Lorentzian lineshapes with power-law tails extending into the tachyon continuum $\omega>ck$. This unusual behavior can be linked to the properties of tachyon poles located beneath $\omega>ck$ branch cuts in the complex $\omega$ plane: as $\gamma$ grows, tachyon poles approach the light cone and hybridize with plasmons. Narrow resonances persisting for $\gamma>\omega$, along with the unusual density and temperature dependence of resonance frequencies, provide clear signatures of the tachyon regime.


title: The Panchromatic Circumgalactic Medium

abstract: Galaxies are surrounded by extended atmospheres, which are often called the circumgalactic medium (CGM) and are the least understood part of galactic ecosystems. The CGM serves as a reservoir of both diffuse, metal-poor gas accreted from the intergalactic medium, and metal-rich gas that is either ejected from galaxies by energetic feedback or stripped from infalling satellites. As such, the CGM is empirically multi-phased and complex in dynamics. Significant progress has been made in the past decade or so in observing the cosmic-ray/B-field, as well as various phases of the CGM. But basic questions remain to be answered. First, what are the energy, mass, and metal contents of the CGM? More specifically, how are they spatially distributed and partitioned in the different components? Moreover, how are they linked to properties of host galaxies and their global clustering and intergalactic medium environments? Lastly, what are the origin, state, and life-cycle of the CGM? This question explores the dynamics of the CGM. Here we illustrate how these questions may be addressed with multi-wavelength observations of the CGM.


title: Dynamical Bonding Driving Mixed Valency in a Metal Boride

abstract: Samarium hexaboride is an anomaly, having many exotic and seemingly mutually incompatible properties. It was proposed to be a mixed-valent semiconductor, and later - a topological Kondo insulator, and yet has a Fermi surface despite being an insulator. We propose a new and unified understanding of SmB$_6$ centered on the hitherto unrecognized dynamical bonding effect: the coexistence of two Sm-B bonding modes within SmB$_6$, corresponding to different oxidation states of the Sm. The mixed valency arises in SmB$_6$ from thermal population of these distinct minima enabled by motion of B. Our model simultaneously explains the thermal valence fluctuations, appearance of magnetic Fermi surface, excess entropy at low temperatures, pressure-induced phase transitions, and related features in Raman spectra and their unexpected dependence on temperature and boron isotope.


title: Restricting Control Flow During Speculative Execution with Venkman

abstract: Side-channel attacks such as Spectre that utilize speculative execution to steal application secrets pose a significant threat to modern computing systems. While program transformations can mitigate some Spectre attacks, more advanced attacks can divert control flow speculatively to bypass these protective instructions, rendering existing defenses useless. In this paper, we present Venkman: a system that employs program transformation to completely thwart Spectre attacks that poison entries in the Branch Target Buffer (BTB) and the Return Stack Buffer (RSB). Venkman transforms code so that all valid targets of a control-flow transfer have an identical alignment in the virtual address space; it further transforms all branches to ensure that all entries added to the BTB and RSB are properly aligned. By transforming all code this way, Venkman ensures that, in any program wanting Spectre defenses, all control-flow transfers, including speculative ones, do not skip over protective instructions Venkman adds to the code segment to mitigate Spectre attacks. Unlike existing defenses, Venkman does not reduce sharing of the BTB and RSB and does not flush these structures, allowing safe sharing and reuse among programs while maintaining strong protection against Spectre attacks. We built a prototype of Venkman on an IBM POWER8 machine. Our evaluation on the SPEC benchmarks and selected applications shows that Venkman increases execution time to 3.47$\times$ on average and increases code size to 1.94$\times$ on average when it is used to ensure that fences are executed to mitigate Spectre attacks. Our evaluation also shows that Spectre-resistant Software Fault Isolation (SFI) built using Venkman incurs a geometric mean of 2.42$\times$ space overhead and 1.68$\times$ performance overhead.


title: Large-scale structures in high Reynolds number rotating Waleffe flow

abstract: We perform direct numerical simulations of rotating turbulent Waleffe flow, the flow between two parallel plates with a sinusoidal streamwise shear driving force, to study the formation of large-scale structures and the mechanisms for momentum transport. We simulate different cyclonic and anti-cyclonic rotations in the range of dimensionless rotation numbers (inverse Rossby numbers) $R_\Omega$ $\in$ $[-0.16, 2.21]$, and fix the Reynolds number to $Re=3.16\times 10^3$, large enough such that the shear transport is almost entirely due to Reynolds stresses and viscous transport is negligible. We find an optimum rotation in anti-cyclonic regime at $R_\Omega=0.63$, where a given streamwise momentum transport in the wall-normal direction is achieved with minimum mean energy of the streamwise flow. We link this optimal transport to the strength of large scale structures, as was done in plane Couette by Brauckmann \& Eckhardt (J. Fluid Mech., 815, 2017). Furthermore, we explore the large-scale structures and their behaviour under spanwise rotation, and find disorganized large structures at $R_\Omega =0$ but highly organized structures in the anti-cyclonic regime, similar to the rolls in rotating plane Couette and turbulent Taylor Couette flow. We compare the large scale structures of plane Couette flow and Waleffe flow, and observe that the streamwise vorticity is localized inside the cores of the rolls. We show that the rolls take energy from the mean flow at long time-scales, and relate these structures to eigenvalues of the streamfunction.


title: Orbital stability of standing waves for the nonlinear Schr\"odinger equation with attractive delta potential and double power repulsive nonlinearity

abstract: In this paper, a nonlinear Schr\"odinger equation with an attractive (focusing) delta potential and a repulsive (defocusing) double power nonlinearity in one spatial dimension is considered. It is shown, via explicit construction, that both standing wave and equilibrium solutions do exist for certain parameter regimes. In addition, it is proved that both types of wave solutions are orbitally stable under the flow of the equation by minimizing the charge/energy functional.


title: Failure-Scenario Maker for Rule-Based Agent using Multi-agent Adversarial Reinforcement Learning and its Application to Autonomous Driving

abstract: We examine the problem of adversarial reinforcement learning for multi-agent domains including a rule-based agent. Rule-based algorithms are required in safety-critical applications for them to work properly in a wide range of situations. Hence, every effort is made to find failure scenarios during the development phase. However, as the software becomes complicated, finding failure cases becomes difficult. Especially in multi-agent domains, such as autonomous driving environments, it is much harder to find useful failure scenarios that help us improve the algorithm. We propose a method for efficiently finding failure scenarios; this method trains the adversarial agents using multi-agent reinforcement learning such that the tested rule-based agent fails. We demonstrate the effectiveness of our proposed method using a simple environment and autonomous driving simulator.


title: Metrics and compactifications of Teichm\"uller spaces of flat tori

abstract: Using the identification of the symmetric space $\mathrm{SL}(n,\mathbb{R})/\mathrm{SO}(n)$ with the Teichm\"uller space of flat $n$-tori of unit volume, we explore several metrics and compactifications of these spaces, drawing inspiration both from Teichm\"uller theory and symmetric spaces. We define and study analogs of the Thurston, Teichm\"uller, and Weil-Petersson metrics. We show the Teichm\"uller metric is a symmetrization of the Thurston metric, which is a polyhedral Finsler metric, and the Weil-Petersson metric is the Riemannian metric of $\mathrm{SL}(n,\mathbb{R})/\mathrm{SO}(n)$ as a symmetric space. We also construct a Thurston-type compactification using measured foliations on $n$-tori, and show that the horofunction compactification with respect to the Thurston metric is isomorphic to it, as well as to a minimal Satake compactification.


title: On a continuum limit of discrete Schr\"odinger operators on square lattice

abstract: The norm resolvent convergence of discrete Schr\"odinger operators to a continuum Schr\"odinger operator in the continuum limit is proved under relatively weak assumptions. This result implies, in particular, the convergence of the spectrum with respect to the Hausdorff distance.


title: A Probabilistic Bitwise Genetic Algorithm for B-Spline based Image Deformation Estimation

abstract: We propose a novel genetic algorithm to solve the image deformation estimation problem by preserving the genetic diversity. As a classical problem, there is always a trade-off between the complexity of deformation models and the difficulty of parameters search in image deformation. 2D cubic B-spline surface is a highly free-form deformation model and is able to handle complex deformations such as fluid image distortions. However, it is challenging to estimate an apposite global solution. To tackle this problem, we develop a genetic operation named probabilistic bitwise operation (PBO) to replace the crossover and mutation operations, which can preserve the diversity during generation iteration and achieve better coverage ratio of the solution space. Furthermore, a selection strategy named annealing selection is proposed to control the convergence. Qualitative and quantitative results on synthetic data show the effectiveness of our method.


title: Unpaired Image Captioning via Scene Graph Alignments

abstract: Most of current image captioning models heavily rely on paired image-caption datasets. However, getting large scale image-caption paired data is labor-intensive and time-consuming. In this paper, we present a scene graph-based approach for unpaired image captioning. Our framework comprises an image scene graph generator, a sentence scene graph generator, a scene graph encoder, and a sentence decoder. Specifically, we first train the scene graph encoder and the sentence decoder on the text modality. To align the scene graphs between images and sentences, we propose an unsupervised feature alignment method that maps the scene graph features from the image to the sentence modality. Experimental results show that our proposed model can generate quite promising results without using any image-caption training pairs, outperforming existing methods by a wide margin.


title: An Exact Auxiliary Variable Gibbs Sampler for a Class of Diffusions

abstract: Stochastic differential equations (SDEs) or diffusions are continuous-valued continuous-time stochastic processes widely used in the applied and mathematical sciences. Simulating paths from these processes is usually an intractable problem, and typically involves time-discretization approximations. We propose an exact Markov chain Monte Carlo sampling algorithm that involves no such time-discretization error. Our sampler is applicable to the problem of prior simulation from an SDE, posterior simulation conditioned on noisy observations, as well as parameter inference given noisy observations. Our work recasts an existing rejection sampling algorithm for a class of diffusions as a latent variable model, and then derives an auxiliary variable Gibbs sampling algorithm that targets the associated joint distribution. At a high level, the resulting algorithm involves two steps: simulating a random grid of times from an inhomogeneous Poisson process, and updating the SDE trajectory conditioned on this grid. Our work allows the vast literature of Monte Carlo sampling algorithms from the Gaussian process literature to be brought to bear to applications involving diffusions. We study our method on synthetic and real datasets, where we demonstrate superior performance over competing methods.


title: Radio follow-up of a candidate gamma-ray transient in the sky localization area of GW170608

abstract: After the identification of a candidate $\gamma$-ray transient in the error region of the binary black hole (BBH) merger GW150914 by the \textit{Fermi} satellite, the question of whether BBH mergers can be associated to electromagnetic counterparts remains highly debated. Here, we present radio follow-up observations of GW170608, a BBH merger that occurred during the second observing run (O2) of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). Our radio follow up focused on a specific field contained in the GW170608 sky localization area, where a candidate high-energy transient was detected by the \textit{Fermi} Large Area Telescope (LAT). We make use of data collected at 1.4\,GHz with the Karl G. Jansky Very Large Array (VLA), as well as with the VLA Low-band Ionosphere and Transient Experiment (VLITE). Our analysis is sensitive to potential radio afterglows with luminosity densities $L_{\rm 1.4\,GHz}\gtrsim 6\times10^{28}$\,erg\,s$^{-1}$\,Hz$^{-1}$. In the most optimistic theoretical models, $\approx 20\%$ of BBH events occurring in massive hosts could be associated with outflows as radio luminous as this. Although we find no evidence for the presence of a radio counterpart associated with the LAT candidate in the GW170608 error region, our analysis demonstrates the feasibility of future radio follow-up observations of well localized BBHs. Comparing our radio upper-limits with theoretical expectations for the radio afterglows potentially associated with jets launched in BBH mergers, we find that for jets of energy $\approx 10^{49}$\,erg seen on-axis, only jet angles $\theta_{jet}\gtrsim 40^{\circ}$ are compatible with the observations.


title: Semantic Alignment: Finding Semantically Consistent Ground-truth for Facial Landmark Detection

abstract: Recently, deep learning based facial landmark detection has achieved great success. Despite this, we notice that the semantic ambiguity greatly degrades the detection performance. Specifically, the semantic ambiguity means that some landmarks (e.g. those evenly distributed along the face contour) do not have clear and accurate definition, causing inconsistent annotations by annotators. Accordingly, these inconsistent annotations, which are usually provided by public databases, commonly work as the ground-truth to supervise network training, leading to the degraded accuracy. To our knowledge, little research has investigated this problem. In this paper, we propose a novel probabilistic model which introduces a latent variable, i.e. the 'real' ground-truth which is semantically consistent, to optimize. This framework couples two parts (1) training landmark detection CNN and (2) searching the 'real' ground-truth. These two parts are alternatively optimized: the searched 'real' ground-truth supervises the CNN training; and the trained CNN assists the searching of 'real' ground-truth. In addition, to recover the unconfidently predicted landmarks due to occlusion and low quality, we propose a global heatmap correction unit (GHCU) to correct outliers by considering the global face shape as a constraint. Extensive experiments on both image-based (300W and AFLW) and video-based (300-VW) databases demonstrate that our method effectively improves the landmark detection accuracy and achieves the state of the art performance.


title: Definite orders with locally free cancellation

abstract: We enumerate all orders in definite quaternion algebras over number fields with the Hermite property; this includes all orders with the cancellation property for locally free modules.


title: Combination of Multiple Global Descriptors for Image Retrieval

abstract: Recent studies in image retrieval task have shown that ensembling different models and combining multiple global descriptors lead to performance improvement. However, training different models for the ensemble is not only difficult but also inefficient with respect to time and memory. In this paper, we propose a novel framework that exploits multiple global descriptors to get an ensemble effect while it can be trained in an end-to-end manner. The proposed framework is flexible and expandable by the global descriptor, CNN backbone, loss, and dataset. Moreover, we investigate the effectiveness of combining multiple global descriptors with quantitative and qualitative analysis. Our extensive experiments show that the combined descriptor outperforms a single global descriptor, as it can utilize different types of feature properties. In the benchmark evaluation, the proposed framework achieves the state-of-the-art performance on the CARS196, CUB200-2011, In-shop Clothes, and Stanford Online Products on image retrieval tasks. Our model implementations and pretrained models are publicly available.


title: Bottom stresses of static packing of granular chains

abstract: We experimentally measure the static stress at the bottom of a granular chains column with a precise and reproducible method. The relation, between the filling mass and the apparent mass converted from the bottom stress, is investigated on various chain lengths. Our measurements reconfirm that the scaling behavior of the stress saturation curves is in accord with the theoretical expectation of the Janssen model. Additionally, the saturation mass is displayed as a nonmonotonic function of the chain length, where a distinguishing transition of the saturation mass is found at the persistence length of the granular chain. We repeat the measurement with another measuring methodology and a silo with different size, respectively, the position of the peak maintains robust. In order to understand the transition of the saturation mass, the friction coefficient and the volume fraction of granular chains are also measured, from which Janssen parameter can be calculated. Finally, we preliminarily measure the bottom stress for two distinct packing structures of long chains, find the effect of the entanglements on the bottom stress, and argue that the entanglements might be responsible for the transition of the saturation mass.


title: Controlling spin-orbit torque by polarization field in multiferroic BiFeO3 based heterostructures

abstract: In the last few years, some ideas of electric manipulations in ferromagnetic heterostructures have been proposed for developing next generation spintronic devices. Among them, the magnetization switching driven by spin-orbit torque (SOT) is being intensely pursued. Especially, how to control the switching current density, which is expected to enrich device functionalities, has aroused much interest among researchers all over the world. In this paper, a novel method to adjust the switching current is proposed, and the BiFeO3 (BFO) based heterostructures with opposite spontaneous polarizations fields show huge changes in both perpendicular magnetic anisotropy and the SOT-induced magnetization switching. The damping-like torques were estimated by using harmonic Hall voltage measurement, and the variation of effective spin Hall angles for the heterostructures with opposite polarizations was calculated to be 272%. At the end of this paper, we have also demonstrated the possible applications of our structure in memory and reconfigurable logic devices.


title: Searching for Beyond the Standard Model Physics with COHERENT Energy and Timing Data

abstract: We search for beyond the Standard Model (BSM) physics by combining COHERENT energy and timing data. Focusing on light, $\lesssim$ GeV mediators, we find the data favor a $\sim 10-1000$ MeV mediator as compared to the Standard Model (SM) best fit at the $\lesssim 2\sigma$ level. The best-fit coupling range is $g \sim 10^{-5}-10^{-3}$. The timing data provides statistical information on the neutrino flavor distributions that is not attainable from the nuclear recoil energies alone. This result accounts for uncertainty in the effective size of the neutron distribution, and highlights the power of including uncertainties on experimental backgrounds, nuclear structure, and BSM physics.


title: Blur Removal via Blurred-Noisy Image Pair

abstract: Complex blur such as the mixup of space-variant and space-invariant blur, which is hard to model mathematically, widely exists in real images. In this paper, we propose a novel image deblurring method that does not need to estimate blur kernels. We utilize a pair of images that can be easily acquired in low-light situations: (1) a blurred image taken with low shutter speed and low ISO noise; and (2) a noisy image captured with high shutter speed and high ISO noise. Slicing the blurred image into patches, we extend the Gaussian mixture model (GMM) to model the underlying intensity distribution of each patch using the corresponding patches in the noisy image. We compute patch correspondences by analyzing the optical flow between the two images. The Expectation Maximization (EM) algorithm is utilized to estimate the parameters of GMM. To preserve sharp features, we add an additional bilateral term to the objective function in the M-step. We eventually add a detail layer to the deblurred image for refinement. Extensive experiments on both synthetic and real-world data demonstrate that our method outperforms state-of-the-art techniques, in terms of robustness, visual quality, and quantitative metrics.


title: On the Weakly Prime-Additive Numbers with Length 4

abstract: In 1992, Erd$\H{o}$s and Hegyv$\'{a}$ri showed that for any prime p, there exist infinitely many length 3 weakly prime-additive numbers divisible by p. In 2018, Fang and Chen showed that for any positive integer m, there exists infinitely many length 3 weakly prime-additive numbers divisible by m if and only if 8 does not divide m. Under the assumption (*) of existence of a prime in certain arithmetic progression with prescribed primitive root, which is true under the Generalized Riemann Hypothesis (GRH), we show for any positive integer m, there exists infinitely many length 4 weakly prime-additive numbers divisible by m. We also present another related result analogous to the length 3 case shown by Fang and Chen.


title: Understanding Exoplanet Atmospheres with UV Observations II: The Far UV and Atmospheric Escape

abstract: Much of the focus of exoplanet atmosphere analysis in the coming decade will be at infrared wavelengths, with the planned launches of the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). However, without being placed in the context of broader wavelength coverage, especially in the optical and ultraviolet, infrared observations produce an incomplete picture of exoplanet atmospheres. Scattering information encoded in blue optical and near-UV (NUV) observations can help determine whether muted spectral features observed in the infrared are due to a hazy/cloudy atmosphere, or a clear atmosphere with a higher mean molecular weight. Farther into the UV, observations can identify atmospheric escape and mass loss from exoplanet atmospheres, providing a greater understanding of the atmospheric evolution of exoplanets, along with composition information from above the cloud deck. In this white paper we focus on the science case for exoplanet observations in the far-UV (FUV); an accompanying white paper led by Jessie Christiansen will focus on the science case in the near-UV (NUV) and blue/optical.


title: Detecting and Gauging Impact on Wikipedia Page Views

abstract: Understanding how various external campaigns or events affect readership on Wikipedia is important to efforts aimed at improving awareness and access to its content. In this paper, we consider how to build time-series models aimed at predicting page views on Wikipedia with the goal of detecting whether there are significant changes to the existing trends. We test these models on two different events: a video campaign aimed at increasing awareness of Hindi Wikipedia in India and the page preview feature roll-out---a means of accessing Wikipedia content without actually visiting the pages---on English and German Wikipedia. Our models effectively estimate the impact of page preview roll-out, but do not detect a significant change following the video campaign in India. We also discuss the utility of other geographies or language editions for predicting page views from a given area on a given language edition.


title: Reinforcement Learning Based Text Style Transfer without Parallel Training Corpus

abstract: Text style transfer rephrases a text from a source style (e.g., informal) to a target style (e.g., formal) while keeping its original meaning. Despite the success existing works have achieved using a parallel corpus for the two styles, transferring text style has proven significantly more challenging when there is no parallel training corpus. In this paper, we address this challenge by using a reinforcement-learning-based generator-evaluator architecture. Our generator employs an attention-based encoder-decoder to transfer a sentence from the source style to the target style. Our evaluator is an adversarially trained style discriminator with semantic and syntactic constraints that score the generated sentence for style, meaning preservation, and fluency. Experimental results on two different style transfer tasks (sentiment transfer and formality transfer) show that our model outperforms state-of-the-art approaches. Furthermore, we perform a manual evaluation that demonstrates the effectiveness of the proposed method using subjective metrics of generated text quality.


title: Robustness of Neural Networks to Parameter Quantization

abstract: Quantization, a commonly used technique to reduce the memory footprint of a neural network for edge computing, entails reducing the precision of the floating-point representation used for the parameters of the network. The impact of such rounding-off errors on the overall performance of the neural network is estimated using testing, which is not exhaustive and thus cannot be used to guarantee the safety of the model. We present a framework based on Satisfiability Modulo Theory (SMT) solvers to quantify the robustness of neural networks to parameter perturbation. To this end, we introduce notions of local and global robustness that capture the deviation in the confidence of class assignments due to parameter quantization. The robustness notions are then cast as instances of SMT problems and solved automatically using solvers, such as dReal. We demonstrate our framework on two simple Multi-Layer Perceptrons (MLP) that perform binary classification on a two-dimensional input. In addition to quantifying the robustness, we also show that Rectified Linear Unit activation results in higher robustness than linear activations for our MLPs.


title: Probabilistic Dense Reconstruction from a Moving Camera

abstract: This paper presents a probabilistic approach for online dense reconstruction using a single monocular camera moving through the environment. Compared to spatial stereo, depth estimation from motion stereo is challenging due to insufficient parallaxes, visual scale changes, pose errors, etc. We utilize both the spatial and temporal correlations of consecutive depth estimates to increase the robustness and accuracy of monocular depth estimation. An online, recursive, probabilistic scheme to compute depth estimates, with corresponding covariances and inlier probability expectations, is proposed in this work. We integrate the obtained depth hypotheses into dense 3D models in an uncertainty-aware way. We show the effectiveness and efficiency of our proposed approach by comparing it with state-of-the-art methods in the TUM RGB-D SLAM and ICL-NUIM dataset. Online indoor and outdoor experiments are also presented for performance demonstration.


title: Virtual User Pairing Non-Orthogonal Multiple Access in Downlink Coordinated Multipoint Transmission

abstract: In this paper, joint transmission coordinated multipoint (JT-CoMP) is exploited by using virtual user pairing non-orthogonal multiple access (VP-NOMA), termed as JT-CoMP VP-NOMA. The technique combines both VP-NOMA for enhancing ergodic sum capacity (ESC) and JT-CoMP for inter-cell interference mitigation. To show the performance gains, ESC of a three-cell scenario is analyzed as a key performance metric. The analytical and simulation results of JT-CoMP VP-NOMA are compared with orthogonal multiple access (OMA), non-orthogonal multiple access (NOMA), and VP-NOMA. It is shown that the proposed JT-CoMP VP-NOMA outperforms the other schemes in the viewpoint of ESC.


title: Document Similarity for Texts of Varying Lengths via Hidden Topics

abstract: Measuring similarity between texts is an important task for several applications. Available approaches to measure document similarity are inadequate for document pairs that have non-comparable lengths, such as a long document and its summary. This is because of the lexical, contextual and the abstraction gaps between a long document of rich details and its concise summary of abstract information. In this paper, we present a document matching approach to bridge this gap, by comparing the texts in a common space of hidden topics. We evaluate the matching algorithm on two matching tasks and find that it consistently and widely outperforms strong baselines. We also highlight the benefits of incorporating domain knowledge to text matching.


title: SciBERT: A Pretrained Language Model for Scientific Text

abstract: Obtaining large-scale annotated data for NLP tasks in the scientific domain is challenging and expensive. We release SciBERT, a pretrained language model based on BERT (Devlin et al., 2018) to address the lack of high-quality, large-scale labeled scientific data. SciBERT leverages unsupervised pretraining on a large multi-domain corpus of scientific publications to improve performance on downstream scientific NLP tasks. We evaluate on a suite of tasks including sequence tagging, sentence classification and dependency parsing, with datasets from a variety of scientific domains. We demonstrate statistically significant improvements over BERT and achieve new state-of-the-art results on several of these tasks. The code and pretrained models are available at https://github.com/allenai/scibert/.


title: Generalized Convolution and Efficient Language Recognition

abstract: Convolution is a broadly useful operation with applications including signal processing, machine learning, probability, optics, polynomial multiplication, and efficient parsing. Usually, however, this operation is understood and implemented in more specialized forms, hiding commonalities and limiting usefulness. This paper formulates convolution in the common algebraic framework of semirings and semimodules and populates that framework with various representation types. One of those types is the grand abstract template and itself generalizes to the free semimodule monad. Other representations serve varied uses and performance trade-offs, with implementations calculated from simple and regular specifications. Of particular interest is Brzozowski's method for regular expression matching. Uncovering the method's essence frees it from syntactic manipulations, while generalizing from boolean to weighted membership (such as multisets and probability distributions) and from sets to n-ary relations. The classic trie data structure then provides an elegant and efficient alternative to syntax. Pleasantly, polynomial arithmetic requires no additional implementation effort, works correctly with a variety of representations, and handles multivariate polynomials and power series with ease. Image convolution also falls out as a special case.


title: Bayesian optimization for tuning and selecting hybrid-density functionals

abstract: The accuracy of some density functional (DF) models, widely used in material science, depends on empirical or free parameters which are commonly tuned using reference physical properties. The optimal value of the free parameters is regularly found using grid search algorithms, which computational complexity scales with the number of points in the grid. In this report, we illustrate that Bayesian optimization (BO), a sample-efficient machine learning algorithm, can efficiently calibrate different density functional models, e.g., hybrid-exchange-correlation and range-separated density functionals. We present that, BO can optimize the free parameters of hybrid-exchange-correlation functionals, with approximately 55 evaluations of the root-mean-square or mean-absolute error functions of the atomization energies and the bond length of the Gaussian-1 (G1) database. We also illustrate that BO can identify, without any prior information, the most appropriate exchange-correlation functional by navigating through the space of density functional models. We optimize and select the free parameters and the exchange-correlation functional form jointly by also minimizing the root-mean-square error function with respect to the atomization energies of the G1 database using BO.


title: Short-term Load Forecasting at Different Aggregation Levels with Predictability Analysis

abstract: Short-term load forecasting (STLF) is essential for the reliable and economic operation of power systems. Though many STLF methods were proposed over the past decades, most of them focused on loads at high aggregation levels only. Thus, low-aggregation load forecast still requires further research and development. Compared with the substation or city level loads, individual loads are typically more volatile and much more challenging to forecast. To further address this issue, this paper first discusses the characteristics of small-and-medium enterprise (SME) and residential loads at different aggregation levels and quantifies their predictability with approximate entropy. Various STLF techniques, from the conventional linear regression to state-of-the-art deep learning, are implemented for a detailed comparative analysis to verify the forecasting performances as well as the predictability using an Irish smart meter dataset. In addition, the paper also investigates how using data processing improves individual-level residential load forecasting with low predictability. Effectiveness of the discussed method is validated with numerical results.


title: Visco- and plasto-elastic fracture of nano-porous polymer sheets

abstract: We study the dependence of the fracture surface energy on the pulling velocity for nano-porous polypropylene (PP) sheets to find two components: the static and dynamic ones. We show that these terms can be interpreted respectively as plastoelastic and viscoelastic components, as has been shown for soft polyethylene (PE) foams in a previous work. Considering significant differences in the pore size, volume fraction and Young's modulus of the present PP and previous PE sheets, the present results suggest a universal physical mechanism for fracture of porous polymer sheets. The simple physical interpretation emerging from the mechanism could be useful for developing tough polymers. Equivalence of Griffith's energy balance in fracture mechanics to a stress criterion is also discussed and demonstrated using the present experimental data.


title: Dynamic Multi Objective Particle Swarm Optimization based on a New Environment Change Detection Strategy

abstract: The dynamic of real-world optimization problems raises new challenges to the traditional particle swarm optimization (PSO). Responding to these challenges, the dynamic optimization has received considerable attention over the past decade. This paper introduces a new dynamic multi-objective optimization based particle swarm optimization (Dynamic-MOPSO).The main idea of this paper is to solve such dynamic problem based on a new environment change detection strategy using the advantage of the particle swarm optimization. In this way, our approach has been developed not just to obtain the optimal solution, but also to have a capability to detect the environment changes. Thereby, DynamicMOPSO ensures the balance between the exploration and the exploitation in dynamic research space. Our approach is tested through the most popularized dynamic benchmark's functions to evaluate its performance as a good method.


title: Torsional Newton-Cartan gravity and strong gravitational fields

abstract: We review how the large $c$ expansion of General Relativity leads to an effective theory in the form of Twistless Torsional Newton-Cartan gravity. We show how this is a strong field expansion around the static sector of General Relativity and illustrate this through two examples.


title: Mask-ShadowGAN: Learning to Remove Shadows from Unpaired Data

abstract: This paper presents a new method for shadow removal using unpaired data, enabling us to avoid tedious annotations and obtain more diverse training samples. However, directly employing adversarial learning and cycle-consistency constraints is insufficient to learn the underlying relationship between the shadow and shadow-free domains, since the mapping between shadow and shadow-free images is not simply one-to-one. To address the problem, we formulate Mask-ShadowGAN, a new deep framework that automatically learns to produce a shadow mask from the input shadow image and then takes the mask to guide the shadow generation via re-formulated cycle-consistency constraints. Particularly, the framework simultaneously learns to produce shadow masks and learns to remove shadows, to maximize the overall performance. Also, we prepared an unpaired dataset for shadow removal and demonstrated the effectiveness of Mask-ShadowGAN on various experiments, even it was trained on unpaired data.


title: Probabilistic Load Forecasting via Point Forecast Feature Integration

abstract: Short-term load forecasting is a critical element of power systems energy management systems. In recent years, probabilistic load forecasting (PLF) has gained increased attention for its ability to provide uncertainty information that helps to improve the reliability and economics of system operation performances. This paper proposes a two-stage probabilistic load forecasting framework by integrating point forecast as a key probabilistic forecasting feature into PLF. In the first stage, all related features are utilized to train a point forecast model and also obtain the feature importance. In the second stage the forecasting model is trained, taking into consideration point forecast features, as well as selected feature subsets. During the testing period of the forecast model, the final probabilistic load forecast results are leveraged to obtain both point forecasting and probabilistic forecasting. Numerical results obtained from ISO New England demand data demonstrate the effectiveness of the proposed approach in the hour-ahead load forecasting, which uses the gradient boosting regression for the point forecasting and quantile regression neural networks for the probabilistic forecasting.


title: Protocols for Packet Quantum Network Intercommunication

abstract: A quantum network, which involves multiple parties pinging each other with quantum messages, could revolutionize communication, computing and basic sciences. The future internet will be a global system of various packet switching quantum and classical networks and we call it \emph{quantum internet}. To build a quantum internet, unified protocols that support the distribution of quantum messages within it are necessary. Intuitively one would extend classical internet protocols to handle quantum messages. However, classical network mechanisms, especially those related to error control and reliable connection, implicitly assume that information can be duplicated, which is not true in the quantum world due to the no-cloning theorem and monogamy of entanglement. In this paper, we investigate and propose protocols for packet quantum network intercommunication. To handle the packet loss problem in transport, we propose a quantum retransmission protocol based on the recursive use of a quantum secret sharing scheme. Other internet protocols are also discussed. In particular, the creation of logical process-to-process connections is accomplished by a quantum version of the three-way handshake protocol.


title: A two-category of Hamiltonian manifolds, and a (1+1+1) field theory

abstract: We define an extended field theory in dimensions $1+1+1$, that takes the form of a `quasi 2-functor' with values in a strict 2-category $\widehat{\mathcal{H}am}$, defined as the `completion of a partial 2-category' $\mathcal{H}am$, notions which we define. Our construction extends Wehrheim and Woodward's Floer Field theory, and is inspired by Manolescu and Woodward's construction of symplectic instanton homology. It can be seen, in dimensions $1+1$, as a real analog of a construction by Moore and Tachikawa. Our construction is motivated by instanton gauge theory in dimensions 3 and 4: we expect to promote $\widehat{\mathcal{H}am}$ to a (sort of) 3-category via equivariant Lagrangian Floer homology, and extend our quasi 2-functor to dimension 4, via equivariant analogues of Donaldson polynomials.


title: Kaonic Nuclear Clusters --- a New Paradigm of Particle and Nuclear Physics ---

abstract: Lambda*=Lambda(1405) plays an essential role in forming kaonic nuclear clusters (KNC). The simplest KNC, K-pp, has the structure Lambda*p=(K-p){I=0}p, where a real kaon migrates between two nucleons, mediating super-strong Lambda*-p$ attraction. Production data of K-pp have been accumulated by DISTO, J-PARC E27 and J-PARC E15 experiments. In the case of K-K-pp, the attractive covalent bond of Lambda*Lambda* is doubly enhanced compared to the Lambda*-p one. As a consequence, a m-th Lambda* multiplet, (Lambda*)m, with m(m-1)/2 bonds becomes more stable than its corresponding neutron aggregate, (n)m, with m=8 ~ 12. This may suggest the possible existence of stable Lambda* matter. The production of K-K-pp by high-energy pp or heavy-ion collisions is awaited as a doorway to so-far unknown Lambda* matter.


title: Controlled inter-state switching between quantized conductance states in resistive devices for multilevel memory

abstract: A detailed understanding of quantization conductance (QC), their correlation with resistive switching phenomena and controlled manipulation of quantized states is crucial for realizing atomic-scale multilevel memory elements. Here, we demonstrate highly stable and reproducible quantized conductance states (QC-states) in Al/Niobium oxide/Pt resistive switching devices. Three levels of control over the QC-states, required for multilevel quantized state memories, like, switching ON to different quantized states, switching OFF from quantized states, and controlled inter-state switching among one QC states to another has been demonstrated by imposing limiting conditions of stop-voltage and current compliance. The well defined multiple QC-states along with a working principle for switching among various states show promise for implementation of multilevel memory devices.


title: Magnetic competition induced colossal magnetoresistance in n-type HgCr2Se4 under high pressures

abstract: The n-type HgCr2Se4 exhibits a sharp semiconductor-to-metal transition (SMT) in resistivity accompanying the ferromagnetic order at TC = 106 K. Here, we investigate the effects of pressure and magnetic field on the concomitant SMT and ferromagnetic order by measuring resistivity, dc and ac magnetic susceptibility, as well as single-crystal neutron diffraction under various pressures up to 8 GPa and magnetic fields up to 8 T. Our results demonstrate that the ferromagnetic metallic ground state of n-type HgCr2Se4 is destabilized and gradually replaced by an antiferromagnetic, most likely a spiral magnetic, and insulating ground state upon the application of high pressure. On the other hand, the application of external magnetic fields can restore the ferromagnetic metallic state again at high pressures, resulting in a colossal magnetoresistance (CMR) as high as ~ 3 * 10^11 % under 5 T and 2 K at 4 GPa. The present study demonstrates that n-type HgCr2Se4 is located at a peculiar critical point where the balance of competion between ferromagnetic and antiferromagnetic interactions can be easily tipped by the external stimuli, providing a new platform for achieving CMR in a single-valent system.


title: Impact of the chromo-electromagnetic field fluctuations on transport coefficients of heavy quarks and shear viscosity to entropy density ratio of quark-gluon plasma

abstract: The chromo-electromagnetic field fluctuations in the quark-gluon plasma (QGP) play an important role as these field fluctuations result energy gain of heavy quarks. We consider these fluctuations and evaluate the transport coefficients, e.g., drag and diffusion coefficients of charm quarks and shear viscosity to entropy density ratio ($\eta/s$) of the QGP. We find a significant effect of such fluctuations on the transport coefficients. These fluctuations cause a reduction of the drag and diffusion coefficients. We also observe that the shear viscosity to entropy density ratio of the QGP is closer to the value obtained in Lattice QCD (LQCD) and functional renormalization group calculations when the effects of such fluctuations are included.


title: G-networks and the optimization of supply chains

abstract: Supply chains are fundamental to the economy of the world and many supply chains focus on perishable items, such as food, or even clothing that is subject to a limited shelf life due to fashion and seasonable effects. G-networks have not been previously applied to this important area. Thus in this paper, we apply G-networks to supply chain systems and investigate an optimal order allocation problem for a N-node supply chain with perishable products that share the same order source of fresh products. The objective is to find an optimal order allocation strategy to minimize the purchase price per object from the perspective of the customers. An analytical solution based on G-networks with batch removal, together with optimization methods are shown to produce the desired results. The results are illustrated by a numerical example with realistic parameters.


title: Suggestion of the DLV dimensionless number system to represent the scaled behavior of structures under impact loads

abstract: A group of dimensionless numbers, termed DLV (Density-Length-Velocity) system, is put forward to represent the scaled behavior of structures under impact loads. It is obtained by means of the Buckingham Pi theorem with an alternative basis. The distinct features of this group of dimensionless numbers are that it relates physical quantities of the impacted structure with essential basis of the Density, the Length and the Velocity, and thus it can represent the scaled influence of material property, geometry characteristic and velocity on the behavior of structures. The newly 15 proposed dimensionless numbers reflect three advantages. (1) The intuitively clear physical significance of these dimensionless numbers, such as the ratios of force intensity, force, moment of inertia to the corresponding dynamic quantities, the Johnson's damage number Dn and Zhao's response number Rn etc. are naturally included. (2) The property of direct matching the dimensionless expression of response equations of dynamic problems with these dimensionless numbers through simple equation analysis; (3) The ability of addressing non-scaling problems for different materials and strain-rate-sensitive as well as the VSG (initial impact Velocity-dynamic flow Stress-impact mass G) system. Four classical impact models are used to verify the direct matching property and the non-scaling addressing ability of the DLV system by equation analysis. The results show that the proposed dimensionless number system is simple, clear and efficient, and we suggest using it to represent the scaled behavior of structures under impact loads.


title: Decomposing information into copying versus transformation

abstract: In many real-world systems, information can be transmitted in two qualitatively different ways: by copying or by transformation. Copying occurs when messages are transmitted without modification, e.g., when an offspring receives an unaltered copy of a gene from its parent. Transformation occurs when messages are modified systematically during transmission, e.g., when mutational biases occur during genetic replication. Standard information-theoretic measures do not distinguish these two modes of information transfer, although they may reflect different mechanisms and have different functional consequences. Starting from a few simple axioms, we derive a decomposition of mutual information into the information transmitted by copying versus the information transmitted by transformation. We begin with a decomposition that applies when the source and destination of the channel have the same set of messages and a notion of message identity exists. We then generalize our decomposition to other kinds of channels, which can involve different source and destination sets and broader notions of similarity. In addition, we show that copy information can be interpreted as the minimal work needed by a physical copying process, which is relevant for understanding the physics of replication. We use the proposed decomposition to explore a model of amino acid substitution rates. Our results apply to any system in which the fidelity of copying, rather than simple predictability, is of critical relevance.

categories: cs.IT cond-mat.stat-mech math.IT stat.ML


title: Revising the Wilks Scoring System for pro RAW Powerlifting

abstract: Purpose: In powerlifting the total result is highly dependent on the athletes bodyweight. Powerlifting is divided to equipped and RAW types. Pro RAW powerlifting competitions use the Wilks scoring system to compare and rank powerlifting results across bodyweights, to choose the winners. The Wilks formulas are 5th order polynomials fitted to equipped powerlifting results from the years 1987 - 1994. The Wilks formulas were meant to be updated every 2 to 5 years, but have never been revised. This study aims to update the Wilks scoring system for use in pro RAW powerlifting. Methods: The 10 highest RAW powerlifting totals per weight class for bodyweights 60 - 175 kg for men, and 44 kg and up for women, were collected from the openpowerlifting.org database. Polynomials were fitted to these data separately for men and women, using polynomial regression. The scores for men and women were calculated using these fits, and the results are assessed. Results: 4th order polynomials were chosen for best describing the dependence of powerlifting total result on bodyweight. The mens scores were normalized to 500 points, and womens to 455 points, on those fitted curves, to bring the top scores to the same level across weight classes of each sex. These scores appeared similar across the different weight classes for both men and women. This even scoring across body weights is here considered the single most important characteristic of a scoring system for pro competitions. Key words: Wilks score, revised prediction curve, polynomial regression, powerlifting


title: Bayesian Estimation Based Parameter Estimation for Composite Load

abstract: Accurate identification of parameters of load models is essential in power system computations, including simulation, prediction, and stability and reliability analysis. Conventional point estimation based composite load modeling approaches suffer from disturbances and noises and provide limited information of the system dynamics. In this work, a statistic (Bayesian Estimation) based distribution estimation approach is proposed for both static (ZIP) and dynamic (Induction Motor) load modeling. When dealing with multiple parameters, Gibbs sampling method is employed. In each iteration, the proposal samples each parameter while keeps others fixed. The proposed method provides a distribution estimation of load models coefficients and is robust to measurement errors.


title: On the mean square displacement in Levy walks

abstract: Many physical and biological processes are modeled by "particles" undergoing L\'evy random walks. A feature of significant interest in these systems is the mean square displacement (MSD) of the particles. Long-time asymptotic approximations of the MSD have been established, via the Tauberian Theorem, for systems in which the distribution of the step durations is asymptotically a power law of infinite variance. We extend these results, using elementary analysis, and obtain closed-form expressions as well as power law bounds for the MSD in equilibrium, and representations of the MSD as sums of super-linear, linear, and sub-linear terms. We show that the super-linear components are determined by the mean and asymptotics of the step durations, but that the linear and sub-linear components (whose size has implications for the accuracy of the asymptotic approximation) depend on the entire distribution function.


title: On Generalizations of the Newton-Raphson-Simpson Method

abstract: We present generalizations of the Newton-Raphson-Simpson method. Specifically, for a positive integer $m$ and the sequence of coefficients of a Taylor series of a function $f(z)$, we define an algorithm we denote by NRS($m$) which is a way to evaluate, in our terminology, a sum of $m$ formal zeros of $f(z)$. We prove that NRS(1) yields the familiar iterations of the Newton-Raphson-Simpson method. We also prove that NRS($m$) is way to evaluate certain $\mathscr{A}$-hypergeometric series defined by Sturmfels. In order to define these algorithms, we make use of combinatorial objects which we call trees with negative vertex degree.


title: A Galactic centre gravitational-wave Messenger

abstract: Our existence in the Universe resulted from a rare combination of circumstances. The same must hold for any highly developed extraterrestrial civilisation, and if they have ever existed in the Milky Way, they would likely be scattered over large distances in space and time. However, all technologically advanced species must be aware of the unique property of the galactic centre: it hosts Sagittarius A* (Sgr A*), the closest supermassive black hole to anyone in the Galaxy. A civilisation with sufficient technical know-how may have placed material in orbit around Sgr A* for research, energy extraction, and communication purposes. In either case, its orbital motion will necessarily be a source of gravitational waves. We show that a Jupiter-mass probe on the retrograde innermost stable circular orbit around Sgr A* emits, depending on the black hole spin, at a frequency of $f_{GW} = 0.63 - 1.07$ mHz and with a power of $P_{GW}=2.7 \times\, 10^{36} - 2.0 \times\, 10^{37}$ erg/s. We discuss that the energy output of a single star is sufficient to stabilise the location of an orbiting probe for a billion years against gravitational wave induced orbital decay. Placing and sustaining a device near Sgr A* is therefore astrophysically possible. Such a probe will emit an unambiguously artificial continuous gravitational wave signal that is observable with LISA-type detectors.


title: On the Theory of Dynamic Graph Regression Problem

abstract: Most of real-world graphs are dynamic, i.e., they change over time by a sequence of update operations. While the regression problem has been studied for static graphs and temporal graphs, it is not investigated for general dynamic graphs. In this paper, we study regression over dynamic graphs. First, we present the notion of update-efficient matrix embedding, that defines conditions sufficient for a matrix embedding to be effectively used for dynamic graph regression (under l2 norm). Then, we show that given a n*m update-efficient matrix embedding (e.g., the adjacency matrix) and after an update operation in the graph, the exact optimal solution of linear regression can be updated in O(nm) time for the revised graph. Moreover, we show that this also holds when the matrix embedding is the Laplacian matrix and the update operations are restricted to edge insertion/deletion. In the end, by conducting experiments over synthetic and real-world graphs, we show the high efficiency of updating the solution of graph regression.