title: Measuring the D/H Ratios of Exoplanets and Brown Dwarfs

abstract: The relative abundance of deuterium and hydrogen is a potent tracer of planet formation and evolution. Jupiter and Saturn have protosolar atmospheric D/H ratios, a relic of substantial gas accretion from the nebula, while the atmospheres of Neptune and Uranus are enhanced in D by accretion of ices into their envelopes. For terrestrial planets, D/H ratios are used to determine the mechanisms of volatile delivery and subsequent atmosphere loss over the lifetime of the planet. Planets and brown dwarfs more massive than 13 M_J quickly fuse their initial D reservoir. Here, we simulate spectra for giant exoplanets and brown dwarfs (2 M_Neptune to ~10 M_Jupiter) from Teff=200--1800 K including both CH3D and HDO to determine the observability of these dominant deuterium isotopologues in mid-infrared thermal emission spectra. Colder objects have stronger molecular features in their spectra, due to the temperature-dependence of molecular cross sections. CH3D is easier to observe than HDO at all temperatures considered, due to the strength of its absorption bands and locations of features at wavelengths with few other strong absorption features. We predict that for nearby cool brown dwarfs, the CH3D abundance can be measured at high signal-to-noise with the James Webb Space Telescope; for objects from 200--800 K closer than 10 pc, a protosolar D/H ratio would be readily observable in 2.5 hours. Moderately young Jupiter-mass planets (100--300 Myr) and young Neptunes (10 Myr) may be discovered with JWST and provide the best targets for detecting deuterium on an exoplanet in the coming decade. Future telescope designs should consider the importance of isotopes for understanding the formation and evolution of planetary atmospheres.


title: Iterative Nonlocal Residual Elasticity

abstract: Motivated by the existing complications of finding solutions of Eringen nonlocal model, an alternative model is developed here. The new formulation of the nonlocal elasticity is centered upon expressing the dynamic equilibrium requirements based on a nonlocal residual stress field. This new nonlocal elasticity is explained from the lattice mechanics and continuum mechanics points of view. Boundary value problems obtained based on the new nonlocal elasticity are solved following a proposed iterative procedure. This iterative procedure is centered upon correcting the solution of the classical field problem for the nonlocal residual field of the elastic domain. Convergence analyses are presented to show the convergence of the iterative procedure to the solution of the nonlocal field problem. The iterative procedure is an integrated part of the proposed nonlocal elasticity. Therefore, the newly developed nonlocal elasticity is give the name iterative nonlocal residual elasticity.


title: Autonomous Urban Localization and Navigation with Limited Information

abstract: Urban environments offer a challenging scenario for autonomous driving. Globally localizing information, such as a GPS signal, can be unreliable due to signal shadowing and multipath errors. Detailed a priori maps of the environment with sufficient information for autonomous navigation typically require driving the area multiple times to collect large amounts of data, substantial post-processing on that data to obtain the map, and then maintaining updates on the map as the environment changes. This paper addresses the issue of autonomous driving in an urban environment by investigating algorithms and an architecture to enable fully functional autonomous driving with limited information. An algorithm to autonomously navigate urban roadways with little to no reliance on an a priori map or GPS is developed. Localization is performed with an extended Kalman filter with odometry, compass, and sparse landmark measurement updates. Navigation is accomplished by a compass-based navigation control law. Key results from Monte Carlo studies show success rates of urban navigation under different environmental conditions. Experiments validate the simulated results and demonstrate that, for given test conditions, an expected range can be found for a given success rate.


title: Distributed Wildfire Surveillance with Autonomous Aircraft using Deep Reinforcement Learning

abstract: Teams of autonomous unmanned aircraft can be used to monitor wildfires, enabling firefighters to make informed decisions. However, controlling multiple autonomous fixed-wing aircraft to maximize forest fire coverage is a complex problem. The state space is high dimensional, the fire propagates stochastically, the sensor information is imperfect, and the aircraft must coordinate with each other to accomplish their mission. This work presents two deep reinforcement learning approaches for training decentralized controllers that accommodate the high dimensionality and uncertainty inherent in the problem. The first approach controls the aircraft using immediate observations of the individual aircraft. The second approach allows aircraft to collaborate on a map of the wildfire's state and maintain a time history of locations visited, which are used as inputs to the controller. Simulation results show that both approaches allow the aircraft to accurately track wildfire expansions and outperform an online receding horizon controller. Additional simulations demonstrate that the approach scales with different numbers of aircraft and generalizes to different wildfire shapes.


title: Neutron stars meet constraints from high and low energy nuclear physics

abstract: A novel equation of state used for analysis of the heavy ion collision experimental data is generalized to also describe the matter inside neutron stars. This approach differs from others by including an induced surface tension caused by the short range repulsion between particles accounted by presence of their hard cores. This new equation of state respects thermodynamic consistency and is free from the usual causality problems. This makes it applicable over a wide range of temperatures and baryon densities, including the ones inside neutron stars. Despite small number of parameters the present equation of state reproduces properties of the saturated nuclear matter and is consistent with the proton flow constraint as well as with astrophysical data obtained from neutron star observations. Accordingly, we found the sets of parameters that agree with the same ones obtained from the heavy ion collision data analysis.


title: Deep clustering: On the link between discriminative models and K-means

abstract: In the context of recent deep clustering studies, discriminative models dominate the literature and report the most competitive performances. These models learn a deep discriminative neural network classifier in which the labels are latent. Typically, they use multinomial logistic regression posteriors and parameter regularization, as is very common in supervised learning. It is generally acknowledged that discriminative objective functions (e.g., those based on the mutual information or the KL divergence) are more flexible than generative approaches (e.g., K-means) in the sense that they make fewer assumptions about the data distributions and, typically, yield much better unsupervised deep learning results. On the surface, several recent discriminative models may seem unrelated to K-means. This study shows that these models are, in fact, equivalent to K-means under mild conditions and common posterior models and parameter regularization. We prove that, for the commonly used logistic regression posteriors, maximizing the $L_2$ regularized mutual information via an approximate alternating direction method (ADM) is equivalent to a soft and regularized K-means loss. Our theoretical analysis not only connects directly several recent state-of-the-art discriminative models to K-means, but also leads to a new soft and regularized deep K-means algorithm, which yields competitive performance on several image clustering benchmarks.


title: Feature Selection using Stochastic Gates

abstract: Feature selection problems have been extensively studied for linear estimation, for instance, Lasso, but less emphasis has been placed on feature selection for non-linear functions. In this study, we propose a method for feature selection in high-dimensional non-linear function estimation problems. The new procedure is based on minimizing the $\ell_0$ norm of the vector of indicator variables that represent if a feature is selected or not. Our approach relies on the continuous relaxation of Bernoulli distributions, which allows our model to learn the parameters of the approximate Bernoulli distributions via gradient descent. This general framework simultaneously minimizes a loss function while selecting relevant features. Furthermore, we provide an information-theoretic justification of incorporating Bernoulli distribution into our approach and demonstrate the potential of the approach on synthetic and real-life applications.


title: The homotopy type of skeleta of the flag complex over a finite vector space

abstract: The aim of this paper is to give a (discrete) Morse theoretic proof of the fact that the $k$-th skeleton of the flag complex $\mathcal{F}$, associated to the lattice of subspaces of a finite dimensional vector space, is homotopy equivalent to a wedge of spheres of dimension $\min\{k,\dim(\mathcal{F})\}$. The tight control provided by Morse theoretic methods allows us to give an explicit formula for the number of spheres appearing in each of these wedge summands.


title: Data-dependent compression of random features for large-scale kernel approximation

abstract: Kernel methods offer the flexibility to learn complex relationships in modern, large data sets while enjoying strong theoretical guarantees on quality. Unfortunately, these methods typically require cubic running time in the data set size, a prohibitive cost in the large-data setting. Random feature maps (RFMs) and the Nystrom method both consider low-rank approximations to the kernel matrix as a potential solution. But, in order to achieve desirable theoretical guarantees, the former may require a prohibitively large number of features J+, and the latter may be prohibitively expensive for high-dimensional problems. We propose to combine the simplicity and generality of RFMs with a data-dependent feature selection scheme to achieve desirable theoretical approximation properties of Nystrom with just O(log J+) features. Our key insight is to begin with a large set of random features, then reduce them to a small number of weighted features in a data-dependent, computationally efficient way, while preserving the statistical guarantees of using the original large set of features. We demonstrate the efficacy of our method with theory and experiments--including on a data set with over 50 million observations. In particular, we show that our method achieves small kernel matrix approximation error and better test set accuracy with provably fewer random features than state-of-the-art methods.


title: Bird Species Classification using Transfer Learning with Multistage Training

abstract: Bird species classification has received more and more attention in the field of computer vision, for its promising applications in biology and environmental studies. Recognizing bird species is difficult due to the challenges of discriminative region localization and fine-grained feature learning. In this paper, we have introduced a Transfer learning based method with multistage training. We have used both Pre-Trained Mask-RCNN and an ensemble model consisting of Inception Nets (InceptionV3 & InceptionResNetV2 ) to get localization and species of the bird from the images respectively. Our final model achieves an F1 score of 0.5567 or 55.67 % on the dataset provided in CVIP 2018 Challenge.


title: Ultrahigh-Energy Cosmic Ray Composition from the Distribution of Arrival Directions

