title: Software Engineering for Fairness: A Case Study with Hyperparameter Optimization

abstract: We assert that it is the ethical duty of software engineers to strive to reduce software discrimination. This paper discusses how that might be done. This is an important topic since machine learning software is increasingly being used to make decisions that affect people's lives. Potentially, the application of that software will result in fairer decisions because (unlike humans) machine learning software is not biased. However, recent results show that the software within many data mining packages exhibits "group discrimination"; i.e. their decisions are inappropriately affected by "protected attributes"(e.g., race, gender, age, etc.). There has been much prior work on validating the fairness of machine-learning models (by recognizing when such software discrimination exists). But after detection, comes mitigation. What steps can ethical software engineers take to reduce discrimination in the software they produce? This paper shows that making \textit{fairness} as a goal during hyperparameter optimization can (a) preserve the predictive power of a model learned from a data miner while also (b) generates fairer results. To the best of our knowledge, this is the first application of hyperparameter optimization as a tool for software engineers to generate fairer software.


title: Combining Parametric and Nonparametric Models for Off-Policy Evaluation

abstract: We consider a model-based approach to perform batch off-policy evaluation in reinforcement learning. Our method takes a mixture-of-experts approach to combine parametric and non-parametric models of the environment such that the final value estimate has the least expected error. We do so by first estimating the local accuracy of each model and then using a planner to select which model to use at every time step as to minimize the return error estimate along entire trajectories. Across a variety of domains, our mixture-based approach outperforms the individual models alone as well as state-of-the-art importance sampling-based estimators.


title: Stable components and layers

abstract: Component graphs $\Gamma_{0}(F)$ are defined for arrays of sets $F$, and in particular for arrays of path components for Vietoris-Rips complexes and Lesnick complexes. The path components of $\Gamma_{0}(F)$ are the {\it stable components} of the array $F$. The stable components for the system of Lesnick complexes $\{ L_{s,k}(X) \}$ for a finite data set $X$ decompose into layers, which are themselves path components of a graph. Combinatorial scoring functions are defined for layers and stable components.


title: Breakdown of electroneutrality in nanopores

abstract: Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to its unique physical properties, and the technological advances that allow us to study them. It is tempting to approach this regime with the tools and knowledge developed for wider microfluidic devices and use continuum models like the Poisson-Nernst-Planck equation. However, it turns out that some of the most basic principles we take for granted in a large system, such as electroneutrality, can breakdown under extreme confinement. We show that in a truly one-dimensional system, interacting with three-dimensional electrostatic interactions, the screening length is exponentially large in ionic spacing, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality is broken, and ionic transport can behave in a completely different way. In this work, we build a theoretical framework for electroneutrality breakdown in a one-dimensional nanopore and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.


title: A UV perspective on mixed anomalies at critical points between bosonic symmetry-protected phases

abstract: Symmetry-protected phases are gapped phases of matter which are distinguished only in the presence of a global symmetry G. These quantum phases lack any symmetry-breaking or topological order, and have short-range entangled ground states. Based on this short-range entanglement property, we give a general argument for the existence of an emergent unitary or anti-unitary $\mathbb{Z}_2$ symmetry at a critical point separating two different bosonic symmetry-protected phases in any dimension. Often, the emergent symmetry group at criticality has a mixed global anomaly. For those phases classified by group cohomology, we identify a criterion for when such a mixed global anomaly is present, and write down representative cocycles for the corresponding anomaly class. We illustrate our results with a series of examples, and make connections to recent results on (2 + 1)-d beyond-Landau critical points.


title: Thermodynamics of symmetric spin--orbital model: One- and two-dimensional cases

abstract: The specific heat and susceptibilities for the two- and one-dimensional spin--orbital models are calculated in the framework of a spherically symmetric self-consistent approach at different temperatures and relations between the parameters of the system. It is shown that even in the absence of the long-range spin and orbital order, the system exhibits the features in the behavior of thermodynamic characteristics, which are typical of those manifesting themselves at phase transitions. Such features are attributed to the quantum entanglement of the coupled spin and orbital degrees of freedom.


title: Simulation Typology and Termination Risks

abstract: The goal of the article is to explore what is the most probable type of simulation in which humanity lives (if any) and how this affects simulation termination risks. We firstly explore the question of what kind of simulation in which humanity is most likely located based on pure theoretical reasoning. We suggest a new patch to the classical simulation argument, showing that we are likely simulated not by our own descendants, but by alien civilizations. Based on this, we provide classification of different possible simulations and we find that simpler, less expensive and one-person-centered simulations, resurrectional simulations, or simulations of the first artificial general intelligence's (AGI's) origin (singularity simulations) should dominate. Also, simulations which simulate the 21st century and global catastrophic risks are probable. We then explore whether the simulation could collapse or be terminated. Most simulations must be terminated after they model the singularity or after they model a global catastrophe before the singularity. Undeniably observed glitches, but not philosophical speculations could result in simulation termination. The simulation could collapse if it is overwhelmed by glitches. The Doomsday Argument in simulations implies termination soon. We conclude that all types of the most probable simulations except resurrectional simulations are prone to termination risks in a relatively short time frame of hundreds of years or less from now.


title: Mather Measures and Ergodic Properties of Kantorovich Operators associated to General Mass Transfers

abstract: We introduce and study the class of linear transfers between probability distributions and the dual class of Kantorovich operators between function spaces. Linear transfers can be seen as an extension of convex lower semi-continuous energies on Wasserstein space, of cost minimizing mass transports, as well as many other couplings between probability measures to which Monge-Kantorovich theory does not readily apply. Basic examples include balayage of measures, martingale transports, optimal Skorokhod embeddings, and the weak mass transports of Talagrand, Marton, Gozlan and others. The class also includes various stochastic mass transports such as the Schr\"odinger bridge associated to a reversible Markov process, and the Arnold-Brenier variational principle for the incompressible Euler equations. We associate to most linear transfers, a critical constant, a corresponding effective linear transfer and additive eigenfunctions to their dual Kantorovich operators, that extend Man\'e's critical value, Aubry-Mather invariant tori, and Fathi's weak KAM solutions for Hamiltonian systems. This amounts to studying the asymptotic properties of the nonlinear Kantorovich operators as opposed to classical ergodic theory, which deals with linear Markov operators. This allows for the extension of Mather theory to other settings such as its stochastic counterpart. We also introduce the class of convex transfers, which includes $p$-powers ($ p \geq 1$) of linear transfers, the logarithmic entropy, the Donsker-Varadhan information, optimal mean field plans, and certain free energies as functions of two probability measures, i.e., where the reference measure is also a variable. Duality formulae for general transfer inequalities follow in a very natural way. This paper is an expanded version of a previously posted but not published work by the authors.


title: Evidence of Systematic Errors in $Spitzer$ Microlens Parallax Measurements

abstract: The microlensing parallax campaign with the $Spitzer$ space telescope aims to measure masses and distances of microlensing events seen towards the Galactic bulge, with a focus on planetary microlensing events. The hope is to measure how the distribution of planets depends on position within the Galaxy. In this paper, we compare 50 microlens parallax measurements from the 2015 $Spitzer$ campaign to three different Galactic models commonly used in microlensing analyses, and we find that $\geq 74\%$ of these events have microlensing parallax values higher than the medians predicted by Galactic models. The Anderson-Darling tests indicate probabilities of $p_{\rm AD} < 6.6 \times 10^{-5}$ for these three Galactic models, while the binomial probability of such a large fraction of large microlensing parallax values is $< 4.6\times 10^{-4}$. Given that many $Spitzer$ light curves show evidence of large correlated errors, we conclude that this discrepancy is probably due to systematic errors in the $Spitzer$ photometry. We find formally acceptable probabilities of $p_{\rm AD} > 0.05$ for subsamples of events with bright source stars ($I_{\rm S} \leq 17.75$) or $Spitzer$ coverage of the light curve peak. This indicates that the systematic errors have a more serious influence on faint events, especially when the light curve peak is not covered by $Spitzer$. We find that multiplying an error bar renormalization factor of 2.2 by the reported error bars on the $Spitzer$ microlensing parallax measurements provides reasonable agreement with all three Galactic models. However, corrections to the uncertainties in the $Spitzer$ photometry itself are a more effective way to address the systematic errors.


title: Broadband Coherent Perfect Absorber with PT-Symmetric 2D-Materials

abstract: We suggest graphene and a two-dimensional (2D) Weyl semimetal (WSM) as 2D materials for the realization of a broadband coherent perfect absorber (CPA) respecting overall PT-symmetry. We also demonstrate the conditions for mutually equal amplitudes and phases of the left and right incoming waves to realize a CPA. 2D materials in our system play the role to enhance the absorption rate of a CPA once the appropriate parameters are inserted in the system. We show that a 2D WSM is more effective than graphene in obtaining the optimal conditions. We display the behavior of each parameter governing the optical system and show that optimal conditions of these parameters give rise to enhancement and possible experimental realization of a broadband CPA-laser.


title: Approximating Orthogonal Matrices with Effective Givens Factorization

abstract: We analyze effective approximation of unitary matrices. In our formulation, a unitary matrix is represented as a product of rotations in two-dimensional subspaces, so-called Givens rotations. Instead of the quadratic dimension dependence when applying a dense matrix, applying such an approximation scales with the number factors, each of which can be implemented efficiently. Consequently, in settings where an approximation is once computed and then applied many times, such a representation becomes advantageous. Although effective Givens factorization is not possible for generic unitary operators, we show that minimizing a sparsity-inducing objective with a coordinate descent algorithm on the unitary group yields good factorizations for structured matrices. Canonical applications of such a setup are orthogonal basis transforms. We demonstrate numerical results of approximating the graph Fourier transform, which is the matrix obtained when diagonalizing a graph Laplacian.


title: Robust Precoding Design for Coarsely Quantized MU-MIMO Under Channel Uncertainties

