title: Analysis of the Herschel DEBRIS Sun-like star sample

abstract: This paper presents a study of circumstellar debris around Sun-like stars using data from the Herschel DEBRIS Key Programme. DEBRIS is an unbiased survey comprising the nearest ~90 stars of each spectral type A-M. Analysis of the 275 F-K stars shows that excess emission from a debris disc was detected around 47 stars, giving a detection rate of 17.1+2.6-2.3 per cent, with lower rates for later spectral types. For each target a blackbody spectrum was fitted to the dust emission to determine its fractional luminosity and temperature. The derived under- lying distribution of fractional luminosity versus blackbody radius in the population showed that most detected discs are concentrated at f ~ 10^-5 and at temperatures corresponding to blackbody radii 7-40 AU, which scales to ~40 AU for realistic dust properties (similar to the current Kuiper belt). Two outlying populations are also evident; five stars have exceptionally bright emission ( f > 5x10^-5), and one has unusually hot dust < 4 AU. The excess emission distributions at all wavelengths were fitted with a steady-state evolution model, showing these are compatible with all stars being born with a narrow belt that then undergoes collisional grinding. However, the model cannot explain the hot dust systems - likely originating in transient events - and bright emission systems - arising potentially from atypically massive discs or recent stirring. The emission from the present-day Kuiper belt is predicted to be close to the median of the population, suggesting that half of stars have either depleted their Kuiper belts (similar to the Solar System), or had a lower planetesimal formation efficiency.


title: Resolution Improvement of the Common Method for Presentating Arbitrary Space Curves Voxel

abstract: The task of voxel resolution for a space curve in video memory of 3D display is set. Furthermore, an approach solution of voxel resolution of arbitrary space curve, given in parametric form, is studied. Numerous numbers of intensive experiments are conducted and interesting results with significant recommendations are presented.


title: Ultrafast heat flow in heterostructures of Au nanoclusters on thin-films: atomic-disorder induced by hot electrons

abstract: We study the ultrafast structural dynamics, in response to electronic excitations, in heterostructures composed of Au$_{923}$ nanoclusters on thin-film substrates with the use of femtosecond electron diffraction. Various forms of atomic motion, such as thermal vibrations, thermal expansion and lattice disordering, manifest as distinct and quantifiable reciprocal-space observables. In photo-excited, supported nanoclusters thermal equilibration proceeds through intrinsic heat flow, between their electrons and their lattice, and extrinsic heat flow between the nanoclusters and their substrate. For an in-depth understanding of this process, we have extended the two-temperature model to the case of 0D/2D heterostructures and used it to describe energy flow among the various subsystems, to quantify interfacial coupling constants, and to elucidate the role of the optical and thermal substrate properties. When lattice heating of Au nanoclusters is dominated by intrinsic heat flow, a reversible disordering of atomic positions occurs, which is absent when heat is injected as hot substrate-phonons. The present analysis indicates that hot electrons can distort the lattice of nanoclusters, even if the lattice temperature is below the equilibrium threshold for surface pre-melting. Based on simple considerations, the effect is interpreted as activation of surface diffusion due to modifications of the potential energy surface at high electronic temperatures. We discuss the implications of such a process in structural changes during surface chemical reactions.


title: The global well-posedness and scattering for the $5$D defocusing conformal invariant NLW with radial initial data in a critical Besov space

abstract: In this paper, we obtain the global well-posedness and scattering for the radial solution to the defocusing conformal invariant nonlinear wave equation with initial data in the critical Besov space $\dot{B}^3_{1,1}\times\dot{B}^2_{1,1}(\mathbb{R}^5)$. This is the five dimensional analogue of \cite{dodson-2016}, which is the first result on the global well-posedness and scattering of the energy subcritical nonlinear wave equation without the uniform boundedness assumption on the critical Sobolev norms employed as a substitute of the missing conservation law with respect to the scaling invariance of the equation. The proof is based on exploiting the structure of the radial solution, developing the Strichartz-type estimates and incorporation of the strategy in \cite{dodson-2016}, where we also avoid a logarithm-type loss by employing the inhomogeneous Strichartz estimates.


title: Left-orderablity for surgeries on $(-2,3,2s+1)$-pretzel knots

abstract: In this paper, we prove that the fundamental group of the manifold obtained by Dehn surgery along a $(-2,3,2s+1)$-pretzel knot ($s\ge 3$) with slope $\frac{p}{q}$ is not left orderable if $\frac{p}{q}\ge 2s+3$, and that it is left orderable if $\frac{p}{q}$ is in a neighborhood of zero depending on $s$.


title: Distributed Synthesis Using Accelerated ADMM

abstract: We propose a convex distributed optimization algorithm for synthesizing robust controllers for large-scale continuous time systems subject to exogenous disturbances. Given a large scale system, instead of solving the larger centralized synthesis task, we decompose the problem into a set of smaller synthesis problems for the local subsystems with a given interconnection topology. Hence, the synthesis problem is constrained to the sparsity pattern dictated by the interconnection topology. To this end, for each subsystem, we solve a local dissipation inequality and then check a small-gain like condition for the overall system. To minimize the effect of disturbances, we consider the $\mathrm{H}_\infty$ synthesis problems. We instantiate the distributed synthesis method using accelerated alternating direction method of multipliers (ADMM) with convergence rate $O(\frac{1}{k^2})$ with $k$ being the number of iterations.


title: Chaos and the Flow Capture Problem: Polluting is Easy, Cleaning is Hard

abstract: Cleaning pollution from a heterogeneous flow environment is far from simple. We consider the flow capture problem, which has flows and sinks in a heterogeneous environment, and investigate the problem of positioning pollutant capture units. We show that arrays of capture units carry a high risk of failure without accounting for environmental heterogeneity and chaos in their placement, design, and operation. Our idealized 2-dimensional models reveal salient features of the problem. Maximum capture efficiency depends on the required capture rate: long term efficiency decreases as the number of capture units increases, whereas short term efficiency increases. If efficiency is important, the capture process should begin as early as feasible. Knowledge of transport controlling flow structures offers predictability for unit placement. We demonstrate two heuristic approaches to near-optimally position capture units.


title: Torsion points on elliptic curves and tame semistable coverings

