title: Charge-stripe order, antiferromagnetism, and spin dynamics in the cuprate-analog nickelate La$_4$Ni$_3$O$_8$

abstract: We report a muon spin rotation ($\mu$SR) study of the cuprate-analog nickelate La$_4$Ni$_3$O$_8$, which undergoes a transition at 105 K to a low-temperature phase with charge-stripe and antiferromagnetic (AFM) order on square planar NiO$_2$ layers. Zero-field $\mu$SR shows that this transition is abrupt and that the AFM configuration is commensurate. Comparison of observed muon precession frequencies with Ni dipolar field calculations yields Ni moments $\lesssim 0.5\mu_B$. Dynamic muon spin relaxation above 105 K suggests critical slowing of Ni spin fluctuations, but is inconsistent with corresponding $^{139}$La NMR results. Critical slowing and an abrupt transition are also observed in the planar cuprate AFM La$_2$CuO$_{4+\delta}$, where they are evidence for weakly interplanar-coupled two-dimensional AFM spin fluctuations, but our $\mu$SR data do not agree quantitatively with theoretical predictions for this scenario.


title: Updated extraction of the APOGEE 1.5273 {\mu}m diffuse interstellar band: a Planck view on the carrier depletion in dense cores

abstract: The latest SDSS/APOGEE data release DR14 has provided an increased number of stellar spectra in the H band and associated stellar models using an innovative algorithm known as The Cannon. We took advantage of these novelties to extract the 15 273 {\AA} near-infrared DIB and to study its link with dust extinction and emission. We modified our automated fitting methods dedicated to hot stars and used in earlier studies with some adaptations motivated by the change from early- or intermediate-type stars to red giants. A new method has also been developed to quantify the upper limits on DIB strengths. We compared our DIB measurements with the stellar extinctions Av from the Starhorse database. We then compared the resulting DIB-extinction ratio with the dust optical depth derived from Planck data, globally and also separately for nearby off-Plane cloud complexes. Our analysis has led to the production of a catalog containing 124 064 new measurements of the 15 273 {\AA} DIB, allowing us to revisit the correlation between DIB strength and dust reddening. The new data reveal clearly that the sky-averaged 15 273 {\AA} DIB strength is linearly correlated with Av over two orders as reported by earlier studies but leveling-off with respect to extinction for highly reddened lines-of-sight behind dense clouds. The comparison with Planck individual optical depths reveals in a conspicuous way this DIB depletion in the dense cores and shows it applies to all off-Plane dense clouds. APOGEE measurements confirm the ubiquity of the 15 273 {\AA} DIB carrier decrease with respect to dust grains in dense cloud cores, in a manner that can be empirically related to the dust optical depth reached in the cloud.


title: On some generalizations of the property $B$-problem of an $n$-uniform hypergraph

abstract: The extremal problem of hypergraph colorings related to Erd\H{o}s--Hajnal property $B$-problem is considered. Let $k$ be a natural number. The problem is to find the value of $m_k(n)$ equal to the minimal number of edges in an $n$-uniform hypergraph not admitting $2$-colorings of the vertex set such that every edge of the hypergraph contains at least $k$ vertices of each color. In this paper we obtain new lower bounds for $m_k(n).$


title: One-shot conclusive multiport quantum dense coding capacities

abstract: We introduce a probabilistic version of the one-shot quantum dense coding protocol in both two- and multiport scenarios, and refer to it as conclusive quantum dense coding. Specifically, we analyze the corresponding capacities of two-qubit, two-qutrit, and three-qubit shared states. We identify cases where Pauli and generalized Pauli operators are not sufficient as encoders to attain the optimal one-shot conclusive quantum dense coding capacities. We find that there is a rich connection between the capacities, and the bipartite and multipartite entanglements of the shared state.


title: A Dimensionality-Reduction Strategy to Compute Shortest Paths in Urban Water Networks

abstract: The efficient computation of shortest paths in complex networks is essential to face new challenges related to critical infrastructures such as a near real-time monitoring and control and the management of big size systems. In particular, using information on the minimum paths in water distribution networks (WDNs) allows to track the diffusion of contaminants and to quantify the resilience and criticality of the system. This is, ultimately, approached by considering dynamically changing path-weights that depend on the flow or on other information available at run-time. These analyses tipically include all the WDN assets but reducing the high degree of physical details with a minimum lost of key information for their performance assessment. This paper proposes a strategy to compute minimum paths that is based on a dimensionality-reduction process. Specifically, the network is partitioned into communities and suitably modified to obtain a reduced complexity representation (e.g., in terms of number of nodes and links). The paper shows how this novel, reduced representation is equivalent to the traditional network on computing the shortest paths. The proposed approach is validated considering two utility networks as case studies. The results show that the proposed method provides the exact solution for the shortest path with a computational-time reduction consistently over 50\% and up to 90\% for some cases. Furthermore, the application of the proposal on WDNs partitioning shows both hydraulic and economic advantages thanks to their monitoring and controlling at near real-time.


title: Output Feedback Controller Synthesis for Negative Imaginary Systems

abstract: This paper presents necessary and sufficient conditions for deriving a strictly proper dynamic controller which satisfies the negative imaginary output feedback control problem. Our synthesis method divides the output feedback control problem into a state feedback problem and a dual output injection problem. Thus, a controller is formed using the solutions to a pair of dual algebraic Riccati equations. Finally, an illustrative example is offered to show how this method may be applied.


title: A Multi Hidden Recurrent Neural Network with a Modified Grey Wolf Optimizer

abstract: Identifying university students' weaknesses results in better learning and can function as an early warning system to enable students to improve. However, the satisfaction level of existing systems is not promising. New and dynamic hybrid systems are needed to imitate this mechanism. A hybrid system (a modified Recurrent Neural Network with an adapted Grey Wolf Optimizer) is used to forecast students' outcomes. This proposed system would improve instruction by the faculty and enhance the students' learning experiences. The results show that a modified recurrent neural network with an adapted Grey Wolf Optimizer has the best accuracy when compared with other models.


title: Primary damage production in the presence of extended defects and growth of vacancy-type dislocation loops in hcp zirconium

abstract: Production rates in long-term predictive radiation damage accumulation models are generally considered independent of the material's microstructure for reactor components. In this study, the effect of pre-existing microstructural elements on primary damage production in alpha-Zr -- and vice-versa -- is assessed by molecular dynamics (MD) simulations. a-type dislocation loops, c-component dislocation loops and a tilt grain boundary (GB) were considered. Primary damage production is reduced in the presence of all these microstructural elements, and clustering behavior is dependent on the microstructure. Collision cascades do not cause a-type loop growth or shrinkage, but they cause c-component loop shrinkage. Cascades in the presence of the GBs produce more vacancies than interstitials. This result, as well as other theoretical, MD and experimental evidence, confirm that vacancy loops will grow in the vacancy supersaturated environment near GBs. Distinct temperature-dependent growth regimes are identified. Also, MD reveals cascade-induced events where a-type vacancy loops are absorbed by GBs. Fe segregation at the loops inhibits this cascade-induced absorption.


title: High Performance Monte Carlo Simulation of Ising Model on TPU Clusters

abstract: Large-scale deep learning benefits from an emerging class of AI accelerators. Some of these accelerators' designs are general enough for compute-intensive applications beyond AI and Cloud TPU is one such example. In this paper, we demonstrate a novel approach using TensorFlow on Cloud TPU to simulate the two-dimensional Ising Model. TensorFlow and Cloud TPU framework enable the simple and readable code to express the complicated distributed algorithm without compromising the performance. Our code implementation fits into a small Jupyter Notebook and fully utilizes Cloud TPU's efficient matrix operation and dedicated high speed inter-chip connection. The performance is highly competitive: it outperforms the best published benchmarks to our knowledge by 60% in single-core and 250% in multi-core with good linear scaling. When compared to Tesla V100 GPU, the single-core performance maintains a ~10% gain. We also demonstrate that using low precision arithmetic---bfloat16---does not compromise the correctness of the simulation results.


title: Classification of commutative pairs of surjective maps of interval, one of which is unimodal

abstract: We introduce here a classification of unimodal maps $[0, 1]\rightarrow [0, 1]$, which commute with piecewise linear surjective maps $[0, 1]\rightarrow [0, 1]$. Remind that if continuous piecewise linear unimodal map $g$ commutes with a non-constant piecewise linear map $\psi$, which is not an iteration of $g$, then $g$ is topologically conjugated with the tent map by piecewise linear conjugacy. We use the obtained classification to illustrate the mentioned fact.


title: Ferroelectric ZrO$_{2}$ monolayers as buffer layers between SrTiO$_{3}$ and Si

abstract: A monolayer of ZrO$_{2}$ has recently been grown on the Si(001) surface and shown to have ferroelectric properties, which signifies the realization of the lowest possible thickness in ferroelectric oxides [M. Dogan et al., Nano Lett., 18 (1) (2018)]. In our previous computational study, we reported on the multiple (meta)stable configurations of ZrO$_{2}$ monolayers on Si, and how switching between a pair of differently polarized configurations may explain the observed ferroelectric behavior of these films [M. Dogan and S. Ismail-Beigi, arXiv:1902.01022 (2019)]. In the current study, we conduct a DFT-based investigation of (i) the effect of oxygen content on the ionic polarization of the oxide, and (ii) the role of zirconia monolayers as buffer layers between silicon and a thicker oxide film that is normally paraelectric on silicon, e.g. SrTiO$_{3}$. We find that (i) total energy-vs-polarization behavior of the monolayers, as well as interface chemistry, is highly dependent on the oxygen content; and (ii) SrTiO$_{3}$/ZrO$_{2}$/Si stacks exhibit multiple (meta)stable configurations and polarization profiles, i.e. zirconia monolayers can induce ferroelectricity in oxides such as SrTiO$_{3}$ when used as a buffer layer. This may enable a robust non-volatile device architecture where the thickness of the gate oxide (here strontium titanate) can be chosen according to the desired properties.


