title: Bayesian Pool-based Active Learning With Abstention Feedbacks

abstract: We study pool-based active learning with abstention feedbacks, where a labeler can abstain from labeling a queried example with some unknown abstention rate. This is an important problem with many useful applications. We take a Bayesian approach to the problem and develop two new greedy algorithms that learn both the classification problem and the unknown abstention rate at the same time. These are achieved by simply incorporating the estimated abstention rate into the greedy criteria. We prove that both of our algorithms have near-optimality guarantees: they respectively achieve a ${(1-\frac{1}{e})}$ constant factor approximation of the optimal expected or worst-case value of a useful utility function. Our experiments show the algorithms perform well in various practical scenarios.


title: New separated polynomial solutions to the Zernike system on the unit disk and interbasis expansion

abstract: The differential equation proposed by Frits Zernike to obtain a basis of polynomial orthogonal solutions on the the unit disk to classify wavefront aberrations in circular pupils, is shown to have a set of new orthonormal solution bases, involving Legendre and Gegenbauer polynomials, in non-orthogonal coordinates close to Cartesian ones. We find the overlaps between the original Zernike basis and a representative of the new set, which turn out to be Clebsch-Gordan coefficients.


title: An explicit symmetric DGLA model of a bi-gon

abstract: We give explicit formulae for a DGLA model of the bi-gon which is symmetric under the geometric symmetries of the cell. This follows the work of Lawrence-Sullivan on the (unique) DGLA model of the interval and its construction uses deeper knowledge of the structure of such models and their localisations for non-simply connected spaces.


title: Dark Energy from pNGB Mediated Dirac Neutrino Condensate

abstract: We consider an extension of the Standard Model that provide an unified description of eV scale neutrino mass and dark energy. An explicit model is presented by augmenting the Standard Model with an $SU(2)_L$ doublet scalar, a singlet scalar and right handed neutrinos where all of them are assumed to be charged under a global $U(1)_X$ symmetry. A light pseudo-Nambu-Goldstone Boson, associated with the spontaneously broken $U(1)_{X}$ symmetry, acts as a mediator of an attractive force leading to a Dirac neutrino condensate, with large correlation length, and a non-zero gap in the right range providing a cosmologically feasible dark energy scenario. The neutrino mass is generated through the usual Dirac seesaw mechanism. Parameter space, reproducing viable dark energy scenario while having neutrino mass in the right ballpark, is presented.


title: Search for the Footprints of New Physics with Laboratory and Cosmic Neutrinos

abstract: Observations of high energy neutrinos, both in the laboratory and from cosmic sources, can be a useful probe in searching for new physics. Such observations can provide sensitive tests of Lorentz invariance violation (LIV), which may be a the result of quantum gravity physics (QG). We review some observationally testable consequences of LIV using effective field theory (EFT) formalism. To do this, one can postulate the existence of additional small LIV terms in free particle Lagrangians, suppressed by powers of the Planck mass. The observational consequences of such terms are then examined. In particular, one can place limits on a class of non-renormalizable, mass dimension five and six Lorentz invariance violating operators that may be the result of QG.


title: FCNC Decays of the Top Quark

abstract: If new physics (e.g. SUSY) does not show up as direct evidence at the LHC, it could still be observable in FCNC processes involving the $t$-quark. We take a close look at the process $t\to c + h/Z$ and show that its branching ratio in the Standard Model is subject to three mechanisms of suppression. To obtain an observable signal, one needs to evade all these mechanisms in a theory beyond the Standard Model. We show that a theory like the cMSSM cannot provide a big enough enhancement. However, in a framework like $R$-parity-violating SUSY, observable signals are a distinct possibility.


title: Transport by Lagrangian Vortices in the Eastern Pacific

abstract: Rotationally coherent Lagrangian vortices (RCLVs) are identified from satellite-derived surface geostrophic velocities in the Eastern Pacific (180$^\circ$-130$^\circ$ W) using the objective (frame-invariant) finite-time Lagrangian-coherent-structure detection method of Haller et al. (2016) based on the Lagrangian-averaged vorticity deviation. RCLVs are identified for 30, 90, and 270 day intervals over the entire satellite dataset, beginning in 1993. In contrast to structures identified using Eulerian eddy-tracking methods, the RCLVs maintain material coherence over the specified time intervals, making them suitable for material transport estimates. Statistics of RCLVs are compared to statistics of eddies identified from sea-surface height (SSH) by Chelton et al. 2011. RCLVs and SSH eddies are found to propagate westward at similar speeds at each latitude, consistent with the Rossby wave dispersion relation. However, RCLVs are uniformly smaller and shorter-lived than SSH eddies. A coherent eddy diffusivity is derived to quantify the contribution of RCLVs to meridional transport; it is found that RCLVs contribute less than 1% to net meridional dispersion and diffusion in this sector, implying that eddy transport of tracers is mostly due to incoherent motions, such as swirling and filamentation outside of the eddy cores, rather than coherent meridional translation of eddies themselves. These findings call into question prior estimates of coherent eddy transport based on Eulerian eddy identification methods.


title: Second-Order Word Embeddings from Nearest Neighbor Topological Features

abstract: We introduce second-order vector representations of words, induced from nearest neighborhood topological features in pre-trained contextual word embeddings. We then analyze the effects of using second-order embeddings as input features in two deep natural language processing models, for named entity recognition and recognizing textual entailment, as well as a linear model for paraphrase recognition. Surprisingly, we find that nearest neighbor information alone is sufficient to capture most of the performance benefits derived from using pre-trained word embeddings. Furthermore, second-order embeddings are able to handle highly heterogeneous data better than first-order representations, though at the cost of some specificity. Additionally, augmenting contextual embeddings with second-order information further improves model performance in some cases. Due to variance in the random initializations of word embeddings, utilizing nearest neighbor features from multiple first-order embedding samples can also contribute to downstream performance gains. Finally, we identify intriguing characteristics of second-order embedding spaces for further research, including much higher density and different semantic interpretations of cosine similarity.


title: Securing Real-Time Internet-of-Things

abstract: Modern embedded and cyber-physical systems are ubiquitous. A large number of critical cyber-physical systems have real-time requirements (e.g., avionics, automobiles, power grids, manufacturing systems, industrial control systems, etc.). Recent developments and new functionality requires real-time embedded devices to be connected to the Internet. This gives rise to the real-time Internet-of-things (RT-IoT) that promises a better user experience through stronger connectivity and efficient use of next-generation embedded devices. However RT- IoT are also increasingly becoming targets for cyber-attacks which is exacerbated by this increased connectivity. This paper gives an introduction to RT-IoT systems, an outlook of current approaches and possible research challenges towards secure RT- IoT frameworks.


title: Distortion and the bridge distance of knots

abstract: We extend techniques due to Pardon to show that there is a lower bound on the distortion of a knot in $\mathbb{R}^3$ proportional to the minimum of the bridge distance and the bridge number of the knot. We also exhibit an infinite family of knots for which the minimum of the bridge distance and the bridge number is unbounded and Pardon's lower bound is constant.


title: Stochastic Replica Voting Machine Prediction of Stable Cubic and Double Perovskite Materials and Binary Alloys