abstract: The sources of ultrahigh-energy cosmic rays (UHECRs) have been difficult to catch. It was recently pointed out that while sources of UHECR protons exhibit anisotropy patterns that become denser and compressed with rising energy, nucleus-emitting-sources give rise to a cepa stratis (onion-like) structure with layers that become more distant from the source position with rising energy. The peculiar shape of the hot spots from nucleus-accelerators is steered by the competition between energy loss during propagation and deflection on the Galactic magnetic field (GMF). Here, we run a full-blown simulation study to accurately characterize the deflections of UHECR nuclei in the GMF. We show that while the cepa stratis structure provides a global description of anisotropy patterns produced by UHECR nuclei en route to Earth, the hot spots are elongated depending on their location in the sky due to the regular structure of the GMF. We demonstrate that with a high-statistics sample at the high-energy-end of the spectrum, like the one to be collected by NASA's POEMMA mission, the energy dependence of the hot-spot contours could become a useful observable to identify the nuclear composition of UHECRs. This new method to determine the nature of the particle species is complementary to those using observables of extensive air showers, and therefore is unaffected by the large systematic uncertainties of hadronic interaction models.


title: Cycle Intersection Graphs and Minimum Decycling Sets of Even Graphs

abstract: We introduce the cycle intersection graph of a graph, an adaptation of the cycle graph of a graph, and use the structure of these graphs to prove an upper bound for the decycling number of all even graphs. This bound is shown to be significantly better when an even graph admits a cycle decomposition in which any two cycles intersect in at most one vertex. Links between the cycle rank of the cycle intersection graph of an even graph and the decycling number of the even graph itself are found. The problem of choosing an ideal cycle decomposition is addressed and is presented as an optimization problem over the space of cycle decompositions of even graphs.


title: Energy-resolved Photoconductivity Mapping in a Monolayer-bilayer WSe2 Lateral Heterostructure

abstract: Vertical and lateral heterostructures of van der Waals materials provide tremendous flexibility for band structure engineering. Since electronic bands are sensitively affected by defects, strain, and interlayer coupling, the edge and heterojunction of these two-dimensional (2D) systems may exhibit novel physical properties, which can be fully revealed only by spatially resolved probes. Here, we report the spatial mapping of photoconductivity in a monolayer-bilayer WSe2 lateral heterostructure under multiple excitation lasers. As the photon energy increases, the light-induced conductivity detected by microwave impedance microscopy first appears along the hetero-interface and bilayer edge, then along the monolayer edge, inside the bilayer area, and finally in the interior of the monolayer region. The sequential emergence of mobile carriers in different sections of the sample is consistent with the theoretical calculation of local energy gaps. Quantitative analysis of the microscopy and transport data also reveals the linear dependence of photoconductivity on the laser intensity and the influence of interlayer coupling on carrier recombination. Combining theoretical modeling, atomic scale imaging, mesoscale impedance microscopy, and device-level characterization, our work suggests an exciting perspective to control the intrinsic band-gap variation in 2D heterostructures down to the few-nanometer regime.


title: Extreme Classification in Log Memory

abstract: We present Merged-Averaged Classifiers via Hashing (MACH) for K-classification with ultra-large values of K. Compared to traditional one-vs-all classifiers that require O(Kd) memory and inference cost, MACH only need O(d log K) (d is dimensionality )memory while only requiring O(K log K + d log K) operation for inference. MACH is a generic K-classification algorithm, with provably theoretical guarantees, which requires O(log K) memory without any assumption on the relationship between classes. MACH uses universal hashing to reduce classification with a large number of classes to few independent classification tasks with small (constant) number of classes. We provide theoretical quantification of discriminability-memory tradeoff. With MACH we can train ODP dataset with 100,000 classes and 400,000 features on a single Titan X GPU, with the classification accuracy of 19.28%, which is the best-reported accuracy on this dataset. Before this work, the best performing baseline is a one-vs-all classifier that requires 40 billion parameters (160 GB model size) and achieves 9% accuracy. In contrast, MACH can achieve 9% accuracy with 480x reduction in the model size (of mere 0.3GB). With MACH, we also demonstrate complete training of fine-grained imagenet dataset (compressed size 104GB), with 21,000 classes, on a single GPU. To the best of our knowledge, this is the first work to demonstrate complete training of these extreme-class datasets on a single Titan X.


title: Efficient Trajectory Optimization for Robot Motion Planning

abstract: Motion planning for multi-jointed robots is challenging. Due to the inherent complexity of the problem, most existing works decompose motion planning as easier subproblems. However, because of the inconsistent performance metrics, only sub-optimal solution can be found by decomposition based approaches. This paper presents an optimal control based approach to address the path planning and trajectory planning subproblems simultaneously. Unlike similar works which either ignore robot dynamics or require long computation time, an efficient numerical method for trajectory optimization is presented in this paper for motion planning involving complicated robot dynamics. The efficiency and effectiveness of the proposed approach is shown by numerical results. Experimental results are used to show the feasibility of the presented planning algorithm.


title: The effects of invasive epibionts on crab-mussel communities: a theoretical approach to understand mussel population decline

abstract: Blue mussels (Mytilus edulis) are an important keystone species that have been declining in the Gulf of Maine. This could be attributed to a variety of complex factors such as indirect effects due to invasion by epibionts, which remains unexplored mathematically. Based on classical optimal foraging theory and anti-fouling defense mechanisms of mussels, we derive an ODE model for crab-mussel interactions in the presence of an invasive epibiont, Didemnum vexillum. The dynamical analysis leads to results on stability, global boundedness and bifurcations of the model. Next, via optimal control methods we predict various ecological outcomes. Our results have key implications for preserving mussel populations in the advent of invasion by non-native epibionts. In particular they help us understand the changing dynamics of local predator-prey communities, due to indirect effects that epibionts confer.


title: Notes on the Sasaki metric

abstract: We survey on the geometry of the tangent bundle of a Riemannian manifold, endowed with the classical metric established by S. Sasaki 60 years ago. Following the results of Sasaki, we try to write and deduce them by different means. Questions of vector fields, mainly those arising from the base, are related as invariants of the classical metric, contact and Hermitian structures. Attention is given to the natural notion of extension or complete lift of a vector field, from the base to the tangent manifold. Few results are original, but finally new equations of the mirror map are considered.


title: Geometric constructions over $\mathbb{C}$ and $\mathbb{F}_2$ for Quantum Information

abstract: In this review paper I present two geometric constructions of distinguished nature, one is over the field of complex numbers $\mathbb{C}$ and the other one is over the two elements field $\mathbb{F}_2$. Both constructions have been employed in the past fifteen years to describe two quantum paradoxes or two resources of quantum information: entanglement of pure multipartite systems on one side and contextuality on the other. Both geometric constructions are linked to representation of semi-simple Lie groups/algebras. To emphasize this aspect one explains on one hand how well-known results in representation theory allows one to see all the classification of entanglement classes of various tripartite quantum systems ($3$ qubits, $3$ fermions, $3$ bosonic qubits...) in a unified picture. On the other hand, one also shows how some weight diagrams of simple Lie groups are encapsulated in the geometry which deals with the commutation relations of the generalized $N$-Pauli group.

categories: quant-ph math-ph math.AG math.CO math.MP


title: Fair Division Minimizing Inequality

abstract: Behavioural economists have shown that people are often averse to inequality and will make choices to avoid unequal outcomes. In this paper, we consider how to allocate indivisible goods fairly so as to minimize inequality. We consider how this interacts with axiomatic properties such as envy-freeness, Pareto efficiency and strategy-proofness. We also consider the computational complexity of computing allocations minimizing inequality. Unfortunately, this is computationally intractable in general so we consider several tractable greedy online mechanisms that minimize inequality. Finally, we run experiments to explore the performance of these methods.


title: Inter-Scanner Harmonization of High Angular Resolution DW-MRI using Null Space Deep Learning

abstract: Diffusion-weighted magnetic resonance imaging (DW-MRI) allows for non-invasive imaging of the local fiber architecture of the human brain at a millimetric scale. Multiple classical approaches have been proposed to detect both single (e.g., tensors) and multiple (e.g., constrained spherical deconvolution, CSD) fiber population orientations per voxel. However, existing techniques generally exhibit low reproducibility across MRI scanners. Herein, we propose a data-driven tech-nique using a neural network design which exploits two categories of data. First, training data were acquired on three squirrel monkey brains using ex-vivo DW-MRI and histology of the brain. Second, repeated scans of human subjects were acquired on two different scanners to augment the learning of the network pro-posed. To use these data, we propose a new network architecture, the null space deep network (NSDN), to simultaneously learn on traditional observed/truth pairs (e.g., MRI-histology voxels) along with repeated observations without a known truth (e.g., scan-rescan MRI). The NSDN was tested on twenty percent of the histology voxels that were kept completely blind to the network. NSDN significantly improved absolute performance relative to histology by 3.87% over CSD and 1.42% over a recently proposed deep neural network approach. More-over, it improved reproducibility on the paired data by 21.19% over CSD and 10.09% over a recently proposed deep approach. Finally, NSDN improved gen-eralizability of the model to a third in vivo human scanner (which was not used in training) by 16.08% over CSD and 10.41% over a recently proposed deep learn-ing approach. This work suggests that data-driven approaches for local fiber re-construction are more reproducible, informative and precise and offers a novel, practical method for determining these models.


title: A Tale of Three Probabilistic Families: Discriminative, Descriptive and Generative Models

abstract: The pattern theory of Grenander is a mathematical framework where patterns are represented by probability models on random variables of algebraic structures. In this paper, we review three families of probability models, namely, the discriminative models, the descriptive models, and the generative models. A discriminative model is in the form of a classifier. It specifies the conditional probability of the class label given the input signal. A descriptive model specifies the probability distribution of the signal, based on an energy function defined on the signal. A generative model assumes that the signal is generated by some latent variables via a transformation. We shall review these models within a common framework and explore their connections. We shall also review the recent developments that take advantage of the high approximation capacities of deep neural networks.