title: Coupled reaction-diffusion equations on adjacent domains

abstract: We consider a reaction-diffusion system for two densities lying in adjacent domains of $\mathbb{R}^N$. We treat two configurations: either a cylinder and its complement, or two half-spaces. Diffusion and reaction heterogeneities for the two densities are considered, and an exchange occurs through the separating boundary. We study the long-time behavior of the solution, and, when it converges to a positive steady state, we prove the existence of an asymptotic speed of propagation in some specific directions. Moreover, we determine how such a speed qualitatively depends with respect to several parameters appearing in the model. In the case $N=2$, we compare such properties to those studied in [6-9] for a model with a line representing a road of fast diffusion at the boundary of a half-plane, which can be seen as a singular limit of the problem studied here.


title: Deriving Iodine-free spectra for high-resolution echelle spectrographs

abstract: We describe a new method to derive clean, iodine-free spectra directly from observations acquired using high-resolution echelle spectrographs equipped with iodine cells. The main motivation to obtain iodine-free spectra is to use portions of the spectrum that are superimposed with the dense forest of iodine absorption lines, in order to retrieve lines that can be used to monitor the magnetic activity of the star, helping to validate candidate planets. In short, we provide a straight-forward methodology to clean the spectra by using the forward model used to derive radial velocities, the Line Spread Function information plus the stellar spectrum without iodine to reconstruct and subtract the iodine spectrum from the observations. We show our results using observations of the star $\tau$ Ceti acquired with the PFS, HIRES and UCLES spectrographs, reaching an iodine-free spectrum correction at the $\sim$1% RMS level. We additionally discuss the limitations and further applications of the method.


title: Fairness in Algorithmic Decision Making: An Excursion Through the Lens of Causality

abstract: As virtually all aspects of our lives are increasingly impacted by algorithmic decision making systems, it is incumbent upon us as a society to ensure such systems do not become instruments of unfair discrimination on the basis of gender, race, ethnicity, religion, etc. We consider the problem of determining whether the decisions made by such systems are discriminatory, through the lens of causal models. We introduce two definitions of group fairness grounded in causality: fair on average causal effect (FACE), and fair on average causal effect on the treated (FACT). We use the Rubin-Neyman potential outcomes framework for the analysis of cause-effect relationships to robustly estimate FACE and FACT. We demonstrate the effectiveness of our proposed approach on synthetic data. Our analyses of two real-world data sets, the Adult income data set from the UCI repository (with gender as the protected attribute), and the NYC Stop and Frisk data set (with race as the protected attribute), show that the evidence of discrimination obtained by FACE and FACT, or lack thereof, is often in agreement with the findings from other studies. We further show that FACT, being somewhat more nuanced compared to FACE, can yield findings of discrimination that differ from those obtained using FACE.


title: Performance Analysis and Enhancements for In-Band Full-Duplex Wireless Local Area Networks

abstract: In-Band Full-Duplex (IBFD) is a technique that enables a wireless node to simultaneously transmit a signal and receive another on the same assigned frequency. Thus, IBFD wireless systems can provide up to twice the channel capacity compared to conventional Half-Duplex (HD) systems. In order to study the feasibility of IBFD networks, reliable models are needed to capture anticipated benefits of IBFD above the physical layer (PHY). In this paper, an accurate analytical model based on Discrete-Time Markov Chain (DTMC) analysis for IEEE 802.11 Distributed Coordination Function (DCF) with IBFD capabilities is proposed. The model captures all parameters necessary to calculate important performance metrics which quantify enhancements introduced as a result of IBFD solutions. Additionally, two frame aggregation schemes for Wireless Local Area Networks (WLANs) with IBFD features are proposed to increase the efficiency of data transmission. Matching analytical and simulation results with less than 1% average errors confirm that the proposed frame aggregation schemes further improve the overall throughput by up to 24% and reduce latency by up to 47% in practical IBFD-WLANs. More importantly, the results assert that IBFD transmission can only reduce latency to a suboptimal point in WLANs, but frame aggregation is necessary to minimize it.


title: How Martian araneiforms get their shapes: morphological analysis and diffusion-limited aggregation model for polar surface erosion

abstract: Araneiforms are radially converging systems of branching troughs exhibiting fractal properties. They are found exclusively in the Southern polar regions of Mars and believed to be result of multiple repetitions of cold CO2 gas jets eruptions. Araneiform troughs get carved by the overpressurized gas rushing underneath a seasonal ice layer towards a newly created opening. Current work is an attempt to quantitatively analyze araneiforms patterns and model their formation mechanism. The dendritic quality of most araneiforms are suggestive that they can be described in terms commonly applied to terrestrial rivers. We have adapted and for the first time applied to Martian araneiforms qualitative morphological analysis typically used for terrestrial rivers. We have shown that the large and well-developed araneiforms (with tributary orders larger than 4) closely follow Horton's law of tributary orders and have bifurcation ratio that falls well inside the range of terrestrial rivers. We have implemented a two-dimensional Diffusion-Limited Aggregation (DLA) model that describes formation of dendrite shapes by mathematical probabilistic means. We compared modeled dendrite shapes to the araneiform shapes observed in the Martian polar regions and evaluated their similarity using the morphological analysis of araneiforms. We showed that DLA model can successfully recreate 2D shapes of different observed araneiforms. Modeling the creation of araneiform patterns with DLA leads to better understanding of seasonal processes that create them.


title: Resource Allocation Mechanism for Media Handling Services in Cloud Multimedia Conferencing

abstract: Multimedia conferencing is the conversational exchange of multimedia content between multiple parties. It has a wide range of applications (e.g., Massively Multiplayer Online Games (MMOGs) and distance learning). Media handling services (e.g., video mixing, transcoding, and compressing) are critical to multimedia conferencing. However, efficient resource usage and scalability still remain important challenges. Unfortunately, the cloud-based approaches proposed so far have several deficiencies in terms of efficiency in resource usage and scaling, while meeting Quality of Service (QoS) requirements. This paper proposes a solution which optimizes resource allocation and scales in terms of the number of participants while guaranteeing QoS. Moreover, our solution composes different media handling services to support the participants' demands. We formulate the resource allocation problem mathematically as an Integer Linear Programming (ILP) problem and design a heuristic for it. We evaluate our proposed solution for different numbers of participants and different participants' geographical distributions. Simulation results show that our resource allocation mechanism can compose the media handling services and allocate the required resources in an optimal manner while honoring the QoS in terms of end-to-end delay.


title: The Semantic Web Rule Language Expressiveness Extensions-A Survey

abstract: The Semantic Web Rule Language (SWRL) is a direct extension of OWL 2 DL with a subset of RuleML, and it is designed to be the rule language of the Semantic Web. This paper explores the state-of-the-art of SWRL's expressiveness extensions proposed over time. As a motivation, the effectiveness of the SWRL/OWL combination in modeling domain facts is discussed while some of the common expressive limitations of the combination are also highlighted. The paper then classifies and presents the relevant language extensions of the SWRL and their added expressive powers to the original SWRL definition. Furthermore, it provides a comparative analysis of the syntax and semantics of the proposed extensions. In conclusion, the decidability requirement and usability of each expressiveness extension are evaluated towards an efficient inclusion into the OWL ontologies.


title: Design, Fabrication and Characterization of FeAl-based Metallic-Intermetallic Laminate (MIL) Composites

abstract: FeAl-based MIL composites of various iron alloys were fabricated with an innovative 'multiple-thin-foil' configuration and 'two-stage reaction' strategy.


title: Skill Acquisition via Automated Multi-Coordinate Cost Balancing

abstract: We propose a learning framework, named Multi-Coordinate Cost Balancing (MCCB), to address the problem of acquiring point-to-point movement skills from demonstrations. MCCB encodes demonstrations simultaneously in multiple differential coordinates that specify local geometric properties. MCCB generates reproductions by solving a convex optimization problem with a multi-coordinate cost function and linear constraints on the reproductions, such as initial, target, and via points. Further, since the relative importance of each coordinate system in the cost function might be unknown for a given skill, MCCB learns optimal weighting factors that balance the cost function. We demonstrate the effectiveness of MCCB via detailed experiments conducted on one handwriting dataset and three complex skill datasets.


title: Radiological images and machine learning: trends, perspectives, and prospects

abstract: The application of machine learning to radiological images is an increasingly active research area that is expected to grow in the next five to ten years. Recent advances in machine learning have the potential to recognize and classify complex patterns from different radiological imaging modalities such as x-rays, computed tomography, magnetic resonance imaging and positron emission tomography imaging. In many applications, machine learning based systems have shown comparable performance to human decision-making. The applications of machine learning are the key ingredients of future clinical decision making and monitoring systems. This review covers the fundamental concepts behind various machine learning techniques and their applications in several radiological imaging areas, such as medical image segmentation, brain function studies and neurological disease diagnosis, as well as computer-aided systems, image registration, and content-based image retrieval systems. Synchronistically, we will briefly discuss current challenges and future directions regarding the application of machine learning in radiological imaging. By giving insight on how take advantage of machine learning powered applications, we expect that clinicians can prevent and diagnose diseases more accurately and efficiently.


title: Weyl points created by a three-dimensional flat band