abstract: A machine learning approach that we term the `Stochastic Replica Voting Machine' (SRVM) algorithm is presented and applied to a binary and a 3-class classification problems in materials science. Here, we employ SRVM to predict candidate compounds capable of forming stable perovskites and double perovskites and further classify binary ($AB$) solids. The results of our binary and ternary classifications compared well to those obtained by SVM and neural network algorithms.


title: Uplift Modeling with Multiple Treatments and General Response Types

abstract: Randomized experiments have been used to assist decision-making in many areas. They help people select the optimal treatment for the test population with certain statistical guarantee. However, subjects can show significant heterogeneity in response to treatments. The problem of customizing treatment assignment based on subject characteristics is known as uplift modeling, differential response analysis, or personalized treatment learning in literature. A key feature for uplift modeling is that the data is unlabeled. It is impossible to know whether the chosen treatment is optimal for an individual subject because response under alternative treatments is unobserved. This presents a challenge to both the training and the evaluation of uplift models. In this paper we describe how to obtain an unbiased estimate of the key performance metric of an uplift model, the expected response. We present a new uplift algorithm which creates a forest of randomized trees. The trees are built with a splitting criterion designed to directly optimize their uplift performance based on the proposed evaluation method. Both the evaluation method and the algorithm apply to arbitrary number of treatments and general response types. Experimental results on synthetic data and industry-provided data show that our algorithm leads to significant performance improvement over other applicable methods.


title: Asymmetric product of left braces and simplicity; new solutions of the Yang-Baxter equation

abstract: The problem of constructing all the non-degenerate involutive set-theoretic solutions of the Yang-Baxter equation recently has been reduced to the problem of describing all the left braces. In particular, the classification of all finite left braces is fundamental in order to describe all finite such solutions of the Yang-Baxter equation. In this paper we continue the study of finite simple left braces with the emphasis on the application of the asymmetric product of left braces in order to construct new classes of simple left braces. We do not only construct new classes but also we interpret all previously known constructions as asymmetric products. Moreover, a construction is given of finite simple left braces with a multiplicative group that is solvable of arbitrary derived length.


title: Asynchronous Coordinate Descent under More Realistic Assumptions

abstract: Asynchronous-parallel algorithms have the potential to vastly speed up algorithms by eliminating costly synchronization. However, our understanding to these algorithms is limited because the current convergence of asynchronous (block) coordinate descent algorithms are based on somewhat unrealistic assumptions. In particular, the age of the shared optimization variables being used to update a block is assumed to be independent of the block being updated. Also, it is assumed that the updates are applied to randomly chosen blocks. In this paper, we argue that these assumptions either fail to hold or will imply less efficient implementations. We then prove the convergence of asynchronous-parallel block coordinate descent under more realistic assumptions, in particular, always without the independence assumption. The analysis permits both the deterministic (essentially) cyclic and random rules for block choices. Because a bound on the asynchronous delays may or may not be available, we establish convergence for both bounded delays and unbounded delays. The analysis also covers nonconvex, weakly convex, and strongly convex functions. We construct Lyapunov functions that directly model both objective progress and delays, so delays are not treated errors or noise. A continuous-time ODE is provided to explain the construction at a high level.


title: A study on exponential-size neighborhoods for the bin packing problem with conflicts

abstract: We propose an iterated local search based on several classes of local and large neighborhoods for the bin packing problem with conflicts. This problem, which combines the characteristics of both bin packing and vertex coloring, arises in various application contexts such as logistics and transportation, timetabling, and resource allocation for cloud computing. We introduce $O(1)$ evaluation procedures for classical local-search moves, polynomial variants of ejection chains and assignment neighborhoods, an adaptive set covering-based neighborhood, and finally a controlled use of 0-cost moves to further diversify the search. The overall method produces solutions of good quality on the classical benchmark instances and scales very well with an increase of problem size. Extensive computational experiments are conducted to measure the respective contribution of each proposed neighborhood. In particular, the 0-cost moves and the large neighborhood based on set covering contribute very significantly to the search. Several research perspectives are open in relation to possible hybridizations with other state-of-the-art mathematical programming heuristics for this problem.


title: Stochastic Ratcheting on a Funneled Energy Landscape is Necessary for Highly Efficient Contractility of Actomyosin Force Dipoles

abstract: Current understanding of how contractility emerges in disordered actomyosin networks of non-muscle cells is still largely based on the intuition derived from earlier works on muscle contractility. This view, however, largely overlooks the free energy gain following passive cross-linker binding, which, even in the absence of active fluctuations, provides a thermodynamic drive towards highly overlapping filamentous states. In this work, we shed light on this phenomenon, showing that passive cross-linkers, when considered in the context of two anti-parallel filaments, generate noticeable contractile forces. However, as binding free energy of cross-linkers is increased, a sharp onset of kinetic arrest follows, greatly diminishing effectiveness of this contractility mechanism, allowing the network to contract only with weakly resisting tensions at its boundary. We have carried out stochastic simulations elucidating this mechanism, followed by a mean-field treatment that predicts how contractile forces asymptotically scale at small and large binding energies, respectively. Furthermore, when considering an active contractile filament pair, based on non-muscle myosin II, we found that the non-processive nature of these motors leads to highly inefficient force generation, due to recoil slippage of the overlap during periods when the motor is dissociated. However, we discovered that passive cross-linkers can serve as a structural ratchet during these unbound motor time spans, resulting in vast force amplification. Our results shed light on the non-equilibrium effects of transiently binding proteins in biological active matter, as observed in the non-muscle actin cytoskeleton, showing that highly efficient contractile force dipoles result from synergy of passive cross-linker and active motor dynamics, via a ratcheting mechanism on a funneled energy landscape.


title: Discovery of Galactic OIV and OV X-ray absorption due to transition temperature gas in the PKS 2155-304 spectrum

abstract: FUV observations have revealed the transition temperature gas (TTG; $\log{T({\mathrm{K}})}$ ~ 5), located in the lower Galactic halo and in the High Velocity Clouds. However, the corresponding X-ray absorption has so far remained mostly undetected. In order to make an improvement in this respect in Galactic X-ray absorption studies, we accumulated very deep (~3 Ms) spectra of the blazar PKS 2155-304 obtained with the spectrometers RGS1, RGS2, LETG/HRC and LETG/ACIS-S and studied the absorption lines due to the intervening Galactic components. The very high quality of the data and the coverage of important wavelengths with at least two independent instruments allowed us to reliably detect ten Galactic lines with better than 99.95% confidence. We discovered significant absorption from blended OIV transitions 1s-2p $^2$S (22.571 \AA ), 1s-2p $^2$P (22.741 \AA ) and 1s-2p $^2$D (22.777 \AA ), and from the OV transition 1s-2p (22.370 \AA ) from TTG at $\log{T({\mathrm{K}})} \thinspace = \thinspace 5.2\pm0.1$. A joint X-ray and FUV analysis indicated that photoionisation is negligible for this component and that the gas is in a cooling transition phase. However, the temperature is high enough that the column density ratio N(OIV)/N(OV) is not significantly different from that in collisional ionisation equilibrium (CIE). Under CIE we obtained $N_{\mathrm{OIV}}$= 3.6$\pm$2.0 $\times 10^{15}$ cm$^{-2}$, corresponding to $N_{\mathrm{H}}$ = 1.0$\pm$0.5 $\times 10^{19} \frac{Z_{\odot}}{Z_{\mathrm{TTG}}}$ cm$^{-2}$.


title: Clinical Intervention Prediction and Understanding using Deep Networks

abstract: Real-time prediction of clinical interventions remains a challenge within intensive care units (ICUs). This task is complicated by data sources that are noisy, sparse, heterogeneous and outcomes that are imbalanced. In this paper, we integrate data from all available ICU sources (vitals, labs, notes, demographics) and focus on learning rich representations of this data to predict onset and weaning of multiple invasive interventions. In particular, we compare both long short-term memory networks (LSTM) and convolutional neural networks (CNN) for prediction of five intervention tasks: invasive ventilation, non-invasive ventilation, vasopressors, colloid boluses, and crystalloid boluses. Our predictions are done in a forward-facing manner to enable "real-time" performance, and predictions are made with a six hour gap time to support clinically actionable planning. We achieve state-of-the-art results on our predictive tasks using deep architectures. We explore the use of feature occlusion to interpret LSTM models, and compare this to the interpretability gained from examining inputs that maximally activate CNN outputs. We show that our models are able to significantly outperform baselines in intervention prediction, and provide insight into model learning, which is crucial for the adoption of such models in practice.


title: The Prediction Advantage: A Universally Meaningful Performance Measure for Classification and Regression

abstract: We introduce the Prediction Advantage (PA), a novel performance measure for prediction functions under any loss function (e.g., classification or regression). The PA is defined as the performance advantage relative to the Bayesian risk restricted to knowing only the distribution of the labels. We derive the PA for well-known loss functions, including 0/1 loss, cross-entropy loss, absolute loss, and squared loss. In the latter case, the PA is identical to the well-known R-squared measure, widely used in statistics. The use of the PA ensures meaningful quantification of prediction performance, which is not guaranteed, for example, when dealing with noisy imbalanced classification problems. We argue that among several known alternative performance measures, PA is the best (and only) quantity ensuring meaningfulness for all noise and imbalance levels.


title: Selective Classification for Deep Neural Networks

abstract: Selective classification techniques (also known as reject option) have not yet been considered in the context of deep neural networks (DNNs). These techniques can potentially significantly improve DNNs prediction performance by trading-off coverage. In this paper we propose a method to construct a selective classifier given a trained neural network. Our method allows a user to set a desired risk level. At test time, the classifier rejects instances as needed, to grant the desired risk (with high probability). Empirical results over CIFAR and ImageNet convincingly demonstrate the viability of our method, which opens up possibilities to operate DNNs in mission-critical applications. For example, using our method an unprecedented 2% error in top-5 ImageNet classification can be guaranteed with probability 99.9%, and almost 60% test coverage.


title: Characteristics of On-time and Late Reward Delivery Projects

abstract: The crowdfunding market size has increased exponentially, reaching tens of billions of dollars and showing the popularity of crowdfunding. However, according to Kickstarter, 35% backers did not receive rewards on time. To maintain the trust between creators and backers, and sustain the crowdfunding business growth, it is crucial to understand how on-time and late reward delivery projects are different. In this paper, we analyze characteristics of on-time and late reward delivery projects, especially, focusing on project descriptions, creator profiles, and activeness and linguistic patterns of creators and backers. Our analysis reveals that the less complicated a project is and more actively a creator responds to backers, the higher an on-time reward delivery probability is. It shows there are significant differences between on-time and late reward delivery projects.


title: Culture, Computation, Morality

abstract: I point to a deep and unjustly ignored relation between culture and computation. I first establish interpretations of Piaget's and Vygotsky's theories of child development with the language of theoretical computer science. Using these interpretations, I argue that the two different possible routes to Piagetian disequilibrium -- a tendency to overaccommodate, and a tendency to overassimilate -- are equivalent to the two distinct cultural tendencies, collectivistism and individualism. I argue that this simple characterization of overaccommodation versus overassimilation provides a satisfying explanation as to why the two cultural tendencies differ in the way they empirically do. All such notions are grounded on a firm mathematical framework for those who prefer the computable, and grounded on my personal history for those who prefer the uncomputable.


title: The Geometry and Topology of Data and Information for Analytics of Processes and Behaviours: Building on Bourdieu and Addressing New Societal Challenges

abstract: We begin by summarizing the relevance and importance of inductive analytics based on the geometry and topology of data and information. Contemporary issues are then discussed. These include how sampling data for representativity is increasingly to be questioned. While we can always avail of analytics from a "bag of tools and techniques", in the application of machine learning and predictive analytics, nonetheless we present the case for Bourdieu and Benz\'ecri-based science of data, as follows. This is to construct bridges between data sources and position-taking, and decision-making. There is summary presentation of a few case studies, illustrating and exemplifying application domains.


title: Interpreting Blackbox Models via Model Extraction

abstract: Interpretability has become incredibly important as machine learning is increasingly used to inform consequential decisions. We propose to construct global explanations of complex, blackbox models in the form of a decision tree approximating the original model---as long as the decision tree is a good approximation, then it mirrors the computation performed by the blackbox model. We devise a novel algorithm for extracting decision tree explanations that actively samples new training points to avoid overfitting. We evaluate our algorithm on a random forest to predict diabetes risk and a learned controller for cart-pole. Compared to several baselines, our decision trees are both substantially more accurate and equally or more interpretable based on a user study. Finally, we describe several insights provided by our interpretations, including a causal issue validated by a physician.


title: Universal terms in the overlap of the ground state of the spin-1/2 XXZ chain with the N\'eel state

abstract: We analyze universal terms that appear in the large system size scaling of the overlap between the N\'eel state and the ground state of the spin-1/2 XXZ chain in the antiferromagnetic regime. In a critical theory, the order one term of the asymptotics of such an overlap may be expressed in terms of $g$-factors, known in the context of boundary conformal field theory. In particular, for the XXZ model in its gapless phase, this term provides access to the Luttinger parameter. In its gapped phase, on the other hand, the order one term simply reflects the symmetry broken nature of the phase. In order to study the large system size scaling of this overlap analytically and to compute the order one term exactly, we use a recently derived finite-size determinant formula and perform an asymptotic expansion. Our analysis confirms the predictions of boundary conformal field theory and enables us to determine the exponent of the leading finite-size correction.


title: HuGaDB: Human Gait Database for Activity Recognition from Wearable Inertial Sensor Networks

abstract: This paper presents a human gait data collection for analysis and activity recognition consisting of continues recordings of combined activities, such as walking, running, taking stairs up and down, sitting down, and so on; and the data recorded are segmented and annotated. Data were collected from a body sensor network consisting of six wearable inertial sensors (accelerometer and gyroscope) located on the right and left thighs, shins, and feet. Additionally, two electromyography sensors were used on the quadriceps (front thigh) to measure muscle activity. This database can be used not only for activity recognition but also for studying how activities are performed and how the parts of the legs move relative to each other. Therefore, the data can be used (a) to perform health-care-related studies, such as in walking rehabilitation or Parkinson's disease recognition, (b) in virtual reality and gaming for simulating humanoid motion, or (c) for humanoid robotics to model humanoid walking. This dataset is the first of its kind which provides data about human gait in great detail. The database is available free of charge https://github.com/romanchereshnev/HuGaDB.


title: Teaching computer code at school

abstract: In today's education systems, there is a deep concern about the importance of teaching code and computer programming in schools. Moving digital learning from a simple use of tools to understanding the processes of the internal functioning of these tools is an old / new debate originated with the digital laboratories of the 1960. Today, it is emerging again under impulse of the large - scale public sphere digitalization and the new constructivist education theories. Teachers and educators discuss not only the viability of code teaching in the classroom, but also the intellectual and cognitive advantages for students. The debate thus takes several orientations and is resourced in the entanglement of arguments and interpretations of any order, technical, educational, cultural, cognitive and psychological. However, that phenomenon which undoubtedly augurs for a profound transformation in the future models of learning and teaching , is predicting a new and almost congenital digital humanism


title: Vehicle Traffic Driven Camera Placement for Better Metropolis Security Surveillance

abstract: Security surveillance is one of the most important issues in smart cities, especially in an era of terrorism. Deploying a number of (video) cameras is a common surveillance approach. Given the never-ending power offered by vehicles to metropolises, exploiting vehicle traffic to design camera placement strategies could potentially facilitate security surveillance. This article constitutes the first effort toward building the linkage between vehicle traffic and security surveillance, which is a critical problem for smart cities. We expect our study could influence the decision making of surveillance camera placement, and foster more research of principled ways of security surveillance beneficial to our physical-world life. Code has been made publicly available.


title: Predictive Analytics for Enhancing Travel Time Estimation in Navigation Apps of Apple, Google, and Microsoft

abstract: The explosive growth of the location-enabled devices coupled with the increasing use of Internet services has led to an increasing awareness of the importance and usage of geospatial information in many applications. The navigation apps (often called Maps), use a variety of available data sources to calculate and predict the travel time as well as several options for routing in public transportation, car or pedestrian modes. This paper evaluates the pedestrian mode of Maps apps in three major smartphone operating systems (Android, iOS and Windows Phone). In the paper, we will show that the Maps apps on iOS, Android and Windows Phone in pedestrian mode, predict travel time without learning from the individual's movement profile. In addition, we will exemplify that those apps suffer from a specific data quality issue which relates to the absence of information about location and type of pedestrian crossings. Finally, we will illustrate learning from movement profile of individuals using various predictive analytics models to improve the accuracy of travel time estimation.


title: Discontinuous Hamiltonian Monte Carlo for discrete parameters and discontinuous likelihoods

abstract: Hamiltonian Monte Carlo has emerged as a standard tool for posterior computation. In this article, we present an extension that can efficiently explore target distributions with discontinuous densities. Our extension in particular enables efficient sampling from ordinal parameters though embedding of probability mass functions into continuous spaces. We motivate our approach through a theory of discontinuous Hamiltonian dynamics and develop a corresponding numerical solver. The proposed solver is the first of its kind, with a remarkable ability to exactly preserve the Hamiltonian. We apply our algorithm to challenging posterior inference problems to demonstrate its wide applicability and competitive performance.


title: Lozi-like maps

abstract: We define a broad class of piecewise smooth plane homeomorphisms which have properties similar to the properties of Lozi maps, including the existence of a hyperbolic attractor. We call those maps Lozi-like. For those maps one can apply our previous results on kneading theory for Lozi maps. We show a strong numerical evidence that there exist Lozi-like maps that have kneading sequences different than those of Lozi maps.


title: Low ionospheric reactions on tropical depressions prior hurricanes

abstract: We study the reactions of the low ionosphere during tropical depressions (TDs) which have been detected before the hurricane appearances in the Atlantic Ocean. We explore 41 TD events using very low frequency (VLF) radio signals emitted by NAA transmitter located in the USA and recorded by VLF receiver located in Belgrade (Serbia). We found VLF signal deviations (caused ionospheric turbulence) in the case of 36 out of 41 TD events (88%). Additionally, we explore 27 TDs which have not been developed in hurricanes and found similar low ionospheric reactions. However, in the sample of 41 TDs which are followed by hurricanes the typical low ionosphere perturbations seem to be more frequent than other TDs.


title: Extensive Globular Cluster Systems Associated with Ultra Diffuse Galaxies in the Coma Cluster

abstract: We present Hubble Space Telescope imaging of two ultra diffuse galaxies (UDGs) with measured stellar velocity dispersions in the Coma cluster. The galaxies, Dragonfly 44 and DFX1, have effective radii of 4.7 kpc and 3.5 kpc and velocity dispersions of $47^{+8}_{-6}$ km/s and $30^{+7}_{-7}$ km/s, respectively. Both galaxies are associated with a striking number of compact objects, tentatively identified as globular clusters: $N_{\rm gc}=74\pm 18$ for Dragonfly 44 and $N_{\rm gc}=62\pm 17$ for DFX1. The number of globular clusters is far higher than expected from the luminosities of the galaxies but is consistent with expectations from the empirical relation between dynamical mass and globular cluster count defined by other galaxies. Combining our data for these two objects with previous HST observations of Coma UDGs we find that UDGs have a factor of $6.9^{+1.0}_{-2.4}$ more globular clusters than other galaxies of the same luminosity, in contrast to a recent study of a similar sample by Amorisco et al. (2017), but consistent with earlier results for individual galaxies. The Harris et al. (2017) relation between globular cluster count and dark matter halo mass implies a median halo mass of $M_{\rm halo}\sim 1.5\times 10^{11}\,{\rm M}_{\odot}$ for the sixteen Coma UDGs that have been observed with HST so far, with the largest and brightest having $M_{\rm halo}\sim 5\times 10^{11}\,{\rm M}_{\odot}$.


title: Cybernetic Health

abstract: Future health ecosystems demand the integration of emerging data technology with an increased focus on preventive medicine. Cybernetics extracts the full potential of data to serve the spectrum of health care, from acute to chronic problems. Building actionable cybernetic navigation tools can greatly empower optimal health decisions, especially by quantifying lifestyle and environmental data. This data to decisions transformation is powered by intuitive event analysis to offer the best semantic abstraction of dynamic living systems. Achieving the goal of preventive health systems in the cybernetic model occurs through the flow of several components. From personalized models we can predict health status using perpetual sensing and data streams. Given these predictions, we give precise recommendations to best suit the prediction for that individual. To enact these recommendations we use persuasive technology in order to deliver and execute targeted interventions.


title: Half-magnetization plateau in a Heisenberg antiferromagnet on a triangular lattice

abstract: We present the phase diagram in a magnetic field of a 2D isotropic Heisenberg antiferromagnet on a triangular lattice. We consider spin-$S$ model with nearest-neighbor ($J_1$) and next-nearest-neighbor ($J_2$) interactions. We focus on the range of $1/8<J_2/J_1<1$, where the ordered states are different from those in the model with only nearest neighbor exchange. A classical ground state in this range has four sublattices and is infinitely degenerate in any field. The actual order is then determined by quantum fluctuations via "order from disorder" phenomenon. We argue that the phase diagram is rich due to competition between competing four-sublattice quantum states which break either $\mathbb{Z}_3$ orientational symmetry or $\mathbb{Z}_4$ sublattice symmetry. At small and high fields, the ground state is a $\mathbb{Z}_3$-breaking canted stripe state, but at intermediate fields the ordered states break $\mathbb{Z}_4$ sublattice symmetry. The most noticeable of such states is "three up, one down" state in which spins in three sublattices are directed along the field and in one sublattice opposite to the field. Such a state breaks no continuous symmetry and has gapped excitations. As the consequence, magnetization has a plateau at exactly one half of the saturation value. We identify gapless states, which border the "three up, one down" state and discuss the transitions between these states and the canted stripe state.


title: An Open-Data Analysis of Heterogeneities in Lung Cancer Premature Mortality Rate and Associated Factors among Toronto Neighborhoods

abstract: In public health, various data are rigorously collected and published with open access. These data reflect the environmental and non-environmental characteristics of heterogeneous neighborhoods in cities. In the present study, we aimed to study the relations between these data and disease risks in heterogeneous neighborhoods. A flexible framework was developed to determine the key factors correlated with diseases and find the most relevant combination of factors to explain observations of diseases through nonlinear analyses. Taking Lung Cancer Premature Mortality Rate (LCPMR) in Toronto as an example, two environmental factors (green space, and industrial pollution) and two non-environmental factors (immigrants, and mental health visits) were identified in the relational analysis of all of the target neighborhoods. To determine the influence of the heterogeneity of the neighborhoods, they were clustered into three different classes. In the most severe class, two additional factors related to dwellings were determined to be involved, which increased the observation's deviance from 48.1% to 80%. The factors determined in this study may assist governments in improving public health policies.


title: Numerical analysis of flexible ultra-thin conformal antenna array

abstract: We present a flexible ultra-thin conformal patch antenna array operating in the Industrial, Scientific, and Medical (ISM) band (2.4-2.4835 GHz) for wireless sensor nodes. The antenna was designed on a flexible polyimide film (Kapton HN) substrate with the size of 95180.1 mm3 and the number of 7 radiating elements. The antenna exhibits S11 < -10 dB at the band of interest with broadband behavior. The measured gain is 4.72 dB in vertical polarization at 2.48 GHz. To obtain the maximum gain of the array, distances of the elements in the array are optimized by using the genetic algorithm (GA) in conjunction with the multilevel fast multipole algorithm (MLFMA). The VSWR is less than 2 over 2.408-2.555 GHz, frequency band. The method of moment (MOM) based on the electric field integral equation (EFIE) is applied to analyze the radiation characteristics of the proposed antenna. The resulting matrix equation obtained by the discretization of the EFIE is solved iteratively using the MLFMA. The radiation patterns and return loss obtained from numerical analysis and simulations show a good impedance matching and a gain enhancement of the proposed antenna.


title: Sideband cooling of small ion Coulomb crystals in a Penning trap

abstract: We have recently demonstrated the laser cooling of a single $^{40}$Ca$^+$ ion to the motional ground state in a Penning trap using the resolved-sideband cooling technique on the electric quadrupole transition S$_{1/2} \leftrightarrow$ D$_{5/2}$. Here we report on the extension of this technique to small ion Coulomb crystals made of two or three $^{40}$Ca$^+$ ions. Efficient cooling of the axial motion is achieved outside the Lamb-Dicke regime on a two-ion string along the magnetic field axis as well as on two- and three-ion planar crystals. Complex sideband cooling sequences are required in order to cool both axial degrees of freedom simultaneously. We measure a mean excitation after cooling of $\bar n_\text{COM}=0.30(4)$ for the centre of mass mode and $\bar n_\text{B}=0.07(3)$ for the breathing mode of the two-ion string with corresponding heating rates of 11(2) s$^{-1}$ and 1(1) s$^{-1}$ at a trap frequency of 162 kHz. The ground state occupation of the axial modes is above 75% for the two-ion planar crystal and the associated heating rates 0.8(5) s$^{-1}$ at a trap frequency of 355 kHz.


title: Ergodicity of Bowen-Margulis measure for the Benoist 3-manifolds

abstract: We study the geodesic flow of a class of 3-manifolds introduced by Benoist which have some hyperbolicity but are non-Riemannian, not CAT(0), and with non-C^1 geodesic flow. The geometries are nonstrictly convex Hilbert geometries in dimension three which admit compact quotient manifolds by discrete groups of projective transformations. We prove the Patterson-Sullivan density is canonical, with applications to counting, and construct explicitly the Bowen-Margulis measure of maximal entropy. The main result of this work is ergodicity of the Bowen-Margulis measure.


title: An effective algorithm for hyperparameter optimization of neural networks

abstract: A major challenge in designing neural network (NN) systems is to determine the best structure and parameters for the network given the data for the machine learning problem at hand. Examples of parameters are the number of layers and nodes, the learning rates, and the dropout rates. Typically, these parameters are chosen based on heuristic rules and manually fine-tuned, which may be very time-consuming, because evaluating the performance of a single parametrization of the NN may require several hours. This paper addresses the problem of choosing appropriate parameters for the NN by formulating it as a box-constrained mathematical optimization problem, and applying a derivative-free optimization tool that automatically and effectively searches the parameter space. The optimization tool employs a radial basis function model of the objective function (the prediction accuracy of the NN) to accelerate the discovery of configurations yielding high accuracy. Candidate configurations explored by the algorithm are trained to a small number of epochs, and only the most promising candidates receive full training. The performance of the proposed methodology is assessed on benchmark sets and in the context of predicting drug-drug interactions, showing promising results. The optimization tool used in this paper is open-source.


title: A Geometric Approach to Dynamical Model-Order Reduction

abstract: Any model order reduced dynamical system that evolves a modal decomposition to approximate the discretized solution of a stochastic PDE can be related to a vector field tangent to the manifold of fixed rank matrices. The Dynamically Orthogonal (DO) approximation is the canonical reduced order model for which the corresponding vector field is the orthogonal projection of the original system dynamics onto the tangent spaces of this manifold. The embedded geometry of the fixed rank matrix manifold is thoroughly analyzed. The curvature of the manifold is characterized and related to the smallest singular value through the study of the Weingarten map. Differentiability results for the orthogonal projection onto embedded manifolds are reviewed and used to derive an explicit dynamical system for tracking the truncated Singular Value Decomposition (SVD) of a time-dependent matrix. It is demonstrated that the error made by the DO approximation remains controlled under the minimal condition that the original solution stays close to the low rank manifold, which translates into an explicit dependence of this error on the gap between singular values. The DO approximation is also justified as the dynamical system that applies instantaneously the SVD truncation to optimally constrain the rank of the reduced solution. Riemannian matrix optimization is investigated in this extrinsic framework to provide algorithms that adaptively update the best low rank approximation of a smoothly varying matrix. The related gradient flow provides a dynamical system that converges to the truncated SVD of an input matrix for almost every initial data.


title: Internet of Things Security Research: A Rehash of Old Ideas or New Intellectual Challenges?

abstract: The Internet of Things (IoT) is a new computing paradigm that spans wearable devices, homes, hospitals, cities, transportation, and critical infrastructure. Building security into this new computing paradigm is a major technical challenge today. However, what are the security problems in IoT that we can solve using existing security principles? And, what are the new problems and challenges in this space that require new security mechanisms? This article summarizes the intellectual similarities and differences between classic information technology security research and IoT security research.


title: Probing the fluctuations of the optical properties in time-resolved spectroscopy

abstract: We show that, in optical pump-probe experiments on bulk samples, the statistical distribution of the intensity of ultrashort light pulses after the interaction with a nonequilibrium complex material can be used to measure the time-dependent noise of the current in the system. We illustrate the general arguments for a photo-excited Peierls material. The transient noise spectroscopy allows to measure to what extent electronic degrees of freedom dynamically obey the fluctuation-dissipation theorem, and how well they thermalize during the coherent lattice vibrations. The proposed statistical measurement developed here provides a new general framework to retrieve dynamical information on the excited distributions in nonequilibrium experiments which could be extended to other degrees of freedom of magnetic or vibrational origin.


title: Designs for estimating the treatment effect in networks with interference

abstract: In this paper we introduce new, easily implementable designs for drawing causal inference from randomized experiments on networks with interference. Inspired by the idea of matching in observational studies, we introduce the notion of considering a treatment assignment as a quasi-coloring" on a graph. Our idea of a perfect quasi-coloring strives to match every treated unit on a given network with a distinct control unit that has identical number of treated and control neighbors. For a wide range of interference functions encountered in applications, we show both by theory and simulations that the classical Neymanian estimator for the direct effect has desirable properties for our designs. This further extends to settings where homophily is present in addition to interference.


title: Data-driven Random Fourier Features using Stein Effect

abstract: Large-scale kernel approximation is an important problem in machine learning research. Approaches using random Fourier features have become increasingly popular [Rahimi and Recht, 2007], where kernel approximation is treated as empirical mean estimation via Monte Carlo (MC) or Quasi-Monte Carlo (QMC) integration [Yang et al., 2014]. A limitation of the current approaches is that all the features receive an equal weight summing to 1. In this paper, we propose a novel shrinkage estimator from "Stein effect", which provides a data-driven weighting strategy for random features and enjoys theoretical justifications in terms of lowering the empirical risk. We further present an efficient randomized algorithm for large-scale applications of the proposed method. Our empirical results on six benchmark data sets demonstrate the advantageous performance of this approach over representative baselines in both kernel approximation and supervised learning tasks.


title: Model-free causal inference of binary experimental data

abstract: For binary experimental data, we discuss randomization-based inferential procedures that do not need to invoke any modeling assumptions. We also introduce methods for likelihood and Bayesian inference based solely on the physical randomization without any hypothetical super population assumptions about the potential outcomes. These estimators have some properties superior to moment-based ones such as only giving estimates in regions of feasible support. Due to the lack of identification of the causal model, we also propose a sensitivity analysis approach which allows for the characterization of the impact of the association between the potential outcomes on statistical inference.


title: Causal inference for social network data

abstract: We describe semiparametric estimation and inference for causal effects using observational data from a single social network. Our asymptotic results are the first to allow for dependence of each observation on a growing number of other units as sample size increases. In addition, while previous methods have implicitly permitted only one of two possible sources of dependence among social network observations, we allow for both dependence due to transmission of information across network ties and for dependence due to latent similarities among nodes sharing ties. We propose new causal effects that are specifically of interest in social network settings, such as interventions on network ties and network structure. We use our methods to reanalyze an influential and controversial study that estimated causal peer effects of obesity using social network data from the Framingham Heart Study; after accounting for network structure we find no evidence for causal peer effects.


title: Quasi-Random Lattice Model for Electrolyte Solutions: State of Art and Future Perspectives

abstract: In this work, the Quasi-Random Lattice (QRL) model is summarized and critically discussed, in order to outline its potentialities and limitations, in perspective of future developments. QRL primarily focuses on the mean activity coefficient of ionic solutions, the model having first been developed in order to provide practical equations, able to involve a minimal number of unknown or unpredictable quantities. QRL at present depends on one adjustable parameter (at given pressure and temperature), experimentally known for many common salts either symmetric or asymmetric, and corresponding to a well-defined concentration, which also sets the upper limit of applicability of the model. For aqueous electrolytes, the concentration-parameter ranges from 1 M to 8 M (about). In the following it will be seen that, although belonging to the class of simplified approaches, QRL can provide very interesting results since its simple parametrisation is more significant, from a theoretical point of view, than so far recognized. A general overview of the QRL theory will first be presented. Then, some preliminary results will be discussed, in particular concerned with volumetric and thermal properties of electrolyte solutions.


title: Convolution estimates and the number of disjoint partitions

abstract: Let $X$ be a finite collection of sets. We count the number of ways a disjoint union of $n-1$ subsets in $X$ is a set in $X$, and estimate this number from above by $|X|^{c(n)}$ where $$ c(n)=\left(1-\frac{(n-1)\ln (n-1)}{n\ln n} \right)^{-1}. $$ This extends the recent result of Kane-Tao, corresponding to the case $n=3$ where $c(3)\approx 1.725$, to an arbitrary finite number of disjoint $n-1$ partitions.


title: Convergence Analysis of Gradient EM for Multi-component Gaussian Mixture

abstract: In this paper, we study convergence properties of the gradient Expectation-Maximization algorithm \cite{lange1995gradient} for Gaussian Mixture Models for general number of clusters and mixing coefficients. We derive the convergence rate depending on the mixing coefficients, minimum and maximum pairwise distances between the true centers and dimensionality and number of components; and obtain a near-optimal local contraction radius. While there have been some recent notable works that derive local convergence rates for EM in the two equal mixture symmetric GMM, in the more general case, the derivations need structurally different and non-trivial arguments. We use recent tools from learning theory and empirical processes to achieve our theoretical results.


title: Properties of the Intergalactic Magnetic Field Constrained by Gamma-ray Observations of Gamma-Ray Bursts

abstract: The magnetic field in intergalactic space gives important information about magnetogenesis in the early universe. The properties of this field can be probed by searching for radiation of secondary e$^+ $e$^-$ pairs created by TeV photons, that produce GeV range radiation by Compton-scattering cosmic microwave background (CMB) photons. The arrival times of the GeV "echo" photons depend strongly on the magnetic field strength and coherence length. A Monte Carlo code that accurately treats pair creation is developed to simulate the spectrum and time-dependence of the echo radiation. The extrapolation of the spectrum of powerful gamma-ray bursts (GRBs) like GRB 130427A to TeV energies is used to demonstrate how the IGMF can be constrained if it falls in the $10^{-21}$ - $10^{-17}$ G range for 1 Mpc coherence length.


title: Neural network closures for nonlinear model order reduction

abstract: Many reduced order models are neither robust with respect to the parameter changes nor cost-effective enough for handling the nonlinear dependence of complex dynamical systems. In this study, we put forth a robust machine learning framework for projection based reduced order modeling of such nonlinear and nonstationary systems. As a demonstration, we focus on a nonlinear advection-diffusion system given by the viscous Burgers equation, which is a prototype setting of more realistic fluid dynamics applications with the same quadratic nonlinearity. In our proposed methodology the effects of truncated modes are modeled using a singe layer feed-forward neural network architecture. The neural network architecture is trained by utilizing both the Bayesian regularization and extreme learning machine approaches, where the latter one is found to be computationally more efficient. A particular effort is devoted to the selection of basis functions considering the proper orthogonal decomposition and Fourier bases. It is shown that the proposed models yield significant improvements in the accuracy over the standard Galerkin projection models with a negligibly small computational overhead and provide reliable predictions with respect to the parameter changes.


title: Experimental data over quantum mechanics simulations for inferring the repulsive exponent of the Lennard-Jones potential in Molecular Dynamics

abstract: The Lennard-Jones (LJ) potential is a cornerstone of Molecular Dynamics (MD) simulations and among the most widely used computational kernels in science. The potential models atomistic attraction and repulsion with century old prescribed parameters ($q=6, \; p=12$, respectively), originally related by a factor of two for simplicity of calculations. We re-examine the value of the repulsion exponent through data driven uncertainty quantification. We perform Hierarchical Bayesian inference on MD simulations of argon using experimental data of the radial distribution function (RDF) for a range of thermodynamic conditions, as well as dimer interaction energies from quantum mechanics simulations. The experimental data suggest a repulsion exponent ($p \approx 6.5$), in contrast to the quantum simulations data that support values closer to the original ($p=12$) exponent. Most notably, we find that predictions of RDF, diffusion coefficient and density of argon are more accurate and robust in producing the correct argon phase around its triple point, when using the values inferred from experimental data over those from quantum mechanics simulations. The present results suggest the need for data driven recalibration of the LJ potential across MD simulations.


title: High voltage assisted mechanical stabilization of single-molecule junctions

abstract: The realization of molecular-based electronic devices depends to a large extent on the ability to mechanically stabilize the involved molecular bonds, while making use of efficient resonant charge transport through the device. Resonant charge transport can induce vibrational instability of molecular bonds, leading to bond rupture under a bias voltage. In this work, we go beyond the wide-band approximation in order to study the phenomenon of vibrational instability in single molecule junctions and show that the energy-dependence of realistic molecule-leads couplings affects the mechanical stability of the junction. We show that the chemical bonds can be stabilized in the resonant transport regime by increasing the bias voltage on the junction. This research provides guidelines for the design of mechanically stable molecular devices operating in the regime of resonant charge transport.


title: Simulating electric field interactions with polar molecules using spectroscopic databases

abstract: Ro-vibrational Stark-associated phenomena of small polyatomic molecules are modelled using extensive spectroscopic data generated as part of the ExoMol project. The external field Hamiltonian is built from the computed ro-vibrational line list of the molecule in question. The Hamiltonian we propose is general and suitable for any polar molecule in the presence of an electric field. By exploiting precomputed data, the often prohibitively expensive computations associated with high accuracy simulations of molecule-field interactions are avoided. Applications to strong terahertz field-induced ro-vibrational dynamics of PH$_3$ and NH$_3$, and spontaneous emission data for optoelectrical Sisyphus cooling of H$_2$CO and CH$_3$Cl are discussed.


title: A Quantum-like Model of Selection Behavior

abstract: In this paper, we introduce a new model of selection behavior under risk that describes an essential cognitive process for comparing values of objects and making a selection decision. This model is constructed by the quantum-like approach that employs the state representation specific to quantum theory, which has the mathematical framework beyond the classical probability theory. We show that our quantum approach can clearly explain the famous examples of anomalies for the expected utility theory, the Ellsberg paradox, the Machina paradox and the disparity between WTA and WTP. Further, we point out that our model mathematically specifies the characteristics of the probability weighting function and the value function, which are basic concepts in the prospect theory.


title: The first Baade-Wesselink analysis of Blazhko RR Lyrae stars: discrepancies between photometrically and spectroscopically determined radius variations

abstract: The simultaneous photometric and spectroscopic observations of the RR Lyrae variables in the globular cluster, M3, published in Jurcsik et al. (2017, Paper I) made it possible to perform Baade-Wesselink (BW) analysis of a large sample of Blazhko stars for the first time. The BW distances of Blazhko stars turned out to be unreliable, as significantly different distances were obtained for the stars of the Blazhko sample and also for the same star in different modulation phases. Even the results of small modulation-amplitude Blazhko stars may be doubtful. This result warns that the application of the BW method to Blazhko stars is not trustworthy. Keeping the distance fixed for each Blazhko star in each modulation phase, a significant difference between the spectroscopic and the photometric radius ($R_{\textrm{sp}}$, $R_{\textrm{ph}}$) variations is detected. The phase and amplitude variations of $R_{\textrm{sp}}$ follow the changes of the light curve during the Blazhko cycle but the $R_{\textrm{ph}}$ curve seems to be not (or only marginally) affected by the modulation. The asynchronous behaviour of $R_{\textrm{sp}}$ and $R_{\textrm{ph}}$ supports the interpretation of the Blazhko effect as a depth-dependent phenomenon, as the spectroscopic radius variation reflects the radial displacement of the line-forming region high in the atmosphere, while the photospheric radius variation is derived from the information of the observed visual-band light emitted mostly by the lower photosphere. The stability of $R_{\textrm{ph}}$ may be interpreted as a strong argument against the non-radial-mode explanation of the Blazhko phenomenon.


title: Effect of finite particle number sampling on baryon number fluctuations

abstract: The effects of finite particle number sampling on the net baryon number cumulants, extracted from fluid dynamical simulations, are studied. The commonly used finite particle number sampling procedure introduces an additional Poissonian (or multinomial if global baryon number conservation is enforced) contribution which increases the extracted moments of the baryon number distribution. If this procedure is applied to a fluctuating fluid dynamics framework one severely overestimates the actual cumulants. We show that the sampling of so called test-particles suppresses the additional contribution to the moments by at least one power of the number of test-particles. We demonstrate this method in a numerical fluid dynamics simulation that includes the effects of spinodal decomposition due to a first order phase transition. Furthermore, in the limit where anti-baryons can be ignored, we derive analytic formulas which capture exactly the effect of particle sampling on the baryon number cumulants. These formulas may be used to test the various numerical particle sampling algorithms.


title: Conscious and controlling elements in combinatorial group testing problems with more defectives

abstract: In combinatorial group testing problems Questioner needs to find a defective element $x\in [n]$ by testing subsets of $[n]$. In [18] the authors introduced a new model, where each element knows the answer for those queries that contain it and each element should be able to identify the defective one. In this article we continue to investigate this kind of models with more defective elements. We also consider related models inspired by secret sharing models, where the elements should share information among them to find out the defectives. Finally the adaptive versions of the different models are also investigated.


title: Critical two-point function for long-range $O(n)$ models below the upper critical dimension

abstract: We consider the $n$-component $|\varphi|^4$ lattice spin model ($n \ge 1$) and the weakly self-avoiding walk ($n=0$) on $\mathbb{Z}^d$, in dimensions $d=1,2,3$. We study long-range models based on the fractional Laplacian, with spin-spin interactions or walk step probabilities decaying with distance $r$ as $r^{-(d+\alpha)}$ with $\alpha \in (0,2)$. The upper critical dimension is $d_c=2\alpha$. For $\epsilon >0$, and $\alpha = \frac 12 (d+\epsilon)$, the dimension $d=d_c-\epsilon$ is below the upper critical dimension. For small $\epsilon$, weak coupling, and all integers $n \ge 0$, we prove that the two-point function at the critical point decays with distance as $r^{-(d-\alpha)}$. This "sticking" of the critical exponent at its mean-field value was first predicted in the physics literature in 1972. Our proof is based on a rigorous renormalisation group method. The treatment of observables differs from that used in recent work on the nearest-neighbour 4-dimensional case, via our use of a cluster expansion.


title: Global existence in the 1D quasilinear parabolic-elliptic chemotaxis system with critical nonlinearity

abstract: The paper should be viewed as complement of an earlier result in [8]. In the paper just mentioned it is shown that 1d case of a quasilinear parabolic-elliptic Keller-Segel system is very special. Namely, unlike in higher dimensions, there is no critical nonlinearity. Indeed, for the nonlinear diffusion of the form 1/u all the solutions, independently on the magnitude of initial mass, stay bounded. However, the argument presented in [8] deals with the Jager-Luckhaus type system. And is very sensitive to this restriction. Namely, the change of variables introduced in [8], being a main step of the method, works only for the Jager-Luckhaus modification. It does not seem to be applicable in the usual version of the parabolic-elliptic Keller-Segel system. The present paper fulfils this gap and deals with the case of the usual parabolic-elliptic version. To handle it we establish a new Lyapunov-like functional (it is related to what was done in [8]), which leads to global existence of the initial-boundary value problem for any initial mass.


title: Pseudogap and Fermi surface in the presence of spin-vortex checkerboard for 1/8-doped lanthanum cuprates

abstract: Lanthanum family of high-temperature cuprate superconductors is known to exhibit both spin and charge electronic modulations around doping level 1/8. We assume that these modulations have the character of two-dimensional spin-vortex checkerboard and investigate whether this assumption is consistent with the Fermi surface and the pseudogap measured by angle-resolved photo-emission spectroscopy. We also explore the possibility of observing quantum oscillations of transport coefficients in such a background. These investigations are based on a model of non-interacting spin-1/2 fermions hopping on a square lattice and coupled through spins to a magnetic field imitating spin-vortex checkerboard. The main results of this article include (i) calculation of Fermi surface containing Fermi arcs at the positions in the Brillouin zone largely consistent with experiments; (ii) identification of factors complicating the observations of quantum oscillations in the presence of spin modulations; and (iii) investigation of the symmetries of the resulting electronic energy bands, which, in particular, indicates that each band is double-degenerate and, in addition, has at least one conical point, where it touches another double-degenerate band. We discuss possible implications these cones may have for the transport properties and the pseudogap.


title: Observing Pulsars with a Phased Array Feed at the Parkes Telescope

abstract: During February 2016, CSIRO Astronomy and Space Science and the Max-Planck-Institute for Radio Astronomy installed, commissioned and carried out science observations with a phased array feed (PAF) receiver system on the 64m diameter Parkes radio telescope. Here we demonstrate that the PAF can be used for pulsar observations and we highlight some unique capabilities. We demonstrate that the pulse profiles obtained using the PAF can be calibrated and that multiple pulsars can be simultaneously observed. Significantly, we find that an intrinsic polarisation leakage of -31dB can be achieved with a PAF beam offset from the centre of the field of view. We discuss the possibilities for using a PAF for future pulsar observations and for searching for fast radio bursts with the Parkes and Effelsberg telescopes.


title: Data Visualization on Day One: Bringing Big Ideas into Intro Stats Early and Often

abstract: In a world awash with data, the ability to think and compute with data has become an important skill for students in many fields. For that reason, inclusion of some level of statistical computing in many introductory-level courses has grown more common in recent years. Existing literature has documented multiple success stories of teaching statistics with R, bolstered by the capabilities of R Markdown. In this article, we present an in-class data visualization activity intended to expose students to R and R Markdown during the first week of an introductory statistics class. The activity begins with a brief lecture on exploratory data analysis in R. Students are then placed in small groups tasked with exploring a new dataset to produce three visualizations that describe particular insights that are not immediately obvious from the data. Upon completion, students will have produced a series of univariate and multivariate visualizations on a real dataset and practiced describing them.


title: Financial Time Series Forecasting: Semantic Analysis Of Economic News

abstract: The paper proposes a method of financial time series forecasting taking into account the semantics of news. For the semantic analysis of financial news the sampling of negative and positive words in economic sense was formed based on Loughran McDonald Master Dictionary. The sampling included the words with high frequency of occurrence in the news of financial markets. For single-root words it has been left only common part that allows covering few words for one request. Neural networks were chosen for modeling and forecasting. To automate the process of extracting information from the economic news a script was developed in the MATLAB Simulink programming environment, which is based on the generated sampling of positive and negative words. Experimental studies with different architectures of neural networks showed a high adequacy of constructed models and confirmed the feasibility of using information from news feeds to predict the stock prices.


title: On factorizations of graphical maps

abstract: We study the categories governing infinity (wheeled) properads. The graphical category, which was already known to be generalized Reedy, is in fact an Eilenberg-Zilber category. A minor alteration to the definition of the wheeled graphical category allows us to show that it is a generalized Reedy category. Finally, we present model structures for Segal properads and Segal wheeled properads.


title: SPH Modeling of Short-crested Waves

abstract: This study investigates short-crested wave breaking over a planar beach by using the mesh-free Smoothed Particle Hydrodynamics model, GPUSPH. The short-crested waves are created by generating intersecting wave trains in a numerical wave basin. We examine the influence of beach slope, incident wave height, and incident wave angle on the generated short-crested waves. Short-crested wave breaking over a steeper beach generates stronger rip currents, and larger circulation cells in front of the beach. Intersecting wave trains with a larger incident wave height drive a more complicated short-crested wave field including isolated breakers and wave amplitude diffraction. Nearshore circulation induced by short-crested wave breaking is greatly influenced by the incident wave angle (or the rip current spacing). There is no secondary circulation cell between the nodal line and the antinodal line if the rip current spacing is narrow. However, there are multiple secondary circulation cells observed when the rip current spacing is relatively large.


title: Detecting the growth of structures in Pure Stellar Disk Models

abstract: We performed a series of 3D N-body simulations where the initial conditions were chosen to get two sets of models; unbarred and barred ones. In this work, we analyze the growth of spirals and bar structures using 1D, and 2D Fourier Transforms FT methods. Spectrograms and diagrams of the amplitude of the Fourier coefficients as a function of time, radius and pitch angle show that the general morphology, of our modeled galaxies, is due to the superposition of structures which have different values of pitch angle and number of arms. Also, we made in barred models a geometric classification of orbits from the bar reference frame showing that the barred potential and the Lagrangian points $L_4$ and $L_5$ catch approximately one-third of the total disk mass.


title: Self-Organized Supercriticality and Oscillations in Networks of Stochastic Spiking Neurons

abstract: Networks of stochastic spiking neurons are interesting models in the area of Theoretical Neuroscience, presenting both continuous and discontinuous phase transitions. Here we study fully connected networks analytically, numerically and by computational simulations. The neurons have dynamic gains that enable the network to converge to a stationary slightly supercritical state (self-organized supercriticality or SOSC) in the presence of the continuous transition. We show that SOSC, which presents power laws for neuronal avalanches plus some large events, is robust as a function of the main parameter of the neuronal gain dynamics. We discuss the possible applications of the idea of SOSC to biological phenomena like epilepsy and dragon king avalanches. We also find that neuronal gains can produce collective oscillations that coexists with neuronal avalanches, with frequencies compatible with characteristic brain rhythms.


title: Deep Multi-instance Networks with Sparse Label Assignment for Whole Mammogram Classification

abstract: Mammogram classification is directly related to computer-aided diagnosis of breast cancer. Traditional methods rely on regions of interest (ROIs) which require great efforts to annotate. Inspired by the success of using deep convolutional features for natural image analysis and multi-instance learning (MIL) for labeling a set of instances/patches, we propose end-to-end trained deep multi-instance networks for mass classification based on whole mammogram without the aforementioned ROIs. We explore three different schemes to construct deep multi-instance networks for whole mammogram classification. Experimental results on the INbreast dataset demonstrate the robustness of proposed networks compared to previous work using segmentation and detection annotations.


title: Safe Model-based Reinforcement Learning with Stability Guarantees

abstract: Reinforcement learning is a powerful paradigm for learning optimal policies from experimental data. However, to find optimal policies, most reinforcement learning algorithms explore all possible actions, which may be harmful for real-world systems. As a consequence, learning algorithms are rarely applied on safety-critical systems in the real world. In this paper, we present a learning algorithm that explicitly considers safety, defined in terms of stability guarantees. Specifically, we extend control-theoretic results on Lyapunov stability verification and show how to use statistical models of the dynamics to obtain high-performance control policies with provable stability certificates. Moreover, under additional regularity assumptions in terms of a Gaussian process prior, we prove that one can effectively and safely collect data in order to learn about the dynamics and thus both improve control performance and expand the safe region of the state space. In our experiments, we show how the resulting algorithm can safely optimize a neural network policy on a simulated inverted pendulum, without the pendulum ever falling down.


title: Path-sum solution of the Weyl Quantum Walk in 3+1 dimensions

abstract: We consider the Weyl quantum walk in 3+1 dimensions, that is a discrete-time walk describing a particle with two internal degrees of freedom moving on a Cayley graph of the group $\mathbb Z^3$, that in an appropriate regime evolves according to Weyl's equation. The Weyl quantum walk was recently derived as the unique unitary evolution on a Cayley graph of $\mathbb Z^3$ that is homogeneous and isotropic. The general solution of the quantum walk evolution is provided here in the position representation, by the analytical expression of the propagator, i.e. transition amplitude from a node of the graph to another node in a finite number of steps. The quantum nature of the walk manifests itself in the interference of the paths on the graph joining the given nodes. The solution is based on the binary encoding of the admissible paths on the graph and on the semigroup structure of the walk transition matrices.


title: Lieb-Robinson bounds, the spectral flow, and stability of the spectral gap for lattice fermion systems

abstract: We prove Lieb-Robinson bounds for a general class of lattice fermion systems. By making use of a suitable conditional expectation onto subalgebras of the CAR algebra, we can apply the Lieb-Robinson bounds much in the same way as for quantum spin systems. We preview how to obtain the spectral flow automorphisms and to prove stability of the spectral gap for frustration-free gapped systems satisfying a Local Topological Quantum Order condition.


title: On beam models and their paraxial approximation

abstract: We derive focused laser pulse solutions to the electromagnetic wave equation in vacuum. After reproducing beam and pulse expressions for the well-known paraxial Gaussian and axicon cases, we apply the method to analyse a laser beam with Lorentzian transverse momentum distribution. Whilst a paraxial approach has some success close to the focal axis and within a Rayleigh range of the focal spot, we find that it incorrectly predicts the transverse fall-off typical of a Lorentzian. Our vector-potential approach is particularly relevant to calculation of quantum electrodynamical processes in weak laser pulse backgrounds.


title: Statistical study of solar dimmings using CoDiT

abstract: We present the results from analyzing the physical and morphological properties of 154 dimmings (transient coronal holes) and the associated flares and coronal mass ejections (CMEs). Each dimming in our 2013 catalog was processed with the semi-automated Coronal Dimming Tracker using Solar Dynamics Observatory AIA 193 A observations and HMI magnetograms. Instead of the typically used difference images, we used our coronal hole detection algorithm to detect transient dark regions "directly" in extreme ultraviolet (EUV) images. This allowed us to study dimmings as the footpoints of CMEs - in contrast with the larger, diffuse dimmings seen in difference images that represent the projected view of the rising, expanding plasma. Studying the footpoint-dimming morphology allowed us to better understand the CME structure in the low corona. While comparing the physical properties of dimmings, flares, and CMEs we were also able to identify relationships between the different parts of this complex eruptive phenomenon. We found that larger dimmings are longer-lived, suggesting that it takes longer to "close down" large open magnetic regions. Also, during their growth phase, smaller dimmings acquire a higher magnetic flux imbalance (i. e., become more unipolar) than larger dimmings. Furthermore, we found that the EUV intensity of dimmings (indicative of local electron density) correlates with how much plasma was removed and how energetic the eruption was. Studying the morphology of dimmings (single, double, fragmented) also helped us identify different configurations of the quasi-open magnetic field.


title: Constraining the Orbit of Supermassive Black Hole Binary 0402+379

abstract: The radio galaxy 0402+379 is believed to host a supermassive black hole binary (SMBHB). The two compact core sources are separated by a projected distance of 7.3 pc, making it the most (spatially) compact resolved SMBHB known. We present new multi-frequency VLBI observations of 0402+379 at 5, 8, 15 and 22 GHz, and combine with previous observations spanning 12 years. A strong frequency dependent core shift is evident, which we use to infer magnetic fields near the jet base. After correcting for these shifts we detect significant relative motion of the two cores at $\beta=v/c=0.0054 \pm 0.0003$ at $PA= -34.4^\circ$. With some assumptions about the orbit, we use this measurement to constrain the orbital period $P\approx 3 \times 10^4$ y and SMBHB mass $M \approx 15 \times 10^9\ M_\odot$. While additional observations are needed to confirm this motion and obtain a precise orbit, this is apparently the first black hole system resolved as a visual binary.


title: Grounded Recurrent Neural Networks

abstract: In this work, we present the Grounded Recurrent Neural Network (GRNN), a recurrent neural network architecture for multi-label prediction which explicitly ties labels to specific dimensions of the recurrent hidden state (we call this process "grounding"). The approach is particularly well-suited for extracting large numbers of concepts from text. We apply the new model to address an important problem in healthcare of understanding what medical concepts are discussed in clinical text. Using a publicly available dataset derived from Intensive Care Units, we learn to label a patient's diagnoses and procedures from their discharge summary. Our evaluation shows a clear advantage to using our proposed architecture over a variety of strong baselines.


title: Nature of Kinetic Scale Turbulence in the Earth's Magnetosheath

abstract: We present a combined observational and theoretical analysis to investigate the nature of plasma turbulence at kinetic scales in the Earth's magnetosheath. In the first decade of the kinetic range, just below the ion gyroscale, the turbulence was found to be similar to that in the upstream solar wind: predominantly anisotropic, low-frequency and kinetic Alfv\'en in nature. A key difference, however, is that the magnetosheath ions are typically much hotter than the electrons, $T_\mathrm{i}\gg T_\mathrm{e}$, which, together with $\beta_\mathrm{i}\sim 1$, leads to a change in behaviour in the second decade, close to electron scales. The turbulence here is characterised by an increased magnetic compressibility, following a mode we term the inertial kinetic Alfv\'en wave, and a steeper spectrum of magnetic fluctuations, consistent with the prediction $E_B(k_\perp)\propto k_\perp^{-11/3}$ that we obtain from a set of nonlinear equations. This regime of plasma turbulence may also be relevant for other astrophysical environments with $T_\mathrm{i}\gg T_\mathrm{e}$, such as the solar corona, hot accretion flows, and regions downstream of collisionless shocks.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.plasm-ph


title: Random ordering formula for sofic and Rokhlin entropy of Gibbs measures

abstract: We prove the explicit formula for sofic and Rokhlin entropy of actions arising from some class of Gibbs measures. It provides a new set of examples with sofic entropy independent of sofic approximations. It is particularilly interresting, since in non-amenable case Rokhlin entropy was computed only in case of Bernoulli actions and for some examples with zero Rokhlin entropy. As an example we show that our formula holds for the supercritical Ising model. We also establish a criterion for uniqueness of Gibbs measure by means of f-invariant.


title: Singularity development and supersymmetry in holography

abstract: We study the effects of supersymmetry on singularity development scenario in holography presented in \cite{Bosch:2017ccw} (BBL). We argue that the singularity persists in a supersymmetric extension of the BBL model. The challenge remains to find a string theory embedding of the singularity mechanism.


title: A Taylor expansion of the square root matrix functional

abstract: This short note provides an explicit description of the Fr\'echet derivatives of the principal square root matrix functional at any order. We present an original formulation that allows to compute sequentially the Fr\'echet derivatives of the matrix square root at any order starting from the first order derivative. A Taylor expansion at any order with an integral remainder term is also provided, yielding the first result of this type for this class of matrix functional.


title: Hashing as Tie-Aware Learning to Rank

abstract: Hashing, or learning binary embeddings of data, is frequently used in nearest neighbor retrieval. In this paper, we develop learning to rank formulations for hashing, aimed at directly optimizing ranking-based evaluation metrics such as Average Precision (AP) and Normalized Discounted Cumulative Gain (NDCG). We first observe that the integer-valued Hamming distance often leads to tied rankings, and propose to use tie-aware versions of AP and NDCG to evaluate hashing for retrieval. Then, to optimize tie-aware ranking metrics, we derive their continuous relaxations, and perform gradient-based optimization with deep neural networks. Our results establish the new state-of-the-art for image retrieval by Hamming ranking in common benchmarks.


title: Simple Pricing Schemes for the Cloud

abstract: The problem of pricing the cloud has attracted much recent attention due to the widespread use of cloud computing and cloud services. From a theoretical perspective, several mechanisms that provide strong efficiency or fairness guarantees and desirable incentive properties have been designed. However, these mechanisms often rely on a rigid model, with several parameters needing to be precisely known in order for the guarantees to hold. In this paper, we consider a stochastic model and show that it is possible to obtain good welfare and revenue guarantees with simple mechanisms that do not make use of the information on some of these parameters. In particular, we prove that a mechanism that sets the same price per time step for jobs of any length achieves at least 50% of the welfare and revenue obtained by a mechanism that can set different prices for jobs of different lengths, and the ratio can be improved if we have more specific knowledge of some parameters. Similarly, a mechanism that sets the same price for all servers even though the servers may receive different kinds of jobs can provide a reasonable welfare and revenue approximation compared to a mechanism that is allowed to set different prices for different servers.


title: Towards Interrogating Discriminative Machine Learning Models

abstract: It is oftentimes impossible to understand how machine learning models reach a decision. While recent research has proposed various technical approaches to provide some clues as to how a learning model makes individual decisions, they cannot provide users with ability to inspect a learning model as a complete entity. In this work, we propose a new technical approach that augments a Bayesian regression mixture model with multiple elastic nets. Using the enhanced mixture model, we extract explanations for a target model through global approximation. To demonstrate the utility of our approach, we evaluate it on different learning models covering the tasks of text mining and image recognition. Our results indicate that the proposed approach not only outperforms the state-of-the-art technique in explaining individual decisions but also provides users with an ability to discover the vulnerabilities of a learning model.


title: Factorization systems on (stable) derivators

abstract: We define triangulated factorization systems on triangulated categories, and prove that a suitable subclass thereof (the normal triangulated torsion theories) corresponds bijectively to $t$-structures on the same category. This result is then placed in the framework of derivators regarding a triangulated category as the base of a stable derivator. More generally, we define derivator factorization systems in the 2-category $\mathrm{PDer}$, describing them as algebras for a suitable strict 2-monad (this result is of independent interest), and prove that a similar characterization still holds true: for a stable derivator $\mathbb D$, a suitable class of derivator factorization systems (the normal derivator torsion theories) correspond bijectively with $t$-structures on the base $\mathbb{D}(\mathbb{1})$ of the derivator. These two result can be regarded as the triangulated- and derivator- analogues, respectively, of the theorem that says that `$t$-structures are normal torsion theories' in the setting of stable $\infty$-categories, showing how the result remains true whatever the chosen model for stable homotopy theory is.