title: A priori error estimates for the optimal control of the integral fractional Laplacian

abstract: We design and analyze solution techniques for a linear-quadratic optimal control problem involving the integral fractional Laplacian. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. We propose two strategies to discretize the fractional optimal control problem: a semidiscrete approach where the control is not discretized - the so-called variational discretization approach - and a fully discrete approach where the control variable is discretized with piecewise constant functions. Both schemes rely on the discretization of the state equation with the finite element space of continuous piecewise polynomials of degree one. We derive a priori error estimates for both solution techniques. We illustrate the theory with two-dimensional numerical tests.


title: Recent Developments in the Analysis of Galaxy Surveys

abstract: These are advanced lecture notes covering recent developments in the methodology used to analyse galaxy surveys. The focus is particularly on direct measurements of the galaxy power spectrum although I also discuss its Fourier transform, the correlation function for comparison. These 2-point statistics, under the assumption that the overdensity field has Gaussian statistics on large-scales, contain the majority of the cosmological signal available from the galaxy distribution. Recent developments in multipole measurements, dealing with systematics, convolving theoretical models with the survey window function, the approximation of covariance matrices, and weighting schemes for measuring evolution with redshift are considered. The focus is on analytic explanation of the issues involved rather than on recent analyses or simulation results.


title: Low Space-Charge Intensity Beams in UMER via Collimation and Solenoid Focusing

abstract: The University of Maryland Electron Ring (UMER) has operated traditionally in the regime of strong space-charge dominated beam transport, but small-current beams are desirable to significantly reduce the direct (incoherent) space-charge tune shift as well as the tune depression. This regime is of interest to model space-charge effects in large proton and ion rings similar to those used in nuclear physics and spallation neutron sources, and also for nonlinear dynamics studies of lattices inspired on the Integrable Optics Test Accelerator (IOTA). We review the definition of space-charge intensity, show a comparison of space-charge parameters in UMER and other machines, and discuss a simple method involving double collimation and solenoid focusing for varying the space-charge intensity of the beam injected into UMER.


title: Coherent diffusion of partial spatial coherence

abstract: Partially coherent light is abundant in many physical systems, and its propagation properties are well understood. Here we extend current theory of propagation of partially coherent light beams to the field of coherent diffusion. Based on a unique four-wave mixing scheme of electro-magnetically induced transparency, an optical speckle pattern is coupled to diffusing atoms in a warm vapor. The spatial coherence propagation properties of light speckles is studied experimentally under diffusion, and is compared to the familiar spatial coherence of speckles under diffraction. An analytic model explaining the results is presented, based on a diffusion analogue of the famous Van Cittert-Zernike theorem.


title: Edge-Insensitive Magnetism and Half Metallicity in Graphene Nanoribbons

abstract: Realizing magnetism in graphene/carbon nanostructures is a decade-long challenge. The magnetic edge state and half metallicity in zigzag graphene nanoribbons are particularly promising [Y.-W. Son, et al., Nature 444, 347 (2006)]. However, its experimental realization has been hindered by the stringent requirement of the mono-hydrogenated zigzag edge. Using first-principle calculations, we predict that free-carrier doping can overcome this challenge and realize ferromagnetism and half-metallicity in narrow graphene nanoribbons of general types of edge structures. This magnetism exists within the density range of gate-doping experiments (~1013 cm-2) and has large spin polarization energy up to 17 meV per carrier, which induces a Zeeman splitting equivalent to an external magnetic field of a few hundred Tesla. Finally, we trace the formation mechanics of this edge-insensitive magnetism to the quantum confinement of the electronic state near the band edge and reveal the scaling law of magnetism versus the ribbon width. Our findings suggest that combining doping with quantum confinement could be a general tool to realize transition-metal-free magnetism in light-element nanostructures.


title: A level-set approach for stochastic optimal control problems under controlled-loss constraints

abstract: We study a family of optimal control problems under a set of controlled-loss constraints holding at different deterministic dates. The characterization of the associated value function by a Hamilton-Jacobi-Bellman equation usually calls for additional strong assumptions on the dynamics of the processes involved and the set of constraints. To treat this problem in absence of those assumptions, we first convert it into a state-constrained stochastic target problem and then apply a level-set approach. With this approach, the state constraints can be managed through an exact penalization technique.


title: X-ray coherent diffraction imaging with an objective lens: towards 3D mapping of thick polycrystals

abstract: We report on a new x-ray imaging method, which combines the high spatial resolution of coherent diffraction imaging with the ability of dark field microscopy to map grains within thick polycrystalline specimens. An x-ray objective serves to isolate a grain and avoid overlap of diffraction spots. Iterative oversampling routines are used to reconstruct the shape and strain field within the grain from the far field intensity pattern. The limitation on resolution caused by the finite numerical aperture of the objective is overcome by the Fourier synthesis of several diffraction patterns. We demonstrate the method by an experimental study of a ~500 nm Pt grain for the two cases of a real and a virtual image plane. In the latter case the spatial resolution is 13 nm rms. Our results confirm that no information on the pupil function of the lens is required and that lens aberrations are not critical.


title: Scaling behavior in the heart rate variability characteristics with age

abstract: In this work we study the characteristics of the heart rate variability (HRV) as a function of age and gender. The analyzed data include previous results reported in the literature. The data obtained in this work expand the range of age studied until now revealing new behaviors not reported before. We analyze some measurements in the time domain,in the frequency domain and nonlinear measurements. We report scaling behaviors and abrupt changes in some measurements. There is also a progressive decrease in the dimensionality of the dynamic system governing the HRV, with the increase in age that is interpreted in terms ofautonomic regulation of cardiac activity.


title: The construction of averaged planetary motion theory by means of computer algebra system Piranha

abstract: The application of computer algebra system Piranha to the investigation of the planetary problem is described in this work. Piranha is an echeloned Poisson series processor authored by F. Biscani from Max Planck Institute for Astronomy in Heidelberg. Using Piranha the averaged semi-analytical motion theory of four-planetary system is constructed up to the second degree of planetary masses. In this work we use the algorithm of the Hamiltonian expansion into the Poisson series in only orbital elements without other variables. The motion equations are obtained analytically in time-averaged elements by Hori-Deprit method. Piranha showed high-performance of analytical manipulations. Different properties of obtained series are discussed. The numerical integration of the motion equations is performed by Everhart method for the Solar system's giant-planets and some exoplanetary systems.


title: Precise characterization of polarizing prisms by optical localization analysis

abstract: We utilize the localization analysis method to precisely determine the light beam positions with the spatial separation beyond the optical diffraction limit. By such a direct spatial measurement, the associated optical setup is built to characterize the beam-shear angle of polarizing prisms that refract the light propagation according to polarization. We use the Nomarski prisms to demonstrate our method. The statistical error is below 0.5 urad, which is two orders of magnitudes smaller than the tolerance value given by most prism manufacturers. Our method provides a stand-alone characterization of prisms with a high accuracy for many quantitative imaging technologies such as the orientation-free DIC microscope and the dual-focus fluorescence correlation spectroscopy.


title: Semigroup expansions for non-selfadjoint Schr\"odinger operators

abstract: A large time expansion for the propagator associated to a semiclassical non-selfadjoint magnetic Schr\"odinger operator is established, in terms of the low lying eigenvalues of the operator.


title: Convolutional Neural Networks and x-vector Embedding for DCASE2018 Acoustic Scene Classification Challenge

abstract: In this paper, the Brno University of Technology (BUT) team submissions for Task 1 (Acoustic Scene Classification, ASC) of the DCASE-2018 challenge are described. Also, the analysis of different methods on the leaderboard set is provided. The proposed approach is a fusion of two different Convolutional Neural Network (CNN) topologies. The first one is the common two-dimensional CNNs which is mainly used in image classification. The second one is a one-dimensional CNN for extracting fixed-length audio segment embeddings, so called x-vectors, which has also been used in speech processing, especially for speaker recognition. In addition to the different topologies, two types of features were tested: log mel-spectrogram and CQT features. Finally, the outputs of different systems are fused using a simple output averaging in the best performing system. Our submissions ranked third among 24 teams in the ASC sub-task A (task1a).


title: Survival prediction using ensemble tumor segmentation and transfer learning

abstract: Segmenting tumors and their subregions is a challenging task as demonstrated by the annual BraTS challenge. Moreover, predicting the survival of the patient using mainly imaging features, while being a desirable outcome to evaluate the treatment of the patient, it is also a difficult task. In this paper, we present a cascaded pipeline to segment the tumor and its subregions and then we use these results and other clinical features together with image features coming from a pretrained VGG-16 network to predict the survival of the patient. Preliminary results with the training and validation dataset show a promising start in terms of segmentation, while the prediction values could be improved with further testing on the feature extraction part of the network.


title: Self-Folding Metasheets: The Optimal Pattern of Strain of Miura-Ori Folded State

abstract: Self-folding origami has emerged as a tool to make functional objects in material science. The common idea is to pattern a sheet with creases and activate them to have the object fold spontaneously into a desired configuration. This article shows that collinear quadrilateral metasheets are able to fold into the Miura-Ori configuration, if we only impose strain on part of their creases. In this study, we define and determine the optimal pattern of strain (OPS) on a collinear quadrilateral metasheet, that is the pattern of minimum "functional" creases with which the self-folding metasheet can fold into Miura-Ori state stably. By comparing the energy evolution along the folding pathway of each possible folded state under OPS, we conclude that the energy predominance of the desired Miura-Ori pathway during the initial period of time accounts for why the OPS works. Furthermore, we measure the projected force of the OPS on the intial flat metasheet and give insights on how to determine the OPS using only local information of the initial flat state.


title: Current Trends and Future Research Directions for Interactive Music