abstract: Following the discovery of topological insulators (TIs), topological Dirac/Weyl semimetal has attracted much recent interest. A prevailing mechanism for the formation of Weyl points is by breaking time-reversal symmetry (TRS) or spatial inversion symmetry of a Dirac point. Here we demonstrate a generic formation mechanism for Weyl points by breaking TRS of a three-dimensional (3D) TI featured with highly degenerate 3D flat bands (FBs). It is in direct contrast to the conventional view that breaking TRS of a 2D/3D TI leads to a Chern insulator exhibiting quantum anomalous Hall effect. Based on a tight-binding model of pyrochlore lattice, we show that this unusual 3D-FB-enabled Weyl state may contain only a minimum of two Weyl points. Furthermore, using first-principles calculations, we identify this Weyl state in a real material Sn$_2$Nb$_2$O$_7$. The main features of the resulting Weyl points are analyzed with respect to symmetry, topological invariant and surface state. Our finding sheds new light on our fundamental understanding of topological physics and significantly extends the scope of Weyl semimetals to attract immediate experimental interest.


title: AutoSlim: Towards One-Shot Architecture Search for Channel Numbers

abstract: We study how to set channel numbers in a neural network to achieve better accuracy under constrained resources (e.g., FLOPs, latency, memory footprint or model size). A simple and one-shot solution, named AutoSlim, is presented. Instead of training many network samples and searching with reinforcement learning, we train a single slimmable network to approximate the network accuracy of different channel configurations. We then iteratively evaluate the trained slimmable model and greedily slim the layer with minimal accuracy drop. By this single pass, we can obtain the optimized channel configurations under different resource constraints. We present experiments with MobileNet v1, MobileNet v2, ResNet-50 and RL-searched MNasNet on ImageNet classification. We show significant improvements over their default channel configurations. We also achieve better accuracy than recent channel pruning methods and neural architecture search methods. Notably, by setting optimized channel numbers, our AutoSlim-MobileNet-v2 at 305M FLOPs achieves 74.2% top-1 accuracy, 2.4% better than default MobileNet-v2 (301M FLOPs), and even 0.2% better than RL-searched MNasNet (317M FLOPs). Our AutoSlim-ResNet-50 at 570M FLOPs, without depthwise convolutions, achieves 1.3% better accuracy than MobileNet-v1 (569M FLOPs). Code and models will be available at: https://github.com/JiahuiYu/slimmable_networks


title: Assessing Simulations of Imperial Dynamics and Conflict in the Ancient World

abstract: The development of models to capture large-scale dynamics in human history is one of the core contributions of cliodynamics. Most often, these models are assessed by their predictive capability on some macro-scale and aggregated measure and compared to manually curated historical data. In this report, we consider the model from Turchin et al. (2013), where the evaluation is done on the prediction of "imperial density": the relative frequency with which a geographical area belonged to large-scale polities over a certain time window. We implement the model and release both code and data for reproducibility. We then assess its behaviour against three historical data sets: the relative size of simulated polities vs historical ones; the spatial correlation of simulated imperial density with historical population density; the spatial correlation of simulated conflict vs historical conflict. At the global level, we show good agreement with population density ($R^2 < 0.75$), and some agreement with historical conflict in Europe ($R^2 < 0.42$). The model instead fails to reproduce the historical shape of individual polities. Finally, we tweak the model to behave greedily by having polities preferentially attacking weaker neighbours. Results significantly degrade, suggesting that random attacks are a key trait of the original model. We conclude by proposing a way forward by matching the probabilistic imperial strength from simulations to inferred networked communities from real settlement data.


title: Speed of sound of pure water to 700 MPa and an equation of state to 2300 MPa

abstract: Sound speeds of pure fluid water are reported between 0.1 and 700 MPa, from 353 K down to the melting curves of the ice phases. The 2 sigma precision and accuracy of the new sound speed measurements are close to 0.02%, with an estimated pressure accuracy of 0.04% up to 700 MPa. Using additional published measurements, a new equation of state is derived extending from 240 to 500 K and from 0.1 to 2300 MPa, covering much of the sub-critical domain of water up to the ice VI - ice VII transition. Analysis of measurements and construction of the equation of state are accomplished with a flexible computational thermodynamic framework based on local basis functions in the form of tensor b splines. Relative to IAPWS-95 (the most comprehensive representation available), improvements in the accuracies of density, sound speed and specific heat are observed above 100 MPa, particularly near the solid-fluid phase boundaries. Accurate thermodynamic properties for aqueous solutions under an increasing range of pressures, temperatures and compositions are needed to address a variety of technical and scientific challenges: the present study provides measurements essential for improving the high-pressure and low temperature representation of water.


title: Spectral Measures of Spiked Random Matrices

abstract: We study two spiked models of random matrices under general frameworks corresponding respectively to additive deformation of random symmetric matrices and multiplicative perturbation of random covariance matrices. In both cases, the limiting spectral measure in the direction of an eigenvector of the perturbation leads to old and new results on the coordinates of eigenvectors.


title: Quantum back-action of variable-strength measurement

abstract: Measuring a quantum system can randomly perturb its state. The strength and nature of this back-action depends on the quantity which is measured. In a partial measurement performed by an ideal apparatus, quantum physics predicts that the system remains in a pure state whose evolution can be tracked perfectly from the measurement record. We demonstrate this property using a superconducting qubit dispersively coupled to a cavity traversed by a microwave signal. The back-action on the qubit state of a single measurement of both signal quadratures is observed and shown to produce a stochastic operation whose action is determined by the measurement result. This accurate monitoring of a qubit state is an essential prerequisite for measurement-based feedback control of quantum systems.


title: Multifractality in random networks with power-law decaying bond strengths

abstract: In this paper we demonstrate numerically that random networks whose adjacency matrices ${\bf A}$ are represented by a diluted version of the Power--Law Banded Random Matrix (PBRM) model have multifractal eigenfunctions. The PBRM model describes one--dimensional samples with random long--range bonds. The bond strengths of the model, which decay as a power--law, are tuned by the parameter $\mu$ as $A_{mn}\propto |m-n|^{-\mu}$; while the sparsity is driven by the average network connectivity $\alpha$: for $\alpha=0$ the vertices in the network are isolated and for $\alpha=1$ the network is fully connected and the PBRM model is recovered. Though it is known that the PBRM model has multifractal eigenfunctions at the critical value $\mu=\mu_c=1$, we clearly show [from the scaling of the relative fluctuation of the participation number $I_2$ as well as the scaling of the probability distribution functions $P(\ln I_2)$] the existence of the critical value $\mu_c\equiv \mu_c(\alpha)$ for $\alpha<1$. Moreover, we characterise the multifractality of the eigenfunctions of our random network model by the use of the corresponding multifractal dimensions $D_q$, that we compute from the finite network-size scaling of the typical eigenfunction participation numbers $\exp\left\langle\ln I_q \right\rangle$.


title: A Posteriori Probabilistic Bounds of Convex Scenario Programs with Validation Tests

abstract: Scenario programs have established themselves as efficient tools towards decision-making under uncertainty. To assess the quality of scenario-based solutions a posteriori, validation tests based on Bernoulli trials have been widely adopted in practice. However, to reach a theoretically reliable judgement of risk, one typically needs to collect massive validation samples. In this work, we propose new a posteriori bounds for convex scenario programs with validation tests, which are dependent on both realizations of support constraints and performance on out-of-sample validation data. The proposed bounds enjoy wide generality in that many existing theoretical results can be incorporated as particular cases. To facilitate practical use, a systematic approach for parameterizing a posteriori probability bounds is also developed, which is shown to possess a variety of desirable properties allowing for easy implementations and clear interpretations. By synthesizing comprehensive information about support constraints and validation tests, improved risk evaluation can be achieved for randomized solutions in comparison with existing a posteriori bounds. Case studies on controller design of aircraft lateral motion are presented to validate the effectiveness of the proposed a posteriori bounds.


title: Scattering amplitudes from finite-volume spectral functions

abstract: A novel proposal is outlined to determine scattering amplitudes from finite-volume spectral functions. The method requires extracting smeared spectral functions from finite-volume Euclidean correlation functions, with a particular complex smearing kernel of width $\epsilon$ which implements the standard $i\epsilon$-prescription. In the $L \to \infty$ limit these smeared spectral functions are therefore equivalent to Minkowskian correlators with a specific time ordering to which a modified LSZ reduction formalism can be applied. The approach is presented for general $m \to n$ scattering amplitudes (above arbitrary inelastic thresholds) for a single-species real scalar field, although generalization to arbitrary spins and multiple coupled channels is likely straightforward. Processes mediated by the single insertion of an external current are also considered. Numerical determination of the finite-volume smeared spectral function is discussed briefly and the interplay between the finite volume, Euclidean signature, and time-ordered $i\epsilon$-prescription is illustrated perturbatively in a toy example.


title: Constraints on non-minimal coupling from quantum cosmology

abstract: Quantum cosmology is investigated in a de Sitter minisuperspace model with a quantized scalar field non-minimally coupled to curvature. Quantum states of the scalar field must satisfy the regularity condition, which requires that the probability of field fluctuations should not increase with their amplitude. We show that this condition imposes constraints on the allowed values of the curvature coupling parameter $\xi$. This is a surprising result, since the field dynamics depends only on the combination $m^2+\xi R$, where $m$ is the field mass and $R = \mathrm{const}$ is the curvature, and does not depend on $\xi$ separately.


title: Roberge-Weiss periodicity, canonical sector and modified Polyakov-loop

abstract: To obtain deeper understanding of QCD properties at finite temperature, we consider the Fourier decomposition of the grand-canonical partition function based on the canonical ensemble method via the imaginary chemical potential. Expectation values are, then, represented by summation over each canonical sector. We point out that the modified Polyakov-loop can play an important role in the canonical ensemble; for example, the Polyakov-loop paradox which is known in the canonical ensemble method can be evaded by considering the quantity. In addition, based on the periodicity issue of the modified Polyakov-loop at finite imaginary chemical potential, we can construct the systematic way to compute the dual quark condensate which has strong unclearness in its foundation in the presence of dynamical quarks so far.