abstract: Recently, multi-user multiple input multiple output (MU-MIMO) systems with low-resolution digital-to-analog converters (DACs) has received considerable attention, owing to the capability of dramatically reducing the hardware cost. Besides, it has been shown that the use of low-resolution DACs enable great reduction in power consumption while maintain the performance loss within acceptable margin, under the assumption of perfect knowledge of channel state information (CSI). In this paper, we investigate the precoding problem for the coarsely quantized MU-MIMO system without such an assumption. The channel uncertainties are modeled to be a random matrix with finite second-order statistics. By leveraging a favorable relation between the multi-bit DACs outputs and the single-bit ones, we first reformulate the original complex precoding problem into a nonconvex binary optimization problem. Then, using the S-procedure lemma, the nonconvex problem is recast into a tractable formulation with convex constraints and finally solved by the semidefinite relaxation (SDR) method. Compared with existing representative methods, the proposed precoder is robust to various channel uncertainties and is able to support a MUMIMO system with higher-order modulations, e.g., 16QAM.


title: Solving the gamma-ray radiative transfer equation for supernovae

abstract: We present a new relativistic radiative-transfer code for $\gamma$-rays of energy less than 5 MeV in supernova (SN) ejecta. This code computes the opacities, the prompt emissivity (i.e. decay), and the scattering emissivity, and solves for the intensity in the co-moving frame. Because of the large expansion velocities of SN ejecta, we ignore redistribution effects associated with thermal motions. The energy deposition is calculated from the energy removed from the radiation field by scattering or photoelectric absorption. This new code yields comparable results to an independent Monte Carlo code. However, both yield non-trivial differences with the results from a pure absorption treatment of $\gamma$-ray transport. A synthetic observer's frame spectrum is also produced from the CMF intensity. At early times when the optical depth to $\gamma$-rays is large, the synthetic spectrum show asymmetric line profiles with redshifted absorption as seen in SN 2014J. This new code is integrated within CMFGEN and allows for an accurate and fast computation of the decay energy deposition in SN ejecta.


title: Lasing with Topological Weyl Semimetal

abstract: Lasing behavior of optically active planar topological Weyl semimetal (TWS) is investigated in view of the Kerr and Faraday rotations. Robust topological character of TWS is revealed by the presence of Weyl nodes and relevant surface conductivities. We focus our attention on the surfaces where no Fermi arcs are formed, and thus Maxwell equations contain topological terms. We explicitly demonstrate that two distinct lasing modes arise because of the presence of effective refractive indices which lead to the birefringence phenomena. Transfer matrix is constructed in such a way that reflection and transmission amplitudes involve $2\times2$ matrix-valued components describing the bimodal character of the TWS laser. We provide associated parameters of the topological laser system yielding the optimal impacts. We reveal that gain values corresponding to the lasing threshold display a quantized behavior, which occurs due to topological character of the system. Our proposal is supported by the corresponding graphical demonstrations. Our observations and predictions suggest a concrete way of forming TWS laser and coherent perfect absorber; and are awaited to be confirmed by an experimental realization based on our computations.


abstract: Proving Unrealizability for Syntax-Guided Synthesis We consider the problem of automatically establishing that a given syntax-guided-synthesis (SyGuS) problem is unrealizable (i.e., has no solution). Existing techniques have quite limited ability to establish unrealizability for general SyGuS instances in which the grammar describing the search space contains infinitely many programs. By encoding the synthesis problem's grammar G as a nondeterministic program P_G, we reduce the unrealizability problem to a reachability problem such that, if a standard program-analysis tool can establish that a certain assertion in P_G always holds, then the synthesis problem is unrealizable. Our method can be used to augment any existing SyGus tool so that it can establish that a successfully synthesized program q is optimal with respect to some syntactic cost -- e.g., q has the fewest possible if-then-else operators. Using known techniques, grammar G can be automatically transformed to generate exactly all programs with lower cost than q -- e.g., fewer conditional expressions. Our algorithm can then be applied to show that the resulting synthesis problem is unrealizable. We implemented the proposed technique in a tool called NOPE. NOPE can prove unrealizability for 59/134 variants of existing linear-integer-arithmetic SyGus benchmarks, whereas all existing SyGus solvers lack the ability to prove that these benchmarks are unrealizable, and time out on them.


title: Preservation of normality by non-oblivious group selection

abstract: We give two different proofs of the fact that non-oblivious selection via regular group sets preserves normality. Non-oblivious here means that whether or not a symbol is selected can depend on the symbol itself. One proof relies on the incompressibility of normal sequences, the other on the use of augmented dynamical systems.


title: A sample-based spectral method approach for solving high-dimensional SPDEs

abstract: Uncertainty quantification appears today as a crucial point in numerous branches of science and engineering. In the past two decades, a growing interest has been devoted to stochastic finite element method (SFEM) for the propagation of uncertainties through physical models governed by stochastic partial differential equations (SPDEs). Despite its success and applications, the SFEM is mainly limited to small-scale and low-dimensional stochastic problems due to the extreme computational cost. In this article, by developing an universal construct of stochastic solution and a general algorithm for linear/nonlinear SFE equation, we explore a new strategy for the solution of high-dimensional stochastic problems, where stochastic problems are transformed into deterministic problems and stochastic algebraic equations. Since computational cost is almost proportional to the stochastic dimensionality of the problem, our method beats the so-called Curse of Dimensionality with great success. Numerical examples, including linear, nonlinear and high-dimensional stochastic problems, are used to demonstrate the method. Results show that our algorithm provides a highly efficient and unified framework for problems involving uncertainties, and is particularly appropriate for high stochastic dimensional problems of practical interest.


title: Multivariate Modeling for Sustainable and Resilient Infrastructure Systems and Communities

abstract: Sustainability and resilience of urban systems are multifaceted concepts, requiring information about multiple system attributes to adequately evaluate and characterize. However, despite the scientific consensus on the multivariate nature of these concepts, many of the existing techniques to model urban sustainability and resilience are unidimensional in nature, focusing on a characterizing a single element of highly interconnected urban systems. We champion a paradigm shift in modeling urban sustainability and resilience, using an integrated approach to simultaneously estimate multiple interconnected (correlated) system attributes of sustainability and resilience as a function of key environmental factors. We present a novel case study and review a few recent studies to illustrate the applicability and benefits of the multivariate approach to modeling urban sustainability and resilience. Our proposed framework can be utilized by infrastructure managers, urban planners, and researchers to conceptualize and assess urban sustainability and resilience more holistically, and to better understand the key factors in advancing the sustainability and resilience of infrastructure systems.


title: A Beurling-Lax-Halmos theorem for spaces with a complete Nevanlinna-Pick factor

abstract: We provide a short argument to establish a Beurling-Lax-Halmos theorem for reproducing kernel Hilbert spaces whose kernel has a complete Nevanlinna-Pick factor. We also record factorization results for pairs of nested invariant subspaces.


title: Elementary numerical methods for double integrals

abstract: Approximations to the integral $\int_a^b\int_c^d f(x,y)\,dy\,dx$ are obtained under the assumption that the partial derivatives of the integrand are in an $L^p$ space, for some $1\leq p\leq\infty$. We assume ${\lVert f_{xy}\rVert}_p$ is bounded (integration over $[a,b]\times[c,d]$), assume ${\lVert f_x(\cdot,c)\rVert}_p$ and ${\lVert f_x(\cdot,d)\rVert}_p$ are bounded (integration over $[a,b]$), and assume ${\lVert f_y(a,\cdot)\rVert}_p$ and ${\lVert f_y(b,\cdot)\rVert}_p$ are bounded (integration over $[c,d]$). The methods are elementary, using only integration by parts and H\"older's inequality. Versions of the trapezoidal rule, composite trapezoidal rule, midpoint rule and composite midpoint rule are given, with error estimates in terms of the above norms.


title: On spectral measures for certain unitary representations of R. Thompson's group F

abstract: The Hilbert space $\mathcal H$ of backward renormalisation of an anyonic quantum spin chain affords a unitary representation of Thompson's group $F$ via local scale transformations. Given a vector in the canonical dense subspace of $\mathcal H$ we show how to calculate the corresponding spectral measure for any element of $F$ and illustrate with some examples. Introducing the "essential part" of an element we show that the spectral measure of any vector in $\mathcal H$ is, apart from possibly finitely many eigenvalues, absolutely continuous with respect to Lebesgue measure. The same considerations and results hold for the Brown-Thompson groups $F_n$ (for which $F=F_2$).


title: 3d dualities from 2d free field correlators: recombination and rank stabilization

abstract: We propose various new 3d N=2 dualities exploiting their recently discovered connection to the duality relations for 2d free field CFT correlators. Most of the dualities involve, as the main building block, a quiver theory with monopole superpotential which enjoys various interesting properties such as being self-dual and reducing, in a suitable real mass deformation, to the familiar T[SU(N)] theory. In particular we propose a duality for the U(N) theory with one adjoint and k+1 fundamental flavors. By iterating some basic dualities we can bring the theory to a stable form which, in turns, allows us to find a dual frame where the rank of the original theory appears as a parameter.


title: The Radio Universe at Low Surface Brightness: Feedback & accretion in the circumgalactic medium

abstract: Massive galaxies at high-z are known to co-evolve with their circumgalactic medium (CGM). If we want to truly understand the role of the CGM in the early evolution of galaxies and galaxy-clusters, we need to fully explore the multi-phase nature of the CGM. We present two novel science cases that utilize low-surface-brightness observations in the radio regime to better understand the CGM around distant galaxies. At the lowest temperatures, observations of widespread molecular gas are providing evidence for the cold baryon cycle that grows massive galaxies. At the highest temperatures, observations of the Sunyaev-Zeldovich Effect are starting to reveal the effect of quasar feedback onto the hot gas in the CGM. We discuss the critical role that radio interferometers with compact configurations in the millimeter regime will play over the next decade in understanding the crucial role of the multi-phase CGM in galaxy evolution.


