title: Most principal permutation classes have nonrational generating functions

abstract: We prove that for any fixed $n$, and for most permutation patterns $q$, the number $\textup{Av}_{n,\ell}(q)$ of $q$-avoiding permutations of length $n$ that consist of $\ell$ skew blocks is a monotone decreasing function of $\ell$. We then show that this implies that for most patterns $q$, the generating function $\sum_{n\geq 0} \textup{Av}_n(q)z^n$ of the sequence $\textup{Av}_n(q)$ of the numbers of $q$-avoiding permutations is not rational. Placing our results in a broader context, we show that for rational power series $F(z)$ and $G(z)$ with nonnegative real coefficients, the relation $F(z)=1/(1-G(z))$ is supercritical, while for most permutation patterns $q$, the corresponding relation is not supercritical.


title: Emergence of leader-follower hierarchy among players in an on-line experiment

abstract: Hierarchical networks are prevalent in nature and society, corresponding to groups of actors - animals, humans or even robots - organised according to a pyramidal structure with decision makers at the top and followers at the bottom. While this phenomenon is seemingly universal, the underlying governing principles are poorly understood. Here we study the emergence of hierarchies in groups of people playing a simple dot guessing game in controlled experiments, lasting for about 40 rounds, conducted over the Internet. During the games, the players had the possibility to look at the answer of a limited number of other players of their choice. This act of asking for advice defines a directed connection between the involved players, and according to our analysis, the initial random configuration of the emerging networks became more structured overt time, showing signs of hierarchy towards the end of the game. In addition, the achieved score of the players appeared to be correlated with their position in the hierarchy. These results indicate that under certain conditions imitation and limited knowledge about the performance of other actors is sufficient for the emergence of hierarchy in a social group.


title: Maximum Entropy Generators for Energy-Based Models

abstract: Maximum likelihood estimation of energy-based models is a challenging problem due to the intractability of the log-likelihood gradient. In this work, we propose learning both the energy function and an amortized approximate sampling mechanism using a neural generator network, which provides an efficient approximation of the log-likelihood gradient. The resulting objective requires maximizing entropy of the generated samples, which we perform using recently proposed nonparametric mutual information estimators. Finally, to stabilize the resulting adversarial game, we use a zero-centered gradient penalty derived as a necessary condition from the score matching literature. The proposed technique can generate sharp images with Inception and FID scores competitive with recent GAN techniques, does not suffer from mode collapse, and is competitive with state-of-the-art anomaly detection techniques.


title: Modal, Truly Counterfactual Communication with On-Chip Demonstration Proposal

abstract: We formalize Salih et al's Counterfactual Communication Protocol (arXiv2018), which allows it not only to be used in with other modes than polarization, but also for interesting extensions (e.g. sending superpositions from Bob to Alice).


title: A priori error estimates for the finite element approximation of Westervelt's quasilinear acoustic wave equation

abstract: We study the spatial discretization of Westervelt's quasilinear strongly damped wave equation by piecewise linear finite elements. Our approach employs the Banach fixed-point theorem combined with a priori analysis of a linear wave model with variable coefficients. Degeneracy of the semi-discrete Westervelt equation is avoided by relying on the inverse estimates for finite element functions and the stability and approximation properties of the interpolation operator. In this way, we obtain optimal convergence rates in $L^2$-based spatial norms for sufficiently small data and mesh size and an appropriate choice of initial approximations. Numerical experiments in a setting of a 1D channel as well as for a focused-ultrasound problem illustrate our theoretical findings.


title: A Unified Analysis of Extra-gradient and Optimistic Gradient Methods for Saddle Point Problems: Proximal Point Approach

abstract: In this paper we consider solving saddle point problems using two variants of Gradient Descent-Ascent algorithms, Extra-gradient (EG) and Optimistic Gradient Descent Ascent (OGDA) methods. We show that both of these algorithms admit a unified analysis as approximations of the classical proximal point method for solving saddle point problems. This viewpoint enables us to develop a new framework for analyzing EG and OGDA for bilinear and strongly convex-strongly concave settings. Moreover, we use the proximal point approximation interpretation to generalize the results for OGDA for a wide range of parameters.


title: Statistics of heat transport across capacitively coupled double quantum dot circuit

abstract: We study heat current and the full statistics of heat fluctuations in a capacitively-coupled double quantum dot system. This work is motivated by recent theoretical studies and experimental works on heat currents in quantum dot circuits. As expected intuitively, within the (static) mean-field approximation, the system at steady-state decouples into two single-dot equilibrium systems with renormalized dot energies, leading to zero average heat flux and fluctuations. This reveals that dynamic correlations induced between electrons on the dots is solely responsible for the heat transport between the two reservoirs. To study heat current fluctuations, we compute steady-state cumulant generating function for heat exchanged between reservoirs using two approaches : Lindblad quantum master equation approach, which is valid for arbitrary coulomb interaction strength but weak system-reservoir coupling strength, and the saddle point approximation for Schwinger-Keldysh coherent state path integral, which is valid for arbitrary system-reservoir coupling strength but weak coulomb interaction strength. Using thus obtained generating functions, we verify steady-state fluctuation theorem for stochastic heat flux and study the average heat current and its fluctuations. We find that the heat current and its fluctuations change non-monotonically with the coulomb interaction strength ($U$) and system-reservoir coupling strength ($\Gamma$) and are suppressed for large values of $U$ and $\Gamma$.

categories: cond-mat.mes-hall cond-mat.other cond-mat.stat-mech quant-ph


title: Finite-Time Stability of Switched and Hybrid Systems with Unstable Modes

abstract: In this work, we study finite-time stability of switched and hybrid systems in the presence of unstable modes. We present sufficient conditions in terms of multiple Lyapunov functions for the origin of the system to be finite time stable. More specifically, we show that even if the value of the Lyapunov function increases in between two switches, i.e., if there are unstable modes in the system, finite-time stability can still be guaranteed if the finite time convergent mode is active long enough. In contrast to earlier work where the Lyapunov functions are required to be decreasing during the continuous flows and non-increasing at the discrete jumps, we allow the Lyapunov functions to increase \emph{both} during the continuous flows and the discrete jumps. As thus, the derived stability results are less conservative compared to the earlier results in the related literature, and in effect allow the hybrid system to have unstable modes. Then, we illustrate how the proposed finite-time stability conditions specialize for a class of switched systems, and present a method on the synthesis of a finite-time stabilizing switching signal for switched linear systems. As a case study, we design a finite-time stable output feedback controller for a linear switched system, in which only one of the modes is controllable and observable. Numerical example demonstrates the efficacy of the proposed methods.


title: Magnetically-Tuned Kondo Effect in a Molecular Double Quantum Dot: Role of the Anisotropic Exchange

abstract: We investigate theoretically and experimentally the singlet-triplet Kondo effect induced by a magnetic field in a molecular junction. Temperature dependent conductance, $G(T)$, is calculated by the numerical renormalization group, showing a strong imprint of the relevant low energy scales, such as the Kondo temperature, exchange and singlet-triplet splitting. We demonstrate the stability of the singlet-triplet Kondo effect against weak spin anisotropy, modeled by an anisotropic exchange. Moderate spin anisotropy manifests itself by lowering the Kondo plateaus, causing the $G(T)$ to deviate from a standard temperature dependence, expected for a spin-half Kondo effect. We propose this scenario as an explanation for anomalous $G(T)$, measured in an organic diradical molecule coupled to gold contacts. We uncover certain new aspects of the singlet-triplet Kondo effect, such as coexistence of spin-polarization on the molecule with Kondo screening and non-perturbative parametric dependence of an effective magnetic field induced by the leads.


title: Global well-posedness for the Phan-Thein-Tanner model in critical Besov spaces without damping

abstract: In this paper, we mainly investigate the Cauchy problem for the Phan-Thein-Tanner (PTT) model. The PPT model can be viewed as a Navier-Stokes equations couple with a nonlinear transport system. This model is derived from network theory for the polymeric fluid. We study about the global well posedness of the PTT model in critical Besov spaces. When the initial data is a small perturbation over around the equilibrium, we prove that the strong solution in critical Besov spaces exists globally.


title: Comparison of post-Newtonian mode amplitudes with numerical relativity simulations of binary black holes

abstract: Gravitational waves from the coalescence of two black holes carry the signature of the strong field dynamics of binary black holes. In this work we have used numerical relativity simulations and post-Newtonian theory to investigate this dynamics. Post-Newtonian theory is a low-velocity expansion that assumes the companion bodies to be point-particles, while numerical relativity treats black holes as extended objects with horizons and fully captures their dynamics. There is a priori no reason for the waveforms computed using these disparate methods to agree with each other, especially at late times when the black holes move close to the speed of light. We find, remarkably, that the leading order amplitudes in post-Newtonian theory agree well with the full general relativity solution for a large set of spherical harmonic modes, even in the most dynamical part of the binary evolution, with only some modes showing distinctly different behavior than that found by numerical relativity simulations. In particular, modes with spherical harmonic indices l = m as well as l = 2, m = 1 are least modified from their dominant post-Newtonian behavior. Understanding the nature of these modes in terms of the post-Newtonian description will aid in formulating better models of the emitted waveforms in the strong field regime of the dynamics.


title: Even simpler modeling of quadruply lensed quasars (and random quartets) using Witt's hyperbola

