title: Velocity Spectra and Model Spectrum in Non-Premixed Jet Flames

abstract: In this contribution, velocity and its dissipation spectra in turbulent non-premixed jet flames are investigated by using two Direct Numerical Simulation (DNS) databases of temporally evolving jets with different bulk jet Reynolds numbers. In the DNS differential diffusion and flame dynamics effects have been accounted for and the flames experience high levels of extinction followed by re-ignition. It turns out that the spectra extracted from different statistically homogeneous planes across the jets and selected times corresponding to various occurring flame regimes (extinction and re-ignition phases) collapse reasonably when normalized by Favre averaged turbulent quantities. Especially (1) in the inertial range, the $\kappa^{-5/3}$ power law is observed with the constant of proportionality of $C_{K}=2.3$. (2) In the dissipation range $\exp(\beta\kappa)$ scaling exists with $\beta=7.2$ (instead of $\beta=5.2$). (3) The location of the peak of the normalized 1D dissipation spectra is in a lower normalized waver number ($\approx 0.08$) compared to the peak of non-reactive model spectrum ($\approx 0.11$). (4) Finally, an adapted 3D dissipation model spectrum proposed which peaks at a normalized waver number $\approx 0.19$ (instead of $\approx 0.26$).


title: Fast Automatic Detection of Geological Boundaries from Multivariate Log Data Using Recurrence

abstract: Manual interpretation of data collected from drill holes for mineral or oil and gas exploration is time-consuming and subjective. Identification of geological boundaries and distinctive rock physical property domains is the first step of interpretation. We introduce a multivariate technique, that can identify geological boundaries from petrophysical or geochemical data. The method is based on time-series techniques that have been adapted to be applicable for detecting transitions in geological spatial data. This method allows for the use of multiple variables in detecting different lithological layers. Additionally, it reconstructs the phase space of a single drill-hole or well to be applicable for further investigations across other holes or wells. The computationally cheap method shows efficiency and accuracy in detecting boundaries between lithological layers, which we demonstrate using examples from mineral exploration boreholes and an offshore gas exploration well.


title: Location Order Recovery in Trails with Low Temporal Resolution

abstract: Researchers who study object movement problems related to topics like traffic flow analysis, patient monitoring, and software operation, need to know the correct order in which objects move. Here, we use the term trail to refer to a series of movements by an object. This paper introduces a new missing data problem that occurs when analyzing trails where there is inadequate temporal resolution on the events. The temporal resolution is inadequate when an object, which can only be in one place at one time, appears in the data to be in two or more locations at once. We refer to this lack of resolution as a broken point. Broken points prevent us from knowing the correct order of movement. We propose a three-phase framework for recovering the location order. Based on the Markov transition network, we are able to find the route with the highest probability. Our results show that this framework can efficiently find the correct location order in trails with low temporal resolution. We also demonstrate that by correcting the location order, the criticality of locations can change significantly.


title: Scanning tunneling microscopy and Raman evidences of silicene nanosheets intercalated into graphite surface at room temperature

abstract: Highly oriented pyrolitic graphite (HOPG) is an inert substrate with a structural honeycomb lattice, well suited for the growth of two-dimensional (2D) silicene layer. It was reported that when Si atoms are deposited on HOPG surface at room temperature, they arrange in two configurations: silicene nanosheets and three dimensional clusters. In this work we demonstrate, by using scanning tunneling microscopy (STM) and Raman spectroscopy, that a third configuration stabilizes in the form of Si 2D nanosheets intercalated below the first top layer of carbon atoms. The Raman spectra reveal a structure located at 538 cm$^{-1}$ which we ascribe to the presence of sp$^2$ Si hybridization. Moreover, the silicon deposition induces several modifications in the graphite D and G Raman modes, which we interpret as an experimental evidence of the intercalation of the silicene nanosheets. The Si atom intercalation at room temperature takes place at the HOPG step edges and it detaches only the outermost graphene layer inducing a strong tensile strain mainly concentrated on the edges of the silicene nanosheets. Theoretical calculations of the structure and energetic viability of the silicene nanosheets, of the strain distribution on the outermost graphene layer and its influence on the Raman resonances support the STM and Raman observations.


title: Evidence for spots on hot stars suggests major revision of stellar physics

abstract: It has long been thought that starspots are not present in the A and B stars because magnetic fields cannot be generated in stars with radiative envelopes. Space observations show that a considerable fraction of these stars vary in light with periods consistent with the expected rotation periods. Here we show that the photometric periods are the same as the rotation periods and that starspots are the likely cause for the light variations. This discovery has wide-ranging implications and suggests that a major revision of the physics of hot stellar envelopes may be required.


title: Hall conductivity of strained $Z_2$ crystals

abstract: We establish topological nature of Hall conductivity of graphene and other $Z_2$ crystals in 2D and 3D in the presence of inhomogeneous perturbations. To this end the lattice Weyl-Wigner formalism is employed. The non-uniform mechanical stress is considered, along with spatially varying magnetic field. The relation of the obtained topological invariant to level counting is clarified.


title: Secondary Inputs for Measuring User Engagement in Immersive VR Education Environments

abstract: This paper presents an experiment to assess the feasibility of using secondary input data as a method of determining user engagement in immersive virtual reality (VR). The work investigates whether secondary data (biosignals) acquired from users are useful as a method of detecting levels of concentration, stress, relaxation etc. in immersive environments, and if they could be used to create an affective feedback loop in immersive VR environments, including educational contexts. A VR Experience was developed in the Unity game engine, with three different levels, each designed to expose the user in one of three different states (relaxation, concentration, stress). While in the VR Experience users physiological responses were measured using ECG and EEG sensors. After the experience users completed questionnaires to establish their perceived state during the levels, and to established the usability of the system. Next a comparison between the reported levels of emotion and the measured signals is presented, which show a strong correspondence between the two measures indicating that biosignals are a useful indicator of emotional state while in VR. Finally we make some recommendations on the practicalities of using biosensors, and design considerations for their incorporation in to a VR system, with particular focus on their integration in to task-based training and educational virtual environments.


title: On the security and privacy of Interac e-Transfers

abstract: Nowadays, the Interac e-Transfer is one of the most important remote payment methods for Canadian consumers. To the best of our knowledge, this paper is the very first to examine the privacy and security of Interac e-Transfers. Experimental results show that the notifications sent to customers via email and SMS contain sensitive private information that can potentially be observed by third parties. Anyone with illegitimate intent can use this information to carry out attacks, including the fraudulent redirection of Standard e-Transfers. Such an attack is shown to be possible at least in an experimental setup but likely also in reality. Recent news articles support this finding. Improvements to overcome these interconnected privacy and security problems are proposed and discussed.


title: Probabilistic Robust Small-Signal Stability Framework using Gaussian Process Learning

abstract: While most power system small-signal stability assessments rely on the reduced Jacobian, which depends non-linearly on the states, uncertain operating points introduce nontrivial hurdles in certifying the system's stability. In this paper, a novel probabilistic robust small-signal stability (PRS) framework is developed for a power system based on Gaussian process (GP) learning. The proposed PRS assessment provides a robust stability certificate for a state subspace, such as that specified by the error bounds of the state estimation, with a given probability. With such a PRS certificate, all inner points of the concerned subspace will be stable with at least the corresponding confidence level. To this end, the behavior of the critical eigenvalue of the reduced Jacobian with state points in a state subspace is learned using GP. The proposed PRS certificate along with the Subspace-based Search and Confidence-based Search mechanisms constitute a holistic framework catering to all scenarios. The proposed framework is a powerful approach to assess the stability under uncertainty because it does not require input uncertainty distributions and other state-specific input-to-output approximations. Further, the critical eigenvalue behavior in a state subspace is analyzed using an upper bound of the eigenvalue variations and their inferences are discussed in detail. The results on three-machine nine-bus WSCC system show that the proposed certificate can find the robust stable state subspace with a given probability.


title: Graph Analysis and Graph Pooling in the Spatial Domain

abstract: The spatial convolution layer which is widely used in the Graph Neural Networks (GNNs) aggregates the feature vector of each node with the feature vectors of its neighboring nodes. The GNN is not aware of the locations of the nodes in the global structure of the graph and when the local structures corresponding to different nodes are similar to each other, the convolution layer maps all those nodes to similar or same feature vectors in the continuous feature space. Therefore, the GNN cannot distinguish two graphs if their difference is not in their local structures. In addition, when the nodes are not labeled/attributed the convolution layers can fail to distinguish even different local structures. In this paper, we propose an effective solution to address this problem of the GNNs. The proposed approach leverages a spatial representation of the graph which makes the neural network aware of the differences between the nodes and also their locations in the graph. The spatial representation which is equivalent to a point-cloud representation of the graph is obtained by a graph embedding method. Using the proposed approach, the local feature extractor of the GNN distinguishes similar local structures in different locations of the graph and the GNN infers the topological structure of the graph from the spatial distribution of the locally extracted feature vectors. Moreover, the spatial representation is utilized to simplify the graph down-sampling problem. A new graph pooling method is proposed and it is shown that the proposed pooling method achieves competitive or better results in comparison with the state-of-the-art methods.


