title: OPE and a low-energy theorem in QCD-like theories

abstract: We verify, both perturbatively and nonperturbatively asymptotically in the ultraviolet (UV), a special case of a low-energy theorem of the NSVZ type in QCD-like theories, recently derived in arXiv:1701.07833, that relates the logarithmic derivative with respect to the gauge coupling, or the logarithmic derivative with respect to the renormalization-group (RG) invariant scale, of an $n$-point correlator of local operators in one side to an $n+1$-point correlator with the insertion of $Tr F^2$ at zero momentum in the other side. Our computation involves the operator product expansion (OPE) of the scalar glueball operator, $Tr F^2$, in massless QCD, worked out perturbatively in arXiv:1209.1516 -- and in its RG-improved form in the present paper -- by means of which we extract both the perturbative divergences and the nonperturbative UV asymptotics in both sides. We also discuss the role of the contact terms in the OPE, both finite and divergent, discovered some years ago in arXiv:1209.1516, in relation to the low-energy theorem. Besides, working the other way around by assuming the low-energy theorem for any 2-point correlator of a multiplicatively renormalizable gauge-invariant operator, we compute in a massless QCD-like theory the corresponding perturbative OPE to the order of $g^2$ and nonperturbative asymptotics. The low-energy theorem has a number of applications: to the renormalization in asymptotically free QCD-like theories, both perturbatively and nonperturbatively in the large-$N$ 't Hooft and Veneziano expansions, and to the way the open/closed string duality may or may not be realized in the would-be solution by canonical string theories for QCD-like theories, both perturbatively and in the 't Hooft large-$N$ expansion. Our computations will also enter further developments based on the low-energy theorem.


title: Hypothesis about Enrichment of Solar System

abstract: Despite significant progress in the understanding of galactic nucleosynthesis and its influence on the solar system neighborhood, challenges remain in the understanding of enrichment of the solar system itself. Based on the detailed review of multi-disciplinary literature, we propose a scenario that an event of nucleogenesis -- not nucleosynthesis (from lower nucleon numbers A to higher A) but nuclear-fission (from higher A to lower A) -- occurred in the inner part of the solar system at one of the stages of its evolution. We propose a feasible mechanism of implementation of such event. The occurrence of such event could help explain the puzzles in yet-unresolved isotopic abundances, certain meteoritic anomalies, as well as peculiarities in the solar system's composition and planetary structure. We also discuss experimental data and available results from existing models (in several relevant sub-fields) that provide support and/or appear consistent with the hypothesis.


title: Normality of Necessary Optimality Conditions for Calculus of Variations Problems with State Constraints

abstract: We consider non-autonomous calculus of variations problems with a state constraint represented by a given closed set. We prove that if the interior of the Clarke tangent cone of the state constraint set is non-empty (this is the constraint qualification that we suggest here), then the necessary optimality conditions apply in the normal form. We establish normality results for (weak) local minimizers and global minimizers, employing two different approaches and invoking slightly diverse assumptions. More precisely, for the local minimizers result, the Lagrangian is supposed to be Lipschitz with respect to the state variable, and just lower semicontinuous in its third variable. On the other hand, the approach for the global minimizers result (which is simpler) requires the Lagrangian to be convex with respect to its third variable, but the Lipschitz constant of the Lagrangian with respect to the state variable might now depend on time.


title: An indirect numerical method for a time-optimal state-constrained control problem in a steady two-dimensional fluid flow

abstract: This article concerns the problem of computing solutions to state-constrained optimal control problems whose trajectory is affected by a flow field. This general mathematical framework is particularly pertinent to the requirements underlying the control of Autonomous Underwater Vehicles in realistic scenarii. The key contribution consists in devising a computational indirect method which becomes effective in the numerical computation of extremals to optimal control problems with state constraints by using the maximum principle in Gamkrelidze's form in which the measure Lagrange multiplier is ensured to be continuous. The specific problem of time-optimal control of an Autonomous Underwater Vehicle in a bounded space set, subject to the effect of a flow field and with bounded actuation, is used to illustrate the proposed approach. The corresponding numerical results are presented and discussed.


title: Engineering $\mathrm{SU}(1,1) \otimes \mathrm{SU}(1,1)$ vibrational states

abstract: We propose an ideal scheme for preparing vibrational $\mathrm{SU(1,1)} \otimes \mathrm{SU(1,1)}$ states in a two-dimensional ion trap using red and blue second sideband resolved driving of two orthogonal vibrational modes. Symmetric and asymmetric driving provide two regimes to realize quantum state engineering of the vibrational modes. In one regime, we show that time evolution synthesizes so-called $\mathrm{SU}(1,1)$ Perelomov coherent states, that is separable squeezed states and their superposition too. The other regime allows engineering of lossless 50/50 $\mathrm{SU}(2)$ beam splitter states that are entangled states. These ideal dynamics are reversible, thus, the non-classical and entangled states produced by our schemes might be used as resources for interferometry.


title: Coming: a Tool for Mining Change Pattern Instances from Git Commits

abstract: Software repositories such as Git have become a relevant source of information for software engineer researcher. For instance, the detection of Commits that fulfill a given criterion (e.g., bugfixing commits) is one of the most frequent tasks done to understand the software evolution. However, to our knowledge, there is not open-source tools that, given a Git repository, returns all the instances of a given change pattern. In this paper we present Coming, a tool that takes an input a Git repository and mines instances of change patterns on each commit. For that, Coming computes fine-grained changes between two consecutive revisions, analyzes those changes to detect if they correspond to an instance of a change pattern (specified by the user using XML), and finally, after analyzing all the commits, it presents a) the frequency of code changes and b) the instances found on each commit. We evaluate Coming on a set of 28 pairs of revisions from Defects4J, finding instances of change patterns that involve If conditions on 26 of them.


title: Teaching the concept of convolution and correlation using Fourier transform

abstract: Convolution operation is indispensable in studying analog optical and digital signal processing. Equally important is the correlation operation. The time domain community often teaches convolution and correlation only with one dimensional time signals. That does not clearly demonstrate the effect of convolution and correlation between two signals. Instead if we consider two dimensional spatial signals, the convolution and correlation operations can be very clearly explained. In this paper, we propose a lecture demonstration of convolution and correlation between two spatial signals using the Fourier transform tool. Both simulation and optical experiments are possible using a variety of object transparencies. The demonstration experiments help to clearly explain the similarity and the difference between convolution and correlation operations. This method of teaching using simulation and hands-on experiments can stimulate the curiosity of the students. The feedback of the students, in my class teaching, has been quite encouraging.


title: MsCGAN: Multi-scale Conditional Generative Adversarial Networks for Person Image Generation

abstract: To synthesize high-quality person images with arbitrary poses is challenging. In this paper, we propose a novel Multi-scale Conditional Generative Adversarial Networks (MsCGAN), aiming to convert the input conditional person image to a synthetic image of any given target pose, whose appearance and the texture are consistent with the input image. MsCGAN is a multi-scale adversarial network consisting of two generators and two discriminators. One generator transforms the conditional person image into a coarse image of the target pose globally, and the other is to enhance the detailed quality of the synthetic person image through a local reinforcement network. The outputs of the two generators are then merged into a synthetic, discriminant and high-resolution image. On the other hand, the synthetic image is downsampled to multiple resolutions as the input to multi-scale discriminator networks. The proposed multi-scale generators and discriminators handling different levels of visual features can benefit to synthesizing high-resolution person images with realistic appearance and texture. Experiments are conducted on the Market-1501 and DeepFashion datasets to evaluate the proposed model, and both qualitative and quantitative results demonstrate the superior performance of the proposed MsCGAN.


title: On Uniform Approximations of Normal Distributions By Jacobi Theta Functions

abstract: In this short note we study uniform approximations to the normal distributions by Jacobi theta functions. We shall show that scaled theta functions approach to a normal distribution exponentially fast.


title: Traces and inverse nodal problems for a class of Sturm-Liouville operators with retarded argument

abstract: In this study, we investigate the traces and solutions of inverse nodal problems of discontinuous Sturm-Liouville operators with retarded argument and with a finite number of transmission conditions.


title: Bayesian Distance Clustering

abstract: Model-based clustering is widely-used in a variety of application areas. However, fundamental concerns remain about robustness. In particular, results can be sensitive to the choice of kernel representing the within-cluster data density. Leveraging on properties of pairwise differences between data points, we propose a class of Bayesian distance clustering methods, which rely on modeling the likelihood of the pairwise distances in place of the original data. Although some information in the data is discarded, we gain substantial robustness to modeling assumptions. The proposed approach represents an appealing middle ground between distance- and model-based clustering, drawing advantages from each of these canonical approaches. We illustrate dramatic gains in the ability to infer clusters that are not well represented by the usual choices of kernel. A simulation study is included to assess performance relative to competitors, and we apply the approach to clustering of brain genome expression data. Keywords: Distance-based clustering; Mixture model; Model-based clustering; Model misspecification; Pairwise distance matrix; Partial likelihood; Robustness.


title: A new characterization of the discrete Sugeno integral

abstract: We introduce a new property of the discrete Sugeno integrals which can be seen as their characterization, too. This property, compatibility with respect to congruences on $[0,1]$, stresses the importance of the Sugeno integrals in multicriteria decision support as well.


title: Some tastings in Morales-Ramis Theory

abstract: In this paper we present a short material concerning to some results in Morales-Ramis theory, which relates two different notions of integrability: Integrability of Hamiltonian Systems through Liouville Arnold Theorem and Integrability of Linear Differential Equations through Differential Galois Theory. As contribution, we obtain the abelian differential Galois group of the variational equation related to a bi-parametric Hamiltonian system,


title: Fairness for Whom? Critically reframing fairness with Nash Welfare Product

