title: Neural Signatures of Motor Skill in the Resting Brain

abstract: Stroke-induced disturbances of large-scale cortical networks are known to be associated with the extent of motor deficits. We argue that identifying brain networks representative of motor behavior in the resting brain would provide significant insights for current neurorehabilitation approaches. Particularly, we aim to investigate the global configuration of brain rhythms and their relation to motor skill, instead of learning performance as broadly studied. We empirically approach this problem by conducting a three-dimensional physical space visuomotor learning experiment during electroencephalographic (EEG) data recordings with thirty-seven healthy participants. We demonstrate that across-subjects variations in average movement smoothness as the quantified measure of subjects' motor skills can be predicted from the global configuration of resting-state EEG alpha-rhythms (8-14 Hz) recorded prior to the experiment. Importantly, this neural signature of motor skill was found to be orthogonal to (independent of) task -- as well as to learning-related changes in alpha-rhythms, which we interpret as an organizing principle of the brain. We argue that disturbances of such configurations in the brain may contribute to motor deficits in stroke, and that reconfiguring stroke patients' brain rhythms by neurofeedback may enhance post-stroke neurorehabilitation.


title: Bolometric correction factors for Active Galactic Nuclei

abstract: The bolometric luminosity of active galactic nuclei (AGN) is difficult to determine and various approximations have been used to calibrate it against different observed properties. Here I combine theoretical calculations of optically thick, geometrically thin accretion disks, and observed X-ray properties of AGN, to provide new bolometric correction factors (\kbol) over a large range of black hole (BH) mass, accretion rate, and spin. This is particularly important in cases where the mass accretion rate cannot be determined from the observed spectral energy distribution, and in cases where luminosity-independent correction factors have been used. Simple powerlaw approximations of \kbol\ are provided for L(5100\AA), L(3000\AA), L(1400\AA), L(2-10 keV) and L(narrow \hb). In all cases the uncertainties are large mostly due to the unknown BH spin. Prior knowledge of the BH mass reduces the uncertainty considerably.


title: Association rule mining and itemset-correlation based variants

abstract: Association rules express implication formed relations among attributes in databases of itemsets. The apriori algorithm is presented, the basis for most association rule mining algorithms. It works by pruning away rules that need not be evaluated based on the user specified minimum support confidence. Additionally, variations of the algorithm are presented that enable it to handle quantitative attributes and to extract rules about generalizations of items, but preserve the downward closure property that enables pruning. Intertransformation of the extensions is proposed for special cases.


title: Analysis of Worst-Case Interference in Underlay Radar-Massive MIMO Spectrum Sharing Scenarios

abstract: In this paper, we consider an underlay radar-massive MIMO spectrum sharing scenario in which massive MIMO base stations (BSs) are allowed to operate outside a circular exclusion zone centered at the radar. Modeling the locations of the massive MIMO BSs as a homogeneous Poisson point process (PPP), we derive an analytical expression for a tight upper bound on the average interference at the radar due to cellular transmissions. The technical novelty is in bounding the worst-case elevation angle for each massive MIMO BS for which we devise a novel construction based on the circumradius distribution of a typical Poisson-Voronoi (PV) cell. While these worst-case elevation angles are correlated for neighboring BSs due to the structure of the PV tessellation, it does not explicitly appear in our analysis because of our focus on the average interference. We also provide an estimate of the nominal average interference by approximating each cell as a circle with area equal to the average area of the typical cell. Using these results, we demonstrate that the gap between the two results remains approximately constant with respect to the exclusion zone radius. Our analysis reveals useful trends in average interference power, as a function of key deployment parameters such as radar/BS antenna heights, number of antenna elements per radar/BS, BS density, and exclusion zone radius.


title: On phase retrieval via matrix completion and the estimation of low rank PSD matrices

abstract: Given underdetermined measurements of a Positive Semi-Definite (PSD) matrix $X$ of known low rank $K$, we present a new algorithm to estimate $X$ based on recent advances in non-convex optimization schemes. We apply this in particular to the phase retrieval problem for Fourier data, which can be formulated as a rank 1 PSD matrix recovery problem. Moreover, we provide theory for how oversampling affects the stability of the lifted inverse problem.


title: BEHRT: Transformer for Electronic Health Records

abstract: Today, despite decades of developments in medicine and the growing interest in precision healthcare, vast majority of diagnoses happen once patients begin to show noticeable signs of illness. Early indication and detection of diseases, however, can provide patients and carers with the chance of early intervention, better disease management, and efficient allocation of healthcare resources. The latest developments in machine learning (more specifically, deep learning) provides a great opportunity to address this unmet need. In this study, we introduce BEHRT: A deep neural sequence transduction model for EHR (electronic health records), capable of multitask prediction and disease trajectory mapping. When trained and evaluated on the data from nearly 1.6 million individuals, BEHRT shows a striking absolute improvement of 8.0-10.8%, in terms of Average Precision Score, compared to the existing state-of-the-art deep EHR models (in terms of average precision, when predicting for the onset of 301 conditions). In addition to its superior prediction power, BEHRT provides a personalised view of disease trajectories through its attention mechanism; its flexible architecture enables it to incorporate multiple heterogeneous concepts (e.g., diagnosis, medication, measurements, and more) to improve the accuracy of its predictions; and its (pre-)training results in disease and patient representations that can help us get a step closer to interpretable predictions.


title: Channel Normalization in Convolutional Neural Network avoids Vanishing Gradients

abstract: Normalization layers are widely used in deep neural networks to stabilize training. In this paper, we consider the training of convolutional neural networks with gradient descent on a single training example. This optimization problem arises in recent approaches for solving inverse problems such as the deep image prior or the deep decoder. We show that for this setup, channel normalization, which centers and normalizes each channel individually, avoids vanishing gradients, whereas, without normalization, gradients vanish which prevents efficient optimization. This effect prevails in deep single-channel linear convolutional networks, and we show that without channel normalization, gradient descent takes at least exponentially many steps to come close to an optimum. Contrary, with channel normalization, the gradients remain bounded, thus avoiding exploding gradients.


title: Adversarial Feature Learning in Brain Interfacing: An Experimental Study on Eliminating Drowsiness Effects

abstract: Across- and within-recording variabilities in electroencephalographic (EEG) activity is a major limitation in EEG-based brain-computer interfaces (BCIs). Specifically, gradual changes in fatigue and vigilance levels during long EEG recording durations and BCI system usage bring along significant fluctuations in BCI performances even when these systems are calibrated daily. We address this in an experimental offline study from EEG-based BCI speller usage data acquired for one hour duration. As the main part of our methodological approach, we propose the concept of adversarial invariant feature learning for BCIs as a regularization approach on recently expanding EEG deep learning architectures, to learn nuisance-invariant discriminative features. We empirically demonstrate the feasibility of adversarial feature learning on eliminating drowsiness effects from event related EEG activity features, by using temporal recording block ordering as the source of drowsiness variability.


title: Dynamics of Dusty Vortices I: Extensions and limitations of the terminal velocity approximation

abstract: Motivated by the stability of dust laden vortices, in this paper we study the terminal velocity approximation equations for a gas coupled to a pressureless dust fluid and present a numerical solver for the equations embedded in the FARGO3D hydrodynamics code. We show that for protoplanetary discs it is possible to use the baricenter velocity in the viscous stress tensor, making it trivial to simulate viscous dusty protoplanetary discs with this model. We also show that the terminal velocity model breaks down around shocks, becoming incompatible with the two fluid model it is derived from. Finally we produce a set of test cases for numerical schemes and demonstrate the performance of our code on these tests. Our implementation embedded in FARGO3D using an unconditionally stable explicit integrator is fast, and exhibits the desired second order spatial convergence for smooth problems.


title: Testing fundamental physics with photon frequency shift

abstract: We propose a high precision satellite experiment to further test Einstein's General Relativity and constrain extended theories of gravity. We consider the frequency shift of a photon radially exchanged between two observers located on Earth and on a satellite in circular orbit in the equatorial plane. In General Relativity there exists a peculiar satellite-distance at which the static contribution to the frequency shift vanishes since the effects induced by pure gravity and special relativity compensate, while it can be non-zero in modified gravities, like in models with screening mechanisms. As an experimental device placed on the satellite we choose a system of hydrogen atoms which can exhibit the $1$s spin-flip transition from the singlet (unaligned proton-electron spins) to the triplet (aligned proton-electron spins) state induced by the absorption of photons at $21.1$cm. The observation of an excited state would indicate that the frequency of the emitted and absorbed photon remains unchanged according to General Relativity. On the contrary, a non-zero frequency shift, as predicted in extended theories of gravity, would prevent the spin-flip transition and the hydrogen atoms from jumping into the excited state. Such a detection would signify a smoking-gun signature of new physics beyond special and general relativity.

categories: gr-qc astro-ph.EP astro-ph.SR hep-th quant-ph


title: Spatial sensitivity analysis for urban land use prediction with physics-constrained conditional generative adversarial networks

abstract: Accurately forecasting urban development and its environmental and climate impacts critically depends on realistic models of the spatial structure of the built environment, and of its dependence on key factors such as population and economic development. Scenario simulation and sensitivity analysis, i.e., predicting how changes in underlying factors at a given location affect urbanization outcomes at other locations, is currently not achievable at a large scale with traditional urban growth models, which are either too simplistic, or depend on detailed locally-collected socioeconomic data that is not available in most places. Here we develop a framework to estimate, purely from globally-available remote-sensing data and without parametric assumptions, the spatial sensitivity of the (\textit{static}) rate of change of urban sprawl to key macroeconomic development indicators. We formulate this spatial regression problem as an image-to-image translation task using conditional generative adversarial networks (GANs), where the gradients necessary for comparative static analysis are provided by the backpropagation algorithm used to train the model. This framework allows to naturally incorporate physical constraints, e.g., the inability to build over water bodies. To validate the spatial structure of model-generated built environment distributions, we use spatial statistics commonly used in urban form analysis. We apply our method to a novel dataset comprising of layers on the built environment, nightlighs measurements (a proxy for economic development and energy use), and population density for the world's most populous 15,000 cities.