title: DPSOM: Deep Probabilistic Clustering with Self-Organizing Maps

abstract: Generating interpretable visualizations from complex data is a common problem in many applications. Two key ingredients for tackling this issue are clustering and representation learning. However, current methods do not yet successfully combine the strengths of these two approaches. Existing representation learning models which rely on latent topological structure such as self-organising maps, exhibit markedly lower clustering performance compared to recent deep clustering methods. To close this performance gap, we (a) present a novel way to fit self-organizing maps with probabilistic cluster assignments (PSOM), (b) propose a new deep architecture for probabilistic clustering (DPSOM) using a VAE, and (c) extend our architecture for time-series clustering (T-DPSOM), which also allows forecasting in the latent space using LSTMs. We show that DPSOM achieves superior clustering performance compared to current deep clustering methods on MNIST/Fashion-MNIST, while maintaining the favourable visualization properties of SOMs. On medical time series, we show that T-DPSOM outperforms baseline methods in time series clustering and time series forecasting, while providing interpretable visualizations of patient state trajectories and uncertainty estimation.


title: Evidence for a pressure-induced phase transition of few-layer graphene to 2D diamond

abstract: We unveil the diamondization mechanism of few-layer graphene compressed in the presence of water, providing robust evidence for the pressure-induced formation of 2D diamond. High-pressure Raman spectroscopy provides evidence of a phase transition occurring in the range of 4-7 GPa for 5-layer graphene and graphite. The pressure-induced phase is partially transparent and indents the silicon substrate. Our combined theoretical and experimental results indicate a gradual top-bottom diamondization mechanism, consistent with the formation of diamondene, a 2D ferromagnetic semiconductor. High-pressure x-ray diffraction on graphene indicates the formation of hexagonal diamond, consistent with the bulk limit of eclipsed-conformed diamondene.


title: Efficient training of energy-based models via spin-glass control

abstract: We introduce a new family of energy-based probabilistic graphical models for efficient unsupervised learning. Its definition is motivated by the control of the spin-glass properties of the Ising model described by the weights of Boltzmann machines. We use it to learn the Bars and Stripes dataset of various sizes and the MNIST dataset, and show how they quickly achieve the performance offered by standard methods for unsupervised learning. Our results indicate that the standard initialization of Boltzmann machines with random weights equivalent to spin-glass models is an unnecessary bottleneck in the process of training. Furthermore, this new family allows for very easy access to low-energy configurations, which points to new, efficient training algorithms. The simplest variant of such algorithms approximates the negative phase of the log-likelihood gradient with no Markov chain Monte Carlo sampling costs at all, and with an accuracy sufficient to achieve good learning and generalization.


title: Engineering generalized Gibbs ensembles with trapped ions

abstract: The concept of generalized Gibbs ensembles (GGEs) has been introduced to describe steady states of integrable models. Recent advances show that GGEs can also be stabilized in nearly integrable quantum systems when driven by external fields and open. Here, we present a weakly dissipative dynamics that drives towards a steady-state GGE and is realistic to implement in systems of trapped ions. We outline the engineering of the desired dissipation by a combination of couplings which can be realized with ion-trap setups and discuss the experimental observables needed to detect a deviation from a thermal state. We present a novel mixed-species motional mode engineering technique in an array of micro-traps and demonstrate the possibility to use sympathetic cooling to construct many-body dissipators. Our work provides a blueprint for experimental observation of GGEs in open systems and opens a new avenue for quantum simulation of driven-dissipative quantum many-body problems.


title: Int-Deep: A Deep Learning Initialized Iterative Method for Nonlinear Problems

abstract: This paper focuses on proposing a deep learning initialized iterative method (Int-Deep) for low-dimensional nonlinear partial differential equations (PDEs). The corresponding framework consists of two phases. In the first phase, an expectation minimization problem formulated from a given nonlinear PDE is approximately resolved with mesh-free deep neural networks to parametrize the solution space. In the second phase, a solution ansatz of the finite element method to solve the given PDE is obtained from the approximate solution in the first phase, and the ansatz can serve as a good initial guess such that Newton's method for solving the nonlinear PDE is able to converge to the ground truth solution with high-accuracy quickly. Systematic theoretical analysis is provided to justify the Int-Deep framework for several classes of problems. Numerical results show that the Int-Deep outperforms existing purely deep learning-based methods or traditional iterative methods (e.g., Newton's method and the Picard iteration method).


title: Pattern engineering of living bacterial colonies using meniscus-driven fluidic channels

abstract: Engineering spatially organized biofilms for creating adaptive and sustainable biomaterials is a forthcoming mission of synthetic biology. Existing technologies of patterning biofilm materials suffer limitations associated with the high technical barrier and the requirements of special equipment. Here we present controlled meniscus-driven fluidics, MeniFluidics; an easily implementable technique for patterning living bacterial populations. We demonstrate multiscale patterning of living-colony and biofilm formation with submillimetre resolution. Relying on fast bacterial spreading in liquid channels, MeniFluidics allows controlled anisotropic bacterial colonies expansion both in space and time. The technique has also been applied for studying collective phenomena in confined bacterial swarming and organizing different fluorescently labelled Bacillus subtilis strains into a converged pattern. We believe that the robustness and low technical barrier of MeniFluidics offer a tool for developing living functional materials, bioengineering and bio-art, and adding to fundamental research of microbial interactions.


title: A New Dynamic Model of a Two-Wheeled Two-Flexible-Beam Inverted Pendulum Robot

abstract: Two-wheeled inverted pendulum robots are designed for self-balancing and they have remarkable advantages. In this paper, a new configuration and consequently dynamic model of one specific robot is presented and its dynamic behavior is analyzed. In this model, two cantilever beams are on the two-wheeled base and they are excited by voltages to the attached piezoelectric actuators. The mathematical model of this system is obtained using the extended Hamilton's Principle. The results show that the governing equations of motion are highly nonlinear and contain several coupled terms. These complex equations are solved numerically and the natural frequencies of the system are extracted. The system is then simulated in both lateral and horizontal plan movements. This paper proposes a new model of two-wheeled two-flexible-beam inverted pendulum Robot and investigates its complex dynamic; however, the derived equations will be validated experimentally and a suitable control strategy will be applied to the system to make it fully automated and more applicable in the future works.


title: Authentication Modeling with Five Generic Processes

abstract: Conceptual modeling is an essential tool in many fields of study, including security specification in information technology systems. As a model, it restricts access to resources and identifies possible threats to the system. We claim that current modeling languages (e.g., Unified Modeling Language, Business Process Model and Notation) lack the notion of genericity, which refers to a limited set of elementary processes. This paper proposes five generic processes for modeling the structural behavior of a system: creating, releasing, transferring, receiving, and processing. The paper demonstrates these processes within the context of public key infrastructure, biometric, and multifactor authentication. The results indicate that the proposed generic processes are sufficient to represent these authentication schemes.


title: Amplification of superkicks in black-hole binaries through orbital eccentricity

abstract: We present new numerical-relativity simulations of eccentric merging black holes with initially antiparallel spins lying in the orbital plane (the so-called \emph{superkick} configuration). Binary eccentricity boosts the recoil of the merger remnant by up to $25\%$. The increase in the energy flux is much more modest, and therefore this kick enhancement is mainly due to asymmetry in the binary dynamics. Our findings might have important consequences for the retention of stellar-mass black holes in star clusters and supermassive black holes in galactic hosts.


title: Non-Linear Supersymmetry and $T\bar T$-like Flows

abstract: The $T\bar{T}$ deformation of a supersymmetric two-dimensional theory preserves the original supersymmetry. Moreover, in several interesting cases the deformed theory possesses additional non-linearly realized supersymmetries. We show this for certain $\mathcal N =(2,2)$ models in two dimensions, where we observe an intriguing similarity with known $\mathcal N=1 $ models in four dimensions. This suggests that higher-dimensional models with non-linearly realized supersymmetries might also be obtained from $T\bar{T}$-like flow equations. We show that in four dimensions this is indeed the case for $\mathcal N=1 $ Born-Infeld theory, as well as for the Goldstino action for spontaneously broken $\mathcal N=1 $ supersymmetry.


title: Three-qubit refrigerator with two-body interactions

abstract: We propose a three-qubit setup for the implementation of a variety of quantum thermal machines where all heat fluxes and work production can be controlled. An important configuration that can be designed is that of an absorption refrigerator, extracting heat from the coldest reservoir without the need of external work supply. Remarkably, we achieve this regime by using only two-body interactions instead of the widely employed three-body interactions. This configuration could be more easily realised in current experimental setups. We model the open-system dynamics with both a global and a local master equation thermodynamic-consistent approach. Finally, we show how this model can be employed as a heat valve, in which by varying the local field of one of the two qubits allows one to control and amplify the heat current between the other qubits.