title: Strong bound between trace distance and Hilbert-Schmidt distance for low-rank states

abstract: The trace distance between two quantum states, $\rho$ and $\sigma$, is an operationally meaningful quantity in quantum information theory. However, in general it is difficult to compute, involving the diagonalization of $\rho - \sigma$. In contrast, the Hilbert-Schmidt distance can be computed without diagonalization, although it is less operationally significant. Here, we relate the trace distance and the Hilbert-Schmidt distance with a bound that is particularly strong when either $\rho$ or $\sigma$ is low rank. Our bound is stronger than the bound one could obtain via the norm equivalence of the Frobenius and trace norms. We also consider bounds that are useful not only for low-rank states but also for low-entropy states. Our results have relevance to quantum information theory, quantum algorithms design, and quantum complexity theory.


title: Some results on the optimal matching problem for the Jacobi model

abstract: We establish some exact asymptotic results for a matching problem with respect to a family of beta distributions. Let $X_1, \ldots, X_n$ be independent random variables with common distribution the symmetric Jacobi measure $d\mu (x) = C_d (1-x^2)^{\frac d2 -1} dx$ with dimension $ d \geq 1$ on $[-1, 1]$, and let $\mu_n = \frac{1}{n} \sum_{i = 1}^{n} \delta_{X_i}$ be the associated empirical measure. We show that $\lim_{n \to \infty} n\E \left[ W_2^2( \mu^n, \mu ) \right] = \sum_{k = 1}^{\infty} \frac{1}{k(k+d-1)}$, where $W_2$ is the quadratic Kantorovich distance with respect to the intrinsic cost $\rho(x, y) = |\arccos(x) - \arccos (y)|$, $(x, y) \in [-1, 1]^2$, associated to the model. When $\mu$ is the product measure of two Jacobi measures with dimensions $d$ and $d'$ respectively, then $\E \left[ W_2^2( \mu^n, \mu ) \right] \approx \frac{\log n}{n}$. In the particular case $d = d' = 1$ (corresponding to the product of arcsine laws), $\lim_{n \to \infty} \frac{n}{\log n} \E \left[ W_2^2( \mu^n, \mu ) \right] = \frac{\pi}{4}$. Similar results do hold for non-symmetric Jacobi distributions. The proofs are based on the recent PDE and mass transportation approach developed by L.~Ambrosio, F.~Stra and D.~Trevisan.


title: Maxima and minima of homogeneous Gaussian random fields over continuous time and uniform grids

abstract: In this paper, for centered homogeneous Gaussian random fields the joint limiting distributions of normalized maxima and minima over continuous time and uniform grids are investigated. It is shown that maxima and minima are asymptotic dependent for strongly dependent homogeneous Gaussian random field with the choice of sparse grid, Pickands' grid or dense grid, while for the weakly dependent Gaussian random field maxima and minima are asymptotically independent.


title: InfoMask: Masked Variational Latent Representation to Localize Chest Disease

abstract: The scarcity of richly annotated medical images is limiting supervised deep learning based solutions to medical image analysis tasks, such as localizing discriminatory radiomic disease signatures. Therefore, it is desirable to leverage unsupervised and weakly supervised models. Most recent weakly supervised localization methods apply attention maps or region proposals in a multiple instance learning formulation. While attention maps can be noisy, leading to erroneously highlighted regions, it is not simple to decide on an optimal window/bag size for multiple instance learning approaches. In this paper, we propose a learned spatial masking mechanism to filter out irrelevant background signals from attention maps. The proposed method minimizes mutual information between a masked variational representation and the input while maximizing the information between the masked representation and class labels. This results in more accurate localization of discriminatory regions. We tested the proposed model on the ChestX-ray8 dataset to localize pneumonia from chest X-ray images without using any pixel-level or bounding-box annotations.


title: Global existence of solutions of initial-boundary value problem for nonlocal parabolic equation with nonlocal boundary condition

abstract: We prove global existence and blow-up of solutions of initial-boundary value problem for nonlinear nonlocal parabolic equation with nonlinear nonlocal boundary condition. Obtained results depend on the behavior of variable coefficients for large values of time.


title: A hint of matter underdensity at low $z$?

abstract: The $\Lambda$CDM cosmological model provides to first approximation a good description of the universe, but various tensions with data, most notably Hubble tension, persist. In this work we confront $\Lambda$CDM with the Pantheon Type Ia supernovae dataset and perform a two-parameter fit of the distance modulus for a running cut-off $z_{\textrm{max}}$. We observe that in a window between $z_{\textrm{max}} \approx 0.1$ and $z_{\textrm{max}} \approx 0.16$ there is a 1 - 2 $\sigma$ discrepancy with the Planck value $\omega_m = 0.315 \pm 0.007$, which points to a potential matter underdensity. For high-energy theorists, the analysis appears to support the de Sitter Swampland conjecture.


title: Tomography of the Cosmic Dawn and Reionization Eras with Multiple Tracers

abstract: The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.


title: Approximate spectral gaps for Markov chains mixing times in high dimensions

abstract: This paper introduces a concept of approximate spectral gap to analyze the mixing time of Markov Chain Monte Carlo (MCMC) algorithms for which the usual spectral gap is degenerate or almost degenerate. We use the idea to analyze a class of MCMC algorithms to sample from mixtures of densities. As an application we study the mixing time of a Gibbs sampler for variable selection in linear regression models. Under some regularity conditions on the signal and the design matrix of the regression problem, we show that for well-chosen initial distributions the mixing time of the Gibbs sampler is polynomial in the dimension of the space.


title: A generalization of convergence actions

abstract: Let a group $G$ act properly discontinuously and cocompactly on a locally compact space $X$. A Hausdorff compact space $Z$ that contains $X$ as an open subspace has the perspectivity property if the action $G\curvearrowright X$ extends to an action $G\curvearrowright Z$, by homeomorphisms, such that for every compact $K\subseteq X$ and every element $u$ of the unique uniform structure compatible with the topology of $Z$, the set $\{gK: g \in G\}$ has finitely many non $u$-small sets. We describe a correspondence between the compact spaces with the perspectivity property with respect to $X$ (and the fixed action of $G$ on it) and the compact spaces with the perspectivity property with respect to $G$ (and the left multiplication on itself). This generalizes a similar result for convergence group actions.


title: Entanglement in a 20-Qubit Superconducting Quantum Computer

abstract: Towards realising larger scale quantum algorithms, the ability to prepare sizeable multi-qubit entangled states with full qubit control is used as a benchmark for quantum technologies. We investigate the extent to which entanglement is found within a prepared graph state on the 20-qubit superconducting quantum computer, IBM Q Poughkeepsie. We prepared a graph state along a path consisting of all twenty qubits within Poughkeepsie and performed full quantum state tomography on all groups of four connected qubits along this path. We determined that each pair of connected qubits was inseparable and hence the prepared state was entangled. Additionally, a genuine multipartite entanglement witness was measured on all qubit subpaths of the graph state and we found genuine multipartite entanglement on chains of up to three qubits.


title: Medical Time Series Classification with Hierarchical Attention-based Temporal Convolutional Networks: A Case Study of Myotonic Dystrophy Diagnosis

abstract: Myotonia, which refers to delayed muscle relaxation after contraction, is the main symptom of myotonic dystrophy patients. We propose a hierarchical attention-based temporal convolutional network (HA-TCN) for myotonic dystrohpy diagnosis from handgrip time series data, and introduce mechanisms that enable model explainability. We compare the performance of the HA-TCN model against that of benchmark TCN models, LSTM models with and without attention mechanisms, and SVM approaches with handcrafted features. In terms of classification accuracy and F1 score, we found all deep learning models have similar levels of performance, and they all outperform SVM. Further, the HA-TCN model outperforms its TCN counterpart with regards to computational efficiency regardless of network depth, and in terms of performance particularly when the number of hidden layers is small. Lastly, HA-TCN models can consistently identify relevant time series segments in the relaxation phase of the handgrip time series, and exhibit increased robustness to noise when compared to attention-based LSTM models.


title: Distributed Algorithms for Fully Personalized PageRank on Large Graphs

abstract: Personalized PageRank (PPR) has enormous applications, such as link prediction and recommendation systems for social networks, which often require the fully PPR to be known. Besides, most of real-life graphs are edge-weighted, e.g., the interaction between users on the Facebook network. However, it is computationally difficult to compute the fully PPR, especially on large graphs, not to mention that most existing approaches do not consider the weights of edges. In particular, the existing approach cannot handle graphs with billion edges on a moderate-size cluster. To address this problem, this paper presents a novel study on the computation of fully edge-weighted PPR on large graphs using the distributed computing framework. Specifically, we employ the Monte Carlo approximation that performs a large number of random walks from each node of the graph, and exploits the parallel pipeline framework to reduce the overall running time of the fully PPR. Based on that, we develop several optimization techniques which (i) alleviate the issue of large nodes that could explode the memory space, (ii) pre-compute short walks for small nodes that largely speedup the computation of random walks, and (iii) optimize the amount of random walks to compute in each pipeline that significantly reduces the overhead. With extensive experiments on a variety of real-life graph datasets, we demonstrate that our solution is several orders of magnitude faster than the state-of-the-arts, and meanwhile, largely outperforms the baseline algorithms in terms of accuracy.


title: Navigation in the Presence of Obstacles for an Agile Autonomous Underwater Vehicle