abstract: Recent studies on disparate impact in machine learning applications have sparked a debate around the concept of fairness along with attempts to formalize its different criteria. Many of these approaches focus on reducing prediction errors while maximizing sole utility of the institution. This work seeks to reconceptualize and critically frame the existing discourse on fairness by underlining the implicit biases embedded in common understandings of fairness in the literature and how they contrast with its corresponding economic and legal definitions. This paper expands the concept of utility and fairness by bringing in concepts from established literature in welfare economics and game theory. We then translate these concepts for the algorithmic prediction domain by defining a formalization of Nash Welfare Product that seeks to expand utility by collapsing that of the institution using the prediction tool and the individual subject to the prediction into one function. We then apply a modulating function that makes the fairness and welfare trade-offs explicit based on designated policy goals and then apply it to a temporal model to take into account the effects of decisions beyond the scope of one-shot predictions. We apply this on a binary classification problem and present results of a multi-epoch simulation based on the UCI Adult Income dataset and a test case analysis of the ProPublica recidivism dataset that show that expanding the concept of utility results in a fairer distribution correcting for the embedded biases in the dataset without sacrificing the classifier accuracy.


title: Optimal Control of Dynamic Bipartite Matching Models

abstract: A dynamic bipartite matching model is given by a bipartite matching graph which determines the possible matchings between the various types of supply and demand items. Both supply and demand items arrive to the system according to a stochastic process. Matched pairs leave the system and the others wait in the queues, which induces a holding cost. We model this problem as a Markov Decision Process and study the discounted cost and the average cost problem. We fully characterize the optimal matching policy for complete matching graphs and for the N -shaped matching graph. In the former case, the optimal policy consists of matching everything and, in the latter case, it prioritizes the matchings in the extreme edges and is of threshold type for the diagonal edge. In addition, for the average cost problem, we compute the optimal threshold value. For more general graphs, we need to consider some assumptions on the cost of the nodes. For complete graphs minus one edge, we provide conditions on the cost of the nodes such that the optimal policy of the N-shaped matching graph extends to this case. For acyclic graphs, we show that, when the cost of the extreme edges is large, the optimal matching policy prioritizes the matchings in the extreme edges. We also study the W-shaped matching graph and, using simulations, we show that there are cases where it is not optimal to prioritize to matchings in the extreme edges.


title: Atomic photoionization cross-sections beyond the electric dipole approximation

abstract: A methodology is developed to compute photoionization cross-sections beyond the electric dipole (BED) approximation from response theory, using Gaussian type orbitals and plane waves for the initial and final states, respectively. The methodology is applied to compute photoionization cross-sections of atoms and ions from the first four rows of the Periodic Table. Analysing the error due to the plane wave description of the photoelectron, we find kinetic energy and concomitant photon energy thresholds above which the plane wave approximation becomes applicable. The correction introduced by going beyond the electric dipole approximation increases with photon energy and depends on the spatial extension of the initial state. In general, the corrections are below 10% for most elements, at a photon energy reaching up to 12 keV.


title: Strong Constraints on Fuzzy Dark Matter from Ultrafaint Dwarf Galaxy Eridanus II

abstract: The fuzzy dark matter (FDM) model treats DM as a bosonic field with astrophysically large de Broglie wavelength. A striking feature of this model is $\mathcal{O}(1)$ fluctuations in the dark matter density on time scales which are shorter than the gravitational timescale. Including for the first time the effect of core oscillations, we demonstrate how such fluctuations lead to heating of star clusters, and thus an increase in their size over time. From the survival of the old star cluster in Eridanus II we infer $m_a\gtrsim 0.6\rightarrow 1\times 10^{-19}\text{ eV}$ within modelling uncertainty if FDM is to compose all of the DM, and derive constraints on the FDM fraction at lower masses. The subhalo mass function in the Milky Way implies $m_a\gtrsim 0.8\times 10^{-21}\text{ eV}$ to successfully form Eridanus II. The window between $10^{-21}\text{ eV}\lesssim m_a\lesssim 10^{-20}\text{ eV}$ is affected by narrow band resonances, and the limited applicability of the diffusion approximation. Some of this window may be consistent with observations of Eridanus II and more detailed investigations are required.


title: New and Simplified Distributed Algorithms for Weighted All Pairs Shortest Paths

abstract: We consider the problem of computing all pairs shortest paths (APSP) and shortest paths for k sources in a weighted graph in the distributed CONGEST model. For graphs with non-negative integer edge weights (including zero weights) we build on a recent pipelined algorithm to obtain $\tilde{O}(\lambda^{1/4}\cdot n^{5/4})$ in graphs with non-negative integer edge-weight at most $\lambda$, and $\tilde{O}(n \cdot \bigtriangleup^{1/3})$ rounds for shortest path distances at most $\bigtriangleup$. Additionally, we simplify some of the procedures in the earlier APSP algorithms for non-negative edge weights in [HNS17,ARKP18]. We also present results for computing h-hop shortest paths and shortest paths from $k$ given sources. In other results, we present a randomized exact APSP algorithm for graphs with arbitrary edge weights that runs in $\tilde{O}(n^{4/3})$ rounds w.h.p. in n, which improves the previous best $\tilde{O}(n^{3/2})$ bound, which is deterministic. We also present an $\tilde{O}(n/\epsilon^2)$-round deterministic $(1+\epsilon)$ approximation algorithm for graphs with non-negative $poly(n)$ integer weights (including zero edge-weights), improving results in [Nanongkai14,LP15] that hold only for positive integer weights.


title: Congruences and the discrete Sugeno integrals on bounded distributive lattices

abstract: We study compatible aggregation functions on a general bounded distributive lattice $L$, where the compatibility is related to the congruences on $L$. As a by-product, a new proof of an earlier result of G. Gr\"atzer is obtained. Moreover, our results yield a new characterization of discrete Sugeno integrals on bounded distributive lattices.


title: Isogenies between K3 Surfaces over $\bar{\mathbb{F}}_p$

abstract: We generalize Mukai and Shafarevich's definitions of isogenies between K3 surfaces over $\mathbb{C}$ to an arbitrary perfect field and describe how to construct isogenous K3 surfaces over $\bar{\mathbb{F}}_p$ by prescribing linear algebraic data when $p$ is large. The main step is to show that isogenies between Kuga-Satake abelian varieties induce isogenies between K3 surfaces, in the context of integral models of Shimura varieties. As a byproduct, we show that every K3 surface of finite height admits a CM lifting under a mild assumption on $p$.


title: Means of infinite sets III

abstract: We study various topics, e.g. accumulation points by a mean, two types of derivative by a mean, two new continuity and a boundedness concepts, we construct new means from old ones, finally we investigate the limit of means.


title: The Rotational Evolution of Young, Binary M Dwarfs

abstract: We have analysed K2 light curves for more than 3,000 low mass stars in the $\sim$8 Myr old Upper Sco association, the $\sim$125 Myr age Pleiades open cluster and the $\sim$700 Myr old Hyades and Praesepe open clusters to determine stellar rotation rates. Many of these K2 targets show two distinct periods, and for the lowest mass stars in these clusters virtually all of these systems with two periods are photometric binaries. The most likely explanation is that we are detecting the rotation periods for both components of these binaries. We explore the evolution of the rotation rate in both components of photometric binaries relative to one another and to non-photometric binary stars. In Upper Sco and the Pleiades, these low mass binary stars have periods that are much shorter on average and much closer to each other than would be true if drawn at random from the M dwarf single stars. In Upper Sco, this difference correlates strongly with the presence or absence of infrared excesses due to primordial circumstellar disks -- the single star population includes many stars with disks, and their rotation periods are distinctively longer on average than their binary star cousins of the same mass. By Praesepe age, the significance of the difference in rotation rate between the single and binary low mass dMs is much less, suggesting that angular momentum loss from winds for fully-convective zero-age main sequence stars erases memory of the rotation rate dichotomy for binary and single very low mass stars at later ages.


title: Optimal cell transport in straight channels and networks

abstract: Flux of rigid or soft particles (such as drops, vesicles, red blood cells, etc.) in a channel is a complex function of particle concentration, which depends on the details of induced dissipation and suspension structure due to hydrodynamic interactions with walls or between neighboring particles. Through two-dimensional and three-dimensional simulations and a simple model that reveals the contribution of the main characteristics of the flowing suspension, we discuss the existence of an optimal volume fraction for cell transport and its dependence on the cell mechanical properties. The example of blood is explored in detail, by adopting the commonly used modeling of red blood cells dynamics. We highlight the complexity of optimization at the level of a network, due to the antagonist evolution of local volume fraction and optimal volume fraction with the channels diameter. In the case of the blood network, the most recent results on the size evolution of vessels along the circulatory network of healthy organs suggest that the red blood cell volume fraction (hematocrit) of healthy subjects is close to optimality, as far as transport only is concerned. However, the hematocrit value of patients suffering from diverse red blood cel pathologies may strongly deviate from optimality.

categories: physics.flu-dyn cond-mat.soft physics.bio-ph q-bio.CB


title: Aspects of asymptotic safety for quantum gravity

abstract: We study fixed points of quantum gravity with renormalisation group methods, and a procedure to remove convergence-limiting poles from the flow. The setup is tested within the $f(R)$ approximation for gravity by solving exact recursive relations up to order $R^{70}$ in the Ricci scalar, combined with resummations and numerical integration. Results include fixed points, scaling exponents, gap in the eigenvalue spectrum, dimensionality of the UV critical surface, fingerprints for weak coupling, and quantum equations of motion. Our findings strengthen the view that ``most of quantum gravity'' is rather weakly coupled. Another novelty are a pair of de Sitter solutions for quantum cosmology, whose occurrence is traced back to the removal of poles. We also address slight disparities of results in the literature, and give bounds on the number of fundamentally free parameters of quantum gravity.


title: Nonreciprocal localization of ultrasound in a viscous medium with asymmetric scatterers