abstract: In this review, it is explained and compared different software and formalisms used in music interaction: sequencers, computer-assisted improvisation, meta- instruments, score-following, asynchronous dataflow languages, synchronous dataflow languages, process calculi, temporal constraints and interactive scores. Formal approaches have the advantage of providing rigorous semantics of the behavior of the model and proving correctness during execution. The main disadvantage of formal approaches is lack of commercial tools.


title: Solitons in Nonlinear Systems and Eigen-states in Quantum Wells

abstract: We study the relations between solitons of nonlinear Schr\"{o}dinger equation described systems and eigen-states of linear Schr\"{o}dinger equation with some quantum wells. Many different non-degenerated solitons are re-derived from the eigen-states in the quantum wells. We show that the vector solitons for coupled system with attractive interactions correspond to the identical eigen-states with the ones of coupled systems with repulsive interactions. The energy eigenvalues of them seem to be different, but they can be reduced to identical ones in the same quantum wells. The non-degenerated solitons for multi-component systems can be used to construct much abundant degenerated solitons in more components coupled cases. On the other hand, we demonstrate soliton solutions in nonlinear systems can be also used to solve the eigen-problems of quantum wells. As an example, we present eigenvalue and eigen-state in a complicated quantum well for which the Hamiltonian belongs to the non-Hermitian Hamiltonian having Parity-Time symmetry. We further present the ground state and the first exited state in an asymmetric quantum double-well from asymmetric solitons. Based on these results, we expect that many nonlinear physical systems can be used to observe the quantum states evolution of quantum wells, such as water wave tank, nonlinear fiber, Bose-Einstein condensate, and even plasma, although some of them are classical physical systems. These relations provide another different way to understand the stability of solitons in nonlinear Schr\"{o}dinger equation described systems, in contrast to the balance between dispersion and nonlinearity.

categories: nlin.PS hep-th math-ph math.MP nlin.SI quant-ph


title: Time discretization schemes for hyperbolic systems on networks by $\epsilon$-expansion

abstract: We consider partial differential equations on networks with a small parameter $\epsilon$, which are hyperbolic for $\epsilon>0$ and parabolic for $\epsilon=0$. With a combination of an $\epsilon$-expansion and Runge-Kutta schemes for constrained systems of parabolic type, we derive a new class of time discretization schemes for hyperbolic systems on networks, which are constrained due to interconnection conditions. For the analysis we consider the coupled system equations as partial differential-algebraic equations based on the variational formulation of the problem. We discuss well-posedness of the resulting systems and estimate the error caused by the $\epsilon$-expansion.


title: Structured decomposition for reversible Boolean functions

abstract: Reversible Boolean function is a one-to-one function which maps $n$-bit input to $n$-bit output. Reversible logic synthesis has been widely studied due to its relationship with low-energy computation as well as quantum computation. In this work, we give a structured decomposition for even reversible Boolean functions (RBF). Specifically, for $n\geq 6$, any even $n$-bit RBF can be decomposed to $7$ blocks of $(n-1)$-bit RBF, where $7$ is a constant independent of $n$; and the positions of those blocks have large degree of freedom. Moreover, if the $(n-1)$-bit RBFs are required to be even as well, we show for $n\geq 10$, $n$-bit RBF can be decomposed to $10$ even $(n-1)$-bit RBFs. For simplicity, we say our decomposition has block depth $7$ and even block depth $10$. Our result improves Selinger's work in block depth model, by reducing the constant from $9$ to $7$; and from $13$ to $10$ when the blocks are limited to be even. We emphasize that our setting is a bit different from Selinger's. In Selinger's constructive proof, each block is one of two specific positions and thus the decomposition has an alternating structure. We relax this restriction and allow each block to act on arbitrary $(n-1)$ bits. This relaxation keeps the block structure and provides more candidates when choosing positions of blocks.


title: Is Neuron Made from Mathematics?

abstract: This paper is to derive a mathematical model for neuron by imposing only a principle of symmetry that two modelers must come up with the same model when one is approaching the problem by modeling the conductances of ion channels and the other by the channel resistances.


title: Fully integrative data analysis of NMR metabolic fingerprints with comprehensive patient data: a case report based on the German Chronic Kidney Disease (GCKD) study

abstract: Omics data facilitate the gain of novel insights into the pathophysiology of diseases and, consequently, their diagnosis, treatment, and prevention. To that end, it is necessary to integrate omics data with other data types such as clinical, phenotypic, and demographic parameters of categorical or continuous nature. Here, we exemplify this data integration issue for a study on chronic kidney disease (CKD), where complex clinical and demographic parameters were assessed together with one-dimensional (1D) 1H NMR metabolic fingerprints. Routine analysis screens for associations of single metabolic features with clinical parameters, which requires confounding variables typically chosen by expert knowledge to be taken into account. This knowledge can be incomplete or unavailable. The results of this article are manifold. We introduce a framework for data integration that intrinsically adjusts for confounding variables. We give its mathematical and algorithmic foundation, provide a state-of-the-art implementation, and give several sanity checks. In particular, we show that the discovered associations remain significant after variable adjustment based on expert knowledge. In contrast, we illustrate that the discovery of associations in routine analysis can be biased by incorrect or incomplete expert knowledge in univariate screening approaches. Finally, we exemplify how our data integration approach reveals important associations between CKD comorbidities and metabolites. Moreover, we evaluate the predictive performance of the estimated models on independent validation data and contrast the results with a naive screening approach.


title: Chebyshev Interpolation for Function in 1D

abstract: This research is concerned with finding the roots of a function in an interval using Chebyshev Interpolation. Numerical results of Chebyshev Interpolation are presented to show that this is a powerful way to simultaneously calculate all the roots in an interval.


title: The Ricci-Bourguignon flow on Heisenberg and quaternion Lie groups

abstract: In this paper, we study the Ricci-Bourguignon flow on higher dimensional classical Heisenberg nilpotent Lie groups and construct a solution of this flow on Heisenberg and quaternion nilpotent Lie groups. In the end, we investigate the deformation of spectrum and length spectrum on compact nilmanifolds obtained of Heisenberg and quaternion nilpotent Lie groups.


title: The seed order

abstract: We study various orders on countably complete ultrafilters on ordinals that coincide and are wellorders under a hypothesis called the Ultrapower Axiom. Our main focus is on the relationship between the Ultrapower Axiom and the linearity of these orders.


title: Warped Hypertime Representations for Long-term Autonomy of Mobile Robots

abstract: This paper presents a novel method for introducing time into discrete and continuous spatial representations used in mobile robotics, by modelling long-term, pseudo-periodic variations caused by human activities. Unlike previous approaches, the proposed method does not treat time and space separately, and its continuous nature respects both the temporal and spatial continuity of the modeled phenomena. The method extends the given spatial model with a set of wrapped dimensions that represent the periodicities of observed changes. By performing clustering over this extended representation, we obtain a model that allows us to predict future states of both discrete and continuous spatial representations. We apply the proposed algorithm to several long-term datasets and show that the method enables a robot to predict future states of repre- sentations with different dimensions. The experiments further show that the method achieves more accurate predictions than the previous state of the art.


title: A maximum-mean-discrepancy goodness-of-fit test for censored data

abstract: We introduce a kernel-based goodness-of-fit test for censored data, where observations may be missing in random time intervals: a common occurrence in clinical trials and industrial life-testing. The test statistic is straightforward to compute, as is the test threshold, and we establish consistency under the null. Unlike earlier approaches such as the Log-rank test, we make no assumptions as to how the data distribution might differ from the null, and our test has power against a very rich class of alternatives. In experiments, our test outperforms competing approaches for periodic and Weibull hazard functions (where risks are time dependent), and does not show the failure modes of tests that rely on user-defined features. Moreover, in cases where classical tests are provably most powerful, our test performs almost as well, while being more general.


title: The asymptotic distance between an ultraflat unimodular polynomial and its conjugate reciprocal

abstract: Let $${\mathcal K}_n := \left\{p_n: p_n(z) = \sum_{k=0}^n{a_k z^k}, \enspace a_k \in {\mathbb C}\,,\enspace |a_k| = 1 \right\}\,.$$ A sequence $(P_n)$ of polynomials $P_n \in {\mathcal K}_n$ is called ultraflat if $(n + 1)^{-1/2}|P_n(e^{it})|$ converge to $1$ uniformly in $t \in {\mathbb R}$. In this paper we prove that $$\frac{1}{2\pi} \int_0^{2\pi}{\left| (P_n - P_n^*)(e^{it}) \right|^q \, dt} \sim \frac{{2}^q \Gamma \left(\frac{q+1}{2} \right)}{\Gamma \left(\frac q2 + 1 \right) \sqrt{\pi}} \,\, n^{q/2}$$ for every ultraflat sequence $(P_n)$ of polynomials $P_n \in {\mathcal K}_n$ and for every $q \in (0,\infty)$, where $P_n^*$ is the conjugate reciprocal polynomial associated with $P_n$, $\Gamma$ is the usual gamma function, and the $\sim$ symbol means that the ratio of the left and right hand sides converges to $1$ as $n \rightarrow \infty$. Another highlight of the paper states that $$\frac{1}{2\pi}\int_0^{2\pi}{\left| (P_n^\prime - P_n^{*\prime})(e^{it}) \right|^2 \, dt} \sim \frac{2n^3}{3}$$ for every ultraflat sequence $(P_n)$ of polynomials $P_n \in {\mathcal K}_n$. We prove a few other new results and reprove some interesting old results as well.


title: The Rudin-Frolik order and the Ultrapower Axiom