abstract: Navigation underwater traditionally is done by keeping a safe distance from obstacles, resulting in "fly-overs" of the area of interest. Movement of an autonomous underwater vehicle (AUV) through a cluttered space, such as a shipwreck or a decorated cave, is an extremely challenging problem that has not been addressed in the past. This paper proposes a novel navigation framework utilizing an enhanced version of Trajopt for fast 3D path-optimization planning for AUVs. A sampling-based correction procedure ensures that the planning is not constrained by local minima, enabling navigation through narrow spaces. Two different modalities are proposed: planning with a known map results in efficient trajectories through cluttered spaces; operating in an unknown environment utilizes the point cloud from the visual features detected to navigate efficiently while avoiding the detected obstacles. The proposed approach is rigorously tested, both on simulation and in-pool experiments, proven to be fast enough to enable safe real-time 3D autonomous navigation for an AUV.


title: Regularized Stochastic Block Model for robust community detection in complex networks

abstract: The stochastic block model is able to generate different network partitions, ranging from traditional assortative communities to disassortative structures. Since the degree-corrected stochastic block model does not specify which mixing pattern is desired, the inference algorithms, which discover the most likely partition of the networks nodes, are likely to get trapped in the local optima of the log-likelihood. Here we introduce a new model constraining nodes' internal degrees ratios in the objective function to stabilize the inference of block models from the observed network data. Given the regularized model, the inference algorithms, such as Markov chain Monte Carlo, reliably finds assortative or disassortive structure as directed by the value of a single parameter. We show experimentally that the inference of our proposed model quickly converges to the desired assortative or disassortative partition while the inference of degree-corrected stochastic block model gets often trapped at the inferior local optimal partitions when the traditional assortative community structure is not strong in the observed networks.


title: ThunderNet: Towards Real-time Generic Object Detection

abstract: Real-time generic object detection on mobile platforms is a crucial but challenging computer vision task. However, previous CNN-based detectors suffer from enormous computational cost, which hinders them from real-time inference in computation-constrained scenarios. In this paper, we investigate the effectiveness of two-stage detectors in real-time generic detection and propose a lightweight two-stage detector named ThunderNet. In the backbone part, we analyze the drawbacks in previous lightweight backbones and present a lightweight backbone designed for object detection. In the detection part, we exploit an extremely efficient RPN and detection head design. To generate more discriminative feature representation, we design two efficient architecture blocks, Context Enhancement Module and Spatial Attention Module. At last, we investigate the balance between the input resolution, the backbone, and the detection head. Compared with lightweight one-stage detectors, ThunderNet achieves superior performance with only 40% of the computational cost on PASCAL VOC and COCO benchmarks. Without bells and whistles, our model runs at 24.1 fps on an ARM-based device. To the best of our knowledge, this is the first real-time detector reported on ARM platforms. Our code and models are available at \url{https://github.com/qinzheng93/ThunderNet}.


title: Note on a solution to domain wall problem with the Lazarides-Shafi mechanism in axion dark matter models

abstract: Axion is a promising candidate of dark matter. After the Peccei-Quinn symmetry breaking, axion strings are formed and attached by domain walls when the temperature of the universe becomes comparable to the QCD scale. Such objects can cause cosmological disasters if they are long-lived. As a solution for it, the Lazarides-Shafi mechanism is often discussed through introduction of a new non-Abelian (gauge) symmetry. We study this mechanism in detail and show configuration of strings and walls. Even if Abelian axion strings with a domain wall number greater than one are formed in the early universe, each of them is split into multiple Alice axion strings due to a repulsive force between the Alice strings even without domain wall. When domain walls are formed as the universe cools down, a single Alice string can be attached by a single wall because a vacuum is connected by a non-Abelian rotation without changing energy. Even if an Abelian axion string attached by domain walls are created due to the Kibble Zurek mechanism at the chiral phase transition, such strings are also similarly split into multiple Alice strings attached by walls in the presence of the domain wall tension. Such walls do not form stable networks since they collapse by the tension of the walls, emitting axions.


title: Euler-Maruyama Approximations for Stochastic McKean-Vlasov Equations with Non-Lipschitz Coefficients

abstract: In this paper we study a type of stochastic McKean-Vlasov equations with non-Lipschitz coefficients. Firstly, by an Euler-Maruyama approximation existence of its weak solutions is proved. And then we observe pathwise uniqueness of its weak solutions. Finally, it is shown that the Euler-Maruyama approximation has an optimal strong convergence rate.


title: Extremal isosystolic metrics with multiple bands of crossing geodesics

abstract: We apply recently developed convex programs to find the minimal-area Riemannian metric on $2n$-sided polygons ($n\geq 3$) with length conditions on curves joining opposite sides. We argue that the Riemannian extremal metric coincides with the conformal extremal metric on the regular $2n$-gon. The hexagon was considered by Calabi. The region covered by the maximal number $n$ of geodesics bands extends over most of the surface and exhibits positive curvature. As $n\to \infty$ the metric, away from the boundary, approaches the well-known round extremal metric on $\mathbb{RP}_2$. We extend Calabi's isosystolic variational principle to the case of regions with more than three bands of systolic geodesics. The extremal metric on $\mathbb{RP}_2$ is a stationary point of this functional applied to a surface with infinite number of systolic bands.


title: Models and Simulations for the Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC)

abstract: We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.


title: Monotonicity of the principal eigenvalue for a linear time-periodic parabolic operator

abstract: We investigate the effect of frequency on the principal eigenvalue of a time-periodic parabolic operator with Dirichlet, Robin or Neumann boundary conditions. The monotonicity and asymptotic behaviors of the principal eigenvalue with respect to the frequency parameter are established. Our results prove a conjecture raised by Hutson, Michaikow and Pol\'{a}\v{c}ik [2001 J. Math. Biol.].


title: Quasinormal Modes of Massive Scalar Field with Nonminimal Coupling in Higher-Dimensional de Sitter Black Hole with Single Rotation

abstract: We analytically investigate the quasinormal modes of the massive scalar field with a nonminimal coupling in the higher-dimensional de Sitter black hole with a single rotation. According to the separated scalar field equation, the boundary conditions of quasinormal modes are well constructed at the outer and cosmological horizons. Then, under near-extremal conditions, where the outer horizon closes to the cosmological horizon, the quasinormal frequencies are obtained and generalized to universal form in the higher-dimensional spacetime. Here, the real part of the frequency includes the scalar field contents, and its imaginary part only depends on the surface gravity at the outer horizon of the black hole.


title: Identities for Bernoulli polynomials related to multiple Tornheim zeta functions

abstract: We show that each member of a doubly infinite sequence of highly nonlinear expressions of Bernoulli polynomials, which can be seen as linear combinations of certain higher-order convolutions, is a multiple of a specific product of linear factors. The special case of Bernoulli numbers has important applications in the study of multiple Tornheim zeta functions. The proof of the main result relies on properties of Eulerian polynomials and higher-order Bernoulli polynomials.


title: Testing the Gravitational Weak Equivalence Principle in the Standard-Model Extension with Binary Pulsars

abstract: The Standard-Model Extension provides a framework to systematically investigate possible violation of the Lorentz symmetry. Concerning gravity, the linearized version was extensively examined. We here cast the first set of experimental bounds on the nonlinear terms in the field equation from the anisotropic cubic curvature couplings. These terms introduce body-dependent accelerations for self-gravitating objects, thus violating the gravitational weak equivalence principle (GWEP). Novel phenomena, that are absent in the linearized gravity, remain experimentally unexplored. We constrain them with precise binary-orbit measurements from pulsar timing, wherein the high density and large compactness of neutron stars are crucial for the test. It is the first study that seeks GWEP-violating signals in a fully anisotropic framework with Lorentz violation.


title: Phase space classification of an Ising Cellular Automaton: the Q2R model

abstract: An exact characterization of the different dynamical behavior that exhibit the space phase of a reversible and conservative cellular automaton, the so called Q2R model, is shown in this paper. Q2R is a cellular automaton which is a dynamical variation of the Ising model in statistical physics and whose space of configurations grows exponentially with the system size. As a consequence of the intrinsic reversibility of the model, the phase space is composed only by configurations that belong to a fixed point or a limit cycle. In this work we classify them in four types accordingly to well differentiated topological characteristics. Three of them, which we call of type S-I, S-II and S-III, share a symmetry property, while the fourth, which we call of type AS, does not. Specifically, we prove that any configuration of Q2R belongs to one of the four previous limit cycles. Moreover, at a combinatorial level, we are able to determine the number of limit cycles for some small periods which are almost always present in the Q2R. Finally, we provide a general overview of the resulting decomposition of the arbitrary size Q2R phase space, in addition, we realize an exhaustive study of a small Ising system 4x4 which is fully analyzed under this new framework.


title: An introduction to a new space debris evolution model -- SOLEM

abstract: SOLEM is the first space debris long-term evolution model of China. This paper describes the principles, components and workflow of the SOLEM. The effects of different mitigation measures based on SOLEM model are analyzed and presented. The limitation of the model is pointed out and its future improvement work-plan is prospected.


title: Finite Time Encryption Schedule in the Presence of an Eavesdropper with Operation Cost

abstract: In this paper, we consider a remote state estimation problem in the presence of an eavesdropper. A smart sensor takes measurement of a discrete linear time-invariant (LTI) process and sends its local state estimate through a wireless network to a remote estimator. An eavesdropper can overhear the sensor transmissions with a certain probability. To enhance the system privacy level, we propose a novel encryption strategy to minimize a linear combination of the expected error covariance at the remote estimator and the negative of the expected error covariance at the eavesdropper, taking into account the cost of the encryption process. We prove the existence of an optimal deterministic and Markovian policy for such an encryption strategy over a finite time horizon. Two situations, namely, with or without knowledge of the eavesdropper estimation error covariance are studied and the optimal schedule is shown to satisfy the threshold-like structure in both cases. Numerical examples are given to illustrate the results.


title: When an attacker meets a cipher-image in 2018: A Year in Review