abstract: Witt (1996) has shown that for an elliptical potential, the four images of a quadruply lensed quasar lie on a rectangular hyperbola that passes through the unlensed quasar position and the center of the potential as well. Wynne and Schechter (2018) have shown that, for the singular isothermal elliptical potential (SIEP), the four images also lie on an `amplitude' ellipse centered on the quasar position with axes parallel to the hyperbola's asymptotes. Witt's hyperbola arises from equating the directions of both sides of the lens equation. The amplitude ellipse derives from equating the magnitudes. One can model any four points as an SIEP in three steps. 1. Find the rectangular hyperbola that passes through the points. 2. Find the aligned ellipse that also passes through them. 3. Find the hyperbola with asymptotes parallel to those of the first that passes through the center of the ellipse and the pair of images closest to each other. The second hyperbola and the ellipse give an SIEP that predicts the positions of the two remaining images where the curves intersect. Pinning the model to the closest pair guarantees a four image model. Such models permit rapid discrimination between gravitationally lensed quasars and random quartets of stars.


title: Learning from Multiple Cities: A Meta-Learning Approach for Spatial-Temporal Prediction

abstract: Spatial-temporal prediction is a fundamental problem for constructing smart city, which is useful for tasks such as traffic control, taxi dispatching, and environmental policy making. Due to data collection mechanism, it is common to see data collection with unbalanced spatial distributions. For example, some cities may release taxi data for multiple years while others only release a few days of data; some regions may have constant water quality data monitored by sensors whereas some regions only have a small collection of water samples. In this paper, we tackle the problem of spatial-temporal prediction for the cities with only a short period of data collection. We aim to utilize the long-period data from other cities via transfer learning. Different from previous studies that transfer knowledge from one single source city to a target city, we are the first to leverage information from multiple cities to increase the stability of transfer. Specifically, our proposed model is designed as a spatial-temporal network with a meta-learning paradigm. The meta-learning paradigm learns a well-generalized initialization of the spatial-temporal network, which can be effectively adapted to target cities. In addition, a pattern-based spatial-temporal memory is designed to distill long-term temporal information (i.e., periodicity). We conduct extensive experiments on two tasks: traffic (taxi and bike) prediction and water quality prediction. The experiments demonstrate the effectiveness of our proposed model over several competitive baseline models.


title: Raking-ratio empirical process with auxiliary information learning

abstract: The raking-ratio method is a statistical and computational method which adjusts the empirical measure to match the true probability of sets of a finite partition. We study the asymptotic behavior of the raking-ratio empirical process indexed by a class of functions when the auxiliary information is given by estimates. We suppose that these estimates result from the learning of the probability of sets of partitions from another sample larger than the sample of the statistician, as in the case of two-stage sampling surveys. Under some metric entropy hypothesis and conditions on the size of the information source sample, we establish the strong approximation of this process and show in this case that the weak convergence is the same as the classical raking-ratio empirical process. We also give possible statistical applications of these results like the strengthening of the $Z$-test and the chi-square goodness of fit test.


title: An Efficient Solver for Cumulative Density Function-based Solutions of Uncertain Kinematic Wave Models

abstract: We develop a numerical framework to implement the cumulative density function (CDF) method for obtaining the probability distribution of the system state described by a kinematic wave model. The approach relies on Monte Carlo Simulations (MCS) of the fine-grained CDF equation of system state, as derived by the CDF method. This fine-grained CDF equation is solved via the method of characteristics. Each method of characteristics solution is far more computationally efficient than the direct solution of the kinematic wave model, and the MCS estimator of the CDF converges relatively quickly. We verify the accuracy and robustness of our procedure via comparison with direct MCS of a particular kinematic wave system, the Saint-Venant equation.


title: The unbalanced reorganization of weaker functional connections induces the altered brain network topology in schizophrenia

abstract: Network neuroscience shed some light on the functional and structural modifications occurring to the brain associated with the phenomenology of schizophrenia. In particular, resting-state functional networks have helped our understanding of the illness by highlighting the global and local alterations within the cerebral organization. We investigated the robustness of the brain functional architecture in forty-four medicated schizophrenic patients and forty healthy comparators through an advanced network analysis of resting-state functional magnetic resonance imaging data. The networks in patients showed more resistance to disconnection than in healthy controls, with an evident discrepancy between the two groups in the node degree distribution computed along a percolation process. Despite a substantial similarity of the basal functional organization between the two groups, the expected hierarchy of healthy brains modular organization is crumbled in schizophrenia, showing a peculiar arrangement of the functional connections, characterized by several topologically equivalent backbones.


title: Mixed-Granularity Human-Swarm Interaction

abstract: We present an augmented reality human-swarm interface that combines two modalities of interaction: environment-oriented and robot-oriented. The environment-oriented modality allows the user to modify the environment (either virtual or physical) to indicate a goal to attain for the robot swarm. The robot-oriented modality makes it possible to select individual robots to reassign them to other tasks to increase performance or remedy failures. Previous research has concluded that environment-oriented interaction might prove more difficult to grasp for untrained users. In this paper, we report a user study which indicates that, at least in collective transport, environment-oriented interaction is more effective than purely robot-oriented interaction, and that the two combined achieve remarkable efficacy.


title: Curvature-Exploiting Acceleration of Elastic Net Computations

abstract: This paper introduces an efficient second-order method for solving the elastic net problem. Its key innovation is a computationally efficient technique for injecting curvature information in the optimization process which admits a strong theoretical performance guarantee. In particular, we show improved run time over popular first-order methods and quantify the speed-up in terms of statistical measures of the data matrix. The improved time complexity is the result of an extensive exploitation of the problem structure and a careful combination of second-order information, variance reduction techniques, and momentum acceleration. Beside theoretical speed-up, experimental results demonstrate great practical performance benefits of curvature information, especially for ill-conditioned data sets.


title: The Tensor Pomeron and Low-x Deep Inelastic Scattering

abstract: The tensor-pomeron model is applied to low-x deep-inelastic lepton-nucleon scattering and photoproduction. We consider c.m. energies in the range 6 - 318 GeV and Q^2 < 50 GeV^2. In addition to the soft tensor pomeron, which has proven quite successful for the description of soft hadronic high-energy reactions, we include a hard tensor pomeron. We also include f_2-reggeon exchange which turns out to be particularly relevant for real-photon-proton scattering at c.m. energies in the range up to 30 GeV. The combination of these exchanges permits a description of the absorption cross sections of real and virtual photons on the proton in the same framework. In particular, a detailed comparison of this two-tensor-pomeron model with the latest HERA data for x < 0.01 is made. Our model gives a very good description of the transition from the small-Q^2 regime where the real or virtual photon behaves hadron-like to the large-Q^2 regime where hard scattering dominates. Our fit allows us, for instance, a determination of the intercepts of the hard pomeron as 1.3008(+73,-84), of the soft pomeron as 1.0935(+76,-64), and of the f_2 reggeon. We find that in photoproduction the hard pomeron does not contribute within the errors of the fit. We show that assuming a vector instead of a tensor character for the pomeron leads to the conclusion that it must decouple in real photoproduction.


title: The self-concordant perceptron is efficient on a sub-family feasibility instances

abstract: Strict linear feasibility or linear separation is usually tackled using efficient approximation/stochastic algorithms (that may even run in sub-linear times in expectation). However, today state of the art for solving exactly/deterministically such instances is to cast them as a linear programming instances. Inversely, this paper introduces a self-concordant perceptron algorithm which tackles directly strict linear feasibility with interior point paradigm. This algorithm has the same worse times complexity than state of the art linear programming algorithms but it complexity can be characterized more precisely eventually proving that it binary complexity is low on a sub-family of linear feasibility.


title: IR-truncated $\mathcal{PT}-$symmetric $ix^3$ model and its asymptotic spectral scaling graph

abstract: The $\mathcal{PT}-$symmetric quantum mechanical $V=ix^3$ model over the real line, $x\in\mathbb{R}$, is infrared (IR) truncated and considered as Sturm-Liouville problem over a finite interval $x\in\left[-L,L\right]\subset\mathbb{R}$. Via WKB and Stokes graph analysis, the location of the complex spectral branches of the $V=ix^3$ model and those of more general $V=-(ix)^{2n+1}$ models over $x\in\left[-L,L\right]\subset\mathbb{R}$ are obtained. The corresponding eigenvalues are mapped onto $L-$invariant asymptotic spectral scaling graphs $\mathcal{R}\subset \mathbb{C}$. These scaling graphs are geometrically invariant and cutoff-independent so that the IR limit $L\to \infty $ can be formally taken. Moreover, an increasing $L$ can be associated with an $\mathcal{R}-$constrained spectral UV$\to$IR renormalization group flow on $\mathcal{R}$. The existence of a scale-invariant $\mathcal{PT}$ symmetry breaking region on each of these graphs allows to conclude that the unbounded eigenvalue sequence of the $ix^3$ Hamiltonian over $x\in\mathbb{R}$ can be considered as tending toward a mapped version of such a $\mathcal{PT}$ symmetry breaking region at spectral infinity. This provides a simple heuristic explanation for the specific eigenfunction properties described in the literature so far and clear complementary evidence that the $\mathcal{PT}-$symmetric $V=-(ix)^{2n+1}$ models over the real line $x\in\mathbb{R}$ are not equivalent to Hermitian models, but that they rather form a separate model class with purely real spectra. Our findings allow us to hypothesize a possible physical interpretation of the non-Rieszian mode behavior as a related mode condensation process.

categories: math-ph cond-mat.stat-mech hep-th math.MP math.SP quant-ph


title: Probing typical black hole microstates

abstract: We investigate the possibility that the geometry dual to a typical AdS black hole microstate corresponds to the extended AdS-Schwarzschild geometry, including a region spacelike to the exterior. We argue that this region can be described by the mirror operators, a set of state-dependent operators in the dual CFT. We probe the geometry of a typical state by considering state-dependent deformations of the CFT Hamiltonian, which have an interpretation as a one-sided analogue of the Gao-Jafferis-Wall traversable wormhole protocol for typical states. We argue that the validity of the conjectured bulk geometry requires that out-of-time-order correlators of simple CFT operators on typical pure states must exhibit the same chaotic effects as thermal correlators at scrambling time. This condition is related to the question of whether the product of operators separated by scrambling time obey the Eigenstate Thermalization Hypothesis. We investigate some of these statements in the SYK model and discuss similarities with state-dependent perturbations of pure states in the SYK model previously considered by Kourkoulou and Maldacena. Finally, we discuss how the mirror operators can be used to implement an analogue of the Hayden-Preskill protocol.


title: WIMPs and stellar-mass primordial black holes are incompatible

abstract: We recently showed that postulated ultracompact minihalos with a steep density profile do not form in realistic simulations with enhanced initial perturbations. In this paper we assume that a small fraction of the dark matter consists of primordial black holes (PBHs) and simulate the formation of structures around them. We find that in this scenario halos with steep density profiles do form, consistent with theoretical predictions. If the rest of the dark matter consists of weakly interacting massive particles (WIMPs), we also show that WIMPs in the dense innermost part of these halos would annihilate and produce a detectable gamma-ray signal. The non-detection of this signal implies that PBHs make up at most one billionth of the dark matter, provided that their mass is greater than one millionth of the mass of the Sun. Similarly, a detection of PBHs would imply that the remaining dark matter could not be WIMPs.


title: Inequalities for some integrals involving modified Lommel functions of the first kind

abstract: In this paper, we obtain inequalities for some integrals involving the modified Lommel function of the first kind $t_{\mu,\nu}(x)$. In most cases, these inequalities are tight in certain limits. We also deduce a tight double inequality, involving the modified Lommel function $t_{\mu,\nu}(x)$, for a generalized hypergeometric function. The inequalities obtained in this paper generalise recent bounds for integrals involving the modified Struve function of the first kind.


title: Relativistic dissipation obeys Chapman-Enskog asymptotics: analytical and numerical evidence as a basis for accurate kinetic simulations

abstract: We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grad's expansion methods. Moreover, we develop a systematic calibration method, connecting the relaxation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison between the analytical results and numerical simulations, shows that the CE method correctly captures dissipative effects, while Grad's method does not. The resulting calibration procedure based on the CE method opens the way to the quantitative kinetic description of dissipative relativistic fluid dynamics under fairly general conditions, namely flows with strongly non-linearities, in non-ideal geometries, across both ultra-relativistic and near-non-relativistic regimes.


title: Coronal Imaging with the Solar UltraViolet Imager

abstract: We investigate the coronal imaging capabilities of the Solar UltraViolet Imager (SUVI) on the Geostationary Operational Environmental Satellite-R series spacecraft. Nominally Sun-pointed, SUVI provides solar images in six Extreme UltraViolet (EUV) wavelengths. On-orbit data indicated that SUVI had sufficient dynamic range and sensitivity to image the corona to the largest heights above the Sun to date while simultaneously imaging the Sun. We undertook a campaign to investigate the existence of the EUV signal well beyond the nominal Sun-centered imaging area of the solar EUV imagers. We off-pointed SUVI line-of-sight by almost one imaging area around the Sun. We present the details of the campaign conducted when the solar cycle is at near the minimum and some results that affirm the EUV presence to beyond three solar radii.


title: Direct $N$-body code on low-power embedded ARM GPUs

abstract: This work arises on the environment of the ExaNeSt project aiming at design and development of an exascale ready supercomputer with low energy consumption profile but able to support the most demanding scientific and technical applications. The ExaNeSt compute unit consists of densely-packed low-power 64-bit ARM processors, embedded within Xilinx FPGA SoCs. SoC boards are heterogeneous architecture where computing power is supplied both by CPUs and GPUs, and are emerging as a possible low-power and low-cost alternative to clusters based on traditional CPUs. A state-of-the-art direct $N$-body code suitable for astrophysical simulations has been re-engineered in order to exploit SoC heterogeneous platforms based on ARM CPUs and embedded GPUs. Performance tests show that embedded GPUs can be effectively used to accelerate real-life scientific calculations, and that are promising also because of their energy efficiency, which is a crucial design in future exascale platforms.


title: A deep XMM-Newton look on the thermally emitting isolated neutron star RX J1605.3+3249

abstract: Previous XMM-Newton observations of the thermally emitting isolated neutron star RX J1605.3+3249 provided a candidate for a shallow periodic signal and evidence of a fast spin down, which suggested a high dipolar magnetic field and an evolution from a magnetar. We obtained a large programme with XMM-Newton to confirm its candidate timing solution, understand the energy-dependent amplitude of the modulation, and investigate the spectral features of the source. We performed extensive high-resolution and broadband periodicity searches in the new observations, using the combined photons of the three EPIC cameras and allowing for moderate changes of pulsed fraction and the optimal energy range for detection. A deep $4\sigma$ upper limit of $1.33(6)\%$ for modulations in the relevant frequency range conservatively rules out the candidate period previously reported. Blind searches revealed no other periodic signal above the $1.5\%$ level $(3\sigma$) in any of the four new observations. While theoretical models fall short at physically describing the complex energy distribution of the source, best-fit X-ray spectral parameters are obtained for a fully or partially ionized neutron star hydrogen atmosphere model with $B=10^{13}$ G, modified by a broad Gaussian absorption line at $385\pm10$ eV. The deep limits from the timing analysis disfavour equally well-fit double temperature blackbody models where both the star surface and small hotspots contribute to the X-ray flux of the source. We identified a low significance ($1\sigma$) temporal trend on the parameters of the source in the analysis of RGS data dating back to 2002, which may be explained by unaccounted calibration issues and spectral model uncertainties. The new dataset also shows no evidence of the previously reported narrow absorption feature at $\sim570$ eV, whose possible transient nature disfavours an atmospheric origin.


title: Learning Disentangled Representations with Reference-Based Variational Autoencoders

abstract: Learning disentangled representations from visual data, where different high-level generative factors are independently encoded, is of importance for many computer vision tasks. Solving this problem, however, typically requires to explicitly label all the factors of interest in training images. To alleviate the annotation cost, we introduce a learning setting which we refer to as "reference-based disentangling". Given a pool of unlabeled images, the goal is to learn a representation where a set of target factors are disentangled from others. The only supervision comes from an auxiliary "reference set" containing images where the factors of interest are constant. In order to address this problem, we propose reference-based variational autoencoders, a novel deep generative model designed to exploit the weak-supervision provided by the reference set. By addressing tasks such as feature learning, conditional image generation or attribute transfer, we validate the ability of the proposed model to learn disentangled representations from this minimal form of supervision.


title: Testing Equality of Autocovariance Operators for Functional Time Series

abstract: We consider strictly stationary stochastic processes of Hilbert space-valued random variables and focus on fully functional tests for the equality of the lag-zero autocovariance operators of several independent functional time series. A moving block bootstrap-based testing procedure is proposed which generates pseudo random elements that satisfy the null hypothesis of interest. It is based on directly bootstrapping the time series of tensor products which overcomes some common difficulties associated with applications of the bootstrap to related testing problems. The suggested methodology can be potentially applied to a broad range of test statistics of the hypotheses of interest. As an example, we establish validity for approximating the distribution under the null of a test statistic based on the Hilbert-Schmidt distance of the corresponding sample lag-zero autocovariance operators, and show consistency under the alternative. As a prerequisite, we prove a central limit theorem for the moving block bootstrap procedure applied to the sample autocovariance operator which is of interest on its own. The finite sample size and power performance of the suggested moving block bootstrap-based testing procedure is illustrated through simulations and an application to a real-life dataset is discussed.


title: Dynamic axion field in the magnetoelectric antiferromagnet chromia

abstract: Chromia is a magnetoelectric insulator whose bulk magnetoelectric susceptibility contains a non-zero pseudoscalar component, like that present in magnetized topological insulators. We reveal the dynamic axion field of chromia by measuring the temperature dependence of its non-linear magnetoelectric response using a lock-in technique with an AC electric field stimulus. The electric field dependence, magnetic field dependence, and temperature behavior of the second harmonic lock-in signal are all in agreement with theoretical expectations. Our results demonstrate that chromia magnetoelectric response couples strongly to its low-energy quantum and thermal spin fluctuations.


title: Distributed Learning of Decentralized Control Policies for Articulated Mobile Robots

abstract: State-of-the-art distributed algorithms for reinforcement learning rely on multiple independent agents, which simultaneously learn in parallel environments while asynchronously updating a common, shared policy. Moreover, decentralized control architectures (e.g., CPGs) can coordinate spatially distributed portions of an articulated robot to achieve system-level objectives. In this work, we investigate the relationship between distributed learning and decentralized control by learning decentralized control policies for the locomotion of articulated robots in challenging environments. To this end, we present an approach that leverages the structure of the asynchronous advantage actor-critic (A3C) algorithm to provide a natural means of learning decentralized control policies on a single articulated robot. Our primary contribution shows individual agents in the A3C algorithm can be defined by independently controlled portions of the robot's body, thus enabling distributed learning on a single robot for efficient hardware implementation. We present results of closed-loop locomotion in unstructured terrains on a snake and a hexapod robot, using decentralized controllers learned offline and online respectively. Preprint of the paper submitted to the IEEE Transactions in Robotics (T-RO) journal in October 2018, and accepted for publication as a regular paper in May 2019.


title: Fluctuation bounds for ergodic averages of amenable groups on uniformly convex Banach spaces

abstract: We study fluctuations of ergodic averages generated by actions of amenable groups. In the setting of an abstract ergodic theorem for locally compact second countable amenable groups acting on uniformly convex Banach spaces, we deduce a highly uniform bound on the number of fluctuations of the ergodic average for a class of F{\o}lner sequences satisfying an analogue of Lindenstrauss's temperedness condition. Equivalently, we deduce a uniform bound on the number of fluctuations over long distances for arbitrary F{\o}lner sequences. As a corollary, these results imply associated bounds for a continuous action of an amenable group on a $\sigma$-finite $L^{p}$ space with $p\in(1,\infty)$.


title: Multiresolution Analysis and Learning for Computational Seismic Interpretation

abstract: We explore the use of multiresolution analysis techniques as texture attributes for seismic image characterization, especially in representing subsurface structures in large migrated seismic data. Namely, we explore the Gaussian pyramid, the discrete wavelet transform, Gabor filters, and the curvelet transform. These techniques are examined in a seismic structure labeling case study on the Netherlands offshore F3 block. In seismic structure labeling, a seismic volume is automatically segmented and classified according to the underlying subsurface structure using texture attributes. Our results show that multiresolution attributes improved the labeling performance compared to using seismic amplitude alone. Moreover, directional multiresolution attributes, such as the curvelet transform, are more effective than the non-directional attributes in distinguishing different subsurface structures in large seismic datasets, and can greatly help the interpretation process.


title: Causal Mediation Analysis Leveraging Multiple Types of Summary Statistics Data

abstract: Summary statistics of genome-wide association studies (GWAS) teach causal relationship between millions of genetic markers and tens and thousands of phenotypes. However, underlying biological mechanisms are yet to be elucidated. We can achieve necessary interpretation of GWAS in a causal mediation framework, looking to establish a sparse set of mediators between genetic and downstream variables, but there are several challenges. Unlike existing methods rely on strong and unrealistic assumptions, we tackle practical challenges within a principled summary-based causal inference framework. We analyzed the proposed methods in extensive simulations generated from real-world genetic data. We demonstrated only our approach can accurately redeem causal genes, even without knowing actual individual-level data, despite the presence of competing non-causal trails.


title: Localizations of Fast Radio Bursts on milliarcsecond scales

abstract: Fast Radio Bursts (FRBs) are transient sources that emit a single radio pulse with a duration of only a few milliseconds. Since the discovery of the first FRB in 2007, tens of similar events have been detected. However, their physical origin remains unclear, and a number of scenarios even larger than the number of known FRBs has been proposed during these years. The presence of repeating bursts in FRB 121102 allowed us to perform a precise localization of the source with the Very Large Array and the European VLBI Network (EVN). Optical observations with Keck, Gemini and HST unveiled the host to be a low-metallicity star-forming dwarf galaxy located at a redshift of 0.193. The EVN results showed that the bursts are co-located (within a projected separation of $< 40$ pc) to a compact and persistent radio source with a size of $< 0.7$ pc inside a star-forming region. This environment resembles the ones where superluminous supernovae (SLSNe) or long gamma-ray bursts are produced. Although the nature of this persistent source and the origin of the bursts remain unknown, scenarios considering a pulsar/magnetar energizing a young SLSN, or a system with a pulsar/magnetar in the vicinity of a massive black hole are the most plausible ones to date. More recent observations have shown that the bursts from FRB 121102 are almost 100% linearly polarized at an unexpectedly high and variable Faraday rotation measure, that has been observed to date only in vicinities of massive black holes. The bursts are thus likely produced from a neutron star in such environment, although the system can still be explained by a young neutron star embedded in a highly magnetized nebula. Upcoming interferometric searches are expected to report tens of these localizations in the coming years, unveil if this source is representative of the whole population or a particular case, and dramatically boosting the field of FRBs.


title: Origin of life's building blocks in Carbon and Nitrogen rich surface hydrothermal vents

abstract: There are two dominant and contrasting classes of origin of life scenarios: those predicting that life emerged in submarine hydrothermal systems, where chemical disequilibrium can provide an energy source for nascent life; and those predicting that life emerged within subaerial environments, where UV catalysis of reactions may occur to form the building blocks of life. Here, we describe a prebiotically plausible environment that draws on the strengths of both scenarios: surface hydrothermal vents. We show how key feedstock molecules for prebiotic chemistry can be produced in abundance in shallow and surficial hydrothermal systems. We calculate the chemistry of volcanic gases feeding these vents over a range of pressures and basalt C/N/O contents. If ultra-reducing carbon-rich nitrogen-rich gases interact with subsurface water at a volcanic vent they result in 1 mM to 1 M concentrations of diacetylene, acetylene, cyanoacetylene, hydrogen cyanide, bisulfite, hydrogen sulfide and soluble iron in vent water. One key feedstock molecule, cyanamide, is not formed in significant quantities within this scenario, suggesting that it may need to be delivered exogenously, or formed from hydrogen cyanide either via organometallic compounds, or by some as yet-unknown chemical synthesis. Given the likely ubiquity of surface hydrothermal vents on young, hot, terrestrial planets, these results identify a prebiotically plausible local geochemical environment, which is also amenable to future lab-based simulation.


title: Anisotropic magnetic excitations and incipient N\'eel order in Ba(Fe$_{1-x}$Mn$_{x}$)$_{2}$As$_{2}$

abstract: It is currently understood that high temperature superconductivity (SC) in the transition metal $(M)$ substituted iron arsenides Ba(Fe$_{1-x}$$M$$_{x}$)$_{2}$As$_{2}$ is promoted by magnetic excitations with wave vectors $(\pi,0)$ or $(0,\pi)$. It is known that while a small amount of Co substitution leads to SC, the same does not occur for Mn for any value of $x$. In this work, magnetic excitations in the iron arsenides Ba(Fe$_{1-x}$Mn$_{x}$)$_{2}$As$_{2}$ ($x=0.0$, $0.007$, $0.009$, $0.08$) are investigated by means of Resonant Inelastic X rays Scattering (RIXS) at the Fe $L_{3}$-edge, for momentum transfer $\boldsymbol{q}$ along the high symmetry Brillouin zone $(\pi,0)$ and $(\pi,\pi)$ directions. It is shown that with increasing Mn content ($x$), the excitations become anisotropic both in dispersion and lineshape. Both effects are detected even for small values of $x$, evidencing a cooperative phenomenon between the Mn impurities, that we ascribe to emerging N\'eel order of the Mn spins. Moreover, for $x=0.08$, the excitations along $\boldsymbol{q}\parallel(\pi,0)$ are strongly damped and nearly non dispersive. This result suggests that phases of arsenides containing local moments at the FeAs layers, as in Mn or Cr substituted phases, do not support high temperature SC due to absence of the appropriate magnetic excitations.


title: Learning Space Partitions for Nearest Neighbor Search

abstract: Space partitions of $\mathbb{R}^d$ underlie a vast and important class of fast nearest neighbor search (NNS) algorithms. Inspired by recent theoretical work on NNS for general metric spaces [Andoni, Naor, Nikolov, Razenshteyn, Waingarten STOC 2018, FOCS 2018], we develop a new framework for building space partitions reducing the problem to balanced graph partitioning followed by supervised classification. We instantiate this general approach with the KaHIP graph partitioner [Sanders, Schulz SEA 2013] and neural networks, respectively, to obtain a new partitioning procedure called Neural Locality-Sensitive Hashing (Neural LSH). On several standard benchmarks for NNS, our experiments show that the partitions obtained by Neural LSH consistently outperform partitions found by quantization-based and tree-based methods as well as classic, data-oblivious LSH.


title: Dating the foundation of Augusta Taurinorum ex sole. The augustean propaganda and the role of Astronomy

abstract: The essay presents the results of a joint study of astronomy and archeology that has allowed to define the foundation date of the city of Turin as a Roman colony, called Iulia Augusta Taurinorum. This multidisciplinary research represents a new reading of the historical-archaeological sources and the use of astronomy according to the Etruscan-Latin gromatica. By taking into consideration the apparent motion of a True Sun, the possible measurement errors, the atmospheric refraction and the elevation of the horizon, and the Julian date in use in astronomy, a numerical program has been elaborated to define the coincidences of the calendars between the azimuth of the main road axis and the course of the sun. For a series of very particular historical and contextual conditions it was therefore possible to trace with sufficient accuracy the day and the year of foundation of the city: January 30, 9 BC, which coincides, not surprisingly, with a particularly important anniversary for Ottaviano Augusto, Emperor from 27 BC to 14 AC. In fact, the astronomical data interfaced with the archaeological, epigraphic and written sources suggest that Turin was born at the end of the Alpine wars, concluded after the acquisition of the Valle di Susa and, more extensively, of the Maritime Alps, Pennine and Graie. Founded as a colony along the road leading to the Gauls, Augusta Taurinorum was inaugurated on the day of the anniversary of the feast of Pax, established by Augustus and celebrated at the Ara Pacis in Rome starting from 9 BC. The large marble altar of Campo Marzio was, as matter of fact, built to celebrate the end of civil wars and the beginning of a new era (confirmed also by the astrological/astronomical ephemerides of that time), which is widely reported in the official Augustean literature.


title: Ballistic and hydrodynamic magnetotransport in narrow channels

abstract: An increasing number of low carrier density materials exhibit a surprisingly large transport mean free path due to inefficient momentum relaxation. Consequently, charge transport in these systems is markedly non-ohmic but rather ballistic or hydrodynamic, features which can be explored by driving current through narrow channels. Using a kinetic equation approach we theoretically investigate how a non-quantizing magnetic field discerns ballistic and hydrodynamic transport, in particular in the spatial dependence of the transverse electric field, $E_y$: We find that $E_y$ is locally enhanced when the flow exhibits a sharp directional anisotropy in the non-equilibrium density. As a consequence, at weak magnetic fields, the curvature of $E_y$ has opposite signs in the ballistic and hydrodynamic regimes. Moreover, we find a robust signature of the onset of non-local correlations in the form of distinctive peaks of the transverse field, which are accessible by local measurements. Our results demonstrate that a purely hydrodynamic approach is insufficient in the Gurzhi regime once a magnetic field is introduced.


title: Human-centric Transfer Learning Explanation via Knowledge Graph [Extended Abstract]

abstract: Transfer learning which aims at utilizing knowledge learned from one problem (source domain) to solve another different but related problem (target domain) has attracted wide research attentions. However, the current transfer learning methods are mostly uninterpretable, especially to people without ML expertise. In this extended abstract, we brief introduce two knowledge graph (KG) based frameworks towards human understandable transfer learning explanation. The first one explains the transferability of features learned by Convolutional Neural Network (CNN) from one domain to another through pre-training and fine-tuning, while the second justifies the model of a target domain predicted by models from multiple source domains in zero-shot learning (ZSL). Both methods utilize KG and its reasoning capability to provide rich and human understandable explanations to the transfer procedure.


title: A tensorized logic programming language for large-scale data

abstract: We introduce a new logic programming language T-PRISM based on tensor embeddings. Our embedding scheme is a modification of the distribution semantics in PRISM, one of the state-of-the-art probabilistic logic programming languages, by replacing distribution functions with multidimensional arrays, i.e., tensors. T-PRISM consists of two parts: logic programming part and numerical computation part. The former provides flexible and interpretable modeling at the level of first order logic, and the latter part provides scalable computation utilizing parallelization and hardware acceleration with GPUs. Combing these two parts provides a remarkably wide range of high-level declarative modeling from symbolic reasoning to deep learning. To embody this programming language, we also introduce a new semantics, termed tensorized semantics, which combines the traditional least model semantics in logic programming with the embeddings of tensors. In T-PRISM, we first derive a set of equations related to tensors from a given program using logical inference, i.e., Prolog execution in a symbolic space and then solve the derived equations in a continuous space by TensorFlow. Using our preliminary implementation of T-PRISM, we have successfully dealt with a wide range of modeling. We have succeeded in dealing with real large-scale data in the declarative modeling. This paper presents a DistMult model for knowledge graphs using the FB15k and WN18 datasets.


title: Enabling Deep All-Sky Searches of Outer Solar System Objects

abstract: A foundational goal of the Large Synoptic Survey Telescope (LSST) is to map the Solar System small body populations that provide key windows into understanding of its formation and evolution. This is especially true of the populations of the Outer Solar System -- objects at the orbit of Neptune $r > 30$AU and beyond. In this whitepaper, we propose a minimal change to the LSST cadence that can greatly enhance LSST's ability to discover faint distant Solar System objects across the entire wide-fast-deep (WFD) survey area. Specifically, we propose that the WFD cadence be constrained so as to deliver least one sequence of $\gtrsim 10$ visits per year taken in a $\sim 10$ day period in any combination of $g, r$, and $i$ bands. Combined with advanced shift-and-stack algorithms (Whidden et al. 2019) this modification would enable a nearly complete census of the outer Solar System to $\sim 25.5$ magnitude, yielding $4-8$x more KBO discoveries than with single-epoch baseline, and enabling rapid identification and follow-up of unusual distant Solar System objects in $\gtrsim 5$x greater volume of space. These increases would enhance the science cases discussed in Schwamb et al. (2018) whitepaper, including probing Neptune's past migration history as well as discovering hypothesized planet(s) beyond the orbit of Neptune (or at least placing significant constraints on their existence).


title: On the K-theory of coordinate axes in affine space

abstract: Let k be a perfect field of characteristic p>0, let A_d be the coordinate ring of the coordinate axes in affine d-space over k, and let I_d be the ideal defining the origin. We evaluate the relative K-groups K_q(A_d,I_d) in terms of p-typical Witt vectors of k. When d=2 the result is due to Hesselholt, and for K_2 it is due to Dennis and Krusemeyer. We also compute the groups K_q(A_d,I_d) in the case where k is an ind-smooth algebra over the rationals, the result being expressed in terms of algebraic de Rham forms. When k is a field of characteristic zero this calculation is due to Geller, Reid and Weibel.


title: A Universal Logic Operator for Interpretable Deep Convolution Networks

abstract: Explaining neural network computation in terms of probabilistic/fuzzy logical operations has attracted much attention due to its simplicity and high interpretability. Different choices of logical operators such as AND, OR and XOR give rise to another dimension for network optimization, and in this paper, we study the open problem of learning a universal logical operator without prescribing to any logical operations manually. Insightful observations along this exploration furnish deep convolution networks with a novel logical interpretation.


title: General Supervision via Probabilistic Transformations

abstract: Different types of training data have led to numerous schemes for supervised classification. Current learning techniques are tailored to one specific scheme and cannot handle general ensembles of training data. This paper presents a unifying framework for supervised classification with general ensembles of training data, and proposes the learning methodology of generalized robust risk minimization (GRRM). The paper shows how current and novel supervision schemes can be addressed under the proposed framework by representing the relationship between examples at test and training via probabilistic transformations. The results show that GRRM can handle different types of training data in a unified manner, and enable new supervision schemes that aggregate general ensembles of training data.


title: Data Interpolations in Deep Generative Models under Non-Simply-Connected Manifold Topology

abstract: Exploiting the deep generative model's remarkable ability of learning the data-manifold structure, some recent researches proposed a geometric data interpolation method based on the geodesic curves on the learned data-manifold. However, this interpolation method often gives poor results due to a topological difference between the model and the dataset. The model defines a family of simply-connected manifolds, whereas the dataset generally contains disconnected regions or holes that make them non-simply-connected. To compensate this difference, we propose a novel density regularizer that make the interpolation path circumvent the holes denoted by low probability density. We confirm that our method gives consistently better interpolation results from the experiments with real-world image datasets.


title: Explainable Failure Predictions with RNN Classifiers based on Time Series Data

abstract: Given key performance indicators collected with fine granularity as time series, our aim is to predict and explain failures in storage environments. Although explainable predictive modeling based on spiky telemetry data is key in many domains, current approaches cannot tackle this problem. Deep learning methods suitable for sequence modeling and learning temporal dependencies, such as RNNs, are effective, but opaque from an explainability perspective. Our approach first extracts the anomalous spikes from time series as events and then builds an RNN classifier with attention mechanisms to embed the irregularity and frequency of these events. A preliminary evaluation on real world storage environments shows that our approach can predict failures within a 3-day prediction window with comparable accuracy as traditional RNN-based classifiers. At the same time it can explain the predictions by returning the key anomalous events which led to those failure predictions.


title: Light-cone reduction of Witten's open string field theory

abstract: We elucidate some exact relations between light-cone and covariant string field theories on the basis of the homological perturbation lemma for $A_{\infty }$. The covariant string field splits into the light-cone string field and trivial excitations of BRST quartets: The latter generates the gauge symmetry and covariance. We first show that the reduction of gauge degrees can be performed by applying the lemma, which gives a refined version of the no-ghost theorem of covariant strings. Then, we demonstrate that after the reduction, gauge-fixed theory can be regarded as a kind of effective field theory and it provides an exact gauge-fixing procedure taking into account interactions. As a result, a novel light-cone string field theory is obtained from Witten's open string field theory.


title: Visualized Insights into the Optimization Landscape of Fully Convolutional Networks

abstract: Many image processing tasks involve image-to-image mapping, which can be addressed well by fully convolutional networks (FCN) without any heavy preprocessing. Although empirically designing and training FCNs can achieve satisfactory results, reasons for the improvement in performance are slightly ambiguous. Our study is to make progress in understanding their generalization abilities through visualizing the optimization landscapes. The visualization of objective functions is obtained by choosing a solution and projecting its vicinity onto a 3D space. We compare three FCN-based networks (two existing models and a new proposed in this paper for comparison) on multiple datasets. It has been observed in practice that the connections from the pre-pooled feature maps to the post-upsampled can achieve better results. We investigate the cause and provide experiments to shows that the skip-layer connections in FCN can promote flat optimization landscape, which is well known to generalize better. Additionally, we explore the relationship between the models generalization ability and loss surface under different batch sizes. Results show that large-batch training makes the model converge to sharp minimizers with chaotic vicinities while small-batch method leads the model to flat minimizers with smooth and nearly convex regions. Our work may contribute to insights and analysis for designing and training FCNs.


title: On Network Science and Mutual Information for Explaining Deep Neural Networks

abstract: In this paper, we present a new approach to interpret deep learning models. By coupling mutual information with network science, we explore how information flows through feedforward networks. We show that efficiently approximating mutual information allows us to create an information measure that quantifies how much information flows between any two neurons of a deep learning model. To that end, we propose NIF, Neural Information Flow, a technique for codifying information flow that exposes deep learning model internals and provides feature attributions.


title: Quantifying Interpretability and Trust in Machine Learning Systems

abstract: Decisions by Machine Learning (ML) models have become ubiquitous. Trusting these decisions requires understanding how algorithms take them. Hence interpretability methods for ML are an active focus of research. A central problem in this context is that both the quality of interpretability methods as well as trust in ML predictions are difficult to measure. Yet evaluations, comparisons and improvements of trust and interpretability require quantifiable measures. Here we propose a quantitative measure for the quality of interpretability methods. Based on that we derive a quantitative measure of trust in ML decisions. Building on previous work we propose to measure intuitive understanding of algorithmic decisions using the information transfer rate at which humans replicate ML model predictions. We provide empirical evidence from crowdsourcing experiments that the proposed metric robustly differentiates interpretability methods. The proposed metric also demonstrates the value of interpretability for ML assisted human decision making: in our experiments providing explanations more than doubled productivity in annotation tasks. However unbiased human judgement is critical for doctors, judges, policy makers and others. Here we derive a trust metric that identifies when human decisions are overly biased towards ML predictions. Our results complement existing qualitative work on trust and interpretability by quantifiable measures that can serve as objectives for further improving methods in this field of research.


title: NIHAO XXI: The emergence of Low Surface Brightness galaxies

abstract: The existence of galaxies with a surface brightness $\mu$ lower than the night sky has been known since three decades. Yet, their formation mechanism and emergence within a $\rm\Lambda CDM$ universe has remained largely undetermined. For the first time, we investigated the origin of Low Surface Brightness (LSB) galaxies with M$_{\star}$$\sim$10$^{9.5-10}$M$_{\odot}$, which we are able to reproduce within hydrodynamical cosmological simulations from the NIHAO suite. The simulated and observed LSBs share similar properties, having large HI reservoir, extended star formation histories and effective radii, low S\'{e}rsic index and slowly rising rotation curves. The formation mechanism of these objects is explored: simulated LSBs form as a result of co-planar co-rotating mergers and aligned accretion of gas at early times, while perpendicular mergers and mis-aligned gas accretion result in higher $\mu$ galaxies by $z$=0. The larger the merger, the stronger the correlation between merger orbital configuration and final $\mu$. While the halo spin parameter is consistently high in simulated LSB galaxies, the impact of halo concentration, feedback-driven gas outflows and merger time only plays a minor-to-no role in determining $\mu$. Interestingly, the formation scenario of such `classical' LSBs differs from the one of less massive, M$_{\star}$$\sim$10$^{7-9}$M$_{\odot}$, Ultra-Diffuse Galaxies, the latter resulting from the effects of SNae driven gas outflows: a M$_{\star}$ of $\sim$10$^9$M$_{\odot}$ thus represents the transition regime between a feedback-dominated to an angular momentum-dominated formation scenario in the LSB realm. Observational predictions are offered regarding spatially resolved star formation rates through LSB discs: these, together with upcoming surveys, can be used to verify the proposed emergence scenario of LSB galaxies.