title: Learning Policies from Self-Play with Policy Gradients and MCTS Value Estimates

abstract: In recent years, state-of-the-art game-playing agents often involve policies that are trained in self-playing processes where Monte Carlo tree search (MCTS) algorithms and trained policies iteratively improve each other. The strongest results have been obtained when policies are trained to mimic the search behaviour of MCTS by minimising a cross-entropy loss. Because MCTS, by design, includes an element of exploration, policies trained in this manner are also likely to exhibit a similar extent of exploration. In this paper, we are interested in learning policies for a project with future goals including the extraction of interpretable strategies, rather than state-of-the-art game-playing performance. For these goals, we argue that such an extent of exploration is undesirable, and we propose a novel objective function for training policies that are not exploratory. We derive a policy gradient expression for maximising this objective function, which can be estimated using MCTS value estimates, rather than MCTS visit counts. We empirically evaluate various properties of resulting policies, in a variety of board games.


title: Pointing task on smart glasses: Comparison of four interaction techniques

abstract: Mobile devices such as smartphones, smartwatches or smart glasses have revolutionized how we interact. We are interested in smart glasses because they have the advantage of providing a simultaneous view of both physical and digital worlds. Despite this potential, pointing task on smart glasses is not really widespread. In this paper, we compared four interaction techniques for selecting targets : (a) the Absolute Head Movement and (b) the Relative Head Movement, where head movement controls the cursor on smart glasses in absolute or relative way, (c) the Absolute Free Hand interaction, where the forefinger control the cursor and (d) the Tactile Surface interaction, where the user controls the cursor via a small touchpad connected to smart glasses. We conducted an experiment with 18 participants. The Tactile Surface and Absolute Head Movement were the most efficient. The Relative Head Movement and Absolute Free Hand interactions were promising and require more exploration for other tasks.


title: Electromagnetic Backgrounds and Potassium-42 Activity in the DEAP-3600 Dark Matter Detector

abstract: The DEAP-3600 experiment is searching for WIMP dark matter with a 3.3 tonne single phase liquid argon (LAr) target, located 2.1 km underground at SNOLAB. The experimental signature of dark matter interactions is keV-scale $^{40}$Ar nuclear recoils (NR) producing 128 nm LAr scintillation photons observed by PMTs. The largest backgrounds in DEAP-3600 are electronic recoils (ER) induced by $\beta$ and $\gamma$-rays originating from internal and external radioactivity in the detector material. A background model of the ER interactions in DEAP-3600 was developed and is described in this work. The model is based on several components which are expected from radioisotopes in the LAr, from ex-situ material assay measurements, and from dedicated independent in-situ analyses. This prior information is used in a Bayesian fit of the ER components to a 247.2 d dataset to model the radioactivity in the surrounding detector materials. While excellent discrimination between ERs and NRs is reached with pulse shape discrimination, utilizing the large difference between fast and slow components of LAr scintillation light, detailed knowledge of the ER background and activity of detector components, sets valuable constraints on other key types of backgrounds in the detector: neutrons and alphas. In addition, the activity of $^{42}$Ar in LAr in DEAP-3600 is determined by measuring the daughter decay of $^{42}$K. This cosmogenically activated trace isotope is a relevant background at higher energies for other rare event searches using atmospheric argon e.g. DarkSide-20k, GERDA or LEGEND. The specific activity of $^{42}$Ar in the atmosphere is found to be $40.4 \pm 5.9$ $\mu$Bq/kg of argon.


title: Multi-task Learning for Multi-modal Emotion Recognition and Sentiment Analysis

abstract: Related tasks often have inter-dependence on each other and perform better when solved in a joint framework. In this paper, we present a deep multi-task learning framework that jointly performs sentiment and emotion analysis both. The multi-modal inputs (i.e., text, acoustic and visual frames) of a video convey diverse and distinctive information, and usually do not have equal contribution in the decision making. We propose a context-level inter-modal attention framework for simultaneously predicting the sentiment and expressed emotions of an utterance. We evaluate our proposed approach on CMU-MOSEI dataset for multi-modal sentiment and emotion analysis. Evaluation results suggest that multi-task learning framework offers improvement over the single-task framework. The proposed approach reports new state-of-the-art performance for both sentiment analysis and emotion analysis.


title: The connection between multiple prices of an Option at a given time with single prices defined at different times: The concept of weak-value in quantum finance

abstract: We introduce a new tool for predicting the evolution of an option for the cases where at some specific time, there is a high-degree of uncertainty for identifying its price. We work over the special case where we can predict the evolution of the system by joining a single price for the Option, defined at some specific time with a pair of prices defined at another instant. This is achieved by describing the evolution of the system through a financial Hamiltonian. The extension to the case of multiple prices at a given instant is straightforward. We also explain how to apply these results in real situations.


title: The paradox of monotone structural QRE

abstract: McKelvey and Palfrey (1995)'s monotone structural Quantal Response Equilibrium theory may be misspecified for the study of monotone behavior.


title: Investigating the origin of the spectral line profiles of the Hot Wolf-Rayet Star WR2

abstract: The hot WN star WR2 (HD6327) has been claimed to have many singular characteristics. To explain its unusually rounded and relatively weak emission line profiles, it has been proposed that WR2 is rotating close to break-up with a magnetically confined wind. Alternatively, the line profiles could be explained by the dilution of WR2's spectrum by that of a companion. In this paper, we present a study of WR2 using near-infrared AO imaging and optical spectroscopy and polarimetry. Our spectra reveal the presence of weak photospheric absorption lines from a ~B2.5-4V companion, which however contributes only ~5-10% to the total light, suggesting that the companion is a background object. Therefore, its flux cannot be causing any significant dilution of the WR star's emission lines. The absence of intrinsic linear continuum polarization from WR2 does not support the proposed fast rotation. Our Stokes V spectrum was not of sufficient quality to test the presence of a moderately strong organized magnetic field but our new modelling indicates that to confine the wind the putative magnetic field must be significantly stronger than was previously suggested sufficiently strong as to make its presence implausible.


title: Curriculum Learning for Domain Adaptation in Neural Machine Translation

abstract: We introduce a curriculum learning approach to adapt generic neural machine translation models to a specific domain. Samples are grouped by their similarities to the domain of interest and each group is fed to the training algorithm with a particular schedule. This approach is simple to implement on top of any neural framework or architecture, and consistently outperforms both unadapted and adapted baselines in experiments with two distinct domains and two language pairs.


title: 3D-VAR for Parametrized Partial Differential Equations: A Certified Reduced Basis Approach

abstract: In this paper, we propose a reduced order approach for 3D variational data assimilation governed by parametrized partial differential equations. In contrast to the classical 3D-VAR formulation that penalizes the measurement error directly, we present a modified formulation that penalizes the experimentally-observable misfit in the measurement space. Furthermore, we include a model correction term that allows to obtain an improved state estimate. We begin by discussing the influence of the measurement space on the amplification of noise and prove a necessary and sufficient condition for the identification of a "good" measurement space. We then propose a certified reduced basis (RB) method for the estimation of the model correction, the state prediction, the adjoint solution and the observable misfit with respect to the true state for real-time and many-query applications. A posteriori bounds are proposed for the error in each of these approximations. Finally, we introduce different approaches for the generation of the reduced basis spaces and the stability-based selection of measurement functionals. The 3D-VAR method and the associated certified reduced basis approximation are tested in a parameter and state estimation problem for a steady-state thermal conduction problem with unknown parameters and unknown Neumann boundary conditions.


title: Ontology-Aware Clinical Abstractive Summarization

abstract: Automatically generating accurate summaries from clinical reports could save a clinician's time, improve summary coverage, and reduce errors. We propose a sequence-to-sequence abstractive summarization model augmented with domain-specific ontological information to enhance content selection and summary generation. We apply our method to a dataset of radiology reports and show that it significantly outperforms the current state-of-the-art on this task in terms of rouge scores. Extensive human evaluation conducted by a radiologist further indicates that this approach yields summaries that are less likely to omit important details, without sacrificing readability or accuracy.


title: Tailoring the spatio-temporal distribution of diffractive focused ultrashort pulses through pulse shaping

abstract: Focusing control of ultrashort pulsed beams is an important research topic, due to its impact to subsequent interaction with matter. In this work, we study the propagation near the focus of ultrashort laser pulses of ~25 fs duration under diffractive focusing. We perform the spatio-spectral and spatio-temporal measurements of their amplitude and phase, complemented by the corresponding simulations. With them, we demonstrate that pulse shaping allows modifying in a controlled way not only the spatio-temporal distribution of the light irradiance in the focal region, but also the way it propagates as well as the frequency distribution within the pulse (temporal chirp). To gain a further intuitive insight, the role of diverse added spectral phase components is analyzed, showing the symmetries that arise for each case. In particular, we compare the effects, similarities and differences of the second and third order dispersion cases.


title: Reconstruction-Aware Imaging System Ranking by use of a Sparsity-Driven Numerical Observer Enabled by Variational Bayesian Inference

abstract: It is widely accepted that optimization of imaging system performance should be guided by task-based measures of image quality (IQ). It has been advocated that imaging hardware or data-acquisition designs should be optimized by use of an ideal observer (IO) that exploits full statistical knowledge of the measurement noise and class of objects to be imaged, without consideration of the reconstruction method. In practice, accurate and tractable models of the complete object statistics are often difficult to determine. Moreover, in imaging systems that employ compressive sensing concepts, imaging hardware and sparse image reconstruction are innately coupled technologies. In this work, a sparsity-driven observer (SDO) that can be employed to optimize hardware by use of a stochastic object model describing object sparsity is described and investigated. The SDO and sparse reconstruction method can therefore be "matched" in the sense that they both utilize the same statistical information regarding the class of objects to be imaged. To efficiently compute the SDO test statistic, computational tools developed recently for variational Bayesian inference with sparse linear models are adopted. The use of the SDO to rank data-acquisition designs in a stylized example as motivated by magnetic resonance imaging (MRI) is demonstrated. This study reveals that the SDO can produce rankings that are consistent with visual assessments of the reconstructed images but different from those produced by use of the traditionally employed Hotelling observer (HO).