title: Networking and processing in optical wireless

abstract: Optical wireless communication (OWC) is a promising technology that can provide high data rates while supporting multiple users. The Optical Wireless (OW) physical layer has been researched extensively, however less work was devoted to multiple access and how the OW front end is connected to the network. In this paper, an OWC system which employs a wavelength division multiple access (WDMA) scheme is studied, for the purpose of supporting multiple users. In addition, a cloud/fog architecture is proposed for the first time for OWC to provide processing capabilities. The cloud/fog-integrated architecture uses visible indoor light to create high data rate connections with potential mobile nodes. These optical wireless nodes are further clustered and used as fog mini servers to provide processing services through the optical wireless channel for other users. Additional fog processing units are located in the room, the building, the campus and at the metro level. Further processing capabilities are provided by remote cloud sites. A mixed-integer linear programming (MILP) model was developed and utilised to optimise resource allocation in the indoor OWC system. A second MILP model was developed to optimise the placement of processing tasks in the different fog and cloud nodes available. The optimisation of tasks placement in the cloud-/fog-integrated architecture was analysed using the MILP models. Multiple scenarios were considered where the mobile node locations were varied in the room and the amount of processing and data rate requested by each optical wireless node is varied. The results help identify the optimum colour and access point to use for communication for a given mobile node location and OWC system configuration, the optimum location to place processing and the impact of the network architecture. Areas for future work are identified.


title: Towards Higher Spectral Efficiency: Rate-2 Full-Diversity Complex Space-Time Block Codes

abstract: The upcoming 5G networks demand high-speed and high spectral-efficiency communications to keep up with the proliferating traffic demands. To this end, Massive multiple-input multiple-output (MIMO) techniques have gained significant traction owing to its ability to achieve these without increasing bandwidth or density of base stations. The preexisting space-time block code (STBC) designs cannot achieve a rate of more than 1 for more than two transmit antennas while preserving the orthogonality and full diversity conditions. In this paper, we present Jagannath codes - a novel complex modulation STBC, that achieves a very high rate of 2 for three and four transmit antennas. The presented designs achieve full diversity and overcome the previously achieved rates with the three and four antenna MIMO systems. We present a detailed account of the code construction of the proposed designs, orthogonality and full diversity analysis, transceiver model and conditional maximum likelihood (ML) decoding. In an effort to showcase the improvement achieved with the presented designs, we compare the rates and delays of some of the known STBCs with the proposed designs. The effective spectral efficiency and coding gain of the presented designs are compared to the Asymmetric Coordinate Interleaved design (ACIOD) and Jafarkhani code. We presented an effective spectral efficiency improvement by a factor of 2 with the proposed Jagannath codes. Owing to the full diversity of the presented designs, we demonstrate significant coding gains (6 dB and 12 dB) with the proposed designs.


title: Numerically exact counting statistics of energy current in the Kondo regime

abstract: We use the inchworm Quantum Monte Carlo method to investigate the full counting statistics of particle and energy currents in a strongly correlated quantum dot. Our method is used to extract the heat fluctuations and entropy production of a quantum thermoelectric device, as well as cumulants of the particle and energy currents. The energy--particle current cross correlations reveal information on the preparation of the system and the interplay of thermal and electric currents. We furthermore demonstrate the signature of a crossover from Coulomb blockade to Kondo physics in the energy current fluctuations, and show how the conventional master equation approach to full counting statistics systematically fails to capture this crossover.


title: Stochastic algorithms with geometric step decay converge linearly on sharp functions

abstract: Stochastic (sub)gradient methods require step size schedule tuning to perform well in practice. Classical tuning strategies decay the step size polynomially and lead to optimal sublinear rates on (strongly) convex problems. An alternative schedule, popular in nonconvex optimization, is called \emph{geometric step decay} and proceeds by halving the step size after every few epochs. In recent work, geometric step decay was shown to improve exponentially upon classical sublinear rates for the class of \emph{sharp} convex functions. In this work, we ask whether geometric step decay similarly improves stochastic algorithms for the class of sharp nonconvex problems. Such losses feature in modern statistical recovery problems and lead to a new challenge not present in the convex setting: the region of convergence is local, so one must bound the probability of escape. Our main result shows that for a large class of stochastic, sharp, nonsmooth, and nonconvex problems a geometric step decay schedule endows well-known algorithms with a local linear rate of convergence to global minimizers. This guarantee applies to the stochastic projected subgradient, proximal point, and prox-linear algorithms. As an application of our main result, we analyze two statistical recovery tasks---phase retrieval and blind deconvolution---and match the best known guarantees under Gaussian measurement models and establish new guarantees under heavy-tailed distributions.


title: On the Equivalence Between Abstract Dialectical Frameworks and Logic Programs

abstract: Abstract Dialectical Frameworks (ADFs) are argumentation frameworks where each node is associated with an acceptance condition. This allows us to model different types of dependencies as supports and attacks. Previous studies provided a translation from Normal Logic Programs (NLPs) to ADFs and proved the stable models semantics for a normal logic program has an equivalent semantics to that of the corresponding ADF. However, these studies failed in identifying a semantics for ADFs equivalent to a three-valued semantics (as partial stable models and well-founded models) for NLPs. In this work, we focus on a fragment of ADFs, called Attacking Dialectical Frameworks (ADF$^+$s), and provide a translation from NLPs to ADF$^+$s robust enough to guarantee the equivalence between partial stable models, well-founded models, regular models, stable models semantics for NLPs and respectively complete models, grounded models, preferred models, stable models for ADFs. In addition, we define a new semantics for ADF$^+$s, called L-stable, and show it is equivalent to the L-stable semantics for NLPs. This paper is under consideration for acceptance in TPLP.


title: High harmonic generation spectroscopy via orbital angular momentum

abstract: We present an experimental technique using orbital angular momentum (OAM) in a fundamental laser field to drive High Harmonic Generation (HHG). The mixing of beams with different OAM allows to generate two laser foci tightly spaced to study the phase and amplitude of HHG produced in diatomic nitrogen. Nitrogen is used as a well studied system to show the quality of OAM based HHG interferometry.


title: A note on collapsibility of acyclic 2-complexes

abstract: We present a Morse-theoretic characterization of collapsibility for 2-dimensional acyclic simplicial complexes by means of the values of normalized optimal combinatorial Morse functions.


title: Cell differentiation: what have we learned in 50 years?

abstract: I revisit two theories of cell differentiation in multicellular organisms published a half-century ago, Stuart Kauffman's global gene regulatory dynamics (GGRD) model and Roy Britten's and Eric Davidson's modular gene regulatory network (MGRN) model, in light of newer knowledge of mechanisms of gene regulation in the metazoans (animals). The two models continue to inform hypotheses and computational studies of differentiation of lineage-adjacent cell types. However, their shared notion (based on bacterial regulatory systems) of gene switches and networks built from them, have constrained progress in understanding the dynamics and evolution of differentiation. Recent work has described unique write-read-rewrite chromatin-based expression encoding in eukaryotes, as well metazoan-specific processes of gene activation and silencing in condensed-phase, enhancer-recruiting regulatory hubs, employing disordered proteins, including transcription factors, with context-dependent identities. These findings suggest an evolutionary scenario in which the origination of differentiation in animals, rather than depending exclusively on adaptive natural selection, emerged as a consequence of a type of multicellularity in which the novel metazoan gene regulatory apparatus was readily mobilized to amplify and exaggerate inherent cell functions of unicellular ancestors. The plausibility of this hypothesis is illustrated by the evolution of the developmental role of Grainyhead-like in the formation of epithelium.


title: Applications of the perturbation formula for Poisson processes to elementary and geometric probability

abstract: The binomial, the negative binomial, the Poisson, the compound Poisson and the Erlang distribution do all admit integral representations with respect to its (continuous) parameter. We use the Margulis-Russo type formulas for Bernoulli and Poisson processes to derive these representations in a unified way and to provide a probabilistic interpretation for the derivatives. By similar variational methods, we obtain apparently new integro-differential identities which the density of a strictly $\alpha$-stable multivariate density satisfies. Then, we extend Crofton's derivative formula known in integral geometry to the case of a Poisson process. Finally we use this extension to give a new probabilistic proof of a version of this formula for binomial point processes.


title: Coupling material and mechanical design processes via computer model calibration

abstract: Computer model calibration typically operates by choosing parameter values in a computer model so that the model output faithfully predicts reality. By using performance targets in place of observed data, we show that calibration techniques can be repurposed to wed engineering and material design, two processes that are traditionally carried out separately. This allows materials to be designed with specific engineering targets in mind while quantifying the associated sources of uncertainty. We demonstrate our proposed approach by "calibrating" material design settings to performance targets for a wind turbine blade.


title: Product of Orthogonal Spheres Parameterization for Disentangled Representation Learning