abstract: We study the structure of the Rudin-Frolik order on countably complete ultrafilters under the assumption that this order is directed. This assumption, called the Ultrapower Axiom, holds in all known canonical inner models. It turns out that assuming the Ultrapower Axiom, much more about the Rudin-Frolik order can be determined. Our main theorem is that under the Ultrapower Axiom, a countably complete ultrafilter has at most finitely many predecessors in the Rudin-Frolik order. In other words, any wellfounded ultrapower (of the universe) is the ultrapower of at most finitely many ultrapowers.


title: Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

abstract: Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems, but need improvements when the problems become harder in the case when (1) the geometry becomes complex (from a few branches to a full pulmonary artery), (2) the model becomes more complex (from fluid-only calculation to coupled fluid-structure interaction calculation), (3) both the fluid and wall models become highly nonlinear, and (4) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from CT (Computed Tomography) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10,000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.


title: The internal relation

abstract: This paper explores various generalizations of the Mitchell order focusing mostly on a generalization called the internal relation. The internal relation lacks the implicit strength requirement in the definition of the Mitchell order, and therefore can fail to be wellfounded. We establish some constraints on the illfoundedness of the internal relation, which leads to a proof of a conjecture of Steel regarding rank-to-rank cardinals.


title: Equivalences between localisations of categories provided by replacements

abstract: We give a characterisation of functors whose induced functor on the level of localisations is an equivalence and where the isomorphism inverse is induced by some kind of replacements such as projective resolutions or cofibrant replacements.


title: On a theorem of Scholze-Weinstein

abstract: Let G be the Tate module of a p-divisble group H over a perfect field k of characteristic p. A theorem of Scholze-Weinstein describes G (and therefore H itself) in terms of the Dieudonne module of H; more precisely, it describes G(C) for "good" semiperfect k-algebras C (which is enough to reconstruct G). In these notes we give a self-contained proof of this theorem and explain the relation with the classical descriptions of the Dieudonne functor from Dieudonne modules to p-divisible groups.


title: Towards geometric Satake correspondence for Kac-Moody algebras -- Cherkis bow varieties and affine Lie algebras of type $A$

abstract: We give a provisional construction of the Kac-Moody Lie algebra module structure on the hyperbolic restriction of the intersection cohomology complex of the Coulomb branch of a framed quiver gauge theory, as a refinement of the conjectural geometric Satake correspondence for Kac-Moody algebras proposed in an earlier paper with Braverman, Finkelberg in 2019. This construction assumes several geometric properties of the Coulomb branch under the torus action. These properties are checked in affine type A, via the identification of the Coulomb branch with a Cherkis bow variety established in a joint work with Takayama.


title: Nucleon electromagnetic form factors at high $Q^2$ from Wilson-clover fermions

abstract: We present results on the nucleon electromagnetic form factors from Lattice QCD at momentum transfer up to about $12$~GeV$^2$. We analyze two gauge ensembles with the Wilson-clover fermion action, a lattice spacing of $a\approx 0.09$~fm and pion masses $m_\pi\approx 170$~MeV and $m_\pi\approx 280$~MeV. In our analysis we employ momentum smearing as well as a set of techniques to investigate excited state effects. Good agreement with experiment and phenomenology is found for the ratios $G_E/G_M$ and $F_2/F_1$, whereas discrepancies are observed for the individual form factors $F_1$ and $F_2$. We discuss various systematics that may affect our calculation.


title: RCF-defense for quantum randomness

abstract: Authors develop new nonparametric methods for verification and monitoring for quantum randomness based on the ranged correlation function (RCF) and a sequence of the ranged amplitudes (SRA). We carried out RCF-analysis of different topology subsamples from raw data of the prototype of a quantum random number generator on homodyne detection. It is shown that in the real system there are weak local regression relations, for which it is possible to introduce a robust criterion of significance, and also precise SRA-identification of long samples statistics is made. The obtained results extend the traditional entropy methods of the useful randomness analysis and open the way for creation of new strict quality quantum standards and defense for physical random numbers generators.


title: Epidemics, production and savings. Why saving is important?

abstract: In this work we combing models of disease dynamics and economic production, and we show the potential implications of this for demonstrating the importance of savings for buffering an economy during the period of an epidemic. Finding an explicit function that relates poverty and the production of a community is an almost impossible task because of the number of variables and parameters that should be taken into account. However, studying the dynamics of an endemic disease in a region that affects its population, and therefore its ability to work, is an honest approach to understanding this function. We propose a model, perhaps the simplest, that couples two dynamics, the dynamics of an endemic disease and the dynamics of a closed economy of products and goods that the community produces in the epidemic period. Some of the results of this study are expected and known in the literature but some others are not. We highlight three of them: the interdependence that exists between health and the product of that economy. The majority of the known results only show the dependence in a single direction. Another of the most important goals in this work is to show the existence of a poverty Malthusian trap in economies with low levels of capital investment due to pandemic disease. And another is to show that there is an optimal level of savings that maximizes the \emph{per capita} product of the economy, which is Machiavellian, but has important implications for market efficiency. The model tells us why saving is important in an aggregate economy despite being affected by an epidemic.


title: Decipherment of Historical Manuscript Images

abstract: European libraries and archives are filled with enciphered manuscripts from the early modern period. These include military and diplomatic correspondence, records of secret societies, private letters, and so on. Although they are enciphered with classical cryptographic algorithms, their contents are unavailable to working historians. We therefore attack the problem of automatically converting cipher manuscript images into plaintext. We develop unsupervised models for character segmentation, character-image clustering, and decipherment of cluster sequences. We experiment with both pipelined and joint models, and we give empirical results for multiple ciphers.


title: Polar Codes with exponentially small error at finite block length

abstract: We show that the entire class of polar codes (up to a natural necessary condition) converge to capacity at block lengths polynomial in the gap to capacity, while simultaneously achieving failure probabilities that are exponentially small in the block length (i.e., decoding fails with probability $\exp(-N^{\Omega(1)})$ for codes of length $N$). Previously this combination was known only for one specific family within the class of polar codes, whereas we establish this whenever the polar code exhibits a condition necessary for any polarization. Our results adapt and strengthen a local analysis of polar codes due to the authors with Nakkiran and Rudra [Proc. STOC 2018]. Their analysis related the time-local behavior of a martingale to its global convergence, and this allowed them to prove that the broad class of polar codes converge to capacity at polynomial block lengths. Their analysis easily adapts to show exponentially small failure probabilities, provided the associated martingale, the ``Arikan martingale'', exhibits a corresponding strong local effect. The main contribution of this work is a much stronger local analysis of the Arikan martingale. This leads to the general result claimed above. In addition to our general result, we also show, for the first time, polar codes that achieve failure probability $\exp(-N^{\beta})$ for any $\beta < 1$ while converging to capacity at block length polynomial in the gap to capacity. Finally we also show that the ``local'' approach can be combined with any analysis of failure probability of an arbitrary polar code to get essentially the same failure probability while achieving block length polynomial in the gap to capacity.


title: Stochasticization of Solutions to the Yang-Baxter Equation

abstract: In this paper we introduce a procedure that, given a solution to the Yang-Baxter equation as input, produces a stochastic (or Markovian) solution to (a possibly dynamical version of) the Yang-Baxter equation. We then apply this "stochasticization procedure" to obtain three new, stochastic solutions to several different forms of the Yang-Baxter equation. The first is a stochastic, elliptic solution to the dynamical Yang-Baxter equation; the second is a stochastic, higher rank solution to the dynamical Yang-Baxter equation; and the third is a stochastic solution to a dynamical variant of the tetrahedron equation.


title: Relative Error of Scaled Poisson Approximation via Stein's Method

abstract: We study the accuracy of a scaled Poisson approximation to the weighted sum of independent Poisson random variables, focusing on in particular the relative error of the tail distribution. A bound on the relative approximation error is established using a modified Stein-Chen method.


title: Detection and Mitigation of Biasing Attacks on Distributed Estimation Networks

abstract: The paper considers a problem of detecting and mitigating biasing attacks on networks of state observers targeting cooperative state estimation algorithms. The problem is cast within the recently developed framework of distributed estimation utilizing the vector dissipativity approach. The paper shows that a network of distributed observers can be endowed with an additional attack detection layer capable of detecting biasing attacks and correcting their effect on estimates produced by the network. An example is provided to illustrate the performance of the proposed distributed attack detector.


title: On Tracking the Physicality of Wi-Fi: A Subspace Approach

abstract: Wi-Fi channel state information (CSI) has emerged as a plausible modality for sensing different human activities as a function of modulations in the wireless signal that travels between wireless devices. Until now, most research has taken a statistical approach and/or purpose-built inference pipeline. Although interesting, these approaches struggle to sustain sensing performances beyond experimental conditions. As such, the full potential of CSI as a general-purpose sensing modality is yet to be realised. We argue a universal approach with well-grounded formalisation is necessary to characterise the relationship between wireless channel modulations (spatial and temporal) and human movement. To this end, we present a formalism for quantifying the changing part of the wireless signal modulated by human motion. Grounded in this formalisation, we then present a new subspace tracking technique to describe the channel statistics in an interpretable way, which succinctly contains the human modulated part of the channel. We characterise the signal and noise subspaces for the case of uncontrolled human movement, and show that these subspaces are dynamic. Our results demonstrate that proposed channel statistics alone can robustly reproduce state-of-the-art application-specific feature engineering baseline, however, across multiple usage scenarios. We expect, our universal channel statistics will yield an effective general-purpose featurisation of wireless channel measurements and will uncover opportunities for applying CSI for a variety of human sensing applications in a robust way.


title: Batch Active Preference-Based Learning of Reward Functions