title: SensorDrop: A Reinforcement Learning Framework for Communication Overhead Reduction on the Edge

abstract: In IoT solutions, it is usually desirable to collect data from a large number of distributed IoT sensors at a central node in the cloud for further processing. One of the main design challenges of such solutions is the high communication overhead between the sensors and the central node (especially for multimedia data). In this paper, we aim to reduce the communication overhead and propose a method that is able to determine which sensors should send their data to the central node and which to drop data. The idea is that some sensors may have data which are correlated with others and some may have data that are not essential for the operation to be performed at the central node. As such decisions are application dependent and may change over time, they should be learned during the operation of the system, for that we propose a method based on Advantage Actor-Critic (A2C) reinforcement learning which gradually learns which sensor's data is cost-effective to be sent to the central node. The proposed approach has been evaluated on a multi-view multi-camera dataset, and we observe a significant reduction in communication overhead with marginal degradation in object classification accuracy.


title: Loops in surfaces and star-fillings

abstract: We discuss a new approach to computing the standard algebraic operations on homotopy classes of loops in surfaces: the homological intersection number, Goldman's Lie bracket, and the author's Lie cobracket. Our approach uses fillings of the surfaces by certain graphs.


title: Bootstrapping Conditional GANs for Video Game Level Generation

abstract: Generative Adversarial Networks (GANs) have shown im-pressive results for image generation. However, GANs facechallenges in generating contents with certain types of con-straints, such as game levels. Specifically, it is difficult togenerate levels that have aesthetic appeal and are playable atthe same time. Additionally, because training data usually islimited, it is challenging to generate unique levels with cur-rent GANs. In this paper, we propose a new GAN architec-ture namedConditional Embedding Self-Attention Genera-tive Adversarial Network(CESAGAN) and a new bootstrap-ping training procedure. The CESAGAN is a modification ofthe self-attention GAN that incorporates an embedding fea-ture vector input to condition the training of the discriminatorand generator. This allows the network to model non-localdependency between game objects, and to count objects. Ad-ditionally, to reduce the number of levels necessary to trainthe GAN, we propose a bootstrapping mechanism in whichplayable generated levels are added to the training set. Theresults demonstrate that the new approach does not only gen-erate a larger number of levels that are playable but also gen-erates fewer duplicate levels compared to a standard GAN.


title: On shock waves from the inhomogeneous Boltzmann equation

abstract: We revisit the problem on the inner structure of shock waves in simple gases modelized by the Boltzmann kinetic equation. In \cite{pomeau1987shock}, a self-similarity approach was proposed for infinite total cross section resulting from a power law interaction, but this self-similar form does not have finite energy. Motivated by the work of Pomeau, Bobylev and Cercignani started the rigorous study of the solutions of the spatial homogeneous Boltzmann equation, focusing on those which do not have finite energy \cite{bobylev2002self,bobylev2003eternal}. In the present work, we provide a correction to the self-similar form, so that the solutions are more physically sound in the sense that the energy is no longer infinite and that the perturbation brought by the shock does not grow at large distances of it on the cold side in the soft potential case.


title: Photoemission Spectrum of Ca2RuO4: Spin Polaron Physics in an S=1 Antiferromagnet with Anisotropies

abstract: We derive an S=1 spin polaron model which describes the motion of a single hole introduced into the S=1 spin antiferromagnetic ground state of Ca2RuO4. We solve the model using the self-consistent Born approximation and show that its hole spectral function qualitatively agrees with the experimentally observed high-binding energy part of the Ca2RuO4 photoemission spectrum. We explain the observed peculiarities of the photoemission spectrum by linking them to two anisotropies present in the employed model---the spin anisotropy and the hopping anisotropy. We verify that these anisotropies, and not the possible differences between the ruthenate (S=1) and the cuprate (S=1/2) spin polaron models, are responsible for the strong qualitative differences between the photoemission spectrum of Ca2RuO4 and of the undoped cuprates.


title: Resurgence and holonomy of the $\phi^{2k}$ model in zero dimension

abstract: We describe the resurgence properties of some partition functions corresponding to field theories in dimension 0. We show that these functions satisfy linear differential equations with polynomial coefficients and then use elementary stability results or holonomic functions to prove resurgence properties, enhancing previously known results on growth estimates for the formal series involved, which had been obtained through a delicate combinatorics.


title: Exploring the atmospheric dynamics of the extreme ultra-hot Jupiter KELT-9b using TESS photometry

abstract: We carry out a phase-curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite (TESS). The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of the secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hot spot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase-curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase-curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $\mu$m transmission spectrum that is featureless at low resolution.


title: Cosmological discordances III: more on measure properties, Large-Scale-Structure constraints, the Hubble constant and Planck

abstract: Consistency between cosmological data sets is essential for ongoing and future cosmological analyses. We first investigate the questions of stability and applicability of some moment-based inconsistency measures to multiple data sets. We show that the recently introduced index of inconsistency (IOI) is numerically stable while it can be applied to multiple data sets. We use an illustrative construction of constraints as well as an example with real data sets (i.e. WMAP versus Planck) to show some limitations of the application of the Karhunen-Loeve decomposition to discordance measures. Second, we perform various consistency analyzes using IOI between multiple current data sets while \textit{working with the entire common parameter spaces}. We find current Large-Scale-Structure (LSS) data sets (Planck CMB lensing, DES lensing-clustering and SDSS RSD) all to be consistent with one another. This is found to be not the case for Planck temperature (TT) versus polarization (TE,EE) data, where moderate inconsistencies are present. Noteworthy, we find a strong inconsistency between joint LSS probes and Planck with IOI=5.27, and a moderate tension between DES and Planck with IOI=3.14. Next, using the IOI metric, we compare the Hubble constant from five independent probes. We confirm previous strong tensions between local measurement (SH0ES) and Planck as well as between H0LiCOW and Planck, but also find new strong tensions between SH0ES measurement and the joint LSS probes with IOI=6.73 (i.e. 3.7-$\sigma$ in 1D) as well as between joint LSS and combined probes SH0ES+H0LiCOW with IOI=8.59 (i.e. 4.1-$\sigma$ in 1D). Whether due to systematic effects in the data sets or problems with the underlying model, sources of these old and new tensions need to be identified and dealt with.


title: Helical structures in layered magnetic superconductors due to indirect exchange interactions mediated by interlayer tunneling

abstract: Motivated by the recent discovery of helical magnetic structure in RbEuFe$_{4}$As$_{4}$, we investigate interlayer ordering of magnetic moments in materials composed of spatially-separated superconducting and ferromagnetically-aligned layers. We consider the interplay between the normal and superconducting indirect exchange interaction mediated by tunneling between the conducting layers. We elaborate a recipe to evaluate the normal interlayer interaction via two-dimensional density of states of an isolated layer and demonstrate that for bands with small fillings, such interaction is typically ferromagnetic and short-range. The nearest-layer interaction is proportional to the ratio of the interlayer hopping and in-plane band width squared. On the other hand, the superconducting contribution always gives antiferromagnetic interaction and may extend over several layers when the interlayer hopping energy exceeds the superconducting gap. The frustration caused by the interplay between the normal and superconducting parts may lead to spiral ground-state magnetic configuration. The four-fold in-plane anisotropy may lock the rotation angle between the moments in the neighboring layers to $90^\circ$, as it was observed in RbEuFe$_{4}$As$_{4}$


title: Strong cosmic censorship in de Sitter space: As strong as ever

abstract: The strong cosmic censorship has recently been put into question for the charged black holes in de Sitter space. We have performed the full non-linear evolution of the massless charged scalar field minimally coupled to the Einstein-Maxwell system in de Sitter space, and found that the non-linear effect can restore the strong cosmic censorship, making it as strong as ever.


title: On counting centralizer subgroups of symmetric groups

abstract: Let $S_{2m}$ be the symmetric group, $h=(1\ 2)(3\ 4)\cdots(2m-1\ 2m)$ and $H=C(h)$. We consider the structure of $gHg^{-1}\cap H$ for any $g\in S_{2m}$. We prove the permutations $g$ which makes $gHg^{-1}\cap H$ have size of polynomial in $m$ have density zero.


title: Partial differential equation regularization for supervised machine learning

abstract: This article is an overview of supervised machine learning problems for regression and classification. Topics include: kernel methods, training by stochastic gradient descent, deep learning architecture, losses for classification, statistical learning theory, and dimension independent generalization bounds. Implicit regularization in deep learning examples are presented, including data augmentation, adversarial training, and additive noise. These methods are reframed as explicit gradient regularization.


title: Velocity Decay Estimates for Boltzmann equation with hard potentials