abstract: Learning representations that can disentangle explanatory attributes underlying the data improves interpretabilty as well as provides control on data generation. Various learning frameworks such as VAEs, GANs and auto-encoders have been used in the literature to learn such representations. Most often, the latent space is constrained to a partitioned representation or structured by a prior to impose disentangling. In this work, we advance the use of a latent representation based on a product space of Orthogonal Spheres PrOSe. The PrOSe model is motivated by the reasoning that latent-variables related to the physics of image-formation can under certain relaxed assumptions lead to spherical-spaces. Orthogonality between the spheres is motivated via physical independence models. Imposing the orthogonal-sphere constraint is much simpler than other complicated physical models, is fairly general and flexible, and extensible beyond the factors used to motivate its development. Under further relaxed assumptions of equal-sized latent blocks per factor, the constraint can be written down in closed form as an ortho-normality term in the loss function. We show that our approach improves the quality of disentanglement significantly. We find consistent improvement in disentanglement compared to several state-of-the-art approaches, across several benchmarks and metrics.


title: An analysis of IQ LINK$^{TM}$

abstract: This is a theoretical and computational strategy exploration of the visually attractive game IQ-Link. Not all games which are visually appealing are worthy of your time as a puzzlist. This analysis gives a would-be addict some idea of what type of game they are falling for.


title: Heavy quark dynamics in a hot magnetized QCD medium

abstract: The heavy quark drag and momentum diffusion have been investigated in a hot magnetized quark-gluon plasma, along the directions parallel and perpendicular to the magnetic field. The analysis is done within the framework of Fokker-Planck dynamics by considering the heavy quark scattering with thermal quarks and gluons at the leading order in the coupling constant. An extended quasiparticle model is adopted to encode the thermal QCD medium interactions in the presence of a magnetic field. Further, the higher Landau level effects on the temperature behaviour of the parallel and perpendicular components of the drag force and diffusion coefficients have studied. It has been observed that both the equation of state and the magnetic field play key roles in the temperature dependence of the heavy quark dynamics.


title: Relational Generalized Few-Shot Learning

abstract: Transferring learned models to novel tasks is a challenging problem, particularly if only very few labeled examples are available. Although this few-shot learning setup has received a lot of attention recently, most proposed methods focus on discriminating novel classes only. Instead, we consider the extended setup of generalized few-shot learning (GFSL), where the model is required to perform classification on the joint label space consisting of both previously seen and novel classes. We propose a graph-based framework that explicitly models relationships between all seen and novel classes in the joint label space. Our model Graph-convolutional Global Prototypical Networks (GcGPN) incorporates these inter-class relations using graph-convolution in order to embed novel class representations into the existing space of previously seen classes in a globally consistent manner. Our approach ensures both fast adaptation and global discrimination, which is the major challenge in GFSL. We demonstrate the benefits of our model on two challenging benchmark datasets.


title: System-Level Development of a User-Integrated Semi-Autonomous Lawn Mowing System: Problem Overview, Basic Requirements, and Proposed Architecture

abstract: This concept paper outlines some recent efforts toward the design and development of user-integrated semi-autonomous home-sized lawn mowing systems from a systems engineering perspective. This is an important and emerging field of study within the robotics and systems engineering communities. The work presented includes a review of current progress on this problem, a discussion of the problem from a systems engineering perspective, a general system architecture developed by the authors, and a preliminary set of design requirements. This work is meant to provide a baseline and motivation for the further development and refinement of these systems within the systems engineering and robotics communities and is relevant to both academic and commercial research.


title: Beyond NP: Quantifying over Answer Sets

abstract: Answer Set Programming (ASP) is a logic programming paradigm featuring a purely declarative language with comparatively high modeling capabilities. Indeed, ASP can model problems in NP in a compact and elegant way. However, modeling problems beyond NP with ASP is known to be complicated, on the one hand, and limited to problems in {\Sigma}^P_2 on the other. Inspired by the way Quantified Boolean Formulas extend SAT formulas to model problems beyond NP, we propose an extension of ASP that introduces quantifiers over stable models of programs. We name the new language ASP with Quantifiers (ASP(Q)). In the paper we identify computational properties of ASP(Q); we highlight its modeling capabilities by reporting natural encodings of several complex problems with applications in artificial intelligence and number theory; and we compare ASP(Q) with related languages. Arguably, ASP(Q) allows one to model problems in the Polynomial Hierarchy in a direct way, providing an elegant expansion of ASP beyond the class NP. Under consideration for acceptance in TPLP.


title: Better Paracoherent Answer Sets with Less Resources

abstract: Answer Set Programming (ASP) is a well-established formalism for logic programming. Problem solving in ASP requires to write an ASP program whose answers sets correspond to solutions. Albeit the non-existence of answer sets for some ASP programs can be considered as a modeling feature, it turns out to be a weakness in many other cases, and especially for query answering. Paracoherent answer set semantics extend the classical semantics of ASP to draw meaningful conclusions also from incoherent programs, with the result of increasing the range of applications of ASP. State of the art implementations of paracoherent ASP adopt the semi-equilibrium semantics, but cannot be lifted straightforwardly to compute efficiently the (better) split semi-equilibrium semantics that discards undesirable semi-equilibrium models. In this paper an efficient evaluation technique for computing a split semi-equilibrium model is presented. An experiment on hard benchmarks shows that better paracoherent answer sets can be computed consuming less computational resources than existing methods. Under consideration for acceptance in TPLP.


title: The Future of Exoplanet Direct Detection

abstract: Diffraction fundamentally limits our ability to image and characterize exoplanets. Current and planned coronagraphic searches for exoplanets are making incredible strides but are fundamentally limited by the inner working angle of a few lambda/D. Some crucial topics, such as demographics of exoplanets within the first 50 Myr and the infrared characterization of terrestrial planets, are beyond the reach of the single aperture angular resolution for the foreseeable future. Interferometry offers some advantages in exoplanet detection and characterization and we explore in this white paper some of the potential scientific breakthroughs possible. We demonstrate here that investments in 'exoplanet interferometry' could open up new possibilities for speckle suppression through spatial coherence, a giant boost in astrometric precision for determining exoplanet orbits, ability to take a census of young giant exoplanets (clusters <50 Myr age), and an unrivaled potential for infrared nulling from space to detect terrestrial planets and search for atmospheric biomarkers. All signs point to an exciting future for exoplanets and interferometers, albeit a promise that will take decades to fulfill.


title: Practical Newton-Type Distributed Learning using Gradient Based Approximations

abstract: We study distributed algorithms for expected loss minimization where the datasets are large and have to be stored on different machines. Often we deal with minimizing the average of a set of convex functions where each function is the empirical risk of the corresponding part of the data. In the distributed setting where the individual data instances can be accessed only on the local machines, there would be a series of rounds of local computations followed by some communication among the machines. Since the cost of the communication is usually higher than the local machine computations, it is important to reduce it as much as possible. However, we should not allow this to make the computation too expensive to become a burden in practice. Using second-order methods could make the algorithms converge faster and decrease the amount of communication needed. There are some successful attempts in developing distributed second-order methods. Although these methods have shown fast convergence, their local computation is expensive and could enjoy more improvement for practical uses. In this study we modify an existing approach, DANE (Distributed Approximate NEwton), in order to improve the computational cost while maintaining the accuracy. We tackle this problem by using iterative methods for solving the local subproblems approximately instead of providing exact solutions for each round of communication. We study how using different iterative methods affect the behavior of the algorithm and try to provide an appropriate tradeoff between the amount of local computation and the required amount of communication. We demonstrate the practicality of our algorithm and compare it to the existing distributed gradient based methods such as SGD.


title: Minimization of visibly pushdown automata is NP-complete

abstract: We show that the minimization of visibly pushdown automata is NP-complete. This result is obtained by introducing immersions, that recognize multiple languages (over a usual, non-visible alphabet) using a common deterministic transition graph, such that each language is associated with an initial state and a set of final states. We show that minimizing immersions is NP-complete, and reduce this problem to the minimization of visibly pushdown automata.


title: Imaging the Key Stages of Planet Formation

abstract: New images of young stars are revolutionizing our understanding of planet formation. ALMA detects large grains in planet-forming disks with few AU scale resolution and scattered light imaging with extreme adaptive optics systems reveal small grains suspended on the disks' flared surfaces. Tantalizing evidence for young exoplanets is emerging from line observations of CO and H-alpha. In this white paper, we explore how even higher angular resolution can extend our understanding of the key stages of planet formation, to resolve accreting circumplanetary disks themselves, and to watch planets forming in situ for the nearest star-forming regions. We focus on infrared observations which are sensitive to thermal emission in the terrestrial planet formation zone and allow access to molecular tracers in warm ro-vibrational states. Successful planet formation theories will not only be able to explain the diverse features seen in disks, but will also be consistent with the rich exoplanet demographics from RV and transit surveys. While we are far from exhausting ground-based techniques, the ultimate combination of high angular resolution and high infrared sensitivity can only be achieved through mid-infrared space interferometry.


title: Classification with the matrix-variate-$t$ distribution

abstract: Matrix-variate distributions can intuitively model the dependence structure of matrix-valued observations that arise in applications with multivariate time series, spatio-temporal or repeated measures. This paper develops an Expectation-Maximization algorithm for discriminant analysis and classification with matrix-variate $t$-distributions. The methodology shows promise on simulated datasets or when applied to the forensic matching of fractured surfaces or the classification of functional Magnetic Resonance, satellite or hand gestures images.


title: Molecules as metric measure spaces with Kato-bounded Ricci curvature

abstract: Set $\Psi:=-\log(\tilde{\Psi})$, with $\tilde{\Psi}>0$ the ground state of an arbitrary molecule with $n$ electrons in the infinite mass limit (neglecting spin/statistics). Let $\Sigma\subset \IR^{3n}$ be the set of singularities of the underlying Coulomb potential. We show that the metric measure space $\IMM$ given by $\IR^{3n}$ with its Euclidean distance and the measure $$ \mu(dx)=e^{-2\Psi(x)}dx $$ has a Bakry-Emery-Ricci tensor which is absolutely bounded by the the function $x\mapsto |x-\Sigma|^{-1}$, which we show to be an element of the Kato class induced by $\IMM$. In addition, it is shown $\IMM$ is stochastically complete, that is, the Brownian motion which is induced by a molecule is nonexplosive, and that the heat semigroup of $\IMM$ has the $L^{\infty}$-to-Lipschitz smoothing property. Our proofs reveal a fundamental connection between the above geometric/probabilistic properties and recently obtained derivative estimates for $e^{\Psi}$ by Fournais/S\o rensen, as well as Aizenman/Simon's Harnack inequality for Schr\"odinger operators. Moreover, our results suggest to study general metric measure spaces having a Ricci curvature which is synthetically bounded from below/above by a function in the underlying Kato class.


title: The Effect of Visual Design in Image Classification

abstract: Financial companies continuously analyze the state of the markets to rethink and adjust their investment strategies. While the analysis is done on the digital form of data, decisions are often made based on graphical representations in white papers or presentation slides. In this study, we examine whether binary decisions are better to be decided based on the numeric or the visual representation of the same data. Using two data sets, a matrix of numerical data with spatial dependencies and financial data describing the state of the S&P index, we compare the results of supervised classification based on the original numerical representation and the visual transformation of the same data. We show that, for these data sets, the visual transformation results in higher predictability skill compared to the original form of the data. We suggest thinking of the visual representation of numeric data, effectively, as a combination of dimensional reduction and feature engineering techniques. In particular, if the visual layout encapsulates the full complexity of the data. In this view, thoughtful visual design can guard against overfitting, or introduce new features -- all of which benefit the learning process, and effectively lead to better recognition of meaningful patterns.


title: Horizontal branch morphology: A new photometric parametrization

abstract: Theory and observations indicate that the distribution of stars along the horizontal branch of Galactic globular clusters mainly depends on the metal content. However, the existence of globular clusters with similar metal content and absolute age but different horizontal branch morphologies, suggests the presence of another parameter affecting the star distribution along the branch. To investigate the variation of the horizontal branch morphology in Galactic globular clusters, we define a new photometric horizontal branch morphology index, overcoming some of the limitations and degeneracies affecting similar indices available in the literature. We took advantage of a sample of 64 Galactic globular clusters, with both space-based imaging data and homogeneous ground-based photometric catalogues in five different bands ($U$,$B$,$V$, $R$, $I$). The new index, $\tau_{HB}$, is defined as the ratio between the areas subtended by the cumulative number distribution in magnitude ($I$) and in colour ($V-I$) of all stars along the horizontal branch. This new index shows a linear trend over the entire range in metallicity (-2.35 $\leq$ [Fe/H] $\leq$ -0.12) covered by our Galactic globular cluster sample. We found a linear relation between $\tau_{HB}$ and absolute cluster ages. We also found a quadratic anti-correlation with [Fe/H], becoming linear when we eliminate the age effect on $\tau_{HB}$ values. Moreover, we identified a subsample of eight clusters that are peculiar according to their $\tau_{HB}$ values. These clusters have bluer horizontal branch morphology when compared to typical ones of similar metallicity. These findings allow us to define them as the 'second parameter' clusters in the sample. A comparison with synthetic horizontal branch models suggests that they cannot be entirely explained with a spread in helium content.


title: MemNet: Memory-Efficiency Guided Neural Architecture Search with Augment-Trim learning

abstract: Recent studies on automatic neural architectures search have demonstrated significant performance, competitive to or even better than hand-crafted neural architectures. However, most of the existing network architecture tend to use residual, parallel structures and concatenation block between shallow and deep features to construct a large network. This requires large amounts of memory for storing both weights and feature maps. This is challenging for mobile and embedded devices since they may not have enough memory to perform inference with the designed large network model. To close this gap, we propose MemNet, an augment-trim learning-based neural network search framework that optimizes not only performance but also memory requirement. Specifically, it employs memory consumption based ranking score which forces an upper bound on memory consumption for navigating the search process. Experiment results show that, as compared to the state-of-the-art efficient designing methods, MemNet can find an architecture which can achieve competitive accuracy and save an average of 24.17% on the total memory needed.


title: Derivation of the two Schwarzians effective action for the Sachdev-Ye-Kitaev spectral form factor

abstract: The Sachdev-Ye-Kitaev model spectral form factor exhibits absence of information loss in the form of a ramp and a plateau, that are typical of random matrix theory. In a large $N$ collective fields description, the ramp was reproduced by Saad, Shenker and Stanford \cite{Saad:2018bqo}, by replica symmetry breaking saddles for a connected component of the analytically continued to real times thermal partition function two point function. We derive a two sides Schwarzians effective action for fluctuations around the ramp critical saddles, by adapting to the two replica system a method by Kitaev and Suh \cite{Kitaev:2017awl} for studying non linear responses to the conformal breaking kinetic operator in regular SYK. Our result confirms \cite{Saad:2018bqo}, where the form of the action was obtained by assuming locality.


title: Efficient Deep Learning Techniques for Multiphase Flow Simulation in Heterogeneous Porous Media

abstract: We present efficient deep learning techniques for approximating flow and transport equations for both single phase and two-phase flow problems. The proposed methods take advantages of the sparsity structures in the underlying discrete systems and can be served as efficient alternatives to the system solvers at the full order. In particular, for the flow problem, we design a network with convolutional and locally connected layers to perform model reductions. Moreover, we employ a custom loss function to impose local mass conservation constraints. This helps to preserve the physical property of velocity solution which we are interested in learning. For the saturation problem, we propose a residual type of network to approximate the dynamics. Our main contribution here is the design of custom sparsely connected layers which take into account the inherent sparse interaction between the input and output. After training, the approximated feed-forward map can be applied iteratively to predict solutions in the long range. Our trained networks, especially in two-phase flow where the maps are nonlinear, show their great potential in accurately approximating the underlying physical system and improvement in computational efficiency. Some numerical experiments are performed and discussed to demonstrate the performance of our proposed techniques.


title: Steady states, squeezing, and entanglement in intracavity triplet down conversion

abstract: Triplet down conversion, the process of converting one high-energy photon into three low-energy photons, may soon be experimentally feasible due to advances in optical resonator technology. We use quantum phase-space techniques to analyse the process of degenerate intracavity triplet down conversion by solving stochastic differential equations within the truncated positive-P representation. The time evolution of both intracavity mode populations are simulated, and the resulting steady-states are examined as a function of the pump intensity. Quantum effects are most pronounced in the region immediately above the semi-classical pumping threshold, where our numerical results differ significantly from semi-classical predictions. Regimes of measurable squeezing and bipartite entanglement are identified from steady-state spectra of the cavity output fields. We validate the truncated positive-P description against Monte Carlo wave function simulations, finding good agreement for low mode populations.


title: Orbital dynamics in realistic galaxy models: NGC 3726, NGC 3877 and NGC 4010

abstract: In the present paper, using a generalization of the Miyamoto and Nagai potential we adjusted the observed rotation curves of three specific spiral galaxies to the analytical circular velocities. The observational data have been taken from a 21 cm-line synthesis imaging survey using the Westerbork Synthesis Radio Telescope, for three particular galaxies in the Ursa Major cluster: NGC 3726, NGC 3877 and NGC 4010. Accordingly, the dynamics of the system is analyzed in terms of the Poincar\'e sections method, finding that for larger values of the angular momentum of the test particle or lower values its total energy the dynamics is mainly regular, while on the opposite cases, the dynamics is mainly chaotic. Our toy model opens the possibility to find chaotic bounded orbits for stars in those particular galaxies.


title: LEGO: Leveraging Experience in Roadmap Generation for Sampling-Based Planning

abstract: We consider the problem of leveraging prior experience to generate roadmaps in sampling-based motion planning. A desirable roadmap is one that is sparse, allowing for fast search, with nodes spread out at key locations such that a low-cost feasible path exists. An increasingly popular approach is to learn a distribution of nodes that would produce such a roadmap. State-of-the-art is to train a conditional variational auto-encoder (CVAE) on the prior dataset with the shortest paths as target input. While this is quite effective on many problems, we show it can fail in the face of complex obstacle configurations or mismatch between training and testing. We present an algorithm LEGO that addresses these issues by training the CVAE with target samples that satisfy two important criteria. Firstly, these samples belong only to bottleneck regions along near-optimal paths that are otherwise difficult-to-sample with a uniform sampler. Secondly, these samples are spread out across diverse regions to maximize the likelihood of a feasible path existing. We formally define these properties and prove performance guarantees for LEGO. We extensively evaluate LEGO on a range of planning problems, including robot arm planning, and report significant gains over heuristics as well as learned baselines.


title: A 2D Parallel Triangle Counting Algorithm for Distributed-Memory Architectures

abstract: Triangle counting is a fundamental graph analytic operation that is used extensively in network science and graph mining. As the size of the graphs that needs to be analyzed continues to grow, there is a requirement in developing scalable algorithms for distributed-memory parallel systems. To this end, we present a distributed-memory triangle counting algorithm, which uses a 2D cyclic decomposition to balance the computations and reduce the communication overheads. The algorithm structures its communication and computational steps such that it reduces its memory overhead and includes key optimizations that leverage the sparsity of the graph and the way the computations are structured. Experiments on synthetic and real-world graphs show that our algorithm obtains an average relative speedup that range between 3.24 and 7.22 out of 10.56 across the datasets using 169 MPI ranks over the performance achieved by 16 MPI ranks. Moreover, we obtain an average speedup of 10.2 times on comparison with previously developed distributed-memory parallel algorithms.


title: Studying transverse momentum distributions with jets at N$^3$LL

abstract: Semi-inclusive deep inelastic scattering (SIDIS) is a promising channel for the extraction of transverse momentum dependent distributions at future colliders. In this context, we recently developed a framework that uses jets (instead of single hadrons) to achieve reduced sensitivity to final-state non-perturbative effects. A suitable non-standard jet definition allows us to apply the factorization formulas valid for hadrons to jets of arbitrary size, by just replacing fragmentation functions with the jet functions we computed. Besides presenting the framework, we will show numerical predictions at N$^3$LL accuracy.


title: Examples of weak amalgamation classes

abstract: We present several examples of hereditary classes of finite structures satisfying the joint embedding property and the weak amalgamation property, but failing the cofinal amalgamation property. These include a continuum-sized family of classes of finite undirected graphs, as well as an example due to Pouzet with countably categorical generic limit.


title: Information-Bottleneck Approach to Salient Region Discovery

abstract: We propose a new method for learning image attention masks in a semi-supervised setting based on the Information Bottleneck principle. Provided with a set of labeled images, the mask generation model is minimizing mutual information between the input and the masked image while maximizing the mutual information between the same masked image and the image label. In contrast with other approaches, our attention model produces a Boolean rather than a continuous mask, entirely concealing the information in masked-out pixels. Using a set of synthetic datasets based on MNIST and CIFAR10 and the SVHN datasets, we demonstrate that our method can successfully attend to features known to define the image class.


title: ZKlaims: Privacy-preserving Attribute-based Credentials using Non-interactive Zero-knowledge Techniques

abstract: In this paper we present ZKlaims: a system that allows users to present attribute-based credentials in a privacy-preserving way. We achieve a zero-knowledge property on the basis of Succinct Non-interactive Arguments of Knowledge (SNARKs). ZKlaims allow users to prove statements on credentials issued by trusted third parties. The credential contents are never revealed to the verifier as part of the proving process. Further, ZKlaims can be presented non-interactively, mitigating the need for interactive proofs between the user and the verifier. This allows ZKlaims to be exchanged via fully decentralized services and storages such as traditional peer-to-peer networks based on distributed hash tables (DHTs) or even blockchains. To show this, we include a performance evaluation of ZKlaims and show how it can be integrated in decentralized identity provider services.


title: Interplay of pulse duration, peak intensity, and particle size in laser-driven electron emission from silica nanospheres

abstract: We present the results of a systematic study of photoelectron emission from gasphase dielectric nanoparticles (SiO2) irradiated by intense 25 fs, 780 nm linearly polarized laser pulses as a function of particle size (20 nm to 750 nm in diameter) and laser intensity. We also introduce an experimental technique to reduce the effects of focal volume averaging. The highest photoelectron energies show a strong size dependence, increasing by a factor of six over the range of particles sizes studied at a fixed intensity. For smaller particle sizes (up to 200 nm), our findings agree well with earlier results obtained with few-cycle, ~4 fs pulses. For large nanoparticles, which exhibit stronger near-field localization due to field-propagation effects, we observe the emission of much more energetic electrons, reaching energies up to ~200 times the ponderomotive energy. This strong deviation in maximum photoelectron energy is attributed to the increase in ionization and charge interaction for many-cycle pulses at similar intensities.


title: Barrier from chaos: operator entanglement dynamics of the reduced density matrix

abstract: It is believed that thermalization drives the reduced density matrix of a subsystem to approach a short-range entangled operator. If the initial state is also short-range entangled, it is possible that the reduced density matrix remains low-entangled throughout thermalization; or there could exist a barrier with high operator entanglement between the initial and thermalized reduced density matrix. In this paper, we study such dynamics in three classes of models: the rational CFTs, the random unitary circuit, and the holographic CFTs, representing systems of increasing quantum chaoticity. We show that in all three classes of models, the operator entanglement (or variant of) exhibits three phases, a linear growth phase, a plateau phase, and a decay phase. The plateau phase characterized by volume-law operator entanglement corresponds to the barrier in operator entanglement. While it is present in all three models, its persistence and exit show interesting distinctions among them. The rational CFTs have the shortest plateau phase, followed by the slowest decay phase; the holographic CFTs mark the opposite end, i.e. having the longest plateau phase followed by a discontinuous drop; and the random unitary circuit shows the intermediate behavior. We discuss the mechanisms underlying these behaviors in operator entanglement barriers, whose persistence might serve as another measure for quantum chaoticity.


title: The $k$-Dimensional Weisfeiler-Leman Algorithm

abstract: In this note, we provide details of the $k$-dimensional Weisfeiler-Leman Algorithm and its analysis from Immerman-Lander (1990). In particular, we present an optimized version of the algorithm that runs in time $O(n^{k+1}\log n)$, where $k$ is fixed (not varying with $n$).


title: A Realistic Roadmap to Formation Flying Space Interferometry

abstract: The ultimate astronomical observatory would be a formation flying interferometer in space, immune to atmospheric turbulence and absorption, free from atmospheric and telescope thermal emission, and reconfigurable to adjust baselines according to the required angular resolution. Imagine the near/mid-infrared sensitivity of the JWST and the far-IR sensitivity of Herschel but with ALMA-level angular resolution, or imagine having the precision control to null host star light across 250m baselines and to detect molecules from the atmospheres of nearby exo-Earths. With no practical engineering limit to the formation's size or number of telescopes in the array, formation flying interferometry will revolutionize astronomy and this White Paper makes the case that it is now time to accelerate investments in this technological area. Here we provide a brief overview of the required technologies needed to allow light to be collected and interfered using separate spacecrafts. We emphasize the emerging role of inexpensive smallSat projects and the excitement for the LISA Gravitational Wave Interferometer to push development of the required engineering building-blocks. We urge the Astro2020 Decadal Survey Committee to highlight the need for a small-scale formation flying space interferometer project to demonstrate end-to-end competency with a timeline for first stellar fringes by the end of the decade.


title: NAT Black Holes

abstract: We study some physical properties of black holes in Null Aether Theory (NAT)--a vector-tensor theory of gravity. We first review the black hole solutions in NAT and then derive the first law of black hole thermodynamics. The temperature of the black holes depends on both the mass and the NAT \textquotedblleft charge" of the black holes. The extreme cases where the temperature vanishes resemble the extreme Reissner-Nordstr\"{o}m black holes. We also discuss the contribution of the NAT charge to the geodesics of massive and massless particles around the NAT black holes.


title: Cooperative Pollution Source Localization and Cleanup with a Bio-inspired Swarm Robot Aggregation

abstract: Using robots for exploration of extreme and hazardous environments has the potential to significantly improve human safety. For example, robotic solutions can be deployed to find the source of a chemical leakage and clean the contaminated area. This paper demonstrates a proof-of-concept bio-inspired exploration method using a swarm robotic system, which is based on a combination of two bio-inspired behaviours: aggregation, and pheromone tracking. The main idea of the work presented is to follow pheromone trails to find the source of a chemical leakage and then carry out a decontamination task by aggregating at the critical zone. Using experiments conducted by a simulated model of a Mona robot, we evaluate the effects of population size and robot speed on the ability of the swarm in a decontamination task. The results indicate the feasibility of deploying robotic swarms in an exploration and cleaning task in an extreme environment.


title: Good Neighbors, Bad Neighbors: The Frequent Network Neighborhood Mapping of the Hippocampus Enlightens Several Structural Factors of the Human Intelligence on a 414-Subject Cohort

abstract: The human connectome has become the very frequent subject of study of brain-scientists, psychologists, and imaging experts in the last decade. With diffusion magnetic resonance imaging techniques, unified with advanced data processing algorithms, today we are able to compute braingraphs with several hundred, anatomically identified nodes and thousands of edges, corresponding to the anatomical connections of the brain. The analysis of these graphs without refined mathematical tools is hopeless. These tools need to address the high error rate of the MRI processing workflow, and need to find structural causes or at least correlations of psychological properties and cerebral connections. Until now, structural connectomics was only rarely able identifying such causes or correlations. In the present work, we study the frequent neighbor sets of the most deeply investigated brain area, the hippocampus. By applying the Frequent Network Neighborhood mapping method, we identified frequent neighbor-sets of the hippocampus, which may influence numerous psychological parameters, including intelligence-related ones. We have found neighbor sets, which have significantly higher frequency in subjects with high-scored Penn Matrix tests, and with low-scored Penn Word Memory tests. Our study utilizes the braingraphs, computed from the imaging data of the Human Connectome Project's 414 subjects, each with 463 anatomically identified nodes.


title: Feller coupling of cycles and Poisson spacings

abstract: Feller (1945) provided a coupling between the counts of cycles of various sizes in a uniform random permutation of $[n]$ and the spacings between successes in a sequence of $n$ independent Bernoulli trials with success probability $1/n$ at the $n$th trial. Arratia, Barbour and Tavar\'e (1992) extended Feller's coupling, to associate cycles of random permutations governed by the Ewens $(\theta)$ distribution with spacings derived from independent Bernoulli trials with success probability $\theta/(n-1+\theta)$ at the $n$th trial, and to conclude that in an infinite sequence of such trials, the numbers of spacings of length $\ell$ are independent Poisson variables with means $\theta/\ell$. Ignatov (1978) first discovered this remarkable result in the uniform case $\theta = 1$, by constructing Bernoulli $(1/n)$ trials as the indicators of record values in a sequence of i.i.d. uniform $[0,1]$ variables. In the present article, the Poisson property of inhomogeneous Bernoulli spacings is explained by a variation of Ignatov's approach for a general $\theta >0$. Moreover, our approach naturally provides random permutations of infinite sets whose cycle counts are exactly given by independent Poisson random variables.


title: Direction Matters: On Influence-Preserving Graph Summarization and Max-cut Principle for Directed Graphs

abstract: Summarizing large-scaled directed graphs into small-scale representations is a useful but less studied problem setting. Conventional clustering approaches, which based on "Min-Cut"-style criteria, compress both the vertices and edges of the graph into the communities, that lead to a loss of directed edge information. On the other hand, compressing the vertices while preserving the directed edge information provides a way to learn the small-scale representation of a directed graph. The reconstruction error, which measures the edge information preserved by the summarized graph, can be used to learn such representation. Compared to the original graphs, the summarized graphs are easier to analyze and are capable of extracting group-level features which is useful for efficient interventions of population behavior. In this paper, we present a model, based on minimizing reconstruction error with non-negative constraints, which relates to a "Max-Cut" criterion that simultaneously identifies the compressed nodes and the directed compressed relations between these nodes. A multiplicative update algorithm with column-wise normalization is proposed. We further provide theoretical results on the identifiability of the model and on the convergence of the proposed algorithms. Experiments are conducted to demonstrate the accuracy and robustness of the proposed method.


title: Voltage Security Analysis of VSC-HVDC Transmission Lines

abstract: Due to generation retirements and growth of renewable energy sources integration, there will be widespread changes in the real and reactive power flow in nowadays power systems. These changes also bring about challenges on system operation and stability. High voltage direct current (HVDC) technology seems to be a promising solution to these new challenges. A study on the impact of the replacement of 500KV AC transmission lines by VSC-HVDC transmission line on system voltage security is conducted. The analysis is based on reactive load margin method and implemented on one of the Dominion Energy planning models. Different control schemes of HVDC are considered. Using k-means clustering method, three representative zones within the network are selected. The results corresponding to them are demonstrated and discussed. It is shown that HVDC lines with P-V control remarkably improve system voltage security, while those with P-PF control scheme have negative effects.


title: Fatou's Theorem for Non-commutative Measures

abstract: A classical theorem of Fatou asserts that the Radon-Nikodym derivative of any finite positive Borel measure, $\mu$, with respect to Lebesgue measure on the complex unit circle, is recovered as the non-tangential limits of its Poisson transform in the complex unit disk. This positive harmonic Poisson transform is the real part of an analytic function whose Taylor coefficients are in fixed proportion to the conjugate moments of $\mu$. Replacing Taylor series in one variable by power series in several non-commuting variables, we show that Fatou's Theorem and related results have natural extensions to the setting of positive harmonic functions in an open unit ball of several non-commuting matrix-variables, and a corresponding class of positive \emph{non-commutative (NC) measures}. Here, an NC measure is any positive linear functional on a certain self-adjoint unital subspace of the Cuntz-Toeplitz algebra, the $C^*-$algebra generated by the left creation operators on the full Fock space.


title: Differentiable Gaussian Process Motion Planning

abstract: Modern trajectory optimization based approaches to motion planning are fast, easy to implement, and effective on a wide range of robotics tasks. However, trajectory optimization algorithms have parameters that are typically set in advance (and rarely discussed in detail). Setting these parameters properly can have a significant impact on the practical performance of the algorithm, sometimes making the difference between finding a feasible plan or failing at the task entirely. We propose a method for leveraging past experience to learn how to automatically adapt the parameters of Gaussian Process Motion Planning (GPMP) algorithms. Specifically, we propose a differentiable extension to the GPMP2 algorithm, so that it can be trained end-to-end from data. We perform several experiments that validate our algorithm and illustrate the benefits of our proposed learning-based approach to motion planning.


title: Few-Layer MoS$_2$/a-Si:H Heterojunction pin-Photodiodes for extended Infrared Detection

abstract: Few-layer molybdenum disulfide (FL-MoS$_2$) films have been integrated into amorphous silicon (a-Si:H) pin photodetectors. To achieve this, vertical a-Si:H photodiodes were grown by plasma-enhanced chemical vapor deposition (PE-CVD) on top of large-scale synthesized and transferred homogeneous FL-MoS$_2$. This novel detector array exhibits long-term stability (more than six month) and outperforms conventional silicon-based pin photodetectors in the infrared range (IR, $\lambda$ = 2120 nm) in terms of sensitivities by up to 50 mAW$^{-1}$. Photodetectivities of up to 2 x 10$^{10}$ Jones and external quantum efficiencies of 3 % are achieved. The detectors further feature the additional functionality of bias-dependent responsivity switching between the different spectral ranges. The realization of such scalable detector arrays is an essential step towards pixelated and wavelength-selective sensors operating in the IR spectral range.


title: Achievable Rate Region for Iterative Multi-User Detection via Low-cost Gaussian Approximation

abstract: We establish a multi-user extrinsic information transfer (EXIT) chart area theorem for the interleave-division multiple-access (IDMA) scheme, a special form of superposition coding, in multiple access channels (MACs). A low-cost multi-user detection (MUD) based on the Gaussian approximation (GA) is assumed. The evolution of mean-square errors (MSE) of the GA-based MUD during iterative processing is studied. We show that the K-dimensional tuples formed by the MSEs of K users constitute a conservative vector field. The achievable rate is a potential function of this conservative field, so it is the integral along any path in the field with value of the integral solely determined by the two path terminals. Optimized codes can be found given the integration paths in the MSE fields by matching EXIT type functions. The above findings imply that i) low-cost GA-based MUD can provide near capacity performance; ii) the sum-rate capacity (region) can be achieved independently of the integration path in the MSE fields; and iii) the integration path can be an extra degree of freedom for code design.


title: Understanding the Political Ideology of Legislators from Social Media Images

abstract: In this paper, we seek to understand how politicians use images to express ideological rhetoric through Facebook images posted by members of the U.S. House and Senate. In the era of social media, politics has become saturated with imagery, a potent and emotionally salient form of political rhetoric which has been used by politicians and political organizations to influence public sentiment and voting behavior for well over a century. To date, however, little is known about how images are used as political rhetoric. Using deep learning techniques to automatically predict Republican or Democratic party affiliation solely from the Facebook photographs of the members of the 114th U.S. Congress, we demonstrate that predicted class probabilities from our model function as an accurate proxy of the political ideology of images along a left-right (liberal-conservative) dimension. After controlling for the gender and race of politicians, our method achieves an accuracy of 59.28% from single photographs and 82.35% when aggregating scores from multiple photographs (up to 150) of the same person. To better understand image content distinguishing liberal from conservative images, we also perform in-depth content analyses of the photographs. Our findings suggest that conservatives tend to use more images supporting status quo political institutions and hierarchy maintenance, featuring individuals from dominant social groups, and displaying greater happiness than liberals.


title: MixConv: Mixed Depthwise Convolutional Kernels

abstract: Depthwise convolution is becoming increasingly popular in modern efficient ConvNets, but its kernel size is often overlooked. In this paper, we systematically study the impact of different kernel sizes, and observe that combining the benefits of multiple kernel sizes can lead to better accuracy and efficiency. Based on this observation, we propose a new mixed depthwise convolution (MixConv), which naturally mixes up multiple kernel sizes in a single convolution. As a simple drop-in replacement of vanilla depthwise convolution, our MixConv improves the accuracy and efficiency for existing MobileNets on both ImageNet classification and COCO object detection. To demonstrate the effectiveness of MixConv, we integrate it into AutoML search space and develop a new family of models, named as MixNets, which outperform previous mobile models including MobileNetV2 [20] (ImageNet top-1 accuracy +4.2%), ShuffleNetV2 [16] (+3.5%), MnasNet [26] (+1.3%), ProxylessNAS [2] (+2.2%), and FBNet [27] (+2.0%). In particular, our MixNet-L achieves a new state-of-the-art 78.9% ImageNet top-1 accuracy under typical mobile settings (<600M FLOPS). Code is at https://github.com/ tensorflow/tpu/tree/master/models/official/mnasnet/mixnet


title: Gapless Dirac surface states in the antiferromagnetic topological insulator MnBi2Te4

abstract: We use high-resolution, tunable angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) calculations to study the electronic properties of single crystals of MnBi2Te4, a material that was predicted to be the first intrinsic antiferromagnetic (AFM) topological insulator. We observe both bulk and surface bands in the electronic spectra, in reasonable agreement with the DFT calculations results. In striking contrast to the earlier literatures showing a full gap opening between two surface band manifolds along (0001) direction, we observed a gapless Dirac cone remain protected in MnBi2Te4 across the AFM transition (TN = 24 K). Our data also reveal the existence of a second Dirac cone closer to the Fermi level, predicted by band structure calculations. Whereas the surface Dirac cones seem to be remarkably insensitive to the AFM ordering, we do observe splitting of the bulk band that develops below the TN . Having a moderately high ordering temperature, MnBi2Te4 provides a unique platform for studying the interplay between topology and magnetic ordering.


title: Agent Modeling as Auxiliary Task for Deep Reinforcement Learning

abstract: In this paper we explore how actor-critic methods in deep reinforcement learning, in particular Asynchronous Advantage Actor-Critic (A3C), can be extended with agent modeling. Inspired by recent works on representation learning and multiagent deep reinforcement learning, we propose two architectures to perform agent modeling: the first one based on parameter sharing, and the second one based on agent policy features. Both architectures aim to learn other agents' policies as auxiliary tasks, besides the standard actor (policy) and critic (values). We performed experiments in both cooperative and competitive domains. The former is a problem of coordinated multiagent object transportation and the latter is a two-player mini version of the Pommerman game. Our results show that the proposed architectures stabilize learning and outperform the standard A3C architecture when learning a best response in terms of expected rewards.


title: Gibbs energy of ices III, V and VI: wholistic thermodynamics and elasticity of the water phase diagram to 2300 MPa

abstract: Gibbs energy representations for ice III, V and VI are reported. These were constructed using new measurements of volumes at high pressure over a range of low temperatures combined with calculated vibrational energies grounded in statistical physics. The collection of representations including ice Ih and water (released as the open source SeaFreeze framework) allow accurate determinations of thermodynamics properties (phase boundaries, density, heat capacity, bulk modulus, thermal expansivity, chemical potentials) and seismic wave velocities over the entire range of conditions encountered in hydrospheres in our solar system (220 - 500K to 2300 MPa). These comprehensive representations allow exploration of the rich spectrum of thermodynamic behavior in the H2O system. Although the results are broadly applicable in science and engineering, their use in habitability analysis in water-rich planetary bodies of our solar system and beyond is particularly relevant.


title: The essential numerical range for unbounded linear operators

abstract: We introduce the concept of essential numerical range $W_{\!e}(T)$ for unbounded Hilbert space operators $T$ and study its fundamental properties including possible equivalent characterizations and perturbation results. Many of the properties known for the bounded case do \emph{not} carry over to the unbounded case, and new interesting phenomena arise which we illustrate by some striking examples. A key feature of the essential numerical range $W_{\!e}(T)$ is that it captures spectral pollution in a unified and minimal way when approximating $T$ by projection methods or domain truncation methods for PDEs.


title: Evaluation of Embeddings of Laboratory Test Codes for Patients at a Cancer Center

abstract: Laboratory test results are an important and generally high dimensional component of a patient's Electronic Health Record (EHR). We train embedding representations (via Word2Vec and GloVe) for LOINC codes of laboratory tests from the EHRs of about 80,000 patients at a cancer center. To include information about lab test outcomes, we also train embeddings on the concatenation of a LOINC code with a symbol indicating normality or abnormality of the result. We observe several clinically meaningful similarities among LOINC embeddings trained over our data. For the embeddings of the concatenation of LOINCs with abnormality codes, we evaluate the performance for mortality prediction tasks and the ability to preserve ordinality properties: i.e. a lab test with normal outcome should be more similar to an abnormal one than to the a very abnormal one.


title: Emergent Kalb-Ramond fields from a dimer model

abstract: The emergence of a Kalb-Ramond field and string charge in the lattice is discussed. The local bosonic model with rotor variables placed on the faces of a cubic lattice is considered. The coupling model consisting of the Maxwell fields and the Kalb-Ramond field is given. This construction naturally incorporates the emerging coupling between both gauge and string fields. In the process, an object that resembles to a D-brane on the lattice is introduced.


title: Steganography Protocols for Quantum Channels

abstract: We study several versions of a quantum steganography problem, in which two legitimate parties attempt to conceal a cypher in a quantum cover transmitted over a quantum channel without arising suspicion from a warden who intercepts the cover. In all our models, we assume that the warden has an inaccurate knowledge of the quantum channel and we formulate several variations of the steganography problem depending on the tasks used as the cover and the cypher task. In particular, when the cover task is classical communication, we show that the cypher task can be classical communication or entanglement sharing; when the cover task is entanglement sharing and the main channel is noiseless, we show that the cypher task can be randomness sharing; when the cover task is quantum communication and the main channel is noiseless, we show that the cypher task can be classical communication. In the latter case, our results improve earlier ones by relaxing the need for a shared key between the transmitter and the receiver and hold under milder assumptions on the cover quantum communication code.


title: Correct-by-Construction Advanced Driver Assistance Systems based on a Cognitive Architecture

abstract: Research into safety in autonomous and semi-autonomous vehicles has, so far, largely been focused on testing and validation through simulation. Due to the fact that failure of these autonomous systems is potentially life-endangering, formal methods arise as a complementary approach. This paper studies the application of formal methods to the verification of a human driver model built using the cognitive architecture ACT-R, and to the design of correct-by-construction Advanced Driver Assistance Systems (ADAS). The novelty lies in the integration of ACT-R in the formal analysis and an abstraction technique that enables finite representation of a large dimensional, continuous system in the form of a Markov process. The situation considered is a multi-lane highway driving scenario and the interactions that arise. The efficacy of the method is illustrated in two case studies with various driving conditions.


title: Fluctuation relations and strong inequalities for thermally isolated systems

abstract: For processes during which a macroscopic system exchanges no heat with its surroundings, the second law of thermodynamics places two lower bounds on the amount of work performed on the system: a weak bound, expressed in terms of a fixed-temperature free energy difference, $W \ge \Delta F_T$ , and a strong bound, given by a fixed-entropy internal energy difference, $W \ge \Delta E_S$ . It is known that statistical inequalities related to the weak bound can be obtained from the nonequilibrium work relation, $\langle\exp (-\beta W)\rangle = \exp(-\beta\Delta F_T)$ . Here we derive an integral fluctuation relation $\langle\exp(-\beta X) \rangle = 1 $ that is constructed specifically for adiabatic processes, and we use this result to obtain inequalities related to the strong bound, $W \ge \Delta E_S$ . We provide both classical and quantum derivations of these results.


title: Bilevel Optimization, Deep Learning and Fractional Laplacian Regularization with Applications in Tomography

abstract: In this work we consider a generalized bilevel optimization framework for solving inverse problems. We introduce fractional Laplacian as a regularizer to improve the reconstruction quality, and compare it with the total variation regularization. We emphasize that the key advantage of using fractional Laplacian as a regularizer is that it leads to a linear operator, as opposed to the total variation regularization which results in a nonlinear degenerate operator. Inspired by residual neural networks, to learn the optimal strength of regularization and the exponent of fractional Laplacian, we develop a dedicated bilevel optimization neural network with a variable depth for a general regularized inverse problem. We also draw some parallels between an activation function in a neural network and regularization. We illustrate how to incorporate various regularizer choices into our proposed network. As an example, we consider tomographic reconstruction as a model problem and show an improvement in reconstruction quality, especially for limited data, via fractional Laplacian regularization. We successfully learn the regularization strength and the fractional exponent via our proposed bilevel optimization neural network. We observe that the fractional Laplacian regularization outperforms total variation regularization. This is specially encouraging, and important, in the case of limited and noisy data.


title: Evaluating the Effect of Geometry and Control on Freeway Merge Bottleneck Capacity

abstract: This paper evaluates the impact of geometric and operational features on freeway capacity at merge bottleneck locations, by analyzing pre-breakdown flow rates and using parametric and non-parametric techniques. These features include number of lanes, free-flow speed or speed limit, length of acceleration lane, and presence or ramp metering. The analysis was carried out on seventeen freeway merge sites across the United States. The number of lanes showed negative relationship with the per-lane average pre-breakdown flow rate. Sites equipped with ramp meters also showed higher pre-breakdown flow rates than unmetered sites. In survival analysis, the number of lanes and presence of ramp meters were the only features found to have a statistically significant impact on the survival (or breakdown) probability. Capacities, defined at the flow rate that corresponds to 15 percent breakdown probability, were 2,048, 1,959, and 1,745 pc/hr/ln at three-lane, four-lane, and five-lane unmetered sites, respectively. Capacities were 2,248 pc/hr/ln at three-lane and 2,132 pc/hr/ln at four-lane metered sites. The effect of acceleration lane length and Free Flow Speed (FFS) on pre-breakdown flow rate and survival probability was not found to be significant.


title: Semi-Supervised Learning by Disentangling and Self-Ensembling Over Stochastic Latent Space

abstract: The success of deep learning in medical imaging is mostly achieved at the cost of a large labeled data set. Semi-supervised learning (SSL) provides a promising solution by leveraging the structure of unlabeled data to improve learning from a small set of labeled data. Self-ensembling is a simple approach used in SSL to encourage consensus among ensemble predictions of unknown labels, improving generalization of the model by making it more insensitive to the latent space. Currently, such an ensemble is obtained by randomization such as dropout regularization and random data augmentation. In this work, we hypothesize -- from the generalization perspective -- that self-ensembling can be improved by exploiting the stochasticity of a disentangled latent space. To this end, we present a stacked SSL model that utilizes unsupervised disentangled representation learning as the stochastic embedding for self-ensembling. We evaluate the presented model for multi-label classification using chest X-ray images, demonstrating its improved performance over related SSL models as well as the interpretability of its disentangled representations.


title: Classical Balayage of Charges and Measures

abstract: We investigate some properties of balayage of charges and measures for subclasses of subharmonic functions and their relationship to the geometry of domain or open set in finite-dimensional Euclidean space where this balayage is considered.


title: Remarks on a limiting case of Hardy type inequalities

abstract: The classical Hardy inequality holds in Sobolev spaces $W_0^{1,p}$ when $1\le p< N$. In the limiting case where $p=N$, it is known that by adding a logarithmic function to the Hardy potential, some inequality which is called the critical Hardy inequality holds in $W_0^{1,N}$. In this note, in order to give an explanation of appearance of the logarithmic function at the potential, we derive the logarithmic function from the classical Hardy inequality with the best constant via some limiting procedure as $p \nearrow N$. And we show that our limiting procedure is also available for the classical Rellich inequality in second order Sobolev spaces $W_0^{2,p}$ with $p \in (1, \frac{N}{2})$ and the Poincar\'e inequality.


title: Nonlinear Reconnection in Magnetized Turbulence

abstract: Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those conditions are met, a new turbulence regime is accessed where reconnection-driven energy dissipation becomes common, rather than the rare feature that it must be when they are not.


title: Asymptotic normality, concentration, and coverage of generalized posteriors