title: Semi-Unsupervised Learning: Clustering and Classifying using Ultra-Sparse Labels

abstract: In semi-supervised learning for classification, it is assumed that every ground truth class of data is present in the small labelled dataset. Many real-world sparsely-labelled datasets are plausibly not of this type. It could easily be the case that some classes of data are found only in the unlabelled dataset -- perhaps the labelling process was biased -- so we do not have any labelled examples to train on for some classes. We call this learning regime $\textit{semi-unsupervised learning}$, an extreme case of semi-supervised learning, where some classes have no labelled exemplars in the training set. First, we outline the pitfalls associated with trying to apply deep generative model (DGM)-based semi-supervised learning algorithms to datasets of this type. We then show how a combination of clustering and semi-supervised learning, using DGMs, can be brought to bear on this problem. We study several different datasets, showing how one can still learn effectively when half of the ground truth classes are entirely unlabelled and the other half are sparsely labelled.


title: Rapid falling of an orbiting moon to its parent planet due to tidal-seismic resonance

abstract: Tidal force plays an important role in the evolution of the planet-moon system. The tidal force of a moon can excite seismic waves in the planet it is orbiting. A tidal-seismic resonance is expected when a tidal force frequency matches a free-oscillation frequency of the planet. Here we show that when the moon is close to the planet, the tidal-seismic resonance can cause large-amplitude seismic waves, which can change the shape of the planet and in turn exert a negative torque on the moon to cause it to fall rapidly toward the planet. We postulate that the tidal-seismic resonance may be an important mechanism which can accelerate planet accretion process. On the other hand, tidal-seismic resonance effect can also be used to interrogate planet interior by long term tracking of the orbital change of the moon.