abstract: A two-dimensional phononic crystal with asymmetric scatterers is used for the study of Anderson localization of sound along one-dimensional disorder produced by random orientation of metallic rods. An exponentially weak transmission of ultrasound is demonstrated for the waves propagating along the direction of disorder. In the perpendicular direction where the scatterers are ordered, sound propagates as extended wave. The {\it PT}-symmetry of the system is broken by dissipative viscous losses and asymmetric shape of the scatterers. Nonreciprocal transmission of sound is observed for both, ordered and disordered, directions. In the localized regime, the nonreciprocity is manifested through different values of localization length for sound propagating in the opposite directions.


title: Nonlinear integro-differential operator regression with neural networks

abstract: This note introduces a regression technique for finding a class of nonlinear integro-differential operators from data. The method parametrizes the spatial operator with neural networks and Fourier transforms such that it can fit a class of nonlinear operators without needing a library of a priori selected operators. We verify that this method can recover the spatial operators in the fractional heat equation and the Kuramoto-Sivashinsky equation from numerical solutions of the equations.

categories: cs.LG physics.comp-ph physics.data-an stat.ML


title: Federated Learning in Distributed Medical Databases: Meta-Analysis of Large-Scale Subcortical Brain Data

abstract: At this moment, databanks worldwide contain brain images of previously unimaginable numbers. Combined with developments in data science, these massive data provide the potential to better understand the genetic underpinnings of brain diseases. However, different datasets, which are stored at different institutions, cannot always be shared directly due to privacy and legal concerns, thus limiting the full exploitation of big data in the study of brain disorders. Here we propose a federated learning framework for securely accessing and meta-analyzing any biomedical data without sharing individual information. We illustrate our framework by investigating brain structural relationships across diseases and clinical cohorts. The framework is first tested on synthetic data and then applied to multi-centric, multi-database studies including ADNI, PPMI, MIRIAD and UK Biobank, showing the potential of the approach for further applications in distributed analysis of multi-centric cohorts


title: HeartBEAT: Heart Beat Estimation through Adaptive Tracking

abstract: In this paper, we propose an algorithm denoted as HeartBEAT that tracks heart rate from wrist-type photoplethysmography (PPG) signals and simultaneously recorded three-axis acceleration data. HeartBEAT contains three major parts: spectrum estimation of PPG signals and acceleration data, elimination of motion artifacts in PPG signals using recursive least Square (RLS) adaptive filters, and auxiliary heuristics. We tested HeartBEAT on the 22 datasets provided in the 2015 IEEE Signal Processing Cup. The first ten datasets were recorded from subjects performing forearm and upper-arm exercises, jumping, or pushing-up. The last twelve datasets were recorded from subjects running on tread mills. The experimental results were compared to the ground truth heart rate, which comes from simultaneously recorded electrocardiogram (ECG) signals. Compared to state-of-the-art algorithms, HeartBEAT not only produces comparable Pearson's correlation and mean absolute error, but also higher Spearman's rho and Kendall's tau.


title: Drastic slowdown of shear waves in unjammed granular suspensions

abstract: We present an experimental investigation of shear elastic wave propagation along the surface of a dense granular suspension. Using an ultrafast ultrasound scanner, we monitor the softening of the shear wave velocity inside the optically opaque medium as the driving amplitude increases. For such nonlinear behavior two regimes are found: in the first regime, we observe a significant shear modulus weakening, but without visible grain rearrangements. In the second regime, there is a clear grain rearrangement accompanied by a modulus decrease up to 88%. A friction model is proposed to describe the interplay between nonlinear elasticity and plasticity, which highlights the crucial effect of contact slipping before contact breaking. Investigation of these nonlinear shear waves may bridge the gap between two disjoint approaches for describing the dynamics near unjamming: linear elastic soft modes and nonlinear collisional shock.


title: Global minima for optimal control of the obstacle problem

abstract: An optimal control problem subject to an elliptic obstacle problem is studied. We obtain a numerical approximation of this problem by discretising the PDE obtained via a Moreau--Yosida type penalisation. For the resulting discrete control problem we provide a condition that allows to decide whether a solution of the necessary first order conditions is a global minimum. In addition we show that the corresponding result can be transferred to the limit problem provided that the above condition holds uniformly in the penalisation and discretisation parameters. Numerical examples with unique global solutions are presented.


title: Statistical Treatment of Inverse Problems Constrained by Differential Equations-Based Models with Stochastic Terms

abstract: This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation (PDE) with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations (DAEs). In both cases we show that the variogram and the hybrid scores show better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.


title: Torque on slowly moving electric or magnetic dipoles in vacuo