abstract: Generalized likelihoods are commonly used to obtain consistent estimators with attractive computational and robustness properties. Formally, any generalized likelihood can be used to define a generalized posterior distribution, but an arbitrarily defined "posterior" cannot be expected to appropriately quantify uncertainty in any meaningful sense. In this article, we provide sufficient conditions under which generalized posteriors exhibit concentration, asymptotic normality (Bernstein-von Mises), an asymptotically correct Laplace approximation, and asymptotically correct frequentist coverage. We apply our results in detail to generalized posteriors for a wide array of generalized likelihoods, including pseudolikelihoods in general, the Gaussian Markov random field pseudolikelihood, the fully observed Boltzmann machine pseudolikelihood, the Ising model pseudolikelihood, the Cox proportional hazards partial likelihood, and a median-based likelihood for robust inference of location. Further, we show how our results can be used to easily establish the asymptotics of standard posteriors for exponential families and generalized linear models. We make no assumption of model correctness so that our results apply with or without misspecification.


title: Political network of central power agents: case of missi dominici

abstract: This study offers several models of social network analysis to examine the organization of central power agents, missi dominici, during the Early Middle Ages. Enriched by statistical analysis, different research hypotheses based on the current historiographical positions have been substantiated. On the one side, the network analysis allowed to highlight the evolution of network structure throughout the studied period and to observe a change in the framework of agents transition between reigns. On the other side, the statistical exploration of the relations between the agents and the places of their assignments confirmed some amplification, with time, of a tendency to recruit the agents among the local aristocracy. Finally, several difficulties related to the analyzing of missing data provided by fragmentary historical records as well as to modeling a complex multimodal political network were mentioned.


title: Incremental and Decremental Fuzzy Bounded Twin Support Vector Machine

abstract: In this paper we present an incremental variant of the Twin Support Vector Machine (TWSVM) called Fuzzy Bounded Twin Support Vector Machine (FBTWSVM) to deal with large datasets and learning from data streams. We combine the TWSVM with a fuzzy membership function, so that each input has a different contribution to each hyperplane in a binary classifier. To solve the pair of quadratic programming problems (QPPs) we use a dual coordinate descent algorithm with a shrinking strategy, and to obtain a robust classification with a fast training we propose the use of a Fourier Gaussian approximation function with our linear FBTWSVM. Inspired by the shrinking technique, the incremental algorithm re-utilizes part of the training method with some heuristics, while the decremental procedure is based on a scored window. The FBTWSVM is also extended for multi-class problems by combining binary classifiers using a Directed Acyclic Graph (DAG) approach. Moreover, we analyzed the theoretical foundations properties of the proposed approach and its extension, and the experimental results on benchmark datasets indicate that the FBTWSVM has a fast training and retraining process while maintaining a robust classification performance.


title: Shape of a recoiling liquid filament

abstract: We study the capillary retraction of a Newtonian semi-infinite liquid filament through analytical methods. We derive a long-time asymptotic-state expansion for the filament profile using a one-dimensional free-surface slender cylindrical flow model based on the three-dimensional axisymmetric Navier-Stokes equations. The analysis identifies three distinct length and time scale regions in the retraction domain: a steady filament section, a growing spherical blob, and an intermediate matching zone. We show that liquid filaments naturally develop travelling capillary waves along their surface and a neck behind the blob. We analytically prove that the wavelength of the capillary waves is approximately 3.63 times the filament's radius at the inviscid limit. Additionally, the waves' asymptotic wavelength, decay length, and the minimum neck size are analysed in terms of the Ohnesorge number. Finally, our findings are compared with previous results from the literature and numerical simulations in Basilisk obtaining a good agreement. This analysis provides a full picture of the recoiling process going beyond the classic result of the velocity of retraction found by Taylor and Culick.


title: Towards Realistic Individual Recourse and Actionable Explanations in Black-Box Decision Making Systems

abstract: Machine learning based decision making systems are increasingly affecting humans. An individual can suffer an undesirable outcome under such decision making systems (e.g. denied credit) irrespective of whether the decision is fair or accurate. Individual recourse pertains to the problem of providing an actionable set of changes a person can undertake in order to improve their outcome. We propose a recourse algorithm that models the underlying data distribution or manifold. We then provide a mechanism to generate the smallest set of changes that will improve an individual's outcome. This mechanism can be easily used to provide recourse for any differentiable machine learning based decision making system. Further, the resulting algorithm is shown to be applicable to both supervised classification and causal decision making systems. Our work attempts to fill gaps in existing fairness literature that have primarily focused on discovering and/or algorithmically enforcing fairness constraints on decision making systems. This work also provides an alternative approach to generating counterfactual explanations.


title: Factors influencing thermal solidification of bent-core trimers

abstract: Bent-core trimers are a simple model system for which the competition between crystallization and glass-formation can be tuned by varying a single parameter:\ the bond angle $\theta_0$. Using molecular dynamics simulations, we examine how varying $\theta_0$ affects their thermal solidification. By examining trends with $\theta_0$, comparing these to trends in trimers' jamming phenomenology, and then focusing on the six $\theta_0$ that are commensurable with close-packed crystalline order, we obtain three key results: \textbf{(i)} the increase in trimers' solidification temperature $T_s(\theta_0)$ as they straighten (as $\theta_0 \to 0^\circ$) is driven by the same gradual loss of \textit{effective} configurational freedom that drives athermal trimers' decreasing $\phi_J(\theta_0)$; \textbf{(ii)} $\theta_0$ that allow formation of both FCC and HCP order crystallize, while $\theta_0$ that only allow formation of HCP order glass-form; \textbf{(iii)} local cluster-level structure at temperatures slightly \textit{above} $T_s(\theta_0)$ is highly predictive of whether trimers will crystallize or glass-form.


title: Hierarchical Transformed Scale Mixtures for Flexible Modeling of Spatial Extremes on Datasets with Many Locations

abstract: Flexible spatial models that allow transitions between tail dependence classes have recently appeared in the literature. However, inference for these models is computationally prohibitive, even in moderate dimensions, due to the necessity of repeatedly evaluating the multivariate Gaussian distribution function. In this work, we attempt to achieve truly high-dimensional inference for extremes of spatial processes, while retaining the desirable flexibility in the tail dependence structure, by modifying an established class of models based on scale mixtures Gaussian processes. We show that the desired extremal dependence properties from the original models are preserved under the modification, and demonstrate that the corresponding Bayesian hierarchical model does not involve the expensive computation of the multivariate Gaussian distribution function. We fit our model to exceedances of a high threshold, and perform coverage analyses and cross-model checks to validate its ability to capture different types of tail characteristics. We use a standard adaptive Metropolis algorithm for model fitting, and further accelerate the computation via parallelization and Rcpp. Lastly, we apply the model to a dataset of a fire threat index on the Great Plains region of the US, which is vulnerable to massively destructive wildfires. We find that the joint tail of the fire threat index exhibits a decaying dependence structure that cannot be captured by limiting extreme value models.


title: Experimental proof of Quantum Zeno-assisted Noise Sensing

abstract: In the ideal quantum Zeno effect, repeated quantum projective measurements can freeze the coherent dynamics of a quantum system. However, in the weak quantum Zeno regime, measurement back-actions can allow the sensing of semi-classical field fluctuations. In this regard, we theoretically show how to combine the controlled manipulation of a quantum two-level system, used as a probe, with a sequence of projective measurements to have direct access to the noise correlation function. We experimentally test the effectiveness of the proposed noise sensing method on a properly engineered Bose-Einstein condensate of $^{87}Rb$ atoms realized on an atom chip. We believe that our quantum Zeno-based approach can open a new path towards novel quantum sensing devices.


title: Noncommutative cosmological models induced by a symplectic formalism coupled to phantom fluids

abstract: In the present letter, we consider homogeneous and isotropic noncommutative cosmological models induced by a symplectic formalism coupled to phantom perfect fluids and a cosmological constant. After computing the field equations, we solve them to find the scalar factor dynamics. We restrict our attention to expansive solutions, that may represent the present expansion of our Universe. Those solutions generate {\it big rip} singularities. We study how the parameters, of those models, modify the time it takes to the scalar factor expands from zero till infinity, at the {\it big rip}.


title: Rapid trial-and-error learning with simulation supports flexible tool use and physical reasoning

abstract: Many animals, and an increasing number of artificial agents, display sophisticated capabilities to perceive and manipulate objects. But human beings remain distinctive in their capacity for flexible, creative tool use -- using objects in new ways to act on the world, achieve a goal, or solve a problem. To study this type of general physical problem solving, we introduce the Virtual Tools game. In this game, people solve a large range of challenging physical puzzles in just a handful of attempts. We propose that the flexibility of human physical problem solving rests on an ability to imagine the effects of hypothesized actions, while the efficiency of human search arises from rich action priors which are updated via observations of the world. We instantiate these components in the "Sample, Simulate, Update" (SSUP) model and show that it captures human performance across 30 levels of the Virtual Tools game. More broadly, this model provides a mechanism for explaining how people condense general physical knowledge into actionable, task-specific plans to achieve flexible and efficient physical problem-solving.


title: Galaxy Clusters Selected via the Sunyaev-Zel'dovich Effect in the SPTpol 100-Square-Degree Survey

abstract: We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $\mu$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass.


title: A note on Levi-Civita functional equation

abstract: In this paper we find the solutions of the functional equation $$f(xy) = g(x)h(y) + \sum_{j=1}^n g_j(x)h_j(y), \;x,y \in M,$$ where $M$ is a monoid, $n\geq 2$, and $g_j$ (for $j=1,...,n$) are linear combinations of at least $2$ distinct nonzero multiplicative functions.