title: Sample Complexity of Estimating the Policy Gradient for Nearly Deterministic Dynamical Systems

abstract: Reinforcement learning is a promising approach to learning robotics controllers. It has recently been shown that algorithms based on finite-difference estimates of the policy gradient are competitive with algorithms based on the policy gradient theorem. We propose a theoretical framework for understanding this phenomenon. Our key insight is that many dynamical systems (especially those of interest in robotics control tasks) are nearly deterministic -- i.e., they can be modeled as a deterministic system with a small stochastic perturbation. We show that for such systems, finite-difference estimates of the policy gradient can have substantially lower variance than estimates based on the policy gradient theorem. Finally, we empirically evaluate our insights in an experiment on the inverted pendulum.


title: The two-loop five-point amplitude in $\mathcal N=8$ supergravity

abstract: We compute the symbol of the two-loop five-point amplitude in $\mathcal N=8$ supergravity. We write an ansatz for the amplitude whose rational prefactors are based on not only 4-dimensional leading singularities, but also $d$-dimensional ones, as the former are insufficient. Our novel $d$-dimensional unitarity-based approach to the systematic construction of an amplitude's rational structures is likely to have broader applications, for example to analogous QCD calculations. We fix parameters in the ansatz by performing numerical integration-by-parts reduction of the known integrand. We find that the two-loop five-point $\mathcal N=8$ supergravity amplitude is uniformly transcendental. We then verify the soft and collinear limits of the amplitude. There is considerable similarity with the corresponding amplitude for $\mathcal N=4$ super-Yang-Mills theory: all the rational prefactors are double copies of the Yang-Mills ones and the transcendental functions overlap to a large degree. As a byproduct, we find new relations between color-ordered loop amplitudes in $\mathcal N=4$ super-Yang-Mills theory.