abstract: In this paper, we study tame Galois coverings of semistable models that arise from torsion points on elliptic curves. These coverings induce Galois morphisms of intersection graphs and we express the decomposition groups of the edges in terms of the reduction type of the elliptic curve. To that end, we first define the reduction type of an elliptic curve $E/K(C)$ on a subgraph of the intersection graph $\Sigma(\mathcal{C})$ of a strongly semistable model $\mathcal{C}$. In particular, we define the notions of good and multiplicative reduction on subgraphs of the intersection graph of $\Sigma(\mathcal{C})$. After this, we show that if an elliptic curve has good reduction on an edge and $\mathrm{char}(k)\nmid{N}$, then the $N$-torsion extension is unramified above that edge, as in the codimension one case. Furthermore, we prove a combinatorial version of a theorem by Serre on transvections for elliptic curves with a non-integral $j$-invariant. Our version states that the Galois representation $\rho_{\ell}:G\rightarrow{\mathrm{SL}_{2}(\mathbb{F}_{\ell})}$ of an elliptic curve with multiplicative reduction at an edge $e$ contains a transvection $\sigma\in{I_{e'/e}}$ for prime numbers $\ell$ that do not divide the slope of the normalized Laplacian of the $j$-invariant.


title: On the relation of Lie algebroids to constrained systems and their BV/BFV formulation

abstract: We observe that a system of irreducible, fiber-linear, first class constraints on T*M is equivalent to the definition of a foliation Lie algebroid over M. The BFV formulation of the constrained system is given by the Hamiltonian lift of the Vaintrob description (E[1],Q) of the Lie algebroid to its cotangent bundle T*E[1]. Affine deformations of the constraints are parametrized by the first Lie algebroid cohomology H^1_Q and lead to irreducible constraints also for much more general Lie algebroids such as Dirac structures; the modified BFV function follows by the addition of a representative of the deformation charge. Adding a Hamiltonian to the system corresponds to a metric g on M. Evolution invariance of the constraint surface introduces a connection nabla on E and one reobtains the compatibility of g with (E,rho,nabla) found previously in the literature. The covariantization of the Hamiltonian to a function on T*E[1] serves as a BFV-Hamiltonian, iff, in addition, this connection is compatible with the Lie algebroid structure, turning (E, rho, [ , ], nabla) into a Cartan-Lie algebroid. The BV formulation of the system is obtained from BFV by a (time-dependent) AKSZ procedure.


title: Network Utility Maximization with Heterogeneous Traffic Flows

abstract: We consider the Network Utility Maximization (NUM) problem for wireless networks in the presence of arbitrary types of flows, including unicast, broadcast, multicast, and anycast traffic. Building upon the recent framework of a universal control policy (UMW), we design a utility optimal cross-layer admission control, routing and scheduling policy, called UMW+. The UMW+ policy takes packet level actions based on a precedence-relaxed virtual network. Using Lyapunov optimization techniques, we show that UMW+ maximizes network utility, while simultaneously keeping the physical queues in the network stable. Extensive simulation results validate the performance of UMW+; demonstrating both optimal utility performance and bounded average queue occupancy. Moreover, we establish a precise one-to-one correspondence between the dynamics of the virtual queues under the UMW+ policy, and the dynamics of the dual variables of an associated offline NUM program, under a subgradient algorithm. This correspondence sheds further insight into our understanding of UMW+.


title: GraphVar 2.0: A user-friendly toolbox for machine learning on functional connectivity measures

abstract: Background: We previously presented GraphVar as a user-friendly MATLAB toolbox for comprehensive graph analyses of functional brain connectivity. Here we introduce a comprehensive extension of the toolbox allowing users to seamlessly explore easily customizable decoding models across functional connectivity measures as well as additional features. New Method: GraphVar 2.0 provides machine learning (ML) model construction, validation and exploration. Machine learning can be performed across any combination of network measures and additional variables, allowing for a flexibility in neuroimaging applications. Results: In addition to previously integrated functionalities, such as network construction and graph-theoretical analyses of brain connectivity with a high-speed general linear model (GLM), users can now perform customizable ML across connectivity matrices, network metrics and additionally imported variables. The new extension also provides parametric and nonparametric testing of classifier and regressor performance, data export, figure generation and high quality export. Comparison with existing methods: Compared to other existing toolboxes, GraphVar 2.0 offers (1) comprehensive customization, (2) an all-in-one user friendly interface, (3) customizable model design and manual hyperparameter entry, (4) interactive results exploration and data export, (5) automated cueing for modelling multiple outcome variables within the same session, (6) an easy to follow introductory review. Conclusions: GraphVar 2.0 allows comprehensive, user-friendly exploration of encoding (GLM) and decoding (ML) modelling approaches on functional connectivity measures making big data neuroscience readily accessible to a broader audience of neuroimaging investigators.


title: The [OIII] profiles of infrared-selected active galactic nuclei: More powerful outflows in the obscured population

abstract: We explore the kinematics of ionized gas via the [O III] $\lambda$5007 emission lines in active galactic nuclei (AGN) selected on the basis of their mid-infrared (IR) emission, and split into obscured and unobscured populations based on their optical-IR colors. After correcting for differences in redshift distributions, we provide composite spectra of spectroscopically and photometrically defined obscured/Type 2 and unobscured/Type 1 AGN from 3500 to 7000 \AA. The IR-selected obscured sources contain a mixture of narrow-lined Type 2 AGN and intermediate sources that have broad H$\alpha$ emission and significantly narrower H$\beta$. Using both [OIII] luminosities and AGN luminosities derived from optical-IR spectral energy distribution fitting, we find evidence for enhanced large-scale obscuration in the obscured sources. In matched bins of luminosity we find that the obscured population typically has broader, more blueshifted \OIII\ emission than in the unobscured sample, suggestive of more powerful AGN-driven outflows. This trend is not seen in spectroscopically classified samples, and is unlikely to be entirely explained by orientation effects. In addition, outflow velocities increase from small to moderate AGN $E(B-V)$ values, before flattening out (as traced by FWHM) and even decreasing (as traced by blueshift). While difficult to fully interpret in a single physical model, due to both the averaging over populations and the spatially-averaged spectra, these results agree with previous findings that simple geometric unification models are insufficient for the IR-selected AGN population, and may fit into an evolutionary model for obscured and unobscured AGN.


title: The Geometry of the Secant Caustic of a Planar Curve

abstract: The secant caustic of a planar curve $M$ is the image of the singular set of the secant map of $M$. We analyse the geometrical properties of the secant caustic of a planar curve, i.e. the number of branches of the secant caustic, the parity of the number of cusps and the number of inflexion points in each branch of this set. In particular, we investigate in detail some of the geometrical properties of the secant caustic of a rosette, i.e. a smooth regular oriented closed curve with non-vanishing curvature.


title: Chinese Text in the Wild

abstract: We introduce Chinese Text in the Wild, a very large dataset of Chinese text in street view images. While optical character recognition (OCR) in document images is well studied and many commercial tools are available, detection and recognition of text in natural images is still a challenging problem, especially for more complicated character sets such as Chinese text. Lack of training data has always been a problem, especially for deep learning methods which require massive training data. In this paper we provide details of a newly created dataset of Chinese text with about 1 million Chinese characters annotated by experts in over 30 thousand street view images. This is a challenging dataset with good diversity. It contains planar text, raised text, text in cities, text in rural areas, text under poor illumination, distant text, partially occluded text, etc. For each character in the dataset, the annotation includes its underlying character, its bounding box, and 6 attributes. The attributes indicate whether it has complex background, whether it is raised, whether it is handwritten or printed, etc. The large size and diversity of this dataset make it suitable for training robust neural networks for various tasks, particularly detection and recognition. We give baseline results using several state-of-the-art networks, including AlexNet, OverFeat, Google Inception and ResNet for character recognition, and YOLOv2 for character detection in images. Overall Google Inception has the best performance on recognition with 80.5% top-1 accuracy, while YOLOv2 achieves an mAP of 71.0% on detection. Dataset, source code and trained models will all be publicly available on the website.


title: From CCR to Levy Processes: An Excursion in Quantum Probability

abstract: This is an expositary article telling a short story made from the leaves of quantum probability with the following ingredients: (i) A special projective, unitary, irreducible and factorizable representation of the euclidean group of a Hilbert space known as the Weyl representation. \item The infinitesimal version of the Weyl representation includes the Heisenberg canonical commutation relations (CCR) of quantum theory. It also yields the three fundamental operator fields known as the creation, conservation and annihilation fields. (ii) The three fundamental fields, with the inclusion of time, lead to quantum stochastic integration and a calculus with an Ito's formula for products of differentials. (iii) Appropriate linear combinations of the fundamental operator processes yield all the L{\'e}vy processes of classical probability theory along with the bonus of Ito's formula for products of their differentials.


title: A Literature Survey on Ontology of Different Computing Platforms in Smart Environments

abstract: Smart environments integrates various types of technologies, including cloud computing, fog computing, and the IoT paradigm. In such environments, it is essential to organize and manage efficiently the broad and complex set of heterogeneous resources. For this reason, resources classification and categorization becomes a vital issue in the control system. In this paper we make an exhaustive literature survey about the various computing systems and architectures which defines any type of ontology in the context of smart environments, considering both, authors that explicitly propose resources categorization and authors that implicitly propose some resources classification as part of their system architecture. As part of this research survey, we have built a table that summarizes all research works considered, and which provides a compact and graphical snapshot of the current classification trends. The goal and primary motivation of this literature survey has been to understand the current state of the art and identify the gaps between the different computing paradigms involved in smart environment scenarios. As a result, we have found that it is essential to consider together several computing paradigms and technologies, and that there is not, yet, any research work that integrates a merged resources classification, taxonomy or ontology required in such heterogeneous scenarios.


title: The accuracy of semi-numerical reionization models in comparison with radiative transfer simulations

abstract: We have developed a modular semi-numerical code that computes the time and spatially dependent ionization of neutral hydrogen (HI), neutral (HeI) and singly ionized helium (HeII) in the intergalactic medium (IGM). The model accounts for recombinations and provides different descriptions for the photoionization rate that are used to calculate the residual HI fraction in ionized regions. We compare different semi-numerical reionization schemes to a radiative transfer (RT) simulation. We use the RT simulation as a benchmark, and find that the semi-numerical approaches produce similar HII and HeII morphologies and power spectra of the HI 21cm signal throughout reionization. As we do not track partial ionization of HeII, the extent of the double ionized helium (HeIII) regions is consistently smaller. In contrast to previous comparison projects, the ionizing emissivity in our semi-numerical scheme is not adjusted to reproduce the redshift evolution of the RT simulation, but directly derived from the RT simulation spectra. Among schemes that identify the ionized regions by the ratio of the number of ionization and absorption events on different spatial smoothing scales, we find those that mark the entire sphere as ionized when the ionization criterion is fulfilled to result in significantly accelerated reionization compared to the RT simulation. Conversely, those that flag only the central cell as ionized yield very similar but slightly delayed redshift evolution of reionization, with up to 20% ionizing photons lost. Despite the overall agreement with the RT simulation, our results suggests that constraining ionizing emissivity sensitive parameters from semi-numerical galaxy formation-reionization models are subject to photon nonconservation.


title: Coarse-grained dynamics of operator and state entanglement

abstract: We give a detailed theory for the leading coarse-grained dynamics of entanglement entropy of states and of operators in generic short-range interacting quantum many-body systems. This includes operators spreading under Heisenberg time evolution, which we find are much less entangled than "typical" operators of the same spatial support. Extending previous conjectures based on random circuit dynamics, we provide evidence that the leading-order entanglement dynamics of a given chaotic system are determined by a function $\mathcal{E}(v)$, which is model-dependent, but which we argue satisfies certain general constraints. In a minimal membrane picture, $\mathcal{E}(v)$ is the "surface tension" of the membrane and is a function of the membrane's orientation $v$ in spacetime. For one-dimensional (1D) systems this surface tension is related by a Legendre transformation to an entanglement entropy growth rate $\Gamma(\partial S/\partial x)$ which depends on the spatial "gradient" of the entanglement entropy $S(x,t)$ across the cut at position $x$. We show how to extract the entanglement growth functions numerically in 1D at infinite temperature using the concept of the operator entanglement of the time evolution operator, and we discuss possible universality of $\mathcal{E}$ at low temperatures. Our theoretical ideas are tested against and informed by numerical results for a quantum-chaotic 1D spin Hamiltonian. These results are relevant to the broad class of chaotic many-particle systems or field theories with spatially local interactions, both in 1D and above.

categories: cond-mat.stat-mech cond-mat.str-el hep-th nlin.CD quant-ph


title: Intermediate-line Emission in AGNs: The Effect of Prescription of the Gas Density

abstract: The requirement of intermediate line component in the recently observed spectra of several AGNs points to possibility of the existence of a physically separate region between broad line region (BLR) and narrow line region (NLR). In this paper we explore the emission from intermediate line region (ILR) by using the photoionization simulations of the gas clouds distributed radially from the AGN center. The gas clouds span distances typical for BLR, ILR and NLR, and the appearance of dust at the sublimation radius is fully taken into account in our model. Single cloud structure is calculated under the assumption of the constant pressure. We show that the slope of the power law cloud density radial profile does not affect the existence of ILR in major types of AGN. We found that the low ionization iron line, Fe~II, appears to be highly sensitive for the presence of dust and therefore becomes potential tracer of dust content in line emitting regions. We show that the use of disk-like cloud density profile computed at the upper part of the accretion disc atmosphere reproduces the observed properties of the line emissivities. In particular, the distance of H${\beta}$ line inferred from our model agrees with that obtained from the reverberation mapping studies in Sy1 galaxy NGC 5548.


title: Verification of Markov Decision Processes with Risk-Sensitive Measures

abstract: We develop a method for computing policies in Markov decision processes with risk-sensitive measures subject to temporal logic constraints. Specifically, we use a particular risk-sensitive measure from cumulative prospect theory, which has been previously adopted in psychology and economics. The nonlinear transformation of the probabilities and utility functions yields a nonlinear programming problem, which makes computation of optimal policies typically challenging. We show that this nonlinear weighting function can be accurately approximated by the difference of two convex functions. This observation enables efficient policy computation using convex-concave programming. We demonstrate the effectiveness of the approach on several scenarios.


title: NETT: Solving Inverse Problems with Deep Neural Networks

abstract: Recovering a function or high-dimensional parameter vector from indirect measurements is a central task in various scientific areas. Several methods for solving such inverse problems are well developed and well understood. Recently, novel algorithms using deep learning and neural networks for inverse problems appeared. While still in their infancy, these techniques show astonishing performance for applications like low-dose CT or various sparse data problems. However, there are few theoretical results for deep learning in inverse problems. In this paper, we establish a complete convergence analysis for the proposed NETT (Network Tikhonov) approach to inverse problems. NETT considers data consistent solutions having small value of a regularizer defined by a trained neural network. We derive well-posedness results and quantitative error estimates, and propose a possible strategy for training the regularizer. Our theoretical results and framework are different from any previous work using neural networks for solving inverse problems. A possible data driven regularizer is proposed. Numerical results are presented for a tomographic sparse data problem, which demonstrate good performance of NETT even for unknowns of different type from the training data. To derive the convergence and convergence rates results we introduce a new framework based on the absolute Bregman distance generalizing the standard Bregman distance from the convex to the non-convex case.


title: Divergent on average directions of Teichmuller geodesic flow

abstract: The set of directions from a quadratic differential that diverge on average under Teichmuller geodesic flow has Hausdorff dimension exactly equal to one-half.


title: Neural Networks Should Be Wide Enough to Learn Disconnected Decision Regions

abstract: In the recent literature the important role of depth in deep learning has been emphasized. In this paper we argue that sufficient width of a feedforward network is equally important by answering the simple question under which conditions the decision regions of a neural network are connected. It turns out that for a class of activation functions including leaky ReLU, neural networks having a pyramidal structure, that is no layer has more hidden units than the input dimension, produce necessarily connected decision regions. This implies that a sufficiently wide hidden layer is necessary to guarantee that the network can produce disconnected decision regions. We discuss the implications of this result for the construction of neural networks, in particular the relation to the problem of adversarial manipulation of classifiers.


title: A computationally universal phase of quantum matter

abstract: We provide the first example of a symmetry protected quantum phase that has universal computational power. Throughout this phase, which lives in spatial dimension two, the ground state is a universal resource for measurement based quantum computation.


title: Synthetic Control Methods and Big Data

abstract: Many macroeconomic policy questions may be assessed in a case study framework, where the time series of a treated unit is compared to a counterfactual constructed from a large pool of control units. I provide a general framework for this setting, tailored to predict the counterfactual by minimizing a tradeoff between underfitting (bias) and overfitting (variance). The framework nests recently proposed structural and reduced form machine learning approaches as special cases. Furthermore, difference-in-differences with matching and the original synthetic control are restrictive cases of the framework, in general not minimizing the bias-variance objective. Using simulation studies I find that machine learning methods outperform traditional methods when the number of potential controls is large or the treated unit is substantially different from the controls. Equipped with a toolbox of approaches, I revisit a study on the effect of economic liberalisation on economic growth. I find effects for several countries where no effect was found in the original study. Furthermore, I inspect how a systematically important bank respond to increasing capital requirements by using a large pool of banks to estimate the counterfactual. Finally, I assess the effect of a changing product price on product sales using a novel scanner dataset.


title: Intelligent Irrigation System Based on Arduino

abstract: This paper explains how to build an intelligent irrigation system using Arduino (a micro- controller) and many devices (humidity, temperature, pressure and water flow sensors). Our irrigation system combines a precise method to determine water balance of soils with an automatic response to water content oscillations. Thus, it is an example of how we can perform better irrigation systems by increasing the precision of measurements but also by automating decisions.


title: A general measure of the impact of priors in Bayesian statistics via Stein's Method

abstract: We propose a measure of the impact of any two choices of prior distributions by quantifying the Wasserstein distance between the respective resulting posterior distributions at any fixed sample size. We illustrate this measure on the normal, Binomial and Poisson models.


title: The Difficulty of Monte Carlo Approximation of Multivariate Monotone Functions

abstract: We study the $L_1$-approximation of $d$-variate monotone functions based on information from $n$ function evaluations. It is known that this problem suffers from the curse of dimensionality in the deterministic setting, that is, the number $n(\varepsilon,d)$ of function evaluations needed in order to approximate an unknown monotone function within a given error threshold $\varepsilon$ grows at least exponentially in $d$. This is not the case in the randomized setting (Monte Carlo setting) where the complexity $n(\varepsilon,d)$ grows exponentially in $\sqrt{d}$ (modulo logarithmic terms) only. An algorithm exhibiting this complexity is presented. Still, the problem remains difficult as best known methods are deterministic if $\varepsilon$ is comparably small, namely $\varepsilon \preceq 1/\sqrt{d}$. This inherent difficulty is confirmed by lower complexity bounds which reveal a joint $(\varepsilon,d)$-dependency and from which we deduce that the problem is not weakly tractable.


title: The TT* deformation at large central charge

abstract: We study Zamolodchikov's TT* deformation of two dimensional quantum field theories in a 't Hooft-like limit, in which we scale the number of degrees of freedom c to infinity and the deformation parameter t to zero, keeping their product t*c fixed (more precisely, we keep energies and distances fixed in units of t*c). In this limit the Hagedorn temperature remains fixed, but other non-local aspects of the theory disappear. We show that in this limit correlation functions may be computed exactly, and they are local in space and polynomials in t. We compute explicitly the deformed three-point functions of the energy-momentum tensor for a TT*-deformed conformal field theory.


title: Model-Based Value Estimation for Efficient Model-Free Reinforcement Learning

abstract: Recent model-free reinforcement learning algorithms have proposed incorporating learned dynamics models as a source of additional data with the intention of reducing sample complexity. Such methods hold the promise of incorporating imagined data coupled with a notion of model uncertainty to accelerate the learning of continuous control tasks. Unfortunately, they rely on heuristics that limit usage of the dynamics model. We present model-based value expansion, which controls for uncertainty in the model by only allowing imagination to fixed depth. By enabling wider use of learned dynamics models within a model-free reinforcement learning algorithm, we improve value estimation, which, in turn, reduces the sample complexity of learning.


title: Fault-tolerant detection of a quantum error

abstract: A critical component of any quantum error-correcting scheme is detection of errors by using an ancilla system. However, errors occurring in the ancilla can propagate onto the logical qubit, irreversibly corrupting the encoded information. We demonstrate a fault-tolerant error-detection scheme that suppresses spreading of ancilla errors by a factor of 5, while maintaining the assignment fidelity. The same method is used to prevent propagation of ancilla excitations, increasing the logical qubit dephasing time by an order of magnitude. Our approach is hardware-efficient, as it uses a single multilevel transmon ancilla and a cavity-encoded logical qubit, whose interaction is engineered in situ by using an off-resonant sideband drive. The results demonstrate that hardware-efficient approaches that exploit system-specific error models can yield advances toward fault-tolerant quantum computation.


title: Nested Algebraic Bethe Ansatz in integrable models: recent results

abstract: This short note summarizes the works done in collaboration between S. Belliard (CEA, Saclay), L. Frappat (LAPTh, Annecy), S. Pakuliak (JINR, Dubna), E. Ragoucy (LAPTh, Annecy), N. Slavnov (Steklov Math. Inst., Moscow) and more recently A. Hutsalyuk (Wuppertal / Moscow) and A. Liashyk (Kiev / Moscow). It presents the construction of Bethe vectors, their scalar products and the form factors of local operator for integrable models based on the (super)algebras $gl_n$, $gl_{m|p}$ or their quantum deformations. It corresponds to two talks given by E.R. and N.S. at \textsl{Correlation functions of quantum integrable systems and beyond}, in honor of Jean-Michel Maillet for his 60's (ENS Lyon, October 2017).


title: Active Matter Class with Second-Order Transition to Quasi-Long-Range Polar Order

abstract: We introduce and study in two dimensions a new class of dry, aligning, active matter that exhibits a direct transition to orientational order, without the phase-separation phenomenology usually observed in this context. Characterized by self-propelled particles with velocity reversals and ferromagnetic alignment of polarities, systems in this class display quasi-long-range polar order with continuously-varying scaling exponents and yet a numerical study of the transition leads to conclude that it does not belong to the Berezinskii-Kosterlitz-Thouless universality class, but is best described as a standard critical point with algebraic divergence of correlations. We rationalize these findings by showing that the interplay between order and density changes the role of defects.


title: Computational International Relations: What Can Programming, Coding and Internet Research Do for the Discipline?

abstract: Computational Social Science emerged as a highly technical and popular discipline in the last few years, owing to the substantial advances in communication technology and daily production of vast quantities of personal data. As per capita data production significantly increased in the last decade, both in terms of its size, bytes, as well as its detail, heart rate monitors, Internet connected appliances, smartphones, social scientists ability to extract meaningful social, political and demographic information from digital data also increased. A vast methodological gap exists in computational international relations, or ComInt, which refers to the use of one or a combination of tools such as data mining, natural language processing, automated text analysis, web scraping, geospatial analysis and machine learning to provide larger and better organized data to test more advanced theories of IR. After providing an overview of the potentials of computational IR and how an IR scholar can establish technical proficiency in computer science, such as starting with Python, R, QGis, ArcGIS or Github, this paper will focus on some of the author's works in providing an idea for IR students on how to think about computational IR. The paper argues that computational methods transcend the methodological schism between qualitative and quantitative approaches and form a solid foundation for building truly multi method research design.


title: RF Energy Harvesting Sensor Networks for Healthcare of Animals: Opportunities and Challenges

abstract: In recent years, the radio frequency (RF) energy harvesting technique is considered as a favorable alternative to supply power for the next-generation wireless sensor networks. Due to the features of energy self-sustainability and long lifetime, the energy harvesting sensor network (RF-EHSN) becomes a promising solution to provide smart healthcare services. Nowadays, many energy harvesting based applications have been developed for monitoring the health status of human beings; how to benefit animals, however, has not yet drawn people's attention. This article explores the potential of applying RF-EHSNs to monitoring the health level of animals. The unique challenges and potential solutions at different layers of an RF-EHSN for animals' healthcare service are studied.


title: Pulling Out All the Tops with Computer Vision and Deep Learning

abstract: We apply computer vision with deep learning -- in the form of a convolutional neural network (CNN) -- to build a highly effective boosted top tagger. Previous work (the "DeepTop" tagger of Kasieczka et al) has shown that a CNN-based top tagger can achieve comparable performance to state-of-the-art conventional top taggers based on high-level inputs. Here, we introduce a number of improvements to the DeepTop tagger, including architecture, training, image preprocessing, sample size and color pixels. Our final CNN top tagger outperforms BDTs based on high-level inputs by a factor of $\sim 2$--3 or more in background rejection, over a wide range of tagging efficiencies and fiducial jet selections. As reference points, we achieve a QCD background rejection factor of 500 (60) at 50\% top tagging efficiency for fully-merged (non-merged) top jets with $p_T$ in the 800--900 GeV (350--450 GeV) range. Our CNN can also be straightforwardly extended to the classification of other types of jets, and the lessons learned here may be useful to others designing their own deep NNs for LHC applications.


title: Martingale decomposition of a $L^2$ space with nonlinear stochastic integrals

abstract: This paper presents a generalization of the Kunita-Watanabe decomposition of a $L^2$ space with nonlinear stochastic integrals where the integrator is a family of continuous martingales bounded in $L^2$. To get the result, a useful relation between the regularity of the martingale family respect to its parameter and the regularity of the integrand in its martingale decomposition is shown.The decomposition presented in the main result is also the solution of an optimization problem in $L^2$. Finally, an example is given where the optimization problem is solved explicitely.


title: Asymptotics of quantum channels: conserved quantities, an adiabatic limit, and matrix product states

abstract: This work derives an analytical formula for the asymptotic state---the quantum state resulting from an infinite number of applications of a general quantum channel on some initial state. For channels admitting multiple fixed or rotating points, conserved quantities---the left fixed/rotating points of the channel---determine the dependence of the asymptotic state on the initial state. The formula stems from a Noether-like theorem stating that, for any channel admitting a full-rank fixed point, conserved quantities commute with that channel's Kraus operators up to a phase. The formula is applied to adiabatic transport of the fixed-point space of channels, revealing cases where the dissipative/spectral gap can close during any segment of the adiabatic path. The formula is also applied to calculate expectation values of noninjective matrix product states (MPS) in the thermodynamic limit, revealing that those expectation values can also be calculated using an MPS with reduced bond dimension and a modified boundary.


title: Maximizing Activity in Ising Networks via the TAP Approximation

abstract: A wide array of complex biological, social, and physical systems have recently been shown to be quantitatively described by Ising models, which lie at the intersection of statistical physics and machine learning. Here, we study the fundamental question of how to optimize the state of a networked Ising system given a budget of external influence. In the continuous setting where one can tune the influence applied to each node, we propose a series of approximate gradient ascent algorithms based on the Plefka expansion, which generalizes the na\"{i}ve mean field and TAP approximations. In the discrete setting where one chooses a small set of influential nodes, the problem is equivalent to the famous influence maximization problem in social networks with an additional stochastic noise term. In this case, we provide sufficient conditions for when the objective is submodular, allowing a greedy algorithm to achieve an approximation ratio of $1-1/e$. Additionally, we compare the Ising-based algorithms with traditional influence maximization algorithms, demonstrating the practical importance of accurately modeling stochastic fluctuations in the system.


title: Predicting Recall Probability to Adaptively Prioritize Study

abstract: Students have a limited time to study and are typically ineffective at allocating study time. Machine-directed study strategies that identify which items need reinforcement and dictate the spacing of repetition have been shown to help students optimize mastery (Mozer & Lindsey 2017). The large volume of research on this matter is typically conducted in constructed experimental settings with fixed instruction, content, and scheduling; in contrast, we aim to develop methods that can address any demographic, subject matter, or study schedule. We show two methods that model item-specific recall probability for use in a discrepancy-reduction instruction strategy. The first model predicts item recall probability using a multiple logistic regression (MLR) model based on previous answer correctness and temporal spacing of study. Prompted by literature suggesting that forgetting is better modeled by the power law than an exponential decay (Wickelgren 1974), we compare the MLR approach with a Recurrent Power Law (RPL) model which adaptively fits a forgetting curve. We then discuss the performance of these models against study datasets comprised of millions of answers and show that the RPL approach is more accurate and flexible than the MLR model. Finally, we give an overview of promising future approaches to knowledge modeling.


title: Cooperation dynamics of generalized reciprocity in state-based social dilemmas

abstract: We introduce a framework for studying social dilemmas in networked societies where individuals follow a simple state-based behavioral mechanism based on generalized reciprocity, which is rooted in the principle "help anyone if helped by someone". Within this general framework, which applies to a wide range of social dilemmas including, among others, public goods, donation and snowdrift games, we study the cooperation dynamics on a variety of complex network examples. By interpreting the studied model through the lenses of nonlinear dynamical systems, we show that cooperation through generalized reciprocity always emerges as the unique attractor in which the overall level of cooperation is maximized, while simultaneously exploitation of the participating individuals is prevented. The analysis elucidates the role of the network structure, here captured by a local centrality measure which uniquely quantifies the propensity of the network structure to cooperation, by dictating the degree of cooperation displayed both at microscopic and macroscopic level. We demonstrate the applicability of the analysis on a practical example by considering an interaction structure that couples a donation process with a public goods game.


title: Approximate Inference for Constructing Astronomical Catalogs from Images

abstract: We present a new, fully generative model for constructing astronomical catalogs from optical telescope image sets. Each pixel intensity is treated as a random variable with parameters that depend on the latent properties of stars and galaxies. These latent properties are themselves modeled as random. We compare two procedures for posterior inference. One procedure is based on Markov chain Monte Carlo (MCMC) while the other is based on variational inference (VI). The MCMC procedure excels at quantifying uncertainty, while the VI procedure is 1000 times faster. On a supercomputer, the VI procedure efficiently uses 665,000 CPU cores to construct an astronomical catalog from 50 terabytes of images in 14.6 minutes, demonstrating the scaling characteristics necessary to construct catalogs for upcoming astronomical surveys.


title: SQL-Rank: A Listwise Approach to Collaborative Ranking

abstract: In this paper, we propose a listwise approach for constructing user-specific rankings in recommendation systems in a collaborative fashion. We contrast the listwise approach to previous pointwise and pairwise approaches, which are based on treating either each rating or each pairwise comparison as an independent instance respectively. By extending the work of (Cao et al. 2007), we cast listwise collaborative ranking as maximum likelihood under a permutation model which applies probability mass to permutations based on a low rank latent score matrix. We present a novel algorithm called SQL-Rank, which can accommodate ties and missing data and can run in linear time. We develop a theoretical framework for analyzing listwise ranking methods based on a novel representation theory for the permutation model. Applying this framework to collaborative ranking, we derive asymptotic statistical rates as the number of users and items grow together. We conclude by demonstrating that our SQL-Rank method often outperforms current state-of-the-art algorithms for implicit feedback such as Weighted-MF and BPR and achieve favorable results when compared to explicit feedback algorithms such as matrix factorization and collaborative ranking.


title: Holomorphic quadratic differentials on graphs and the chromatic polynomial

abstract: We study "holomorphic quadratic differentials" on graphs. We relate them to the reactive power in an LC circuit, and also to the chromatic polynomial of a graph. Specifically, we show that the chromatic polynomial $\chi$ of a graph $G$, at negative integer values, can be evaluated as the degree of a certain rational mapping, arising from the defining equations for a holomorphic quadratic differential. This allows us to give an explicit integral expression for $\chi(-k)$.


title: Separators and Adjustment Sets in Causal Graphs: Complete Criteria and an Algorithmic Framework

abstract: Principled reasoning about the identifiability of causal effects from non-experimental data is an important application of graphical causal models. This paper focuses on effects that are identifiable by covariate adjustment, a commonly used estimation approach. We present an algorithmic framework for efficiently testing, constructing, and enumerating $m$-separators in ancestral graphs (AGs), a class of graphical causal models that can represent uncertainty about the presence of latent confounders. Furthermore, we prove a reduction from causal effect identification by covariate adjustment to $m$-separation in a subgraph for directed acyclic graphs (DAGs) and maximal ancestral graphs (MAGs). Jointly, these results yield constructive criteria that characterize all adjustment sets as well as all minimal and minimum adjustment sets for identification of a desired causal effect with multivariate exposures and outcomes in the presence of latent confounding. Our results extend several existing solutions for special cases of these problems. Our efficient algorithms allowed us to empirically quantify the identifiability gap between covariate adjustment and the do-calculus in random DAGs and MAGs, covering a wide range of scenarios. Implementations of our algorithms are provided in the R package dagitty.


title: Redundancy allocation in finite-length nested codes for nonvolatile memories

abstract: In this paper, we investigate the optimum way to allocate redundancy of finite-length nested codes for modern nonvolatile memories suffering from both permanent defects and transient errors (erasures or random errors). A nested coding approach such as partitioned codes can handle both permanent defects and transient errors by using two parts of redundancy: 1) redundancy to deal with permanent defects and 2) redundancy for transient errors. We consider two different channel models of the binary defect and erasure channel (BDEC) and the binary defect and symmetric channel (BDSC). The transient errors of the BDEC are erasures and the BDSC's transient errors are modeled by the binary symmetric channel, respectively. Asymptotically, the probability of recovery failure can converge to zero if the capacity region conditions of nested codes are satisfied. However, the probability of recovery failure of finite-length nested codes can be significantly variable for different redundancy allocations even though they all satisfy the capacity region conditions. Hence, we formulate the redundancy allocation problem of finite-length nested codes to minimize the recovery failure probability. We derive the upper bounds on the probability of recovery failure and use them to estimate the optimal redundancy allocation. Numerical results show that our estimated redundancy allocation matches well the optimal redundancy allocation.


