title: The mass fraction of halo stars contributed by the disruption of globular clusters in the E-MOSAICS simulations

abstract: Globular clusters (GCs) have been posited, alongside dwarf galaxies, as significant contributors to the field stellar population of the Galactic halo. In order to quantify their contribution, we examine the fraction of halo stars formed in stellar clusters in the suite of 25 present-day Milky Way-mass cosmological zoom simulations from the E-MOSAICS project. We find that a median of $2.3$ and $0.3$ per cent of the mass in halo field stars formed in clusters and GCs, defined as clusters more massive than $5\times 10^3$ and $10^5~M_{\odot}$, respectively, with the $25$-$75$th percentiles spanning $1.9$-$3.0$ and $0.2$-$0.5$ per cent being caused by differences in the assembly histories of the host galaxies. Under the extreme assumption that no stellar cluster survives to the present day, the mass fractions increase to a median of $5.9$ and $1.8$ per cent. These small fractions indicate that the disruption of GCs plays a sub-dominant role in the build-up of the stellar halo. We also determine the contributed halo mass fraction that would present signatures of light-element abundance variations considered to be unique to GCs, and find that clusters and GCs would contribute a median of $1.1$ and $0.2$ per cent, respectively. We estimate the contributed fraction of GC stars to the Milky Way halo, based on recent surveys, and find upper limits of $2$-$5$ per cent (significantly lower than previous estimates), suggesting that models other than those invoking strong mass loss are required to describe the formation of chemically enriched stellar populations in GCs.


title: Not all stars form in clusters -- $Gaia$-DR2 uncovers the origin of OB associations

abstract: Historically, it has often been asserted that most stars form in compact clusters. In this scenario, present-day gravitationally-unbound OB associations are the result of the expansion of initially gravitationally-bound star clusters. However, this paradigm is inconsistent with recent results, both theoretical and observational, that instead favour a hierarchical picture of star formation in which stars are formed across a continuous distribution of gas densities and most OB associations never were bound clusters. Instead they are formed in-situ as the low-density side of this distribution, rather than as the remnants of expanding clusters. We utilise the second $Gaia$ data release to quantify the degree to which OB associations are undergoing expansion and, therefore, whether OB associations are the product of expanding clusters, or whether they were born in-situ, as the large-scale, globally-unbound associations that we see today. We find that the observed kinematic properties of associations are consistent with highly substructured velocity fields and additionally require some degree of localised expansion from sub-clusters within the association. While most present-day OB associations do exhibit low levels of expansion, there is no significant correlation between radial velocity and radius. Therefore, the large-scale structure of associations is not set by the expansion of clusters, rather it is a relic of the molecular gas cloud from which the association was formed. This finding is inconsistent with a monolithic model of association formation and instead favours a hierarchical model, in which OB associations form in-situ, following the fractal structure of the gas from which they form.


title: The Stellar Population of Metal-Poor Galaxies at z $\approx$ 0.8 and the Evolution of the Mass-Metallicity Relation

abstract: We present results from deep Spitzer/Infrared Array Camera (IRAC) observations of 28 metal-poor, strongly star-forming galaxies selected from the DEEP2 Galaxy Survey. By modelling infrared and optical photometry, we derive stellar masses and other stellar properties. We determine that these metal-poor galaxies have low stellar masses, $M_{\star}$ $\approx10^{8.1}$-$10^{9.5}$ $M_{\odot}$. Combined with the Balmer-derived star formation rates (SFRs), these galaxies have average inverse SFR/$M_{\star}$ of $\approx$100 Myr. The evolution of stellar mass-gas metallicity relation to $z\approx0.8$ is measured by combining the modelled masses with previously obtained spectroscopic measurements of metallicity from [O III] $\lambda$4363 detections. Here, we include measurements for 79 galaxies from the Metal Abundances across Cosmic Time Survey. Our mass-metallicity relation is lower at a given stellar mass than at $z=0.1$ by 0.27 dex. This demonstrates a strong evolution in the mass-metallicity relation, $(1+z)^{-1.45^{+0.61}_{-0.76}}$. We find that the shape of the $z\approx0.8$ mass-metallicity relation, a steep rise in metallicity at low stellar masses, transitioning to a plateau at higher masses, is consistent with $z\sim0.1$ studies. We also compare the evolution in metallicity between $z\approx0.8$ and $z\sim0.1$ against recent strong-line diagnostic studies at intermediate redshifts and find good agreement. Specifically, we find that lower mass galaxies ($4\times10^8$ $M_{\odot}$) built up their metal content 1.6 times more rapidly than higher mass galaxies ($10^{10}$ $M_{\odot}$). Finally, we examine whether the mass-metallicity relation has a secondary dependence on SFR, and statistically concluded that there is no strong secondary dependence for $z\approx0.8$ low-mass galaxies.


title: Stellar winds pump the heart of the Milky Way

abstract: The central super-massive black hole of the Milky Way, Sgr A*, accretes at a very low rate making it a very underluminous galactic nucleus. Despite the tens of Wolf-Rayet stars present within the inner parsec supplying ${\sim}10^{-3}\rm\ M_{\odot}\ yr^{-1}$ in stellar winds, only a negligible fraction of this material ($<10^{-4}$) ends up being accreted onto Sgr A*. The recent discovery of cold gas (${\sim}10^4\rm\ K$) in its vicinity raised questions about how such material could settle in the hostile (${\sim}10^7\rm\ K$) environment near Sgr A*. In this work we show that the system of mass-losing stars blowing winds can naturally account for both the hot, inefficient accretion flow, as well as the formation of a cold disk-like structure. We run hydrodynamical simulations using the grid-based code Ramses starting as early in the past as possible to observe the state of the system at the present time. Our results show that the system reaches a quasi-steady state in about ${\sim}500\rm\ yr$ with material being captured at a rate of ${\sim}10^{-6}\rm\ M_{\odot}\ yr^{-1}$ at scales of ${\sim}10^{-4}\rm\ pc$, consistent with the observations and previous models. However, on longer timescales ($\gtrsim3000\rm\ yr$) the material accumulates close to the black hole in the form of a disk. Considering the duration of the Wolf-Rayet phase (${\sim}10^5\rm\ yr$), we conclude that this scenario likely has already happened, and could be responsible for the more active past of Sgr A*, and/or its current outflow. We argue that the hypothesis of the mass-losing stars being the main regulator of the activity of the black hole deserves further consideration.


title: Accessing the axion via compact object binaries