abstract: Data generation and labeling are usually an expensive part of learning for robotics. While active learning methods are commonly used to tackle the former problem, preference-based learning is a concept that attempts to solve the latter by querying users with preference questions. In this paper, we will develop a new algorithm, batch active preference-based learning, that enables efficient learning of reward functions using as few data samples as possible while still having short query generation times. We introduce several approximations to the batch active learning problem, and provide theoretical guarantees for the convergence of our algorithms. Finally, we present our experimental results for a variety of robotics tasks in simulation. Our results suggest that our batch active learning algorithm requires only a few queries that are computed in a short amount of time. We then showcase our algorithm in a study to learn human users' preferences.


title: Multi-Institutional Deep Learning Modeling Without Sharing Patient Data: A Feasibility Study on Brain Tumor Segmentation

abstract: Deep learning models for semantic segmentation of images require large amounts of data. In the medical imaging domain, acquiring sufficient data is a significant challenge. Labeling medical image data requires expert knowledge. Collaboration between institutions could address this challenge, but sharing medical data to a centralized location faces various legal, privacy, technical, and data-ownership challenges, especially among international institutions. In this study, we introduce the first use of federated learning for multi-institutional collaboration, enabling deep learning modeling without sharing patient data. Our quantitative results demonstrate that the performance of federated semantic segmentation models (Dice=0.852) on multimodal brain scans is similar to that of models trained by sharing data (Dice=0.862). We compare federated learning with two alternative collaborative learning methods and find that they fail to match the performance of federated learning.


title: Interband electron pairing for superconductivity from the breakdown of the Born-Oppenheimer approximation

abstract: The origin of interband electron pairing responsible for enhancing superconductivity and the factors controlling its strength were examined. We show that the interband electron pairing is a natural consequence of breaking down the Born-Oppenheimer approximation during the electron-phonon interactions. Its strength is determined by the pair-state excitations around the Fermi surfaces that take place to form a superconducting state. Fermi surfaces favorable for the pairing were found and its implications were discussed.


title: Extension groups between atoms in abelian categories

abstract: We introduce the extension groups between atoms in an abelian category. For a locally noetherian Grothendieck category, the localizing subcategories closed under injective envelopes are characterized in terms of those extension groups. We also introduce the virtual duals of the extension groups between atoms to measure the global dimension of the category. A new topological property of atom spectra is revealed and it is used to relate the projective dimensions of atoms with the Krull-Gabriel dimensions. As a byproduct of the topological observation, we show that there exists a spectral space that is not homeomorphic to the atom spectrum of any abelian category.


title: Galactic Panspermia

abstract: We present an analytic model to estimate the total number of rocky or icy objects that could be captured by planetary systems within the Milky Way galaxy and result in panspermia should they harbor life. We estimate the capture rate of objects ejected from planetary systems over the entire phase space as well as time. Our final expression for the capture rate depends upon the velocity dispersion as well as the characteristic biological survival time and the size of the captured object. We further take into account the number of stars that an interstellar object traverses, as well as the scale height and length of the Milky Way's disk. The likelihood of Galactic panspermia is strongly dependent upon the survival lifetime of the putative organisms as well as the velocity of the transporter. Velocities between $10-100$ km s$^{-1}$ result in the highest probabilities. However, given large enough survival lifetimes, even hypervelocity objects traveling at over 1000 km s$^{-1}$ have a significant chance of capture, thereby increasing the likelihood of panspermia. Thus, we show that panspermia is not exclusively relegated to solar-system sized scales, and the entire Milky Way could potentially be exchanging biotic components across vast distances.


title: A Simple Java Code Generator for ACL2 Based on a Deep Embedding of ACL2 in Java

abstract: AIJ (ACL2 In Java) is a deep embedding in Java of an executable, side-effect-free, non-stobj-accessing subset of the ACL2 language without guards. ATJ (ACL2 To Java) is a simple Java code generator that turns ACL2 functions into AIJ representations that are evaluated by the AIJ interpreter. AIJ and ATJ enable possibly verified ACL2 code to run as, and interoperate with, Java code, without much of the ACL2 framework or any of the Lisp runtime. The current speed of the resulting Java code may be adequate to some applications.


title: Formalising Filesystems in the ACL2 Theorem Prover: an Application to FAT32

abstract: In this work, we present an approach towards constructing executable specifications of existing filesystems and verifying their functional properties in a theorem proving environment. We detail an application of this approach to the FAT32 filesystem. We also detail the methodology used to build up this type of executable specification through a series of models which incrementally add features of the target filesystem. This methodology has the benefit of allowing the verification effort to start from simple models which encapsulate features common to many filesystems and which are thus suitable for reuse.


title: Trapezoidal Generalization over Linear Constraints

abstract: We are developing a model-based fuzzing framework that employs mathematical models of system behavior to guide the fuzzing process. Whereas traditional fuzzing frameworks generate tests randomly, a model-based framework can deduce tests from a behavioral model using a constraint solver. Because the state space being explored by the fuzzer is often large, the rapid generation of test vectors is crucial. The need to generate tests quickly, however, is antithetical to the use of a constraint solver. Our solution to this problem is to use the constraint solver to generate an initial solution, to generalize that solution relative to the system model, and then to perform rapid, repeated, randomized sampling of the generalized solution space to generate fuzzing tests. Crucial to the success of this endeavor is a generalization procedure with reasonable size and performance costs that produces generalized solution spaces that can be sampled efficiently. This paper describes a generalization technique for logical formulae expressed in terms of Boolean combinations of linear constraints that meets the unique performance requirements of model-based fuzzing. The technique represents generalizations using trapezoidal solution sets consisting of ordered, hierarchical conjunctions of linear constraints that are more expressive than simple intervals but are more efficient to manipulate and sample than generic polytopes. Supporting materials contain an ACL2 proof that verifies the correctness of a low-level implementation of the generalization algorithm against a specification of generalization correctness. Finally a post-processing procedure is described that results in a restricted trapezoidal solution that can be sampled (solved) rapidly and efficiently without backtracking, even for integer domains. While informal correctness arguments are provided, a formal proof of the correctness of the restriction algorithm remains as future work.


title: Incremental SAT Library Integration Using Abstract Stobjs

abstract: We describe an effort to soundly use off-the-shelf incremental SAT solvers within ACL2 by modeling the behavior of a SAT solver library as an abstract stobj. The interface allows ACL2 programs to use incremental SAT solvers, and the abstract stobj model allows us to reason about the behavior of an incremental SAT library so as to show that algorithms implemented using it are correct, as long as the library is bug-free.


title: Using ACL2 in the Design of Efficient, Verifiable Data Structures for High-Assurance Systems

abstract: Verification of algorithms and data structures utilized in modern autonomous and semi-autonomous vehicles for land, sea, air, and space presents a significant challenge. Autonomy algorithms, e.g., route planning, pattern matching, and inference, are based on complex data structures such as directed graphs and algebraic data types. Proof techniques for these data structures exist, but are oriented to unbounded, functional realizations, which are not typically efficient in either space or time. Autonomous systems designers, on the other hand, generally limit the space and time allocations for any given function, and require that algorithms deliver results within a finite time, or suffer a watchdog timeout. Furthermore, high-assurance design rules frown on dynamic memory allocation, preferring simple array-based data structure implementations. In order to provide efficient implementations of high-level data structures used in autonomous systems with the high assurance needed for accreditation, we have developed a verifying compilation technique that supports the "natural" functional proof style, but yet applies to more efficient data structure implementations. Our toolchain features code generation to mainstream programming languages, as well as GPU-based and hardware-based realizations. We base the Intermediate Verification Language for our toolchain upon higher-order logic; however, we have used ACL2 to develop our efficient yet verifiable data structure design. ACL2 is particularly well-suited for this work, with its sophisticated libraries for reasoning about aggregate data structures of arbitrary size, efficient execution of formal specifications, as well as its support for "single-threaded objects" -- functional datatypes with imperative "under the hood" implementations. In this paper, we detail our high-assurance data structure design approach, including examples in ACL2 of common algebraic data types implemented using this design approach, proofs of correctness for those data types carried out in ACL2, as well as sample ACL2 implementations of relevant algorithms utilizing these efficient, high-assurance data structures.


title: Adding 32-bit Mode to the ACL2 Model of the x86 ISA

abstract: The ACL2 model of the x86 Instruction Set Architecture was built for the 64-bit mode of operation of the processor. This paper reports on our work to extend the model with support for 32-bit mode, recounting the salient aspects of this activity and identifying the ones that required the most work.


title: The Fundamental Theorem of Algebra in ACL2

abstract: We report on a verification of the Fundamental Theorem of Algebra in ACL2(r). The proof consists of four parts. First, continuity for both complex-valued and real-valued functions of complex numbers is defined, and it is shown that continuous functions from the complex to the real numbers achieve a minimum value over a closed square region. An important case of continuous real-valued, complex functions results from taking the traditional complex norm of a continuous complex function. We think of these continuous functions as having only one (complex) argument, but in ACL2(r) they appear as functions of two arguments. The extra argument is a "context", which is uninterpreted. For example, it could be other arguments that are held fixed, as in an exponential function which has a base and an exponent, either of which could be held fixed. Second, it is shown that complex polynomials are continuous, so the norm of a complex polynomial is a continuous real-valued function and it achieves its minimum over an arbitrary square region centered at the origin. This part of the proof benefits from the introduction of the "context" argument, and it illustrates an innovation that simplifies the proofs of classical properties with unbound parameters. Third, we derive lower and upper bounds on the norm of non-constant polynomials for inputs that are sufficiently far away from the origin. This means that a sufficiently large square can be found to guarantee that it contains the global minimum of the norm of the polynomial. Fourth, it is shown that if a given number is not a root of a non-constant polynomial, then it cannot be the global minimum. Finally, these results are combined to show that the global minimum must be a root of the polynomial. This result is part of a larger effort in the formalization of complex polynomials in ACL2(r).