abstract: The torque on a moving electric or magnetic dipole in slow motion is deduced using the Lorentz transformation of the fields to first order in v/c. It is shown that the obtained equations are independent of the model adopted for the dipole, whether it is of Amperian or Gilbertian type, thus showing the complete validity of the Amp\`ere equivalence principle even in dynamical conditions. The torque is made of three terms: beside that due to the direct torque on the dipole there are two more terms: one due to the torque on the associated perpendicular dual-dipole caused by motion, while the other is the inertial torque due to the displacement of the dipole which carries with it the field linear momentum, or the hidden momentum.


title: EdgeSpeechNets: Highly Efficient Deep Neural Networks for Speech Recognition on the Edge

abstract: Despite showing state-of-the-art performance, deep learning for speech recognition remains challenging to deploy in on-device edge scenarios such as mobile and other consumer devices. Recently, there have been greater efforts in the design of small, low-footprint deep neural networks (DNNs) that are more appropriate for edge devices, with much of the focus on design principles for hand-crafting efficient network architectures. In this study, we explore a human-machine collaborative design strategy for building low-footprint DNN architectures for speech recognition through a marriage of human-driven principled network design prototyping and machine-driven design exploration. The efficacy of this design strategy is demonstrated through the design of a family of highly-efficient DNNs (nicknamed EdgeSpeechNets) for limited-vocabulary speech recognition. Experimental results using the Google Speech Commands dataset for limited-vocabulary speech recognition showed that EdgeSpeechNets have higher accuracies than state-of-the-art DNNs (with the best EdgeSpeechNet achieving ~97% accuracy), while achieving significantly smaller network sizes (as much as 7.8x smaller) and lower computational cost (as much as 36x fewer multiply-add operations, 10x lower prediction latency, and 16x smaller memory footprint on a Motorola Moto E phone), making them very well-suited for on-device edge voice interface applications.

categories: eess.AS cs.LG cs.NE cs.SD eess.SP stat.ML


title: $2\times2$ Hypergeometric operators with diagonal eigenvalues

abstract: In this work we classify all the order-two Hypergeometric operators $D$, symmetric with respect to some $2\times 2$ irreducible matrix-weight $W$ such that $DP_n=P_n\left(\begin{smallmatrix} \lambda_n&0\\0&\mu_n \end{smallmatrix} \right)$ with no repetition among the eigenvalues $\{\lambda_n,\mu_n\}_{n\in\mathbb N_0}$, where $\{P_n\}_{n\in\mathbb N_0}$ is the (unique) sequence of monic orthogonal polynomials with respect to $W$. We obtain, in a very explicit way, a three parameter family of such operators and weights. We also give the corresponding monic orthongonal polynomials, their three term recurrence relation and their squared matrix-norms.


title: Dissipative Dynamics of a Single Polymer in Solution: A Lowe-Andersen Approach

abstract: We study the equilibrium dynamics of a single polymer chain under good solvent condition. Special emphasis is laid on varying the drag force experienced by the chain while it moves. To this end we model the solvent in a mesoscopic manner by employing the Lowe-Andersen approach of dissipative particle dynamics which is known to reproduce hydrodynamic effects. Our approach captures the correct static behavior in equilibrium. Regarding the dynamics, we investigate the scaling of the self-diffusion coefficient $D$ with respect to the length of the polymer $N$, yielding results that are compatible with the Zimm scaling $D \sim N^{-3/5}$.


title: Long time behavior of a mean-field model of interacting neurons

abstract: We study the long time behavior of the solution to some McKean-Vlasov stochastic differential equation (SDE) driven by a Poisson process. In neuroscience, this SDE models the asymptotic dynamic of the membrane potential of a spiking neuron in a large network. We prove that for a small enough interaction parameter, any solution converges to the unique (in this case) invariant measure. To this aim, we first obtain global bounds on the jump rate and derive a Volterra type integral equation satisfied by this rate. We then replace temporary the interaction part of the equation by a deterministic external quantity (we call it the external current). For constant current, we obtain the convergence to the invariant measure. Using a perturbation method, we extend this result to more general external currents. Finally, we prove the result for the non-linear McKean-Vlasov equation.


title: Ferromagnetic nodal-line metal in monolayer {\em h}-InC

abstract: Based on first-principles calculations, we predict a new two-dimensional ferromagnetic material that exhibits exotic Fermi surface topology. We show that monolayer hexagonal indium carbide ({\em h}-InC) is thermodynamically and dynamically stable, and it energetically favors the ferromagnetic ordering of spins. The perfectly planar geometry in two dimensions, together with ferromagnetism, gives rise to a unique opportunity to encounter intriguing electronic properties, captured in the Fermi surface and band topology. We show that multiple nodal lines coexist in momentum space, accompanied by the electron and hole pockets that touch each other linearly at the nodal lines. Inclusion of spin-orbit coupling enriches the magnetic and electronic properties of {\em h}-InC. Spin-orbit coupling leads to an easy-plane type magnetocrystalline anisotropy, and the nodal lines can be tuned into topological nodal points, contingent upon the magnetization direction. Symmetry analysis and a tight-binding model are provided to explain the nodal structure of the bands. Our findings suggest {\em h}-InC as a new venue for supporting carbon-based magnetism and exotic band topology in two dimensions.


title: Nonparametric Bayesian Lomax delegate racing for survival analysis with competing risks

abstract: We propose Lomax delegate racing (LDR) to explicitly model the mechanism of survival under competing risks and to interpret how the covariates accelerate or decelerate the time to event. LDR explains non-monotonic covariate effects by racing a potentially infinite number of sub-risks, and consequently relaxes the ubiquitous proportional-hazards assumption which may be too restrictive. Moreover, LDR is naturally able to model not only censoring, but also missing event times or event types. For inference, we develop a Gibbs sampler under data augmentation for moderately sized data, along with a stochastic gradient descent maximum a posteriori inference algorithm for big data applications. Illustrative experiments are provided on both synthetic and real datasets, and comparison with various benchmark algorithms for survival analysis with competing risks demonstrates distinguished performance of LDR.


title: Deep Person Re-identification for Probabilistic Data Association in Multiple Pedestrian Tracking

abstract: We present a data association method for vision-based multiple pedestrian tracking, using deep convolutional features to distinguish between different people based on their appearances. These re-identification (re-ID) features are learned such that they are invariant to transformations such as rotation, translation, and changes in the background, allowing consistent identification of a pedestrian moving through a scene. We incorporate re-ID features into a general data association likelihood model for multiple person tracking, experimentally validate this model by using it to perform tracking in two evaluation video sequences, and examine the performance improvements gained as compared to several baseline approaches. Our results demonstrate that using deep person re-ID for data association greatly improves tracking robustness to challenges such as occlusions and path crossings.


title: Differential Galois Theory and Isomonodromic Deformations

abstract: We present a geometric setting for the differential Galois theory of $G$-invariant connections with parameters. As an application of some classical results on differential algebraic groups and Lie algebra bundles, we see that the Galois group of a connection with parameters with simple structural group $G$ is determined by its isomonodromic deformations. This allows us to compute the Galois groups with parameters of the general Fuchsian special linear system and of Gauss hypergeometric equation.


title: Weak Semi-Markov CRFs for NP Chunking in Informal Text

abstract: This paper introduces a new annotated corpus based on an existing informal text corpus: the NUS SMS Corpus (Chen and Kan, 2013). The new corpus includes 76,490 noun phrases from 26,500 SMS messages, annotated by university students. We then explored several graphical models, including a novel variant of the semi-Markov conditional random fields (semi-CRF) for the task of noun phrase chunking. We demonstrated through empirical evaluations on the new dataset that the new variant yielded similar accuracy but ran in significantly lower running time compared to the conventional semi-CRF.


title: Ultrasonic tracking of a sinking ball in a vibrated dense granular suspension

abstract: Observing and understanding the movement of an intruder through opaque dense suspensions such as quicksand remains a practical and conceptual challenge. Here we use an ultrasonic probe to investigate the dynamics of a steel ball sinking in a 3D dense glass bead packing saturated by water. We show that the frictional model developed for dry granular media can be used to describe the ball motion induced by horizontal vibration. From this rheology we infer the static friction coefficient and effective viscosity that decrease when increasing the vibration intensity. Our main finding is that the vibration-induced reduction of the yield stress and increase of the sinking depth are presumably due to induced slipping at the grain contacts but without visible plastic rearrangements of grains, in contrast to dry granular packings. To explain these results, we propose a mechanism of acoustic lubrication that reduces the inter-particle friction and leads to a decrease of the yield stress. This scenario is different from the mechanism of liquefaction usually invoked in loosely packed quicksands where the vibration-induced compaction increases the pore pressure and decreases the confining pressure on the solid skeleton, thus reducing the granular resistance to external load.


title: On the possibility of homogeneous nucleation of water droplets and/or ice crystals during the bounces of a cavitation bubble

abstract: Acoustic cavitation is known to trigger ice nucleation in supercooled water. Several competing and still debatable mechanisms have been proposed in the literature and are related to the pressure field in the vicinity of the bubble at the end of its collapse. Numerical simulations of the bubble dynamics show that during the bubble expansions in the bounces following the main collapse, the bubble core temperature reaches values far below 0 {\deg}C for time periods of about 500 ns. The water vapour present in the bubble during these time intervals explores the liquid and the solid region of the phase diagrams before going back to the vapour region. On the base of approximate nucleation kinetics calculations, we examine to what extent liquid droplets could nucleate homogenously in the bubble core during these excursions. We also discuss the possibility that the nucleated clusters reach the bubble wall and trigger ice nucleation in the surrounding liquid if the latter is supercooled.


title: Deformations of algebraic schemes via Reedy-Palamodov cofibrant resolutions

abstract: Let $X$ be a Noetherian separated and finite dimensional scheme over a field $\mathbb{K}$ of characteristic zero. The goal of this paper is to study deformations of $X$ over a differential graded local Artin $\mathbb{K}$-algebra by using local Tate-Quillen resolutions, i.e., the algebraic analog of the Palamodov's resolvent of a complex space. The above goal is achieved by describing the DG-Lie algebra controlling deformation theory of a diagram of differential graded commutative algebras, indexed by a direct Reedy category.


title: Leveraging Elastic instabilities for Amplified Performance: spine-inspired high-speed and high-force soft robots

abstract: Soft machines typically exhibit slow locomotion speed and low manipulation strength because of intrinsic limitations of soft materials. Here, we present a generic design principle that harnesses mechanical instability for a variety of spine-inspired fast and strong soft machines. Unlike most current soft robots that are designed as inherently and unimodally stable, our design leverages tunable snap-through bistability to fully explore the ability of soft robots to rapidly store and release energy within tens of milliseconds. We demonstrate this generic design principle with three high-performance soft machines: High-speed cheetah-like galloping crawlers with locomotion speeds of 2.68 body length/s, high-speed underwater swimmers (0.78 body length/s), and tunable low-to-high-force soft grippers with over 1 to 103 stiffness modulation (maximum load capacity is 11.4 kg). Our study establishes a new generic design paradigm of next-generation high-performance soft robots that are applicable for multifunctionality, different actuation methods, and materials at multiscales.


title: Finite Volume Simulation Framework for Die Casting with Uncertainty Quantification

abstract: The present paper describes the development of a novel and comprehensive computational framework to simulate solidification problems in materials processing, specifically casting processes. Heat transfer, solidification and fluid flow due to natural convection are modeled. Empirical relations are used to estimate the microstructure parameters and mechanical properties. The fractional step algorithm is modified to deal with the numerical aspects of solidification by suitably altering the coefficients in the discretized equation to simulate selectively only in the liquid and mushy zones. This brings significant computational speed up as the simulation proceeds. Complex domains are represented by unstructured hexahedral elements. The algebraic multigrid method, blended with a Krylov subspace solver is used to accelerate convergence. State of the art uncertainty quantification technique is included in the framework to incorporate the effects of stochastic variations in the input parameters. Rigorous validation is presented using published experimental results of a solidification problem.


title: Security-Reliability Tradeoff for Distributed Antenna Systems in Heterogeneous Cellular Networks

abstract: In this paper, we investigate physical-layer security for a spectrum-sharing heterogeneous cellular network comprised of a macro cell and a small cell, where a passive eavesdropper is assumed to tap the transmissions of both the macro cell and small cell. In the macro cell, a macro base station (MBS) equipped with multiple distributed antennas sends its confidential information to a macro user (MU) through an opportunistic transmit antenna. Meanwhile, in the small cell, a small base station (SBS) transmits to a small user (SU) over the same spectrum used by MBS. We propose an interference-canceled opportunistic antenna selection (IC-OAS) scheme to enhance physical-layer security for the heterogeneous network. To be specific, when MBS sends its confidential message to MU through an opportunistic distributed antenna, a special signal is artificially designed and emitted at MBS to ensure that the received interference at MU from SBS is canceled out. For comparison, the conventional interference-limited opportunistic antenna selection (IL-OAS) is considered as a benchmark. We characterize the security-reliability tradeoff (SRT) for the proposed IC-OAS and conventional IL-OAS schemes in terms of deriving their closed-form expressions of intercept probability and outage probability. Numerical results show that compared with the conventional IL-OAS, the proposed IC-OAS scheme not only brings SRT benefits to the macro cell, but also has the potential of improving the SRT of small cell by increasing the number of distributed antennas. Additionally, by jointly taking into account the macro cell and small cell, an overall SRT of the proposed IC-OAS scheme is shown to be significantly better than that of the conventional IL-OAS approach in terms of a sum intercept probability versus sum outage probability.


title: Monatomic Co, CoO$_2$, and CoO$_3$ Nanowires on Ir(100) and Pt(100) surfaces: Formation, Structure, and Energetics

abstract: In this study we investigate the structural and chemical changes of monatomic CoO$_2$ chains grown self-organized on the Ir(100) surface [P. Ferstl et al., PRL 117, 2016, 046101] and on Pt(100) under reducing and oxidizing conditions. By a combination of quantitative low-energy electron diffraction, scanning tunnelling microscopy, and density functional theory we show that the cobalt oxide wires are completely reduced by H$_2$ at temperatures above 320 K and a 3x1 ordered Ir$_2$Co or Pt$_2$Co surface alloy is formed. Depending on temperature the surface alloy on Ir(100) is either hydrogen covered (T < 400 K) or clean and eventually undergoes an irreversible order-disorder transition at about 570 K. The Pt$_2$Co surface alloy disorders with the desorption of hydrogen, whereby Co submerges into subsurface sites. Vice versa, applying stronger oxidants than O$_2$ such as NO$_2$ leads to the formation of CoO3 chains on Ir(100) in a 3x1 superstructure. On Pt(100) such a CoO$_3$ phase could not be prepared so far, which however, is due to the UHV conditions of our experiments. As revealed by theory this phase will become stable in a regime of higher pressure. In general, the structures can be reversibly switched on both surfaces using the respective agents O$_2$, NO$_2$ and H$_2$.


title: Supervising strong learners by amplifying weak experts

abstract: Many real world learning tasks involve complex or hard-to-specify objectives, and using an easier-to-specify proxy can lead to poor performance or misaligned behavior. One solution is to have humans provide a training signal by demonstrating or judging performance, but this approach fails if the task is too complicated for a human to directly evaluate. We propose Iterated Amplification, an alternative training strategy which progressively builds up a training signal for difficult problems by combining solutions to easier subproblems. Iterated Amplification is closely related to Expert Iteration (Anthony et al., 2017; Silver et al., 2017), except that it uses no external reward function. We present results in algorithmic environments, showing that Iterated Amplification can efficiently learn complex behaviors.


title: Probing the IGMF with the next generation Cherenkov telescopes

abstract: Intergalactic space is believed to contain non-zero magnetic fields (the Intergalactic Magnetic Field: IGMF) which at scales of Mpc would have intensities below $10^{-9}$ G. Very high energy (VHE $>$100 GeV) gamma rays coming from blazars can produce e$^+$e$^-$ pairs when interacting with the Extragalactic Background Light (EBL) and the Cosmic Microwave Background, generating an electromagnetic cascade of Mpc scale. The IGMF may produce a detectable broadening of the emission beam that could lead to important constrains both on the IGMF intensity and its coherence length. Using the Monte Carlo-based Elmag code, we simulate the electromagnetic cascade corresponding to two detected TeV sources: PKS 2155-304 visible from the South and H1426+428 visible from the North. Assuming an EBL model and intrinsic spectral properties of the sources we obtain the spectral and angular distribution of photons when they arrive at Earth. We include the response of the next generation Cherenkov telescopes by using simplified models for CTA (Cherenkov Telescope Array)-south and CTA-north based on a full simulation of each array performance. Combining the instrument properties with the simulated source fluxes, we calculate the telescope point spread function for null and non-null IGMF intensities and develop a method to test the statistical feasibility of detecting IGMF imprints by comparing the resulting angular distributions. Our results show that for the analysed source PKS 2155-304 corresponding to the southern site, CTA should be able to detect IGMF with intensities stronger than 10$^{-14.5}$G within an observation time of $\sim$100 hours.


title: Conceptual Organization is Revealed by Consumer Activity Patterns

abstract: Meaning may arise from an element's role or interactions within a larger system. For example, hitting nails is more central to people's concept of a hammer than its particular material composition or other intrinsic features. Likewise, the importance of a web page may result from its links with other pages rather than solely from its content. One example of meaning arising from extrinsic relationships are approaches that extract the meaning of word concepts from co-occurrence patterns in large, text corpora. The success of these methods suggest that human activity patterns may reveal conceptual organization. However, texts do not directly reflect human activity, but instead serve a communicative function and are usually highly curated or edited to suit an audience. Here, we apply methods devised for text to a data source that directly reflects thousands of individuals' activity patterns, namely supermarket purchases. Using product co-occurrence data from nearly 1.3m shopping baskets, we trained a topic model to learn 25 high-level concepts (or "topics"). These topics were found to be comprehensible and coherent by both retail experts and consumers. Topics ranged from specific (e.g., ingredients for a stir-fry) to general (e.g., cooking from scratch). Topics tended to be goal-directed and situational, consistent with the notion that human conceptual knowledge is tailored to support action. Individual differences in the topics sampled predicted basic demographic characteristics. These results suggest that human activity patterns reveal conceptual organization and may give rise to it.


title: Leveraging Product as an Activation Function in Deep Networks

abstract: Product unit neural networks (PUNNs) are powerful representational models with a strong theoretical basis, but have proven to be difficult to train with gradient-based optimizers. We present windowed product unit neural networks (WPUNNs), a simple method of leveraging product as a nonlinearity in a neural network. Windowing the product tames the complex gradient surface and enables WPUNNs to learn effectively, solving the problems faced by PUNNs. WPUNNs use product layers between traditional sum layers, capturing the representational power of product units and using the product itself as a nonlinearity. We find the result that this method works as well as traditional nonlinearities like ReLU on the MNIST dataset. We demonstrate that WPUNNs can also generalize gated units in recurrent neural networks, yielding results comparable to LSTM networks.


title: Learning to Recognize Discontiguous Entities

abstract: This paper focuses on the study of recognizing discontiguous entities. Motivated by a previous work, we propose to use a novel hypergraph representation to jointly encode discontiguous entities of unbounded length, which can overlap with one another. To compare with existing approaches, we first formally introduce the notion of model ambiguity, which defines the difficulty level of interpreting the outputs of a model, and then formally analyze the theoretical advantages of our model over previous existing approaches based on linear-chain CRFs. Our empirical results also show that our model is able to achieve significantly better results when evaluated on standard data with many discontiguous entities.


title: High-Density Qubit Wiring: Pin-Chip Bonding for Fully Vertical Interconnects

abstract: Large-scale quantum computers with more than $10^5$ qubits will likely be built within the next decade. Trapped ions, semiconductor devices, and superconducting qubits among other physical implementations are still confined in the realm of medium-scale quantum integration ($\sim 100$ qubits); however, they show promise toward large-scale quantum integration. Building large-scale quantum processing units will require truly scalable control and measurement classical coprocessors as well as suitable wiring methods. In this blue paper, we introduce a fully vertical interconnect that will make it possible to address $\sim 10^5$ superconducting qubits fabricated on a single silicon or sapphire chip: Pin-chip bonding. This method permits signal transmission from DC to $\sim 10$ GHz, both at room temperature and at cryogenic temperatures down to $\sim 10$ mK. At temperatures below $\sim 1$ K, the on-chip wiring contact resistance is close to zero and all signal lines are in the superconducting state. High-density wiring is achieved by means of a fully vertical interconnect that interfaces the qubit array with a network of rectangular coaxial ribbon cables. Pin-chip bonding is fully compatible with classical high-density test board applications as well as with other qubit implementations.


title: Forecasting Time Series with VARMA Recursions on Graphs

abstract: Graph-based techniques emerged as a choice to deal with the dimensionality issues in modeling multivariate time series. However, there is yet no complete understanding of how the underlying structure could be exploited to ease this task. This work provides contributions in this direction by considering the forecasting of a process evolving over a graph. We make use of the (approximate) time-vertex stationarity assumption, i.e., timevarying graph signals whose first and second order statistical moments are invariant over time and correlated to a known graph topology. The latter is combined with VAR and VARMA models to tackle the dimensionality issues present in predicting the temporal evolution of multivariate time series. We find out that by projecting the data to the graph spectral domain: (i) the multivariate model estimation reduces to that of fitting a number of uncorrelated univariate ARMA models and (ii) an optimal low-rank data representation can be exploited so as to further reduce the estimation costs. In the case that the multivariate process can be observed at a subset of nodes, the proposed models extend naturally to Kalman filtering on graphs allowing for optimal tracking. Numerical experiments with both synthetic and real data validate the proposed approach and highlight its benefits over state-of-the-art alternatives.


title: Strategy and guidelines for the calibration of the ITER radial neutron camera

abstract: A calibration procedure is proposed for the ITER Radial Neutron Camera. No in-vessel calibration using external neutron sources is required: instead, it is proposed to rely on embedded sources, reference ITER pulses and cross-calibration with ITER fission chambers and activation system coupled to Monte Carlo simulations of radiation transport for the validation of the RNC calibration and for its tracking during ITER lifetime.


title: Principal component analysis of sunspot cycle shape

abstract: We study the shape of sunspot cycles using the Wolf sunspot numbers and group sunspot numbers of solar cycles 1-23. We determine the most typical "model" cycles and the most asymmetric cycles, and test the validity of the two Waldmeier rules: the anti-correlation between cycle height and the length of its ascending phase (rule 1), and between cycle height and the length of the preceding cycle (rule 2). We applied the principal component analysis to sunspot cycles and studied the first two components, which describe the average cycle shape and cycle asymmetry, respectively. We also calculated their autocorrelation in order to study their recurrence properties. The best model cycles for Wolf numbers are SC12, SC14, and SC16, the successive even cycles from a long period of rather low overall solar activity. We find that the model cycles in eight different analyses using both sunspot series are almost exclusively even cycles. Correspondingly, the most asymmetric cycles are odd cycles. We find that both Waldmeier rules are valid for the whole Wolf number series of 23 cycles. Waldmeier rule 2 is also valid for group number series although its significance is weaker. Waldmeier rule 1 is not significant for the original group number series, but becomes significant for the proxy series. For separate centuries, Waldmeier rules are not always valid for Wolf numbers and very rarely for group numbers. Our results also offer a new interpretation for the Gnevyshev gap. In addition to being a local depression of solar activity, the Gnevyshev gap is a separatrix that divides cycles into two parts whose relative intensities determine the cycle asymmetry. The Gnevyshev gap is the zero value time of PC2, located approximately 33-42% into the cycle after its start.


title: Reverse Quantum Annealing Approach to Portfolio Optimization Problems

abstract: We investigate a hybrid quantum-classical solution method to the mean-variance portfolio optimization problems. Starting from real financial data statistics and following the principles of the Modern Portfolio Theory, we generate parametrized samples of portfolio optimization problems that can be related to quadratic binary optimization forms programmable in the analog D-Wave Quantum Annealer 2000Q. The instances are also solvable by an industry-established Genetic Algorithm approach, which we use as a classical benchmark. We investigate several options to run the quantum computation optimally, ultimately discovering that the best results in terms of expected time-to-solution as a function of number of variables for the hardest instances set are obtained by seeding the quantum annealer with a solution candidate found by a greedy local search and then performing a reverse annealing protocol. The optimized reverse annealing protocol is found to be more than 100 times faster than the corresponding forward quantum annealing on average.


title: Monotonic Distributive Semilattices

abstract: In the study of algebras related to non-classical logics, (distributive) semilattices are always present in the background. For example, the algebraic semantic of the $\{\rightarrow,\wedge,\top\}$-fragment of intuitionistic logic is the variety of implicative meet-semilattices \cite{CelaniImplicative} \cite{ChajdaHalasKuhr}. In this paper we introduce and study the class of distributive meet-semilattices endowed with a monotonic modal operator $m$. We study the representation theory of these algebras using the theory of canonical extensions and we give a topological duality for them. Also, we show how our new duality extends to some particular subclasses.


title: Constraints on quintessence scalar field models using cosmological observations

abstract: We consider a varieties of quintessence scalar field models in a homogeneous and isotropic geometry of the universe with zero spatial curvature aiming to provide stringent constraints using a series of cosmological data sets, namely, the cosmic microwave background observations (CMB), baryon acoustic oscillations (BAO), joint light curve analysis (JLA) from supernovae type Ia, redshift space distortions (RSD), and the cosmic chronometers (CC). From the qualitative evolution of the models, we find all of them are able to execute a fine transition from the past decelerating phase to the presently accelerating expansion where in addition, the equation of state of the scalar field (also the effective equation of state) might be close to that of the $\Lambda$CDM cosmology depending on its free parameters. From the observational analyses, we find that the scalar field parameters are unconstrained irrespective of all the observational datasets. In fact, we find that the quintessence scalar field models are pretty much determined by the CMB observations since any of the external datasets such as BAO, JLA, RSD, CC does not add any constraining power to CMB. Additionally, we observe a strong negative correlation between the parameters $H_0$ (present value of the Hubble parameter), $\Omega_{m0}$ (density parameter for the matter sector, i.e., cold dark matter plus baryons) exists, while no correlation between $H_0$, and $\sigma_8$ (amplitude of the matter fluctuation) are not correlated. We also comment that the present models are unable to reconcile the tension on $H_0$. Finally, we conclude our work with the Bayesian analyses which report that the non-interacting $\Lambda$CDM model is preferred over all the quintessence scalar field models.


title: Memory formation in matter

abstract: Memory formation in matter is a theme of broad intellectual relevance; it sits at the interdisciplinary crossroads of physics, biology, chemistry, and computer science. Memory connotes the ability to encode, access, and erase signatures of past history in the state of a system. Once the system has completely relaxed to thermal equilibrium, it is no longer able to recall aspects of its evolution. Memory of initial conditions or previous training protocols will be lost. Thus many forms of memory are intrinsically tied to far-from-equilibrium behavior and to transient response to a perturbation. This general behavior arises in diverse contexts in condensed matter physics and materials: phase change memory, shape memory, echoes, memory effects in glasses, return-point memory in disordered magnets, as well as related contexts in computer science. Yet, as opposed to the situation in biology, there is currently no common categorization and description of the memory behavior that appears to be prevalent throughout condensed-matter systems. Here we focus on material memories. We will describe the basic phenomenology of a few of the known behaviors that can be understood as constituting a memory. We hope that this will be a guide towards developing the unifying conceptual underpinnings for a broad understanding of memory effects that appear in materials.


title: Remote sensing to reduce the effects of spatial autocorrelation on design-based inference for forest inventory using systematic samples

abstract: Systematic sampling is often used to select plot locations for forest inventory estimation. However, it is not possible to derive a design-unbiased variance estimator for a systematic sample using one random start. As a result, many forest inventory analysts resort to applying variance estimators that are design-unbiased following simple random sampling to their systematic samples even though this typically leads to conservative estimates of error. We explore the influence of spatial autocorrelation on variance estimation when systematic sampling is employed using repeated sampling. We generate a sequence of 1000 synthetic populations with increasing spatial autocorrelation and repeatedly sample from each to examine how the performance of estimators change as spatial autocorrelation changes. We also repeatedly sample from a tree census plot in Harvard Forest, Massachusetts and examine the performance of similar estimators. Results indicate that applying variance estimators that are unbiased following simple random sampling to systematic samples from populations exhibiting stronger spatial autocorrelation tend to be more conservative. We also find that incorporating ancillary wall-to-wall covariates, e.g., remote sensing data, using generalized regression estimators can reduce variance over-estimation by explaining some or all of the spatial structure in the population.


title: Measurements in Interstellar Space of Galactic Cosmic Ray Isotopes of Li, Be, B and N, Ne Nuclei Between 40-160 MeV/nuc by the CRS Instrument on Voyager 1

abstract: In this paper we report a study of the isotopic composition of Li, Be, B and N, Ne nuclei from a 5 year time period beyond the heliopause using the CRS instruments on Voyager. By comparing the isotopic ratios, 15N/14N and 22Ne/20Ne outside the heliosphere as measured at Voyager, and which are found to be significantly lower than those measured at the same energy inside the heliosphere, we have provided strong evidence that cosmic rays of this energy have lost as much as 200 MeV/nuc or more in the solar modulation process. This is in accordance with the so called force field description of this overall modulation by Gleeson and Axford. The measurements at Voyager confirm that the unusual 14N and 22Ne cosmic ray source abundances relative to solar abundances made earlier inside the heliosphere extend to the lower energies not accessible from near Earth measurements. The low energy Li, Be and B nuclei, which are believed to be purely secondary nuclei, are found to have a (previously unobservable) peak in the differential intensity spectrum at ~100 MeV/nuc. This is in agreement with propagation predictions. The intensities of these nuclei are ~10-20% higher than those predicted in a propagation model with a matter path length lambda = 9 g/cm2 at these low energies. The isotopic composition of Li, Be and B nuclei is also consistent with that expected from propagation through interstellar matter.


title: Hypocoercivity for a BGK model for gas mixtures

abstract: We consider a kinetic model for a two component gas mixture without chemical reactions. Our goal is to study hypocoercivity for the linearized BGK model for gas mixtures in continuous phase space. By constructing an entropy functional, we can prove exponential relaxation to equilibrium with explicit rates. Our strategy is based on the entropy and spectral methods adapting Lyapunov's direct method as presented in [Achleitner-Arnold-Carlen, 2016] for the one species linearized BGK model. In comparison to the one species case, we start with two partial differential equations, one for each species. These equations are coupled due to interspecies interactions, which requires additional estimates on these interspecies terms.


title: A Modern Take on the Bias-Variance Tradeoff in Neural Networks

abstract: The bias-variance tradeoff tells us that as model complexity increases, bias falls and variances increases, leading to a U-shaped test error curve. However, recent empirical results with over-parameterized neural networks are marked by a striking absence of the classic U-shaped test error curve: test error keeps decreasing in wider networks. This suggests that there might not be a bias-variance tradeoff in neural networks with respect to network width, unlike was originally claimed by, e.g., Geman et al. (1992). Motivated by the shaky evidence used to support this claim in neural networks, we measure bias and variance in the modern setting. We find that both bias and variance can decrease as the number of parameters grows. To better understand this, we introduce a new decomposition of the variance to disentangle the effects of optimization and data sampling. We also provide theoretical analysis in a simplified setting that is consistent with our empirical findings.


title: Big and nef classes, Futaki Invariant and resolutions of cubic threefolds

abstract: In this note we revisit and extend few classical and recent results on the definition and use of the Futaki invariant in connection with the existence problem for Kaehler constant scalar curvature metrics on polarized algebraic manifolds, especially in the case of resolution of singularities. The general inspiration behind this work is no doubt the beautiful 1992 paper by Ding and Tian which contains the germs of a huge amount of the successive developments in this fundamental problem, and it is a great pleasure to dedicate this to Professor G. Tian on the occasion of his birthday!


title: Transition probabilities for infinite two-sided loop-erased random walks

abstract: The infinite two-sided loop-erased random walk (LERW) is a measure on infinite self-avoiding walks that can be viewed as giving the law of the `middle part' of an infinite LERW loop going through 0 and infinity. In this note we derive expressions for transition probabilities for this model in dimensions two and up. In the plane, the formula can be further expressed in terms of a Laplacian with signed weights acting on certain discrete harmonic functions at the tips of the walk, and taking a determinant. The discrete harmonic functions are closely related to a discrete version of the complex square-root.


