title: Reflected Light Phase Curves in the TESS Era

abstract: The reflected light signal from a planet throughout its orbit is a powerful probe of a planet's atmospheric properties. There are a number of planets that are amenable to reflected light phase curve studies with present and future space-based instrumentation and here we assess our ability to characterize these worlds. Using simulated TESS populations we identify the Nine, a set of archetypal exoplanets with the potential to be bright in reflected light, because of their radii and proximity to their star, while still being cool enough to have minimal thermal contamination at optical wavelengths. For each planet we compute albedo spectra for several cloud and atmosphere assumptions (e.g. thermochemical equilibrium, solar composition). We find that in the TESS bandpass the estimated contrast at optical wavelengths is typically <10 ppm except for the brightest, largest, or closest in planets with the highest lofted clouds where contrast can reach a few tens of ppm. Meanwhile, in a bluer bandpass (0.3--0.5 microns) the estimated contrast can be as high as 150 ppm but typically 10-50 ppm. In the temperature range of interest, planets with the highest, most extensive cloud decks are generally darker at bluer wavelengths than cloudless planets because of the low single scattering albedos of their primary condensate constituents. Our models suggest that Neptune-sized planets with relatively low insolation and small semi-major axes are the most conducive to reflected light phase curve studies in TESS.


title: Emergence of Object Segmentation in Perturbed Generative Models

abstract: We introduce a novel framework to build a model that can learn how to segment objects from a collection of images without any human annotation. Our method builds on the observation that the location of object segments can be perturbed locally relative to a given background without affecting the realism of a scene. Our approach is to first train a generative model of a layered scene. The layered representation consists of a background image, a foreground image and the mask of the foreground. A composite image is then obtained by overlaying the masked foreground image onto the background. The generative model is trained in an adversarial fashion against a discriminator, which forces the generative model to produce realistic composite images. To force the generator to learn a representation where the foreground layer corresponds to an object, we perturb the output of the generative model by introducing a random shift of both the foreground image and mask relative to the background. Because the generator is unaware of the shift before computing its output, it must produce layered representations that are realistic for any such random perturbation. Finally, we learn to segment an image by defining an autoencoder consisting of an encoder, which we train, and the pre-trained generator as the decoder, which we freeze. The encoder maps an image to a feature vector, which is fed as input to the generator to give a composite image matching the original input image. Because the generator outputs an explicit layered representation of the scene, the encoder learns to detect and segment objects. We demonstrate this framework on real images of several object categories.


title: CCM property and Lyubeznik numbers under Groebner deformations

abstract: In this note we draw some interesting consequences of the recent results on squarefree Groebner degenerations obtained by Conca and the second author.


title: Graph Learning Network: A Structure Learning Algorithm

abstract: Recently, graph neural networks (GNNs) have proved to be suitable in tasks on unstructured data. Particularly in tasks as community detection, node classification, and link prediction. However, most GNN models still operate with static relationships. We propose the Graph Learning Network (GLN), a simple yet effective process to learn node embeddings and structure prediction functions. Our model uses graph convolutions to propose expected node features, and predict the best structure based on them. We repeat these steps recursively to enhance the prediction and the embeddings.


title: Evaluating structure learning algorithms with a balanced scoring function

abstract: Several structure learning algorithms have been proposed towards discovering causal or Bayesian Network (BN) graphs. The validity of these algorithms tends to be evaluated by assessing the relationship between the learnt and the ground truth graph. However, there is no agreed scoring metric to determine this relationship. Moreover, this paper shows that some of the commonly used metrics tend to be biased in favour of graphs that minimise edges. While graphs that are less complex are desirable, some of the metrics favour underfitted graphs, thereby encouraging limited propagation of evidence. This paper proposes the Balanced Scoring Function (BSF) that eliminates this bias by adjusting the reward function based on the difficulty of discovering an edge, or no edge, proportional to their occurrence rate in the ground truth graph. The BSF score can be used in conjunction with other traditional metrics to provide an alternative and unbiased assessment about the capability of a structure learning algorithm in discovering causal or BN graphs.


title: Structured Monte Carlo Sampling for Nonisotropic Distributions via Determinantal Point Processes

abstract: We propose a new class of structured methods for Monte Carlo (MC) sampling, called DPPMC, designed for high-dimensional nonisotropic distributions where samples are correlated to reduce the variance of the estimator via determinantal point processes. We successfully apply DPPMCs to problems involving nonisotropic distributions arising in guided evolution strategy (GES) methods for RL, CMA-ES techniques and trust region algorithms for blackbox optimization, improving state-of-the-art in all these settings. In particular, we show that DPPMCs drastically improve exploration profiles of the existing evolution strategy algorithms. We further confirm our results, analyzing random feature map estimators for Gaussian mixture kernels. We provide theoretical justification of our empirical results, showing a connection between DPPMCs and structured orthogonal MC methods for isotropic distributions.


title: Magnetic phases of orbital bipartite optical lattices

abstract: In the Hamburg cold atom experiment with orbital states in an optical lattice, $s$- and $p$-orbital atomic states hybridize between neighbouring sites. In this work we show how this alternation of sites hosting $s$- and $p$-orbital states gives rise to a plethora of different magnetic phases, quantum and classical. We focus on phases whose properties derive from frustration originating from a competition between nearest and next nearest neighbouring exchange interactions. The physics of the Mott insulating phase with unit filling is described by an effective spin-1/2 Hamiltonian showing great similarities with the $J_1$-$J_2$ model. Based on the knowledge of the $J_1$-$J_2$ model, together with numerical simulations, we discuss the possibility of realising a quantum spin liquid phase in the present optical lattice system. In the superfluid regime we consider the parameter regime where the $s$-orbital states can be adiabatically eliminated to give an effective model for the $p$-orbital atoms. At the mean-field level we derive a generalized classical $XY$ model, and show that it may support maximum frustration. When quantum fluctuations can be disregarded, the ground state is expected to be a spin glass. Even with quantum fluctuations present it has been debated whether a spin liquid may persist at the point of full frustration.


title: Picky Partners: The Pairing of Component Masses in Binary Black Hole Mergers

abstract: We examine the relationship between individual black hole (BH) masses in merging binary black hole (BBH) systems. Analyzing the ten BBH detections from LIGO/Virgo's first two observing runs, we find that the masses of the component BHs comprising each binary are unlikely to be randomly drawn from the same underlying distribution. Instead, the two BHs of a given binary prefer to be of comparable mass. We show that it is $\sim 5$ times more likely that the component BHs in a given binary are always equal (to within 5\%) than that they are randomly paired. If we assume that the probability of a merger between two BHs scales with the mass ratio $q$ as $q^\beta$, so that $\beta=0$ corresponds to random pairings, we find $\beta>0$ is favored at credibility 0.987. By modeling the mass distribution, we find that the median mass ratio is $q_{50\%} = 0.91^{+0.05}_{-0.17}$ at 90\% credibility. While the pairing between BHs depends on their mass ratio, we find no evidence that it depends on the total mass of the system: it is $\sim 6$ times more likely that the pairing depends purely on the mass ratio than on the total mass. We predict that 99\% of BBHs detected by LIGO/Virgo will have mass ratios $q > 0.5$. We conclude that merging black holes do not form random pairings; instead they are selective about their partners, preferring to mate with black holes of a similar mass. The details of these selective pairings provide insight into the underlying formation channels of merging binaries.


title: Gate Controlled Anomalous Phase Shift in Al/InAs Josephson Junctions

abstract: In a standard Josephson junction the current is zero when the phase difference between the superconducting leads is zero. This condition is protected by parity and time-reversal symmetries. However, the combined presence of spin-orbit coupling and magnetic field breaks these symmetries and can lead to a finite supercurrent even when the phase difference is zero. This is the so called anomalous Josephson effect -- the hallmark effect of superconducting spintronics --and can be characterized by the corresponding anomalous phase shift ($\phi_0$). We report the observation of a tunable anomalous Josephson effect in InAs/Al Josephson junctions measured via a superconducting quantum interference device (SQUID). By gate controlling the density of InAs we are able to tune the spin-orbit coupling of the Josephson junction by more than one order of magnitude. This gives us the ability to tune $\phi_0$, and opens several new opportunities for superconducting spintronics, and new possibilities for realizing and characterizing topological superconductivity.


title: Strain engineering in Ge/GeSn core/shell nanowires

abstract: Strain engineering in Sn-rich group IV semiconductors is a key enabling factor to exploit the direct band gap at mid-infrared wavelengths. Here, we investigate the effect of strain on the growth of GeSn alloys in a Ge/GeSn core/shell nanowire geometry. Incorporation of Sn content in the 10-20 at.% range is achieved with Ge core diameters ranging from 50nm to 100nm. While the smaller cores lead to the formation of a regular and homogeneous GeSn shell, larger cores lead to the formation of multi-faceted sidewalls and broadened segregation domains, inducing the nucleation of defects. This behavior is rationalized in terms of the different residual strain, as obtained by realistic finite element method simulations. The extended analysis of the strain relaxation as a function of core and shell sizes, in comparison with the conventional planar geometry, provides a deeper understanding of the role of strain in the epitaxy of metastable GeSn semiconductors.


title: Efficient and tunable Aharonov-Bohm quantum heat engine