title: A Chromaticity Analysis and PSF Subtraction Techniques for SCExAO/CHARIS data

abstract: We present an analysis of instrument performance using new observations taken with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS) instrument and the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system. In a correlation analysis of our datasets (which use the broadband mode covering J through K band in a single spectrum), we find that chromaticity in the SCExAO/CHARIS system is generally worse than temporal stability. We also develop a point spread function (PSF) subtraction pipeline optimized for the CHARIS broadband mode, including a forward modelling-based exoplanet algorithmic throughput correction scheme. We then present contrast curves using this newly developed pipeline. An analogous subtraction of the same datasets using only the H band slices yields the same final contrasts as the full JHK sequences; this result is consistent with our chromaticity analysis, illustrating that PSF subtraction using spectral differential imaging (SDI) in this broadband mode is generally not more effective than SDI in the individual J, H, or K bands. In the future, the data processing framework and analysis developed in this paper will be important to consider for additional SCExAO/CHARIS broadband observations and other ExAO instruments which plan to implement a similar integral field spectrograph broadband mode.


title: Performance Analysis of Non-DC-Biased OFDM

abstract: The performance analysis of a novel optical modulation scheme is presented in this paper. The basic concept is to transmit signs of modulated optical orthogonal frequency division multiplexing (O-OFDM) symbols and absolute values of the symbols separately by two information carrying units: 1) indices of two light emitting diodes (LED) transmitters that represent positive and negative signs separately; and 2) optical intensity symbols that carry the absolute values of signals. The new approach, referred as to non-DC-biased OFDM (NDC-OFDM), uses the optical spatial modulation (OSM) technique to eliminate the effect of the clipping distortion in DC-biased optical OFDM (DCO-OFDM). In addition, it can achieve similar advantages as the conventional unipolar modulation scheme, asymmetrically clipped optical OFDM (ACO-OFDM), without using additional subcarriers. In this paper, the analytical BER performance is compared with the Monte Carlo result in order to prove the reliability of the new method. Moreover, the practical BER performance of NDC-OFDM with DCO-OFDM and ACO-OFDM is compared for different constellation sizes to verify the improvement of NDC-OFDM on the spectral and power efficiencies.


title: An Analysis of a Fishing Model with Nonlinear Harvesting Function

abstract: In this study, considering the importance of how to exploit renewable natural resources, we analyze a fishing model with nonlinear harvesting function in which the players at the equilibrium point do a static game with complete information that, according to the calculations, will cause a waste of energy for both players and so the selection of cooperative strategies along with the agreement between the players is the result of this research.


title: Counterfactual Off-Policy Evaluation with Gumbel-Max Structural Causal Models

abstract: We introduce an off-policy evaluation procedure for highlighting episodes where applying a reinforcement learned (RL) policy is likely to have produced a substantially different outcome than the observed policy. In particular, we introduce a class of structural causal models (SCMs) for generating counterfactual trajectories in finite partially observable Markov Decision Processes (POMDPs). We see this as a useful procedure for off-policy "debugging" in high-risk settings (e.g., healthcare); by decomposing the expected difference in reward between the RL and observed policy into specific episodes, we can identify episodes where the counterfactual difference in reward is most dramatic. This in turn can be used to facilitate review of specific episodes by domain experts. We demonstrate the utility of this procedure with a synthetic environment of sepsis management.


title: A Rough Super-Brownian Motion

abstract: We study the scaling limit of a branching random walk in static random environment in dimension $d=1,2$ and show that it is given by a super-Brownian motion in a white noise potential. In dimension $1$ we characterize the limit as the unique weak solution to the stochastic PDE: \[\partial_t \mu = (\Delta {+} \xi) \mu {+} \sqrt{2\nu \mu} \tilde{\xi}\] for independent space white noise $\xi$ and space-time white noise $\tilde{\xi}$. In dimension $2$ the study requires paracontrolled theory and the limit process is described via a martingale problem. In both dimensions we prove persistence of this rough version of the super-Brownian motion.


title: The fundamental metallicity relation emerges from the local anti-correlation between star formation rate and gas-phase metallicity existing in disk galaxies

abstract: The fundamental metallicity relation (FMR) states that galaxies of the same stellar mass but larger star formation rate (SFR) tend to have smaller gas-phase metallicity (<Zg>). It is thought to be fundamental because it naturally arises from the stochastic feeding of star-formation from external metal-poor gas accretion, a process extremely elusive to observe but essential according the cosmological simulations of galaxy formation. In this letter, we show how the FMR emerges from the local anti-correlation between SFR surface density and Zg recently observed to exist in disk galaxies. We analytically derive the global FMR from the local law, and then show that both relations agree quantitatively when considering the star-forming galaxies of the MaNGA survey. Thus, understanding the FMR becomes equivalent to understanding the origin of the anti-correlation between SFR and metallicity followed by the set of star-forming regions of any typical galaxy. The correspondence between local and global laws is not specific of the FMR, so that a number of local relations should exist associated with known global relations.


title: Discriminative Sleep Patterns of Alzheimer's Disease via Tensor Factorization

abstract: Sleep change is commonly reported in Alzheimer's disease (AD) patients and their brain wave studies show decrease in dreaming and non-dreaming stages. Although sleep disturbance is generally considered as a consequence of AD, it might also be a risk factor of AD as new biological evidence shows. Leveraging National Sleep Research Resource (NSRR), we built a unique cohort of 83 cases and 331 controls with clinical variables and EEG signals. Supervised tensor factorization method was applied for this temporal dataset to extract discriminative sleep patterns. Among the 30 patterns extracted, we identified 5 significant patterns (4 patterns for AD likely and 1 pattern for normal ones) and their visual patterns provide interesting linkage to sleep with repeated wakefulness, insomnia, epileptic seizure, and etc. This study is preliminary but findings are interesting, which is a first step to provide quantifiable evidences to measure sleep as a risk factor of AD.


title: Minimax estimation of smooth optimal transport maps

abstract: Brenier's theorem is a cornerstone of optimal transport that guarantees the existence of an optimal transport map $T$ between two probability distributions $P$ and $Q$ over $\mathbb{R}^d$ under certain regularity conditions. The main goal of this work is to establish the minimax estimation rates for such a transport map from data sampled from $P$ and $Q$ under additional smoothness assumptions on $T$. To achieve this goal, we develop an estimator based on the minimization of an empirical version of the semi-dual optimal transport problem, restricted to truncated wavelet expansions. This estimator is shown to achieve near minimax optimality using new stability arguments for the semi-dual and a complementary minimax lower bound. Furthermore, we provide numerical experiments on synthetic data supporting our theoretical findings and highlighting the practical benefits of smoothness regularization. These are the first minimax estimation rates for transport maps in general dimension.


title: Modeling the Protoplanetary Disks of Two Brown Dwarfs in the Taurus Molecular Cloud

abstract: Measuring the properties of protoplanetary disks around brown dwarfs is central to understanding the formation of brown dwarfs and their planetary companions. We present modeling of CFHT Tau 4 and 2M0444, two brown dwarfs with protoplanetary disks in the Taurus Molecular Cloud. By combining modeling of the spectral energy distributions and ALMA images, we obtain disk radii and masses for these objects; these parameters can be used to constrain brown dwarf formation and planet formation, respectively. We find that the disk around CFHT Tau 4 has a mass of 0.42 M$_{Jup}$ and a radius of 80 au; we find 2M0444's disk to have a mass of 2.05 M$_{Jup}$ and a radius of 100 au. These radii are more consistent with those predicted by theoretical simulations of brown dwarf formation via undisturbed condensation from a mass reservoir than those predicted by ejection from the formation region. Furthermore, the disk mass of 2M0444 suggests that planet formation may be possible in this disk, although the disk of 2M0444 is likely not massive enough to form planets. The disk properties measured here provide constraints to theoretical models of brown dwarf formation and the formation of their planetary companions.


title: Combining Representation Learning with Tensor Factorization for Risk Factor Analysis - an application to Epilepsy and Alzheimer's disease

abstract: Existing studies consider Alzheimer's disease (AD) a comorbidity of epilepsy, but also recognize epilepsy to occur more frequently in patients with AD than those without. The goal of this paper is to understand the relationship between epilepsy and AD by studying causal relations among subgroups of epilepsy patients. We develop an approach combining representation learning with tensor factorization to provide an in-depth analysis of the risk factors among epilepsy patients for AD. An epilepsy-AD cohort of ~600,000 patients were extracted from Cerner Health Facts data (50M patients). Our experimental results not only suggested a causal relationship between epilepsy and later onset of AD ( p = 1.92e-51), but also identified five epilepsy subgroups with distinct phenotypic patterns leading to AD. While such findings are preliminary, the proposed method combining representation learning with tensor factorization seems to be an effective approach for risk factor analysis.


title: Effect of Strain on Charge Density Wave Order in the Holstein Model

abstract: We investigate charge ordering in the Holstein model in the presence of anisotropic hopping, $t_x, t_y=1-\delta, 1 + \delta$, as a model of the effect of strain on charge density wave (CDW) materials. Using Quantum Monte Carlo simulations, we show that the CDW transition temperature is relatively insensitive to moderate anisotropy $\delta \lesssim 0.3$, but begins to decrease more rapidly at $\delta \gtrsim 0.4$. However, the density correlations, as well as the kinetic energies parallel and perpendicular to the compressional axis, change significantly for moderate $\delta$. Accompanying mean-field theory calculations show a similar qualitative structure, with the transition temperature relatively constant at small $\delta$ and a more rapid decrease for larger strains. We also obtain the density of states $N(\omega)$, which provides clear signal of the charge ordering transition at large strain, where finite size scaling of the charge structure factor is extremely difficult because of the small value of the order parameter.