title: Theoretical evidence of spin-orbital-entangled $J_{\mathbf{eff}}$=1/2 state in the 3$d$ transition metal oxide CuAl$_2$O$_4$

abstract: Transition metal oxides exhibit various competing phases and exotic phenomena depending on how their reaction to the rich degeneracy of the $d$-orbital. Large spin-orbit coupling (SOC) reduces this degeneracy in a unique way by providing a spin-orbital-entangled ground state for 4$d$ and 5$d$ transition metal compounds. In particular, the spin-orbital-entangled Kramers doublet, known as the $J_{\mathbf{eff}}$=1/2 pseudospin, appears in layered iridates and $\alpha$-RuCl$_3$, manifesting a relativistic Mott insulating phase. Such entanglement, however, seems barely attainable in 3$d$ transition metal oxides, where the SOC is small and the orbital angular momentum is easily quenched. From experimental and theoretical evidence, here we report on the CuAl$_2$O$_4$ spinel as the first example of a $J_{\mathbf{eff}}$=1/2 Mott insulator in 3$d$ transition metal compounds. Based on the experimental study, including synthesis of the cubic CuAl$_2$O$_4$ single crystal, density functional theory and dynamical mean field theory calculations reveal that the $J_{\mathbf{eff}}$=1/2 state survives the competition with an orbital-momentum-quenched $S$=1/2 state. The electron-addition spectra probing unoccupied states are well described by the $j_{\mathbf{eff}}$=1/2 hole state, whereas electron-removal spectra have a rich multiplet structure. The fully relativistic entity found in CuAl$_2$O$_4$ provides new insight into the untapped regime where the spin-orbital-entangled Kramers pair coexists with strong electron correlation.