abstract: We propose a quantum heat engine based on an Aharonov-Bohm interferometer in a two-terminal geometry, and investigate its thermoelectric performances in the linear response regime. Sizeable thermopower (up to $\sim 0.3\,\text{mV}$/K) as well as $ZT$ values largely exceeding unity can be achieved by simply adjusting parameters of the setup and temperature bias across the interferometer leading to thermal efficiency at maximum power approaching $30\%$ of the Carnot limit. This is close to the optimal efficiency at maximum power achievable for a two-terminal heat engine. Changing the magnetic flux, the asymmetry of the structure, a side-gate bias voltage through a capacitively-coupled electrode and the transmission of the T-junctions connecting the AB ring to the contacts allows to finely tune the operation of the quantum heat engine. The exploration of the parameters' space demonstrates that the high performances of the Aharonov-Bohm two-terminal device as a quantum heat engine are stable over a wide range of temperatures and length imbalances, promising towards experimental realization.


title: Regret Bounds for Thompson Sampling in Episodic Restless Bandit Problems

abstract: Restless bandit problems are instances of non-stationary multi-armed bandits. These problems have been studied well from the optimization perspective, where the goal is to efficiently find a near-optimal policy when system parameters are known. However, very few papers adopt a learning perspective, where the parameters are unknown. In this paper, we analyze the performance of Thompson sampling in episodic restless bandits with unknown parameters. We consider a general policy map to define our competitor and prove an $\tilde{\mathcal{O}}(\sqrt{T})$ Bayesian regret bound. Our competitor is flexible enough to represent various benchmarks including the best fixed action policy, the optimal policy, the Whittle index policy, or the myopic policy. We also present empirical results that support our theoretical findings.


title: End-to-end capacities of a quantum communication network

abstract: In quantum mechanics, a fundamental law prevents quantum communications to simultaneously achieve high rates and long distances. This limitation is well known for point-to-point protocols, where two parties are directly connected by a quantum channel, but not yet fully understood in protocols with quantum repeaters. Here we solve this problem bounding the ultimate rates for transmitting quantum information, entanglement and secret keys via quantum repeaters. We derive single-letter upper bounds for the end-to-end capacities achievable by the most general (adaptive) protocols of quantum and private communication, from a single repeater chain to an arbitrarily-complex quantum network, where systems may be routed through single or multiple paths. We analytically establish these capacities under fundamental noise models, including bosonic loss which is the most important for optical communications. In this way, our results provide the ultimate benchmarks for testing the optimal performance of repeater-assisted quantum communications.

categories: quant-ph cond-mat.other math-ph math.MP physics.optics


title: On the uniqueness of solutions to quadratic BSDEs with non-convex generators and unbounded terminal conditions

abstract: We prove a uniqueness result of the unbounded solution for a quadratic backward stochastic differential equation whose terminal condition is unbounded and whose generator $g$ may be non-Lipschitz continuous in the state variable $y$, non-convex (non-concave) in the state variable $z$, and instead satisfies a strictly quadratic condition and an additional assumption. The key observation is that if the generator is strictly quadratic, then the quadratic variation of the first component of the solution admits an exponential moment. Typically, a Lipschitz perturbation of some convex (concave) function satisfies the additional assumption mentioned above. This generalizes some results obtained in \cite{BriandHu2006PTRF} and \cite{BriandHu2008PTRF}.


title: The (Non-)Utility of Structural Features in BiLSTM-based Dependency Parsers

abstract: Classical non-neural dependency parsers put considerable effort on the design of feature functions. Especially, they benefit from information coming from structural features, such as features drawn from neighboring tokens in the dependency tree. In contrast, their BiLSTM-based successors achieve state-of-the-art performance without explicit information about the structural context. In this paper we aim to answer the question: How much structural context are the BiLSTM representations able to capture implicitly? We show that features drawn from partial subtrees become redundant when the BiLSTMs are used. We provide a deep insight into information flow in transition- and graph-based neural architectures to demonstrate where the implicit information comes from when the parsers make their decisions. Finally, with model ablations we demonstrate that the structural context is not only present in the models, but it significantly influences their performance.


title: Bounds for multi-end communication over quantum networks

abstract: Quantum and private communications are affected by a fundamental limitation which severely restricts the optimal rates that are achievable by two distant parties. To overcome this problem, one needs to introduce quantum repeaters and, more generally, quantum communication networks. Within a quantum network, other problems and features may appear when we move from the basic unicast setting of single-sender/single-receiver to more complex multi-end scenarios, where multiple senders and multiple receivers simultaneously use the network to communicate. Assuming various configurations, including multiple-unicast, multicast, and multiple-multicast communication, we bound the optimal rates for transmitting quantum information, distributing entanglement, or generating secret keys in quantum networks connected by arbitrary quantum channels. These bounds cannot be surpassed by the most general adaptive protocols of quantum network communication.


title: Entropic Regularisation of Robust Optimal Transport

abstract: Grogan et al [11,12] have recently proposed a solution to colour transfer by minimising the Euclidean distance L2 between two probability density functions capturing the colour distributions of two images (palette and target). It was shown to be very competitive to alternative solutions based on Optimal Transport for colour transfer. We show that in fact Grogan et al's formulation can also be understood as a new robust Optimal Transport based framework with entropy regularisation over marginals.


title: Nonvolatile Spintronic Memory Cells for Neural Networks

abstract: A new spintronic nonvolatile memory cell analogous to 1T DRAM with non-destructive read is proposed. The cells can be used as neural computing units. A dual-circuit neural network architecture is proposed to leverage these devices against the complex operations involved in convolutional networks. Simulations based on HSPICE and Matlab were performed to study the performance of this architecture when classifying images as well as the effect of varying the size and stability of the nanomagnets. The spintronic cells outperform a purely charge-based implementation of the same network, consuming about 100 pJ total per image processed.


title: Resolvent-based study of compressibility effects on supersonic turbulent boundary layers

abstract: The resolvent formulation of McKeon & Sharma (2010) is applied to supersonic turbulent boundary layers to study the validity of Morkovin's hypothesis, which postulates that high-speed turbulence structures in zero pressure-gradient turbulent boundary layers remain largely the same as its incompressible counterpart. Supersonic zero-pressure-gradient turbulent boundary layers with adiabatic wall boundary conditions at Mach numbers ranging from 2 to 4 are considered. Resolvent analysis highlights two distinct regions of the supersonic turbulent boundary layer in the wave parameter space: the relatively supersonic region and the relatively subsonic region. In the relatively supersonic region, where the flow is supersonic relative to the freestream, resolvent modes display structures consistent with Mach wave radiation that are absent in the incompressible regime. In the relatively subsonic region, we show that the low-rank approximation of the resolvent operator is an effective approximation of the full system and that the response modes predicted by the model exhibit universal and geometrically self-similar behaviour via a transformation given by the semi-local scaling. Moreover, with the semi-local scaling, we show that the resolvent modes follow the same scaling law as their incompressible counterparts in this region, which has implications for modelling and the prediction of turbulent high-speed wall-bounded flows. We also show that the thermodynamic variables exhibit similar mode shapes to the streamwise velocity modes, supporting the strong Reynolds analogy. Finally, we demonstrate that the principal resolvent modes can be used to capture the energy distribution between momentum and thermodynamic fluctuations.


title: What Makes Training Multi-Modal Classification Networks Hard?

abstract: Consider end-to-end training of a multi-modal vs. a single-modal network on a task with multiple input modalities: the multi-modal network receives more information, so it should match or outperform its single-modal counterpart. In our experiments, however, we observe the opposite: the best single-modal network always outperforms the multi-modal network. This observation is consistent across different combinations of modalities and on different tasks and benchmarks. This paper identifies two main causes for this performance drop: first, multi-modal networks are often prone to overfitting due to increased capacity. Second, different modalities overfit and generalize at different rates, so training them jointly with a single optimization strategy is sub-optimal. We address these two problems with a technique we call Gradient Blending, which computes an optimal blend of modalities based on their overfitting behavior. We demonstrate that Gradient Blending outperforms widely-used baselines for avoiding overfitting and achieves state-of-the-art accuracy on various tasks including human action recognition, ego-centric action recognition, and acoustic event detection.


title: Dark Energy Survey Year 1 Results: Wide field mass maps via forward fitting in harmonic space

abstract: We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser-Squires method, the technique is then applied to data from the Dark Energy Survey Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations.


title: A survey of Object Classification and Detection based on 2D/3D data

abstract: Recently, by using deep neural network based algorithms, object classification, detection and semantic segmentation solutions are significantly improved. However, one challenge for 2D image-based systems is that they cannot provide accurate 3D location information. This is critical for location sensitive applications such as autonomous driving and robot navigation. On the other hand, 3D methods, such as RGB-D and RGB-LiDAR based systems, can provide solutions that significantly improve the RGB only approaches. That is why this is an interesting research area for both industry and academia. Compared with 2D image-based systems, 3D-based systems are more complicated due to the following five reasons: 1) Data representation itself is more complicated. 3D images can be represented by point clouds, meshes, volumes. 2D images have pixel grid representations. 2) The computation and memory resource requirement is higher as an extra dimension is added. 3) Different distribution of the objects and difference in scene areas between indoor and outdoor make one unified framework hard to achieve. 4) 3D data, especially for the outdoor scenario, is sparse compared with the dense 2D images which makes the detection task more challenging. Finally, large size labelled datasets, which are extremely important for supervised based algorithms, are still under construction compared with well-built 2D datasets such as ImageNet. Based on challenges listed above, the described systems are organized by application scenarios, data representation methods and main tasks addressed. At the same time, critical 2D based systems which greatly influence the 3D ones are also introduced to show the connection between them.


title: Mean-dependent nonstationary spatial models