title: A Near-Optimal Separation Principle for Nonlinear Stochastic Systems Arising in Robotic Path Planning and Control

abstract: We consider nonlinear stochastic systems that arise in path planning and control of mobile robots. As is typical of almost all nonlinear stochastic systems, the optimally solving problem is intractable. We provide a design approach which yields a tractable design that is quantifiably near-optimal. We exhibit a "separation" principle under a small noise assumption consisting of the optimal open-loop design of nominal trajectory followed by an optimal feedback law to track this trajectory, which is different from the usual effort of separating estimation from control. As a corollary, we obtain a trajectory-optimized linear quadratic regulator design for stochastic nonlinear systems with Gaussian noise.


title: Variable stars in the northern Galactic plane from KISOGP

abstract: We have conducted a large-scale survey of the northern plane using Kiso Wide Field Camera attached to Schmidt telescope at Kiso observatory. The KISOGP (KWFC Intensive Survey of the Galactic Plane) project have made 40-70 epoch observations in I band of about 320 sq. degrees for 5 years starting in 2012. The limiting magnitude is around 16.5 in I. In the data analysis so far, we detected a couple of thousands of variable stars including approximately 100 Cepheids and more than 700 Miras. Roughly 90 percent of them were not previously reported as variable stars, indicating that there are still many relatively bright variables to be found in the Galactic plane.