title: Real Vector Spaces and the Cauchy-Schwarz Inequality in ACL2(r)

abstract: We present a mechanical proof of the Cauchy-Schwarz inequality in ACL2(r) and a formalisation of the necessary mathematics to undertake such a proof. This includes the formalisation of $\mathbb{R}^n$ as an inner product space. We also provide an application of Cauchy-Schwarz by formalising $\mathbb R^n$ as a metric space and exhibiting continuity for some simple functions $\mathbb R^n\to\mathbb R$. The Cauchy-Schwarz inequality relates the magnitude of a vector to its projection (or inner product) with another: \[|\langle u,v\rangle| \leq \|u\| \|v\|\] with equality iff the vectors are linearly dependent. It finds frequent use in many branches of mathematics including linear algebra, real analysis, functional analysis, probability, etc. Indeed, the inequality is considered to be among "The Hundred Greatest Theorems" and is listed in the "Formalizing 100 Theorems" project. To the best of our knowledge, our formalisation is the first published proof using ACL2(r) or any other first-order theorem prover.


title: Convex Functions in ACL2(r)

abstract: This paper builds upon our prior formalisation of R^n in ACL2(r) by presenting a set of theorems for reasoning about convex functions. This is a demonstration of the higher-dimensional analytical reasoning possible in our metric space formalisation of R^n. Among the introduced theorems is a set of equivalent conditions for convex functions with Lipschitz continuous gradients from Yurii Nesterov's classic text on convex optimisation. To the best of our knowledge a full proof of the theorem has yet to be published in a single piece of literature. We also explore "proof engineering" issues, such as how to state Nesterov's theorem in a manner that is both clear and useful.


title: Smtlink 2.0

abstract: Smtlink is an extension of ACL2 with Satisfiability Modulo Theories (SMT) solvers. We presented an earlier version at ACL2'2015. Smtlink 2.0 makes major improvements over the initial version with respect to soundness, extensibility, ease-of-use, and the range of types and associated theory-solvers supported. Most theorems that one would want to prove using an SMT solver must first be translated to use only the primitive operations supported by the SMT solver -- this translation includes function expansion and type inference. Smtlink 2.0 performs this translation using a sequence of steps performed by verified clause processors and computed hints. These steps are ensured to be sound. The final transliteration from ACL2 to Z3's Python interface requires a trusted clause processor. This is a great improvement in soundness and extensibility over the original Smtlink which was implemented as a single, monolithic, trusted clause processor. Smtlink 2.0 provides support for FTY defprod, deflist, defalist, and defoption types by using Z3's arrays and user-defined data types. We have identified common usage patterns and simplified the configuration and hint information needed to use Smtlink.


title: Hint Orchestration Using ACL2's Simplifier

abstract: This paper describes a strategy for providing hints during an ACL2 proof, implemented in a utility called use-termhint. An extra literal is added to the goal clause and simplified along with the rest of the goal until it is stable under simplification, after which the simplified literal is examined and a hint extracted from it. This simple technique supports some commonly desirable yet elusive features. It supports providing different hints to different cases of a case split, as well as binding variables so as to avoid repeating multiply referenced subterms. Since terms used in these hints are simplified in the same way as the rest of the goal, this strategy is also more robust against changes in the rewriting normal form than hints in which terms from the goal are written out explicitly.


title: Reply to Comment on "Gravitational Waves From Ultra- Short Period Exoplanets"

abstract: Wong et al. (2018) recently performed an encouraging criticism to our paper "Gravitational waves from ultra-short period exoplanets" (Cunha, Silva, Lima 2018) exploring the potentialities of a subset of exoplanets with extremely short periods (less than 80 min) as a possible scientific target to the planned space-based LISA observatory. Here we call attention to some subtleties and limitations underlying the basic criticism which in our view were not properly stressed in their comment. Particularly, simple estimates show that a sphere encircling the Earth with a radius of 250 pc may accommodate a population $ \sim 10^{4}$ ultra-short period exoplanets with characteristic strain of the same order or higher than the ones analyzed in our paper. This means that the question related to the gravitational wave pattern of ultra-short period exoplanets may be surpassed near future by the LISA instrument with new and more definitive data.


title: Least Squares Normalized Cross Correlation

abstract: Direct methods are widely used for alignment of models to images, due to their accuracy, since they minimize errors in the domain of measurement noise. They have leveraged least squares minimizations, for simple, efficient, variational optimization, since the seminal 1981 work of Lucas & Kanade, and normalized cross correlation (NCC), for robustness to intensity variations, since at least 1972. Despite the complementary benefits of these two well known methods, they have not been effectively combined to address local variations in intensity. Many ad-hoc NCC frameworks, sub-optimal least squares methods and image transformation approaches have thus been proposed instead, each with their own limitations. This work shows that a least squares optimization of NCC without approximation is not only possible, but straightforward and efficient. A robust, locally normalized formulation is introduced to mitigate local intensity variations and partial occlusions. Finally, sparse features with oriented patches are proposed for further efficiency. The resulting framework is simple to implement, computationally efficient and robust to local intensity variations. It is evaluated on the image alignment problem, showing improvements in both convergence rate and computation time over existing lighting invariant methods.


title: The Andoni--Krauthgamer--Razenshteyn characterization of sketchable norms fails for sketchable metrics

abstract: Andoni, Krauthgamer and Razenshteyn (AKR) proved (STOC 2015) that a finite-dimensional normed space $(X,\|\cdot\|_X)$ admits a $O(1)$ sketching algorithm (namely, with $O(1)$ sketch size and $O(1)$ approximation) if and only if for every $\varepsilon\in (0,1)$ there exist $\alpha\geqslant 1$ and an embedding $f:X\to \ell_{1-\varepsilon}$ such that $\|x-y\|_X\leqslant \|f(x)-f(y)\|_{1-\varepsilon}\leqslant \alpha \|x-y\|_X$ for all $x,y\in X$. The "if part" of this theorem follows from a sketching algorithm of Indyk (FOCS 2000). The contribution of AKR is therefore to demonstrate that the mere availability of a sketching algorithm implies the existence of the aforementioned geometric realization. Indyk's algorithm shows that the "if part" of the AKR characterization holds true for any metric space whatsoever, i.e., the existence of an embedding as above implies sketchability even when $X$ is not a normed space. Due to this, a natural question that AKR posed was whether the assumption that the underlying space is a normed space is needed for their characterization of sketchability. We resolve this question by proving that for arbitrarily large $n\in \mathbb{N}$ there is an $n$-point metric space $(M(n),d_{M(n)})$ which is $O(1)$-sketchable yet for every $\varepsilon\in (0,\frac12)$, if $\alpha(n)\geqslant 1$ and $f_n:M(n)\to \ell_{1-\varepsilon}$ are such that $d_{M(n)}(x,y)\leqslant \|f_n(x)-f_n(y)\|_{1-\varepsilon}\leqslant \alpha(n) d_{M(n)}(x,y)$ for all $x,y\in M(n)$, then necessarily $\lim_{n\to \infty} \alpha(n)= \infty$.


title: Compact and lightweight 1.5 {\mu}m lidar with a multi-mode fiber coupling free-running InGaAs/InP single-photon detector

abstract: We present a compact and lightweight 1.5 {\mu}m lidar using a free-running single-photon detector (SPD) based on a multi-mode fiber (MMF) coupling InGaAs/InP negative feedback avalanche diode. The ultimate light detection sensitivity of SPD highly reduces the power requirement of laser, whilst the enhanced collection efficiency due to MMF coupling significantly reduces the volume and weight of telescopes. We develop a specific algorithm for the corrections of errors caused by the SPD and erbium-doped fiber amplifier to extract accurate backscattering signals. We also perform a comparison between single-mode fiber (SMF) coupling and MMF coupling in the lidar receiver, and the results show that the collection efficiency with MMF coupling is five times higher than SMF coupling. In order to validate the functionality, we use the lidar system for the application of cloud detection. The lidar system exhibits the ability to detect both the cloud base height and the thickness of multi-layer clouds to an altitude of 12 km with a temporal resolution of 1 s and a spatial resolution of 15 m. Due to the advantages of compactness and lightweight, our lidar system can be installed on unmanned aerial vehicles for wide applications in practice.


title: Test particles in relativistic resistive magnetohydrodynamics

abstract: The Black Hole Accretion Code (BHAC) has recently been extended with the ability to evolve charged test particles according to the Lorentz force within resistive relativistic magnetohydrodynamics simulations. We apply this method to evolve particles in a reconnecting current sheet that forms due to the coalescence of two magnetic flux tubes in 2D Minkowski spacetime. This is the first analysis of charged test particle evolution in resistive relativistic magnetohydrodynamics simulations. The energy distributions of an ensemble of 100.000 electrons are analyzed, as well as the acceleration of particles in the plasmoids that form in the reconnection layer. The effect of the Lundquist number, magnetization, and plasma-$\beta$ on the particle energy distribution is explored for a range of astrophysically relevant parameters. We find that electrons accelerate to non-thermal energies in the thin current sheets in all cases. We find two separate acceleration regimes: An exponential increase of the Lorentz factor during the island coalescence where the acceleration depends linearly on the resistivity and a nonlinear phase with high variability. These results are relevant for determining energy distributions and acceleration sites obtaining radiation maps in large-scale magnetohydrodynamics simulations of black hole accretion disks and jets.


title: Particle orbits at the magnetopause: Kelvin-Helmholtz induced trapping