title: Infinite All-Layers Simple Foldability

abstract: We study the problem of deciding whether a crease pattern can be folded by simple folds (folding along one line at a time) under the infinite all-layers model introduced by [Akitaya et al., 2017], in which each simple fold is defined by an infinite line and must fold all layers of paper that intersect this line. This model is motivated by folding in manufacturing such as sheet-metal bending. We improve on [Arkin et al., 2004] by giving a deterministic $O(n)$-time algorithm to decide simple foldability of 1D crease patterns in the all-layers model. Then we extend this 1D result to 2D, showing that simple foldability in this model can be decided in linear time for unassigned axis-aligned orthogonal crease patterns on axis-aligned 2D orthogonal paper. On the other hand, we show that simple foldability is strongly NP-complete if a subset of the creases have a mountain-valley assignment, even for an axis-aligned rectangle of paper.


title: Learning Neurosymbolic Generative Models via Program Synthesis

abstract: Significant strides have been made toward designing better generative models in recent years. Despite this progress, however, state-of-the-art approaches are still largely unable to capture complex global structure in data. For example, images of buildings typically contain spatial patterns such as windows repeating at regular intervals; state-of-the-art generative methods can't easily reproduce these structures. We propose to address this problem by incorporating programs representing global structure into the generative model---e.g., a 2D for-loop may represent a configuration of windows. Furthermore, we propose a framework for learning these models by leveraging program synthesis to generate training data. On both synthetic and real-world data, we demonstrate that our approach is substantially better than the state-of-the-art at both generating and completing images that contain global structure.