abstract: We establish pointwise polynomial decay estimates in velocity space for the spatially inhomogeneous Boltzmann equation without cutoff, in the case of hard potentials ($\gamma +2s > 2$), under the assumption that the mass, energy, and entropy densities are bounded above, and the mass density is bounded below. These estimates are self-generating, i.e. they do not require corresponding decay assumptions on the initial data. Our results extend the recent work of Imbert-Mouhot-Silvestre (arXiv:1804.06135), which addressed the case of moderately soft potentials ($\gamma + 2s \in [0,2]$).


title: Invariant embeddings of unimodular random planar graphs

abstract: Consider an ergodic unimodular random one-ended planar graph $\G$ of finite expected degree. We prove that it has an isometry-invariant locally finite embedding in the Euclidean plane if and only if it is invariantly amenable. By "locally finite" we mean that any bounded open set intersects finitely many embedded edges. In particular, there exist invariant embeddings in the Euclidean plane for the Uniform Infinite Planar Triangulation and for the critical Augmented Galton-Watson Tree conditioned to survive. Roughly speaking, a unimodular embedding of $\G$ is one that is jointly unimodular with $\G$ when viewed as a decoration. We show that $\G$ has a unimodular embedding in the hyperbolic plane if it is invariantly nonamenable, and it has a unimodular embedding in the Euclidean plane if and only if it is invariantly amenable. Similar claims hold for representations by tilings instead of embeddings.


title: Fair Division of Goods in the Shadow of Market Values

abstract: Inheritances, divorces or liquidations of companies require common assets to be divided among the entitled parties. Legal methods usually consider the market value of goods, while fair division theory takes into account the parties' preferences expressed as utilities. I combine the two practices to define a procedure that optimally allocates divisible goods with market values to people with easily elicited preferences. Imposing exact equality in the bundles' monetary values may produce unacceptable solutions. I drop the tight requirement and suggest a procedure in which the differences in the monetary values are explained in terms of satisfaction per monetary share as perceived by the agents. A robustness study shows the consequences of misspecification in the model parameters.


title: A Stochastic-Robust Approach to Hierarchical Generation-Transmission Expansion Planning

abstract: In this paper, the differences between an integrated and hierarchical generation and transmission expansion planning approaches are first described. Then, the hierarchical approach is described in detail. In general terms, in this scheme the investment decisions are made in steps, instead of an overall optimization scheme. This paper proposes a stochastic hierarchical generation-transmission expansion planning methodology based on a three-step procedure, as follows: 1. In this step, an integrated expansion planning problem of generation and interregional interconnections is solved; 2. Taking the optimal expansion plan of (generation and interconnections) into account, a production costing simulation with the detailed network representation is performed without monitoring circuit flow limits (except the ones in the interconnections, which are monitored); this simulation produces a set of optimal dispatch scenarios (vectors of bus loads and generation); 3. A Transmission Expansion Planning (TEP) model is then applied to determine the least-cost transmission expansion plan that is robust with respect to all operation scenarios of step 2, using an enhanced Benders decomposition scheme that: (a) incorporates a subset of the operation scenarios in the investment module; and (b) presents a "warm-up" step with a "greedy" algorithm that produces a (good) feasible solution and an initial set of feasibility cuts. The application of the hierarchical planning scheme is illustrated with a realistic multi-country generation and transmission planning case study of the Central America's electricity market.


title: Kilonovae

abstract: The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergo rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a ``kilonova', which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ~hours-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.


title: Inference of a mesoscopic population model from population spike trains