title: The Future of Ad Blocking: An Analytical Framework and New Techniques

abstract: We present a systematic study of ad blocking - and the associated "arms race" - as a security problem. We model ad blocking as a state space with four states and six state transitions, which correspond to techniques that can be deployed by either publishers or ad blockers. We argue that this is a complete model of the system. We propose several new ad blocking techniques, including ones that borrow ideas from rootkits to prevent detection by anti-ad blocking scripts. Another technique uses the insight that ads must be recognizable by humans to comply with laws and industry self-regulation. We have built prototype implementations of three of these techniques, successfully blocking ads and evading detection. We systematically evaluate our proposed techniques, along with existing ones, in terms of security, practicality, and legality. We characterize the order of growth of the development effort required to create/maintain ad blockers as a function of the growth of the web. Based on our state-space model, our new techniques, and this systematization, we offer insights into the likely "end game" of the arms race. We challenge the widespread assumption that the arms race will escalate indefinitely, and instead identify a combination of evolving technical and legal factors that will determine the outcome.


title: Slavnov Identities for Infinite Derivative Theory of Gravity

abstract: In this paper, we shall consider an infinite-derivative theory of gravity, with a view to making it renormalisable. First, we derive the modified superficial degree of divergence. Next, we establish that the theory is invariant under BRS (Becchi-Rouet-Stora) transformations. Consequently, we show that the corresponding Slavnov identities hold for the theory. Finally, we summarise our results and discuss renormalisability.