abstract: This paper aims to review the encountered technical contradictions when an attacker meets the cipher-images encrypted by the image encryption schemes (algorithms) proposed in 2018 from the viewpoint of an image cryptanalyst. The most representative works among them are selected and classified according to their essential structures. Almost all image cryptanalysis works published in 2018 are surveyed due to their small number. The challenging problems on design and analysis of image encryption schemes are summarized to receive the attentions of both designers and attackers (cryptanalysts) of image encryption schemes, which may promote solving scenario-oriented image security problems with new technologies.


title: Connecting Program Synthesis and Reachability: Automatic Program Repair using Test-Input Generation

abstract: We prove that certain formulations of program synthesis and reachability are equivalent. Specifically, our constructive proof shows the reductions between the template-based synthesis problem, which generates a program in a pre-specified form, and the reachability problem, which decides the reachability of a program location. This establishes a link between the two research fields and allows for the transfer of techniques and results between them. To demonstrate the equivalence, we develop a program repair prototype using reachability tools. We transform a buggy program and its required specification into a specific program containing a location reachable only when the original program can be repaired, and then apply an off-the-shelf test-input generation tool on the transformed program to find test values to reach the desired location. Those test values correspond to repairs for the original programm. Preliminary results suggest that our approach compares favorably to other repair methods.


title: New analytical tools for HDG in elasticity, with applications to elastodynamics

abstract: We present some new analytical tools for the error analysis of hybridizable discontinuous Galerkin (HDG) method for linear elasticity. These tools allow us to analyze more variants of HDG method using the projection-based approach, which renders the error analysis simple and concise. The key result is a tailored projection for the Lehrenfeld-Sch\"{o}berl type HDG (HDG+ for simplicity) methods. By using the projection we recover the error estimates of HDG+ for steady-state and time-harmonic elasticity in a simpler analysis. We also present a semi-discrete (in space) HDG+ method for transient elastic waves and prove it is uniformly-in-time optimal convergent by using the projection-based error analysis. Numerical experiments supporting our analysis are presented at the end.


title: Comment on "An X-ray chimney extending hundreds of parsecs above and below the Galactic Centre" (2019, Nature, 567, 34)

abstract: A recent article "An X-ray chimney extending hundreds of parsecs above and below the Galactic Centre" (2019, Nature, 567, 34) reported the detection of chimney-like X-ray-emitting features above and below the Galactic Center from XMM-Newton observations. We note here that these features were already reported by our Suzaku papers: Nakashima et al. (2013, ApJ, 773, 20, arXiv:1310.4236) for the southern feature and Nakashima et al. (2019, ApJ, in press, arXiv:1903.02571) for the northern feature. In particular, Nakashima et al. (2013) show that the ionization state of the southern feature is not in collisional ionization equilibrium and is most likely in a recombining or over-ionized state, which suggests its origin in the Galactic Center about 0.1 Myr ago.


title: SymInfer: Inferring Program Invariants using Symbolic States

abstract: We introduce a new technique for inferring program invariants that uses symbolic states generated by symbolic execution. Symbolic states, which consist of path conditions and constraints on local variables, are a compact description of sets of concrete program states and they can be used for both invariant inference and invariant verification. Our technique uses a counterexample-based algorithm that creates concrete states from symbolic states, infers candidate invariants from concrete states, and then verifies or refutes candidate invariants using symbolic states. The refutation case produces concrete counterexamples that prevent spurious results and allow the technique to obtain more precise invariants. This process stops when the algorithm reaches a stable set of invariants. We present SymInfer, a tool that implements these ideas to automatically generate invariants at arbitrary locations in a Java program. The tool obtains symbolic states from Symbolic PathFinder and uses existing algorithms to infer complex (potentially nonlinear) numerical invariants. Our preliminary results show that SymInfer is effective in using symbolic states to generate precise and useful invariants for proving program safety and analyzing program runtime complexity. We also show that SymInfer outperforms existing invariant generation systems.


title: Development of a coronal mass ejection arrival time forecasting system using interplanetary scintillation observations

abstract: Coronal mass ejections (CMEs) cause disturbances in the environment of the Earth when they arrive at the Earth. However, the prediction of the arrival of CMEs still remains a challenge. We have developed an interplanetary scintillation (IPS) estimation system based on a global magnetohydrodynamic (MHD) simulation of the inner heliosphere to predict the arrival time of CMEs. In this system, the initial speed of a CME is roughly derived from white light coronagraph observations. Then, the propagation of the CME is calculated by a global MHD simulation. The IPS response is estimated by the three-dimensional density distribution of the inner heliosphere derived from the MHD simulation. The simulated IPS response is compared with the actual IPS observations made by the Institute for Space-Earth Environmental Research, Nagoya University, and shows good agreement with that observed. We demonstrated how the simulation system works using a halo CME event generated by a X9.3 flare observed on September 5, 2017. We find that the CME simulation that best estimates the IPS observation can more accurately predict the time of arrival of the CME at the Earth. These results suggest that the accuracy of the CME arrival time can be improved if our current MHD simulations include IPS data.


title: An Improved Approach for Semantic Graph Composition with CCG

abstract: This paper builds on previous work using Combinatory Categorial Grammar (CCG) to derive a transparent syntax-semantics interface for Abstract Meaning Representation (AMR) parsing. We define new semantics for the CCG combinators that is better suited to deriving AMR graphs. In particular, we define relation-wise alternatives for the application and composition combinators: these require that the two constituents being combined overlap in one AMR relation. We also provide a new semantics for type raising, which is necessary for certain constructions. Using these mechanisms, we suggest an analysis of eventive nouns, which present a challenge for deriving AMR graphs. Our theoretical analysis will facilitate future work on robust and transparent AMR parsing using CCG.


title: A Large-Scale Multi-Length Headline Corpus for Analyzing Length-Constrained Headline Generation Model Evaluation

abstract: Browsing news articles on multiple devices is now possible. The lengths of news article headlines have precise upper bounds, dictated by the size of the display of the relevant device or interface. Therefore, controlling the length of headlines is essential when applying the task of headline generation to news production. However, because there is no corpus of headlines of multiple lengths for a given article, previous research on controlling output length in headline generation has not discussed whether the system outputs could be adequately evaluated without multiple references of different lengths. In this paper, we introduce two corpora, which are Japanese News Corpus (JNC) and JApanese MUlti-Length Headline Corpus (JAMUL), to confirm the validity of previous evaluation settings. The JNC provides common supervision data for headline generation. The JAMUL is a large-scale evaluation dataset for headlines of three different lengths composed by professional editors. We report new findings on these corpora; for example, although the longest length reference summary can appropriately evaluate the existing methods controlling output length, this evaluation setting has several problems.


title: Constrained deformations of positive scalar curvature metrics

abstract: We present a series of results concerning the interplay between the scalar curvature of a manifold and the mean curvature of its boundary. In particular, we give a complete topological characterization of those compact 3-manifolds that support Riemannian metrics of positive scalar curvature and mean-convex boundary and, in any such case, we prove that the associated moduli space of metrics is path-connected. The methods we employ are flexible enough to allow the construction of continuous paths of positive scalar curvature metrics with minimal boundary, and to derive similar conclusions in that context as well. Our work relies on a combination of earlier fundamental contributions by Gromov-Lawson and Schoen-Yau, on the smoothing procedure designed by Miao, and on the interplay of Perelman's Ricci flow with surgery and conformal deformation techniques introduced by Cod\'a Marques in dealing with the closed case.


title: Dimensionality and Anisotropicity Dependence of Radiative Recombination in Nanostructured Phosphorene

abstract: The interplay between dimensionality and anisotropicity leads to intriguing optoelectronic properties and exciton dynamics in low dimensional semiconductors. In this study we use nanostructured phosphorene as a prototypical example to unfold such complex physics and develop a general first-principles framework to study exciton dynamics in low dimensional systems. Specifically we derived the radiative lifetime and light emission intensity from 2D to 0D systems based on many-body perturbation theory, and investigated the dimensionality and anisotropicity effects on radiative recombination lifetime both at 0 K and finite temperature, as well as polarization and angle dependence of emitted light. We show that the radiative lifetime at 0 K increases by an order of $10^3$ with the lowering of one dimension (i.e. from 2D to 1D nanoribbons or from 1D to 0D quantum dots). We also show that obtaining the radiative lifetime at finite temperature requires accurate exciton dispersion beyond the effective mass approximation. Finally, we demonstrate that monolayer phosphorene and its nanostructures always emit linearly polarized light consistent with experimental observations, different from in-plane isotropic 2D materials like MoS2 and h-BN that can emit light with arbitrary polarization, which may have important implications for quantum information applications.


title: How to pick the domain randomization parameters for sim-to-real transfer of reinforcement learning policies?

abstract: Recently, reinforcement learning (RL) algorithms have demonstrated remarkable success in learning complicated behaviors from minimally processed input. However, most of this success is limited to simulation. While there are promising successes in applying RL algorithms directly on real systems, their performance on more complex systems remains bottle-necked by the relative data inefficiency of RL algorithms. Domain randomization is a promising direction of research that has demonstrated impressive results using RL algorithms to control real robots. At a high level, domain randomization works by training a policy on a distribution of environmental conditions in simulation. If the environments are diverse enough, then the policy trained on this distribution will plausibly generalize to the real world. A human-specified design choice in domain randomization is the form and parameters of the distribution of simulated environments. It is unclear how to the best pick the form and parameters of this distribution and prior work uses hand-tuned distributions. This extended abstract demonstrates that the choice of the distribution plays a major role in the performance of the trained policies in the real world and that the parameter of this distribution can be optimized to maximize the performance of the trained policies in the real world


title: Atrial Fibrillation Detection Using Deep Features and Convolutional Networks