abstract: To understand how rich dynamics emerge in neural populations, we require models exhibiting a wide range of activity patterns while remaining interpretable in terms of connectivity and single-neuron dynamics. However, it has been challenging to fit such mechanistic spiking networks at the single neuron scale to empirical population data. To close this gap, we propose to fit such data at a meso scale, using a mechanistic but low-dimensional and hence statistically tractable model. The mesoscopic representation is obtained by approximating a population of neurons as multiple homogeneous `pools' of neurons, and modelling the dynamics of the aggregate population activity within each pool. We derive the likelihood of both single-neuron and connectivity parameters given this activity, which can then be used to either optimize parameters by gradient ascent on the log-likelihood, or to perform Bayesian inference using Markov Chain Monte Carlo (MCMC) sampling. We illustrate this approach using a model of generalized integrate-and-fire neurons for which mesoscopic dynamics have been previously derived, and show that both single-neuron and connectivity parameters can be recovered from simulated data. In particular, our inference method extracts posterior correlations between model parameters, which define parameter subsets able to reproduce the data. We compute the Bayesian posterior for combinations of parameters using MCMC sampling and investigate how the approximations inherent to a mesoscopic population model impact the accuracy of the inferred single-neuron parameters.


title: Beyond Linearization: On Quadratic and Higher-Order Approximation of Wide Neural Networks

abstract: Recent theoretical work has established connections between over-parametrized neural networks and linearized models governed by he Neural Tangent Kernels (NTKs). NTK theory leads to concrete convergence and generalization results, yet the empirical performance of neural networks are observed to exceed their linearized models, suggesting insufficiency of this theory. Towards closing this gap, we investigate the training of over-parametrized neural networks that are beyond the NTK regime yet still governed by the Taylor expansion of the network. We bring forward the idea of \emph{randomizing} the neural networks, which allows them to escape their NTK and couple with quadratic models. We show that the optimization landscape of randomized two-layer networks are nice and amenable to escaping-saddle algorithms. We prove concrete generalization and expressivity results on these randomized networks, which lead to sample complexity bounds (of learning certain simple functions) that match the NTK and can in addition be better by a dimension factor when mild distributional assumptions are present. We demonstrate that our randomization technique can be generalized systematically beyond the quadratic case, by using it to find networks that are coupled with higher-order terms in their Taylor series.


title: Optimal Performance of a Three-level Quantum Refrigerator

abstract: We study the optimal performance of a three-level quantum refrigerator using two different objective functions: cooling power and $\chi$-function. For both cases, we obtain general expressions for the coefficient of performance (COP) and derive its well-known lower and upper bounds for the limiting cases when the ratio of system-bath coupling constants at the hot and cold contacts approaches infinity and zero, respectively. We also show that the cooling power is optimizable only in the local region with respect to one control frequency, while $\chi$-function can be optimized globally with respect to two control frequencies. Additionally, we show that in the low-temperatures regime, our model of refrigerator can be mapped to Feynman's ratchet and pawl model, a classical mesoscopic heat engine. In the parameter regime where both cooling power and $\chi$-function can be optimized, we compare the cooling power of the quantum refrigerator at maximum $\chi$-function with the optimum cooling power.


title: Supersymmetry and Hodge theory on Sasakian and Vaisman manifolds

abstract: Sasakian manifolds are odd-dimensional counterpart to Kahler manifolds. They can be defined as contact manifolds equipped with an invariant Kahler structure on their symplectic cone. The quotient of this cone by the homothety action is a complex manifold called Vaisman. We study harmonic forms and Hodge decomposition on Vaisman and Sasakian manifolds. We construct a Lie superalgebra associated to a Sasakian manifold in the same way as the Kahler supersymmetry algebra is associated to a Kahler manifold. We use this construction to produce a self-contained, coordinate-free proof of the results by Tachibana, Kashiwada and Sato on the decomposition of harmonic forms and cohomology of Sasakian and Vaisman manifolds. In the last section, we compute the supersymmetry algebra of Sasakian manifolds explicitly.


title: The Atmospheric Circulation of Ultra-hot Jupiters

abstract: Recent observations of ultra-hot Jupiters with dayside temperatures in excess of $2500~\mathrm{K}$ have found evidence for new physical processes at play in their atmospheres. In this work, we investigate the effects of the dissociation of molecular hydrogen and recombination of atomic hydrogen on the atmospheric circulation of ultra-hot Jupiters. To do so, we incorporate these effects into a general circulation model (GCM) for hot Jupiter atmospheres, and run a large suite of models varying the incident stellar flux, rotation period, and strength of frictional drag. We find that including hydrogen dissociation and recombination reduces the fractional day-to-night temperature contrast of ultra-hot Jupiter atmospheres and causes the speed of the equatorial jet to decrease in simulations with fixed rotation. This is because the large energy input required for hydrogen dissociation cools the dayside of the planet, and the energy released due to hydrogen recombination warms the nightside. The resulting decrease in the day-to-night temperature contrast reduces the day-to-night pressure gradient that drives the circulation, resulting in weaker wind speeds. The results from our GCM experiments qualitatively agree with previous theory which found that the fractional day-night temperature contrast of ultra-hot Jupiters should decrease with increasing equilibrium temperature due to hydrogen dissociation and recombination. Lastly, we compute full-phase light curves from our suite of GCM experiments, finding that the reduced day-to-night temperature contrast in ultra-hot Jupiter atmospheres causes a smaller phase curve amplitude. The reduction in phase curve amplitude due to hydrogen dissociation and recombination could explain the relatively small phase curve amplitudes of observed ultra-hot Jupiters.


title: A Pseudo-Likelihood Approach to Linear Regression with Partially Shuffled Data

abstract: Recently, there has been significant interest in linear regression in the situation where predictors and responses are not observed in matching pairs corresponding to the same statistical unit as a consequence of separate data collection and uncertainty in data integration. Mismatched pairs can considerably impact the model fit and disrupt the estimation of regression parameters. In this paper, we present a method to adjust for such mismatches under ``partial shuffling" in which a sufficiently large fraction of (predictors, response)-pairs are observed in their correct correspondence. The proposed approach is based on a pseudo-likelihood in which each term takes the form of a two-component mixture density. Expectation-Maximization schemes are proposed for optimization, which (i) scale favorably in the number of samples, and (ii) achieve excellent statistical performance relative to an oracle that has access to the correct pairings as certified by simulations and case studies. In particular, the proposed approach can tolerate considerably larger fraction of mismatches than existing approaches, and enables estimation of the noise level as well as the fraction of mismatches. Inference for the resulting estimator (standard errors, confidence intervals) can be based on established theory for composite likelihood estimation. Along the way, we also propose a statistical test for the presence of mismatches and establish its consistency under suitable conditions.


title: Verification of Neural Network Behaviour: Formal Guarantees for Power System Applications

abstract: This paper presents for the first time, to our knowledge, a framework for verifying neural network behavior in power system applications. Up to this moment, neural networks have been applied in power systems as a black-box; this has presented a major barrier for their adoption in practice. Developing a rigorous framework based on mixed integer linear programming, our methods can determine the range of inputs that neural networks classify as safe or unsafe, and are able to systematically identify adversarial examples. Such methods have the potential to build the missing trust of power system operators on neural networks, and unlock a series of new applications in power systems. This paper presents the framework, methods to assess and improve neural network robustness in power systems, and addresses concerns related to scalability and accuracy. We demonstrate our methods on the IEEE 9-bus, 14-bus, and 162-bus systems, treating both N-1 security and small-signal stability.


title: Minimax Bounds for Distributed Logistic Regression

abstract: We consider a distributed logistic regression problem where labeled data pairs $(X_i,Y_i)\in \mathbb{R}^d\times\{-1,1\}$ for $i=1,\ldots,n$ are distributed across multiple machines in a network and must be communicated to a centralized estimator using at most $k$ bits per labeled pair. We assume that the data $X_i$ come independently from some distribution $P_X$, and that the distribution of $Y_i$ conditioned on $X_i$ follows a logistic model with some parameter $\theta\in\mathbb{R}^d$. By using a Fisher information argument, we give minimax lower bounds for estimating $\theta$ under different assumptions on the tail of the distribution $P_X$. We consider both $\ell^2$ and logistic losses, and show that for the logistic loss our sub-Gaussian lower bound is order-optimal and cannot be improved.


title: The isomorphic Kottman constant of a Banach space

abstract: We show that the Kottman constant $K(\cdot)$, together with its symmetric and finite variations, is continuous with respect to the Kadets metric, and they are log-convex, hence continuous, with respect to the interpolation parameter in a complex interpolation schema. Moreover, we show that $K(X)\cdot K(X^*)\geqslant 2$ for every infinite-dimensional Banach space $X$. We also consider the isomorphic Kottman constant (defined as the infimum of the Kottman constants taken over all renormings of the space) and solve the main problem left open in [CaGoPa17], namely that the isomorphic Kottman constant of a twisted-sum space is the maximum of the constants of the respective summands. Consequently, the Kalton--Peck space may be renormed to have Kottman's constant arbitrarily close to $\sqrt{2}$. For other classical parameters, such as the Whitley and the James constants, we prove the continuity with respect to the Kadets metric.


title: Large degrees in scale-free inhomogeneous random graphs

abstract: We consider a class of scale-free inhomogeneous random graphs, which includes some long-range percolation models. We study the maximum degree in such graphs in a growing observation window and show that its limiting distribution is Frechet. We achieve this by proving convergence of the underlying point process of the degrees to a certain Poisson process. Estimating the index of the power-law tail for the typical degree distribution is an important question in statistics. We prove consistency of the Hill estimator for the inverse of the tail exponent of the typical degree distribution.


title: Response of a quantum disordered spin system to a local periodic drive

abstract: We consider a one-dimensional spin chain system with quenched disorder and in the presence of a local periodic drive. We study the time evolution of the system in the Floquet basis and evaluate the fidelity susceptibility, which is a measure of how a given state changes under a small perturbation, of states to a weak periodic drive. We demonstrate that the statistical properties of the fidelity susceptibility over different disorder realizations can be used to identify two phases of the system: (1) the many-body localized phase, in which the susceptibility exhibits long tails while its average value decreases rapidly as disorder increases; and (2) the ergodic phase, in which the susceptibility distribution is narrow and its average value weakly depends on disorder. This distinction in the average value of the susceptibility between the two phases develops readily for systems with ten or more spins. Therefore, recently built networks of qubits subject to a local drive can simulate dynamics of a system in the many-body localization regime. We also show that the spin accumulation speed is correlated with the fidelity susceptibility and can also be used to distinguish the two phases.


title: On the Stable Resolution Limit of Total Variation Regularization for Spike Deconvolution

abstract: The stability of spike deconvolution, which aims at recovering point sources from their convolution with a point spread function (PSF), is known to be related to the separation between those sources. When the observations are noisy, it is critical to ensure support stability, where the deconvolution does not lead to spurious, or oppositely, missing estimates of the point sources. In this paper, we study the resolution limit of stably recovering the support of two closely located point sources using the Beurling-LASSO estimator, which is a convex optimization approach based on total variation regularization. We establish a sufficient separation criterion between the sources, depending only on the PSF, above which the Beurling-LASSO estimator is guaranteed to return a stable estimate of the point sources, with the same number of estimated elements as that of the ground truth. Our result highlights the impact of PSF on the resolution limit in the noisy setting, which was not evident in previous studies of the noiseless setting. Towards the end, we show that the same resolution limit applies to resolving two close-located sources in conjunction of other well-separated sources.


title: Strong planar subsystem symmetry-protected topological phases and their dual fracton orders

abstract: We classify subsystem symmetry-protected topological (SSPT) phases in $3+1$D protected by planar subsystem symmetries, which are dual to abelian fracton topological orders. We distinguish between weak SSPTs, which can be constructed by stacking $2+1$D SPTs, and strong SSPTs, which cannot. We identify signatures of strong phases, and show by explicit construction that such phases exist. A classification of strong phases is presented for an arbitrary finite abelian group. Finally, we show that fracton orders realizable via $p$-string condensation are dual to weak SSPTs, while strong SSPTs do not admit such a realization.

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


title: Uncomputability of Phase Diagrams

abstract: The phase diagram of a material is of central importance to describe the properties and behaviour of a condensed matter system. We prove that the general task of determining the quantum phase diagram of a many-body Hamiltonian is uncomputable, by explicitly constructing a one-parameter family of Hamiltonians for which this is the case. This work builds off recent results from Cubitt et al. and Bausch et al., proving undecidability of the spectral gap problem. However, in all previous constructions, the Hamiltonian was necessarily a discontinuous function of its parameters, making it difficult to derive rigorous implications for phase diagrams or related condensed matter questions. Our main technical contribution is to prove undecidability of the spectral gap for a continuous, single-parameter family of translationally invariant, nearest-neighbour spin-lattice Hamiltonians on a 2D square lattice: $H(\varphi)$ where $\varphi\in \mathbb R$. As well as implying uncomputablity of phase diagrams, our result also proves that undecidability can hold for a set of positive measure of a Hamiltonian's parameter space, whereas previous results only implied undecidability on a zero measure set.


title: Elastic multipole method for describing linear deformation of infinite 2D solid structures with circular holes and inclusions

abstract: Elastic materials with holes and inclusions are important in a large variety of contexts ranging from construction material to biological membranes. More recently, they have also been exploited in mechanical metamaterials, where the geometry of highly deformable structures is responsible for their unusual properties, such as negative Poisson's ratio, mechanical cloaking, and tunable phononic band gaps. Understanding how such structures deform in response to applied external loads is thus crucial for designing novel mechanical metamaterials. Here we present a method for predicting the linear response of infinite 2D solid structures with circular holes and inclusions by employing analogies with electrostatics. Just like an external electric field induces polarization (dipoles, quadrupoles and other multipoles) of conductive and dielectric objects, external stress induces elastic multipoles inside holes and inclusions. Stresses generated by these induced elastic multipoles then lead to interactions between holes and inclusions, which induce additional polarization and thus additional deformation of holes and inclusions. We present a method that expands the induced polarization in a series of elastic multipoles, which systematically takes into account the interactions of inclusions and holes with the external stress field and also between them. The results of our method show good agreement with both linear finite element simulations and experiments.


title: The Tutte polynomial and toric Nakajima quiver varieties

abstract: For a quiver $Q$, we take $\mathcal{M}$ an associated toric Nakajima quiver variety and $\Gamma$ the underlying graph. In this article, we give a direct relation between a specialisation of the Tutte polynomial of $\Gamma$, the Kac polynomial of $Q$ and the Poincar\'e polynomial of $\mathcal{M}$. We do this by giving a cell decomposition of $\mathcal{M}$ indexed by spanning trees of $\Gamma$ and `geometrising' the deletion and contraction operators on graphs. These relations have been previously established by Sturmfels-Hausel and (Crawley-Boovey)-Van den Bergh, however the methods here are more hands-on.