title: Microstructured Stark Decelerator

abstract: A Stark decelerator is an effective tool for controlling motional degrees of freedom of polar molecules. Due to technical limitations, many of the current Stark decelerators focus on molecules in low-field-seeking quantum states and are built based on a fixed electrode size and spacing (type-A architecture). Here, we present an alternative method based on a different architecture, so-called type-B, with microstructured electrodes and a simpler electric field pulse timing with the prospect of producing cold, quasi-CW molecular beams. We demonstrate with simulations the feasibility of this method and show that a 2 cm decelerator consisting of 100 stages can bring molecules from 200 m/s to a near standstill in about 150 microseconds. Subsequently, the molecules can be reaccelerated to 200 m/s in a phase-stable manner. Our simulations focus on the longitudinal transmission leaving the problem of a transverse confinement for a future study. This work points to a possible way of integrating polar molecule accelerators with solid-state electronics.


title: The observed $\Omega_{c}^{0}$ resonances as pentaquark states

abstract: In present work, we investigate the spectrum of several low-lying $sscq\bar{q}$ pentaquark configurations employing the constituent quark model, within which the hyperfine interaction between quarks is taken to be mediated by Goldstone boson exchange. Our numerical results show that four $sscq\bar{q}$ configurations with $J^{P}=1/2^{-}$ or $J^{P}=3/2^{-}$ lie at energies very close to the recently observed five $\Omega_{c}^{0}$ states by LHCb collaboration, this indicates that the $sscq\bar{q}$ pentaquark configurations may form sizable components of the observed $\Omega_{c}^{0}$ resonances.