abstract: Atrial fibrillation is a cardiac arrhythmia that affects an estimated 33.5 million people globally and is the potential cause of 1 in 3 strokes in people over the age of 60. Detection and diagnosis of atrial fibrillation (AFIB) is done noninvasively in the clinical environment through the evaluation of electrocardiograms (ECGs). Early research into automated methods for the detection of AFIB in ECG signals focused on traditional bio-medical signal analysis to extract important features for use in statistical classification models. Artificial intelligence models have more recently been used that employ convolutional and/or recurrent network architectures. In this work, significant time and frequency domain characteristics of the ECG signal are extracted by applying the short-time Fourier trans-form and then visually representing the information in a spectrogram. Two different classification approaches were investigated that utilized deep features in the spectrograms construct-ed from ECG segments. The first approach used a pretrained DenseNet model to extract features that were then classified using Support Vector Machines, and the second approach used the spectrograms as direct input into a convolutional network. Both approaches were evaluated against the MIT-BIH AFIB dataset, where the convolutional network approach achieved a classification accuracy of 93.16%. While these results do not surpass established automated atrial fibrillation detection methods, they are promising and warrant further investigation given they did not require any noise prefiltering, hand-crafted features, nor a reliance on beat detection.


title: Replica symmetry breaking in random matrix model: a toy model of wormhole networks

abstract: We study the replica symmetry breaking (RSB) in the Gaussian Unitary Ensemble (GUE) random matrix model. We find that the RSB occurs at the transition temperature $T_{\text{RSB}}\sim N^{-2/3}$ in the large $N$ limit. We argue that this transition originates from a landscape of almost degenerate local minima of free energy labeled by Young diagrams. We also discuss a possible implication of our findings for the multiply-coupled Sachdev-Ye-Kitaev models and their holographic dual Jackiw-Teitelboim gravity from the viewpoint of ER=EPR conjecture.


title: What you get is what you see: Decomposing Epistemic Planning using Functional STRIPS

abstract: Epistemic planning --- planning with knowledge and belief --- is essential in many multi-agent and human-agent interaction domains. Most state-of-the-art epistemic planners solve this problem by compiling to propositional classical planning, for example, generating all possible knowledge atoms, or compiling epistemic formula to normal forms. However, these methods become computationally infeasible as problems grow. In this paper, we decompose epistemic planning by delegating reasoning about epistemic formula to an external solver. We do this by modelling the problem using \emph{functional STRIPS}, which is more expressive than standard STRIPS and supports the use of external, black-box functions within action models. Exploiting recent work that demonstrates the relationship between what an agent `sees' and what it knows, we allow modellers to provide new implementations of externals functions. These define what agents see in their environment, allowing new epistemic logics to be defined without changing the planner. As a result, it increases the capability and flexibility of the epistemic model itself, and avoids the exponential pre-compilation step. We ran evaluations on well-known epistemic planning benchmarks to compare with an existing state-of-the-art planner, and on new scenarios based on different external functions. The results show that our planner scales significantly better than the state-of-the-art planner against which we compared, and can express problems more succinctly.


title: Beam energy dependence of (anti-)deuteron production in Au+Au collisions at RHIC

abstract: We report the energy dependence of mid-rapidity (anti-)deuteron production in Au+Au collisions at $\sqrt{s_\text{NN}} =\ $7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV, measured by the STAR experiment at RHIC. The yield of deuterons is found to be well described by the thermal model. The collision energy, centrality, and transverse momentum dependence of the coalescence parameter $B_2$ are discussed. We find that the values of $B_2$ for anti-deuterons are systematically lower than those for deuterons, indicating that the correlation volume of anti-baryons is larger than that of baryons at $\sqrt{s_\text{NN}}$ from 19.6 to 39 GeV. In addition, values of $B_2$ are found to vary with collision energy and show a broad minimum around $\sqrt{s_\text{NN}}=\ $20 to 40 GeV, which might imply a change of the equation of state of the medium in these collisions.


title: BubbleNets: Learning to Select the Guidance Frame in Video Object Segmentation by Deep Sorting Frames

abstract: Semi-supervised video object segmentation has made significant progress on real and challenging videos in recent years. The current paradigm for segmentation methods and benchmark datasets is to segment objects in video provided a single annotation in the first frame. However, we find that segmentation performance across the entire video varies dramatically when selecting an alternative frame for annotation. This paper address the problem of learning to suggest the single best frame across the video for user annotation-this is, in fact, never the first frame of video. We achieve this by introducing BubbleNets, a novel deep sorting network that learns to select frames using a performance-based loss function that enables the conversion of expansive amounts of training examples from already existing datasets. Using BubbleNets, we are able to achieve an 11% relative improvement in segmentation performance on the DAVIS benchmark without any changes to the underlying method of segmentation.


title: Wasserstein Dependency Measure for Representation Learning

abstract: Mutual information maximization has emerged as a powerful learning objective for unsupervised representation learning obtaining state-of-the-art performance in applications such as object recognition, speech recognition, and reinforcement learning. However, such approaches are fundamentally limited since a tight lower bound of mutual information requires sample size exponential in the mutual information. This limits the applicability of these approaches for prediction tasks with high mutual information, such as in video understanding or reinforcement learning. In these settings, such techniques are prone to overfit, both in theory and in practice, and capture only a few of the relevant factors of variation. This leads to incomplete representations that are not optimal for downstream tasks. In this work, we empirically demonstrate that mutual information-based representation learning approaches do fail to learn complete representations on a number of designed and real-world tasks. To mitigate these problems we introduce the Wasserstein dependency measure, which learns more complete representations by using the Wasserstein distance instead of the KL divergence in the mutual information estimator. We show that a practical approximation to this theoretically motivated solution, constructed using Lipschitz constraint techniques from the GAN literature, achieves substantially improved results on tasks where incomplete representations are a major challenge.


title: Technical Report: Optimal Control of Piecwise-smooth Control Systems via Singular Perturbations

abstract: This paper investigates optimal control problems formulated over a class of piecewise-smooth vector fields. Instead of optimizing over the discontinuous system directly, we instead formulate optimal control problems over a family of regularizations which are obtained by "smoothing out" the discontinuity in the original system. It is shown that the smooth problems can be used to obtain accurate derivative information about the non-smooth problem, under standard regularity conditions. We then indicate how the regularizations can be used to consistently approximate the non-smooth optimal control problem in the sense of Polak. The utility of these smoothing techniques is demonstrated in an in-depth example, where we utilize recently developed reduced-order modeling techniques from the dynamic walking community to generate motion plans across contact sequences for a 18-DOF model of a lower-body exoskeleton.


title: Anywhere Decoding: Low-Overhead Uplink Interference Management for Wireless Networks

abstract: Inter-cell interference (ICI) is one of the major performance-limiting factors in the context of modern cellular systems. To tackle ICI, coordinated multi-point (CoMP) schemes have been proposed as a key technology for next-generation mobile communication systems. Although CoMP schemes offer promising theoretical gains, their performance could degrade significantly because of practical issues such as limited backhaul. To address this issue, we explore a novel uplink interference management scheme called anywhere decoding, which requires exchanging just a few bits of information per coding interval among the base stations (BSs). In spite of the low overhead of anywhere decoding, we observe considerable gains in the outage probability performance of cell-edge users, compared to no cooperation between BSs. Additionally, asymptotic results of the outage probability for high-SNR regimes demonstrate that anywhere decoding schemes achieve full spatial diversity through multiple decoding opportunities, and they are within 1.5 dB of full cooperation.


title: A dataset for resolving referring expressions in spoken dialogue via contextual query rewrites (CQR)

abstract: We present Contextual Query Rewrite (CQR) a dataset for multi-domain task-oriented spoken dialogue systems that is an extension of the Stanford dialog corpus (Eric et al., 2017a). While previous approaches have addressed the issue of diverse schemas by learning candidate transformations (Naik et al., 2018), we instead model the reference resolution task as a user query reformulation task, where the dialog state is serialized into a natural language query that can be executed by the downstream spoken language understanding system. In this paper, we describe our methodology for creating the query reformulation extension to the dialog corpus, and present an initial set of experiments to establish a baseline for the CQR task. We have released the corpus to the public [1] to support further research in this area.


title: The initiation temperatures in nanothermite reactions

abstract: An assumption on general regularities and chemical mechanisms of solid-state reactions in nanofilms and nanothermite mixtures is made. It is demonstrated that the moving force of all Al-based nanothermite reactions is the synthesis of the Al2O3 phase with a high negative enthalpy of formation. According to this assumption, all Al-based nanothermite reactions should have the same initiation temperature Tin. Indeed, an analysis of literature data shows that the synthesis of the Al2O3 phase in Al/Fe2O3, Al/CuO, Al/Co3O4, Al/MoO3, Al/Bi2O3, and Al/NiO nanothermite mixtures starts at the same temperature Tin ~ 510{\deg}C. We also demonstrate the same initiation temperatures ~250{\deg}C, ~300{\deg}C, and ~180{\deg}C for Zr-, Mg- and In-based nanothermite reactions, respectively. It is predicted that nanothermite reactions based on other fuels have their own initiation temperatures.


title: A Fast Free-viewpoint Video Synthesis Algorithm for Sports Scenes

abstract: In this paper, we report on a parallel freeviewpoint video synthesis algorithm that can efficiently reconstruct a high-quality 3D scene representation of sports scenes. The proposed method focuses on a scene that is captured by multiple synchronized cameras featuring wide-baselines. The following strategies are introduced to accelerate the production of a free-viewpoint video taking the improvement of visual quality into account: (1) a sparse point cloud is reconstructed using a volumetric visual hull approach, and an exact 3D ROI is found for each object using an efficient connected components labeling algorithm. Next, the reconstruction of a dense point cloud is accelerated by implementing visual hull only in the ROIs; (2) an accurate polyhedral surface mesh is built by estimating the exact intersections between grid cells and the visual hull; (3) the appearance of the reconstructed presentation is reproduced in a view-dependent manner that respectively renders the non-occluded and occluded region with the nearest camera and its neighboring cameras. The production for volleyball and judo sequences demonstrates the effectiveness of our method in terms of both execution time and visual quality.


title: Breakdown of agreement between non-relativistic and relativistic quantum dynamical predictions for the periodically kicked rotor

abstract: As an extension of our recent work (Results in Physics 12, 147-152 (2019)), we show, both theoretically and numerically, that the breakdown of agreement between the non-relativistic and relativistic quantum dynamical predictions in the non-relativistic regime also occurs for an electron subjected to a time-dependent force.


title: Successive-Cancellation Decoding of Linear Source Code

abstract: This paper investigates the error probability of several decoding methods for a source code with decoder side information, where the decoding methods are: 1) symbol-wise maximum a posteriori decoding, 2) successive-cancellation decoding, and 3) stochastic successive-cancellation decoding. The proof of the effectiveness of a decoding method is reduced to that for an arbitrary decoding method, where `effective' means that the error probability goes to zero as $n$ goes to infinity. Furthermore, we revisit the polar source code showing that stochastic successive-cancellation decoding, as well as successive-cancellation decoding, is effective for this code.