title: Operational relevance of resource theories of quantum measurements

abstract: For any resource theory it is essential to identify tasks for which resource objects offer advantage over free objects. We show that this identification can always be accomplished for resource theories of quantum measurements in which free objects form a convex subset of measurements on a given Hilbert space. To this aim we prove that every resource measurement offers advantage for some quantum state discrimination task. Moreover, we give an operational interpretation of robustness, which quantifies the minimal amount of noise that must be added to a measurement to make it free. Specifically, we show that this geometric quantity is related to the maximal relative advantage that a resource measurement offers in a class of minimal-error state discrimination problems. Finally, we apply our results to two classes of free measurements: incoherent measurements (measurements that are diagonal in the fixed basis) and separable measurements (measurements whose effects are separable operators). For both of these scenarios we find, in the asymptotic setting in which the dimension or the number of particles increase to infinity, the maximal relative advantage that resource measurements offer for state discrimination tasks.


title: F1/10: An Open-Source Autonomous Cyber-Physical Platform

abstract: In 2005 DARPA labeled the realization of viable autonomous vehicles (AVs) a grand challenge; a short time later the idea became a moonshot that could change the automotive industry. Today, the question of safety stands between reality and solved. Given the right platform the CPS community is poised to offer unique insights. However, testing the limits of safety and performance on real vehicles is costly and hazardous. The use of such vehicles is also outside the reach of most researchers and students. In this paper, we present F1/10: an open-source, affordable, and high-performance 1/10 scale autonomous vehicle testbed. The F1/10 testbed carries a full suite of sensors, perception, planning, control, and networking software stacks that are similar to full scale solutions. We demonstrate key examples of the research enabled by the F1/10 testbed, and how the platform can be used to augment research and education in autonomous systems, making autonomy more accessible.


title: Algorithms for Fairness in Sequential Decision Making

abstract: It has recently been shown that if feedback effects of decisions are ignored, then imposing fairness constraints such as demographic parity or equality of opportunity can actually exacerbate unfairness. We propose to address this challenge by modeling feedback effects as Markov decision processes (MDPs). First, we propose analogs of fairness properties for the MDP setting. Second, we propose algorithms for learning fair decision-making policies for MDPs. Finally, we demonstrate the need to account for dynamical effects using simulations on a loan applicant MDP.


title: Magnetic Moment of Leptons

abstract: We show that the magnetic moment of leptons is significantly modified in thermal background as compared to the corresponding vacuum value. We compare the magnetic moment of all different leptons near nucleosynthesis. It is shown that the significance of thermal corrections depends on the temperature of the universe and the respective lepton mass. In the early universe, particle mass was growing quadratically with temperature which affects the corresponding value of magnetic moment. Intrinsic magnetic moment is inversely proportional to mass whereas thermal corrections to the neutrino dipole moment is proportional to neutrino mass for the Dirac type neutrino in the minimal standard model. Therefore the effect of temperature is not the same for charged leptons and neutrinos.


title: End-to-End Optimized Transmission over Dispersive Intensity-Modulated Channels Using Bidirectional Recurrent Neural Networks

abstract: We propose an autoencoding sequence-based transceiver for communication over dispersive channels with intensity modulation and direct detection (IM/DD), designed as a bidirectional deep recurrent neural network (BRNN). The receiver uses a sliding window technique to allow for efficient data stream estimation. We find that this sliding window BRNN (SBRNN), based on end-to-end deep learning of the communication system, achieves a significant bit-error-rate reduction at all examined distances in comparison to previous block-based autoencoders implemented as feed-forward neural networks (FFNNs), leading to an increase of the transmission distance. We also compare the end-to-end SBRNN with a state-of-the-art IM/DD solution based on two level pulse amplitude modulation with an FFNN receiver, simultaneously processing multiple received symbols and approximating nonlinear Volterra equalization. Our results show that the SBRNN outperforms such systems at both 42 and 84\,Gb/s, while training fewer parameters. Our novel SBRNN design aims at tailoring the end-to-end deep learning-based systems for communication over nonlinear channels with memory, such as the optical IM/DD fiber channel.