title: False Discovery and Its Control in Low Rank Estimation

abstract: Models specified by low-rank matrices are ubiquitous in contemporary applications. In many of these problem domains, the row/column space structure of a low-rank matrix carries information about some underlying phenomenon, and it is of interest in inferential settings to evaluate the extent to which the row/column spaces of an estimated low-rank matrix signify discoveries about the phenomenon. However, in contrast to variable selection, we lack a formal framework to assess true/false discoveries in low-rank estimation; in particular, the key source of difficulty is that the standard notion of a discovery is a discrete one that is ill-suited to the smooth structure underlying low-rank matrices. We address this challenge via a geometric reformulation of the concept of a discovery, which then enables a natural definition in the low-rank case. We describe and analyze a generalization of the Stability Selection method of Meinshausen and B\"uhlmann to control for false discoveries in low-rank estimation, and we demonstrate its utility compared to previous approaches via numerical experiments.


title: Template-Based Image Reconstruction from Sparse Tomographic Data

abstract: We propose a variational regularisation approach for the problem of template-based image reconstruction from indirect, noisy measurements as given, for instance, in X-ray computed tomography. An image is reconstructed from such measurements by deforming a given template image. The image registration is directly incorporated into the variational regularisation approach in the form of a partial differential equation that models the registration as either mass- or intensity-preserving transport from the template to the unknown reconstruction. We provide theoretical results for the proposed variational regularisation for both cases. In particular, we prove existence of a minimiser, stability with respect to the data, and convergence for vanishing noise when either of the abovementioned equations is imposed and more general distance functions are used. Numerically, we solve the problem by extending existing Lagrangian methods and propose a multilevel approach that is applicable whenever a suitable downsampling procedure for the operator and the measured data can be provided. Finally, we demonstrate the performance of our method for template-based image reconstruction from highly undersampled and noisy Radon transform data. We compare results for mass- and intensity-preserving image registration, various regularisation functionals, and different distance functions. Our results show that very reasonable reconstructions can be obtained when only few measurements are available and demonstrate that the use of a normalised cross correlation-based distance is advantageous when the image intensities between the template and the unknown image differ substantially.