title: Surges of collective human activity emerge from simple pairwise correlations

abstract: Human populations exhibit complex behaviors---characterized by long-range correlations and surges in activity---across a range of social, political, and technological contexts. Yet it remains unclear where these collective behaviors come from, or if there even exists a set of unifying principles. Indeed, existing explanations typically rely on context-specific mechanisms, such as traffic jams driven by work schedules or spikes in online traffic induced by significant events. However, analogies with statistical mechanics suggest a more general mechanism: that collective patterns can emerge organically from fine-scale interactions within a population. Here, across four different modes of human activity, we show that the simplest correlations in a population---those between pairs of individuals---can yield accurate quantitative predictions for the large-scale behavior of the entire population. To quantify the minimal consequences of pairwise correlations, we employ the principle of maximum entropy, making our description equivalent to an Ising model whose interactions and external fields are notably calculated from past observations of population activity. In addition to providing accurate quantitative predictions, we show that the topology of learned Ising interactions resembles the network of inter-human communication within a population. Together, these results demonstrate that fine-scale correlations can be used to predict large-scale social behaviors, a perspective that has critical implications for modeling and resource allocation in human populations.


title: Integrating Human-Provided Information Into Belief State Representation Using Dynamic Factorization

abstract: In partially observed environments, it can be useful for a human to provide the robot with declarative information that represents probabilistic relational constraints on properties of objects in the world, augmenting the robot's sensory observations. For instance, a robot tasked with a search-and-rescue mission may be informed by the human that two victims are probably in the same room. An important question arises: how should we represent the robot's internal knowledge so that this information is correctly processed and combined with raw sensory information? In this paper, we provide an efficient belief state representation that dynamically selects an appropriate factoring, combining aspects of the belief when they are correlated through information and separating them when they are not. This strategy works in open domains, in which the set of possible objects is not known in advance, and provides significant improvements in inference time over a static factoring, leading to more efficient planning for complex partially observed tasks. We validate our approach experimentally in two open-domain planning problems: a 2D discrete gridworld task and a 3D continuous cooking task. A supplementary video can be found at http://tinyurl.com/chitnis-iros-18.