title: Improving Limited Angle CT Reconstruction with a Robust GAN Prior

abstract: Limited angle CT reconstruction is an under-determined linear inverse problem that requires appropriate regularization techniques to be solved. In this work we study how pre-trained generative adversarial networks (GANs) can be used to clean noisy, highly artifact laden reconstructions from conventional techniques, by effectively projecting onto the inferred image manifold. In particular, we use a robust version of the popularly used GAN prior for inverse problems, based on a recent technique called corruption mimicking, that significantly improves the reconstruction quality. The proposed approach operates in the image space directly, as a result of which it does not need to be trained or require access to the measurement model, is scanner agnostic, and can work over a wide range of sensing scenarios.


title: A Function Space View of Bounded Norm Infinite Width ReLU Nets: The Multivariate Case

abstract: A key element of understanding the efficacy of overparameterized neural networks is characterizing how they represent functions as the number of weights in the network approaches infinity. In this paper, we characterize the norm required to realize a function $f:\mathbb{R}^d\rightarrow\mathbb{R}$ as a single hidden-layer ReLU network with an unbounded number of units (infinite width), but where the Euclidean norm of the weights is bounded, including precisely characterizing which functions can be realized with finite norm. This was settled for univariate univariate functions in Savarese et al. (2019), where it was shown that the required norm is determined by the L1-norm of the second derivative of the function. We extend the characterization to multivariate functions (i.e., networks with d input units), relating the required norm to the L1-norm of the Radon transform of a (d+1)/2-power Laplacian of the function. This characterization allows us to show that all functions in Sobolev spaces $W^{s,1}(\mathbb{R})$, $s\geq d+1$, can be represented with bounded norm, to calculate the required norm for several specific functions, and to obtain a depth separation result. These results have important implications for understanding generalization performance and the distinction between neural networks and more traditional kernel learning.


title: Self-supervised learning for autonomous vehicles perception: A conciliation between analytical and learning methods

abstract: Nowadays, supervised deep learning techniques yield the best state-of-the-art prediction performances for a wide variety of computer vision tasks. However, such supervised techniques generally require a large amount of manually labeled training data. In the context of autonomous vehicles perception, this requirement is critical, as the distribution of sensor data can continuously change and include several unexpected variations. It turns out that a category of learning techniques, referred to as self-supervised learning (SSL), consists of replacing the manual labeling effort by an automatic labeling process. Thanks to their ability to learn on the application time and in varying environments, state-of-the-art SSL techniques provide a valid alternative to supervised learning for a variety of different tasks, including long-range traversable area segmentation, moving obstacle instance segmentation, long-term moving obstacle tracking, or depth map prediction. In this tutorial-style article, we present an overview and a general formalization of the concept of self-supervised learning (SSL) for autonomous vehicles perception. This formalization provides helpful guidelines for developing novel frameworks based on generic SSL principles. Moreover, it enables to point out significant challenges in the design of future SSL systems.


title: Fast simulation of detector effects in Rivet

abstract: We describe the design and implementation of detector-bias emulation in the Rivet MC event analysis system. Implemented using C++ efficiency and kinematic smearing functors, it allows detector effects to be specified within an analysis routine, customised to the exact phase-space and reconstruction working points of the analysis. A set of standard detector functions for the physics objects of Runs 1 and 2 of the ATLAS and CMS experiments is also provided. Finally, as jet substructure is an important class of physics observable usually considered to require an explicit detector simulation, we demonstrate that a smearing approach, tuned to available substructure data and implemented in Rivet, can accurately reproduce jet-structure biases observed by ATLAS.


title: Particle and entropy production in the Running Vacuum Universe

abstract: We study particle production and the corresponding entropy increase in the context of cosmology with dynamical vacuum. We focus on the particular form that has been called "running vacuum model" (RVM), which is known to furnish a successful description of the overall current observations at a competitive level with the concordance $\Lambda$CDM model. It also provides an elegant global explanation of the cosmic history from a non-singular initial state in the very early universe up to our days and further into the final de Sitter era. The model has no horizon problem and provides an alternative explanation for the early inflation and its graceful exit, as well as a powerful mechanism for generating the large entropy of the current universe. The energy-momentum tensor of matter is generally non-conserved in such context owing to particle creation or annihilation. We analyze general thermodynamical aspects of particle and entropy production in the RVM. We first study the entropy of particles in the comoving volume during the early universe and late universe. Then, in order to obtain a more physical interpretation, we pay attention to the entropy contribution from the cosmological apparent horizon, its interior and its surface. On combining the inner volume entropy with the entropy on the horizon, we elucidate with detailed calculations whether the evolution of the entropy of the RVM universe satisfies the Generalized Second Law of Thermodynamics. We find it is so and we prove that the essential reason for it is the existence of a positive cosmological constant.


title: Decoupling of Lattice and Orbital Degrees of Freedom in an Iron-Pnictide Superconductor

abstract: The interplay of structural and electronic phases in iron-based superconductors is a central theme in the search for the superconducting pairing mechanism. While electronic nematicity, defined as the breaking of four-fold symmetry triggered by electronic degrees of freedom, is competing with superconductivity, the effect of purely structural orthorhombic order is unexplored. Here, using x-ray diffraction (XRD), we reveal a new structural orthorhombic phase with an exceptionally high onset temperature ($T_\mathrm{ort} \sim 250$ K), which coexists with superconductivity ($T_\mathrm{c} = 25$ K), in an electron-doped iron-pnictide superconductor far from the underdoped region. Furthermore, our angle-resolved photoemission spectroscopy (ARPES) measurements demonstrate the absence of electronic nematic order as the driving mechanism, in contrast to other underdoped iron pnictides where nematicity is commonly found. Our results establish a new, high temperature phase in the phase diagram of iron-pnictide superconductors and impose strong constraints for the modeling of their superconducting pairing mechanism.


title: Stochastic Multistage Co-optimization of Generation and Transmission Expansion Planning

abstract: This paper describes a generation and transmission (G&T) expansion planning tool based on Benders decomposition and multistage stochastic optimization: (i) A MIP-based "investment module" determines a trial expansion plan; (ii) an SDDP-based "operation module" calculates the expected operation costs for the trial plan and Benders cuts to the investment module based on the expected marginal costs of the capacity constraints at the optimal solution. The integrated G&T planning approach is illustrated with a realistic planning study for the Bolivian power system.


title: The Explosion of Helium Stars Evolved With Mass Loss

abstract: Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored.


title: APEX: Adaptive Ext4 File System for Enhanced Data Recoverability in Edge Devices

abstract: Recently Edge Computing paradigm has gained significant popularity both in industry and academia. With its increased usage in real-life scenarios, security, privacy and integrity of data in such environments have become critical. Malicious deletion of mission-critical data due to ransomware, trojans and viruses has been a huge menace and recovering such lost data is an active field of research. As most of Edge computing devices have compute and storage limitations, difficult constraints arise in providing an optimal scheme for data protection. These devices mostly use Linux/Unix based operating systems. Hence, this work focuses on extending the Ext4 file system to APEX (Adaptive Ext4): a file system based on novel on-the-fly learning model that provides an Adaptive Recover-ability Aware file allocation platform for efficient post-deletion data recovery and therefore maintaining data integrity. Our recovery model and its lightweight implementation allow significant improvement in recover-ability of lost data with lower compute, space, time, and cost overheads compared to other methods. We demonstrate the effectiveness of APEX through a case study of overwriting surveillance videos by CryPy malware on Raspberry-Pi based Edge deployment and show 678% and 32% higher recovery than Ext4 and current state-of-the-art File Systems. We also evaluate the overhead characteristics and experimentally show that they are lower than other related works.


title: Low-frequency variability in massive stars: Core generation or surface phenomenon?

abstract: Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.


title: Constraining Star Formation Histories of Blue Galaxies using the Scatter between Stellar Mass and Halo Mass