title: Detecting cities in aerial night-time images by learning structural invariants using single reference augmentation

abstract: This paper examines, if it is possible to learn structural invariants of city images by using only a single reference picture when producing transformations along the variants in the dataset. Previous work explored the problem of learning from only a few examples and showed that data augmentation techniques benefit performance and generalization for machine learning approaches. First a principal component analysis in conjunction with a Fourier transform is trained on a single reference augmentation training dataset using the city images. Secondly a convolutional neural network is trained on a similar dataset with more samples. The findings are that the convolutional neural network is capable of finding images of the same category whereas the applied principal component analysis in conjunction with a Fourier transform failed to solve this task.


title: Immutable Autobiography of Smart Cars

abstract: The popularity of smart cars is increasing around the world as they offer a wide range of services and conveniences.These smart cars are equipped with a variety of sensors generating a large amount of data, many of which are sensitive. Besides, there are multiple parties involved in a lifespan of a smart car ,such as manufacturers, car owners, government agencies, and third-party service providers who also produce data about the vehicle. In addition to managing and sharing data amongst these entities in a secure and privacy-friendly way which is a great challenge itself, there exists a trust deficit about some types of data as they remain under the custody of the car owner(e.g. satellite navigation and mileage data) and can easily be manipulated. In this paper, we propose a blockchain supported architecture enabling the owner of a smart car to create an immutable record of every data, called the auto biography of a car, generated within its lifespan. We also explain how the trust about this record is guaranteed by the immutability characteristic of the blockchain. Furthermore, the paper describes how the proposed architecture enables a secure and privacy-friendly sharing of smart car data between different parties in a secure yet privacy-friendly manner.


title: Equivalence of magnetic field and particle dilution in the strange metal state of CeCoIn$_5$

abstract: The Bardeen-Cooper-Schrieffer mechanism for superconductivity is a triumph of the theory of many-body systems. Implicit in its formulation is the existence of long-lived (quasi)particles, originating from the electronic building blocks of the materials, which interact to form Cooper pairs that move coherently in lock-step. The challenge of unconventional superconductors is that it is not only unclear what the nature of the interactions are, but whether the familiar quasi-particles that form a superconducting condensate even exist. In this work, we reveal, by the study of applied magnetic field in electronically diluted materials, that the metallic properties of the unconventional superconductor CeCoIn$_5$ are determined by the degree of quantum entanglement that (Kondo) hybridizes local and itinerant electrons. This work suggests that the properties of the strange metallic state are a reflection of the disentanglement of the many-body state into the underlying electronic building blocks of the system itself.


title: Transport signatures of surface states in a Weyl semimetal: evidence of field driven Fermi arc interferometry

abstract: A signature property of Weyl semimetals is the existence of topologically protected surface states - arcs in momentum space that connect Weyl points in the bulk. However, the presence of bulks states makes detection of surface contributions to the transport challenging. Here we present a magnetoresistance study of high-quality samples of the prototypical Weyl semimetal, TaAs. By measuring the Shubnikov de Haas effect, we reveal the presence of a two-dimensional cyclotron orbit. This orbit is quantitatively consistent with the interference of coherent quasiparticles traversing two distinct Fermi arcs on the [001] crystallographic surface. The observation of this effect suggests that high magnetic fields can be used to study not only the transport properties of Fermi arcs, but also the interference of their quantum mechanical wavefunctions.


title: Mid-to-far infrared tunable perfect absorption by a sub - \lambda/100 nanofilm in a fractal phasor resonant cavity

abstract: Integrating an absorbing thin film into a resonant cavity is the most practical way to achieve perfect absorption of light at a selected wavelength in the mid-to-far infrared, as required to target blackbody radiation or molecular fingerprints. The cavity is designed to resonate and enable perfect absorption in the film at the chosen wavelength \lambda. However, in current state-of-the-art designs, a still large absorbing film thickness (\lambda/50) is needed and tuning the perfect absorption wavelength over a broad range requires changing the cavity materials. Here, we introduce a new resonant cavity concept to achieve perfect absorption of infrared light in much thinner and thus really nanoscale films, with a broad wavelength tunability using a single set of cavity materials. It requires a nanofilm with giant refractive index and small extinction coefficient (found in emerging semi-metals, semi-conductors and topological insulators) backed by a transparent spacer and a metal mirror. The nanofilm acts both as absorber and multiple reflector for the internal cavity waves, which after escaping follow a fractal phasor trajectory. This enables a totally destructive optical interference for a nanofilm thickness more than 2 orders of magnitude smaller than \lambda. With this remarkable effect, we demonstrate angle-insensitive perfect absorption in sub - \lambda/100 bismuth nanofilms, at a wavelength tunable from 3 to 20 \mu m.


title: A Survey on The Ideal Structure of Leavitt Path Algebras

abstract: There is an extensive recent literature on the graded, non-graded, prime, primitive, maximal ideals of Leavitt path algebras. In this introductory level survey, we will be giving an overview of different types of ideals and the correspondence between the lattice of ideals and the lattice of hereditary and saturated subsets of the graph over which the Leavitt path algebra is constructed.


title: Mainumby: un Ayudante para la Traducci\'on Castellano-Guaran\'i

abstract: A wide range of applications play an important role in the daily work of the modern human translator. However, the computational tools designed to aid in the process of translation only benefit translation from or to a small minority of the 7,000 languages of the world, those that we may call "privileged languages". As for those translators who work with the remaining languages, the marginalized languages in the digital world, they cannot benefit from the tools that are speeding up the production of translation in the privileged languages. We may ask whether it is possible to bridge the gap between what is available for these languages and for the marginalized ones. This paper proposes a framework for computer-assisted translation into marginalized languages and its implementation in a web application for Spanish-Guarani translation. The proposed system is based on a new theory for phrase-level translation in contexts where adequate bilingual corpora are not available: Translation by Generalized Segments (referred to as Minimal Dependency Translation in previous work).


title: Corrosion of Si, C, and SiC in molten salt

abstract: Corrosion of Si, C, and SiC in fluoride salt has been studied by ab initio molecular dynamics. The standard dissolution potential for Si is found to be smaller (easier to corrode) than that of C. The dissolved Si attracts F- ions and forms SiF62-, whereas the dissolved C species forms neutral CF4 molecules. A swapping mechanism is identified for the initial corrosion stage, where Si first comes to the surface and then is dissolved, leaving behind chain- and ring-like C structures. A strategy to suppress SiC corrosion is also discussed based on Be doping, including avoiding Be2C formation.


title: Equilibrium measures and equilibrium potentials in the Born-Infeld model