title: Simplified Weak Galerkin and New Finite Difference Schemes for the Stokes Equation

abstract: This article presents a simplified formulation for the weak Galerkin finite element method for the Stokes equation without using the degrees of freedom associated with the unknowns in the interior of each element as formulated in the original weak Galerkin algorithm. The simplified formulation preserves the important mass conservation property locally on each element and allows the use of general polygonal partitions. A particular application of the simplified weak Galerkin on rectangular partitions yields a new class of 5- and 7-point finite difference schemes for the Stokes equation. An explicit formula is presented for the computation of the element stiffness matrices on arbitrary polygonal elements. Error estimates of optimal order are established for the simplified weak Galerkin finite element method in the H^1 and L^2 norms. Furthermore, a superconvergence of order O(h^{1.5}) is established on rectangular partitions for the velocity approximation in the H^1 norm at cell centers, and a similar superconvergence is derived for the pressure approximation in the L^2 norm at cell centers. Some numerical results are reported to confirm the convergence and superconvergence theory.


title: Symmetry in optics and photonics: a group theory approach

abstract: Group theory (GT) provides a rigorous framework for studying symmetries in various disciplines in physics ranging from quantum field theories and the standard model to fluid mechanics and chaos theory. To date, the application of such a powerful tool in optical physics remains limited. Over the past few years however, several quantum-inspired symmetry principles (such as parity-time invariance and supersymmetry) have been introduced in optics and photonics for the first time. Despite the intense activities in these new research directions, only few works utilized the power of group theory. Motivated by this status quo, here we present a brief overview of the application of GT in optics, deliberately choosing examples that illustrate the power of this tool in both continuous and discrete setups. We hope that this review will stimulate further research that exploits the full potential of GT for investigating various symmetry paradigms in optics, eventually leading to new photonic devices.