abstract: We present constraints on the timescale of star formation variability and the correlation between star formation and host halo accretion histories of star-forming (SF) central galaxies from the scatter of the stellar-to-halo mass relation (SHMR). SF galaxies are found to have a tight relationship between their star formation rates and stellar masses on the so-called "star-forming sequence" (SFS), which characterizes their SF histories and $M_*$ growths. Meanwhile, observed SHMR constraints connect $M_*$ growth to halo accretion history. Combining these observed trends with a cosmological $N$-body simulation, we present flexible models that track the SF, $M_*$, and halo accretion histories of SF central galaxies at $z<1$ while reproducing the observed stellar mass function and SFS in SDSS DR7. We find that the scatter in SHMR at $M_h=10^{12}M_\odot$, $\sigma_{M_*|M_h=10^{12}M_\odot}$, is sensitive to the timescale of star formation variability, $t_{\rm duty}$, and the correlation coefficient, $r$, between SF and halo accretion histories: shorter $t_{\rm duty}$ and higher $r$ both result in tighter $\sigma_{M_*|M_h=10^{12}M_\odot}$. To reproduce a constant 0.2 dex scatter over $z=1-0$, our models require $t_{\rm duty}\leq1.5$ Gyr for $r=0.99$ or $r>0.6$ for $t_{\rm duty}=0.1$ Gyr. For $r\sim0.6$, as found in the literature, $t_{\rm duty}<0.2$ Gyr is necessary. Meanwhile, to reproduce the tightening of $\sigma_{M_*|M_h=10^{12}M_\odot}=0.35$ to 0.2 dex from $z=1-0$ in hydrodynamical simulations, our models require $t_{\rm duty}=0.1$ Gyr for $r>0.5$. Although, the lack of consensus on $\sigma_{M_*|M_h=10^{12}M_\odot}$ at $M_h=10^{12}M_\odot$ and at $z=1$ from observations and galaxy formation models remains the main bottleneck in precisely constraining $r$ and $t_{\rm duty}$, we demonstrate that SHMR can be used to constrain the SF and host halo accretion histories of SF central galaxies.


title: "Lagrangian Disks" in M-theory

abstract: While the study of bordered (pseudo-)holomorphic curves with boundary on Lagrangian submanifolds has a long history, a similar problem that involves (special) Lagrangian submanifolds with boundary on complex surfaces appears to be largely overlooked in both physics and math literature. We relate this problem to geometry of coassociative submanifolds in $G_2$ holonomy spaces and to $Spin(7)$ metrics on 8-manifolds with $T^2$ fibrations. As an application to physics, we propose a large class of brane models in type IIA string theory that generalize brane brick models on the one hand and 2d theories $T[M_4]$ on the other.


title: Searching for Solar Siblings in APOGEE and $Gaia$ DR2 with N-body Simulations

abstract: We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67.


title: A few StePS forward in unveiling the complexity of galaxy evolution: light-weighted stellar ages of intermediate redshift galaxies with WEAVE

abstract: The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. We produce spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assess our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise levels. Then, we analyze such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ~0.1 Gyr for stellar populations older than ~1.5 Gyr, pushing up to an age of ~1 Gyr for stellar populations older than ~5 Gyr. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.


title: Festina Lente: EFT Constraints from Charged Black Hole Evaporation in de Sitter

abstract: In the Swampland philosophy of constraining EFTs from black hole mechanics we study charged black hole evaporation in de Sitter space. We establish how the black hole mass and charge change over time due to both Hawking radiation and Schwinger pair production as a function of the masses and charges of the elementary particles in the theory. We find a lower bound on the mass of charged particles by demanding that large charged black holes evaporate back to empty de Sitter space, in accordance with the thermal picture of the de Sitter static patch. This bound is satisfied by the charged spectrum of the Standard Model. We discuss phenomenological implications for the cosmological hierarchy problem and inflation. Enforcing the thermal picture also leads to a heuristic remnant argument for the Weak Gravity Conjecture in de Sitter space, where the usual kinematic arguments do not work. We also comment on a possible relation between WGC and universal bounds on equilibration times. All in all, charged black holes in de Sitter should make haste to evaporate, but they should not rush it.


title: Refined Bounds on MeV-scale Thermal Dark Sectors from BBN and the CMB

abstract: New light states thermally coupled to the Standard Model plasma alter the expansion history of the Universe and impact the synthesis of the primordial light elements. In this work, we carry out an exhaustive and precise analysis of the implications of MeV-scale BSM particles in Big Bang Nucleosynthesis (BBN) and for Cosmic Microwave Background (CMB) observations. We find that, BBN observations set a lower bound on the thermal dark matter mass of $m_\chi > 0.4\,\text{MeV}$ at $2\sigma$. This bound is independent of the spin and number of internal degrees of freedom of the particle, of the annihilation being s-wave or p-wave, and of the annihilation final state. Furthermore, we show that current BBN plus CMB observations constrain purely electrophilic and neutrinophilic BSM species to have a mass, $m_\chi > 3.7\,\text{MeV}$ at $2\sigma$. We explore the reach of future BBN measurements and show that upcoming CMB missions should improve the bounds on light BSM thermal states to $m_\chi > (10-15)\,\text{MeV}$. Finally, we demonstrate that very light BSM species thermally coupled to the SM plasma are highly disfavoured by current cosmological observations.


title: Quantized Coulomb Branches, Monopole Bubbling and Wall-Crossing Phenomena in 3d $\mathcal{N}=4$ Theories

abstract: To study the quantized Coulomb branch of 3d $\mathcal{N}=4$ unitary SQCD theories, we propose a new method to compute correlators of monopole and Casimir operators that are inserted in the $\mathbb{R}\times\mathbb{R}^2_\epsilon$ Omega background. This method combines results from supersymmetric localization with inputs from the brane realisation of the correlators in type IIB string theory. The main challenge is the computation of the partition functions of certain Super-Matrix-Models (SMMs), which appear in the contribution of monopole bubbling sectors and are realised as the theory living on the D1 strings in the brane construction. We find that the non-commutativity arising in the monopole operator insertions is related to a wall-crossing phenomenon in the FI parameter space of the SMM. We illustrate our method in various examples and we provide explicit results for arbitrary correlators of non-bubbling bare monopole operators. We also discuss the realisation of the non-commutative product as a Moyal (star) product and use it to successfully test our results.


title: On the stability of bimetric structure formation

abstract: Bimetric gravity can reproduce the accelerated expansion of the Universe, without a cosmological constant. However, the stability of these solutions to linear perturbations has been questioned, suggesting exponential growth of structure in this approximation. We present a simple model of structure formation, for which an analytical solution is derived. The solution is well-behaved, showing that there is no physical instability with respect to these perturbations. The model can yield a growth of structure exhibiting measurable differences from $\Lambda$CDM.


title: A Dynamical Mass of $70 \pm 5$ Jupiter Masses for Gliese 229B, the First T Dwarf

abstract: We combine Keck/HIRES radial velocities, imaging with HiCIAO/Subaru and the Hubble Space Telescope, and absolute astrometry from Hipparcos and Gaia to measure a dynamical mass of $70 \pm 5$ Jupiter masses for the brown dwarf companion to Gl 229. Gl 229B was the first imaged brown dwarf to show clear signs of methane in its atmosphere. Cooling models have been used to estimate a mass in the range of 20 - 55 Jupiter masses, much lower than our measured value. We argue that our high dynamical mass is unlikely to be due to perturbations from additional unseen companions or to Gl 229B being itself a binary, and we find no evidence of a previously claimed radial velocity planet around Gl 229A. Future Gaia data releases will confirm the reliability of the absolute astrometry, though the data pass all quality checks in both Hipparcos and Gaia. Our dynamical mass implies a very old age for Gl 229, in some tension with kinematic and activity age indicators, and/or shortcomings in brown dwarf cooling models. Gl 229B joins a small but growing list of T dwarfs with masses approaching the minimum mass for core hydrogen ignition.


title: First star-forming structures in fuzzy cosmic filaments

abstract: In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.


title: Chiral instabilities & the onset of chiral turbulence in QED plasmas

abstract: We present a first principles study of chiral plasma instabilities and the onset of chiral turbulence in QED plasmas far from equilibrium. By performing classical-statistical lattice simulations of the microscopic theory, we show that the generation of strong helical magnetic fields from a helicity imbalance in the fermion sector proceeds via three distinct phases. During the initial linear instability regime the helicity imbalance of the fermion sector causes an exponential growth(damping) of magnetic field modes with right(left) handed polarization, for which we extract the characteristic growth (damping) rates. Secondary growth of unstable modes accelerates the helicity transfer from fermions to gauge fields and ultimately leads to the emergence of a self-similar scaling regime characteristic of decaying turbulence, where magnetic helicity is efficiently transferred to macroscopic length scales. Within this turbulent regime the evolution of magnetic helicity spectrum can be described by an infrared power-spectrum with spectral exponent $\kappa$ and dynamical scaling exponents $\alpha,\beta$, which we determine from our simulations.


title: Hinge states in a system of coupled Rashba layers