title: Nonparametric Inference under B-bits Quantization

abstract: Statistical inference based on lossy or incomplete samples is often needed in research areas such as signal/image processing, medical image storage, remote sensing, signal transmission. In this paper, we propose a nonparametric testing procedure based on samples quantized to $B$ bits through a computationally efficient algorithm. Under mild technical conditions, we establish the asymptotic properties of the proposed test statistic and investigate how the testing power changes as $B$ increases. In particular, we show that if $B$ exceeds a certain threshold, the proposed nonparametric testing procedure achieves the classical minimax rate of testing (Shang and Cheng, 2015) for spline models. We further extend our theoretical investigations to a nonparametric linearity test and an adaptive nonparametric test, expanding the applicability of the proposed methods. Extensive simulation studies {together with a real-data analysis} are used to demonstrate the validity and effectiveness of the proposed tests.


title: Width Provably Matters in Optimization for Deep Linear Neural Networks

abstract: We prove that for an $L$-layer fully-connected linear neural network, if the width of every hidden layer is $\tilde\Omega (L \cdot r \cdot d_{\mathrm{out}} \cdot \kappa^3 )$, where $r$ and $\kappa$ are the rank and the condition number of the input data, and $d_{\mathrm{out}}$ is the output dimension, then gradient descent with Gaussian random initialization converges to a global minimum at a linear rate. The number of iterations to find an $\epsilon$-suboptimal solution is $O(\kappa \log(\frac{1}{\epsilon}))$. Our polynomial upper bound on the total running time for wide deep linear networks and the $\exp\left(\Omega\left(L\right)\right)$ lower bound for narrow deep linear neural networks [Shamir, 2018] together demonstrate that wide layers are necessary for optimizing deep models.


title: Theoretically Principled Trade-off between Robustness and Accuracy

abstract: We identify a trade-off between robustness and accuracy that serves as a guiding principle in the design of defenses against adversarial examples. Although this problem has been widely studied empirically, much remains unknown concerning the theory underlying this trade-off. In this work, we decompose the prediction error for adversarial examples (robust error) as the sum of the natural (classification) error and boundary error, and provide a differentiable upper bound using the theory of classification-calibrated loss, which is shown to be the tightest possible upper bound uniform over all probability distributions and measurable predictors. Inspired by our theoretical analysis, we also design a new defense method, TRADES, to trade adversarial robustness off against accuracy. Our proposed algorithm performs well experimentally in real-world datasets. The methodology is the foundation of our entry to the NeurIPS 2018 Adversarial Vision Challenge in which we won the 1st place out of ~2,000 submissions, surpassing the runner-up approach by $11.41\%$ in terms of mean $\ell_2$ perturbation distance.


title: Three candidate globular clusters discovered in the Galactic bulge

abstract: This work reports the discovery of three new globular clusters (GCs) towards the Galactic bulge - Camargo 1107, 1108, and 1109. The discovery was made using the WISE, 2MASS, VVV, and Gaia-DR2 photometry. The new findings are old (12.0-13.5 Gyr) and metal-poor GCs ([Fe/H] < -1.5 dex) located in the bulge area close to the Milky Way (MW) mid-plane. Although the old ages and low metallicities suggest that the newly discovered GCs are likely associated with the inner halo the possibility of these clusters being part of a primordial bulge GC subpopulation cannot be ruled out. Camargo 1107, for instance, presents a metallicity of [Fe/H] = -2.2 +/- 0.4 dex and an age of 13.5 +/- 2 Gyr, which may suggest that this cluster formed just after the Big Bang in the very early Universe. The discovery of GCs such as the new findings is crucial to built a coherent picture of the inner Galaxy. It is likely that at least a few more dozens of GCs are still to be discovered in the bulge.


title: Deterministic 2-Dimensional Temperature-1 Tile Assembly Systems Cannot Compute

abstract: We consider non cooperative binding in so called `temperature 1', in deterministic (here called {\it confluent}) tile self-assembly systems (1-TAS) and prove the standing conjecture that such systems do not have universal computational power. We call a TAS whose maximal assemblies contain at least one ultimately periodic assembly path {\it para-periodic}. We observe that a confluent 1-TAS has at most one maximal producible assembly, $\alpha_{max}$, that can be considered a union of path assemblies, and we show that such a system is always para-periodic. This result is obtained through a superposition and a combination of two paths that produce a new path with desired properties, a technique that we call \emph{co-grow} of two paths. Moreover we provide a characterization of an $\alpha_{max}$ of a confluent 1-TAS as one of two possible cases, so called, a grid or a disjoint union of combs. To a given $\alpha_{max}$ we can associate a finite labeled graph, called \emph{quipu}, such that the union of all labels of paths in the quipu equals $\alpha_{max}$, therefore giving a finite description for $\alpha_{max}$. This finite description implies that $\alpha_{max}$ is a union of semi-affine subsets of $\mathbb{Z}^2$ and since such a finite description can be algorithmicly generated from any 1-TAS, 1-TAS cannot have universal computational power.


title: Learning Interpretable Models with Causal Guarantees

abstract: Machine learning has shown much promise in helping improve the quality of medical, legal, and financial decision-making. In these applications, machine learning models must satisfy two important criteria: (i) they must be causal, since the goal is typically to predict individual treatment effects, and (ii) they must be interpretable, so that human decision makers can validate and trust the model predictions. There has recently been much progress along each direction independently, yet the state-of-the-art approaches are fundamentally incompatible. We propose a framework for learning interpretable models from observational data that can be used to predict individual treatment effects (ITEs). In particular, our framework converts any supervised learning algorithm into an algorithm for estimating ITEs. Furthermore, we prove an error bound on the treatment effects predicted by our model. Finally, in an experiment on real-world data, we show that the models trained using our framework significantly outperform a number of baselines.


title: Using Erasure Feedback for Online Timely Updating with an Energy Harvesting Sensor

abstract: A real-time status updating system is considered, in which an energy harvesting sensor is acquiring measurements regarding some physical phenomenon and sending them to a destination through an erasure channel. The setting is online, in which energy arrives in units according to a Poisson process with unit rate, with arrival times being revealed causally over time. Energy is saved in a unit-sized battery. The sensor is notified by the destination of whether updates were erased via feedback. Updates need to reach the destination successfully in a timely fashion, namely, such that the long term average age of information, defined as the time elapsed since the latest successful update has reached the destination, is minimized. First, it is shown that the optimal status update policy has a renewal structure: successful update times should constitute a renewal process. Then, threshold-greedy policies are investigated: a new update is transmitted, following a successful one, only if the age of information grows above a certain threshold; and if it is erased, then all subsequent update attempts are greedily scheduled whenever energy is available. The optimal threshold-greedy policy is then analytically derived.


title: Entropic repulsion for the occupation-time field of random interlacements conditioned on disconnection

abstract: We investigate percolation of the vacant set of random interlacements on $\mathbb{Z}^d$, $d\geq 3$, in the strongly percolative regime. We consider the event that the interlacement set at level $u$ disconnects the discrete blow-up of a compact set $A\subseteq \mathbb{R}^d$ from the boundary of an enclosing box. We derive asymptotic large deviation upper bounds on the probability that the local averages of the occupation times deviate from a specific function depending on the harmonic potential of $A$, when disconnection occurs. If certain critical levels coincide, which is plausible but open at the moment, these bounds imply that conditionally on disconnection, the occupation-time profile undergoes an entropic push governed by a specific function depending on $A$. Similar entropic repulsion phenomena conditioned on disconnection by level-sets of the discrete Gaussian free field on $\mathbb{Z}^d$, $d \geq 3$, have been obtained by the authors in arxiv:1808.09947. Our proofs rely crucially on the `solidification estimates' developed in arXiv:1706.07229 by A.-S. Sznitman and the second author.


title: Forecasting Transformative AI: An Expert Survey

abstract: Transformative AI technologies have the potential to reshape critical aspects of society in the near future. However, in order to properly prepare policy initiatives for the arrival of such technologies accurate forecasts and timelines are necessary. A survey was administered to attendees of three AI conferences during the summer of 2018 (ICML, IJCAI and the HLAI conference). The survey included questions for estimating AI capabilities over the next decade, questions for forecasting five scenarios of transformative AI and questions concerning the impact of computational resources in AI research. Respondents indicated a median of 21.5% of human tasks (i.e., all tasks that humans are currently paid to do) can be feasibly automated now, and that this figure would rise to 40% in 5 years and 60% in 10 years. Median forecasts indicated a 50% probability of AI systems being capable of automating 90% of current human tasks in 25 years and 99% of current human tasks in 50 years. The conference of attendance was found to have a statistically significant impact on all forecasts, with attendees of HLAI providing more optimistic timelines with less uncertainty. These findings suggest that AI experts expect major advances in AI technology to continue over the next decade to a degree that will likely have profound transformative impacts on society.


title: The GstLAL Search Analysis Methods for Compact Binary Mergers in Advanced LIGO's Second and Advanced Virgo's First Observing Runs

abstract: After their successful first observing run (September 12, 2015 - January 12, 2016), the Advanced LIGO detectors were upgraded to increase their sensitivity for the second observing run (November 30, 2016 - August 26, 2017). The Advanced Virgo detector joined the second observing run on August 1, 2017. We discuss the updates that happened during this period in the GstLAL-based inspiral pipeline, which is used to detect gravitational waves from the coalescence of compact binaries both in low latency and an offline configuration. These updates include deployment of a zero-latency whitening filter to reduce the over-all latency of the pipeline by up to 32 seconds, incorporation of the Virgo data stream in the analysis, introduction of a single-detector search to analyze data from the periods when only one of the detectors is running, addition of new parameters to the likelihood ratio ranking statistic, increase in the parameter space of the search, and introduction of a template mass-dependent glitch-excision thresholding method.