title: Geometric random graphs and Rado sets of continuous functions

abstract: We prove the existence of Rado sets in the Banach space of continuous functions on [0,1]. A countable dense set S is Rado if with probability 1, the infinite geometric random graph on S, formed by probabilistically making adjacent elements of S that are within unit distance of each other, is unique up to isomorphism. We show that for a suitable measure which we construct, almost all countable dense sets in the subspaces of piecewise linear functions and of polynomials are Rado. Moreover, all graphs arising from such sets are of a unique isomorphism type. For the subspace of Brownian motion paths, almost all countable subsets are Rado (for a suitable measure) and the resulting graphs are of a unique isomorphism type. We show that the graph arising from piecewise linear functions and polynomials is not isomorphic to the graph arising from Brownian motion paths. Moreover, these graphs are non-isomorphic to graphs arising from Rado sets in $\mathbb{R}^n$, or the sequence spaces $c$ and $c_0$.


title: On generalized Walsh bases

abstract: This paper continues the study of orthonormal bases (ONB) of $L^2[0,1]$ introduced in \cite{DPS14} by means of Cuntz algebra $\mathcal{O}_N$ representations on $L^2[0,1]$. For $N=2$, one obtains the classic Walsh system. We show that the ONB property holds precisely because the $\mathcal{O}_N$ representations are irreducible. We prove an uncertainty principle related to these bases. As an application to discrete signal processing we find a fast generalized transform and compare this generalized transform with the classic one with respect to compression and sparse signal recovery.


title: Improving Sentiment Analysis in Arabic Using Word Representation

abstract: The complexities of Arabic language in morphology, orthography and dialects makes sentiment analysis for Arabic more challenging. Also, text feature extraction from short messages like tweets, in order to gauge the sentiment, makes this task even more difficult. In recent years, deep neural networks were often employed and showed very good results in sentiment classification and natural language processing applications. Word embedding, or word distributing approach, is a current and powerful tool to capture together the closest words from a contextual text. In this paper, we describe how we construct Word2Vec models from a large Arabic corpus obtained from ten newspapers in different Arab countries. By applying different machine learning algorithms and convolutional neural networks with different text feature selections, we report improved accuracy of sentiment classification (91%-95%) on our publicly available Arabic language health sentiment dataset [1]


title: A social Network Analysis of the Operations Research/Industrial Engineering Faculty Hiring Network

abstract: We study the U.S. Operations Research/Industrial-Systems Engineering (ORIE) faculty hiring network, consisting of 1,179 faculty origin and destination data together with attribute data from 83 ORIE departments. A social network analysis of faculty hires can reveal important patterns in an academic field, such as the existence of a hierarchy or sociological aspects such as the presence of communities of departments. We first statistically test for the existence of a linear hierarchy in the network and for its steepness. We find a near linear hierarchical order of the departments, proposing a new index for hiring networks, which we contrast with other indicators of hierarchy, including published rankings. A single index is not capable to capture the full structure of a complex network, however, so we next fit a latent exponential random graph model (ERGM) to the network, which is able to reproduce its main observed characteristics: high incidence of self-hiring, skewed out-degree distribution, low density and clustering. Finally, we use the latent variables in the ERGM to simplify the network to one where faculty hires take place among three groups of departments. We contrast our findings with those reported for other related disciplines, Computer Science and Business.


title: Efficient sampling of reversible cross-linking polymers: Self-assembly of single-chain polymeric nanoparticles