abstract: Nonstationarity is a major challenge in analyzing spatial data. For example, daily precipitation measurements may have increased variability and decreased spatial smoothness in areas with high mean rainfall. Common nonstationary covariance models introduce parameters specific to each location, giving a highly-parameterized model which is difficult to fit. We develop a nonstationary spatial model that uses the mean to determine the covariance in a region, resulting in a far simpler, albeit more specialized, model. We explore inferential and predictive properties of the model under various simulated data situations. We show that this model in certain circumstances improves predictions compared to a standard stationary spatial model. We further propose a computationally efficient approximation that has comparable predictive accuracy. We also develop a test for nonstationary data and show it reliably identifies nonstationarity. We apply these methods to daily precipitation in Puerto Rico.


title: The Global Magneto-Ionic Medium Survey: Polarimetry of the Southern Sky from 300 to 480 MHz

abstract: Much data on the Galactic polarized radio emission has been gathered in the last five decades. All-sky surveys have been made, but only in narrow, widely spaced frequency bands, and the data are inadequate for the characterization of Faraday rotation, the main determinant of the appearance of the polarized radio sky at decimetre wavelengths. We describe a survey of the polarized radio emission from the Southern sky, aiming to characterize the magneto-ionic medium, particularly the strength and configuration of the magnetic field. This work is part of the Global Magneto-Ionic Medium Survey (GMIMS). We have designed and built a feed and receiver covering the band 300 to 900 MHz for the CSIRO Parkes 64-m Telescope. We have surveyed the entire sky between declinations -90 and +20 degrees. We present data covering 300 to 480 MHz with angular resolution 81' to 45'. The survey intensity scale is absolutely calibrated, based on measurements of resistors at known temperatures and on an assumed flux density and spectral index for Taurus A. Data are presented as brightness temperatures. We have applied Rotation Measure Synthesis to the data to obtain a Faraday depth cube of resolution 5.9 radians per metre squared, sensitivity of 60 mK of polarized intensity, and angular resolution 1.35 degrees. The data presented in this paper are available at the Canadian Astronomy Data Centre.


title: Learning Representations by Humans, for Humans

abstract: When machine predictors can achieve higher performance than the human decision-makers they support, improving the performance of human decision-makers is often conflated with improving machine accuracy. Here we propose a framework to directly support human decision-making, in which the role of machines is to reframe problems rather than to prescribe actions through prediction. Inspired by the success of representation learning in improving performance of machine predictors, our framework learns human-facing representations optimized for human performance. This "Mind Composed with Machine" framework incorporates a human decision-making model directly into the representation learning paradigm and is trained with a novel human-in-the-loop training procedure. We empirically demonstrate the successful application of the framework to various tasks and representational forms.


title: The role of bot squads in the political propaganda on Twitter

abstract: Social Media are nowadays the privileged channel for information spreading and news checking. Unexpectedly for most of the users, automated accounts, also known as social bots, contribute more and more to this process of news spreading. Using Twitter as a benchmark, we consider the traffic exchanged, over one month of observation, on a specific topic, namely the migration flux from Northern Africa to Italy. We measure the significant traffic of tweets only, by implementing an entropy-based null model that discounts the activity of users and the virality of tweets. Results show that social bots play a central role in the exchange of significant content. Indeed, not only the strongest hubs have a number of bots among their followers higher than expected, but furthermore a group of them, that can be assigned to the same political tendency, share a common set of bots as followers. The retwitting activity of such automated accounts amplifies the presence on the platform of the hubs' messages.


title: Choosing Transfer Languages for Cross-Lingual Learning

abstract: Cross-lingual transfer, where a high-resource transfer language is used to improve the accuracy of a low-resource task language, is now an invaluable tool for improving performance of natural language processing (NLP) on low-resource languages. However, given a particular task language, it is not clear which language to transfer from, and the standard strategy is to select languages based on ad hoc criteria, usually the intuition of the experimenter. Since a large number of features contribute to the success of cross-lingual transfer (including phylogenetic similarity, typological properties, lexical overlap, or size of available data), even the most enlightened experimenter rarely considers all these factors for the particular task at hand. In this paper, we consider this task of automatically selecting optimal transfer languages as a ranking problem, and build models that consider the aforementioned features to perform this prediction. In experiments on representative NLP tasks, we demonstrate that our model predicts good transfer languages much better than ad hoc baselines considering single features in isolation, and glean insights on what features are most informative for each different NLP tasks, which may inform future ad hoc selection even without use of our method. Code, data, and pre-trained models are available at https://github.com/neulab/langrank


title: Quantum Gravity from Timelike Liouville theory

abstract: A proper definition of the path integral of quantum gravity has been a long-standing puzzle because the Weyl factor of the Euclidean metric has a wrong-sign kinetic term. We propose a definition of two-dimensional Liouville quantum gravity with cosmological constant using conformal bootstrap for the timelike Liouville theory coupled to supercritical matter. We prove a no-ghost theorem for the states in the BRST cohomology. We show that the four-point function constructed by gluing the timelike Liouville three-point functions is well defined and crossing symmetric (numerically) for external Liouville energies corresponding to \textit{all} physical states in the BRST cohomology with the choice of the Ribault-Santachiara contour for the internal energy.


title: A decomposition of the Jacobian of a Humbert-Edge curve

abstract: A \textit{Humbert-Edge curve of type} $n$ is a non-degenerate smooth complete intersection of $n-1$ diagonal quadrics. Such a curve has an interesting geometry since it has a natural action of the group $(\mathbb{Z}/2\mathbb{Z})^n$. We present here a decomposition of its Jacobian variety as a product of Prym-Tyurin varieties, and we compute the kernel of the corresponding isogeny.


title: Slow quenches in Chern insulator ribbons

abstract: We investigate slow quenches in Chern insulators in ribbon geometry. We consider the Qi-Wu-Zhang model and slowly ramp the parameters (large time of the quench $\tau$) from a non-topological (Chern number = 0) to a topological regime (Chern number $\ne$ 0). In contrast to the Haldane model considered in [L. Privitera and G. E. Santoro, Phys. Rev. B 93, 241406(R) (2016)] earlier, the in-gap state degeneracy point is pinned to an inversion symmetric momentum, which changes the behavior drastically. The density of excitations in the in-gap states scales with the quench time as $\tau^{-1/2}$ as the ramp becomes slow, and the Kibble-Zurek mechanism applies. Despite the slower scaling of the density of in-gap excitations with $\tau$, the Hall conductance after the quench deviates from that of the ground state of the final Hamiltonian by an amount that drops as $\tau^{-1}$.


title: Vector-Valued Graph Trend Filtering with Non-Convex Penalties

abstract: This work studies the denoising of piecewise smooth graph signals that exhibit inhomogeneous levels of smoothness over a graph, where the value at each node can be vector-valued. We extend the graph trend filtering framework to denoising vector-valued graph signals with a family of non-convex regularizers, which exhibit superior recovery performance over existing convex regularizers. Using an oracle inequality, we establish the statistical error rates of first-order stationary points of the proposed non-convex method for generic graphs. Furthermore, we present an ADMM-based algorithm to solve the proposed method and establish its convergence. Numerical experiments are conducted on both synthetic and real-world data for denoising, support recovery, event detection, and semi-supervised classification.


title: An [O III] search for extended emission around AGN with H I mapping: a distant cloud ionized by Mkn 1

abstract: Motivated by the discovery of large-scale ionized clouds around AGN host galaxies, and particularly the large fraction of those which are consistent with photoionized gaseous tidal debris, we have searched for [O III] emission over wide fields around a set of Seyfert galaxies previously mapped in H I, many of which show extended gas features. The detection threshold was designed to reach mean emission-line surface brightness 10 times fainter than seen in such AGN-ionized clouds as Hanny's Voorwerp, so that similar structures at larger distances (and ages) could be detected. Of 24 Seyfert galaxies, we find one extended emission feature, a discrete cloud projected 12 kpc from the center of Mkn 1 and spanning a transverse extent of 8 kpc. Optical spectroscopy of the Mkn 1 cloud confirms its redshift association with the Mkn 1- NGC 451 galaxy pair, shows it to closely match the kinematics of nearby H I, and reveals emission-line ratios requiring photoionization by the AGN at roughly the direct observed luminosity of the nucleus. Given the small fraction of H I features with detected [O III] emission, we constrain the typical opening angle of ionization cones in Seyfert galaxies to be of order 20 deg, if active episodes are long compared to the light-travel times involved. An appendix presents a derivation of an analytical expression for the probability of intersection of a cone with randomly oriented arcs, approximating the geometry of H I clouds and tails exposed to ionization cones. For the entire sample, the full opening angle of bicones must be <20 deg if the AGN are continuously bright for scales longer than the light-travel times across the H I structures. Since many ionization cones are observed to be much broader than this, our low detection fraction may add to evidence for the ubiquity of strong variations in AGN luminosity on scales 10,000-100,000 years.


title: Tweeting your Destiny: Profiling Users in the Twitter Landscape around an Online Game

abstract: Social media has become a major communication channel for communities centered around video games. Consequently, social media offers a rich data source to study online communities and the discussions evolving around games. Towards this end, we explore a large-scale dataset consisting of over 1 million tweets related to the online multiplayer shooter Destiny and spanning a time period of about 14 months using unsupervised clustering and topic modelling. Furthermore, we correlate Twitter activity of over 3,000 players with their playtime. Our results contribute to the understanding of online player communities by identifying distinct player groups with respect to their Twitter characteristics, describing subgroups within the Destiny community, and uncovering broad topics of community interest.


title: A New Approach to Lossy Network Compression of a Tuple of Correlated Multivariate Gaussian RVs