title: Cherenkov Detectors Fast Simulation Using Neural Networks

abstract: We propose a way to simulate Cherenkov detector response using a generative adversarial neural network to bypass low-level details. This network is trained to reproduce high level features of the simulated detector events based on input observables of incident particles. This allows the dramatic increase of simulation speed. We demonstrate that this approach provides simulation precision which is consistent with the baseline and discuss possible implications of these results.


title: Step Change Improvement in ADMET Prediction with PotentialNet Deep Featurization

abstract: The Absorption, Distribution, Metabolism, Elimination, and Toxicity (ADMET) properties of drug candidates are estimated to account for up to 50% of all clinical trial failures. Predicting ADMET properties has therefore been of great interest to the cheminformatics and medicinal chemistry communities in recent decades. Traditional cheminformatics approaches, whether the learner is a random forest or a deep neural network, leverage fixed fingerprint feature representations of molecules. In contrast, in this paper, we learn the features most relevant to each chemical task at hand by representing each molecule explicitly as a graph, where each node is an atom and each edge is a bond. By applying graph convolutions to this explicit molecular representation, we achieve, to our knowledge, unprecedented accuracy in prediction of ADMET properties. By challenging our methodology with rigorous cross-validation procedures and prospective analyses, we show that deep featurization better enables molecular predictors to not only interpolate but also extrapolate to new regions of chemical space.


title: Universality class of the percolation in two-dimensional lattices with distortion

abstract: Mitra et al. [Phys. Rev. E 99 (2019) 012117] proposed a new percolation model that includes distortion in the square lattice and concluded that it may belong to the same universality class as the ordinary percolation. But the conclusion is questionable since their results of critical exponents are not consistent. In this paper, we reexamined the new model with high precision in the square, triangular, and honeycomb lattices by using the Newman-Ziff algorithm. Through the finite-size scaling, we obtained the percolation threshold of the infinite-size lattice and critical exponents ($\nu$ and $\beta$). Our results of the critical exponents are the same as those of the classical percolation within error bars, and the percolation in distorted lattices is confirmed to belong to the universality class of the classical percolation in two dimensions.


title: Hierarchical Pooling Structure for Weakly Labeled Sound Event Detection

abstract: Sound event detection with weakly labeled data is considered as a problem of multi-instance learning. And the choice of pooling function is the key to solving this problem. In this paper, we proposed a hierarchical pooling structure to improve the performance of weakly labeled sound event detection system. Proposed pooling structure has made remarkable improvements on three types of pooling function without adding any parameters. Moreover, our system has achieved competitive performance on Task 4 of Detection and Classification of Acoustic Scenes and Events (DCASE) 2017 Challenge using hierarchical pooling structure.


title: Coupling the Dirac and Einstein equations through geometry

abstract: We show that the exterior algebra bundle over a curved spacetime can be used as framework in which both the Dirac and the Einstein equations can be obtained. These equations and their coupling follow from the variational principle applied to a Lagrangian constructed from natural geometric invariants. We also briefly indicate how other forces can potentially be incorporated within this geometric framework.


title: The Prasad conjectures for U(2), SO(4) and Sp(4)

abstract: We will give a proof to the Prasad conjectures for $U_2$, $SO_4$ and $Sp_4$ over a quadratic field extension.


title: Magnitude homology of geodesic metric spaces with an upper curvature bound

abstract: In this article, we study the magnitude homology of geodesic metric spaces of curvature $\leq \kappa$, especially ${\rm CAT}(\kappa)$ spaces. We will show that the magnitude homology $MH^{l}_{n}(X)$ of such a meric space $X$ vanishes for small $l$ and all $n > 0$. Conseqently, we can compute a total $\mathbb{Z}$-degree magnitude homology for small $l$ for the shperes $\mathbb{S}^{n}$, the Euclid spaces $\mathbb{E}^{n}$, the hyperbolic spaces $\mathbb{H}^{n}$, and real projective spaces $\mathbb{RP}^{n}$ with the standard metric. We also show that an existence of closed geodesic in a metric space guarantees the non-triviality of magnitude homology.


title: Separation of time-scales for the seed bank diffusion and its jump-diffusion limit

abstract: We investigate the scaling limit of the seed bank diffusion when reproduction and migration (to and from the seed bank) happen on different time-scales. More precisely, we consider the case when migration is `slow' and reproduction is `standard' (in the original time-scale) and then switch to a new, accelerated time-scale, where migration is `standard' and reproduction is `fast'. This is motivated by models for bacterial dormancy, where periods of quiescence can be orders of magnitude larger than reproductive times, and where it is expected to find non-trivial degenerate genealogies on the evolutionary time-scale. However, the above scaling regime is not only interesting from a biological perspective, but also from a mathematical point of view, since it provides a prototypical example where the expected scaling limit of a continuous diffusion should (and will be) a jump-diffusion. For this situation, standard convergence results often seem to fail in multiple ways. For example, since the set of continuous paths from a closed subset of the c\`adl\`ag paths in each of the classical Skorohod topologies $J_1, J_2, M_1$ and $M_2$, none of them can be employed for tightness on path-space. Further, a na\"ive direct rescaling of the Markov generator corresponding to the continuous diffusion immediately leads to a blow-up of the diffusion coefficient. Still, one can identify a well-defined limit via duality in a surprisingly non-technical way. Indeed, we show that a certain duality relation is in some sense stable under passage to the limit and allows an identification of the limit, avoiding all technicalities related to the blow-up in the classical generator. The result then boils down to a convergence criterion for time-continuous Markov chains in a separation of time-scales regime, which is of independent interest.


title: Source-Frequency Phase-Referencing Observation of AGNs with KaVA Using Simultaneous Dual-Frequency Receiving

abstract: The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.


title: Estimating effective population size changes from preferentially sampled genetic sequences

abstract: Coalescent theory combined with statistical modeling allows us to estimate effective population size fluctuations from molecular sequences of individuals sampled from a population of interest. When sequences are sampled serially through time and the distribution of the sampling times depends on the effective population size, explicit statistical modeling of sampling times improves population size estimation. Previous work assumed that the genealogy relating sampled sequences is known and modeled sampling times as an inhomogeneous Poisson process with log-intensity equal to a linear function of the log-transformed effective population size. We improve this approach in two ways. First, we extend the method to allow for joint Bayesian estimation of the genealogy, effective population size trajectory, and other model parameters. Next, we improve the sampling time model by incorporating additional sources of information in the form of time-varying covariates. We validate our new modeling framework using a simulation study and apply our new methodology to analyses of population dynamics of seasonal influenza and to the recent Ebola virus outbreak in West Africa.


title: The characteristic variety for Feigin and Odesskii's elliptic algebras

abstract: This paper examines an algebraic variety that controls an important part of the structure and representation theory of the algebra $Q_{n,k}(E,\tau)$ introduced by Feigin and Odesskii. The $Q_{n,k}(E,\tau)$'s are a family of quadratic algebras depending on a pair of coprime integers $n>k\ge 1$, an elliptic curve $E$, and a point $\tau\in E$. It is already known that the structure and representation theory of $Q_{n,1}(E,\tau)$ is controlled by the geometry associated to $E$ embedded as a degree $n$ normal curve in the projective space $\mathbb P^{n-1}$, and by the way in which the translation automorphism $z\mapsto z+\tau$ interacts with that geometry. For $k\ge 2$ a similar phenomenon occurs: $(E,\tau)$ is replaced by $(X_{n/k},\sigma)$ where $X_{n/k}\subseteq\mathbb P^{n-1}$ is the characteristic variety of the title and $\sigma$ is an automorphism of it that is determined by the negative continued fraction for $\frac{n}{k}$. There is a surjective morphism $\Phi:E^g \to X_{n/k}$ where $g$ is the length of that continued fraction. The main result in this paper is that $X_{n/k}$ is a quotient of $E^g$ by the action of an explicit finite group. We also prove some assertions made by Feigin and Odesskii. The morphism $\Phi$ is the natural one associated to a particular invertible sheaf $\mathcal L_{n/k}$ on $E^g$. The generalized Fourier-Mukai transform associated to $\mathcal L_{n/k}$ sends the set of isomorphism classes of degree-zero invertible $\mathcal O_E$-modules to the set of isomorphism classes of indecomposable locally free $\mathcal O_E$-modules of rank $k$ and degree $n$. Thus $X_{n/k}$ has an importance independent of the role it plays in relation to $Q_{n,k}(E,\tau)$. The backward $\sigma$-orbit of each point on $X_{n/k}$ determines a point module for $Q_{n,k}(E,\tau)$.