abstract: Black holes in binaries with other compact objects can provide natural venues for indirect detection of axions or other ultralight fields. The superradiant instability associated with a rapidly spinning black hole leads to the creation of an axion cloud which carries energy and angular momentum from the black hole. This cloud will then decay via gravitational wave emission. We show that the energy lost as a result of this process tends toward an outspiraling of the binary orbit. A given binary system is sensitive to a narrow range of axion masses, determined by the mass of the black hole. Pulsar-black hole binaries, once detected in the electromagnetic band, will allow high-precision measurements of black hole mass loss via timing measurements of the companion pulsar. This avenue of investigation is particularly promising in light of the recent preliminary announcements of two candidate black hole-neutron star mergers by LIGO/VIRGO (#S190814bv and #S190426c). We demonstrate that for such a binary system with a typical millisecond pulsar and a 3-solar-mass black hole, axions with masses between $2.7 \times 10^{-12}$ eV and $3.2 \times 10^{-12}$ eV are detectable. Recent gravitational wave observations by LIGO/VIRGO of binary black hole mergers imply that, for these binaries, gravitational radiation from the rotating quadrupole moment is a dominant effect, causing an inspiraling orbit. With some reasonable assumptions about the period of the binary when it formed and the spins of the black holes, these observations rule out possible axion masses between $3 \times 10^{-13}$ eV and $6 \times 10^{-13}$ eV. Future binary black hole observations, for example by LISA, are expected to provide more robust bounds. In some circumstances, neutron stars may also undergo superradiant instabilities, and binary pulsars could be used to exclude axions with certain masses and matter couplings.


title: New class of generalized coupling theories

abstract: We propose a new class of gravity theories which are characterized by a nontrivial coupling between the gravitational metric and matter mediated by an auxiliary rank-2 tensor. The actions generating the field equations are constructed so that these theories are equivalent to general relativity in a vacuum, and only differ from general relativity theory within a matter distribution. We analyze in detail one of the simplest realizations of these generalized coupling theories. We show that in this case the propagation speed of gravitational radiation in matter is different from its value in vacuum and that this can be used to weakly constrain the (single) additional parameter of the theory. An analysis of the evolution of homogeneous and isotropic spacetimes in the same framework shows that there exist cosmic histories with both an inflationary phase and a dark era characterized by a different expansion rate.


title: The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

abstract: We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates.


title: Meridional flows in the disk around a young star

abstract: Protoplanetary disks are known to posses a stunning variety of substructure in the distribution of their mm~sized grains, predominantly seen as rings and gaps (Andrews et al. 2018), which are frequently interpreted as due to the shepherding of large grains by either hidden, still-forming planets within the disk (Zhang et al. 2018) or (magneto-)hydrodynamic instabilities (Flock et al. 2015). The velocity structure of the gas offers a unique probe of both the underlying mechanisms driving the evolution of the disk, the presence of embedded planets and characterising the transportation of material within the disk, such as following planet-building material from volatile-rich regions to the chemically-inert midplane, or detailing the required removal of angular momentum. Here we present the radial profiles of the three velocity components of gas in upper disk layers in the disk of HD 163296 as traced by 12CO molecular emission. These velocities reveal significant flows from the disk surface towards the midplane of disk at the radial locations of gaps argued to be opened by embedded planets (Isella et al. 2016, 2018, Teague et al. 2018, Pinte et al. 2018), bearing striking resemblance to meridional flows, long predicted to occur during the early stages of planet formation (Szulagyi et al. 2014, Morbidelli et al. 2014, Fung & Chiang 2016, Dong et al. 2019). In addition, a persistent radial outflow is seen at the outer edge of the disk, potentially the base of a wind associated with previously detected extended emission (Klaassen et al. 2013).


title: Design and Characterization of Q-enhanced Silicon Nitride Racetrack Micro-Resonators

abstract: Q-enhanced racetrack micro-resonators for the silicon nitride photonics integration platform have been designed, fabricated and characterized. The proposed geometries permit to mitigate the impact of radiation loss at curved waveguides, one of the major limitations of silicon nitride circuits, therefore providing an increase of the intrinsic Q factor of micro-resonators when compared with the conventional structures with the same bent radii. The schemes put forward in this work permit a reduction of the size of the devices that has a direct impact on the integration scale in this platform. When used in the curved sections of waveguides routing optical signals within an integrated photonic circuit, these geometries provide a reduction of the radiation loss and permit the use of smaller bent radii and to increase the circuit integration density.


title: Automated Transient Detection with Shapelet Analysis in Image-subtracted Data

abstract: We present a method for characterizing image-subtracted objects based on shapelet analysis to identify transient events in ground-based time-domain surveys. We decompose the image-subtracted objects onto a set of discrete Zernike polynomials and use their resulting coefficients to compare them to other point-like objects. We derive a norm in this Zernike space that we use to score transients for their point-like nature and show that it is a powerful comparator for distinguishing image artifacts, or residuals, from true astrophysical transients. Our method allows for a fast and automated way of scanning overcrowded, wide-field telescope images with minimal human interaction and we reduce the large set of unresolved artifacts left unidentified in subtracted observational images. We evaluate the performance of our method using archival intermediate Palomar Transient Factory and Dark Energy Camera survey images. However, our technique allows flexible implementation for a variety of different instruments and data sets. This technique shows a reduction in image subtraction artifacts by 99.95% for surveys extending up to hundreds of square degrees and has strong potential for automated transient identification in electromagnetic follow-up programs triggered by the Laser Interferometer Gravitational Wave Observatory-Virgo Scientific Collaboration.


title: Electron Transport in Gaseous Detectors with a Python-based Monte Carlo Simulation Code

abstract: Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (`swarm') properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.


title: Half-mass radii of quiescent and star-forming galaxies evolve slowly from 0 < z < 2.5: implications for galaxy assembly histories

abstract: We use high-resolution, multi-band imaging of ~16,500 galaxies in the CANDELS fields at 0 < z < 2.5 to study the evolution of color gradients and half-mass radii over cosmic time. We find that galaxy color gradients at fixed mass evolve rapidly between z~2.5 and z~1, but remain roughly constant below z~1. This result implies that the sizes of both star-forming and quiescent galaxies increase much more slowly than previous studies found using half-light radii. The half-mass radius evolution of quiescent galaxies is fully consistent with a model which uses observed minor merger rates to predict the increase in sizes due to the accretion of small galaxies. Progenitor bias may still contribute to the growth of quiescent galaxies, particularly if we assume a slower timescale for the minor merger growth model. The slower half-mass radius evolution of star-forming galaxies is in tension with cosmological simulations and semi-analytic galaxy models. Further detailed, consistent comparisons with simulations are required to place these results in context.


title: How a minimal learning agent can infer the existence of unobserved variables in a complex environment

abstract: According to a mainstream position in contemporary cognitive science and philosophy, the use of abstract compositional concepts is both a necessary and a sufficient condition for the presence of genuine thought. In this article, we show how the ability to develop and utilise abstract conceptual structures can be achieved by a particular kind of learning agents. More specifically, we provide and motivate a concrete operational definition of what it means for these agents to be in possession of abstract concepts, before presenting an explicit example of a minimal architecture that supports this capability. We then proceed to demonstrate how the existence of abstract conceptual structures can be operationally useful in the process of employing previously acquired knowledge in the face of new experiences, thereby vindicating the natural conjecture that the cognitive functions of abstraction and generalisation are closely related. Keywords: concept formation, projective simulation, reinforcement learning, transparent artificial intelligence, theory formation, explainable artificial intelligence (XAI)


title: Auxiliary Space Preconditioning of Finite Element Equations Using a Nonconforming Interior Penalty Reformulation and Static Condensation

abstract: We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between neighboring elements. The resulting bilinear form, and a Schur complement (reduced) version of it, are utilized in a number of auxiliary space preconditioners for the original conforming finite element discretization problem. These preconditioners are analyzed on the fine scale and their performance is illustrated on model second order scalar elliptic problems discretized with high order elements.


title: The cosmic evolution of the stellar mass$-$velocity dispersion relation of early-type galaxies

abstract: We study the evolution of the observed correlation between central stellar velocity dispersion $\sigma_\mathrm{e}$ and stellar mass $M_*$ of massive ($M_*\gtrsim 3\times 10^{10}\,\mathrm{M_\odot}$) early-type galaxies (ETGs) out to redshift $z\approx 2.5$, exploiting a Bayesian hierarchical inference formalism. Collecting ETGs from state-of-the-art literature samples, we build a $fiducial$ sample ($0\lesssim z\lesssim 1$), which is obtained with homogeneous selection criteria, but also a less homogeneous $extended$ sample ($0\lesssim z\lesssim 2.5$). Based on the fiducial sample, we find that the $M_*$-$\sigma_\mathrm{e}$ relation is well represented by $\sigma_\mathrm{e}\propto M_*^{\beta}(1+z)^{\zeta}$, with $\beta\simeq 0.18$ independent of redshift and $\zeta\simeq 0.4$ (at given $M_*$, $\sigma_\mathrm{e}$ decreases for decreasing $z$, for instance by a factor of $\approx1.3$ from $z=1$ to $z=0$). When the slope $\beta$ is allowed to evolve, we find it increasing with redshift: $\beta(z)\simeq 0.16+0.26\log(1+z)$ describes the data as well as constant $\beta\simeq 0.18$. The intrinsic scatter of the $M_*$-$\sigma_\mathrm{e}$ relation is $\simeq0.08$ dex in $\sigma_\mathrm{e}$ at given $M_*$, independent of redshift. Our results suggest that, on average, the velocity dispersion of $individual$ massive ($M_*\gtrsim 3\times 10^{11}\,M_\odot$) ETGs decreases with time while they evolve from $z\approx 1$ to $z\approx 0$. The analysis of the extended sample leads to results similar to that of the fiducial sample, with slightly stronger redshift dependence of the normalisation ($\zeta\simeq 0.5$) and weaker redshift dependence of the slope (${\rm d} \beta/{\rm d} \log (1+z)\simeq 0.18$) when $\beta$ varies with time. At $z=2$ ETGs with $M_*\approx 10^{11}\,M_\odot$ have, on average, $\approx1.7$ higher $\sigma_\mathrm{e}$ than ETGs of similar stellar mass at $z=0$.


title: Autonomous Aerial Cinematography In Unstructured Environments With Learned Artistic Decision-Making

abstract: Aerial cinematography is revolutionizing industries that require live and dynamic camera viewpoints such as entertainment, sports, and security. However, safely piloting a drone while filming a moving target in the presence of obstacles is immensely taxing, often requiring multiple expert human operators. Hence, there is demand for an autonomous cinematographer that can reason about both geometry and scene context in real-time. Existing approaches do not address all aspects of this problem; they either require high-precision motion-capture systems or GPS tags to localize targets, rely on prior maps of the environment, plan for short time horizons, or only follow artistic guidelines specified before flight. In this work, we address the problem in its entirety and propose a complete system for real-time aerial cinematography that for the first time combines: (1) vision-based target estimation; (2) 3D signed-distance mapping for occlusion estimation; (3) efficient trajectory optimization for long time-horizon camera motion; and (4) learning-based artistic shot selection. We extensively evaluate our system both in simulation and in field experiments by filming dynamic targets moving through unstructured environments. Our results indicate that our system can operate reliably in the real world without restrictive assumptions. We also provide in-depth analysis and discussions for each module, with the hope that our design tradeoffs can generalize to other related applications. Videos of the complete system can be found at: https://youtu.be/ookhHnqmlaU.


title: Local and blowing-up solutions for an integro-differential diffusion equation and system

abstract: In the present paper initial problems for the semilinear integro-differential diffusion equation and system are considered. The analogue of Duhamel principle for the linear integro-differential diffusion equation is proved. The results on existence of local mild solutions and Fujita-type critical exponents to the semilinear integro-differential diffusion equation and system are presented.


title: The Renyi Gaussian Process: Towards Improved Generalization

abstract: We introduce an alternative closed form lower bound on the Gaussian process ($\mathcal{GP}$) likelihood based on the R\'enyi $\alpha$-divergence. This new lower bound can be viewed as a convex combination of the Nystr\"om approximation and the exact $\mathcal{GP}$. The key advantage of this bound, is its capability to control and tune the enforced regularization on the model and thus is a generalization of the traditional variational $\mathcal{GP}$ regression. From a theoretical perspective, we provide the convergence rate and risk bound for inference using our proposed approach. Experiments on real data show that the proposed algorithm may be able to deliver improvement over several $\mathcal{GP}$ inference methods.


title: Discussion of "The Blessings of Multiple Causes" by Wang and Blei

abstract: This commentary has two goals. We first critically review the deconfounder method and point out its advantages and limitations. We then briefly consider three possible ways to address some of the limitations of the deconfounder method.


title: A remark on an obstacle problem with lower regularity

abstract: We construct a monotone quantity for the classical obstacle problem with non-smooth obstacle, and show that the blow-ups are homogeneous functions of degree $\alpha<2$.


title: A note on the extremal noncentral sections of the cross-polytope

abstract: We find minimal and maximal length of intersections of lines at a fixed distance to the origin with the cross-polytope. We also find maximal volume noncentral sections of the cross-polypote by hyperplanes which are at a fixed large distance to the origin and minimal volume sections by symmetric slabs of a large fixed width. This parallels recent results about noncentral sections of the cube due to Moody, Stone, Zach and Zvavitch.


title: Gate field effects on the topological insulator BiSbTeSe2 interface

abstract: Interfaces between two topological insulators are of fundamental interest in condensed matter physics. Inspired by experimental efforts, we study interfacial processes between two slabs of BiSbTeSe2 (BSTS) via first principles calculations. Topological surface states are absent for the BSTS interface at its equilibrium separation, but our calculations show that they appear if the inter-slab distance is greater than 6 Ang. More importantly, we find that topological interface states can be preserved by inserting two or more layers of hexagonal boron nitride between the two BSTS slabs. In experiments, the electric current tunneling through the interface is insensitive to back gate voltage when the bias voltage is small. Using a first-principles based method that allows us to simulate gate field, we show that at low bias the extra charge induced by a gate voltage resides on the surface that is closest to the gate electrode, leaving the interface almost undoped. This provides clues to understand the origin of the observed insensitivity of transport properties to back voltage at low bias. Our study resolves a few questions raised in experiment, which does not yet offer a clear correlation between microscopic physics and transport data. We provide a road map for the design of vertical tunneling junctions involving the interface between two topological insulators.


title: Reduced-Order Modeling of Deep Neural Networks

abstract: We introduce a new method for speeding up the inference of deep neural networks. It is somewhat inspired by the reduced-order modeling techniques for dynamical systems.The cornerstone of the proposed method is the maximum volume algorithm. We demonstrate efficiency on neural networks pre-trained on different datasets. We show that in many practical cases it is possible to replace convolutional layers with much smaller fully-connected layers with a relatively small drop in accuracy.


title: Adaptive Exploration in Linear Contextual Bandit

abstract: Contextual bandits serve as a fundamental model for many sequential decision making tasks. The most popular theoretically justified approaches are based on the optimism principle. While these algorithms can be practical, they are known to be suboptimal asymptotically. On the other hand, existing asymptotically optimal algorithms for this problem do not exploit the linear structure in an optimal way and suffer from lower-order terms that dominate the regret in all practically interesting regimes. We start to bridge the gap by designing an algorithm that is asymptotically optimal and has good finite-time empirical performance. At the same time, we make connections to the recent literature on when exploration-free methods are effective. Indeed, if the distribution of contexts is well behaved, then our algorithm acts mostly greedily and enjoys sub-logarithmic regret. Furthermore, our approach is adaptive in the sense that it automatically detects the nice case. Numerical results demonstrate significant regret reductions by our method relative to several baselines.


title: Regularity for the two phase singular perturbation problems

abstract: We prove that an a priori BMO gradient estimate for the two phase singular perturbation problem implies Lipschitz regularity for the limits. This problem arises in the mathematical theory of combustion where the reaction-diffusion is modelled by the $p$-Laplacian. A key tool in our approach is the weak energy identity. Our method proves a natural and intrinsic characterization of the free boundary points and can be applied to more general classes of solutions.


title: A simple approach for vortex core visualization

abstract: We propose a method to visualize vortex cores based on manipulation of the pressure field produced by isolated vortices in incompressible flow. Under ideal conditions the function $D=2|\nabla P|/\nabla^2 P$ yields an approximate distance to vortex centerlines. As opposed to local methods to identify coherent structures, isosurfaces of $D$ produce a field of vortex tubes equidistant to the vortex core center which, ideally, are independent of vortex intensity or size. In contrast to other line-vortex identification methods, which typically rely on algorithms to detect vortex core lines and frequently need complex implementations, the proposed method can be computed from the local Eulerian velocity and pressure fields as easily as vortex identification methods such as the $Q$ and $\lambda_2$ criteria. $2|\nabla P|/\nabla^2 P$ results in the exact distance to the core center for a Rankine vortex and is in general valid for the region of a vortex where there is pure rotation, yielding an approximation to the distance farther from the core in other simple one-dimensional vortex models. The methodology performs well in all tests we attempted, though limitations are presented and discussed. The method is demonstrated for a canonical Burgers vortex, a Bodewadt vortex, homogeneous isotropic turbulent flow, the wake of a propeller, a heaving plate, a turning containership and the airwake of a surface combatant. The proposed method helps to better visualize vortical flow fields by displaying vortex cores, complementing methods like $Q$ and $\lambda_2$ which display vortical volumes.


title: High energy harmonic maps and degeneration of minimal surfaces

abstract: Let $S$ be a closed surface of genus $g \geq 2$ and let $\rho$ be a maximal $\mathrm{PSL}(2, \mathbb{R}) \times \mathrm{PSL}(2, \mathbb{R})$ surface group representation. By a result of Schoen, there is a unique $\rho$-equivariant minimal surface $\widetilde{\Sigma}$ in $\mathbb{H}^{2} \times \mathbb{H}^{2}$. We study the induced metrics on these minimal surfaces and prove the limits are precisely mixed structures. In the second half of the paper, we provide a geometric interpretation: the minimal surfaces $\widetilde{\Sigma}$ degenerate to the core of a product of two $\mathbb{R}$-trees. As a consequence, we obtain a compactification of the space of maximal representations of $\pi_{1}(S)$ into $\mathrm{PSL}(2, \mathbb{R}) \times \mathrm{PSL}(2, \mathbb{R})$.


title: Answering Complex Open-domain Questions Through Iterative Query Generation

abstract: It is challenging for current one-step retrieve-and-read question answering (QA) systems to answer questions like "Which novel by the author of 'Armada' will be adapted as a feature film by Steven Spielberg?" because the question seldom contains retrievable clues about the missing entity (here, the author). Answering such a question requires multi-hop reasoning where one must gather information about the missing entity (or facts) to proceed with further reasoning. We present GoldEn (Gold Entity) Retriever, which iterates between reading context and retrieving more supporting documents to answer open-domain multi-hop questions. Instead of using opaque and computationally expensive neural retrieval models, GoldEn Retriever generates natural language search queries given the question and available context, and leverages off-the-shelf information retrieval systems to query for missing entities. This allows GoldEn Retriever to scale up efficiently for open-domain multi-hop reasoning while maintaining interpretability. We evaluate GoldEn Retriever on the recently proposed open-domain multi-hop QA dataset, HotpotQA, and demonstrate that it outperforms the best previously published model despite not using pretrained language models such as BERT.


title: Quasi-one-dimensional approximation for Bose-Einstein condensates transversely trapped by a funnel potential

abstract: Starting from the standard three-dimensional (3D) Gross-Pitaevskii equation (GPE) and using a variational approximation, we derive an effective one-dimensional nonpolynomial Schr\"odinger equation (1D-NPSE) governing the axial dynamics of atomic Bose-Einstein condensates (BECs) under the action of a singular but physically relevant funnel-shaped transverse trap, i.e., an attractive 2D potential $\sim-1/r$ (where $r$ is the radial coordinate in the transverse plane), in combination with the repulsive self-interaction. Wave functions of the trapped BEC are regular, in spite of the potential's singularity. The model applies to a condensate of particles (small molecules) carrying a permanent electric dipole moment in the field of a uniformly charged axial thread, as well as to a quantum gas of magnetic atoms pulled by an axial electric current. By means of numerical simulations, we verify that the effective 1D-NPSE provides accurate static and dynamical results, in comparison to the full 3D GPE, for both repulsive and attractive signs of the intrinsic nonlinearity.


title: Trans-pars-planar illumination enables 200{\deg} ultra-wide field pediatric fundus photography for easy examination of the retina

abstract: This study is to test the feasibility of using trans-pars-planar illumination for ultrawide field pediatric fundus photography. Fundus examination of the peripheral retina is essential for clinical management of pediatric eye diseases. However, current pediatric fundus cameras with traditional trans-pupillary illumination provide a limited field of view (FOV), making it difficult to access the peripheral retina adequately for a comprehensive assessment of eye conditions. Here, we report the first demonstration of trans-pars-planar illumination in ultra-wide field pediatric fundus photography. For proof-of-concept validation, all off-the-shelf optical components were selected to construct a lab prototype pediatric camera (PedCam). By freeing the entire pupil for imaging purpose only, the trans-pars-planar illumination enables a 200o FOV in a snapshot fundus image, allowing easy visualization of both the central and peripheral retina up to the ora serrata. A low-cost, easy-to-use ultra-wide field PedCam provides a unique opportunity to foster affordable telemedicine in rural and underserved areas.


title: Constrained Bayesian Optimization with Max-Value Entropy Search

abstract: Bayesian optimization (BO) is a model-based approach to sequentially optimize expensive black-box functions, such as the validation error of a deep neural network with respect to its hyperparameters. In many real-world scenarios, the optimization is further subject to a priori unknown constraints. For example, training a deep network configuration may fail with an out-of-memory error when the model is too large. In this work, we focus on a general formulation of Gaussian process-based BO with continuous or binary constraints. We propose constrained Max-value Entropy Search (cMES), a novel information theoretic-based acquisition function implementing this formulation. We also revisit the validity of the factorized approximation adopted for rapid computation of the MES acquisition function, showing empirically that this leads to inaccurate results. On an extensive set of real-world constrained hyperparameter optimization problems we show that cMES compares favourably to prior work, while being simpler to implement and faster than other constrained extensions of Entropy Search.


title: The NAI Suite -- Drafting and Reasoning over Legal Texts

abstract: A prototype for automated reasoning over legal texts, called NAI, is presented. As an input, NAI accepts formalized logical representations of such legal texts that can be created and curated using an integrated annotation interface. The prototype supports automated reasoning over the given text representation and multiple quality assurance procedures. The pragmatics of the NAI suite as well its feasibility in practical applications is studied on a fragment of the Smoking Prohibition (Children in Motor Vehicles) (Scotland) Act 2016 of the Scottish Parliament.


title: Spontaneous Rotational Symmetry Breaking in a Kramers Two-Level System

abstract: Here, I develop a model for a two-level system that respects the time-reversal symmetry of the atom Hamiltonian and the Kramers theorem. The two-level system is formed by two Kramers pairs of excited and ground states. It is shown that due to the spin-orbit interaction it is in general impossible to find a basis of atomic states for which the crossed transition dipole moment vanishes. The parametric electric polarizability of the Kramers two-level system for a definite ground-state is generically nonreciprocal. I apply the developed formalism to study Casimir-Polder forces and torques when the two-level system is placed nearby either a reciprocal or a nonreciprocal substrate. In particular, I investigate the stable equilibrium orientation of the two-level system when both the atom and the reciprocal substrate have symmetry of revolution about some axis. Surprisingly, it is found that when chiral-type dipole transitions are dominant the stable ground state is not the one in which the symmetry axes of the atom and substrate are aligned. The reason is that the rotational symmetry may be spontaneously broken by the quantum vacuum fluctuations, so that the ground state has less symmetry than the system itself.


title: Three-Nucleon Force Effects in the FSI configuration of the d(n, nn)p Breakup Reaction

abstract: We investigated three-nucleon (3N) force effects in the final state interaction (FSI) configuration of the d(n, nn)p breakup reaction at the incoming nucleon energy E n = 200 MeV. Although 3N force effects for the elastic nucleon-deuteron scattering cross section at comparable energies are located predominantly in the region of intermediate and backward angles, the corresponding 3N force effects for the integrated FSI configuration breakup cross section are found also at forward scattering angles.


title: Kinematics: A Clean Diagnostic for Separating Supernova Remnants from HII Regions in Nearby Galaxies

abstract: Many more supernova remnants (SNRs) are now known in external galaxies than in the Milky Way. Most of these SNRs have been identified using narrow-band imaging, separating SNRs from HII regions on the basis of [SII]:H-alpha ratios that are elevated compared to HII regions. However, the boundary between SNRs and HII regions is not always distinct, especially at low surface brightness. Here we explore velocity structure as a possible criterion for separating SNRs from HII regions, using a sample of well-studied SNRs in the Large Magellanic Cloud (LMC) as well as a small number of SNRs in the galaxy M83. We find, perhaps not surprisingly, that even at large diameters, SNRs exhibit velocity broadening sufficient to readily distinguish them from HII regions. We thus suggest that the purity of most extragalactic samples would be greatly improved through spectroscopic observations with a velocity resolution of order 50~km/s$.


title: Erbium-doped aluminophosphosilicate all-fiber laser operating at 1584 nm

abstract: We report on an ytterbium-free erbium-doped aluminophosphosilicate all-fiber laser, producing an output power of 25 W at a wavelength of 1584 nm with a slope efficiency of 30% with respect to the 976 nm absorbed pump power. The simple cavity design proposed takes advantage of fiber Bragg gratings written directly in the gain fiber. The single-mode erbium-doped aluminophosphosilicate fiber was fabricated in-house and was doped with 0.06 mol.% of Er2O3, 1.77 mol.% of Al2O3 and 1.04 mol.% of P2O5. The incorporation of aluminium and phosphorus into the fiber core allowed for an increased concentration of erbium without inducing significant clustering, while keeping the numerical aperture low to ensure a single-mode laser operation.


title: Spectral analysis of Fermi-LAT gamma-ray bursts with known redshift and their potential use as cosmological standard candles

abstract: Long duration Gamma-Ray Bursts (LGRBs) may serve as standard candles to constrain cosmological parameters by probing the Hubble diagram well beyond the range of redshift currently accessible using type-Ia supernovae. The standardization of GRBs is based on phenomenological relations between two or more parameters found from spectral modeling, of which one is strongly dependent on the cosmological model. The Amati relation links the source-frame energy ${E_{\mathrm{i,p}}}$ at which the prompt gamma-ray spectral energy distribution $\nu F_\nu$ peaks, and the isotropic-equivalent bolometric energy ${E_{\mathrm{iso}}}$ emitted during the prompt phase. We performed spectral analysis of 26 GRBs with known redshift that have been detected by the Fermi-Large Area Telescope (LAT) during its nine years of operations from July 2008 to September 2017, thus extending the computation of ${E_{\mathrm{iso}}}$ to the 100 MeV range. Multiple components are required to fit the spectra of a number of GRBs. We found that the Amati relation is satisfied by the 25 LGRBs, with best fit parameters similar to previous studies that used data from different satellite experiments, while the only short GRB with known redshift is an outlier. Using the Amati relation we extend the Hubble diagram to redshift 4.35 and constrain the Hubble constant and dark-energy density in the $\Lambda$CDM model, with Fermi-LAT GRBs alone and together with another sample of 94 GRBs and with the latest Supernovae type-Ia data. Our results are consistent with the currently acceptable ranges of those cosmological parameters within errors.


title: Superconducting nanowire single-photon detectors fabricated from atomic-layer-deposited NbN

abstract: High-quality ultra-thin films of niobium nitride (NbN) are developed by plasma-enhanced atomic layer deposition (PEALD) technique. Superconducting nanowire single-photon detectors (SNSPDs) patterned from this material exhibit high switching currents and saturated internal efficiencies over a broad bias range at 1550 nm telecommunication wavelength. Statistical analyses on hundreds of fabricated devices show near-unity throughput yield due to exceptional homogeneity of the films. The ALD-NbN material represents an ideal superconducting material for fabricating large single-photon detector arrays combining high efficiency, low jitter, low dark counts.


title: Parameter-free velocity-dependent one-scale model for domain walls

abstract: We develop a parameter-free velocity-dependent one-scale model for the evolution of the characteristic length $L$ and root-mean-square velocity $\sigma_v$ of standard domain wall networks in homogeneous and isotropic cosmologies. We compare the frictionless scaling solutions predicted by our model, in the context of cosmological models having a power law evolution of the scale factor $a$ as a function of the cosmic time $t$ ($a \propto t^\lambda$, $0< \lambda < 1$), with the corresponding results obtained using field theory numerical simulations. We show that they agree well (within a few $\%$) for root-mean-square velocities $\sigma_v$ smaller than $0.2 \, c$ ($\lambda \ge 0.9$), where $c$ is the speed of light in vacuum, but significant discrepancies occur for larger values of $\sigma_v$ (smaller values of $\lambda$). We identify problems with the determination of $L$ and $\sigma_v$ from numerical field theory simulations which might potentially be responsible for these discrepancies.


title: Transfer Learning for Algorithm Recommendation

abstract: Meta-Learning is a subarea of Machine Learning that aims to take advantage of prior knowledge to learn faster and with fewer data [1]. There are different scenarios where meta-learning can be applied, and one of the most common is algorithm recommendation, where previous experience on applying machine learning algorithms for several datasets can be used to learn which algorithm, from a set of options, would be more suitable for a new dataset [2]. Perhaps the most popular form of meta-learning is transfer learning, which consists of transferring knowledge acquired by a machine learning algorithm in a previous learning task to increase its performance faster in another and similar task [3]. Transfer Learning has been widely applied in a variety of complex tasks such as image classification, machine translation and, speech recognition, achieving remarkable results [4,5,6,7,8]. Although transfer learning is very used in traditional or base-learning, it is still unknown if it is useful in a meta-learning setup. For that purpose, in this paper, we investigate the effects of transferring knowledge in the meta-level instead of base-level. Thus, we train a neural network on meta-datasets related to algorithm recommendation, and then using transfer learning, we reuse the knowledge learned by the neural network in other similar datasets from the same domain, to verify how transferable is the acquired meta-knowledge.


title: Transients following white dwarfs mergers

abstract: Mergers of white dwarfs (WDs) may lead to a variety of transient astrophysical events, SNIa being one possible outcome. Lyutikov & Toonen (2017, 2019) argued that mergers of WDs result, under various parameter regimes, in unusual central engine-powered supernova and a type of short Gamma Ray Bursts that show extended emission tails. Observations by Gvaramadze et al. (2019) of the central star and the nebula J005311 match to the details the model of Lyutikov & Toonen (2017, 2019) for the immediate product of a merger of a heavy ONeMg WD with CO WD (age, luminosity, stellar size, hydrogen deficiency and chemical composition).


title: Subsistence of ice-covered lakes during the Hesperian at Gale crater, Mars

abstract: Sedimentary deposits characterized by the Mars Science Laboratory Curiosity rover provide evidence that Gale crater, Mars intermittently hosted a fluvio-lacustrine environment during the Hesperian. However, estimates of the CO2 content of the atmosphere at the time the sediments in Gale crater were deposited are far less than needed by any climate model to maintain temperatures warm enough for sustained open water lake conditions due to the low solar energy input available at that time. We incrementally test the existence of open water conditions using energy balance calculations for the global, regional, and seasonal temperatures, and we assess if the preservation of liquid water was possible under perennial ice covers. We found scenarios where lacustrine conditions are preserved in a cold climate, where the resupply of water by the inflow of rivers and high precipitation rates are substituted by an abutting glacier. For equatorial temperatures as low as 240K-255K, the ice thickness ranges from 3-10 m, a value comparable to the range of those for the perennially ice-covered lakes in Antarctica (3-6 m). The ice-covered lake hypothesis is a compelling way to decouple the mineralogy and the climate by limiting the gas exchanges between the sediment and the CO2 atmosphere, and it eliminates the requirement for global mean temperatures above the freezing point. Not only do ice-covered lakes provide a baseline for exploring the range of possible lake scenarios for Gale crater during the Hesperian that is fully consistent with climate studies, but also they might have been ideal environments to sustain life on Mars.


title: Dynamically Aggregating Diverse Information

abstract: An agent has access to multiple information sources, each of which provides information about a different attribute of an unknown state. Information is acquired continuously -- where the agent chooses both which sources to sample from, and also how to allocate attention across them -- until an endogenously chosen time, at which point a decision is taken. We provide an exact characterization of the optimal information acquisition strategy under weak conditions on the agent's prior belief about the different attributes. We then apply this characterization to derive new results regarding: (1) endogenous information acquisition for binary choice, (2) strategic information provision by biased news sources, and (3) the dynamic consequences of attention manipulation.


title: On Harmonic Approximations of Inharmonic Signals

abstract: In this work, we present the misspecified Gaussian Cram\'er-Rao lower bound for the parameters of a harmonic signal, or pitch, when signal measurements are collected from an almost, but not quite, harmonic model. For the asymptotic case of large sample sizes, we present a closed-form expression for the bound corresponding to the pseduo-true fundamental frequency. Using simulation studies, it is shown that the bound is sharp and is attained by maximum likelihood estimators derived under the misspecified harmonic assumption. It is shown that misspecified harmonic models achieve a lower mean squared error than correctly specified unstructured models for moderately inharmonic signals. Examining voices from a speech database, we conclude that human speech belongs to this class of signals, verifying that the use of a harmonic model for voiced speech is preferable.


title: Bayesian variable selection in hierarchical difference-in-differences models

abstract: A popular method for estimating a causal treatment effect with observational data is the difference-in-differences (DiD) model. In this work, we consider an extension of the classical DiD setting to the hierarchical context in which data cannot be matched at the most granular level (e.g., individual-level differences are unobservable). We propose a Bayesian hierarchical difference-in-differences (HDiD) model which estimates the treatment effect by regressing the treatment on a latent variable representing the mean change in group-level outcome. We present theoretical and empirical results showing that an HDiD model that fails to adjust for a particular class of confounding variables, or confounding with the baseline (pre-treatment) outcomes, biases the treatment effect estimate. We propose and implement various approaches to perform variable selection using a structured Bayesian spike-and-slab model in the HDiD context. Our proposed methods leverage the temporal structure within the DiD context to select those covariates that lead to unbiased and efficient estimation of the causal treatment effect. We evaluate the methods' properties through theoretical results and simulation, and we use them to assess the impact of primary care redesign of clinics in Minnesota on the management of diabetes outcomes from 2008 to 2017.


title: Games of Incomplete Information Played By Statisticians

abstract: Players are statistical learners who learn about payoffs from data. They may interpret the same data differently, but have common knowledge of a class of learning procedures. I propose a metric for the analyst's "confidence" in a strategic prediction, based on the probability that the prediction is consistent with the realized data. The main results characterize the analyst's confidence in a given prediction as the quantity of data grows large, and provide bounds for small datasets. The approach generates new predictions, e.g. that speculative trade is more likely given high-dimensional data, and that coordination is less likely given noisy data.


title: Subluminal cosmological bounce beyond Horndeski

abstract: We address the issue of potential superluminal propagation of gravitational waves in backgrounds neighboring the previously suggested bounce [arXiv:1807.08361] in beyond Horndeski theory. We find that the bouncing solution lies right at the boundary of the region where the gravitational waves propagate at speed exceeding that of light, i.e. that solution suffers superluminality problem. We suggest a novel version of a completely stable bouncing model where both scalar and tensor perturbations remain safely subluminal not only on the solution itself but also in its neighbourhood. The model remains free of superluminality when extra matter in the form of radiation or, more generally, ideal fluid with equation of state parameter $w\leq 1/3$ (and also somewhat higher) is added. Superluminality reappears when extra matter is added whose sound velocity is equal or close to 1 in flat space; an example is scalar field minimally coupled to metric. The latter property is characteristic of all beyond Horndeski cosmologies; we briefly discuss its significance.


title: Privacy Preserving Count Statistics

abstract: The ability to preserve user privacy and anonymity is important. One of the safest ways to maintain privacy is to avoid storing personally identifiable information (PII), which poses a challenge for maintaining useful user statistics. Probabilistic counting has been used to find the cardinality of a multiset when precise counting is too resource intensive. In this paper, probabilistic counting is used as an anonymization technique that provides a reliable estimate of the number of unique users. We extend previous work in probabilistic counting by considering its use for preserving user anonymity, developing application guidelines and including hash collisions in the estimate. Our work complements previous method by attempting to explore the causes of the deviation of uncorrected estimate from the real value. The experimental results show that if the proper register size is used, collision compensation provides estimates are as good as, if not better than, the original probabilistic counting. We develop a new anonymity metric to precisely quantify the degree of anonymity the algorithm provides.


title: Free Division Rings of Fractions of Crossed Products of Groups With Conradian Left-Orders

abstract: Let $D$ be a division ring of fractions of a crossed product $F[G,\eta,\alpha]$ where $F$ is a skew field and $G$ is a group with Conradian left-order $\leq$. For $D$ we introduce the notion of freeness with respect to $\leq$ and show that $D$ is free in this sense if and only if $D$ can canonically be embedded into the endomorphism ring of the right $F$-vector space $F((G))$ of all formal power series in $G$ over $F$ with respect to $\leq$. From this we obtain that all division rings of fractions of $F[G,\eta,\alpha]$ which are free with respect to at least one Conradian left-order of $G$ are isomorphic and that they are free with respect to any Conradian left-order of $G$. Moreover, $F[G,\eta,\alpha]$ possesses a division ring of fraction which is free in this sense if and only if the rational closure of $F[G,\eta,\alpha]$ in the endomorphism ring of the corresponding right $F$-vector space $F((G))$ is a skew field.


title: Measuring the Completeness of Theories

abstract: We use machine learning to provide a tractable measure of the amount of predictable variation in the data that a theory captures, which we call its "completeness." We apply this measure to three problems: assigning certain equivalents to lotteries, initial play in games, and human generation of random sequences. We discover considerable variation in the completeness of existing models, which sheds light on whether to focus on developing better models with the same features or instead to look for new features that will improve predictions. We also illustrate how and why completeness varies with the experiments considered, which highlights the role played in choosing which experiments to run.


title: Quadruplex detection in human cells

abstract: DNA and RNA quadruplexes are extensively studied for the key roles they are suspected to play in the cellular regulation networks at both genomic and transcriptomic levels. The reliable detection of quadruplexes in cells was and remains a challenging task. Here, we describe the various strategies that have been implemented over the past years to visualize functionally relevant DNA and RNA quadruplexes in human cells, from immunodetection studies to the design and use of quadruplex-specific turn-on fluorescent probes.


title: Slingshot droplet generator for clogging-free and tunable production of small uniform drops in air

abstract: Submillimeter air-borne drops of well defined size are notoriously difficult to produce in the lab. Here, we demonstrate that any simple droplet dripper, in which millimetric droplets are created by letting them drip under gravity from a hollow needle, is readily turned into a versatile droplet generator capable of producing drops smaller than the needle diameter. This is achieved by vibrating the needle base at a frequency close to a flexing resonance of the needle tube. Through high-speed imaging and detailed image analysis we show that the droplets flung off one-by-one from the needle tip are uniform in size. The produced drop size is mainly set by the driving frequency and is practically independent of needle geometry or driving amplitude, making the generator robust and easily tunable. Like for a simple dripper, the drop ejection rate is independently controlled through the flow-rate. Since the needle diameter can be much larger than the produced drop size, clogging is not an issue. By using needles with different resonant frequencies we could produce droplets with radii anywhere between 100 microns to 1 mm with a high reproducibility (the production of smaller drops is possible in principle). We propose a mass-spring model for the vibrated pendant drop that captures all our observations.


title: The Sarkisov program on log surfaces

abstract: We show that the Sarkisov program holds for $\mathbb{Q}$-factorial log surfaces and log canonical surfaces over any algebraically closed field.


title: Robustness of the Kondo effect in a quantum dot coupled to Majorana zero modes

abstract: The prospect of using semiconductor quantum dots as an experimental tool to distinguish Majorana zero modes (MZMs) from other zero-energy excitations such as Kondo resonances has brought up the fundamental question of whether topological superconductivity and the Kondo effect can coexist in these systems. Here, we study the Kondo effect in a quantum dot coupled to a metallic contact and to a pair of MZMs. We consider a situation in which the MBS are spin polarized in opposite directions. By using numerical renormalization-group calculations and scaling analysis of the renormalization group equations, we show that the Kondo effect takes place at low temperatures, regardless the coupling to the MZMs. Interestingly, we find that the Kondo singlet essentially decouples from the MZMs such that the residual impurity entropy can show local non-Fermi liquid properties characteristic of the single Majorana excitations. This offers the possibility of tuning between Fermi-liquid and non-Fermi-liquid regimes simply by changing the quantum dot-MZM couplings.


title: Jet quenching in the hadron gas: an exploratory study

abstract: The suppression of high momentum particles in heavy-ion collisions in comparison to elementary reactions is one of the main indications for the formation of a quark-gluon plasma. In recent studies, full jets are being reconstructed and substructure observables are gaining importance in assessing the medium modifications of hard probes. In this work, the effect of the late stage hadronic interactions are explored within the hadronic transport approach SMASH (Simulating Many Accelerated Strongly-interacting Hadrons). High momentum particles are incorporated in a radially expanding hadron gas to analyse the corresponding angular distributions, also refered to as `jet shape' observables. We find that the full hadron gas can be approximated with a pion gas with constant elastic cross-sections of 100 mb. In addition, the temperature and probe energy dependence of diffusion coefficients $\tilde{q}$ and $\tilde{e}$ quantifying the transverse and parallel momentum transfers are extracted. The species dependence and the importance of different interaction types are investigated. Parametrizations are presented that can be employed in future jet quenching calculations to include the effect of the hadronic phase.


title: Techno-economic analysis of solar PV power-to-heat-to-power storage and trigeneration in the residential sector

abstract: This article assesses whether it is profitable to store solar PV electricity in the form of heat and convert it back to electricity on demand. The impact of a number of technical and economic parameters on the profitability of a self-consumption residential system located in Madrid is assessed. The proposed solution comprises two kinds of heat stores: a low- or medium-grade heat store for domestic hot water and space heating, and a high-grade heat store for combined heat and power generation. Two cases are considered where the energy that is wasted during the conversion of heat into electricity is employed to satisfy either the heating demand, or both heating and cooling demands by using a thermally-driven heat pump. We compare these solutions against a reference case that relies on the consumption of grid electricity and natural gas and uses an electrically-driven heat pump for cooling. The results show that, under relatively favourable conditions, the proposed solution that uses an electrically-driven heat pump could provide electricity savings in the range of 70-90% with a payback period of 12-15 years, plus an additional 10-20% reduction in the fuel consumption. Shorter payback periods, lower than 10 years, could be attained by using a highly efficient thermally driven heat pump, at the expense of increasing the fuel consumption and the greenhouse gas emissions. Hybridising this solution with solar thermal heating could enable significant savings on the global emissions, whilst keeping a high amount of savings in grid electricity (> 70 %) and a reasonably short payback period (< 12 years).


title: End-to-end particle and event identification at the Large Hadron Collider with CMS Open Data

abstract: From particle identification to the discovery of the Higgs boson, deep learning algorithms have become an increasingly important tool for data analysis at the Large Hadron Collider (LHC). We present an innovative end-to-end deep learning approach for jet identification at the Compact Muon Solenoid (CMS) experiment at the LHC. The method combines deep neural networks with low-level detector information, such as calorimeter energy deposits and tracking information, to build a discriminator to identify different particle species. Using two physics examples as references: electron vs. photon discrimination and quark vs. gluon discrimination, we demonstrate the performance of the end-to-end approach on simulated events with full detector geometry as available in the CMS Open Data. We also offer insights into the importance of the information extracted from various sub-detectors and describe how end-to-end techniques can be extended to event-level classification using information from the whole CMS detector.


title: SGD Learns One-Layer Networks in WGANs

abstract: Generative adversarial networks (GANs) are a widely used framework for learning generative models. Wasserstein GANs (WGANs), one of the most successful variants of GANs, require solving a minmax optimization problem to global optimality, but are in practice successfully trained using stochastic gradient descent-ascent. In this paper, we show that, when the generator is a one-layer network, stochastic gradient descent-ascent converges to a global solution with polynomial time and sample complexity.


title: Deep learning predictions of galaxy merger stage and the importance of observational realism

abstract: Machine learning is becoming a popular tool to quantify galaxy morphologies and identify mergers. However, this technique relies on using an appropriate set of training data to be successful. By combining hydrodynamical simulations, synthetic observations and convolutional neural networks (CNNs), we quantitatively assess how realistic simulated galaxy images must be in order to reliably classify mergers. Specifically, we compare the performance of CNNs trained with two types of galaxy images, stellar maps and dust-inclusive radiatively transferred images, each with three levels of observational realism: (1) no observational effects (idealized images), (2) realistic sky and point spread function (semi-realistic images), (3) insertion into a real sky image (fully realistic images). We find that networks trained on either idealized or semi-real images have poor performance when applied to survey-realistic images. In contrast, networks trained on fully realistic images achieve 87.1% classification performance. Importantly, the level of realism in the training images is much more important than whether the images included radiative transfer, or simply used the stellar maps (87.1% compared to 79.6% accuracy, respectively). Therefore, one can avoid the large computational and storage cost of running radiative transfer with a relatively modest compromise in classification performance. Making photometry-based networks insensitive to colour incurs a very mild penalty to performance with survey-realistic data (86.0% with r-only compared to 87.1% with gri). This result demonstrates that while colour can be exploited by colour-sensitive networks, it is not necessary to achieve high accuracy and so can be avoided if desired. We provide the public release of our statistical observational realism suite, RealSim, as a companion to this paper.


title: Newton transformations and motivic invariants at infinity of plane curves

abstract: In this article we give an expression of the motivic Milnor fiber at infinity and the motivic nearby cycles at infinity of a polynomial $f$ in two variables with coefficients in an algebraic closed field of characteristic zero. This expression is given in terms of some motives associated to the faces of the Newton polygons appearing in the Newton algorithm at infinity of $f$ without any condition of convenience or non degeneracy. In the complex setting, we compute the Euler characteristic of the generic fiber of $f$ in terms of the area of the surfaces associated to faces of the Newton polygons. Furthermore, if $f$ has isolated singularities, we compute similarly the classical invariants at infinity $\lambda_{c}(f)$ which measures the non equisingularity at infinity of the fibers of $f$ in $\mathbb P^2$, and we prove the equality between the topological and the motivic bifurcation sets and give an algorithm to compute them.


title: Braiding Majorana Fermions and Creating Quantum Logic Gates with Vortices on a Periodic Pinning Structure

abstract: We show how vortices that support Majorana fermions when placed on a periodic pinning array can be used for vortex exchange and independent braiding by performing a series of specific moves with a probe tip. Using these braiding operations, we demonstrate realizations of a Hadamard and a CNOT gate. We specifically consider the first matching field at which there is one vortex per pinning site, and we show that there are two basic dynamic operations, move and exchange, from which basic braiding operations can be constructed in order to create specific logic gates. The periodic pinning array permits both control of the world lines of the vortices and freedom for vortex manipulation using a set of specific moves of the probe during which the probe tip strength and height remain unchanged. We measure the robustness of the different moves against thermal effects and show that the three different operations produce distinct force signatures on the moving tip.


title: A class of prime fusion categories of dimension $2^N$

abstract: We study a class of strictly weakly integral fusion categories $\mathfrak{I}_{N, \zeta}$, where $N \geq 1$ is a natural number and $\zeta$ is a $2^N$th root of unity, that we call $N$-Ising fusion categories. An $N$-Ising fusion category has Frobenius-Perron dimension $2^{N+1}$ and is a graded extension of a pointed fusion category of rank 2 by the cyclic group of order $\mathbb Z_{2^N}$. We show that every braided $N$-Ising fusion category is prime and also that there exists a slightly degenerate $N$-Ising braided fusion category for all $N > 2$. We also prove a structure result for braided extensions of a rank 2 pointed fusion category in terms of braided $N$-Ising fusion categories.


title: Cherenkov Light in Liquid Scintillator at the NOvA Experiment

abstract: NOvA is a long-baseline neutrino oscillation experiment with two functionally identical liquid scintillator tracking detectors, i.e. the Near Detector (ND) and the Far Detector (FD). One of NOvA's physics goal is to measure neutrino oscillation parameters by studying $\nu_e$ appearance and $\nu_{\mu}$ disappearance with the Neutrinos at the Main Injector (NuMI) beam at Fermi National Accelerator Laboratory. An accurate light model is a prerequisite for precise charged particle energy estimation in the detectors. Particle energy is needed in event reconstruction and classification, both of which are critical to constraining neutrino oscillation parameters. In this paper, I will explain the details of the data-driven tuning of the Cherenkov model and the impact of the new scintillator model


title: Representation stability for sequences of 0-Hecke modules

abstract: We define a new category analogous to ${\bf FI}$ for the $0$-Hecke algebra $H_n(0)$ called the $0$-Hecke category, $\mathcal{H}$, indexing sequences of representations of $H_n(0)$ as $n$ varies under suitable compatibility conditions. We establish a new type of representation stability in this setting and prove it is implied by being a finitely generated $\mathcal{H}$-module. We then provide examples of $\mathcal{H}$-modules and discuss further desirable properties these modules possess.


title: Hybrid Thermal Yagi-Uda Nanoantennas for Directional and Narrow Band Long-wavelength IR Radiation Sources

abstract: We investigate the possibility of spatially and spectrally controlling the thermal infrared emission by exploitation of the Yagi-Uda antenna design. Hybrid antennas composed of both SiC and Au rods are considered and the contributions of emission from all the elements, at a given equilibrium temperature, are taken into account. We show that the detrimental effect due to thermal emission from the not ideal parasitic elements drastically affect the performances of conventional thermal Au antennas in the 12 micrometers wavelength range. Nevertheless, our results show that the hybrid approach allows the development of efficient narrow-band and high directivity sources. The possibility of exploiting the Yagi-Uda design both in transmission and in reception modes, may open the way to the realization of miniaturized, efficient, robust and cheap sensor devices for mass-market applications.


title: Compact Network Training for Person ReID

abstract: The task of person re-identification (ReID) has attracted growing attention in recent years leading to improved performance, albeit with little focus on real-world applications. Most SotA methods are based on heavy pre-trained models, e.g. ResNet50 (~25M parameters), which makes them less practical and more tedious to explore architecture modifications. In this study, we focus on a small-sized randomly initialized model that enables us to easily introduce architecture and training modifications suitable for person ReID. The outcomes of our study are a compact network and a fitting training regime. We show the robustness of the network by outperforming the SotA on both Market1501 and DukeMTMC. Furthermore, we show the representation power of our ReID network via SotA results on a different task of multi-object tracking.


title: Hierarchical Interactive Architecture based on Coalition Formation for Neighborhood System Decision Making

abstract: In recent years, significant efforts for improving the technical and the economic performance of smart grids have been applied with the presence of different players making decisions in these grids. This paper proposes a bi-level hierarchical structure for designing and planning distributed energy resources in home microgrids (H-MGs). Such a small-scale grid inside the market environment permits energy exchange among distributed energy resources and home microgrids through a pool market. In this paper, each home microgrid's planner maximizes the performance of distributed energy resources while exchanging energy with other H-MGs. The problem is formulated as a high-level problem decomposed into a set of low-level market clearing problems. The global optimal performance in terms of energy cost is met for a market structure (H-MGs, consumers, and retailers) at a Nash optimum point for the formulated scheduling game considering local and general constraints on the spot market. In general, the upper-level structure is based on H-MG generation competition for maximizing their individual and/or group income in the process of forming coalition with other H-MGs. In addition, the functionality of the lower-level of the hierarchical structure is governed by a market clearing based on the price response method by all the DERs enabling H-MGs to change the spot market strategic behavior. Prices are obtained as dual variables of power balance equations. This paper will investigate a set of bilevel games, including nonlinear programming problems solved through a complementary problem method. Using the binary theory, the bilevel hierarchical structure will be replanned using a nonlinear problem. Results prove that the proposed structure can increase the players's profit.


title: Remark on the synergy between the heat kernel techniques and the parity anomaly

abstract: In this paper, we demonstrate that not only the heat kernel techniques are useful for computation of the parity anomaly, but also the parity anomaly turns out to be a powerful mean in studying the heat kernel. We show that the gravitational parity anomaly on 4D manifolds with boundaries can be calculated using the general structure of the heat kernel coefficient $a_5$ for mixed boundary conditions, keeping all the weights of various geometric invariants as unknown numbers. The symmetry properties of the $\eta$-invariant allow to fix all the relevant unknowns. As a byproduct of this calculation we get an efficient and independent crosscheck (and confirmation) of the correction of the general structure of $a_5$ for mixed boundary conditions, previously suggested in Ref. [59].


title: Methods for Preparing Quantum Gases of Lithium

abstract: Lithium is an important element in atomic quantum gas experiments because its interactions are highly tunable, due to broad Feshbach resonances and zero-crossings, and because it has two stable isotopes, $^6$Li, a fermion, and $^7$Li, a boson. Although lithium has special value for these reasons, it also presents experimental challenges. In this article, we review some of the methods that have been developed or adapted to confront these challenges, including beam and vapor sources, Zeeman slowers, sub-Doppler laser cooling, laser sources at 671 nm, and all-optical methods for trapping and cooling. Additionally, we provide spectral diagrams of both $^6$Li and $^7$Li, and present plots of Feshbach resonances for both isotopes.


title: MUTE: Data-Similarity Driven Multi-hot Target Encoding for Neural Network Design

abstract: Target encoding is an effective technique to deliver better performance for conventional machine learning methods, and recently, for deep neural networks as well. However, the existing target encoding approaches require significant increase in the learning capacity, thus demand higher computation power and more training data. In this paper, we present a novel and efficient target encoding scheme, MUTE to improve both generalizability and robustness of a target model by understanding the inter-class characteristics of a target dataset. By extracting the confusion level between the target classes in a dataset, MUTE strategically optimizes the Hamming distances among target encoding. Such optimized target encoding offers higher classification strength for neural network models with negligible computation overhead and without increasing the model size. When MUTE is applied to the popular image classification networks and datasets, our experimental results show that MUTE offers better generalization and defense against the noises and adversarial attacks over the existing solutions.


title: A new fitting concept for the robust determination of S\'ersic model parameters

abstract: The S\'ersic law (SL) offers a versatile functional form for the structural characterization of galaxies near and far. Whereas applying it to galaxies with a genuine SL luminosity distribution yields a robust determination of the S\'ersic exponent eta and effective surface brightness $\mu_{\rm eff}$, this is not necessarily the case for galaxies whose surface brightness profiles (SBPs) appreciably deviate from the SL (eg, early-type galaxies with a depleted core and nucleated dwarf ellipticals, or most late-type galaxies-LTGs). In this general case of "imperfect" SL profiles, the best-fitting solution may significantly depend on the radius (or surface brightness) interval fit and corrections for point spread function (PSF) convolution effects. Such uncertainties may then affect, in a non-easily predictable manner, automated structural studies of galaxies. We present a fitting concept (iFIT) that permits a robust determination of the equivalent SL model for the general case of galaxies with imperfect SL profiles. iFIT has been extensively tested on synthetic data with a S\'ersic index 0.3<${\eta}$<4.2 and an effective radius 1<$\rm{R}_{eff}$ (arcs)<20. Applied to non PSF-convolved data, iFIT can infer the S\'ersic exponent eta with an absolute error of <0.2 even for shallow SBPs. As for PSF-degraded data, iFIT can recover the input SL model parameters with a satisfactorily accuracy almost over the entire considered parameter space as long as FWHM(PSF)<$\rm{R}_{eff}$. Tests indicate that iFIT shows little sensitivity on PSF corrections and the SBP limiting surface brightness, and that subtraction of the best-fitting SL model in two different bands yields a good match to the observed radial color profile. The publicly available iFIT offers an efficient tool for the non-supervised structural characterization of large galaxy samples, as those expected to become available with Euclid and LSST.


title: Investigation of the energy levels and the structure of the different states of the 24Mg nucleus

abstract: In this work, the 24Mg nucleus is considered in the cluster model by solving the Schrodinger and Klein- Gordon equations from the Nikiforov- Uvarov (NU) method. A local potential is used for these two equations that is compatible with the Hafstad-Teller potential. By substituting this potential in the Schrodinger and Klein- Gordon equations, the energy levels and wave functions are obtained from NU method. We obtain a first order equation in terms of E for the Schrodinger equation and a fourth-order equation in terms of ER for the Klein-Gordon equation. Therefore, the obtained equation from the Klein-Gordon equation has both the real solutions and imaginary solutions that we consider only the real solutions for this equation. For more accuracy, the spin-orbit and tensor potentials are added to the central potential as the perturbed terms and the first-order correction value of the energy levels is calculated. These energy levels are used for 24Mg nucleus and the good results are obtained in agreement with the experimental data. Of course, one can see that the calculated results from the Schrodinger and Klein-Gordon equations, i.e. non-relativistic and relativistic respectively, are nearly the same. By using the Coulomb repulsive potential, it is shown that the structure of the ground state and the first excited state of this nucleus is as octahedral and the configuration of the alpha particles for the second excited state is as pentahedral, i.e. as 12C +12C . Finally, we obtain the critical point for the calculation of maximum energy of this nucleus in the cluster model.


title: SCALPEL3: a scalable open-source library for healthcare claims databases

abstract: This article introduces SCALPEL3, a scalable open-source framework for studies involving Large Observational Databases (LODs). Its design eases medical observational studies thanks to abstractions allowing concept extraction, high-level cohort manipulation, and production of data formats compatible with machine learning libraries. SCALPEL3 has successfully been used on the SNDS database (see Tuppin et al. (2017)), a huge healthcare claims database that handles the reimbursement of almost all French citizens. SCALPEL3 focuses on scalability, easy interactive analysis and helpers for data flow analysis to accelerate studies performed on LODs. It consists of three open-source libraries based on Apache Spark. SCALPEL-Flattening allows denormalization of the LOD (only SNDS for now) by joining tables sequentially in a big table. SCALPEL-Extraction provides fast concept extraction from a big table such as the one produced by SCALPEL-Flattening. Finally, SCALPEL-Analysis allows interactive cohort manipulations, monitoring statistics of cohort flows and building datasets to be used with machine learning libraries. The first two provide a Scala API while the last one provides a Python API that can be used in an interactive environment. Our code is available on GitHub. SCALPEL3 allowed to extract successfully complex concepts for studies such as Morel et al (2017) or studies with 14.5 million patients observed over three years (corresponding to more than 15 billion healthcare events and roughly 15 TeraBytes of data) in less than 49 minutes on a small 15 nodes HDFS cluster. SCALPEL3 provides a sharp interactive control of data processing through legible code, which helps to build studies with full reproducibility, leading to improved maintainability and audit of studies performed on LODs.


title: Edge Probes of Topological Order

abstract: According to the bulk-edge correspondence principle, the physics of the gapless edge in the quantum Hall effect determines topological order in the gapped bulk. As the bulk is less accessible, the last two decades saw the emergence of several experimental techniques that invoke the study of the compressible edge. We review the properties of the edge, and describe several experimental techniques that include shot noise and thermal noise measurements, interferometry, and energy (thermal) transport at the edge. We pay special attention to the filling factor 5/2 in the first excited Landau level (in two-dimensional electron gas in GaAs), where experimental evidence of a non-abelian topological order was found. A brief discussion is devoted to recent interferometry experiments that uncovered unexpected physics in the integer quantum Hall effect. The chapter also addresses the theory of edge states, for systems with abelian and non-abelian topological orders.


title: Analyzing Large Receptive Field Convolutional Networks for Distant Speech Recognition

abstract: Despite significant efforts over the last few years to build a robust automatic speech recognition (ASR) system for different acoustic settings, the performance of the current state-of-the-art technologies significantly degrades in noisy reverberant environments. Convolutional Neural Networks (CNNs) have been successfully used to achieve substantial improvements in many speech processing applications including distant speech recognition (DSR). However, standard CNN architectures were not efficient in capturing long-term speech dynamics, which are essential in the design of a robust DSR system. In the present study, we address this issue by investigating variants of large receptive field CNNs (LRF-CNNs) which include deeply recursive networks, dilated convolutional neural networks, and stacked hourglass networks. To compare the efficacy of the aforementioned architectures with the standard CNN for Wall Street Journal (WSJ) corpus, we use a hybrid DNN-HMM based speech recognition system. We extend the study to evaluate the system performances for distant speech simulated using realistic room impulse responses (RIRs). Our experiments show that with fixed number of parameters across all architectures, the large receptive field networks show consistent improvements over the standard CNNs for distant speech. Amongst the explored LRF-CNNs, stacked hourglass network has shown improvements with a 8.9% relative reduction in word error rate (WER) and 10.7% relative improvement in frame accuracy compared to the standard CNNs for distant simulated speech signals.


title: Wasserstein GANs for MR Imaging: from Paired to Unpaired Training

abstract: Lack of ground-truth MR images impedes the common supervised training of neural networks for image reconstruction. To cope with this challenge, this paper leverages unpaired adversarial training for reconstruction networks, where the inputs are undersampled k-space and naively reconstructed images from one dataset, and the labels are high-quality images from another dataset. The reconstruction networks consist of a generator which suppresses the input image artifacts, and a discriminator using a pool of (unpaired) labels to adjust the reconstruction quality. The generator is an unrolled neural network -- a cascade of convolutional and data consistency layers. The discriminator is also a multilayer CNN that plays the role of a critic scoring the quality of reconstructed images based on the Wasserstein distance. Our experiments with knee MRI datasets demonstrate that the proposed unpaired training enables diagnostic-quality reconstruction when high-quality image labels are not available for the input types of interest, or when the amount of labels is small. In addition, our adversarial training scheme can achieve better image quality (as rated by expert radiologists) compared with the paired training schemes with pixel-wise loss.


title: DYTurbo: Fast predictions for Drell-Yan processes

abstract: Drell-Yan lepton pair production processes are extremely important for Standard Model (SM) precision tests and for beyond the SM searches at hadron colliders. Fast and accurate predictions are essential to enable the best use of the precision measurements of these processes; they are used for parton density fits, for the extraction of fundamental parameters of the SM, and for the estimation of background processes in searches. This paper describes a new numerical program, DYTurbo, for the calculation of the QCD transverse-momentum resummation of Drell-Yan cross sections up to next-to-next-to-leading logarithmic accuracy combined with the fixed-order results at next-to-next-to-leading order ($\mathcal{O}(\alpha_{\mathrm{S}}^2)$), including the full kinematical dependence of the decaying lepton pair with the corresponding spin correlations and the finite-width effects. The DYTurbo program is an improved reimplementation of the DYqT, DYqT and DYNNLO programs, which provides fast and numerically precise predictions through the factorisation of the cross section into production and decay variables, and the usage of quadrature rules based on interpolating functions for the integration over kinematic variables.


title: The effect of improved high-energy muon cross-sections

abstract: A precise simulation of muons with energies above a TeV is crucial for neutrino telescopes or cosmic ray experiments. To further increase the precision of these simulations, improved cross-section calculations are needed. At these energies, the interaction probability is dominated by bremsstrahlung for large energy losses and electron-positron pair production for small energy losses. Improved analytical calculations for these processes were developed with more precise screening functions of the target atom as well as higher order corrections reducing the systematic uncertainties below the percent level. The new calculations are already implemented in the new version of the lepton propagator PROPOSAL, which was designed to be highly performant for the propagation through large volumes of media using interpolation tables and to do systematic studies with its multiple available cross-section calculations. The new calculations of the differential cross-section result in a maximum deviation of 3 percent to the current standard. Their effects on the muon simulation with the resulting flux and energy loss distribution are presented.


title: Incongruences for modular forms and applications to partition functions

abstract: The study of arithmetic properties of coefficients of modular forms $f(\tau) = \sum a(n)q^n$ has a rich history, including deep results regarding congruences in arithmetic progressions. Recently, work of C.-S. Radu, S. Ahlgren, B. Kim, N. Andersen, and S. L\"{o}brich have employed the $q$-expansion principle of P. Deligne and M. Rapoport in order to determine more about where these congruences can occur. Here, we extend the method to give additional results for a large class of modular forms. We also give analogous results for generalized Frobenius partitions and the two mock theta functions $f(q)$ and $\omega(q).$


title: Adaptive Low-Rank Approximations for Operator Equations: Accuracy Control and Computational Complexity

abstract: The challenge of mastering computational tasks of enormous size tends to frequently override questioning the quality of the numerical outcome in terms of accuracy. By this we do not mean the accuracy within the discrete setting, which itself may also be far from evident for ill-conditioned problems or when iterative solvers are involved. By accuracy-controlled computation we mean the deviation of the numerical approximation from the exact solution of an underlying continuous problem in a relevant metric, which has been the initiating interest in the first place. Can the accuracy of a numerical result be rigorously certified - a question that is particularly important in the context of uncertainty quantification, when many possible sources of uncertainties interact. This is the guiding question throughout this article, which reviews recent developments of low-rank approximation methods for problems in high spatial dimensions. In particular, we highlight the role of adaptivity when dealing with such strongly nonlinear methods that integrate in a natural way issues of discrete and continuous accuracy.


title: Cosmic Ray isotopes measured by AMS02

abstract: The measurement of light isotope components in cosmic rays (CR) like3He,2H (D)and6Li are some of the most valuable tools for understanding the propagation of CR in the galaxy and constrain the models that describe it. In this work a new preliminary measurement of3He/4He,6Li/7Li and D/p flux ratios are presented, each obtained from the data of the AMS-02 experiment, along with the measurement of the isotopic component fluxes.


title: Hybrid FFT algorithm for fast demagnetization field calculations on non-equidistant magnetic layers

abstract: In micromagnetic simulations, the demagnetization field is by far the computationally most expensive field component and often a limiting factor in large multilayer systems. We present an exact method to calculate the demagnetization field of magnetic layers with arbitrary thicknesses. In this approach we combine the widely used fast-Fourier-transform based circular convolution method with an explicit convolution using a generalized form of the Newell formulas. We implement the method both for central processors and graphics processors and find that significant speedups for irregular multilayer geometries can be achieved. Using this method we optimize the geometry of a magnetic random-access memory cell by varying a single specific layer thickness and simulate a hysteresis curve to determine the resulting switching field.


title: Alleviating Bottlenecks for DNN Execution on GPUs via Opportunistic Computing

abstract: Edge computing and IoT applications are severely constrained by limited hardware resource. This makes memory consuming DNN frameworks not applicable to edge computing. Simple algorithms such as direct convolution are finding their way in embedded machine learning. As one of the most widely used platforms for DNN acceleration, GPUs face the bottleneck of on-chip bandwidth. This work introduces a GPU DNN execution architecture that targets on relieving the on-chip bandwidth bottleneck by reducing data movement through opportunistic computing. We first investigate data access patterns in the hardware view rather than the software view. Then we propose two opportunistic computing techniques to predictably perform computation when data is available with the help of assistant warps. By moving computation to data, our techniques are able to significantly reduce data movement and relieve the DNN execution bottleneck. Our evaluation results show that the proposed technique can improve DNN application performance as much as 55%.


title: Variable Metric Proximal Gradient Method with Diagonal Barzilai-Borwein Stepsize

abstract: Variable metric proximal gradient (VM-PG) is a widely used class of convex optimization method. Lately, there has been a lot of research on the theoretical guarantees of VM-PG with different metric selections. However, most such metric selections are dependent on (an expensive) Hessian, or limited to scalar stepsizes like the Barzilai-Borwein (BB) stepsize with lots of safeguarding. Instead, in this paper we propose an adaptive metric selection strategy called the diagonal Barzilai-Borwein (BB) stepsize. The proposed diagonal selection better captures the local geometry of the problem while keeping per-step computation cost similar to the scalar BB stepsize i.e. $O(n)$. Under this metric selection for VM-PG, the theoretical convergence is analyzed. Our empirical studies illustrate the improved convergence results under the proposed diagonal BB stepsize, specifically for ill-conditioned machine learning problems for both synthetic and real-world datasets.


title: The Gran Telescopio Canarias OSIRIS Broad Band First Data Release

abstract: We present the first release of GTC OSIRIS Broad Band data archive. This is an effort conducted in the framework of the Spanish Virtual Observatory to help optimize science from the Gran Telescopio Canarias Archive. Data Release 1 includes 6 788 broad-band images in the Sloan griz filters obtained between April 2009 and January 2014 and the associated catalogue with roughly 6.23 million detections of more than 630 000 unique sources. The catalogue contains standard PSF and Kron aperture photometry with a mean accuracy better than 0.09 and 0.15 mag, respectively. The relative astrometric residuals are always better than 30 mas and better than 15 mas in most cases. The absolute astrometric uncertainty of the catalogue is of 0.12 arcsec. In this paper we describe the procedure followed to build the image archive and the associated catalogue, as well as the quality tests carried out for validation. To illustrate some of the scientific potential of the catalogue, we also provide two examples of its scientific exploitation: discovery and identification of asteroids and cool dwarfs.


title: Health Monitoring in Smart Homes Utilizing Internet of Things

abstract: In recent years the concept of the Internet of Things (IoT) has evolved to connect commercial gadgets together with the medical field to facilitate an unprecedented range of accessibility. The development of medical devices connected to internet of things has been praised for the potential of alleviating the strain on the modern healthcare system by giving users the opportunity to reside in the home during treatment or recovery. With the IoT becoming more prevalent and available at a commercial level, there exists room for integration into emerging, intelligent environments such as smart homes. When used in tandem with conventional healthcare, the IoT offers a vast range of custom-tailored treatment options. This paper studies recent state-of-the-art research on the field of IoT for health monitoring and smart homes, examines several potential use-cases of blending the technology, and proposes integration with an existing smart home testbed for further study. Challenges of adoption and future research on the topic are also discussed.