abstract: The classical Gray and Wyner source coding for a simple network for sources that generate a tuple of multivariate, correlated Gaussian random variables $(Y_1,Y_2)$ is re-examined using the geometric approach of Gaussian random variables, and the weak stochastic realization of correlated Gaussian random variables. New results are: (1) The formulation, methods and algorithms to parametrize all random variables $W : \Omega \rightarrow {\mathbb R}^n $ which make the two components of the tuple $(Y_1,Y_2)$ conditionally independent, according to the weak stochastic realization of $(Y_1, Y_2)$. (2) The transformation of random variables $(Y_1,Y_2)$ via non-singular transformations $(S_1,S_2)$, into their canonical variable form. (3) A formula for Wyner's lossy common information for joint decoding with mean-square error distortions. (4) The methods are shown to be of fundamental importance to the parametrization of the lossy rate region of the Gray and Wyner source coding problem, and the calculation of the smallest common message rate $R_0$ on the Gray and Wyner source problem, when the sum rate $R_0+R_1+R_2$ is arbitrary close to the joint rate distortion function $R_{Y_1, Y_2}(\Delta_1, \Delta_2)$ of joint decoding. The methods and algorithms may be applicable to other problems of multi-user communication, such as, the multiple access channel, etc. The discussion is largely self-contained and proceeds from first principles; basic concepts of weak stochastic realization theory of multivariate correlated Gaussian random variables are reviewed, while certain results are developed to meet the requirement of results (1)-(4).


title: Inference in latent factor regression with clusterable features

abstract: Regression models, in which the observed features $X \in \R^p$ and the response $Y \in \R$ depend, jointly, on a lower dimensional, unobserved, latent vector $Z \in \R^K$, with $K< p$, are popular in a large array of applications, and mainly used for predicting a response from correlated features. In contrast, methodology and theory for inference on the regression coefficient $\beta$ relating $Y$ to $Z$ are scarce, since typically the un-observable factor $Z$ is hard to interpret. Furthermore, the determination of the asymptotic variance of an estimator of $\beta$ is a long-standing problem, with solutions known only in a few particular cases. To address some of these outstanding questions, we develop inferential tools for $\beta$ in a class of factor regression models in which the observed features are signed mixtures of the latent factors. The model specifications are practically desirable, in a large array of applications, render interpretability to the components of $Z$, and are sufficient for parameter identifiability. Without assuming that the number of latent factors $K$ or the structure of the mixture is known in advance, we construct computationally efficient estimators of $\beta$, along with estimators of other important model parameters. We benchmark the rate of convergence of $\beta$ by first establishing its $\ell_2$-norm minimax lower bound, and show that our proposed estimator is minimax-rate adaptive. Our main contribution is the provision of a unified analysis of the component-wise Gaussian asymptotic distribution of $\wh \beta$ and, especially, the derivation of a closed form expression of its asymptotic variance, together with consistent variance estimators. The resulting inferential tools can be used when both $K$ and $p$ are independent of the sample size $n$, and when both, or either, $p$ and $K$ vary with $n$, while allowing for $p > n$.


title: NMR relaxation in the spin-1 Heisenberg chain

abstract: We consider the isotropic $S=1$ Heisenberg chain with a finite Haldane gap $\Delta$ and use state-of-the-art numerical techniques to investigate its dynamical properties at finite temperature, focusing on the nuclear spin-lattice relaxation rate $1/T_1$ measured in nuclear magnetic resonance (NMR) experiments for instance. In particular, we analyze the contributions from modes with momenta close to $q\approx 0$ and $q\approx \pi$ as a function of temperature. At high-temperature, we observe spin diffusion with a non-trivial exponent. At low-temperature, we argue that a simple activated behavior $1/T_1 \propto\exp(-\Delta/T)$ can only be observed at temperatures much smaller than the gap $\Delta$.


title: Leveraging Latent Features for Local Explanations

abstract: As the application of deep neural networks proliferates in numerous areas such as medical imaging, video surveillance, and self driving cars, the need for explaining the decisions of these models has become a hot research topic, both at the global and local level. Locally, most explanation methods have focused on identifying relevance of features, limiting the types of explanations possible. In this paper, we investigate a new direction by leveraging latent features to generate contrastive explanations; predictions are explained not only by highlighting aspects that are in themselves sufficient to justify the classification, but also by new aspects which if added will change the classification. The key contribution of this paper lies in how we add features to rich data in a formal yet humanly interpretable way that leads to meaningful results. Our new definition of "addition" uses latent features to move beyond the limitations of previous explanations and resolve an open question laid out in Dhurandhar, et. al. (2018), which creates local contrastive explanations but is limited to simple datasets such as grayscale images. The strength of our approach in creating intuitive explanations that are also quantitatively superior to other methods is demonstrated on three diverse image datasets (skin lesions, faces, and fashion apparel). A user study with 200 participants further exemplifies the benefits of contrastive information, which can be viewed as complementary to other state-of-the-art interpretability methods.


title: Parton Energy Loss and the Generalized Jet Transport Coefficient

abstract: We revisit radiative parton energy loss in deeply inelastic scattering (DIS) off a large nucleus within the perturbative QCD approach. We calculate the gluon radiation spectra induced by double parton scattering in DIS without collinear expansion in the transverse momentum of initial gluons as in the original high-twist approach. The final radiative gluon spectrum can be expressed in terms of the convolution of hard partonic parts and unintegrated or transverse momentum dependent (TMD) quark-gluon correlations. The TMD quark-gluon correlation can be factorized approximately as a product of initial quark distribution and TMD gluon distribution which can be used to define the generalized or TMD jet transport coefficient. Under the static scattering center and soft radiative gluon approximation, we recover the result by Gylassy-Levai-Vitev (GLV) in the first order of the opacity expansion. The difference as a result of the soft radiative gluon approximation is investigated numerically under the static scattering center approximation.


title: Building trust for continuous variable quantum states

abstract: In this work we develop new methods for the characterisation of continuous variable quantum states using heterodyne measurement in both the trusted and untrusted settings. First, building on quantum state tomography with heterodyne detection, we introduce a reliable method for continuous variable quantum state certification, which directly yields the elements of the density matrix of the state considered with analytical confidence intervals. This method neither needs mathematical reconstruction of the data nor discrete binning of the sample space, and uses a single Gaussian measurement setting. Second, beyond quantum state tomography and without its identical copies assumption, we promote our reliable tomography method to a general efficient protocol for verifying continuous variable pure quantum states with Gaussian measurements against fully malicious adversaries, i.e., making no assumptions whatsoever on the state generated by the adversary. These results are obtained using a new analytical estimator for the expected value of any operator acting on a continuous variable quantum state with bounded support over the Fock basis, computed with samples from heterodyne detection of the state.


title: Fallout: Reading Kernel Writes From User Space

abstract: Recently, out-of-order execution, an important performance optimization in modern high-end processors, has been revealed to pose a significant security threat, allowing information leaks across security domains. In particular, the Meltdown attack leaks information from the operating system kernel to user space, completely eroding the security of the system. To address this and similar attacks, without incurring the performance costs of software countermeasures, Intel includes hardware-based defenses in its recent Coffee Lake R processors. In this work, we show that the recent hardware defenses are not sufficient. Specifically, we present Fallout, a new transient execution attack that leaks information from a previously unexplored microarchitectural component called the store buffer. We show how unprivileged user processes can exploit Fallout to reconstruct privileged information recently written by the kernel. We further show how Fallout can be used to bypass kernel address space randomization. Finally, we identify and explore microcode assists as a hitherto ignored cause of transient execution. Fallout affects all processor generations we have tested. However, we notice a worrying regression, where the newer Coffee Lake R processors are more vulnerable to Fallout than older generations.


title: Spatial Evolutionary Generative Adversarial Networks

abstract: Generative adversary networks (GANs) suffer from training pathologies such as instability and mode collapse. These pathologies mainly arise from a lack of diversity in their adversarial interactions. Evolutionary generative adversarial networks apply the principles of evolutionary computation to mitigate these problems. We hybridize two of these approaches that promote training diversity. One, E-GAN, at each batch, injects mutation diversity by training the (replicated) generator with three independent objective functions then selecting the resulting best performing generator for the next batch. The other, Lipizzaner, injects population diversity by training a two-dimensional grid of GANs with a distributed evolutionary algorithm that includes neighbor exchanges of additional training adversaries, performance based selection and population-based hyper-parameter tuning. We propose to combine mutation and population approaches to diversity improvement. We contribute a superior evolutionary GANs training method, Mustangs, that eliminates the single loss function used across Lipizzaner's grid. Instead, each training round, a loss function is selected with equal probability, from among the three E-GAN uses. Experimental analyses on standard benchmarks, MNIST and CelebA, demonstrate that Mustangs provides a statistically faster training method resulting in more accurate networks.


title: The convexity package for Hamiltonian actions on conformal symplectic manifolds

abstract: Consider a Hamiltonian action of a compact connected Lie group on a conformal symplectic manifold. We prove a convexity theorem for the moment map under the assumption that the action is of Lee type, which establishes an analog of Kirwan's convexity theorem in conformal symplectic geometry.


title: Which group algebras cannot be made zero by imposing a single non-monomial relation?

abstract: For which groups $G$ is it true that for all fields $k$, every non-monomial element of the group algebra $k\,G$ generates a proper $2$-sided ideal? The only groups for which we know this are the torsion-free abelian groups. We would like to know whether it also holds for all free groups. It is shown that the above property fails for wide classes of groups: for every group $G$ that contains an element $g\neq 1$ whose image in $G/[g,G]$ has finite order (in particular, every group containing a $g\neq 1$ that itself has finite order, or that satisfies $g\in [g,G])$; and for every group containing an element $g$ which commutes with a conjugate $hgh^{-1}\neq g$ (in particular, for every nonabelian solvable group). Results are obtained on closure properties of the class of groups satisfying the stated condition. Many further questions are raised; in particular, a plausible Freiheitssatz for group algebras of free groups is noted.