abstract: The Kelvin-Helmholtz instability (KHI) is a known mechanism for penetration of solar wind matter into the magnetosphere. Using three-dimensional, resistive magnetohydrodynamic simulations, the double mid-latitude reconnection (DMLR) process was shown to efficiently exchange solar wind matter into the magnetosphere, through mixing and reconnection. Here, we compute test particle orbits through DMLR configurations. In the instantaneous electromagnetic fields, charged particle trajectories are integrated using the guiding centre approximation. The mechanisms involved in the electron particle orbits and their kinetic energy evolutions are studied in detail, to identify specific signatures of the DMLR through particle characteristics. The charged particle orbits are influenced mainly by magnetic curvature drifts. We identify complex, temporarily trapped, trajectories where the combined electric field and (reconnected) magnetic field variations realize local cavities where particles gain energy before escaping. By comparing the orbits in strongly deformed fields due to the KHI development, with the textbook mirror-drift orbits resulting from our initial configuration, we identify effects due to current sheets formed in the DMLR process. We do this in various representative stages during the DMLR development.


title: Analysis of Maximal Topologies Achieving Optimal DoF and DoF $\frac{1}{n}$ in Topological Interference Management

abstract: Topological interference management (TIM) can obtain degrees of freedom (DoF) gains with no channel state information at the transmitters (CSIT) except topological information of network in the interference channel. It was shown that TIM achieves the optimal symmetric DoF when internal conflict does not exist among messages. However, it is difficult to assure whether a specific topology can achieve the optimal DoF without scrutinizing internal conflict, which requires lots of works. Also, it is hard to design a specific optimal topology directly from the conventional condition for the optimal DoF. With these problems in mind, we propose a method to derive maximal topology directly in TIM, named as alliance construction in K-user interference channel. That is, it is proved that a topology is maximal if and only if it is derived from alliance construction. We translate a topology design by alliance construction in message graph into topology matrix and propose conditions for maximal topology matrix (MTM). Moreover, we propose a generalized alliance construction that derives a topology achieving DoF 1/n for n>=3 by generalizing sub-alliances. A topology matrix can also be used to analyze maximality of topology with DoF 1/n.


title: Prediction-based compensation for gate on/off latency during respiratory-gated radiotherapy

abstract: During respiratory-gated radiotherapy (RGRT), gate on and off latencies cause deviations of gating windows, possibly leading to delivery of low- and high-dose radiations to tumors and normal tissues, respectively. Currently, there are no RGRT systems that have definite tools to compensate for the delays. To address the problem, we propose a framework consisting of two steps: 1) multi-step-ahead prediction and 2) prediction-based gating. For each step, we have devised a specific algorithm to accomplish the task. Numerical experiments were performed using respiratory signals of a phantom and ten volunteers, and our prediction-based RGRT system exhibited superior performance in more than a few signal samples. In some, however, signal prediction and prediction-based gating did not work well, maybe due to signal irregularity and/or baseline drift. The proposed approach has potential applicability in RGRT, and further studies are needed to verify and refine the constituent algorithms.


title: Complementary-Label Learning for Arbitrary Losses and Models

abstract: In contrast to the standard classification paradigm where the true class is given to each training pattern, complementary-label learning only uses training patterns each equipped with a complementary label, which only specifies one of the classes that the pattern does not belong to. The goal of this paper is to derive a novel framework of complementary-label learning with an unbiased estimator of the classification risk, for arbitrary losses and models---all existing methods have failed to achieve this goal. Not only is this beneficial for the learning stage, it also makes model/hyper-parameter selection (through cross-validation) possible without the need of any ordinarily labeled validation data, while using any linear/non-linear models or convex/non-convex loss functions. We further improve the risk estimator by a non-negative correction and gradient ascent trick, and demonstrate its superiority through experiments.


title: Three-loop planar master integrals for heavy-to-light form factors

abstract: We calculate analytically the three-loop planar master integrals relevant for heavy-to-light form factors using the method of differential equations. After choosing a proper canonical basis, the boundary conditions are easy to be determined, and the solution of differential equations is greatly simplified. The results for seventy-one master integrals at general kinematics are all expressed in terms of harmonic polylogarithms.


title: Task Runtime Prediction in Scientific Workflows Using an Online Incremental Learning Approach

abstract: Many algorithms in workflow scheduling and resource provisioning rely on the performance estimation of tasks to produce a scheduling plan. A profiler that is capable of modeling the execution of tasks and predicting their runtime accurately, therefore, becomes an essential part of any Workflow Management System (WMS). With the emergence of multi-tenant Workflow as a Service (WaaS) platforms that use clouds for deploying scientific workflows, task runtime prediction becomes more challenging because it requires the processing of a significant amount of data in a near real-time scenario while dealing with the performance variability of cloud resources. Hence, relying on methods such as profiling tasks' execution data using basic statistical description (e.g., mean, standard deviation) or batch offline regression techniques to estimate the runtime may not be suitable for such environments. In this paper, we propose an online incremental learning approach to predict the runtime of tasks in scientific workflows in clouds. To improve the performance of the predictions, we harness fine-grained resources monitoring data in the form of time-series records of CPU utilization, memory usage, and I/O activities that are reflecting the unique characteristics of a task's execution. We compare our solution to a state-of-the-art approach that exploits the resources monitoring data based on regression machine learning technique. From our experiments, the proposed strategy improves the performance, in terms of the error, up to 29.89%, compared to the state-of-the-art solutions.


title: Phase-mediated magnon chaos-order transition in cavity optomagnonics

abstract: Magnon as a quantized spin wave has attracted extensive attentions in various fields of physics, such as magnon spintronics, microwave photonics, and cavity quantum electrodynamics. Here, we explore theoretically magnon chaos-order transition in cavity optomagnonics, which is still remains largely unexplored in this emerging field. We find that the evolution of the magnon experiences the transition from order to period-doubling bifurcation and finally enter chaos by adjusting the microwave driving power. Different from the normal chaos, the magnon chaos-order transition proposed here is phase-mediated. Beyond their fundamental scientific significance, our results will contribute to the comprehending of nonlinear phenomena and chaos in optomagnonical system, and may find applications in chaos-based secure communication.


title: Assignment Mechanisms under Distributional Constraints

abstract: We study the assignment problem of objects to agents with heterogeneous preferences under distributional constraints. Each agent is associated with a publicly known type and has a private ordinal ranking over objects. We are interested in assigning as many agents as possible. Our first contribution is a generalization of the well-known and widely used serial dictatorship. Our mechanism maintains several desirable properties of serial dictatorship, including strategyproofness, Pareto efficiency, and computational tractability while satisfying the distributional constraints with a small error. We also propose a generalization of the probabilistic serial algorithm, which finds an ordinally efficient and envy-free assignment, and also satisfies the distributional constraints with a small error. We show, however, that no ordinally efficient and envy-free mechanism is also weakly strategyproof. Both of our algorithms assign at least the same number of students as the optimum fractional assignment.


title: L\'{e}vy-walk-like Langevin dynamics

abstract: Continuous time random walks and Langevin equations are two classes of stochastic models for describing the dynamics of particles in the natural world. While some of the processes can be conveniently characterized by both of them, more often one model has significant advantages (or has to be used) compared with the other one. In this paper, we consider the weakly damped Langevin system coupled with a new subordinator|$\alpha$-dependent subordinator with $1<\alpha<2$. We pay attention to the diffusion behaviour of the stochastic process described by this coupled Langevin system, and find the super-ballistic diffusion phenomena for the system with an unconfined potential on velocity but sub-ballistic superdiffusion phenomenon with a confined potential, which is like L\'{e}vy walk for long times. One can further note that the two-point distribution of inverse subordinator affects mean square displacement of this coupled weakly damped Langevin system in essential.


title: Axial Seamount: Periodic tidal loading reveals stress dependence of the earthquake size distribution (b value)

abstract: Earthquake size-frequency distributions commonly follow a power law, with the b value often used to quantify the relative proportion of small and large events. Laboratory experiments have found that the b value of microfractures decreases with increasing stress. Studies have inferred that this relationship also holds for earthquakes based on observations of earthquake b values varying systematically with faulting style, depth, and for subduction zone earthquakes, plate age. However, these studies are limited by small sample sizes despite aggregating events over large regions, which precludes the ability to control for other variables that might also affect earthquake b values such as rock heterogeneity and fault roughness. Our natural experiment in a unique seafloor laboratory on Axial Seamount involves analyzing the size-frequency distribution of ~60,000 microearthquakes which delineate a ring-fault system in a 25 km3 block of crust that experiences periodic tidal loading of +/-20 kPa. We find that above a threshold stress amplitude, b value is inversely correlated with tidal stress. The earthquake b value varies by ~0.09 per kPa change in Coulomb stress. Our results support the potential use of b values to estimate small stress variations in the Earth's crust.


title: GraphMP: I/O-Efficient Big Graph Analytics on a Single Commodity Machine

abstract: Recent studies showed that single-machine graph processing systems can be as highly competitive as cluster-based approaches on large-scale problems. While several out-of-core graph processing systems and computation models have been proposed, the high disk I/O overhead could significantly reduce performance in many practical cases. In this paper, we propose GraphMP to tackle big graph analytics on a single machine. GraphMP achieves low disk I/O overhead with three techniques. First, we design a vertex-centric sliding window (VSW) computation model to avoid reading and writing vertices on disk. Second, we propose a selective scheduling method to skip loading and processing unnecessary edge shards on disk. Third, we use a compressed edge cache mechanism to fully utilize the available memory of a machine to reduce the amount of disk accesses for edges. Extensive evaluations have shown that GraphMP could outperform existing single-machine out-of-core systems such as GraphChi, X-Stream and GridGraph by up to 51, and can be as highly competitive as distributed graph engines like Pregel+, PowerGraph and Chaos.