abstract: We consider a system of stacked tunnel-coupled two-dimensional electron- and hole-gas layers with Rashba spin-orbit interactions subjected to a staggered Zeeman field. The interplay of different intra-layer tunnel couplings results in a phase transition to a topological insulator phase in three dimensions hosting gapless surface states. The staggered Zeeman field further enriches the topological phase diagram by generating a second-order topological insulator phase hosting gapless hinge states. The emergence of the topological phases is proven analytically in the regime of small Zeeman field and confirmed by numerical simulations in the non-perturbative region of the phase diagram. The topological phases are stable against external perturbations and disorder.


title: RedMaPPer: Evolution and Mass Dependence of the Conditional Luminosity Functions of Red Galaxies in Galaxy Clusters

abstract: We characterize the luminosity distribution, halo mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 RedMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both central and satellite galaxies and the mass of their host halos, and show that this model is well-fit by the data. Using a larger galaxy cluster sample than previously employed in the literature, we find that the luminosities of central galaxies scale as $\langle \log L \rangle \propto A_L \log (M_{200b})$, with $A_L=0.39\pm0.04$, and that the scatter of the central--galaxy luminosity at fixed $M_{200b}$ ( $\sigma_{\log L|M}$) is $0.23 ^{+0.05}_{-0.04}$ dex, with the error bar including systematics due to miscentering of the cluster finder, photometry, and photometric redshift estimation. Our data prefers a positive correlation between the luminosity of central galaxies and the observed richness of clusters at a fixed halo mass, with an effective correlation coefficient $d_{\rm{eff}}=0.36^{+0.17}_{-0.16}$. The characteristic luminosity of satellites becomes dimmer from $z=0.3$ to $z=0.1$ by $\sim 20\%$ after accounting for passive evolution. We estimate the fraction of galaxy clusters where the brightest galaxy is not the central to be $P_{\rm{BNC}} \sim 20\%$. We discuss implications of these findings in the context of galaxy evolution and the galaxy--halo connection.


title: Gravitational Absorption Lines

abstract: We consider the gravitational analogue of Lyman-alpha absorption lines in astronomical spectroscopy. If Einstein gravity with minimally coupled matter is valid up to the Planck scale, quantum bound states absorb gravitons of a specific frequency with Planckian cross section, $\sigma_{\text{abs}} \approx l_p^2$. Consequently, one can show that gravitational absorption by bound states is inefficient in ordinary gravity. If observed, gravitational absorption lines would therefore constitute a powerful smoking gun of new exotic astrophysical bound states (near extremal bound states) or new gravitational physics, as well as give direct evidence of the quantized nature of the gravitational field. We provide, as an example of new gravitational physics near the Planck scale, a non-minimal coupling of the matter fields which breaks the equivalence principle on-shell. We lay out a model in which absorption lines in the primordial gravitational wave spectrum are produced as a consequence of this coupling.


title: Vertex algebras of CohFT-type

abstract: Representations of certain vertex algebras, here called of CohFT-type, can be used to construct vector bundles of coinvariants and conformal blocks on moduli spaces of stable curves [DGT2]. We show that such bundles define semisimple cohomological field theories. As an application, we give an expression for their total Chern character in terms of the fusion rules, following the approach and computation in [MOPPZ] for bundles given by integrable modules over affine Lie algebras. It follows that the Chern classes are tautological. Examples and open problems are discussed.


title: Efficient Quantum Walk Circuits for Metropolis-Hastings Algorithm

abstract: We present a detailed circuit implementation of Szegedy's quantization of the Metropolis-Hastings walk. This quantum walk is usually defined with respect to an oracle. We find that a direct implementation of this oracle requires costly arithmetic operations and thus reformulate the quantum walk in a way that circumvents the implementation of that specific oracle and which closely follows the classical Metropolis-Hastings walk. We also present heuristic quantum algorithms that use the quantum walk in the context of discrete optimization problems and numerically study their performances. Our numerical results indicate polynomial quantum speedups in heuristic settings.


title: The inner view of NGC 1052 using multiple X-ray observations

abstract: In this paper, we make a multi-epoch analysis of NGC 1052, one of the prototypical LLAGN, using XMM-Newton, Suzaku and NuSTAR observations, taken from 2001 to 2017. This is the first time that results from NuSTAR observations are reported for NGC 1052. On the technical aspects, we found a wavelength-dependent calibration issue between simultaneous XMM-Newton and NuSTAR spectra. It is described by a change on the photon index of $\rm{ \Gamma_{NuSTAR}- \Gamma_{XMM-Newton}=0.17\pm0.04}$. We use ancillary Chandra data to decontaminate the nuclear spectrum from circumnuclear contributors. We find that two baseline models can fit the broad (0.5-50 keV) X-ray spectrum of the source. One consists of a power-law like continuum which is absorbed by a uniform absorber, reflection from neutral material, and a separate power-law component in the soft band. The second model presents a clumpy absorber. The reflection component is still present, but not the soft band power-law. Instead, absorption by a warm absorber is necessary to fit the spectra. This is the first time that a reflection component is established in this object, thanks to high energy data from NuSTAR. This component is constant in flux and shape, supporting the idea that is produced away from the central source (probably the torus). We find flux, spectral slope and absorption variations on timescales of months to years. We also find that a patchy-absober can explain the behaviour of this source better as it is $\sim$ 200 times more likely than the uniform absober while it yields to smaller intrinsic variations.


title: Merger induced clump formation in distant infrared luminous starburst galaxies

abstract: While the formation of stellar clumps in distant galaxies is usually attributed to gravitational violent disk instabilities, we show here that major mergers also represent a competitive mechanism to form bright clumps. Using ~0.1'' resolution ACS F814W images in the entire COSMOS field, we measure the fraction of clumpy emission in 109 main sequence (MS) and 79 Herschel-detected starbursts (off-MS) galaxies at 0.5 < z < 0.9, representative of normal versus merger induced star-forming activity, respectively. We additionally identify merger samples from visual inspection and from Gini-M20 morphological parameters. Regardless of the merger criteria adopted, the clumpiness distribution of merging systems is different from that of normal isolated disks at > 99.5 % confidence level, with the former reaching higher clumpiness values, up to 20 % of the total galaxy emission. We confirm the merger induced clumpiness enhancement with novel hydrodynamical simulations of colliding galaxies with gas fractions typical of z ~ 0.7. Multi-wavelength images of three starbursts in the CANDELS field support the young nature of clumps, which are likely merger products rather than older pre-existing structures. Finally, for a subset of 19 starbursts with existing near-IR rest frame spectroscopy, we find that the clumpiness is mildly anti-correlated with the merger phase, decreasing towards final coalescence. Our result can explain recent ALMA detections of clumps in hyperluminous high-z starbursts, while normal objects are smooth. This work raises a question on the role of mergers on the origin of clumps in high redshift galaxies in general.


title: Symmetries for a High Level Neural Decoder on the Toric Code

abstract: Surface codes are a promising method of quantum error correction and the basis of many proposed quantum computation implementations. However, their efficient decoding is still not fully explored. Recently, approaches based on machine learning techniques have been proposed by Torlai and Melko as well as Varsamopoulos et al. In these approaches, a so called high level decoder is used to post-correct an underlying decoder by correcting logical errors. A significant problem is that these methods require large amounts of training data even for relatively small code distances. The above-mentioned methods were tested on the rotated surface code which encodes one logical qubit. Here, we show that they are viable even for the toric surface code which encodes two logical qubits. Furthermore, we explain how symmetries of the toric code can be exploited to reduce the amount of training data that is required to obtain good decoding results. Finally, we compare different underlying decoders and show that the accuracy of high level decoding noticeably depends on the quality of the underlying decoder in the realistic case of imperfect training.


title: Harnessing the Power of Infinitely Wide Deep Nets on Small-data Tasks

abstract: Recent research shows that the following two models are equivalent: (a) infinitely wide neural networks (NNs) trained under l2 loss by gradient descent with infinitesimally small learning rate (b) kernel regression with respect to so-called Neural Tangent Kernels (NTKs) (Jacot et al., 2018). An efficient algorithm to compute the NTK, as well as its convolutional counterparts, appears in Arora et al. (2019a), which allowed studying performance of infinitely wide nets on datasets like CIFAR-10. However, super-quadratic running time of kernel methods makes them best suited for small-data tasks. We report results suggesting neural tangent kernels perform strongly on low-data tasks. 1. On a standard testbed of classification/regression tasks from the UCI database, NTK SVM beats the previous gold standard, Random Forests (RF), and also the corresponding finite nets. 2. On CIFAR-10 with 10 - 640 training samples, Convolutional NTK consistently beats ResNet-34 by 1% - 3%. 3. On VOC07 testbed for few-shot image classification tasks on ImageNet with transfer learning (Goyal et al., 2019), replacing the linear SVM currently used with a Convolutional NTK SVM consistently improves performance. 4. Comparing the performance of NTK with the finite-width net it was derived from, NTK behavior starts at lower net widths than suggested by theoretical analysis(Arora et al., 2019a). NTK's efficacy may trace to lower variance of output.