title: Robust measurement of innovation performances in Europe with a hierarchy of interacting composite indicators

abstract: For long time the measurement of innovation has been in the forefront of policy makers' and researchers' agenda worldwide. Therefore, there is an ongoing debate about which indicators should be used to measure innovation. Recent approaches have favoured the use of composite innovation indicators. However, there is no consensus about the appropriate methodology to aggregate the varying dimensions of innovation into a single summary indicator. One of the best known examples of composite innovation indicators is the European Innovation Scoreboard (EIS). It is a relevant tool for benchmarking innovation in Europe. Still, the EIS lacks a proper scheme for weighting the included indicators according to their relative importance. In this context, we propose an appraisal methodology permitting to take into consideration the interaction of criteria and robustness concerns related to the elicitation of the weights assigned to the elementary indicators. With this aim, we apply the hierarchical-SMAA-Choquet integral approach. This integrated multicriteria decision making (MCDM) method helps the users to rank and benchmark countries' innovation performance taking into account the importance and interaction of criteria assigned by themselves, rather than equal weights or weights exogenously fixed by external experts.


title: Response of active Brownian particles to boundary driving

abstract: We computationally study the behavior of underdamped active Brownian particles in a sheared channel geometry. Due to their underdamped dynamics, the particles carry momentum a characteristic distance away from the boundary before it is dissipated into the substrate. We correlate this distance with the persistence of particle trajectories, determined jointly by their friction and self-propulsion. Within this characteristic length, we observe new and counterintuitive phenomena stemming from the interplay of activity, interparticle interactions, and the boundary driving. Depending on values of friction and self-propulsion, interparticle interactions can either aid or hinder momentum transport. More dramatically, in certain cases we observe a flow reversal near the wall, which we correlate with an induced polarization of the particle self-propulsion directions. We rationalize these results in terms of a simple kinetic picture of particle trajectories.


title: Heterogeneous causal effects with imperfect compliance: a Bayesian machine learning approach

abstract: This paper introduces an innovative Bayesian machine learning algorithm to draw interpretable inference on heterogeneous causal effects in the presence of imperfect compliance (e.g., under an irregular assignment mechanism). We show, through Monte Carlo simulations, that the proposed Bayesian Causal Forest with Instrumental Variable (BCF-IV) methodology outperforms other machine learning techniques tailored for causal inference in discovering and estimating the heterogeneous causal effects while controlling for the familywise error rate (or - less stringently - for the false discovery rate) at leaves' level. BCF-IV sheds a light on the heterogeneity of causal effects in instrumental variable scenarios and, in turn, provides the policy-makers with a relevant tool for targeted policies. Its empirical application evaluates the effects of additional funding on students' performances. The results indicate that BCF-IV could be used to enhance the effectiveness of school funding on students' performance.


title: Dynamic Traffic Scene Classification with Space-Time Coherence

abstract: This paper examines the problem of dynamic traffic scene classification under space-time variations in viewpoint that arise from video captured on-board a moving vehicle. Solutions to this problem are important for realization of effective driving assistance technologies required to interpret or predict road user behavior. Currently, dynamic traffic scene classification has not been adequately addressed due to a lack of benchmark datasets that consider spatiotemporal evolution of traffic scenes resulting from a vehicle's ego-motion. This paper has three main contributions. First, an annotated dataset is released to enable dynamic scene classification that includes 80 hours of diverse high quality driving video data clips collected in the San Francisco Bay area. The dataset includes temporal annotations for road places, road types, weather, and road surface conditions. Second, we introduce novel and baseline algorithms that utilize semantic context and temporal nature of the dataset for dynamic classification of road scenes. Finally, we showcase algorithms and experimental results that highlight how extracted features from scene classification serve as strong priors and help with tactical driver behavior understanding. The results show significant improvement from previously reported driving behavior detection baselines in the literature.


title: Remarks on WDC sets

abstract: We study WDC sets, which form a substantial generalization of sets with positive reach and still admit the definition of curvature measures. Main results concern WDC sets $A\subset \mathbb{R}^2$. We prove that, for such $A$, the distance function $d_A= {\rm dist}(\cdot,A)$ is a `DC aura' for $A$, which implies that each locally WDC set in $\mathbb{R}^2$ is a WDC set. An another consequence is that compact WDC subsets of $\mathbb{R}^2$ form a Borel subset of the space of all compact sets.


title: A cooling anomaly of high-mass white dwarfs

abstract: Recently, the power of Gaia data has revealed an enhancement of high-mass white dwarfs (WDs) on the Hertzsprung--Russell diagram, called the Q branch. This branch is located at the high-mass end of the recently identified crystallization branch. Investigating its properties, we find that the number density and velocity distribution on the Q branch cannot be explained by the cooling delay of crystallization alone, suggesting the existence of an extra cooling delay. To quantify this delay, we statistically compare two age indicators -- the dynamical age inferred from transverse velocity, and the photometric isochrone age -- for more than one thousand high-mass WDs (1.08--1.23 $M_\odot$) selected from Gaia Data Release 2. We show that about 6 % of the high-mass WDs must experience an 8 Gyr extra cooling delay on the Q branch, in addition to the crystallization and merger delays. This cooling anomaly is a challenge for WD cooling models. We point out that $^{22}$Ne settling in C/O-core WDs could account for this extra cooling delay.


title: New Solutions to the Reflection Equation with Braces

abstract: The reflection equation of Cherednik is a counterpart to the celebrated Yang-Baxter equation, with importance in the theory of integrable systems. We obtain several new solutions of the reflection equation using braces building on the work of Smoktunowicz, Vendramin and Weston. In particular, we show that every brace yields several simple solutions and that a natural class of solutions is a near-ring. We find more solutions for factorizable rings. Some of our solutions apply to the original parameter-dependent equation.


abstract: Much of the recent work on learning molecular representations has been based on Graph Convolution Networks (GCN). These models rely on local aggregation operations and can therefore miss higher-order graph properties. To remedy this, we propose Path-Augmented Graph Transformer Networks (PAGTN) that are explicitly built on longer-range dependencies in graph-structured data. Specifically, we use path features in molecular graphs to create global attention layers. We compare our PAGTN model against the GCN model and show that our model consistently outperforms GCNs on molecular property prediction datasets including quantum chemistry (QM7, QM8, QM9), physical chemistry (ESOL, Lipophilictiy) and biochemistry (BACE, BBBP).


title: Geolocating Political Events in Text

abstract: This work introduces a general method for automatically finding the locations where political events in text occurred. Using a novel set of 8,000 labeled sentences, I create a method to link automatically extracted events and locations in text. The model achieves human level performance on the annotation task and outperforms previous event geolocation systems. It can be applied to most event extraction systems across geographic contexts. I formalize the event--location linking task, describe the neural network model, describe the potential uses of such a system in political science, and demonstrate a workflow to answer an open question on the role of conventional military offensives in causing civilian casualties in the Syrian civil war.


title: Nonlinear analysis of the fluid-solid transition in a model for ordered biological tissues

abstract: The rheology of biological tissues is important for their function, and we would like to better understand how single cells control global tissue properties such as tissue fluidity. A confluent tissue can fluidize when cells diffuse by executing a series of cell rearrangements, or T1 transitions. In a disordered 2D vertex model, the tissue fluidizes when the T1 energy barriers disappear as the target shape index approaches a critical value ($s^*_{0} \sim 3.81$), and the shear modulus describing the linear response also vanishes at this same critical point. However, the ordered ground states of 2D vertex models become linearly unstable at a lower value of the target shape index (3.72) [1,2]. We investigate whether the ground states of the 2D vertex model are fluid-like or solid-like between 3.72 and 3.81 $-$ does the "equation of state" for these systems have two branches, like glassy particulate matter, or only one? Using four-cell and many-cell numerical simulations, we demonstrate that for a hexagonal ground state, T1 energy barriers only vanish at $\sim 3.81$, indicating that ordered systems have the same critical point as disordered systems. We also develop a simple geometric argument that correctly predicts how non-linear stabilization disappears at $s^*_{0}$ in ordered systems.


title: A Topological Characterization of the Middle Perversity Intersection Complex for Arbitrary Complex Algebraic Varieties

abstract: For an arbitrary complex algebraic variety which is not necessarily pure dimensional, the intersection complex can be defined as the direct sum of the Deligne-Goresky-MacPherson intersection complexes of each irreducible component. We give two axiomatic topological characterizations of the middle perversity direct sum intersection complex, one stratification dependent and the other stratification independent. To accomplish this, we show that this direct sum intersection complex can be constructed using Deligne's construction in the more general context of topologically stratified spaces. A consequence of these characterizations is the invariance of this direct sum intersection complex under homeomorphisms.


title: The Fundamental Solution to One-Dimensional Degenerate Diffusion Equation, I

abstract: In this work we adopt a combination of probabilistic approach and analytic methods to study the fundamental solutions to variations of the Wright-Fisher equation in one dimension. To be specific, we consider a diffusion equation on $\left(0,\infty\right)$ whose diffusion coefficient vanishes at the boundary 0, equipped with the Cauchy initial data and the Dirichlet boundary condition. One type of diffusion operator that has been extensively studied is the one whose diffusion coefficient vanishes linearly at 0. Our main goal is to extend the study to cases when the diffusion coefficient has a general order of degeneracy. We primarily focus on the fundamental solution to such a degenerate diffusion equation. In particular, we study the regularity properties of the fundamental solution near 0, and investigate how the order of degeneracy of the diffusion operator and the Dirichlet boundary condition jointly affect these properties. We also provide estimates for the fundamental solution and its derivatives near 0.


title: An adaptive nearest neighbor rule for classification

abstract: We introduce a variant of the $k$-nearest neighbor classifier in which $k$ is chosen adaptively for each query, rather than supplied as a parameter. The choice of $k$ depends on properties of each neighborhood, and therefore may significantly vary between different points. (For example, the algorithm will use larger $k$ for predicting the labels of points in noisy regions.) We provide theory and experiments that demonstrate that the algorithm performs comparably to, and sometimes better than, $k$-NN with an optimal choice of $k$. In particular, we derive bounds on the convergence rates of our classifier that depend on a local quantity we call the `advantage' which is significantly weaker than the Lipschitz conditions used in previous convergence rate proofs. These generalization bounds hinge on a variant of the seminal Uniform Convergence Theorem due to Vapnik and Chervonenkis; this variant concerns conditional probabilities and may be of independent interest.