abstract: In this paper, we consider the electrostatic Born-Infeld model \begin{equation*} \tag{$\mathcal{BI}$} \left\{ \begin{array}{rcll} -\operatorname{div}\left(\displaystyle\frac{\nabla \phi}{\sqrt{1-|\nabla \phi|^2}}\right)&=& \rho & \hbox{in }\mathbb{R}^N, \\[6mm] \displaystyle\lim_{|x|\to \infty}\phi(x)&=& 0 \end{array} \right. \end{equation*} where $\rho$ is a charge distribution on the boundary of a bounded domain $\Omega\subset \mathbb{R}^N$. We are interested in its equilibrium measures, i.e. charge distributions which minimize the electrostatic energy of the corresponding potential among all possible distributions with fixed total charge. We prove existence of equilibrium measures and we show that the corresponding equilibrium potential is unique and constant in $\overline \Omega$. Furthermore, for smooth domains, we obtain the uniqueness of the equilibrium measure, we give its precise expression, and we verify that the equilibrium potential solves ($\mathcal{BI}$). Finally we characterize balls in $\mathbb{R}^N$ as the unique sets among all bounded $C^{2,\alpha}$-domains $\Omega$ for which the equilibrium distribution is a constant multiple of the surface measure on $\partial\Omega$. The same results are obtained also for Taylor approximations of the electrostatic energy.


title: An Exploration of Dropout with RNNs for Natural Language Inference

abstract: Dropout is a crucial regularization technique for the Recurrent Neural Network (RNN) models of Natural Language Inference (NLI). However, dropout has not been evaluated for the effectiveness at different layers and dropout rates in NLI models. In this paper, we propose a novel RNN model for NLI and empirically evaluate the effect of applying dropout at different layers in the model. We also investigate the impact of varying dropout rates at these layers. Our empirical evaluation on a large (Stanford Natural Language Inference (SNLI)) and a small (SciTail) dataset suggest that dropout at each feed-forward connection severely affects the model accuracy at increasing dropout rates. We also show that regularizing the embedding layer is efficient for SNLI whereas regularizing the recurrent layer improves the accuracy for SciTail. Our model achieved an accuracy 86.14% on the SNLI dataset and 77.05% on SciTail.


title: Level Set Methods for Stochastic Discontinuity Detection in Nonlinear Problems

abstract: Stochastic physical problems governed by nonlinear conservation laws are challenging due to solution discontinuities in stochastic and physical space. In this paper, we present a level set method to track discontinuities in stochastic space by solving a Hamilton-Jacobi equation. By introducing a speed function that vanishes at discontinuities, the iso-zero of the level set problem coincide with the discontinuities of the conservation law. The level set problem is solved on a sequence of successively finer grids in stochastic space. The method is adaptive in the sense that costly evaluations of the conservation law of interest are only performed in the vicinity of the discontinuities during the refinement stage. In regions of stochastic space where the solution is smooth, a surrogate method replaces expensive evaluations of the conservation law. The proposed method is tested in conjunction with different sets of localized orthogonal basis functions on simplex elements, as well as frames based on piecewise polynomials conforming to the level set function. The performance of the proposed method is compared to existing adaptive multi-element generalized polynomial chaos methods.


title: Matrix models for topological strings: exact results in the planar limit

abstract: We study the large N expansion of a family of matrix models related to topological strings on toric Calabi-Yau threefolds. These matrix models compute spectral observables of underlying operators obtained by quantizing the mirror curves. They have the form of a deformed O(2) matrix model, with a specific non-polynomial potential involving the Faddeev quantum dilogarithm. Their planar limit is studied using a particular conformal mapping depending on two parameters, from which several universal results can be obtained. As expected, the spectral curves controlling the planar limit of the matrix models are the mirror curves themselves, which in our cases have genus 1. Our results encompass all those toric geometries with genus $1$ mirror where an explicit one-cut matrix integral is known: local $P^2$, local $F_0$, local $F_2$, and degenerations of the resolved $C^3/Z_5$, the resolved $C^3/Z_6$ and the resolved $Y^{3,0}$ geometries amongst others.


title: A Stacked Autoencoder Neural Network based Automated Feature Extraction Method for Anomaly detection in On-line Condition Monitoring

abstract: Condition monitoring is one of the routine tasks in all major process industries. The mechanical parts such as a motor, gear, bearings are the major components of a process industry and any fault in them may cause a total shutdown of the whole process, which may result in serious losses. Therefore, it is very crucial to predict any approaching defects before its occurrence. Several methods exist for this purpose and many research are being carried out for better and efficient models. However, most of them are based on the processing of raw sensor signals, which is tedious and expensive. Recently, there has been an increase in the feature based condition monitoring, where only the useful features are extracted from the raw signals and interpreted for the prediction of the fault. Most of these are handcrafted features, where these are manually obtained based on the nature of the raw data. This of course requires the prior knowledge of the nature of data and related processes. This limits the feature extraction process. However, recent development in the autoencoder based feature extraction method provides an alternative to the traditional handcrafted approaches; however, they have mostly been confined in the area of image and audio processing. In this work, we have developed an automated feature extraction method for on-line condition monitoring based on the stack of the traditional autoencoder and an on-line sequential extreme learning machine(OSELM) network. The performance of this method is comparable to that of the traditional feature extraction approaches. The method can achieve 100% detection accuracy for determining the bearing health states of NASA bearing dataset. The simple design of this method is promising for the easy hardware implementation of Internet of Things(IoT) based prognostics solutions.


title: On tradeoffs between treatment time and plan quality of volumetric-modulated arc therapy with sliding-window delivery

abstract: The purpose of this study is to give an exact formulation of optimization of volumetric-modulated arc therapy (VMAT) with sliding-window delivery, and to investigate the plan quality effects of decreasing the number of sliding-window sweeps made on the 360-degree arc for a faster VMAT treatment. In light of the exact formulation, we interpret an algorithm previously suggested in the literature as a heuristic method for solving this optimization problem. By first making a generalization, we suggest a modification of this algorithm for better handling of plans with fewer sweeps. In a numerical study involving one prostate and one lung case, plans with varying treatment times and number of sweeps are generated. It is observed that, as the treatment time restrictions become tighter, fewer sweeps may lead to better plan quality. Performance of the original and the modified version of the algorithm is evaluated in parallel. Applying the modified version results in better objective function values and less dose discrepancies between optimized and accurate dose, and the advantages are pronounced with decreasing number of sweeps.


title: A database linking piano and orchestral MIDI scores with application to automatic projective orchestration

abstract: This article introduces the Projective Orchestral Database (POD), a collection of MIDI scores composed of pairs linking piano scores to their corresponding orchestrations. To the best of our knowledge, this is the first database of its kind, which performs piano or orchestral prediction, but more importantly which tries to learn the correlations between piano and orchestral scores. Hence, we also introduce the projective orchestration task, which consists in learning how to perform the automatic orchestration of a piano score. We show how this task can be addressed using learning methods and also provide methodological guidelines in order to properly use this database.


title: CNN inference acceleration using dictionary of centroids

abstract: It is well known that multiplication operations in convolutional layers of common CNNs consume a lot of time during inference stage. In this article we present a flexible method to decrease both computational complexity of convolutional layers in inference as well as amount of space to store them. The method is based on centroid filter quantization and outperforms approaches based on tensor decomposition by a large margin. We performed comparative analysis of the proposed method and series of CP tensor decomposition on ImageNet benchmark and found that our method provide almost 2.9 times better computational gain. Despite the simplicity of our method it cannot be applied directly in inference stage in modern frameworks, but could be useful for cases calculation flow could be changed, e.g. for CNN-chip designers.


title: Exact solution for a black hole embedded in a nonstatic dust-filled universe

abstract: An exact solution of the Lema\^{i}tre--Tolman--Bondi class is investigated as a possible model of the Schwarzschild-like black hole embedded in a non-static dust-filled universe for the three types of spatial curvature. The solution is obtained in comoving coordinates by means of the mass function method. It is shown that the central part of space contains a Schwarzschild-like black hole. The R-T-structure of the resulting spacetime is built. It is shown that the solution includes both the Schwarzschild and Friedmann solutions as its natural limits. The geodesic equations for test particles are analyzed. The particle observable velocities are found. The trajectories of the test particles are built from the point of view of both comoving and distant observers. For the distant observer, the results coincide with the Schwarzschild picture within a second-order accuracy near the symmetry center.


title: CoHSI III: Long proteins and implications for protein evolution

abstract: The length distribution of proteins measured in amino acids follows the CoHSI (Conservation of Hartley-Shannon Information) probability distribution. In previous papers we have verified various predictions of this using the Uniprot database but here we explore a novel predicted relationship between the longest proteins and evolutionary time. We demonstrate from both theory and experiment that the longest protein and the total number of proteins are intimately related by Information Theory and we give a simple formula for this. We stress that no evolutionary explanation is necessary; it is an intrinsic property of a CoHSI system. While the CoHSI distribution favors the appearance of proteins with fewer than 750 amino acids (characteristic of most functional proteins or their constituent domains) its intrinsic asymptotic power-law also favors the appearance of unusually long proteins; we predict that there are as yet undiscovered proteins longer than 45,000 amino acids. In so doing, we draw an analogy between the process of protein folding driven by favorable pathways (or funnels) through the energy landscape of protein conformations, and the preferential information pathways through which CoHSI exerts its constraints in discrete systems. Finally, we show that CoHSI predicts the recent appearance in evolutionary time of the longest proteins, specifically in eukaryotes because of their richer unique alphabet of amino acids, and by merging with independent phylogenetic data, we confirm a predicted consistent relationship between the longest proteins and documented and potential undocumented mass extinctions.


title: Autonomous Functional Locomotion in a Tendon-Driven Limb via Limited Experience

abstract: Robots will become ubiquitously useful only when they can use few attempts to teach themselves to perform different tasks, even with complex bodies and in dynamical environments. Vertebrates, in fact, successfully use trial-and-error to learn multiple tasks in spite of their intricate tendon-driven anatomies. Roboticists find such tendon-driven systems particularly hard to control because they are simultaneously nonlinear, under-determined (many tendon tensions combine to produce few net joint torques), and over-determined (few joint rotations define how many tendons need to be reeled-in/payed-out). We demonstrate---for the first time in simulation and in hardware---how a model-free approach allows few-shot autonomous learning to produce effective locomotion in a 3-tendon/2-joint tendon-driven leg. Initially, an artificial neural network fed by sparsely sampled data collected using motor babbling creates an inverse map from limb kinematics to motor activations, which is analogous to juvenile vertebrates playing during development. Thereafter, iterative reward-driven exploration of candidate motor activations simultaneously refines the inverse map and finds a functional locomotor limit-cycle autonomously. This biologically-inspired algorithm, which we call G2P (General to Particular), enables versatile adaptation of robots to changes in the target task, mechanics of their bodies, and environment. Moreover, this work empowers future studies of few-shot autonomous learning in biological systems, which is the foundation of their enviable functional versatility.