title: Visual binary stars with partially missing data: Introducing multiple imputation in astrometric analysis

abstract: Partial measurements of relative position are a relatively common event during the observation of visual binary stars. However, these observations are typically discarded when estimating the orbit of a visual pair. In this article we present a novel framework to characterize the orbits from a Bayesian standpoint, including partial observations of relative position as an input for the estimation of orbital parameters. Our aim is to formally incorporate the information contained in those partial measurements in a systematic way into the final inference. In the statistical literature, an imputation is defined as the replacement of a missing quantity with a plausible value. To compute posterior distributions of orbital parameters with partial observations, we propose a technique based on Markov chain Monte Carlo with multiple imputation. We present the methodology and test the algorithm with both synthetic and real observations, studying the effect of incorporating partial measurements in the parameter estimation. Our results suggest that the inclusion of partial measurements into the characterization of visual binaries may lead to a reduction in the uncertainty associated to each orbital element, in terms of a decrease in dispersion measures (such as the interquartile range) of the posterior distribution of relevant orbital parameters. The extent to which the uncertainty decreases after the incorporation of new data (either complete or partial) depends on how informative those newly-incorporated measurements are. Quantifying the information contained in each measurement remains an open issue.


title: Supervised Learning of the Next-Best-View for 3D Object Reconstruction

abstract: Motivated by the advances in 3D sensing technology and the spreading of low-cost robotic platforms, 3D object reconstruction has become a common task in many areas. Nevertheless, the selection of the optimal sensor pose that maximizes the reconstructed surface is a problem that remains open. It is known in the literature as the next-best-view planning problem. In this paper, we propose a novel next-best-view planning scheme based on supervised deep learning. The scheme contains an algorithm for automatic generation of datasets and an original three-dimensional convolutional neural network (3D-CNN) used to learn the next-best-view. Unlike previous work where the problem is addressed as a search, the trained 3D-CNN directly predicts the sensor pose. We present a comparison of the proposed network against a similar net, and we present several experiments of the reconstruction of unknown objects validating the effectiveness of the proposed scheme.


title: Evidence for crisis-induced intermittency during geomagnetic superchron transitions

abstract: The geomagnetic field's dipole undergoes polarity reversals in irregular time intervals. Particularly long periods (of the order of $10^7$yrs) without reversals, named superchrons, have occurred at least three times in history. We provide observational evidence for high non-Gaussianity in the vicinity of a transition to and from a geomagnetic superchron, consisting of a sharp increase in high-order moments (skewness and kurtosis) of the dipole's distribution. Such increase in the moments is a universal feature of crisis-induced intermittency in low-dimensional dynamical systems undergoing global bifurcations. This suggests temporal variation of the underlying parameters of the physical system. Through a low dimensional system that models the geomagnetic reversals we show that the increase in the high-order moments during transitions to geomagnetic superchrons is caused by the progressive destruction of global periodic orbits exhibiting both polarities as the system approaches a merging bifurcation. We argue that the non-gaussianity in this system is caused by the redistribution of the attractor around local cycles as global ones are destroyed.


title: Using Bursty Announcements for Detecting BGP Routing Anomalies

abstract: Despite the robust structure of the Internet, it is still susceptible to disruptive routing updates that prevent network traffic from reaching its destination. Our research shows that BGP announcements that are associated with disruptive updates tend to occur in groups of relatively high frequency, followed by periods of infrequent activity. We hypothesize that we may use these bursty characteristics to detect anomalous routing incidents. In this work, we use manually verified ground truth metadata and volume of announcements as a baseline measure, and propose a burstiness measure that detects prior anomalous incidents with high recall and better precision than the volume baseline. We quantify the burstiness of inter-arrival times around the date and times of four large-scale incidents: the Indosat hijacking event in April 2014, the Telecom Malaysia leak in June 2015, the Bharti Airtel Ltd. hijack in November 2015, and the MainOne leak in November 2018; and three smaller scale incidents that led to traffic interception: the Belarusian traffic direction in February 2013, the Icelandic traffic direction in July 2013, and the Russian telecom that hijacked financial services in April 2017. Our method leverages the burstiness of disruptive update messages to detect these incidents. We describe limitations, open challenges, and how this method can be used for routing anomaly detection.


title: Contractive projections and real positive maps on operator algebras