title: From Halfspace M-depth to Multiple-output Expectile Regression

abstract: Despite the renewed interest in the Newey and Powell (1987) concept of expectiles in fields such as econometrics, risk management, and extreme value theory, expectile regression---or, more generally, M-quantile regression---unfortunately remains limited to single-output problems. To improve on this, we introduce hyperplane-valued multivariate M-quantiles that show strong advantages, for instance in terms of equivariance, over the various point-valued multivariate M-quantiles available in the literature. Like their competitors, our multivariate M-quantiles are directional in nature and provide centrality regions when all directions are considered. These regions define a new statistical depth, the halfspace M-depth, whose deepest point, in the expectile case, is the mean vector. Remarkably, the halfspace M-depth can alternatively be obtained by substituting, in the celebrated Tukey (1975) halfspace depth, M-quantile outlyingness for standard quantile outlyingness, which supports a posteriori the claim that our multivariate M-quantile concept is the natural one. We investigate thoroughly the properties of the proposed multivariate M-quantiles, of halfspace M-depth, and of the corresponding regions. Since our original motivation was to define multiple-output expectile regression methods, we further focus on the expectile case. We show in particular that expectile depth is smoother than the Tukey depth and enjoys interesting monotonicity properties that are extremely promising for computational purposes. Unlike their quantile analogs, the proposed multivariate expectiles also satisfy the coherency axioms of multivariate risk measures. Finally, we show that our multivariate expectiles indeed allow performing multiple-output expectile regression, which is illustrated on simulated and real data.


title: Physical mechanisms involved in the formation and operation of memory devices based on a monolayer of gold nanoparticles-polythiophene hybrid materials