abstract: We present a new simulation technique to study systems of polymers functionalized by reactive sites that bind/unbind forming reversible linkages. Functionalized polymers feature self-assembly and responsive properties that are unmatched by systems lacking selective interactions. The scales at which the functional properties of these materials emerge are difficult to model, especially in the reversible regime where such properties result from many binding/unbinding events. This difficulty is related to large entropic barriers associated with the formation of intra-molecular loops. In this work we present a simulation scheme that sidesteps configurational costs by dedicated Monte Carlo moves capable of binding/unbinding reactive sites in a single step. Cross-linking reactions are implemented by trial moves that reconstruct chain sections attempting, at the same time, a dimerization reaction between pairs of reactive sites. The model is parametrized by the reaction equilibrium constant of the reactive species free in solution. This quantity can be obtained by means of experiments or atomistic/quantum simulations. We use the proposed methodology to study self-assembly of single--chain polymeric nanoparticles, starting from flexible precursors carrying regularly or randomly distributed reactive sites. During a single run, almost all pairs of reactive monomers interact at least once. We focus on understanding differences in the morphology of chain nanoparticles when linkages are reversible as compared to the well studied case of irreversible reactions. Intriguingly, we find that the size of regularly functionalsized chains, in good solvent conditions, is non-monotonous as a function of the degree of functionalization. We clarify how this result follows from excluded volume interactions and is peculiar of reversible linkages and regular functionalizations.


title: SalientDSO: Bringing Attention to Direct Sparse Odometry

abstract: Although cluttered indoor scenes have a lot of useful high-level semantic information which can be used for mapping and localization, most Visual Odometry (VO) algorithms rely on the usage of geometric features such as points, lines and planes. Lately, driven by this idea, the joint optimization of semantic labels and obtaining odometry has gained popularity in the robotics community. The joint optimization is good for accurate results but is generally very slow. At the same time, in the vision community, direct and sparse approaches for VO have stricken the right balance between speed and accuracy. We merge the successes of these two communities and present a way to incorporate semantic information in the form of visual saliency to Direct Sparse Odometry - a highly successful direct sparse VO algorithm. We also present a framework to filter the visual saliency based on scene parsing. Our framework, SalientDSO, relies on the widely successful deep learning based approaches for visual saliency and scene parsing which drives the feature selection for obtaining highly-accurate and robust VO even in the presence of as few as 40 point features per frame. We provide extensive quantitative evaluation of SalientDSO on the ICL-NUIM and TUM monoVO datasets and show that we outperform DSO and ORB-SLAM - two very popular state-of-the-art approaches in the literature. We also collect and publicly release a CVL-UMD dataset which contains two indoor cluttered sequences on which we show qualitative evaluations. To our knowledge this is the first paper to use visual saliency and scene parsing to drive the feature selection in direct VO.


title: Real-Time Detection of Hybrid and Stealthy Cyber-Attacks in Smart Grid

abstract: For a safe and reliable operation of the smart grid, timely detection of cyber-attacks is of critical importance. Moreover, considering smarter and more capable attackers, robust detection mechanisms are needed against a diverse range of cyber-attacks. With these purposes, we propose a robust online detection algorithm for (possibly combined) false data injection (FDI) and jamming attacks, that also provides online estimates of the unknown and time-varying attack parameters and recovered state estimates. Further, considering smarter attackers that are capable of designing stealthy attacks to prevent the detection or to increase the detection delay of the proposed algorithm, we propose additional countermeasures. Numerical studies illustrate the quick and reliable response of the proposed detection mechanisms against hybrid and stealthy cyber-attacks.


title: Modal approach to the controllability problem of distributed parameter systems with damping

abstract: This paper is devoted to the controllability analysis of a class of linear control systems in a Hilbert space. It is proposed to use the minimum energy controls of a reduced lumped parameter system for solving the infinite dimensional steering problem approximately. Sufficient conditions of the approximate controllability are formulated for a modal representation of a flexible structure with small damping.


title: Ring loss: Convex Feature Normalization for Face Recognition

abstract: We motivate and present Ring loss, a simple and elegant feature normalization approach for deep networks designed to augment standard loss functions such as Softmax. We argue that deep feature normalization is an important aspect of supervised classification problems where we require the model to represent each class in a multi-class problem equally well. The direct approach to feature normalization through the hard normalization operation results in a non-convex formulation. Instead, Ring loss applies soft normalization, where it gradually learns to constrain the norm to the scaled unit circle while preserving convexity leading to more robust features. We apply Ring loss to large-scale face recognition problems and present results on LFW, the challenging protocols of IJB-A Janus, Janus CS3 (a superset of IJB-A Janus), Celebrity Frontal-Profile (CFP) and MegaFace with 1 million distractors. Ring loss outperforms strong baselines, matches state-of-the-art performance on IJB-A Janus and outperforms all other results on the challenging Janus CS3 thereby achieving state-of-the-art. We also outperform strong baselines in handling extremely low resolution face matching.


title: Male pelvic synthetic CT generation from T1-weighted MRI using 2D and 3D convolutional neural networks

abstract: To achieve magnetic resonance (MR)-only radiotherapy, a method needs to be employed to estimate a synthetic CT (sCT) for generating electron density maps and patient positioning reference images. We investigated 2D and 3D convolutional neural network (CNN) methods to generate a male pelvic sCT using a T1-weighted MR image. A retrospective study was performed using CTs and T1-weighted MR images of 20 prostate cancer patients. The proposed 2D CNN model, which contained 27 convolutional layers, was modified from the SegNet for better performance. 3D version of the CNN model was also developed. Both CNN models were trained from scratch to map intensities of T1-weighted MR images to CT Hounsfield Unit (HU) values. Each sCT was generated in a five-fold-cross-validation framework and compared with the corresponding CT using voxel-wise mean absolute error (MAE), and dice similarity coefficient (DSC), recall, and precision for bony structures. Wilcoxon signed-rank tests were performed to evaluate the differences between the both models. The MAE averaged across all patients were 40.5 $\pm$ 5.4 HU and 37.6 $\pm$ 5.1 HU for the 2D and 3D CNN models, respectively. The DSC, recall, and precision of the bony structures were 0.81 $\pm$ 0.04, 0.85 $\pm$ 0.04, and 0.77 $\pm$ 0.09 for the 2D CNN model, and 0.82 $\pm$ 0.04, 0.84 $\pm$ 0.04, and 0.80 $\pm$ 0.08 for the 3D CNN model, respectively. P values of the Wilcoxon signed-rank tests were less than 0.05 except for recall, which was 0.6. The 2D and 3D CNN models generated accurate pelvic sCTs for the 20 patients using T1-weighted MR images. The evaluation metrics and statistical tests indicated that the 3D model was able to generate sCTs with better MAE, bone DSC, and bone precision. The accuracy of the dose calculation and patient positioning using generated sCTs will be tested and compared for the two models in the future.


title: Development of GEM Detectors at Hampton University

abstract: Two GEM telescopes, each consisting of three 10x10 cm$^2$ triple-GEM chambers were built, tested and operated by the Hampton University group. The GEMs are read out with APV25 frontend chips and FPGA based digitizing electronics developed by INFN Rome. The telescopes were used for the luminosity monitoring system at the OLYMPUS experiment at DESY in Germany, with positron and electron beams at 2 GeV. The GEM elements have been recycled to serve in another two applications: Three GEM elements are used to track beam particles in the MUSE experiment at PSI in Switzerland. A set of four elements has been configured as a prototype tracker for phase 1a of the DarkLight experiment at the Low-Energy Recirculator Facility (LERF) at Jefferson Lab in Newport News, USA, in a first test run in summer 2016. The Hampton group is responsible for the DarkLight phase-I lepton tracker in preparation. Further efforts are ongoing to optimize the data acquisition speed for GEM operations in MUSE and DarkLight. An overview of the group's GEM detector related activities will be given.


title: Materials data validation and imputation with an artificial neural network

abstract: We apply an artificial neural network to model and verify material properties. The neural network algorithm has a unique capability to handle incomplete data sets in both training and predicting, so it can regard properties as inputs allowing it to exploit both composition-property and property-property correlations to enhance the quality of predictions, and can also handle a graphical data as a single entity. The framework is tested with different validation schemes, and then applied to materials case studies of alloys and polymers. The algorithm found twenty errors in a commercial materials database that were confirmed against primary data sources.


title: A Characterization of Boundary Representations of Positive Matrices in the Hardy Space via the Abel Product

abstract: Spectral measures give rise to a natural harmonic analysis on the unit disc via a boundary representation of a positive matrix arising from a spectrum of the measure. We consider in this paper the reverse: for a positive matrix in the Hardy space of the unit disc we consider which measures, if any, yield a boundary representation of the positive matrix. We prove a characterization of those representing measures via a matrix identity by introducing a new operator product called the Abel Product.


title: Quantum annealing versus classical machine learning applied to a simplified computational biology problem

abstract: Transcription factors regulate gene expression, but how these proteins recognize and specifically bind to their DNA targets is still debated. Machine learning models are effective means to reveal interaction mechanisms. Here we studied the ability of a quantum machine learning approach to predict binding specificity. Using simplified datasets of a small number of DNA sequences derived from actual binding affinity experiments, we trained a commercially available quantum annealer to classify and rank transcription factor binding. The results were compared to state-of-the-art classical approaches for the same simplified datasets, including simulated annealing, simulated quantum annealing, multiple linear regression, LASSO, and extreme gradient boosting. Despite technological limitations, we find a slight advantage in classification performance and nearly equal ranking performance using the quantum annealer for these fairly small training data sets. Thus, we propose that quantum annealing might be an effective method to implement machine learning for certain computational biology problems.


title: Structure of Anomalies of 4d SCFTs from M5-branes, and Anomaly Inflow