title: Joint Rate Control and Power Allocation for Non-Orthogonal Multiple Access Systems

abstract: This paper investigates the optimal resource allocation of a downlink non-orthogonal multiple access (NOMA) system consisting of one base station and multiple users. Unlike existing short-term NOMA designs that focused on the resource allocation for only the current transmission timeslot, we aim to maximize a long-term network utility by jointly optimizing the data rate control at the network layer and the power allocation among multiple users at the physical layer, subject to practical constraints on both the short-term and long-term power consumptions. To solve this problem, we leverage the recently-developed Lyapunov optimization framework to convert the original long-term optimization problem into a series of online rate control and power allocation problems in each timeslot. The power allocation problem, however, is shown to be non-convex in nature and thus cannot be solved with a standard method. However, we explore two structures of the optimal solution and develop a dynamic programming based power allocation algorithm, which can derive a globally optimal solution, with a polynomial computational complexity. Extensive simulation results are provided to evaluate the performance of the proposed joint rate control and power allocation framework for NOMA systems, which demonstrate that the proposed NOMA design can significantly outperform multiple benchmark schemes, including orthogonal multiple access (OMA) schemes with optimal power allocation and NOMA schemes with non-optimal power allocation, in terms of average throughput and data delay.


title: Vortices enhance diffusion in dense granular flows