abstract: Understanding the physical and chemical mechanisms occurring during the forming process and operation of an organic resistive memory device is a major issue for better performances. Various mechanisms were suggested in vertically stacked memory structures, but the analysis remains indirect and needs destructive characterization (e.g. cross-section to access the organic layers sandwiched between electrodes). Here, we report a study on a planar, monolayer thick, hybrid nanoparticle/molecule device (10 nm gold nanoparticles embedded in an electro-generated poly(2-thienyl-3,4-(ethylenedioxy)thiophene) layer), combining, in situ, on the same device, physical (scanning electron microscope, physico-chemical (thermogravimetry and mass spectroscopy, Raman spectroscopy) and electrical (temperature dependent current-voltage) characterizations. We demonstrate that the forming process causes an increase in the gold particle size, almost 4 times larger than the starting nanoparticles, and that the organic layer undergoes a significant chemical rearrangement from a sp3 to sp2 amorphous carbon material. Temperature dependent electrical characterizations of this nonvolatile memory confirm that the charge transport mechanism in the device is consistent with a trap-filled space charge limited current in the off state, the sp2 amorphous carbon material containing many electrically active defects.


title: Evaluation of pilot jobs for Apache Spark applications on HPC clusters

abstract: Big Data has become prominent throughout many scientific fields and, as a result, scientific communities have sought out Big Data frameworks to accelerate the processing of their increasingly data-intensive pipelines. However, while scientific communities typically rely on High-Performance Computing (HPC) clusters for the parallelization of their pipelines, many popular Big Data frameworks such as Hadoop and Apache Spark were primarily designed to be executed on dedicated commodity infrastructures. This paper evaluates the benefits of pilot jobs over traditional batch submission for Apache Spark on HPC clusters. Surprisingly, our results show that the speed-up provided by pilot jobs over batch scheduling is moderate to inexistent (0.98 on average) despite the presence of long queuing times. In addition, pilot jobs provide an extra layer of scheduling that complexifies debugging and deployment. We conclude that traditional batch scheduling should remain the default strategy to deploy Apache Spark applications on HPC clusters.


title: Matrix-Free Preconditioning in Online Learning

abstract: We provide an online convex optimization algorithm with regret that interpolates between the regret of an algorithm using an optimal preconditioning matrix and one using a diagonal preconditioning matrix. Our regret bound is never worse than that obtained by diagonal preconditioning, and in certain setting even surpasses that of algorithms with full-matrix preconditioning. Importantly, our algorithm runs in the same time and space complexity as online gradient descent. Along the way we incorporate new techniques that mildly streamline and improve logarithmic factors in prior regret analyses. We conclude by benchmarking our algorithm on synthetic data and deep learning tasks.


title: Higgs Parity Grand Unification

abstract: The vanishing of the Higgs quartic coupling of the Standard Model at high energies may be explained by spontaneous breaking of Higgs Parity. Taking Higgs Parity to originate from the Left-Right symmetry of the $SO(10)$ gauge group, leads to a new scheme for precision gauge coupling unification that is consistent with proton decay. We compute the relevant running of couplings and threshold corrections to allow a precise correlation among Standard Model parameters. The scheme has a built-in solution for obtaining a realistic value for $m_b/m_\tau$, which further improves the precision from gauge coupling unification, allowing the QCD coupling constant to be predicted to the level of 1 % or, alternatively, the top quark mass to 0.2 %. Future measurements of these parameters may significantly constrain the detailed structure of the theory. We also study an $SO(10)$ embedding of quark and lepton masses, showing how large neutrino mixing is compatible with small quark mixing, and predict a normal neutrino mass hierarchy. The strong CP problem may be explained by combining Higgs Parity with space-time parity.


title: Extending Monocular Visual Odometry to Stereo Camera Systems by Scale Optimization

abstract: This paper proposes a novel approach for extending monocular visual odometry to a stereo camera system. The proposed method uses an additional camera to accurately estimate and optimize the scale of the monocular visual odometry, rather than triangulating 3D points from stereo matching. Specifically, the 3D points generated by the monocular visual odometry are projected onto the other camera of the stereo pair, and the scale is recovered and optimized by directly minimizing the photometric error. It is computationally efficient, adding minimal overhead to the stereo vision system compared to straightforward stereo matching, and is robust to repetitive texture. Additionally, direct scale optimization enables stereo visual odometry to be purely based on the direct method. Extensive evaluation on public datasets (e.g., KITTI), and outdoor environments (both terrestrial and underwater) demonstrates the accuracy and efficiency of a stereo visual odometry approach extended by scale optimization, and its robustness in environments with challenging textures.


title: Variational Diffusion Autoencoders with Random Walk Sampling

abstract: Variational autoencoders (VAEs) and generative adversarial networks (GANs) enjoy an intuitive connection to manifold learning: in training the decoder/generator is optimized to approximate a homeomorphism between the data distribution and the sampling space. This is a construction that strives to define the data manifold. A major obstacle to VAEs and GANs, however, is choosing a suitable prior that matches the data topology. Well-known consequences of poorly picked priors are posterior and mode collapse. To our knowledge, no existing method sidesteps this user choice. Conversely, $\textit{diffusion maps}$ automatically infer the data topology and enjoy a rigorous connection to manifold learning, but do not scale easily or provide the inverse homeomorphism (i.e. decoder/generator). We propose a method that combines these approaches into a generative model that inherits the asymptotic guarantees of $\textit{diffusion maps}$ while preserving the scalability of deep models. We prove approximation theoretic results for the dimension dependence of our proposed method. Finally, we demonstrate the effectiveness of our method with various real and synthetic datasets.


title: DBFFT: A displacement based FFT approach for non-linear homogenization of the mechanical behavior

abstract: Most of the FFT methods available for homogenization of the mechanical response use the strain/deformation gradient as unknown, imposing their compatibility using Green's functions or projection operators. This implies the allocation of redundant information and, when the method is based in solving a linear equation, the rank-deficiency of the resulting system. In this work we propose a fast, robust and memory-efficient FFT homogenization framework in which the displacement field on the Fourier space is the unknown: the displacement based FFT (DBFFT) algorithm. The framework allows any general non-linear constitutive behavior for the phases and direct strain, stress and mixed control of the macroscopic load. In the linear case, the method results in a linear system defined in terms of linear operators in the Fourier space and that does not require a reference medium. The system has an associated full rank Hermitian matrix and can be solved using iterative Krylov solvers and allows the use of preconditioners. A preconditioner is proposed to improve the efficiency of the system resolution. Finally, some numerical examples including elastic, hyperelastic and viscoplastic materials are solved to check the accuracy and efficiency of the method. The computational cost reduction respect the Galerkin-FFT was around 30%.


title: Prioritized Sequence Experience Replay

abstract: Experience replay is widely used in deep reinforcement learning algorithms and allows agents to remember and learn from experiences from the past. In an effort to learn more efficiently, researchers proposed prioritized experience replay (PER) which samples important transitions more frequently. In this paper, we propose Prioritized Sequence Experience Replay (PSER) a framework for prioritizing sequences of experience in an attempt to both learn more efficiently and to obtain better performance. We compare the performance of PER and PSER sampling techniques in a tabular Q-learning environment and in DQN on the Atari 2600 benchmark. We prove theoretically that PSER is guaranteed to converge faster than PER and empirically show PSER substantially improves upon PER.


title: The influence of gravity on the Boltzmann entropy of a closed universe

abstract: This contribution inquires into Clausius' proposal that "the entropy of the world tends to a maximum.'" The question is raised whether the entropy of "the world" actually does have a maximum; and if the answer is "Yes!," what such states of maximum entropy look like, and if the answer is "No!," what this could entail for the fate of the universe. Following R. Penrose, "the world" is modelled by a closed Friedman--Lemaitre type universe, in which a three-dimensional spherical "space" is filled with "matter" consisting of $N$ point particles, their large-scale distribution being influenced by their own gravity. As "entropy of matter" the Boltzmann entropy for a (semi-)classical macrostate, and Boltzmann's ergodic ensemble formulation of it for an isolated thermal equilibrium state, are studied. Since the notion of a Boltzmann entropy is not restricted to classical non-relativistic physics, the inquiry will take into account quantum theory as well as relativity theory; we also consider black hole entropy. Model universes having a maximum entropy state and those which don't will be encountered. It is found that the answer to our maximum entropy question is not at all straightforward at the general-relativistic level. In particular, it is shown that the increase in Bekenstein--Hawking entropy of general-relativistic black holes does not always compensate for the Boltzmann entropy of a piece of matter swallowed by a black hole.

categories: physics.hist-ph gr-qc math-ph math.MP physics.class-ph


title: Fairness and Missing Values

abstract: The causes underlying unfair decision making are complex, being internalised in different ways by decision makers, other actors dealing with data and models, and ultimately by the individuals being affected by these decisions. One frequent manifestation of all these latent causes arises in the form of missing values: protected groups are more reluctant to give information that could be used against them, delicate information for some groups can be erased by human operators, or data acquisition may simply be less complete and systematic for minority groups. As a result, missing values and bias in data are two phenomena that are tightly coupled. However, most recent techniques, libraries and experimental results dealing with fairness in machine learning have simply ignored missing data. In this paper, we claim that fairness research should not miss the opportunity to deal properly with missing data. To support this claim, (1) we analyse the sources of missing data and bias, and we map the common causes, (2) we find that rows containing missing values are usually fairer than the rest, which should not be treated as the uncomfortable ugly data that different techniques and libraries get rid of at the first occasion, and (3) we study the trade-off between performance and fairness when the rows with missing values are used (either because the technique deals with them directly or by imputation methods). We end the paper with a series of recommended procedures about what to do with missing data when aiming for fair decision making.


title: Zeroth-Order Stochastic Alternating Direction Method of Multipliers for Nonconvex Nonsmooth Optimization

abstract: Alternating direction method of multipliers (ADMM) is a popular optimization tool for the composite and constrained problems in machine learning. However, in many machine learning problems such as black-box attacks and bandit feedback, ADMM could fail because the explicit gradients of these problems are difficult or infeasible to obtain. Zeroth-order (gradient-free) methods can effectively solve these problems due to that the objective function values are only required in the optimization. Recently, though there exist a few zeroth-order ADMM methods, they build on the convexity of objective function. Clearly, these existing zeroth-order methods are limited in many applications. In the paper, thus, we propose a class of fast zeroth-order stochastic ADMM methods (i.e., ZO-SVRG-ADMM and ZO-SAGA-ADMM) for solving nonconvex problems with multiple nonsmooth penalties, based on the coordinate smoothing gradient estimator. Moreover, we prove that both the ZO-SVRG-ADMM and ZO-SAGA-ADMM have convergence rate of $O(1/T)$, where $T$ denotes the number of iterations. In particular, our methods not only reach the best convergence rate $O(1/T)$ for the nonconvex optimization, but also are able to effectively solve many complex machine learning problems with multiple regularized penalties and constraints. Finally, we conduct the experiments of black-box binary classification and structured adversarial attack on black-box deep neural network to validate the efficiency of our algorithms.


title: Recursive Sketches for Modular Deep Learning

abstract: We present a mechanism to compute a sketch (succinct summary) of how a complex modular deep network processes its inputs. The sketch summarizes essential information about the inputs and outputs of the network and can be used to quickly identify key components and summary statistics of the inputs. Furthermore, the sketch is recursive and can be unrolled to identify sub-components of these components and so forth, capturing a potentially complicated DAG structure. These sketches erase gracefully; even if we erase a fraction of the sketch at random, the remainder still retains the `high-weight' information present in the original sketch. The sketches can also be organized in a repository to implicitly form a `knowledge graph'; it is possible to quickly retrieve sketches in the repository that are related to a sketch of interest; arranged in this fashion, the sketches can also be used to learn emerging concepts by looking for new clusters in sketch space. Finally, in the scenario where we want to learn a ground truth deep network, we show that augmenting input/output pairs with these sketches can theoretically make it easier to do so.


title: Protecting quantum entanglement from leakage and qubit errors via repetitive parity measurements

abstract: Protecting quantum information from errors is essential for large-scale quantum computation. Quantum error correction (QEC) encodes information in entangled states of many qubits, and performs parity measurements to identify errors without destroying the encoded information. However, traditional QEC cannot handle leakage from the qubit computational space. Leakage affects leading experimental platforms, based on trapped ions and superconducting circuits, which use effective qubits within many-level physical systems. We investigate how two-transmon entangled states evolve under repeated parity measurements, and demonstrate the use of hidden Markov models to detect leakage using only the record of parity measurement outcomes required for QEC. We show the stabilization of Bell states over up to 26 parity measurements by mitigating leakage using postselection, and correcting qubit errors using Pauli-frame transformations. Our leakage identification method is computationally efficient and thus compatible with real-time leakage tracking and correction in larger quantum processors.


title: On a class of generalized solutions to equations describing incompressible viscous fluids

abstract: We consider a class of viscous fluids with a general monotone dependence of the viscous stress on the symmetric velocity gradient. We introduce the concept of dissipative solution to the associated initial boundary value problem inspired by the measure-valued solutions for the inviscid (Euler) system. We show the existence as well as the weak-strong uniqueness property in the class of dissipative solutions. Finally, the dissipative solution enjoying certain extra regularity coincides with a strong solution of the same problem.


title: Accelerating Min-Max Optimization with Application to Minimal Bounding Sphere

abstract: We study the min-max optimization problem where each function contributing to the max operation is strongly-convex and smooth with bounded gradient in the search domain. By smoothing the max operator, we show the ability to achieve an arbitrarily small positive optimality gap of $\delta$ in $\tilde{O}(1/\sqrt{\delta})$ computational complexity (up to logarithmic factors) as opposed to the state-of-the-art strong-convexity computational requirement of $O(1/\delta)$. We apply this important result to the well-known minimal bounding sphere problem and demonstrate that we can achieve a $(1+\varepsilon)$-approximation of the minimal bounding sphere, i.e. identify an hypersphere enclosing a total of $n$ given points in the $d$ dimensional unbounded space $\mathbb{R}^d$ with a radius at most $(1+\varepsilon)$ times the actual minimal bounding sphere radius for an arbitrarily small positive $\varepsilon$, in $\tilde{O}(n d /\sqrt{\varepsilon})$ computational time as opposed to the state-of-the-art approach of core-set methodology, which needs $O(n d /\varepsilon)$ computational time.


title: Sub-Turing Islands in the Wild

abstract: Recently, there has been growing debate as to whether or not static analysis can be truly sound. In spite of this concern, research on techniques seeking to at least partially answer undecidable questions has a long history. However, little attention has been given to the more empirical question of how often an exact solution might be given to a question despite the question being, at least in theory, undecidable. This paper investigates this issue by exploring sub-Turing islands -- regions of code for which a question of interest is decidable. We define such islands and then consider how to identify them. We implemented Cook, a prototype for finding sub-Turing islands and applied it to a corpus of 1100 Android applications, containing over 2 million methods. Results reveal that 55\% of the all methods are sub-Turing. Our results also provide empirical, scientific evidence for the scalability of sub-Turing island identification. Sub-Turing identification has many downstream applications, because islands are so amenable to static analysis. We illustrate two downstream uses of the analysis. In the first, we found that over 37\% of the verification conditions associated with runtime exceptions fell within sub-Turing islands and thus are statically decidable. A second use of our analysis is during code review where it provides guidance to developers. The sub-Turing islands from our study turns out to contain significantly fewer bugs than `theswamp' (non sub-Turing methods). The greater bug density in the swamp is unsurprising; the fact that bugs remain prevalent in islands is, however, surprising: these are bugs whose repair can be fully automated.


title: Single-Boundary Control of the Two-Phase Stefan System

abstract: This paper presents the control design of the two-phase Stefan problem. The two-phase Stefan problem is a representative model of liquid-solid phase transition by describing the time evolutions of the temperature profile which is divided by subdomains of liquid and solid phases as the liquid-solid moving interface position. The mathematical formulation is given by two diffusion partial differential equations (PDEs) defined on a time-varying spatial domain described by an ordinary differential equation (ODE) driven by the Neumann boundary values of both PDE states, resulting in a nonlinear coupled PDE-ODE-PDE system. We design a state feedback control law by means of energy-shaping to stabilize the interface position to a desired setpoint by using single boundary heat input. We prove that the closed-loop system under the control law ensures some conditions for model validity and the global exponential stability estimate is shown in $L_2$ norm. Furthermore, the robustness of the closed-loop stability with respect to the uncertainties of the physical parameters is shown. Numerical simulation is provided to illustrate the good performance of the proposed control law in comparison to the control design for the one-phase Stefan problem.


title: Antiparticles as particles: QED$_{2+1}$ and composite Fermions at $\nu = \frac{1}{2}$

abstract: We present a continuum operator map that inverts the role of particles and antiparticles in three dimensional QED. This is accomplished by the attachment of specific holomorphic Wilson lines to Dirac fermions. We show that this nonlocal map provides a continuum realization of Son's composite fermions at $\nu = \frac{1}{2}$. The inversion of Landau levels and the Dirac-cone--composite-Fermion duality is explicitly demonstrated for the case of slowly varying magnetic fields. The role of Maxwell-terms as well as the connection of this construction to a gauge-invariant formulation of 3D gauge theories is also elaborated upon.


title: Training Data Subset Search with Ensemble Active Learning

abstract: Deep Neural Networks (DNNs) often rely on very large datasets for training. Given the large size of such datasets, it is conceivable that they contain certain samples that either do not contribute or negatively impact the DNN's optimization. Modifying the training distribution in a way that excludes such samples could provide an effective solution to both improve performance and reduce training time. In this paper, we propose to scale up ensemble Active Learning (AL) methods to perform acquisition at a large scale (10k to 500k samples at a time). We do this with ensembles of hundreds of models, obtained at a minimal computational cost by reusing intermediate training checkpoints. This allows us to automatically and efficiently perform a training data subset search for large labeled datasets. We observe that our approach obtains favorable subsets of training data, which can be used to train more accurate DNNs than training with the entire dataset. We perform an extensive experimental study of this phenomenon on three image classification benchmarks (CIFAR-10, CIFAR-100 and ImageNet), as well as an internal object detection benchmark for prototyping perception models for autonomous driving. Unlike existing studies, our experiments on object detection are at the scale required for production-ready autonomous driving systems. We provide insights on the impact of different initialization schemes, acquisition functions and ensemble configurations at this scale. Our results provide strong empirical evidence that optimizing the training data distribution can provide significant benefits on large scale vision tasks.


title: Imaging the Molecular Interstellar Medium in a Gravitationally Lensed Star-forming Galaxy at z=5.7

abstract: Aims: We present and study spatially resolved imaging obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of multiple $^{12}$CO($J=$6$-$5, 8$-$7 and 9$-$8) and two H$_2$O(2$_{02}-$1$_{11}$ and 2$_{11}-$2$_{02}$) emission lines and cold dust continuum toward the gravitationally lensed dusty star forming galaxy SPT0346-52 at z=$5.656$. Methods: Using a visibility-domain source-plane reconstruction we probe the structure and dynamics of the different components of the interstellar medium (ISM) in this galaxy down to scales of 1 kpc in the source plane. Results: Measurements of the intrinsic sizes of the different CO emission lines indicate that the higher J transitions trace more compact regions in the galaxy. Similarly, we find smaller dust continuum intrinsic sizes with decreasing wavelength, based on observations at rest-frame 130, 300 and 450$\mu$m. The source shows significant velocity structure, and clear asymmetry where an elongated structure is observed in the source plane with significant variations in their reconstructed sizes. This could be attributed to a compact merger or turbulent disk rotation. The differences in velocity structure through the different line tracers, however, hint at the former scenario in agreement with previous [CII] line imaging results. Measurements of the CO line ratios and magnifications yield significant variations as a function of velocity, suggesting that modeling of the ISM using integrated values could be misinterpreted. Modeling of the ISM in SPT0346-52 based on delensed fluxes indicate a highly dense and warm medium, qualitatively similar to that observed in high redshift quasar hosts.


title: Nucleation and growth bottleneck in the conductivity recovery dynamics of nickelate thin films

abstract: We investigate THz conductivity dynamics in NdNiO$_3$ and EuNiO$_3$ thin films following a photoinduced thermal quench into the metallic state and reveal a clear contrast between first- and second-order dynamics. While in EuNiO$_3$ the conductivity recovers exponentially, in NdNiO$_3$ the recovery is non-exponential and slower than a simple thermal model. Crucially, it is consistent with first-order dynamics and well-described by a 2d Avrami model, with supercooling leading to metastable phase coexistence. The large transients seen in our films are promising for fast electronic (and magnetic) switching applications.


title: Melting of charge order in the low-temperature state of an electronic ferroelectric

abstract: Strong electronic interactions can drive a system into a state with a symmetry breaking. Lattice frustration or competing interactions tend to prevent a symmetry breaking, leading to quantum disordered phases. In spin systems frustration can produce a spin liquid state. Frustration of a charge degree of freedom also can result in various exotic states, however, experimental data on these effects is scarce. In this work we demonstrate how a charge ordered ferroelectric looses the order on cooling to low temperatures using an example of a Mott insulator on a weakly anisotropic triangular lattice $\kappa$-(BEDT-TTF)$_2$Hg(SCN)$_2$Cl. Typically, a low temperature ordered state is a ground state of a system, and the demonstrated re-entrant behavior is unique. Raman scattering spectroscopy finds that this material enters an insulating ferroelectric `dipole solid' state at $T=30~K$, but below $T=15~K$ the order melts, while preserving the insulating energy gap. The resulting phase diagram is relevant to other quantum paraelectric materials.


title: Adapting Text Embeddings for Causal Inference

abstract: Does adding a theorem to a paper affect its chance of acceptance? Does labeling a post with the author's gender affect the post popularity? This paper develops a method to estimate such causal effects from observational text data, adjusting for confounding features of the text such as the subject or writing quality. We assume that the text suffices for causal adjustment but that, in practice, it is prohibitively high-dimensional. To address this challenge, we develop causally sufficient embeddings, low-dimensional document representations that preserve sufficient information for causal identification and allow for efficient estimation of causal effects. Causally sufficient embeddings combine two ideas. The first is supervised dimensionality reduction: causal adjustment requires only the aspects of text that are predictive of both the treatment and outcome. The second is efficient language modeling: representations of text are designed to dispose of linguistically irrelevant information, and this information is also causally irrelevant. Our method adapts language models (specifically, word embeddings and topic models) to learn document embeddings that are able to predict both treatment and outcome. We study causally sufficient embeddings with semi-synthetic datasets and find that they improve causal estimation over related embedding methods. We illustrate the methods by answering the two motivating questions---the effect of a theorem on paper acceptance and the effect of a gender label on post popularity. Code and data available at https://github.com/vveitch/causal-text-embeddings-tf2}{github.com/vveitch/causal-text-embeddings-tf2


title: Multinuclear Absolute MR Thermometry

abstract: Non-invasive measurement of absolute temperature is important for proper characterization of various pathologies and for evaluation of thermal dose during interventional procedures. The proton magnetic resonance (MR) frequency shift method can be used to map relative temperature changes in tissues; however, spatiotemporal variations in the main magnetic field and the lack of an internal frequency reference at each location challenge the determination of absolute temperature. Here, we introduce a novel multinuclear method for absolute MR thermometry, based on the fact that the proton and sodium nuclei exhibit a unique and distinct characteristic frequency dependence with temperature. A one-to-one mapping between the precession frequency difference of the two nuclei and absolute temperature is demonstrated. Proof-of-concept experiments were conducted in aqueous solutions with different NaCl concentrations, in agarose gel samples, and in freshly excised ex vivo mouse tissues. One-dimensional chemical shift imaging was also conducted at two steady-state temperatures, demonstrating excellent agreement with infrared measurements.


title: Disorder and interaction in chiral chains: Majoranas vs complex fermions

abstract: We study the low-energy physics of a chain of Majorana fermions in the presence of interaction and disorder, emphasizing the difference between Majoranas and conventional (complex) fermions. While in the non-interacting limit both models are equivalent (in particular, belong to the same symmetry class BDI and flow towards the same infinite-randomness critical fixed point), their behavior differs drastically once interaction is added. Our density-matrix renormalization group calculations show that the complex-fermion chain remains at the non-interacting fixed point. On the other hand, the Majorana fermion chain experiences a spontaneous symmetry breaking and localizes for repulsive interaction. To explain the instability of the critical Majorana chain with respect to a combined effect of interaction and disorder, we consider interaction as perturbation to the infinite-randomness fixed point and calculate numerically two-wavefunction correlation functions that enter interaction matrix elements. The numerical results supported by analytical arguments exhibit a rich structure of critical eigenstate correlations. This allows us to identify a relevant interaction operator that drives the Majorana chain away from the infinite randomness fixed point. For the case of complex fermions, the interaction is irrelevant.


title: Fair Decision Making using Privacy-Protected Data

abstract: Data collected about individuals is regularly used to make decisions that impact those same individuals. We consider settings where sensitive personal data is used to decide who will receive resources or benefits. While it is well known that there is a tradeoff between protecting privacy and the accuracy of decisions, we initiate a first-of-its-kind study into the impact of formally private mechanisms (based on differential privacy) on fair and equitable decision-making. We empirically investigate novel tradeoffs on two real-world decisions made using U.S. Census data (allocation of federal funds and assignment of voting rights benefits) as well as a classic apportionment problem. Our results show that if decisions are made using an $\epsilon$-differentially private version of the data, under strict privacy constraints (smaller $\epsilon$), the noise added to achieve privacy may disproportionately impact some groups over others. We propose novel measures of fairness in the context of randomized differentially private algorithms and identify a range of causes of outcome disparities.