abstract: We study the 't Hooft anomalies of four-dimensional superconformal field theories that arise from M5-branes wrapped on a punctured Riemann surface. In general there are two independent contributions to the anomalies. There is a bulk term obtained by integrating the anomaly polynomial of the world-volume theory on the M5-branes over the Riemann surface; this contribution knows about the punctures only through its dependence on the Euler characteristic of the surface. The second set of contributions comes from local data at the punctures; these terms are independent from the bulk data of the surface. Using anomaly inflow in M-theory, we describe the general structure of the anomalies for cases when the four-dimensional theories preserve N=2 supersymmetry. We additionally discuss the anomalies corresponding to (p,q) punctures in N=1 theories.


title: Low-mass eclipsing binaries in the WFCAM Transit Survey: the persistence of the M-dwarf radius inflation problem

abstract: We present the characterization of 5 new short-period low-mass eclipsing binaries from the WFCAM Transit Survey. The analysis was performed by using the photometric WFCAM J-mag data and additional low- and intermediate-resolution spectroscopic data to obtain both orbital and physical properties of the studied sample. The light curves and the measured radial velocity curves were modeled simultaneously with the JKTEBOP code, with Markov chain Monte Carlo simulations for the error estimates. The best-model fit have revealed that the investigated detached binaries are in very close orbits, with orbital separations of $2.9 \leq a \leq 6.7$ $R_{\odot}$ and short periods of $0.59 \leq P_{\rm orb} \leq 1.72$ d, approximately. We have derived stellar masses between $0.24$ and $0.72$ $M_{\odot}$ and radii ranging from $0.42$ to $0.67$ $R_{\odot}$. The great majority of the LMEBs in our sample has an estimated radius far from the predicted values according to evolutionary models. The components with derived masses of $M < 0.6$ $M_{\odot}$ present a radius inflation of $\sim$$9\%$ or more. This general behavior follows the trend of inflation for partially-radiative stars proposed previously. These systems add to the increasing sample of low-mass stellar radii that are not well-reproduced by stellar models. They further highlight the need to understand the magnetic activity and physical state of small stars. Missions like TESS will provide many such systems to perform high-precision radius measurements to tightly constrain low-mass stellar evolution models.


title: A novel approach for fusion of heterogeneous sources of data

abstract: With advancements in sensor technology, a heterogeneous set of data, containing samples of scalar, waveform signal, image, or even structured point cloud are becoming increasingly popular. Developing a statistical model, representing the behavior of the underlying system based upon such a heterogeneous set of data can be used in monitoring, control, and optimization of the system. Unfortunately, available methods only focus on the scalar and curve data and do not provide a general framework that can integrate different sources of data to construct a model. This paper poses the problem of estimating a process output, measured by a scalar, curve, an image, or a point cloud by a set of heterogeneous process variables such as scalar process setting, sensor readings, and images. We introduce a general approach in which each set of input data (predictor) as well as the output measurements are represented by tensors. We formulate a linear regression model between the input and output tensors and estimate the parameters by minimizing a least square loss function. In order to avoid overfitting and to reduce the number of parameters to be estimated, we decompose the model parameters using several bases, spanning the input and output spaces. Next, we learn both the bases and their spanning coefficients when minimizing the loss function using an alternating least square (ALS) algorithm. We show that such a minimization has a closed-form solution in each iteration and can be computed very efficiently. Through several simulation and case studies, we evaluate the performance of the proposed method. The results reveal the advantage of the proposed method over some benchmarks in the literature in terms of the mean square prediction error.


title: Stochastic multi-symplectic Runge-Kutta methods for stochastic Hamiltonian PDEs

abstract: In this paper, we consider stochastic Runge-Kutta methods for stochastic Hamiltonian partial differential equations and present some sufficient conditions for multisymplecticity of stochastic Runge-Kutta methods of stochastic Hamiltonian partial differential equations. Particularly, we apply these ideas to stochastic Maxwell equations with multiplicative noise, possessing the stochastic multi-symplectic conservation law and energy conservation law. Theoretical analysis shows that the methods can preserve both the discrete stochastic multi-symplectic conservation law and discrete energy conservation law almost surely.


title: Homological approximations for Profinite and Pro-$p$ limit groups

abstract: We study homological approximations of the profinite completion of a limit group (see Thm.~A) and obtain the analogous of Bridson and Howie's Theorem for the profinite completion of a non-abelian limit group (see Thm.~B).


title: Entanglement concentration and purification of two-mode squeezed microwave photons in circuit QED

abstract: We present a theoretical proposal for a physical implementation of entanglement concentration and purification protocols for two-mode squeezed microwave photons in circuit quantum electrodynamics (QED). First, we give the description of the cross-Kerr effect induced between two resonators in circuit QED. Then we use the cross-Kerr media to design the effective quantum nondemolition (QND) measurement on microwave-photon number. By using the QND measurement, the parties in quantum communication can accomplish the entanglement concentration and purification of nonlocal two-mode squeezed microwave photons. We discuss the feasibility of our schemes by giving the detailed parameters which can be realized with current experimental technology. Our work can improve some practical applications in continuous-variable microwave-based quantum information processing.


title: Invariants of deformations of quotient surface singularities

abstract: We find all $P$-resolutions of quotient surface singularities (especially, tetrahedral, octahedral, and icosahedral singularities) together with their dual graphs, which reproduces Jan Steven's list [Manuscripta Math. 1993] of the numbers of $P$-resolutions of each singularities. We then compute the dimensions and Milnor numbers of the corresponding irreducible components of the reduced base spaces of versal deformations of each singularities. Furthermore we realize Milnor fibers as complements of certain divisors (depending only on the singularities) in rational surfaces via the minimal model program for 3-folds. Then we compare Milnor fibers with minimal symplectic fillings, where the latter are classified by Bhupal and Ono [Nagoya Math. J. 2012]. As an application, we show that there are 6 pairs of entries in the list of Bhupal and Ono [Nagoya Math. J. 2012] such that two entries in each pairs represent diffeomorphic minimal symplectic fillings.


title: The PPT square conjecture holds generically for some classes of independent states

abstract: Let $|\psi\rangle\langle \psi|$ be a random pure state on $\mathbb{C}^{d^2}\otimes \mathbb{C}^s$, where $\psi$ is a random unit vector uniformly distributed on the sphere in $\mathbb{C}^{d^2}\otimes \mathbb{C}^s$. Let $\rho_1$ be random induced states $\rho_1=Tr_{\mathbb{C}^s}(|\psi\rangle\langle \psi |)$ whose distribution is $\mu_{d^2,s}$; and let $\rho_2$ be random induced states following the same distribution $\mu_{d^2,s}$ independent from $\rho_1$. Let $\rho$ be a random state induced by the entanglement swapping of $\rho_1$ and $\rho_2$. We show that the empirical spectrum of $\rho- {1\mkern -4mu{\rm l}}/d^2$ converges almost surely to the Marcenko-Pastur law with parameter $c^2$ as $d\rightarrow \infty$ and $s/d \rightarrow c$. As an application, we prove that the state $\rho$ is separable generically if $\rho_1, \rho_2$ are PPT entangled.


title: Learning Longer-term Dependencies in RNNs with Auxiliary Losses

abstract: Despite recent advances in training recurrent neural networks (RNNs), capturing long-term dependencies in sequences remains a fundamental challenge. Most approaches use backpropagation through time (BPTT), which is difficult to scale to very long sequences. This paper proposes a simple method that improves the ability to capture long term dependencies in RNNs by adding an unsupervised auxiliary loss to the original objective. This auxiliary loss forces RNNs to either reconstruct previous events or predict next events in a sequence, making truncated backpropagation feasible for long sequences and also improving full BPTT. We evaluate our method on a variety of settings, including pixel-by-pixel image classification with sequence lengths up to 16\,000, and a real document classification benchmark. Our results highlight good performance and resource efficiency of this approach over competitive baselines, including other recurrent models and a comparable sized Transformer. Further analyses reveal beneficial effects of the auxiliary loss on optimization and regularization, as well as extreme cases where there is little to no backpropagation.


title: Critical behavior of magnetization in URhAl:Quasi-two-dimensional Ising system with long-range interactions

abstract: The critical behavior of dc magnetization in the uranium ferromagnet URhAl with the hexagonal ZrNiAl-type crystal structure has been studied around the ferromagnetic transition temperature T_C. The critical exponent beta for the temperature dependence of the spontaneous magnetization below T_C, gamma for the magnetic susceptibility, and delta for the magnetic isotherm at T_C have been obtained with a modified Arrott plot, a Kouvel-Fisher plot, the critical isotherm analysis and the scaling analysis. We have determined the critical exponents as beta = 0.287 +- 0.005, gamma = 1.47 +- 0.02, and delta = 6.08 +- 0.04 by the scaling analysis and the critical isotherm analysis. These critical exponents satisfy the Widom scaling law delta=1+gamma/beta. URhAl has strong uniaxial magnetic anisotropy, similar to its isostructural UCoAl that has been regarded as a three-dimensional (3D) Ising system in previous studies. However, the universality class of the critical phenomenon in URhAl does not belong to the 3D Ising model (beta = 0.325, gamma = 1.241, and delta = 4.82) with short-range exchange interactions between magnetic moments. The determined exponents can be explained with the results of the renormalization group approach for a two-dimensional (2D) Ising system coupled with long-range interactions decaying as J(r)~r^-(d+sigma) with sigma = 1.44. We suggest that the strong hybridization between the uranium 5f and rhodium 4d electrons in the U-Rh_I layer in the hexagonal crystal structure is a source of the low dimensional magnetic property. The present result is contrary to current understandings of the physical properties in a series of isostructural UTX uranium ferromagnets (T: transition metals, X: p-block elements) based on the 3D Ising model.


title: A Generic Top-N Recommendation Framework For Trading-off Accuracy, Novelty, and Coverage