abstract: This Letter introduces unexpected diffusion properties in dense granular flows, and shows that they result from the development of partially jammed clusters of grains, or granular vortices. Transverse diffusion coefficients $D$ and average vortex sizes $\ell$ are systematically measured in simulated plane shear flows at differing internal numbers $I$ revealing (i) a strong deviation from the expected scaling $D\propto d^2 \dot \gamma$ involving the grain size $d$ and shear rate $\dot \gamma$ and (ii) an increase in average vortex size $\ell$ at low $I$, following $\ell\propto dI^{-\frac{1}{2}}$ but limited by the system size. A general scaling $D\propto \ell d \dot \gamma $ is introduced that captures all the measurements and highlights the key role of vortex size. This leads to establishing a scaling for the diffusivity in dense granular flow as $D\propto d^2 \sqrt{\dot \gamma/ t_i}$ involving the geometric average of shear time $1/\dot\gamma$ and inertial time $t_i$ as the relevant time scale. Analysis of grain trajectories further evidence that this diffusion process arises from a vortex-driven random walk.


title: A Gromov hyperbolic metric vs the hyperbolic and other related metrics

abstract: We mainly consider two metrics: a Gromov hyperbolic metric and a scale invariant Cassinian metric. We compare these two metrics and obtain their relationship with certain well-known hyperbolic-type metrics, leading to several inclusion relations between the associated metric balls.


title: Transient transition from free carrier metallic state to exciton insulating state in GaAs by ultrafast photoexcitation

abstract: We present systematic studies of the transient dynamics of GaAs by ultrafast time-resolved optical reflectivity. In photo excited non-equilibrium states, we found a sign reverse in transient reflectivity spectra $\Delta R/R$ (t $>$ 0), from positive around room temperature to negative at cryogenic temperatures. The former corresponds to a transient free carrier metallic state, while the latter is attributed to an exciton insulating state, in which the transient electronic properties is mostly dominated by excitons, resulting in a transient metal-insulator transition (MIT). Two transition temperatures (T$_1$ and T$_2$) are well identified by analysing the intensity change of the time-resolved optical spectra. We found that photoexcited MIT starts emerging at T$_1$ as high as $\sim$ 230 K, in terms of a negative dip feature at 0.4 ps, and becomes stabilized below T$_2$ associated with a negative constant after 40 ps in spectra. Our results address a phase diagram that provides a framework for MIT through temperature and photoexcitation, and shed light on the understanding of light-semiconductor interaction and exciton physics.


title: Flexible Cache-Aided Networks with Backhauling

abstract: Caching at the edge is a promising technique to cope with the increasing data demand in wireless networks. This paper analyzes the performance of cellular networks consisting of a tier macro-cell wireless backhaul nodes overlaid with a tier of cache-aided small cells. We consider both static and dynamic association policies for content delivery to the user terminals and analyze their performance. In particular, we derive closed-form expressions for the area spectral efficiency and the energy efficiency, which are used to optimize relevant design parameters such as the density of cache-aided small cells and the storage size. By means of this approach, we are able to draw useful design insights for the deployment of highly performing cache-aided tiered networks.


title: Optimized contraction scheme for tensor-network states

abstract: In the tensor-network framework, the expectation values of two-dimensional quantum states are evaluated by contracting a double-layer tensor network constructed from initial and final tensor-network states. The computational cost of carrying out this contraction is generally very high, which limits the largest bond dimension of tensor-network states that can be accurately studied to a relatively small value. We propose an optimized contraction scheme to solve this problem by mapping the double-layer tensor network onto an intersected single-layer tensor network. This reduces greatly the bond dimensions of local tensors to be contracted and improves dramatically the efficiency and accuracy of the evaluation of expectation values of tensor-network states. It almost doubles the largest bond dimension of tensor-network states whose physical properties can be efficiently and reliably calculated, and it extends significantly the application scope of tensor-network methods.

categories: cond-mat.str-el cond-mat.stat-mech physics.comp-ph quant-ph