title: Scheme dependence of asymptotically free solutions

abstract: Recent studies have provided evidence for the existence of new asymptotically free trajectories in non-Abelian particle models without asymptotic symmetry in the high-energy limit. We extend these results to a general ${\rm SU}(N_{\rm L})\times {\rm SU}(N_{\rm c})$ Higgs-Yukawa model that includes the non-Abelian sector of the standard model, finding further confirmation for such scenarios for a wide class of regularizations that account for threshold behavior persisting to highest energies. We construct these asymptotically free trajectories within conventional $\overline{\text{MS}}$ schemes and systematic weak coupling expansions. The existence of these solutions is argued to be a scheme-independent phenomenon, as demonstrated for mass-dependent schemes based on general momentum-space infrared regularizations. A change of scheme induces a map of the theory's coupling space onto itself, which in the present case also translates into a reparametrization of the space of asymptotically free solutions.


title: Exploring the shell model of high-velocity features of Type Ia supernovae using TARDIS

abstract: We explore the possible nature of high-velocity features in Type~Ia supernovae by presenting synthetic spectra generated from hydrodynamic models of interaction between the supernova and a compact circumstellar shell. We use TARDIS to compute the spectra and compare model spectra to data from SN~2011fe at 2, 5, and 9 days after the explosion. We apply abundance models to the shells that consist of either a hydrogen, helium, or carbon-oxygen substrate with overall solar abundance of metals and depletion or enhancement of calcium abundance. We vary the calcium abundances in the shell and the ejecta to study the effect on the photospheric and high-velocity components of the calcium near-infrared triplet. The substrate leaves no imprint upon the spectra, but helium substrates are discouraged by demanding excessive calcium abundances. We find that we can approximately reproduce the blended high-velocity and photospheric velocity features at 2 and 5 days. The 9 day spectrum shows a distinct high-velocity component of the line. We are unable to reproduce this specific feature with the current models. We also explore helium-shell double detonation models, finding they tend to give no high-velocity features or excessively strong features. A very carefully chosen helium-shell mass might give a reasonable reproduction of the observed calcium features. We discuss the implications of our models for the nature of high-velocity features and their significance to the evolution and explosion of SN~Ia.


title: Memory-free dynamics for the TAP equations of Ising models with arbitrary rotation invariant ensembles of random coupling matrices

abstract: We propose an iterative algorithm for solving the Thouless-Anderson-Palmer (TAP) equations of Ising models with arbitrary rotation invariant (random) coupling matrices. In the thermodynamic limit, we prove by means of the dynamical functional method that the proposed algorithm converges when the so-called de Almeida Thouless (AT) criterion is fulfilled. Moreover, we give exact analytical expressions for the rate of the convergence.


title: Fine-Grained Analysis of Optimization and Generalization for Overparameterized Two-Layer Neural Networks

abstract: Recent works have cast some light on the mystery of why deep nets fit any data and generalize despite being very overparametrized. This paper analyzes training and generalization for a simple 2-layer ReLU net with random initialization, and provides the following improvements over recent works: (i) Using a tighter characterization of training speed than recent papers, an explanation for why training a neural net with random labels leads to slower training, as originally observed in [Zhang et al. ICLR'17]. (ii) Generalization bound independent of network size, using a data-dependent complexity measure. Our measure distinguishes clearly between random labels and true labels on MNIST and CIFAR, as shown by experiments. Moreover, recent papers require sample complexity to increase (slowly) with the size, while our sample complexity is completely independent of the network size. (iii) Learnability of a broad class of smooth functions by 2-layer ReLU nets trained via gradient descent. The key idea is to track dynamics of training and generalization via properties of a related kernel.


title: Graph heat mixture model learning

abstract: Graph inference methods have recently attracted a great interest from the scientific community, due to the large value they bring in data interpretation and analysis. However, most of the available state-of-the-art methods focus on scenarios where all available data can be explained through the same graph, or groups corresponding to each graph are known a priori. In this paper, we argue that this is not always realistic and we introduce a generative model for mixed signals following a heat diffusion process on multiple graphs. We propose an expectation-maximisation algorithm that can successfully separate signals into corresponding groups, and infer multiple graphs that govern their behaviour. We demonstrate the benefits of our method on both synthetic and real data.


title: Galaxy shape measurement synergies between LSST and Euclid

abstract: We demonstrate that a joint analysis of LSST-like ground-based imaging with Euclid-like space-based imaging leads to increased precision and accuracy in galaxy shape measurements. At galaxy magnitudes of $i \sim 24.5$, a combined survey analysis increases the effective galaxy number density for cosmic shear studies by $\sim 50$ percent in comparison to an analysis of each survey alone. Using a realistic distribution of galaxy sizes, ellipticities and magnitudes down to $i = 25.2$, we simulate LSST-like and Euclid-like images of over one million isolated galaxies. We compare the precision and accuracy of the recovered galaxy ellipticities for four different analyses: LSST-only, Euclid-only, a simultaneous joint-pixel analysis of the two surveys, and a simple catalogue-level survey combination. In the faint and small-galaxy regime, where neither survey excels alone, we find a $\sim 20$ percent increase in the precision of galaxy shape measurement when we adopt a joint-pixel analysis, compared to a catalogue-level combination. As the statistical power of cosmic shear is dominated by intrinsic ellipticity noise, however, this improvement in shape measurement noise only leads to a $\sim 5$ percent improvement in the effective number density of galaxies for lensing studies. We view this as the minimum improvement that should be expected from a joint-pixel analysis over a less accurate catalogue-level combination, as the former will also improve the capability of LSST to de-blend close objects.


title: Heavy ion jet physics studies using precision jet substructure and quark-gluon jet classification

abstract: In the first talk I discuss the usefulness of jet grooming for testing jet quenching mechanisms, and I present a calculation of soft-drop jet mass distribution in proton-proton and heavy ion collisions. In the second talk I discuss the comprehensive studies of quark and gluon jet substructure as a way to search for qualitatively new jet modification signatures. I emphasize the importance of precise understanding of jet substructure as well as powerful analysis tools in future heavy ion jet physics studies, with direct applications to hard probes using heavy flavors and hadronically-decaying electroweak bosons.


title: On the cosmological constant problem

abstract: An additional variation of the Einstein-Hilbert action with respect to the Planck mass provides a constraint on the average Ricci scalar that prevents vacuum energy from gravitating. Consideration of the evolution of the inhomogeneous matter distribution in the Universe with evaluation of the averaging constraint on disconnected matter cells that ultimately form isolated gravitationally bound structures yields a backreaction effect that self-consistently produces the cosmological constant of the background. A uniform prior on our location in the formation of these isolated structures implies a mean expectation for the present cosmological constant energy density parameter of $\Omega_{\Lambda}=0.704$, giving rise to a late-time acceleration of the cosmic expansion and a coincident current energy density of matter.


title: Organic fiducial inference

abstract: A substantial generalisation is put forward of the theory of subjective fiducial inference as it was outlined in earlier papers. In particular, this theory is extended to deal with cases where the data are discrete or categorical rather than continuous, and cases where there was important pre-data knowledge about some or all of the model parameters. The system for directly expressing and then handling this pre-data knowledge, which is via what are referred to as global and local pre-data functions for the parameters concerned, is distinct from that which involves attempting to directly represent this knowledge in the form of a prior distribution function over these parameters, and then using Bayes' theorem. In this regard, the individual attributes of what are identified as three separate types of fiducial argument, namely the strong, moderate and weak fiducial arguments, form an integral part of the theory that is developed. Various practical examples of the application of this theory are presented, including examples involving binomial, Poisson and multinomial data. The fiducial distribution functions for the parameters of the models in these examples are interpreted in terms of a generalised definition of subjective probability that was set out previously.


title: Zonal Jet Creation from Secondary Instability of Drift Waves for Plasma Edge Turbulence

abstract: A new strategy is presented to explain the creation and persistence of zonal flows widely observed in plasma edge turbulence. The core physics in the edge regime of the magnetic-fusion tokamaks can be described qualitatively by the one-state modified Hasegawa-Mima (MHM) model, which creates enhanced zonal flows and more physically relevant features in comparison with the familiar Charney-Hasegawa-Mima (CHM) model for both plasma and geophysical flows. The generation mechanism of zonal jets is displayed from the secondary instability analysis via nonlinear interactions with a background base state. Strong exponential growth in the zonal modes is induced due to a non-zonal drift wave base state in the MHM model, while stabilizing damping effect is shown with a zonal flow base state. Together with the selective decay effect from the dissipation, the secondary instability offers a complete characterization of the convergence process to the purely zonal structure. Direct numerical simulations with and without dissipation are carried out to confirm the instability theory. It shows clearly the emergence of a dominant zonal flow from pure non-zonal drift waves with small perturbation in the initial configuration. In comparison, the CHM model does not create instability in the zonal modes and usually converges to homogeneous turbulence.


title: The Freyd-Mitchell Embedding Theorem

abstract: Given a small abelian category $\mathcal{A}$, the Freyd-Mitchell embedding theorem states the existence of a ring $R$ and an exact full embedding $\mathcal{A} \rightarrow R$-Mod. This theorem is useful as it allows one to prove general results about abelian categories within the context of $R$-modules. The goal of this report is to flesh out the proof of the embedding theorem. We shall follow closely the material and approach presented in Freyd (1964). This means we will encounter such concepts as projective generators, injective cogenerators, the Yoneda embedding, injective envelopes, Grothendieck categories, subcategories of mono objects and subcategories of absolutely pure objects.