abstract: We study contractive projections, isometries, and real positive maps on algebras of operators on a Hilbert space. For example we find generalizations and variants of certain classical results on contractive projections on C*-algebras and JB-algebras due to Choi, Effros, St{\o}rmer, Friedman and Russo, and others. In fact most of our arguments generalize to contractive `real positive' projections on Jordan operator algebras, that is on a norm-closed space A of operators on a Hilbert space which are closedunder the Jordan product. We also prove many new general results on real positive maps which are foundational to the study of such maps, and of interest in their own right. We also prove a new Banach-Stone type theorem for isometries between operator algebras or Jordan operator algebras. An application of this is given to the characterization of symmetric real positive projections.


title: The Voronoi Cell in a saturated Circle Packing and an elementary proof of Thue's theorem

abstract: The famous Kepler conjecture has a less spectacular, two-dimensional equivalent: The theorem of Thue states that the densest circle packing in the Euclidean plane has a hexagonal structure. A common proof uses Voronoi cells and analyzes their area applying Jensen's inequality on convex functions to receive a local estimate which is globally valid. Based on the concept of Voronoi cells, we will introduce a new tessellation into so-called L-triangles which can be related to fundamental parallelograms of lattice circle packings. Therefore a globally disordered circle packing can be reduced to locally ordered configurations: We will show how the theorem of Lagrange on lattice circle packings can be applied to non-lattice circle packings. Thus we receive a new proof of Thue's theorem.


title: Magnetic phase diagram of the infinite-U Hubbard model with nearest- and next nearest-neighbor hoppings

abstract: We study the infinite-U Hubbard model on ladders of 2, 4 and 6 legs with nearest (t) and next-nearest (t') neighbor hoppings by means of the density-matrix renormalization group algorithm. In particular, we analyze the stability of the Nagaoka state for several values of t' when we vary the electron density $(\rho)$ from half-filling to the low-density limit. We build the two-dimensional phase diagram, where the fully spin-polarized and paramagnetic states prevail. We find that the inclusion of a non-frustrating next nearest neighbor hopping stabilizes the fully spin-polarized phase up until |t'/t|=0.5. Surprisingly, for this value of t', the ground state is fully spin-polarized for almost any electron density 1 $\gtrsim \rho \gtrsim$ 0, connecting the Nagaoka state to itinerant ferromagnetism at low density. Also, we find that the previously found checkerboard insulator phase at t'=0 and $\rho$=0.75 is unstable against t'.


title: Planted Hitting Set Recovery in Hypergraphs

abstract: In various application areas, networked data is collected by measuring interactions involving some specific set of core nodes. This results in a network dataset containing the core nodes along with a potentially much larger set of fringe nodes that all have at least one interaction with a core node. In many settings, this type of data arises for structures that are richer than graphs, because they involve the interactions of larger sets; for example, the core nodes might be a set of individuals under surveillance, where we observe the attendees of meetings involving at least one of the core individuals. We model such scenarios using hypergraphs, and we study the problem of core recovery: if we observe the hypergraph but not the labels of core and fringe nodes, can we recover the "planted" set of core nodes in the hypergraph? We provide a theoretical framework for analyzing the recovery of such a set of core nodes and use our theory to develop a practical and scalable algorithm for core recovery. The crux of our analysis and algorithm is that the core nodes are a hitting set of the hypergraph, meaning that every hyperedge has at least one node in the set of core nodes. We demonstrate the efficacy of our algorithm on a number of real-world datasets, outperforming competitive baselines derived from network centrality and core-periphery measures.


title: A Learning based Branch and Bound for Maximum Common Subgraph Problems

abstract: Branch-and-bound (BnB) algorithms are widely used to solve combinatorial problems, and the performance crucially depends on its branching heuristic.In this work, we consider a typical problem of maximum common subgraph (MCS), and propose a branching heuristic inspired from reinforcement learning with a goal of reaching a tree leaf as early as possible to greatly reduce the search tree size.Extensive experiments show that our method is beneficial and outperforms current best BnB algorithm for the MCS.


title: Efficient computation of mean reverting portfolios using cyclical coordinate descent

abstract: The econometric challenge of finding sparse mean reverting portfolios based on a subset of a large number of assets is well known. Many current state-of-the-art approaches fall into the field of co-integration theory, where the problem is phrased in terms of an eigenvector problem with sparsity constraint. Although a number of approximate solutions have been proposed to solve this NP-hard problem, all are based on relatively simple models and are limited in their scalability. In this paper we leverage information obtained from a heterogeneous simultaneous graphical dynamic linear model (H-SGDLM) and propose a novel formulation of the mean reversion problem, which is phrased in terms of a quasi-convex minimisation with a normalisation constraint. This new formulation allows us to employ a cyclical coordinate descent algorithm for efficiently computing an exact sparse solution, even in a large universe of assets, while the use of H-SGDLM data allows us to easily control the required level of sparsity. We demonstrate the flexibility, speed and scalability of the proposed approach on S\&P$500$, FX and ETF futures data.


title: The Penetration Rate Effect of Connected and Automated Vehicles in Mixed Traffic Routing

abstract: We study the problem of routing Connected and Automated Vehicles (CAVs) in the presence of mixed traffic (coexistence of regular vehicles and CAVs). In this setting, we assume that all CAVs belong to the same fleet, and can be routed using a centralized controller. The routing objective is to minimize a given overall fleet traveling cost (travel time or energy consumption). We assume that regular vehicles (non-CAVs) choose their routing decisions selfishly to minimize their traveling time. We propose an algorithm that deals with the routing interaction between CAVs and regular uncontrolled vehicles. We investigate the effect of assigning system-centric routes under different penetration rates (fractions) of CAVs. To validate our method, we apply the proposed routing algorithms to the Braess Network and to a sub-network of the Eastern Massachusetts (EMA) transportation network using actual traffic data provided by the Boston Region Metropolitan Planning Organization. The results suggest that collaborative routing decisions of CAVs improve not only the cost of CAVs, but also that of the non-CAVs. Furthermore, even a small CAV penetration rate can ease congestion for the entire network.


title: Task-Driven Data Verification via Gradient Descent

abstract: We introduce a novel algorithm for the detection of possible sample corruption such as mislabeled samples in a training dataset given a small clean validation set. We use a set of inclusion variables which determine whether or not any element of the noisy training set should be included in the training of a network. We compute these inclusion variables by optimizing the performance of the network on the clean validation set via "gradient descent on gradient descent" based learning. The inclusion variables as well as the network trained in such a way form the basis of our methods, which we call Corruption Detection via Gradient Descent (CDGD). This algorithm can be applied to any supervised machine learning task and is not limited to classification problems. We provide a quantitative comparison of these methods on synthetic and real world datasets.


title: Faraday instability on a sphere: Floquet analysis

abstract: Standing waves appear at the surface of a spherical viscous liquid drop subjected to radial parametric oscillation. This is the spherical analogue of the Faraday instability. Modifying the Kumar & Tuckerman (1994) planar solution to a spherical interface, we linearize the governing equations about the state of rest and solve the resulting equations by using a spherical harmonic decomposition for the angular dependence, spherical Bessel functions for the radial dependence, and a Floquet form for the temporal dependence. Although the inviscid problem can, like the planar case, be mapped exactly onto the Mathieu equation, the spherical geometry introduces additional terms into the analysis. The dependence of the threshold on viscosity is studied and scaling laws are found. It is shown that the spherical thresholds are similar to the planar infinite-depth thresholds, even for small wavenumbers for which the curvature is high. A representative time-dependent Floquet mode is displayed.


title: A note on time analyticity for ancient solutions of the heat equation

abstract: It is well known that generic solutions of the heat equation are not analytic in time in general. Here it is proven that ancient solutions with exponential growth are analytic in time in ${\M} \times (-\infty, 0]$. Here $\M=\R^n$ or is a manifold with Ricci curvature bounded from below. Consequently a necessary and sufficient condition is found on the solvability of backward heat equation in the class of functions with exponential growth.


title: Fast Wiener filtering of CMB maps with Neural Networks

abstract: We show how a neural network can be trained to Wiener filter masked CMB maps to high accuracy. We propose an innovative neural network architecture, the WienerNet, which guarantees linearity in the data map. Our method does not require Wiener filtered training data, but rather learns Wiener filtering from tailored loss functions which are mathematically guaranteed to be minimized by the exact solution. Once trained, the neural network Wiener filter is extremely fast, about a factor of 1000 faster than the standard conjugate gradient method. Wiener filtering is the computational bottleneck in many optimal CMB analyses, including power spectrum estimation, lensing and non-Gaussianities, and our method could potentially be used to speed them up by orders of magnitude with minimal loss of optimality. The method should also be useful to analyze other statistical fields in cosmology.


title: Singularity Properties of Graph Varieties

abstract: For a graph $G=(V,E)$, and a symplectic vector space $(W, \left<\cdot,\cdot\right>)$, we define a variety $X(G,W)$ consisting of all functions $w:V\to W$ satisfying $\left<w(u), w(v)\right> = 0$ for any edge $\{u,v\}$ in $G$. We study the singularities of this varieties.


title: Vignetted-aperture correction for spectral cameras with integrated thin-film Fabry-P\'erot filters

abstract: Spectral cameras with integrated thin-film Fabry-P\'erot filters have become increasingly important in many applications. These applications often require the detection of spectral features at specific wavelengths or to quantify small variations in the spectrum. This can be challenging since thin-film filters are sensitive to the angle of incidence of the light. In prior work we modeled and corrected for the distribution of incident angles for an ideal finite aperture. Many real lenses however experience vignetting. Therefore in this article we generalize our model to the more common case of a vignetted aperture, which changes the distribution of incident angles. We propose a practical method to estimate the model parameters and correct undesired shifts in measured spectra. This is experimentally validated for a lens mounted on a visible to near-infrared spectral camera.


title: Consensus-based Interpretable Deep Neural Networks with Application to Mortality Prediction

abstract: Deep neural networks have achieved remarkable success in various challenging tasks. However, the black-box nature of such networks is not acceptable to critical applications, such as healthcare. In particular, the existence of adversarial examples and their overgeneralization to irrelevant, out-of-distribution inputs with high confidence makes it difficult, if not impossible, to explain decisions by such networks. In this paper, we analyze the underlying mechanism of generalization of deep neural networks and propose an ($n$, $k$) consensus algorithm which is insensitive to adversarial examples and can reliably reject out-of-distribution samples. Furthermore, the consensus algorithm is able to improve classification accuracy by using multiple trained deep neural networks. To handle the complexity of deep neural networks, we cluster linear approximations of individual models and identify highly correlated clusters among different models to capture feature importance robustly, resulting in improved interpretability. Motivated by the importance of building accurate and interpretable prediction models for healthcare, our experimental results on an ICU dataset show the effectiveness of our algorithm in enhancing both the prediction accuracy and the interpretability of deep neural network models on one-year patient mortality prediction. In particular, while the proposed method maintains similar interpretability as conventional shallow models such as logistic regression, it improves the prediction accuracy significantly.


title: Primordial scalar gravitational waves produced at the QCD phase transition due to the trace anomaly

abstract: Relying only on the standard model of elementary particles and gravity, we study the details of a new source of gravitational waves whose origin is in quantum physics. Namely, it is well known that massless fields in curved backgrounds suffer from the so-called "trace anomaly". This anomaly can be cast in terms of new scalar degrees of freedom which take account of macroscopic effects of quantum matter in gravitational fields. The linearized effective action for these fields describes scalar (as opposed to transverse) gravitational waves, which are absent in Einstein's theory. Since these new degrees of freedom couple directly to the gauge field scalars in QCD, the epoch of the QCD phase transition in early universe is a possible source of primordial cosmological gravitational radiation. While the anomaly is most likely fully unsuppressed at the QCD densities (temperature is much higher than the u and d quark masses), just to be careful we introduced the window function which cuts-off very low frequencies where the anomaly effect might be suppressed. We then calculated the characteristic strain of the properly adjusted gravitational waves signal today. The region of the parameter space with no window function gives a stronger signal, and both the strain and the frequencies fall within the sensitivity of the near future gravitational wave experiments (e.g. LISA and The Big Bang Observer). The possibility that one can study quantum physics with gravitational wave astronomy even in principle is exciting, and will be of value for future endeavors in this field.


title: Exploring Interactions with Voice-Controlled TV

abstract: Intelligent agents such as Alexa, Siri, and Google Assistant are now built into streaming TV systems, allowing people to use voice input to navigate the increasingly complex set of apps available on a TV. However, these systems typically support a narrow range of control- and search-oriented commands, and do not support deeper recommendation or exploration queries. To learn about how people interact with a recommendation-oriented voice-controlled TV, we use research through design methods to explore an early prototype movie recommendation system where the only input modality is voice. We describe in-depth qualitative research sessions with 11 participants. We contribute implications for designers of voice-controlled TV: mitigating the drawbacks of voice-only interactions, navigating the tension between expressiveness and efficiency, and building voice-driven recommendation interfaces that facilitate exploration.


title: Influence of transverse field on the spin-3/2 Blume-Capel model on rectangular lattice

abstract: Transverse field effect on thermodynamic properties of the spin-3/2 Blume-Capel model on rectangular lattice in which the interactions in perpendicular directions differ in signs is studied within the mean field approximation. Phase diagrams in the (transverse field, temperature) plane are constructed for various values of single-ion anisotropy.


title: Reconstructing high-dimensional Hilbert-valued functions via compressed sensing

abstract: We present and analyze a novel sparse polynomial technique for approximating high-dimensional Hilbert-valued functions, with application to parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our theoretical framework treats the function approximation problem as a joint sparse recovery problem, where the set of jointly sparse vectors is possibly infinite. To achieve the simultaneous reconstruction of Hilbert-valued functions in both parametric domain and Hilbert space, we propose a novel mixed-norm based $\ell_1$ regularization method that exploits both energy and sparsity. Our approach requires extensions of concepts such as the restricted isometry and null space properties, allowing us to prove recovery guarantees for sparse Hilbert-valued function reconstruction. We complement the enclosed theory with an algorithm for Hilbert-valued recovery, based on standard forward-backward algorithm, meanwhile establishing its strong convergence in the considered infinite-dimensional setting. Finally, we demonstrate the minimal sample complexity requirements of our approach, relative to other popular methods, with numerical experiments approximating the solutions of high-dimensional parameterized elliptic PDEs.


title: On Structured Lyapunov Functions and Dissipativity in Interconnected LTI Systems

abstract: In this paper we study connections between structured storage or Lyapunov functions of a class of interconnected systems (dynamical networks) and dissipativity properties of the individual systems. We prove that if a dynamical network, composed as a set of linear time invariant (LTI) systems interconnected over an acyclic graph, admits an additive quadratic Lyapunov function, then the individual systems in the network are dissipative with respect to a (nonempty) set of interconnection neutral supply functions. Each supply function from this set is defined on a single interconnection link in the network. Specific characterizations of neutral supply functions are presented which imply robustness of network stability/dissiptivity to removal of interconnection links.


title: New perspectives on superfluidity in resonantly--driven polariton fluids

abstract: In this paper we discuss, within the Gross--Pitaevskii framework, superfluidity, soliton nucleation, and instabilities in a non-equilibrium polariton fluid injected by a spatially localized and continuous-wave coherent pump and flowing against a defect located outside the pump spot. In contrast to equilibrium condensates, the steady-state solutions of the driven-dissipative equations in this specific geometry hardly show a clean superfluid flow around the defect and rather feature a crossover from shallow to deep soliton-like perturbation. This is explained in terms of the properties of one-dimensional flows, in particular their weak dependence on the pump parameters and their rapid transition to a super-sonic regime under the effect of the quantum pressure; such a highly nonlinear behaviour calls for quantitative experimental tests of the underlying Gross--Pitaevskii equation. The role of disorder and of a incoherent reservoir in inducing non-stationary behaviours with moving vortices is also highlighted.


title: Single-photon-level light storage in cold atoms using the Autler-Townes splitting protocol

abstract: Broadband spin-photon interfaces for long-lived storage of photonic quantum states are key elements for quantum information technologies. Yet, reliable operation of such memories in the quantum regime is challenging due to photonic noise arising from technical and/or fundamental limitations in the storage-and-recall processes controlled by strong electromagnetic fields. Here, we experimentally implement a single-photon-level spin-wave memory in a laser-cooled Rubidium gas, based on the recently proposed Autler-Townes splitting (ATS) protocol. We demonstrate storage of 20-ns-long laser pulses, each containing an average of 0.1 photons, for 200 ns with an efficiency of $12.5\%$ and signal-to-noise ratio above 30. Notably, the robustness of ATS spin-wave memory against motional dephasing allows for an all-spatial filtering of the control-field noise, yielding an ultra-low unconditional noise probability of $3.3\times10^{-4}$, without the complexity of spectral filtering. These results highlight that broadband ATS memory in ultracold atoms is a preeminent option for storing quantum light.

categories: quant-ph physics.app-ph physics.atom-ph physics.optics


title: Variational Regret Bounds for Reinforcement Learning

abstract: We consider undiscounted reinforcement learning in Markov decision processes (MDPs) where both the reward functions and the state-transition probabilities may vary (gradually or abruptly) over time. For this problem setting, we propose an algorithm and provide performance guarantees for the regret evaluated against the optimal non-stationary policy. The upper bound on the regret is given in terms of the total variation in the MDP. This is the first variational regret bound for the general reinforcement learning setting.


title: Hyperbolic Nodal Band Structures and Knot Invariants

abstract: We extend the list of known band structure topologies to include a large family of hyperbolic nodal links and knots, occurring both in conventional Hermitian systems where their stability relies on discrete symmetries, and in the dissipative non-Hermitian realm where the knotted nodal lines are generic and thus stable towards any small perturbation. We show that these nodal structures, taking the forms of Turk's head knots, appear in both continuum- and lattice models with relatively short-ranged hopping that is within experimental reach. To determine the topology of the nodal structures, we devise an efficient algorithm for computing the Alexander polynomial, linking numbers and higher order Milnor invariants based on an approximate and well controlled parameterisation of the knot.

categories: cond-mat.mes-hall math-ph math.AT math.MP physics.optics quant-ph


title: Alternative Decohering Histories in Quantum Mechanics

abstract: We continue our efforts to understand, within the framework of the quantum mechanics of the universe as a whole, the quasiclassical realm of familiar experience as a feature emergent from the Hamiltonian of the elementary particles and the initial condition of the universe. Quantum mechanics assigns probabilities to exhaustive sets of alternative decoherent histories of the universe. We introduce and define the notion of strong decoherence. We replace the notion of maximal sets of alternative decohering histories by defining the more useful concept of "full" sets of alternative strongly decohering histories. These full sets fall into equivalence classes each of which is characterized by a basis in Hilbert space. Finally we describe our continuing efforts to find measures of classicality --- measures that could be applied to such full sets of alternative strongly decohering so as to characterize a quasiclassical realm.


title: Synchronization to big-data: nudging the Navier-Stokes equations for data assimilation of turbulent flows

abstract: Nudging is an important data assimilation technique where partial field measurements are used to control the evolution of a dynamical system and/or to reconstruct the entire phase-space configuration of the supplied flow. Here, we apply it to the toughest problem in fluid dynamics: three dimensional homogeneous and isotropic turbulence. By doing numerical experiments we perform a systematic assessment of how well the technique reconstructs large- and small-scales features of the flow with respect to the quantity and the quality/type of data supplied to it. The types of data used are: (i) field values on a fixed number of spatial locations (Eulerian nudging), (ii) Fourier coefficients of the fields on a fixed range of wavenumbers (Fourier nudging), or (iii) field values along a set of moving probes inside the flow (Lagrangian nudging). We present state-of-the-art quantitative measurements of the scale-by-scale {\it transition to synchronization} and a detailed discussion of the probability distribution function of the reconstruction error, by comparing the nudged field and the {\it truth} point-by-point. Furthermore, we show that for more complex flow configurations, like the case of anisotropic rotating turbulence, the presence of cyclonic and anticyclonic structures leads to unexpectedly better performances of the algorithm. We discuss potential further applications of nudging to a series of applied flow configurations, including the problem of field-reconstruction in thermal Rayleigh-B\'enard convection and in magnetohydrodynamics (MHD), and to the determination of optimal parametrisation for small-scale turbulent modeling. Our study fixes the standard requirements for future applications of nudging to complex turbulent flows.

categories: physics.flu-dyn physics.ao-ph physics.data-an physics.geo-ph


title: From Brain Imaging to Graph Analysis: a study on ADNI's patient cohort

abstract: In this paper, we studied the association between the change of structural brain volumes to the potential development of Alzheimer's disease (AD). Using a simple abstraction technique, we converted regional cortical and subcortical volume differences over two time points for each study subject into a graph. We then obtained substructures of interest using a graph decomposition algorithm in order to extract pivotal nodes via multi-view feature selection. Intensive experiments using robust classification frameworks were conducted to evaluate the performance of using the brain substructures obtained under different thresholds. The results indicated that compact substructures acquired by examining the differences between patient groups were sufficient to discriminate between AD and healthy controls with an area under the receiver operating curve of 0.72.


title: Optimizing Ground-based Observations of O2 in Earth Analogs

abstract: We present the result of calculations to optimize the search for molecular oxygen (O2) in Earth analogs transiting around nearby, low-mass stars using ground-based, high-resolution, Doppler shift techniques. We investigate a series of parameters, namely spectral resolution, wavelength coverage of the observations, and sky coordinates and systemic velocity of the exoplanetary systems, to find the values that optimize detectability of O2. We find that increasing the spectral resolution of observations to R = 300,000 - 400,000 from the typical R ~ 100,000, more than doubles the average depth of O2 lines in planets with atmospheres similar to Earth's. Resolutions higher than about 500,000 do not produce significant gains in the depths of the O2 lines. We confirm that observations in the O2 A-band are the most efficient except for M9V host stars, for which observations in the O2 NIR-band are more efficient. Combining observations in the O2 A, B, and NIR -bands can reduce the number of transits needed to produce a detection of O2 by about 1/3 in the case of white noise limited observations. However, that advantage disappears in the presence of typical levels of red noise. Therefore, combining observations in more than one band produces no significant gains versus observing only in the A-band, unless red-noise can be significantly reduced. Blending between the exoplanet's O2 lines and telluric O2 lines is a known problem. We find that problem can be alleviated by increasing the resolution of the observations, and by giving preference to targets near the ecliptic.


title: Convergence and Synchronization in Networks of Piecewise-Smooth Systems via Distributed Discontinuous Coupling

abstract: Complex networks are a successful framework to describe collective behaviour in many applications, but a notable gap remains in the current literature, that of proving asymptotic convergence in networks of piecewise-smooth systems. Indeed, a wide variety of physical systems display discontinuous dynamics that change abruptly, including dry friction mechanical oscillators, electrical power converters, and biological neurons. In this paper, we study how to enforce global asymptotic state-synchronization in these networks. Specifically, we propose the addition of a distributed discontinuous coupling action to the commonly used diffusive coupling protocol. Moreover, we provide analytical estimates of the thresholds on the coupling gains required for convergence, and highlight the importance of a new connectivity measure, which we named minimum density. The theoretical results are illustrated by a set of representative examples.


title: Constraints and Horizons for de Sitter with Extra Dimensions

abstract: In order for spacetimes with static extra dimensions to have 4-dimensional de Sitter expansion they must have at least positive curvature, warping sourced by the 4-d expansion, or violate the null energy condition everywhere in the extra dimensions. We show how this constraint arises from the null Raychaudhuri equation, and that it is independent of the matter content, the Einstein equations, and is true point-by-point in the extra dimensions (not integrated), setting it apart from other no-go theorems in the literature. We present two previously known examples -- a Freund-Rubin compactification with bulk cosmological constant, and a Randall-Sundrum model where the mismatch between the brane and bulk cosmological constants sources non-trivial warping -- which evade the constraint, and discuss the implications for flux compactifications. We also show that any spacetime with static compact extra dimensions and 4-dimensional de Sitter expansion has an apparent horizon and antitrapped region in the extra dimensions, which may have interesting implications for the dynamics of bulk fields in cosmology.


title: Nonlinear Semi-Parametric Models for Survival Analysis

abstract: Semi-parametric survival analysis methods like the Cox Proportional Hazards (CPH) regression (Cox, 1972) are a popular approach for survival analysis. These methods involve fitting of the log-proportional hazard as a function of the covariates and are convenient as they do not require estimation of the baseline hazard rate. Recent approaches have involved learning non-linear representations of the input covariates and demonstrate improved performance. In this paper we argue against such deep parameterizations for survival analysis and experimentally demonstrate that more interpretable semi-parametric models inspired from mixtures of experts perform equally well or in some cases better than such overly parameterized deep models.


title: Gray s decomposition on doubly warped product manifolds and applications

abstract: A. Gray presented an interesting $O\left( n\right) $ invariant decomposition of the covariant derivative of the Ricci tensor. Manifolds whose Ricci tensor satisfies the defining property of each orthogonal class are called Einstein-like manifolds. In the present paper, we answered the following question: Under what condition(s), does a factor manifold $M_{i},i=1,2$ of a doubly warped product manifold $M=_{f_{2}}M_{1}\times _{f_{1}}M_{2}$ lie in the same Einstein-like class of $M$? By imposing sufficient and necessary conditions on the warping functions, an inheritance property of each class is proved. As an application, Einstein-like doubly warped product space-times of type $\mathcal{A},$ $\mathcal{B}$ or $\mathcal{P}$ are considered.


title: On Borel subalgebras of quantum groups

abstract: For a quantum group, we study those right coideal subalgebras, for which all irreducible representations are one-dimensional. If a right coideal subalgebra is maximal with this property, then we call it a Borel subalgebra. Besides the positive part of the quantum group and its reflections, we find new unfamiliar Borel subalgebras, for example, ones containing copies of the quantum Weyl algebra. Given a Borel subalgebra, we study its induced (Verma-)modules and prove among others that they have all irreducible finite-dimensional modules as quotients. We give two structural conjectures involving the associated graded right coideal subalgebra, which we prove in certain cases. In particular, they predict the shape of all triangular Borel subalgebras. As examples, we determine all Borel subalgebras of $U_q(\mathfrak{sl}_2)$ and $U_q(\mathfrak{sl}_3)$ and discuss the induced modules.


title: Graph-Theoretic Stability Conditions for Metzler Matrices and Monotone Systems

abstract: This paper studies the graph-theoretic conditions for stability of positive monotone systems. Using concepts from input-to-state stability and network small-gain theory, we first establish necessary and sufficient conditions for the stability of linear positive systems described by Metzler matrices. Specifically, we derive two sets of stability conditions based on two forms of input-to-state stability gains for Metzler systems, namely max-interconnection gains and sum-interconnection gains. Based on the max-interconnection gains, we show that the cyclic small-gain theorem becomes necessary and sufficient for the stability of Metzler systems; based on the sum-interconnection gains, we obtain novel graph-theoretic conditions for the stability of Metzler systems. All these conditions highlight the role of cycles in the interconnection graph and unveil how the structural properties of the graph affect stability. Finally, we extend our results to the nonlinear monotone system and obtain similar sufficient conditions for global asymptotic stability.


title: Average Transition Conditions for Electromagnetic Fields at a Metascreen of Vanishing Thickness

abstract: Using a dipole interaction model, we derive generalized sheet transition conditions (GSTCs) for electromagnetic fields at the surface of a metascreen consisting of an array of arbitrarily shaped apertures in a perfectly conducting screen of zero thickness. Use of the GSTCs permits modeling of structures containing perforated surfaces much more rapidly than is possible with full-wave numerical simulations. These conditions require that the period of the array be smaller than about a third of a wavelength in the surrounding media, and generalize many results previously presented in the literature. They are validated by comparison with results of finite-element modeling, and show excellent agreement when conditions of their derivation are satisfied.


title: Orbital stability in the Solar System for arbitrary inclinations and eccentricities: planetary perturbations versus resonances

abstract: Applying the technique of dynamical maps we study the orbital stability of test particles in the Solar System in the space (a,e,i) defined by 0.1<a<38 au, 0<e<0.9 and 0<i<180 identifying the unstable and stable regions. We find stable niches where small bodies can survive even for very high eccentricities. Mean motion resonances play a fundamental role providing stability against the planetary perturbations specially for high inclination orbits. A stability stripe around i=150 is present all along the Solar System. We found that the population of objects with semimajor axes between 10 and 30 au is evolving inside a highly unstable region according to our maps. For the inner Solar System we found that the region between the Hildas and Jupiter is more stable for high eccentricity orbits than for low eccentricity ones.


title: Average Transition Conditions for Electromagnetic Fields at a Metascreen of Nonzero Thickness

abstract: Using a dipole interaction model, we derive generalized sheet transition conditions (GSTCs) for electromagnetic fields at the surface of a metascreen consisting of an array of arbitrarily shaped apertures in a perfectly conducting screen of nonzero thickness. The simple analytical formulas obtained are validated through comparison with full-wave numerical simulations.


title: An Algebraic Classification of Exceptional EFTs Part II: Supersymmetry

abstract: We present a novel approach to classify supersymmetric effective field theories (EFTs) whose scattering amplitudes exhibit enhanced soft limits. These enhancements arise due to non-linearly realised symmetries on the Goldstone modes of such EFTs and we classify the algebras that these symmetries can form. Our main focus is on so-called exceptional algebras which lead to field-dependent transformation rules and EFTs with the maximum possible soft enhancement at a given derivative power counting. We adapt existing techniques for Poincar\'{e} invariant theories to the supersymmetric case, and introduce superspace inverse Higgs constraints as a method of reducing the number of Goldstone modes while maintaining all symmetries. Restricting to the case of a single Goldstone supermultiplet in four dimensions, we classify the exceptional algebras and EFTs for a chiral, Maxwell or real linear supermultiplet. Moreover, we show how our algebraic approach allows one to read off the soft weights of the different component fields from superspace inverse Higgs trees, which are the algebraic cousin of the on-shell soft data one provides to soft bootstrap EFTs using on-shell recursion. Our Lie-superalgebraic approach extends the results of on-shell methods and provides a complementary perspective on non-linear realisations.


title: $Spitzer$ Parallax of OGLE-2018-BLG-0596: A Low-mass-ratio Planet around an M-dwarf

abstract: We report the discovery of a $Spitzer$ microlensing planet OGLE-2018-BLG-0596Lb, with preferred planet-host mass ratio $q \sim 2\times10^{-4}$. The planetary signal, which is characterized by a short $(\sim 1~{\rm day})$ "bump" on the rising side of the lensing light curve, was densely covered by ground-based surveys. We find that the signal can be explained by a bright source that fully envelops the planetary caustic, i.e., a "Hollywood" geometry. Combined with the source proper motion measured from $Gaia$, the $Spitzer$ satellite parallax measurement makes it possible to precisely constrain the lens physical parameters. The preferred solution, in which the planet perturbs the minor image due to lensing by the host, yields a Uranus-mass planet with a mass of $M_{\rm p} = 13.9\pm1.6~M_{\oplus}$ orbiting a mid M-dwarf with a mass of $M_{\rm h} = 0.23\pm0.03~M_{\odot}$. There is also a second possible solution that is substantially disfavored but cannot be ruled out, for which the planet perturbs the major image. The latter solution yields $M_{\rm p} = 1.2\pm0.2~M_{\oplus}$ and $M_{\rm h} = 0.15\pm0.02~M_{\odot}$. By combining the microlensing and $Gaia$ data together with a Galactic model, we find in either case that the lens lies on the near side of the Galactic bulge at a distance $D_{\rm L} \sim 6\pm1~{\rm kpc}$. Future adaptive optics observations may decisively resolve the major image/minor image degeneracy.


title: On the Convergence Rate of Variants of the Conjugate Gradient Algorithm in Finite Precision Arithmetic

abstract: We consider three mathematically equivalent variants of the conjugate gradient (CG) algorithm and how they perform in finite precision arithmetic. It was shown in [{\em Behavior of slightly perturbed Lanczos and conjugate-gradient recurrences}, Lin.~Alg.~Appl., 113 (1989), pp.~7-63] that under certain conditions the convergence of a slightly perturbed CG computation is like that of exact CG for a matrix with many eigenvalues distributed throughout tiny intervals about the eigenvalues of the given matrix, the size of the intervals being determined by how closely these conditions are satisfied. We determine to what extent each of these variants satisfies the desired conditions, using a set of test problems and show that there is significant correlation between how well these conditions are satisfied and how well the finite precision computation converges before reaching its ultimately attainable accuracy. We show that for problems where the width of the intervals containing the eigenvalues of the associated exact CG matrix makes a significant difference in the behavior of exact CG, the different CG variants behave differently in finite precision arithmetic. For problems where the interval width makes little difference or where the convergence of exact CG is essentially governed by the upper bound based on the square root of the condition number of the matrix, the different CG variants converge similarly in finite precision arithmetic until the ultimate level of accuracy is achieved, although this ultimate level of accuracy may be different for the different variants. This points to the need for testing new CG variants on problems that are especially sensitive to rounding errors.


title: Multiwavelength radio observations of a Brightest Cluster Galaxy at z=1.71: Detection of a modest Active Galactic Nucleus and evidence for extended star formation

abstract: We present deep, multiwavelength radio observations of SpARCS104922.6+564032.5, a z = 1.71 galaxy cluster with a starbusting core. Observations were made with the Karl G. Jansky Very Large Array (JVLA) in 3 bands: 1-2 GHz, 4-8 GHz and 8-12 GHz. We detect a radio source coincident with the Brightest Cluster Galaxy (BCG) that has a spectral index of {\alpha}=0.44\pm 0.29 and is indicative of emission from an Active Galactic Nucleus. The radio luminosity is consistent with the average luminosity of the lower redshift BCG sample, but the flux densities are 6{\sigma} below the predicted values of the star-forming Spectral Energy Distribution based on far infrared data. Our new fit fails to simultaneously describe the far infrared and radio fluxes. This, coupled with the fact that no other bright source is detected in the vicinity of the BCG implies that the star formation region, traced by the infrared emission, is extended or clumpy and not located directly within the BCG. Thus, we suggest that the star-forming core might not be driven by a single major wet merger, but rather by several smaller galaxies stripped of their gas or by a displaced cooling flow, although more data are needed to confirm any of those scenarios.