abstract: Standard collaborative filtering approaches for top-N recommendation are biased toward popular items. As a result, they recommend items that users are likely aware of and under-represent long-tail items. This is inadequate, both for consumers who prefer novel items and because concentrating on popular items poorly covers the item space, whereas high item space coverage increases providers' revenue. We present an approach that relies on historical rating data to learn user long-tail novelty preferences. We integrate these preferences into a generic re-ranking framework that customizes balance between accuracy and coverage. We empirically validate that our proposedframework increases the novelty of recommendations. Furthermore, by promoting long-tail items to the right group of users, we significantly increase the system's coverage while scalably maintaining accuracy. Our framework also enables personalization of existing non-personalized algorithms, making them competitive with existing personalized algorithms in key performance metrics, including accuracy and coverage.


title: Non-Hamiltonian Kelvin wave generation on vortices in Bose-Einstein condensates

abstract: Ultra-cold quantum turbulence is expected to decay through a cascade of Kelvin waves. These helical excitations couple vorticity to the quantum fluid causing long wavelength phonon fluctuations in a Bose-Einstein condensate. This interaction is hypothesized to be the route to relaxation for turbulent tangles in quantum hydrodynamics. The local induction approximation is the lowest order approximation to the Biot-Savart velocity field induced by a vortex line and, because of its integrability, is thought to prohibit energy transfer by Kelvin waves. Using the Biot-Savart description, we derive a generalization to the local induction approximation which predicts that regions of large curvature can reconfigure themselves as Kelvin wave packets. While this generalization preserves the arclength metric, a quantity conserved under the Eulerian flow of vortex lines, it also introduces a non-Hamiltonian structure on the geometric properties of the vortex line. It is this non-Hamiltonian evolution of curvature and torsion which provides a resolution to the missing Kelvin wave motion. In this work, we derive corrections to the local induction approximation in powers of curvature and state them for utilization in vortex filament methods. Using the Hasimoto transformation, we arrive at a nonlinear integro-differential equation which reduces to a modified nonlinear Schr\"odinger type evolution of the curvature and torsion on the vortex line. We show that this modification seeks to disperse localized curvature profiles. At the same time, the non-Hamiltonian break in integrability bolsters the deforming curvature profile and simulations show that this dynamic results in Kelvin wave propagation along the dispersive vortex medium.

categories: cond-mat.quant-gas cond-mat.stat-mech math-ph math.MP


title: Prospects for Detecting light bosons at the FCC-ee and CEPC

abstract: We look at the prospects for detecting light bosons, $X$, at proposed Z factories assuming a production of $10^{12}$ Z bosons. Such a large yield is within the design goals of future FCC-ee and CEPC colliders. Specifically we look at the cases where $X$ is either a singlet scalar which mixes with the standard model Higgs or a vector boson with mass $1\lesssim M_X \lesssim 80$ GeV. We find that several channels are particularly promising for discovery prospects. In particular $Z\rightarrow f \bar{f} X$ and $Z \rightarrow V_Q X$ gives a promising signal above a very clean standard model background. We also discuss several channels that have too large a background to be useful.


title: Deep Learning for Causal Inference

abstract: In this paper, we propose deep learning techniques for econometrics, specifically for causal inference and for estimating individual as well as average treatment effects. The contribution of this paper is twofold: 1. For generalized neighbor matching to estimate individual and average treatment effects, we analyze the use of autoencoders for dimensionality reduction while maintaining the local neighborhood structure among the data points in the embedding space. This deep learning based technique is shown to perform better than simple k nearest neighbor matching for estimating treatment effects, especially when the data points have several features/covariates but reside in a low dimensional manifold in high dimensional space. We also observe better performance than manifold learning methods for neighbor matching. 2. Propensity score matching is one specific and popular way to perform matching in order to estimate average and individual treatment effects. We propose the use of deep neural networks (DNNs) for propensity score matching, and present a network called PropensityNet for this. This is a generalization of the logistic regression technique traditionally used to estimate propensity scores and we show empirically that DNNs perform better than logistic regression at propensity score matching. Code for both methods will be made available shortly on Github at: https://github.com/vikas84bf


title: Cavity-free quantum optomechanical cooling by atom-modulated radiation

abstract: We theoretically study the radiation-induced interaction between the mechanical motion of an oscillating mirror and a remotely trapped atomic cloud. When illuminated by continuous-wave radiation, the mirror motion will induce red and blue sideband radiation, which respectively increases and reduces motional excitation. We find that by suitably driving $\Lambda$-level atoms, the mirror correlation with a specific radiation sideband could be converted from the outgoing to the incoming radiation. Such process allows us to manipulate the heating and cooling effects. Particularly, we develop an optomechanical cooling strategy that can mutually cancel the heating effect of the outgoing and incoming radiations, thus the motional ground state is attainable by net cooling. Without the necessity of cavity installation or perfect alignment, our proposal complements other efforts in quantum cooling of macroscopic objects.


title: A Simple Nearly-Optimal Restart Scheme For Speeding-Up First Order Methods

abstract: We present a simple scheme for restarting first-order methods for convex optimization problems. Restarts are made based only on achieving specified decreases in objective values, the specified amounts being the same for all optimization problems. Unlike existing restart schemes, the scheme makes no attempt to learn parameter values characterizing the structure of an optimization problem, nor does it require any special information that would not be available in practice (unless the first-order method chosen to be employed in the scheme itself requires special information). As immediate corollaries to the main theorems, we show that when some well-known first-order methods are employed in the scheme, the resulting complexity bounds are nearly optimal for particular -- yet quite general -- classes of problems.


title: A Global Information Based Adaptive Threshold for Grouping Large Scale Global Optimization Problems

abstract: By taking the idea of divide-and-conquer, cooperative coevolution (CC) provides a powerful architecture for large scale global optimization (LSGO) problems, but its efficiency relies highly on the decomposition strategy. It has been shown that differential grouping (DG) performs well on decomposing LSGO problems by effectively detecting the interaction among decision variables. However, its decomposition accuracy depends highly on the threshold. To improve the decomposition accuracy of DG, a global information based adaptive threshold setting algorithm (GIAT) is proposed in this paper. On the one hand, by reducing the sensitivity of the indicator in DG to the roundoff error and the magnitude of contribution weight of subcomponent, we proposed a new indicator for two variables which is much more sensitive to their interaction. On the other hand, instead of setting the threshold only based on one pair of variables, the threshold is generated from the interaction information for all pair of variables. By conducting the experiments on two sets of LSGO benchmark functions, the correctness and robustness of this new indicator and GIAT were verified.


title: Drift and its mediation in terrestrial orbits

abstract: The slow deformation of terrestrial orbits in the medium range, subject to lunisolar resonances, is well approximated by a family of Hamiltonian flow with $2.5$ degree-of-freedom. The action variables of the system may experience chaotic variations and large drift that we may quantify. Using variational chaos indicators, we compute high-resolution portraits of the action space. Such refined meshes allow to reveal the existence of tori and structures filling chaotic regions. Our elaborate computations allow us to isolate precise initial conditions near specific zones of interest and study their asymptotic behaviour in time. Borrowing classical techniques of phase- space visualisation, we highlight how the drift is mediated by the complement of the numerically detected KAM tori.


title: Structural and Thermal Stability of Graphyne and Graphdiyne Nanoscroll Structures

abstract: Graphynes and graphdiynes are generic names for families of two-dimensional carbon allotropes, where acetylenic groups connect benzenoid-like hexagonal rings, with the co-existence of sp and sp2 hybridized carbon atoms. The main differences between graphynes and graphdiynes are the number of acetylenic groups (one and two for graphynes and graphdiynes, respectively). Similarly to graphene nanoscrolls, graphyne and graphdiynes nanoscrolls are nanosized membranes rolled up into papyrus-like structures. In this work we investigated through fully atomistic reactive molecular dynamics simulations the structural and thermal (up to 1000K) stability of alpha,beta,gamma-graphyne and alpha,beta,gamma-graphdiyne scrolls. Our results show that stable nanoscrolls can be formed for all the structures investigated here, although they are less stable than corresponding graphene scrolls. This can be explained as a consequence of the higher graphyne/graphdiyne structural porosity in relation to graphene, which results in decreased {\pi}-{\pi} stacking interactions.


title: Blue straggler stars beyond the Milky Way: a non-segregated population in the Large Magellanic Cloud cluster NGC 2213

abstract: Using the high-resolution observations obtained by the Hubble Space Telescope, we analyzed the blue straggler stars (BSSs) in the Large Magellanic Cloud cluster NGC 2213. We found that the radial distribution of BSSs is consistent with that of the normal giant stars in NGC 2213, showing no evidence of mass segregation. However, an analytic calculation carried out for these BSSs shows that they are already dynamically old, because the estimated half-mass relaxation time for these BSSs is significantly shorter than the isochronal age of the cluster. We also performed direct N-body simulations for a NGC 2213-like cluster to understand the dynamical processes that lead to this none-segregated radial distribution of BSSs. Our numerical simulation shows that the presence of black hole subsystems inside the cluster centre can significantly affect the dynamical evolution of BSSs. The combined effects of the delayed segregation, binary disruption and exchange interactions of BSS progenitor binaries may result in this none-segregated radial distribution of BSSs in NGC 2213.


title: Autoencoding topology

abstract: The problem of learning a manifold structure on a dataset is framed in terms of a generative model, to which we use ideas behind autoencoders (namely adversarial/Wasserstein autoencoders) to fit deep neural networks. From a machine learning perspective, the resulting structure, an atlas of a manifold, may be viewed as a combination of dimensionality reduction and "fuzzy" clustering.


title: GOODS-ALMA: 1.1 mm galaxy survey - I. Source catalogue and optically dark galaxies

abstract: We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes