title: A New Method of Modelling Tuneable Lasers with Functional Composition

abstract: A new nonlinear model is proposed for tuneable lasers. Using the generalized nonlinear Schr\"odinger equation as a starting point, expressions for the transformations undergone by the pulse are derived for each of the five components (gain, loss, dispersion, modulation, and nonlinearity) within the laser cavity. These transformations are then composed to give the overall effect of one trip around the cavity.


title: Local roughness exponent in the nonlinear molecular-beam-epitaxy universality class in one-dimension

abstract: We report local roughness exponents, $\alpha_{\text{loc}}$, for three interface growth models in one dimension which are believed to belong the non-linear molecular-beam-epitaxy (nMBE) universality class represented by the Villain-Lais-Das Sarma (VLDS) stochastic equation. We applied an optimum detrended fluctuation analysis (ODFA) [Luis et al., Phys. Rev. E 95, 042801 (2017)] and compared the outcomes with standard detrending methods. We observe in all investigated models that ODFA outperforms the standard methods providing exponents in the narrow interval $\alpha_{\text{loc}}\in[0.96,0.98]$ consistent with renormalization group predictions for the VLDS equation. In particular, these exponent values are calculated for the Clarke-Vvdensky and Das Sarma-Tamborenea models characterized by very strong corrections to the scaling, for which large deviations of these values had been reported. Our results strongly support the absence of anomalous scaling in the nMBE universality class and the existence of corrections in the form $\alpha_{\text{loc}}=1-\epsilon$ of the one-loop renormalization group analysis of the VLDS equation.


title: Kernel Transformer Networks for Compact Spherical Convolution

abstract: Ideally, 360{\deg} imagery could inherit the deep convolutional neural networks (CNNs) already trained with great success on perspective projection images. However, existing methods to transfer CNNs from perspective to spherical images introduce significant computational costs and/or degradations in accuracy. In this work, we present the Kernel Transformer Network (KTN). KTNs efficiently transfer convolution kernels from perspective images to the equirectangular projection of 360{\deg} images. Given a source CNN for perspective images as input, the KTN produces a function parameterized by a polar angle and kernel as output. Given a novel 360{\deg} image, that function in turn can compute convolutions for arbitrary layers and kernels as would the source CNN on the corresponding tangent plane projections. Distinct from all existing methods, KTNs allow model transfer: the same model can be applied to different source CNNs with the same base architecture. This enables application to multiple recognition tasks without re-training the KTN. Validating our approach with multiple source CNNs and datasets, we show that KTNs improve the state of the art for spherical convolution. KTNs successfully preserve the source CNN's accuracy, while offering transferability, scalability to typical image resolutions, and, in many cases, a substantially lower memory footprint.


title: Complex poles and spectral function of Yang-Mills theory

abstract: We derive general relationships between the number of complex poles of a propagator and the sign of the spectral function originating from the branch cut in the Minkowski region under some assumptions on the asymptotic behaviors of the propagator. We apply this relation to the mass-deformed Yang-Mills model with one-loop quantum corrections, which is identified with a low-energy effective theory of the Yang-Mills theory, to show that the gluon propagator in this model has a pair of complex conjugate poles or "tachyonic" poles of multiplicity two, in accordance with the fact that the gluon field has a negative spectral function, while the ghost propagator has at most one "unphysical" pole. Finally, we discuss implications of these results for gluon confinement and other non-perturbative aspects of the Yang-Mills theory.


title: Covariant holographic entanglement negativity for disjoint intervals in $AdS_3/CFT_2$

abstract: We advance a construction for the covariant holographic entanglement negativity for time dependent mixed states of disjoint intervals in $(1+1)$ dimensional conformal field theories ($CFT_{1+1}$) dual to bulk non static $AdS_3$ geometries. Application of our proposal to such mixed states in a $CFT_{1+1}$ dual to bulk non extremal and extremal rotating BTZ black holes exactly reproduces the replica technique results in the large central charge limit. We also investigate the time dependent holographic entanglement negativity for such mixed states in a $CFT_{1+1}$ dual to a bulk Vaidya-$AdS_3$ geometry in the context of their thermalization involving bulk black hole formation.


title: Classical Channel Gravity in the Newtonian Limit

abstract: We present a minimal model for the quantum evolution of matter under the influence of classical gravity in the Newtonian limit. Based on a continuous measurement-feedback channel that acts simultaneously on all constituent masses of a given quantum system, the model scales and applies consistently to arbitrary mass densities, and it recovers the classical Newton force between macroscopic masses. The concomitant loss of coherence is set by a model parameter, does not depend on mass, and can thus be confined to unobservable time scales for micro- and macroscopic systems alike. The model can be probed in high-precision matter-wave interferometry, and ultimately tested in recently proposed optomechanical quantum gravity experiments.


title: Ground-based detection of an extended helium atmosphere in the Saturn-mass exoplanet WASP-69b

abstract: Hot gas giant exoplanets can lose part of their atmosphere due to strong stellar irradiation, affecting their physical and chemical evolution. Studies of atmospheric escape from exoplanets have mostly relied on space-based observations of the hydrogen Lyman-{\alpha} line in the far ultraviolet which is strongly affected by interstellar absorption. Using ground-based high-resolution spectroscopy we detect excess absorption in the helium triplet at 1083 nm during the transit of the Saturn-mass exoplanet WASP-69b, at a signal-to-noise ratio of 18. We measure line blue shifts of several km/s and post transit absorption, which we interpret as the escape of part of the atmosphere trailing behind the planet in comet-like form. [Additional notes by authors: Furthermore, we provide upper limits for helium signals in the atmospheres of the exoplanets HD 209458b, KELT-9b, and GJ 436b. We investigate the host stars of all planets with detected helium signals and those of the three planets we derive upper limits for. In each case we calculate the X-ray and extreme ultraviolet flux received by these planets. We find that helium is detected in the atmospheres of planets (orbiting the more active stars and) receiving the larger amount of irradiation from their host stars.]


title: Black Holes, Heavy States, Phase Shift and Anomalous Dimensions

abstract: We compute the phase shift of a highly energetic particle traveling in the background of an asymptotically AdS black hole. In the dual CFT, the phase shift is related to a four point function in the Regge limit. The black hole mass is translated to the ratio between the conformal dimension of a heavy operator and the central charge. This ratio serves as a useful expansion parameter; its power measures the number of stress tensors appearing in the intermediate channel. We compute the leading term in the phase shift in a holographic CFT of arbitrary dimensionality using Conformal Regge Theory and observe complete agreement with the gravity result. In a two-dimensional CFT with a large central charge the heavy-heavy-light-light Virasoro vacuum block reproduces the gravity phase shift to all orders in the expansion parameter. We show that the leading order phase shift is related to the anomalous dimensions of certain double trace operators and verify this agreement using known results for the latter. We also perform a separate gravity calculation of these anomalous dimensions to second order in the expansion parameter and compare with the phase shift expansion.


title: Variable selection in high-dimensional linear model with possibly asymmetric or heavy-tailed errors

abstract: We consider the problem of automatic variable selection in a linear model with asymmetric or heavy-tailed errors when the number of explanatory variables diverges with the sample size. For this high-dimensional model, the penalized least square method is not appropriate and the quantile framework makes the inference more difficult because to the non differentiability of the loss function. We propose and study an estimation method by penalizing the expectile process with an adaptive LASSO penalty. Two cases are considered: the number of model parameters is smaller and afterwards larger than the sample size, the two cases being distinct by the adaptive penalties considered. For each case we give the rate convergence and establish the oracle properties of the adaptive LASSO expectile estimator. The proposed estimators are evaluated through Monte Carlo simulations and compared with the adaptive LASSO quantile estimator. We applied also our estimation method to real data in genetics when the number of parameters is greater than the sample size.


title: Trade-off between sensitivity and selectivity in olfactory receptor neuron

abstract: It was observed before that due to convergence in the olfactory system a possible amplification can be as large as the degree of convergence. This is in the case when a single impulse from the converging inputs is enough to trigger the secondary neuron. On the other hand, if a number of impulses are required for triggering, a gain in discriminating ability may be obtained along with decrease in sensitivity gained due to the convergence. We discuss this trade-off in terms of concrete estimates using olfactory sensory neuron and the set of its receptor proteins as an example of system with convergence.


title: Deep Variational Transfer: Transfer Learning through Semi-supervised Deep Generative Models

abstract: In real-world applications, it is often expensive and time-consuming to obtain labeled examples. In such cases, knowledge transfer from related domains, where labels are abundant, could greatly reduce the need for extensive labeling efforts. In this scenario, transfer learning comes in hand. In this paper, we propose Deep Variational Transfer (DVT), a variational autoencoder that transfers knowledge across domains using a shared latent Gaussian mixture model. Thanks to the combination of a semi-supervised ELBO and parameters sharing across domains, we are able to simultaneously: (i) align all supervised examples of the same class into the same latent Gaussian Mixture component, independently from their domain; (ii) predict the class of unsupervised examples from different domains and use them to better model the occurring shifts. We perform tests on MNIST and USPS digits datasets, showing DVT's ability to perform transfer learning across heterogeneous datasets. Additionally, we present DVT's top classification performances on the MNIST semi-supervised learning challenge. We further validate DVT on a astronomical datasets. DVT achieves states-of-the-art classification performances, transferring knowledge across real stars surveys datasets, EROS, MACHO and HiTS, . In the worst performance, we double the achieved F1-score for rare classes. These experiments show DVT's ability to tackle all major challenges posed by transfer learning: different covariate distributions, different and highly imbalanced class distributions and different feature spaces.


title: unVEil the darknesS of The gAlactic buLgE (VESTALE)

abstract: The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.


title: Taking the Scenic Route: Automatic Exploration for Videogames

abstract: Machine playtesting tools and game moment search engines require exposure to the diversity of a game's state space if they are to report on or index the most interesting moments of possible play. Meanwhile, mobile app distribution services would like to quickly determine if a freshly-uploaded game is fit to be published. Having access to a semantic map of reachable states in the game would enable efficient inference in these applications. However, human gameplay data is expensive to acquire relative to the coverage of a game that it provides. We show that off-the-shelf automatic exploration strategies can explore with an effectiveness comparable to human gameplay on the same timescale. We contribute generic methods for quantifying exploration quality as a function of time and demonstrate our metric on several elementary techniques and human players on a collection of commercial games sampled from multiple game platforms (from Atari 2600 to Nintendo 64). Emphasizing the diversity of states reached and the semantic map extracted, this work makes productive contrast with the focus on finding a behavior policy or optimizing game score used in most automatic game playing research.


title: Polarization diversity phase modulator for measuring frequency-bin entanglement of a biphoton frequency comb in a depolarized channel

abstract: Phase modulation has emerged as a technique to create and manipulate high-dimensional frequency-bin entanglement. A necessary step to extending this technique to depolarized channels, such as those in a quantum networking environment, is the ability to perform phase modulation independent of photon polarization. This also necessary to harness hypertanglement in the polarization and frequency degrees of freedom for operations like Bell state discrimination. However, practical phase modulators are generally sensitive to the polarization of light and this makes them unsuited to such applications. We overcome this limitation by implementing a polarization diversity scheme to measure frequency-bin entanglement in arbitrarily polarized photon pairs.


title: Induced graphs of uniform spanning forests

abstract: Given a subgraph $H$ of a graph $G$, the induced graph of $H$ is the largest subgraph of $G$ whose vertex set is the same as that of $H$. Our paper concerns the induced graphs of the components of $\operatorname{WSF}(G)$, the wired spanning forest on $G$, and, to a lesser extent, $\operatorname{FSF}(G)$, the free uniform spanning forest. We show that the induced graph of each component of $\operatorname{WSF}(\mathbb Z^d$) is almost surely recurrent when $d\ge 8$. Moreover, the effective resistance between two points on the ray of the tree to infinity within a component grows linearly when $d\ge9$. For any vertex-transitive graph $G$, we establish the following resampling property: Given a vertex $o$ in $G$, let $\mathcal T_o$ be the component of $\operatorname{WSF}(G)$ containing $o$ and $\overline{\mathcal{T}_o}$ be its induced graph. Conditioned on $\overline{\mathcal{T}_o}$, the tree $\mathcal T_o$ is distributed as $\operatorname{WSF}(\overline{\mathcal{T}_o})$. For any graph $G$, we also show that if $\mathcal T_o$ is the component of $\operatorname{FSF}(G)$ containing $o$ and $\overline{\mathcal{T}_o}$ is its induced graph, then conditioned on $\overline{\mathcal{T}_o}$, the tree $\mathcal T_o$ is distributed as $\operatorname{FSF}(\overline{\mathcal{T}_o})$.


title: Backdooring Convolutional Neural Networks via Targeted Weight Perturbations

abstract: We present a new type of backdoor attack that exploits a vulnerability of convolutional neural networks (CNNs) that has been previously unstudied. In particular, we examine the application of facial recognition. Deep learning techniques are at the top of the game for facial recognition, which means they have now been implemented in many production-level systems. Alarmingly, unlike other commercial technologies such as operating systems and network devices, deep learning-based facial recognition algorithms are not presently designed with security requirements or audited for security vulnerabilities before deployment. Given how young the technology is and how abstract many of the internal workings of these algorithms are, neural network-based facial recognition systems are prime targets for security breaches. As more and more of our personal information begins to be guarded by facial recognition (e.g., the iPhone X), exploring the security vulnerabilities of these systems from a penetration testing standpoint is crucial. Along these lines, we describe a general methodology for backdooring CNNs via targeted weight perturbations. Using a five-layer CNN and ResNet-50 as case studies, we show that an attacker is able to significantly increase the chance that inputs they supply will be falsely accepted by a CNN while simultaneously preserving the error rates for legitimate enrolled classes.


title: Ab-initio approach to the stability and the structural, electronic and magnetic properties of the (001) Znfe2O4 surface terminations

abstract: We present a Density Functional Theory (DFT) based study of the structural and magnetic properties of the (001) surface of the semiconducting oxide ZnFe2O4 (spinel structure). The calculations were performed using the DFT based ab initio plane wave and pseudopotential method as implemented in the Quantum Espresso code. The all electron Full-potential linearized-augmented-plane-wave method (FP-LAPW) was also employed to check the accuracy of plane wave method. In both calculations the DFT+U methodology was employed and different (001) surface terminations of ZnFe2O4 were studied: We find that the surface terminated in Zn is the stable one. For all the (001) surface terminations our calculations predict that the Zn-Fe cationic inversion (antisites), which are defects in bulk ZnFe2O4, becomes stable and an integral part of the surface. Also, a ferrimagnetic behavior is predicted for the case of antisites in the superficial layer. Our results for different properties of the surface of ZnFe2O4 are compared with those obtained in bulk samples and those reported in the literature.


title: Spin-states in multiphoton pair production for circularly polarized light

abstract: Scalar and fermionic particle pair production in rotating electric fields is investigated in the nonperturbative multiphoton regime. Angular momentum distribution functions in above-threshold pair production processes are calculated numerically within quantum kinetic theory and discussed on the basis of a photon absorption model. The particle spectra can be understood if the spin states of the particle-antiparticle pair are taken into account.

categories: hep-ph physics.atom-ph physics.plasm-ph quant-ph


title: Playing with and against Hedge

abstract: Hedge has been proposed as an adaptive scheme, which guides an agent's decision in resource selection and distribution problems that can be modeled as a multi-armed bandit full information game. Such problems are encountered in the areas of computer and communication networks, e.g. network path selection, load distribution, network interdiction, and also in problems in the area of transportation. We study Hedge under the assumption that the total loss that can be suffered by the player in each round is upper bounded. In this paper, we study the worst performance of Hedge.


title: Lattice formulation of axion inflation. Application to preheating

abstract: We present a lattice formulation of an interaction ${\phi\over\Lambda} F\tilde F$ between an axion and some $U(1)$ gauge sector with the following properties: it reproduces the continuum theory up to $\mathcal{O}(dx_\mu^2)$ corrections, it preserves exact gauge invariance and shift symmetry on the lattice, and it is suitable for self-consistent expansion of the Universe. The lattice equations of motion can no longer be solved by explicit methods, but we propose an implicit method to overcome this difficulty, which preserves the relevant system constraints down to arbitrary (tunable) precision. As a first application we study, in a comoving grid in $(3+1)$ dimensions, the last efolds of axion-inflation with quadratic potential and the preheating stage following afterwards. We fully account for the inhomogeneity and non-linearity of the system, including the gauge field contribution to the expansion rate of the Universe and its backreaction into the axion dynamics. We characterize in detail, as a function of the coupling, the energy transfer from the axion to the gauge field. Two coupling regimes are identified, sub- and super-critical, depending on whether the final energy fraction stored in the gauge field is below or above $\sim 50\%$ of the total energy. The Universe is very efficiently reheated for super-critical couplings, rapidly entering in a radiation dominated stage. Our results on preheating confirm previously published results.


title: The integral trace form as a complete invariant for real $S_n$ number fields

abstract: In the past the first named author has studied to what extent the integral trace can characterize a number field beyond what the discriminant does. The cases of cyclic number fields and non-totally real fields are more or less settled, concluding that for such fields the integral trace does not always characterize the field. In this paper we show that the integral trace is a complete invariant for degree $n$, $S_n$ real number number fields that satisfy certain ramification bound. Among the real $S_n$ fields that our results cover, there are those of square free different ideal. Moreover, for such fields we find an explicit description of the isometry group of the integral trace.


title: Bounds for the number of cohomological automorphic representations of $GL_3/\mathbb{Q}$ in the weight aspect

abstract: We prove a power saving over the trivial bound for the number of cohomological cuspidal automorphic representations of fixed level and growing weight on $GL_3/\mathbb{Q}$, by adapting the methods of our earlier paper on $GL_2$.


title: Young Stars and their Variability with LSST

abstract: Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.


title: A Complementary Third Law for Black Hole Thermodynamics

abstract: There are some examples in the literature, in which despite the fact that the underlying theory or model does not impose a lower bound on the size of black holes, the final temperature under Hawking evaporation is nevertheless finite and nonzero. We show that under some loose conditions, the black hole is necessarily an effective remnant, in the sense that its evaporation time is infinite. That is, the final state that there is nonzero finite temperature despite having no black hole remaining cannot be realized. We discuss the limitations, subtleties, and the implications of this result, which is reminiscent of the third law of black hole thermodynamics, but with the roles of temperature and size interchanged. We therefore refer to our result as the "complementary third law" for black hole thermodynamics.


title: The Diverse Science Return from a Wide-Area Survey of the Galactic Plane

abstract: The overwhelming majority of objects visible to LSST lie within the Galactic Plane. Though many previous surveys have avoided this region for fear of stellar crowding, LSST's spatial resolution combined with its state-of-the-art Difference Image Analysis mean that it can conduct a high cadence survey of most of the Galaxy for the first time. Here we outline the many areas of science that would greatly benefit from an LSST survey that included the Galactic Plane, Magellanic Clouds and Bulge at a cadence of 2-3 d. Particular highlights include measuring the mass spectrum of black holes, and mapping the population of exoplanets in the Galaxy in relation to variations in star forming environments. But the same survey data will provide a goldmine for a wide range of science, and we explore possible survey strategies which maximize the scientific return for a number of fields including young stellar objects, cataclysmic variables and Neptune Trojans.


title: Entanglement Suppression and Emergent Symmetries of Strong Interactions

abstract: Entanglement suppression in the strong interaction $S$-matrix is shown to be correlated with approximate spin-flavor symmetries that are observed in low-energy baryon interactions, the Wigner $SU(4)$ symmetry for two flavors and an $SU(16)$ symmetry for three flavors. We conjecture that dynamical entanglement suppression is a property of the strong interactions in the infrared, giving rise to these emergent symmetries and providing powerful constraints on the nature of nuclear and hypernuclear forces in dense matter.


title: Mapping the Periphery and Variability of the Magellanic Clouds

abstract: We propose a combination of a modified Wide-Fast-Deep survey, a mini-survey of the South Celestial Pole, and a Deep Drilling-style survey to produce a 3-D map of the Magellanic System and to provide a detailed census of the transient and variable populations in the Clouds. We support modifying the Wide-Fast-Deep survey to cover the declination range $-72.25\deg<{\rm Dec}<12.4\deg$ and the Galactic latitude range $|b|>15\deg$, as proposed in a separate white paper. We additionally propose a mini-survey covering the 950$\deg^2$ with ${\rm Dec} < -72.25$ in $ugriz$ to the standard LSST single-exposure depth and with 40 visits per filter per field. Finally, we propose a mini-survey covering $\sim100 \deg^2$ of the main bodies of the Clouds with twelve total pointings, 2000 total visits per field, and shorter exposure time.


title: Local convergence of large random triangulations coupled with an Ising model

abstract: We prove the existence of the local weak limit of the measure obtained by sampling random triangulations of size $n$ decorated by an Ising configuration with a weight proportional to the energy of this configuration. To do so, we establish the algebraicity and the asymptotic behaviour of the partition functions of triangulations with spins for any boundary condition. In particular, we show that these partition functions all have the same phase transition at the same critical temperature. Some properties of the limiting object -- called the Infinite Ising Planar Triangulation -- are derived, including the recurrence of the simple random walk at the critical temperature.


title: An atomically thin oxide layer on the elemental superconductor Ta(001) surface

abstract: Recently the oxygen-reconstructed tantalum surface Ta(001)-p(3$\times$3)-O has experienced considerable attention due its use as a potential platform for Majorana physics in adatom chains. Experimental studies using scanning tunneling microscopy and spectroscopy found rich atomic and electronic structures already for the clean Ta(001)-O surface, which we combine here with $ab~initio$ methods. We discover two metastable superstructures at the root of the different topographic patterns, discuss its emergence during annealing, and identify the electronic properties. The latter is determined as the sole origin for the contrast reversal seen at positive bias. The observed effects are essentially connected to the two distinct oxygen states appearing on the surface in different geometries. The second superstructure was found in simulations by introducing oxygen vacancies, what was also observed in tantalum pentoxide systems. Additionally we study the charge distribution on the oxidized surface and underline its importance for the adsorption process of polarizable atoms and molecules.


title: A Cadence to Reduce Aliasing in LSST

abstract: Regular sampling in the time domain results in aliasing in the frequency domain that complicates the accurate determination of the periods of astrophysical variables. We propose to actively break the regularity of this sampling by providing an additional consideration for the scheduler that weights fields according to when observations will contribute the least to aliasing. The current aliases for each field can be computed during daytime from the history of observations. We can then calculate the times when additional observations would worsen or alleviate these aliases for different fields. The scheduler should give preference to observation epochs that lessen the effect of aliasing, while still meeting all other cadence requirements.


title: Continuous Cadence Acquisition of the LSST Deep Drilling Fields

abstract: To extend LSST's coverage of the transient and variable sky down to minute timescales, we propose that observations of the Deep Drilling Fields are acquired in sequences of continuous exposures each lasting 2--4 hours. This will allow LSST to resolve rapid stellar variability such as short-period pulsations, exoplanet transits, ultracompact binary systems, and flare morphologies, while still achieving the desired co-added depths for the selected fields. The greater number of observations of each Deep Drilling Field pushes these mini-surveys deep in terms of both sensitivity to low-amplitude variability and co-added depth. Saving the individual 15-second exposures will yield an effective Nyquist limit of $\approx0.031$ Hz (32 seconds). Resolved short-period variability of targets in these fields will aid the interpretation of sparse observations of a greater number of variables in the main survey. If this cadence strategy conflicts with the science goals of individual Deep Drilling Fields, at least a subset of the additional observations of each field should be obtained continuously. This strategy should also be considered for the proposed Galactic Plane mini survey, which will observe a greater number of stellar variables and transients.


title: Large Synoptic Survey Telescope White Paper; The Case for Matching U-band on Deep Drilling Fields

abstract: U-band observations with the LSST have yet to be fully optimized in cadence. The straw man survey design is a simple coverage of the medium-deep-fast survey. Here we argue that deep coverage of the four deep drilling fields (XMM-LSS, ECDFS, ELAIS-S1 and COSMOS) has a much higher scientific return, given that these are also the target of the Southern Hemisphere's Square Kilometer Array Pathfinder, the MeerKAT specifically, deep radio observations.


title: Extending ROOT through Modules

abstract: The ROOT software framework is foundational for the HEP ecosystem, providing capabilities such as IO, a C++ interpreter, GUI, and math libraries. It uses object-oriented concepts and build-time components to layer between them. We believe additional layering formalisms will benefit ROOT and its users. We present the modularization strategy for ROOT which aims to formalize the description of existing source components, making available the dependencies and other metadata externally from the build system, and allow post-install additions of functionality in the runtime environment. components can then be grouped into packages, installable from external repositories to deliver post-install step of missing packages. This provides a mechanism for the wider software ecosystem to interact with a minimalistic install. Reducing intra-component dependencies improves maintainability and code hygiene. We believe helping maintain the smallest "base install" possible will help embedding use cases. The modularization effort draws inspiration from the Java, Python, and Swift ecosystems. Keeping aligned with the modern C++, this strategy relies on forthcoming features such as C++ modules. We hope formalizing the component layer will provide simpler ROOT installs, improve extensibility, and decrease the complexity of embedding in other ecosystems


title: Presto-Color: A Photometric Survey Cadence for Explosive Physics & Fast Transients

abstract: We identify minimal observing cadence requirements that enable photometric astronomical surveys to detect and recognize fast and explosive transients and fast transient features. Observations in two different filters within a short time window (e.g., g-and-i, or r-and-z, within < 0.5 hr) and a repeat of one of those filters with a longer time window (e.g., > 1.5 hr) are desirable for this purpose. Such an observing strategy delivers both the color and light curve evolution of transients on the same night. This allows the identification and initial characterization of fast transient -- or fast features of longer timescale transients -- such as rapidly declining supernovae, kilonovae, and the signatures of SN ejecta interacting with binary companion stars or circumstellar material. Some of these extragalactic transients are intrinsically rare and generally all hard to find, thus upcoming surveys like the Large Synoptic Survey Telescope (LSST) could dramatically improve our understanding of their origin and properties. We colloquially refer to such a strategy implementation for the LSST as the Presto-Color strategy (rapid-color). This cadence's minimal requirements allow for overall optimization of a survey for other science goals.


title: Vortices of $SO(2)$ gauged Skyrmions in $2+1$ dimensions

abstract: Vortices the $SO(2)$ gauged planar Skyrme model, with a) only Maxwell, b) only Chern-Simons, and c) both Maxwell and Chern-Simons dynamics are studied systematically. In cases a) and b), where both models feature a single parameter $\lambda$ (the coupling of the potential term), the dependence of the energy on $\lambda$ is analysed. It is shown that the plots of the energy $vs.$ $\lambda$ feature discontinuities and branches. In case c), the emphasis is on the evolution of the topological charge, taking non-integer values. Throughout, the properties studied are contrasted with those of the corresponding Abelian Higgs models.


title: A Higher Cadence Subsurvey Located in the Galactic Plane

abstract: Presently, the Galactic plane receives relatively few observations compared to most of the LSST footprint. While this may address static science, the plane will also represent the highest density of variable Galactic sources. The proper characterization of variability of these sources will benefit greatly from observations at a higher cadence.


title: Continuous Performance Benchmarking Framework for ROOT

abstract: Foundational software libraries such as ROOT are under intense pressure to avoid software regression, including performance regressions. Continuous performance benchmarking, as a part of continuous integration and other code quality testing, is an industry best-practice to understand how the performance of a software product evolves over time. We present a framework, built from industry best practices and tools, to help to understand ROOT code performance and monitor the efficiency of the code for a several processor architectures. It additionally allows historical performance measurements for ROOT I/O, vectorization and parallelization sub-systems.


title: A simple approach to construct confidence bands for a regression function with incomplete data

abstract: A long-standing problem in the construction of asymptotically correct confidence bands for a regression function $m(x)=E[Y|X=x]$, where $Y$ is the response variable influenced by the covariate $X$, involves the situation where $Y$ values may be missing at random, and where the selection probability, the density function $f(x)$ of $X$, and the conditional variance of $Y$ given $X$ are all completely unknown. This can be particularly more complicated in nonparametric situations. In this paper we propose a new kernel-type regression estimator and study the limiting distribution of the properly normalized versions of the maximal deviation of the proposed estimator from the true regression curve. The resulting limiting distribution will be used to construct uniform confidence bands for the underlying regression curve with asymptotically correct coverages. The focus of the current paper is on the case where $X\in \mathbb{R}$. We also perform numerical studies to assess the finite-sample performance of the proposed method. In this paper, both mechanics and the theoretical validity of our methods are discussed.


title: Blazars and Fast Radio Bursts with LSST

abstract: The aim of this white paper is to discuss the observing strategies for the LSST Wide-Fast-Deep that would improve the study of blazars (emission variability, census, environment) and Fast Radio Bursts (FRBs). For blazars, these include the adoption of: i) a reference filter to allow reconstruction of a well-sampled light curve not affected by colour changes effects; ii) two snapshots/visit with different exposure times to avoid saturation during flaring states; iii) a rolling cadence to get better-sampled light curves at least in some time intervals. We also address the potential importance of Target of Opportunity (ToO) observations of blazar neutrino sources, and the advantages of a Minisurvey with a star trail cadence (see white paper by David Thomas et al.) for both the blazar science and the detection of possible very fast optical counterparts of FRBs.


title: Radial Acceleration and Tully-Fisher Relations in Conformal Gravity

abstract: In 2016 McGaugh, Lelli and Schombert established a universal Radial Acceleration Relation for centripetal accelerations in spiral galaxies. Their work showed a strong correlation between observed centripetal accelerations and those predicted by luminous Newtonian matter alone. Through the use of the fitting function that they introduced, mass discrepancies in spiral galaxies can be constrained in a uniform manner that is completely determined by the baryons in the galaxies. Here we present a new empirical plot of the observed centripetal accelerations and the luminous Newtonian expectations, which more than doubles the number of observed data points considered by McGaugh et al. while retaining the Radial Acceleration Relation. If this relation is not to be due to dark matter, it would then have to be due to an alternate gravitational theory that departs from Newtonian gravity in some way. In this paper we show how the candidate alternate conformal gravity theory can provide a natural description of the Radial Acceleration Relation, without any need for dark matter or its free halo parameters. We discuss how the empirical Tully-Fisher relation follows as a consequence of conformal gravity.


title: The 11 years of low-activity of the magnetar XTE J1810-197

abstract: In 2003, the magnetar XTE J1810-197 started an outburst that lasted until early 2007. In the following 11 years, the source stayed in a quiescent/low activity phase. XTE J1810-197 is one of the closest magnetars, hence its X-ray properties can be studied in detail even in quiescence and an extended monitoring has been carried out to study its long term timing and spectral evolution. Here, we report the results of new X-ray observations, taken between September 2017 and April 2018, with XMM-Newton, Chandra and Nicer. We derived a phase-connected timing solution yielding a frequency derivative of -9.26(6)x10^-14 Hz s-1. This value is consistent with that measured between 2009 and 2011, indicating that the pulsar spin-down rate remained quite stable during the long quiescent period. A spectral analysis of all the X-ray observations taken between 2009 and 2018 does not reveal significant spectral and/or flux variability. The spectrum of XTE J1810-197 can be described by the sum of two thermal components with temperatures of 0.15 and 0.3 keV, plus a power law component with photon index 0.6. We also found evidence for an absorption line at ~1.2 keV and width of 0.1 keV. Thanks to the long exposure time of the summed XMM-Newton observations, we could also carry out a phase-resolved spectral analysis for this source in quiescence. This showed that the flux modulation can be mainly ascribed to the the warmer of the two thermal components, whose flux varies by ~45 per cent along the pulse phase.


title: Pushing Automated Abundance Derivations Into the Cool Star Regime: A Test Using Three G and Two K Stars in Praesepe

abstract: We present the results of an abundance analysis of three G and two K dwarfs in the Praesepe open cluster based on high-resolution, moderate signal-to-noise-ratio spectra obtained with the ARC 3.5-m Telescope at Apache Point Observatory. Using a Principle Component Analysis and the BACCHUS automated spectral analysis code, we determined stellar parameters and abundances of up to 24 elements for each of our targets, which range in temperature from 6000 to 4600 K. The average derived iron abundance for the three G stars is 0.17+/-0.07 dex, consistent with the 0.12+/-0.04 dex derived by Boesgaard et al. (2013) for their sample of 11 solar-type Praesepe members, which included these G stars. To investigate the efficacy of using automated routines to derive the abundances of cooler main-sequence stars, we compared the abundances of the K dwarfs to those of the G dwarfs. Our abundances are consistent to <=0.1~dex for 13 of the 18 elements we report for all five of the stars, providing more evidence that G and K stars in a given open cluster are chemically homogeneous. The median difference between the mean G and K stars abundances is 0.08+/-0.05 dex, despite serious challenges with the noisier data for the fainter K dwarfs. Our results are encouraging for chemical tagging, as they indicate that it may be possible to use automated abundance determination techniques to identify chemically related main-sequence stars across larger temperature ranges than are usually considered in these experiments.


title: Kernelization of Packing Problems

abstract: Kernelization algorithms are polynomial-time reductions from a problem to itself that guarantee their output to have a size not exceeding some bound. For example, d-Set Matching for integers d>2 is the problem of finding a matching of size at least k in a given d-uniform hypergraph and has kernels with O(k^d) edges. Bodlaender et al. [JCSS 2009], Fortnow and Santhanam [JCSS 2011], Dell and Van Melkebeek [JACM 2014] developed a framework for proving lower bounds on the kernel size for certain problems, under the complexity-theoretic hypothesis that coNP is not contained in NP/poly. Under the same hypothesis, we show tight lower bounds for the kernelization of d-Set Matching and other packing problems. Our bounds are tight for d-Set Matching: It does not have kernels with O(k^{d-{\epsilon}}) edges for any {\epsilon}>0 unless the hypothesis fails. By reduction, this transfers to a bound of O(k^{d-1-{\epsilon}}) for the problem of finding k vertex-disjoint cliques of size d in standard graphs. Obtaining tight bounds for graph packing problems is challenging: We make first progress in this direction by showing non-trivial kernels with O(k^2.5) edges for the problem of finding k vertex-disjoint paths of three edges each. If the paths have d edges each, we improve the straightforward O(k^{d+1}) kernel to a uniform polynomial kernel where the exponent of the kernel size is independent of k. Most of our lower bound proofs follow a general scheme that we discover: To exclude kernels of size O(k^{d-{\epsilon}}) for a problem in d-uniform hypergraphs, one should reduce from a carefully chosen d-partite problem that is still NP-hard. As an illustration, we apply this scheme to the vertex cover problem, which allows us to replace the number-theoretical construction by Dell and Van Melkebeek [JACM 2014] with shorter elementary arguments.


title: Scheduling a Human Channel

abstract: We consider a system where a human operator processes a sequence of tasks that are similar in nature under a total time constraint. In these systems, the performance of the operator depends on its past utilization. This is akin to $\textit{state-dependent}$ channels where the past actions of the transmitter affects the future quality of the channel (also known as $\textit{action-dependent}$ or $\textit{use-dependent}$ channels). For $\textit{human channels}$, a well-known psychological phenomena, known as $\textit{Yerkes-Dodson law}$, states that a human operator performs worse when he/she is over-utilized or under-utilized. Over such a $\textit{use-dependent}$ human channel, we consider the problem of maximizing a utility function, which is monotonically increasing and concave in the time allocated for each task, under explicit minimum and maximum $\textit{utilization}$ constraints. We show that the optimal solution is to keep the utilization ratio of the operator as high as possible, and to process all the tasks. We prove that the optimal policy consists of two major strategies: utilize the operator without resting until reaching the maximum allowable utilization ratio, and then alternate between working and resting the operator each time reaching the maximum allowable utilization at the end of work-period. We show that even though the tasks are similar in difficulty, the time allocated for the tasks can be different depending on the strategy in which a task is processed; however, the tasks processed in the same strategy are processed equally.

categories: cs.IT cs.NI cs.SY eess.SP math.IT math.OC


title: On top Fourier coefficients of certain automorphic representations of GLn

abstract: In this paper, we study top Fourier coefficients of certain automorphic representations of $\mathrm{GL}_n(\mathbb{A})$. In particular, we prove a conjecture of Jiang on top Fourier coefficients of isobaric automorphic representations of $\mathrm{GL}_n(\mathbb{A})$ of form $$ \Delta(\tau_1, b_1) \boxplus \Delta(\tau_2, b_2) \boxplus \cdots \boxplus \Delta(\tau_r, b_r)\,, $$ where $\Delta(\tau_i,b_i)$'s are Speh representations in the discrete spectrum of $\mathrm{GL}_{a_ib_i}(\mathbb{A})$ with $\tau_i$'s being unitary cuspidal representations of $\mathrm{GL}_{a_i}(\mathbb{A})$, and $n = \sum_{i=1}^r a_ib_i$. Endoscopic lifting images of the discrete spectrum of classical groups form a special class of such representations. The result of this paper will facilitate the study of automorphic forms of classical groups occurring in the discrete spectrum.


title: Generation and sampling of quantum states of light in a silicon chip

abstract: Implementing large instances of quantum algorithms requires the processing of many quantum information carriers in a hardware platform that supports the integration of different components. While established semiconductor fabrication processes can integrate many photonic components, the generation and algorithmic processing of many photons has been a bottleneck in integrated photonics. Here we report the on-chip generation and processing of quantum states of light with up to eight photons in quantum sampling algorithms. Switching between different optical pumping regimes, we implement the Scattershot, Gaussian and standard boson sampling protocols in the same silicon chip, which integrates linear and nonlinear photonic circuitry. We use these results to benchmark a quantum algorithm for calculating molecular vibronic spectra. Our techniques can be readily scaled for the on-chip implementation of specialised quantum algorithms with tens of photons, pointing the way to efficiency advantages over conventional computers.


title: On completely regular codes of covering radius 1 in the halved hypercubes

abstract: We consider constructions of covering-radius-1 completely regular codes, or, equivalently, equitable 2-partitions (regular 2-partitions, perfect 2-colorings), of halved n-cubes. Keywords: completely regular code, equitable partition, regular partition, partition design, perfect coloring, halved hypercube.


title: On generation of node distributions for meshless PDE discretizations

abstract: In this paper we present an algorithm that is able to generate locally regular node layouts with spatially variable nodal density for interiors of arbitrary domains in two, three and higher dimensions. It is demonstrated that the generated node distributions are suitable to use in the RBF-FD method, which is demonstrated by solving thermo-fluid problem in 2D and 3D. Additionally, local minimal spacing guarantees are proven for both uniform and variable nodal densities. The presented algorithm has time complexity $O(N)$ to generate $N$ nodes with constant nodal spacing and $O(N \log N)$ to generate variably spaced nodes. Comparison with existing algorithms is performed in terms of node quality, time complexity, execution time and PDE solution accuracy.


title: A strategy for LSST to unveil a population of kilonovae without gravitational-wave triggers

abstract: We present a cadence optimization strategy to unveil a large population of kilonovae using optical imaging alone. These transients are generated during binary neutron star and potentially neutron star-black hole mergers and are electromagnetic counterparts to gravitational-wave signals detectable in nearby events with Advanced LIGO, Advanced Virgo, and other interferometers that will come online in the near future. Discovering a large population of kilonovae will allow us to determine how heavy element production varies with the intrinsic parameters of the merger and across cosmic time. The rate of binary neutron star mergers is still uncertain, but only few (less than 15) events with associated kilonovae may be detectable per year within the horizon of next-generation ground-based interferometers. The rapid evolution (hours to days) at optical/infrared wavelengths, relatively low luminosity, and the low volumetric rate of kilonovae makes their discovery difficult, especially during blind surveys of the sky. We propose future large surveys to adopt a rolling cadence in which g-i observations are taken nightly for blocks of 10 consecutive nights. With the current baseline2018a cadence designed for the Large Synoptic Survey Telescope (LSST), less than 7.5 poorly-sampled kilonovae are expected to be detected in both the Wide Fast Deep (WFD) and Deep Drilling Fields (DDF) surveys per year, under optimistic assumptions on their rate, duration, and luminosity. We estimate the proposed strategy to return up to about 272 GW170817-like kilonovae throughout the LSST WFD survey, discovered independently from gravitational-wave triggers.


title: A reciprocal branching problem for automorphic representations and global Vogan packets

abstract: Let $G$ be a group and $H$ be a subgroup of $G$. The classical branching rule (or symmetry breaking) asks: For an irreducible representation $\pi$ of $G$, determine the occurrence of an irreducible representation $\sigma$ of $H$ in the restriction of $\pi$ to $H$. The reciprocal branching problem of this classical branching problem is to ask: For an irreducible representation $\sigma$ of $H$, find an irreducible representation $\pi$ of $G$ such that $\sigma$ occurs in the restriction of $\pi$ to $H$. For automorphic representations of classical groups, the branching problem has been addressed by the well-known global Gan-Gross-Prasad conjecture. In this paper, we investigate the reciprocal branching problem for automorphic representations of special orthogonal groups using the twisted automorphic descent method as developed in [JZ15]. The method may be applied to other classical groups as well.


title: Transfer learning for vision-based tactile sensing

abstract: Due to the complexity of modeling the elastic properties of materials, the use of machine learning algorithms is continuously increasing for tactile sensing applications. Recent advances in deep neural networks applied to computer vision make vision-based tactile sensors very appealing for their high-resolution and low cost. A soft optical tactile sensor that is scalable to large surfaces with arbitrary shape is discussed in this paper. A supervised learning algorithm trains a model that is able to reconstruct the normal force distribution on the sensor's surface, purely from the images recorded by an internal camera. In order to reduce the training times and the need for large datasets, a calibration procedure is proposed to transfer the acquired knowledge across multiple sensors while maintaining satisfactory performance.


title: Flat bands and dynamical localization of spin-orbit coupled Bose-Einstein condensates

abstract: Flat bands and dynamical localization of binary mixtures of Bose-Einstein condensates, with spin-orbit coupling subjected to a deep optical lattice which is shaking in time and to a periodic time modulation of the Zeeman field, are investigated. In contrast with usual dynamical localization in the absence of spin-orbit coupling, we find that to fully suppress the tunneling in the system the optical lattice shaking is not enough, and a proper tuning of the spin-orbit term, achievable via the Zeeman field modulation, is also required. This leads to a sequence of Zeeman parameter values where energy bands become flat, the tunneling in the system is suppressed, and the dynamical localization phenomenon occurs. Exact wave functions at the dynamical localization points show that the binary mixture localizes on a dimer with the two components occupying different sites. This type of localization occurs in exact form also for the ground state of the system at the dynamical localization points in the presence of nonlinearity and remains valid, although in approximate form, for a wide range of the Zeeman parameter around these points. The possibility of observing the above phenomena in real experiments is also briefly discussed.


title: Effective approach to lepton observables: the seesaw case

abstract: In the absence of direct evidence of new physics, any ultraviolet theory can be reduced to its specific set of low-energy effective operators. As a case study, we derive the effective field theory for the seesaw extension of the Standard Model, with sterile neutrinos of mass $M>m_W$. We systematically compute all Wilson coefficients generated at one loop. Hence, it becomes straightforward to (i) identify the seesaw parameters compatible with the smallness of neutrino masses; (ii) compute precision lepton observables, which may be sensitive to scales as large as $M\sim 10^3$ TeV; and (iii) establish sharp correlations among those observables. We find that the flavour-conserving Wilson coefficients set an upper bound on the flavour-violating ones. The low-energy limits on $\mu\to e$ and $\tau\to e,\mu$ transitions suppress flavour violation in $Z$ and Higgs decays, as well as electric dipole moments, far beyond the experimental reach. The precision measurements of $G_F$, $m_W$, and $Z$ partial decay widths set more stringent bounds than present and future limits on $\tau\to e,\mu$ transitions. We also present a general spurion analysis, to compare the seesaw with different models, thus assessing the discriminating potential of the effective approach.


title: Infrared massive gluon propagator from a BRST-invariant Gribov horizon in a family of covariant gauges

abstract: A distinctive feature observed in lattice simulations of confining nonabelian gauge theories, such as Quantum Chromodynamics, is the presence of a dynamical mass for the gauge field in the low energy regime of the theory. In the Gribov-Zwanziger framework in the Landau gauge, such mass is a consequence of the generation of the dimension two condensates $\langle A_\mu^aA_\mu^a\rangle$ and $\langle \bar\varphi_\mu^{ab}\varphi_\mu^{ab}-\bar\omega_\mu^{ab}\omega_\mu^{ab}\rangle$, where $A$ is the gluon field and the fields $\bar\varphi$, $\varphi$, $\bar\omega$, and $\omega$ are Zwanziger's auxiliary fields. In this work, we show that, in the recently developed BRST-invariant version of the Refined Gribov-Zwanziger theory, these condensates can be introduced in a BRST-invariant way for a family of $R_\xi$ gauges. Their values are explicitly computed to first order and turn out to be independent of the gauge parameters contained in the gauge-fixing condition, as expected from the BRST invariance of the formulation. This fact supports the possibility of a gauge-parameter independent nonzero infrared gluon mass, whose value is the same as the one in the Landau gauge.


title: Transverse Velocities with the Moving Lens Effect

abstract: Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.


title: LinKS: Discovering galaxy-scale strong lenses in the Kilo-Degree Survey using Convolutional Neural Networks

abstract: We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention.


title: Bound States of Pseudo-Dirac Dark Matter

abstract: We study the bound-state spectrum in a simple model of pseudo-Dirac dark matter, and examine how the rate of bound-state formation through radiative capture compares to Sommerfeld-enhanced annihilation. We use this model as an example to delineate the new features induced by the presence of a mass splitting between the dark matter and a nearly-degenerate partner, compared to the case where only a single dark-matter-like state is present. We provide a simple analytic prescription for estimating the spectrum of bound states in systems containing a mass splitting, which in turn allows characterization of the resonances due to near-zero-energy bound states, and validate this estimate both for pseudo-Dirac dark matter and for the more complex case of wino dark matter. We demonstrate that for pseudo-Dirac dark matter the capture rate into deeply bound states is, to a good approximation, simply related to the Sommerfeld enhancement factor.


title: Variational Saccading: Efficient Inference for Large Resolution Images

abstract: Image classification with deep neural networks is typically restricted to images of small dimensionality such as 224 x 244 in Resnet models [24]. This limitation excludes the 4000 x 3000 dimensional images that are taken by modern smartphone cameras and smart devices. In this work, we aim to mitigate the prohibitive inferential and memory costs of operating in such large dimensional spaces. To sample from the high-resolution original input distribution, we propose using a smaller proxy distribution to learn the co-ordinates that correspond to regions of interest in the high-dimensional space. We introduce a new principled variational lower bound that captures the relationship of the proxy distribution's posterior and the original image's co-ordinate space in a way that maximizes the conditional classification likelihood. We empirically demonstrate on one synthetic benchmark and one real world large resolution DSLR camera image dataset that our method produces comparable results with ~10x faster inference and lower memory consumption than a model that utilizes the entire original input distribution. Finally, we experiment with a more complex setting using mini-maps from Starcraft II [56] to infer the number of characters in a complex 3d-rendered scene. Even in such complicated scenes our model provides strong localization: a feature missing from traditional classification models.


title: The criticality of multiplicative processes

abstract: Keeping in view applications to numerical simulations of the evolution of a nuclear reactor core around criticality, we use a general mathematical framework for describing the evolutions of multiplicative processes (processes involving particle creation) both in particle generations and in time. This framework allows us to obtain, within a same formalism, two corresponding estimates of the multiplication factor which describes the growth of particle numbers at large times. We obtain the relative positions of both estimates with respect to each other and to criticality. These relations may show particularly useful when simulating in a realistic way the monitoring of nuclear cores in subcritical states, such as is the case for Accelerator Driven Systems (ADS). More generally, this study applies to various multiplicative processes which can be found in nature.


title: Ferroelectric polarization rotation in order-disorder-type LiNbO3 thin films

abstract: The direction of ferroelectric polarization is prescribed by the symmetry of the crystal structure. Therefore, rotation of the polarization direction is largely limited, despite the opportunity it offers in understanding important dielectric phenomena such as piezoelectric response near the morphotropic phase boundaries and practical applications such as ferroelectric memory. In this study, we report the observation of continuous rotation of ferroelectric polarization in order-disorder type LiNbO3 thin films. The spontaneous polarization could be tilted from an out-of-plane to an in-plane direction in the thin film by controlling the Li vacancy concentration within the hexagonal lattice framework. Partial inclusion of monoclinic-like phase is attributed to the breaking of macroscopic inversion symmetry along different directions and the emergence of ferroelectric polarization along the in-plane direction.


title: Use Dimensionality Reduction and SVM Methods to Increase the Penetration Rate of Computer Networks

abstract: In the world today computer networks have a very important position and most of the urban and national infrastructure as well as organizations are managed by computer networks, therefore, the security of these systems against the planned attacks is of great importance. Therefore, researchers have been trying to find these vulnerabilities so that after identifying ways to penetrate the system, they will provide system protection through preventive or countermeasures. SVM is one of the major algorithms for intrusion detection. In this research, we studied a variety of malware and methods of intrusion detection, provide an efficient method for detecting attacks and utilizing dimension reduction.Thus, we will be able to detect attacks by carefully combining these two algorithms and pre-processes that are performed before the two on the input data. The main question raised is how we can identify attacks on computer networks with the above-mentioned method. In anomalies diagnostic method, by identifying behavior as a normal behavior for the user, the host, or the whole system, any deviation from this behavior is considered as an abnormal behavior, which can be a potential occurrence of an attack. The network intrusion detection system is used by anomaly detection method that uses the SVM algorithm for classification and SVD to reduce the size. Steps of the proposed method include pre-processing of the data set, feature selection, support vector machine, and evaluation.The NSL-KDD data set has been used to teach and test the proposed model. In this study, we inferred the intrusion detection using the SVM algorithm for classification and SVD for diminishing dimensions with no classification algorithm.Also the KNN algorithm has been compared in situations with and without diminishing dimensions,the results have shown that the proposed method has a better performance than comparable methods.


title: Approximate Calculation of Tukey's Depth and Median With High-dimensional Data

abstract: We present a new fast approximate algorithm for Tukey (halfspace) depth level sets and its implementation-ABCDepth. Given a $d$-dimensional data set for any $d\geq 1$, the algorithm is based on a representation of level sets as intersections of balls in $\mathbb{R}^d$. Our approach does not need calculations of projections of sample points to directions. This novel idea enables calculations of approximate level sets in very high dimensions with complexity which is linear in $d$, which provides a great advantage over all other approximate algorithms. Using different versions of this algorithm we demonstrate approximate calculations of the deepest set of points ("Tukey median") and Tukey's depth of a sample point or out-of-sample point, all with a linear in $d$ complexity. An additional theoretical advantage of this approach is that the data points are not assumed to be in "general position". Examples with real and synthetic data show that the executing time of the algorithm in all mentioned versions in high dimensions is much smaller than the time of other implemented algorithms. Also, our algorithms can be used with thousands of multidimensional observations.


title: Elasticity tetrads, mixed axial-gravitational anomalies, and 3+1d quantum Hall effect

abstract: For two-dimensional topological insulators, the integer and intrinsic (without external magnetic field) quantum Hall effect is described by the gauge anomalous (2+1)-dimensional [2+1d] Chern-Simons (CS) response for the background gauge potential of the electromagnetic U(1) field. The Hall conductance is given by the quantized prefactor of the CS term, which is a momentum-space topological invariant. Here, we show that three-dimensional crystalline topological insulators with no other symmetries are described by a topological (3+1)-dimensional [3+1d] mixed CS term. In addition to the electromagnetic U(1) gauge field, this term contains elasticity tetrad fields $E^{\ a}_{\mu}({\bf r},t) = \partial_{\mu}X^a(\mathbf{r},t)$ which are gradients of crystalline U(1) phase fields $X^a(\mathbf{r},t)$ and describe the deformations of the crystal. For a crystal in three spatial dimensions $a=1,2,3$ and the mixed axial-gravitational response contains three parameters protected by crystalline symmetries: the weak momentum-space topological invariants. The response of the Hall conductance to the deformations of the crystal is quantized in terms of these invariants. In the presence of dislocations, the anomalous 3+1d CS term describes the Callan-Harvey anomaly inflow mechanism. The response can be extended to all odd spatial dimensions. The elasticity tetrads, being the gradients of the lattice U(1) fields, have canonical dimension of inverse length. Similarly, if such tetrad fields enter general relativity, the metric becomes dimensionful, but the physical parameters, such as Newton's constant, the cosmological constant, and masses of particles, become dimensionless.


title: Magnetic Braking and Damping of Differential Rotation in Massive Stars

abstract: Fragmentation of highly differentially rotating massive stars that undergo collapse has been suggested as a possible channel for binary black hole formation. Such a scenario could explain the formation of the new population of massive black holes detected by the LIGO/VIRGO gravitational wave laser interferometers. We probe that scenario by performing general relativistic magnetohydrodynamic simulations of differentially rotating massive stars supported by thermal radiation pressure plus a gas pressure perturbation. The stars are initially threaded by a dynamically weak, poloidal magnetic field confined to the stellar interior. We find that magnetic braking and turbulent viscous damping via magnetic winding and the magnetorotational instability in the bulk of the star redistribute angular momentum, damp differential rotation and induce the formation of a massive and nearly uniformly rotating inner core surrounded by a Keplerian envelope. The core + disk configuration evolves on a secular timescale and remains in quasi-stationary equilibrium until the termination of our simulations. Our results suggest that the high degree of differential rotation required for $m=2$ seed density perturbations to trigger gas fragmentation and binary black hole formation is likely to be suppressed during the normal lifetime of the star prior to evolving to the point of dynamical instability to collapse. Other cataclysmic events, such as stellar mergers leading to collapse, may therefore be necessary to reestablish sufficient differential rotation and density perturbations to drive nonaxisymmetric modes leading to binary black hole formation.


title: Adiabatic two-qubit gates in capacitively coupled quantum dot hybrid qubits

abstract: The ability to tune qubits to flat points in their energy dispersions ("sweet spots") is an important tool for mitigating the effects of charge noise and dephasing in solid-state devices. However, the number of derivatives that must be simultaneously set to zero grows exponentially with the number of coupled qubits, making the task untenable for as few as two qubits. This is a particular problem for adiabatic gates, due to their slower speeds. Here, we propose an adiabatic two-qubit gate for quantum dot hybrid qubits, based on the tunable, electrostatic coupling between distinct charge configurations. We confirm the absence of a conventional sweet spot, but show that controlled-Z (CZ) gates can nonetheless be optimized to have fidelities of $\sim$99% for a typical level of quasistatic charge noise ($\sigma_\varepsilon$$\simeq$1 $\mu$eV). We then develop the concept of a dynamical sweet spot (DSS), for which the time-averaged energy derivatives are set to zero, and identify a simple pulse sequence that achieves an approximate DSS for a CZ gate, with a 5$\times$ improvement in the fidelity. We observe that the results depend on the number of tunable parameters in the pulse sequence, and speculate that a more elaborate sequence could potentially attain a true DSS.


title: Phenomenology of TeV-scale scalar Leptoquarks in the EFT

abstract: We examine new aspects of leptoquark (LQ) phenomenology using effective field theory (EFT). We construct a complete set of leading effective operators involving SU(2) singlets scalar LQ and the SM fields up to dimension six. We show that, while the renormalizable LQ-lepton-quark interaction Lagrangian can address the persistent hints for physics beyond the Standard Model in the B-decays $\bar B \to D^{(*)} \tau \bar\nu$, $\bar B \to \bar K \ell^+ \ell^-$ and in the measured anomalous magnetic moment of the muon, the LQ higher dimensional effective operators may lead to new interesting effects associated with lepton number violation. These include the generation of one-loop sub-eV Majorana neutrino masses, mediation of neutrinoless double-$\beta$ decay and novel LQ collider signals. For the latter, we focus on 3rd generation LQ ($\phi_3$) in a framework with an approximate $Z_3$ generation symmetry, and show that one class of the dimension five LQ operators may give rise to a striking asymmetric same-charge $\phi_3 \phi_3$ pair-production signal, which leads to low background same-sign leptons signals at the LHC. For example, with $M_{\phi_3} \sim 1$ TeV and a new physics scale of $\Lambda \sim 5$ TeV, we expect about $5000$ positively charged $\tau^+ \tau^+$ events via $pp \to \phi_3 \phi_3 \to \tau^+ \tau^+ + 2 \cdot j_b$ ($j_b$=b-jet) at the 13 TeV LHC with an integrated luminosity of 300 fb$^{-1}$. It is interesting to note that, in the LQ EFT framework, the expected same-sign lepton signals have a rate which is several times larger than the QCD LQ-mediated opposite-sign leptons signals, $gg,q \bar q \to \phi_3 \phi_3^* \to \ell^+ \ell^- +X$. We also consider the same-sign charged lepton signals in the LQ EFT framework at higher energy hadron colliders such as a 27 TeV HE-LHC and a 100 TeV FCC-hh.


title: The non-linear X-ray/UV relation in active galactic nuclei: contribution of instrumental effects on the X-ray variability

abstract: We have recently demonstrated that the non-linear relation between ultraviolet and X-ray luminosity in quasars is very tight (with an intrinsic dispersion of ~0.2 dex), once contaminants (e.g. dust reddening, X-ray absorption), variability, and differences in the active galactic nuclei (AGN) physical properties are taken into account. This relation has thus the great potential to advance our understanding in both supermassive black hole accretion physics and observational cosmology, by targeting a single class of objects. Here we focus on the various contributions to the observed X-ray variability in a homogenous sample of 791 quasars selected from SDSS-DR7 with X-ray data from the 3XMM-DR7 source catalogue. The 250 quasars in this cleaned data set with at least two X-ray observations typically vary with a standard deviation of fractional variation of 15-30% on timescales of weeks/years. Yet, when the count rates are computed at progressively smaller off-axis values, the same quantity is reduced to roughly 10-25%. This suggests that, when estimating variability indicators, part of the quoted variability amplitude could be due to instrumental/calibration issues rather than true variations in the quasar emission.


title: High-energy monitoring of NGC 4593 II. Broadband spectral analysis: testing the two-corona model

abstract: It is widely believed that the primary X-ray emission of AGN is due to the Comptonisation of optical-UV photons from a hot electron corona, while the origin of the 'soft-excess' is still uncertain and matter of debate. A second Comptonisation component, called warm corona, was therefore proposed to account for the soft-excess, and found in agreement with the optical-UV to X-ray emission of a sample of Seyfert galaxies. In this context, we exploit the broadband XMM-Newton and NuSTAR simultaneous observations of the Seyfert galaxy NGC 4593 to further test the so called "two corona model". The NGC 4593 spectra are well reproduced by the model, from the optical/UV to the hard X-rays. Moreover, the data reveal a significant correlation between the hot and the warm corona parameters during our monitoring campaign.


title: Modelling Metropolitan-area Ambulance Mobility under Blue Light Conditions

abstract: Actions taken immediately following a life-threatening personal health incident are critical for the survival of the sufferer. The timely arrival of specialist ambulance crew in particular often makes the difference between life and death. As a consequence, it is critical that emergency ambulance services achieve short response times. This objective sets a considerable challenge to ambulance services worldwide, especially in metropolitan areas where the density of incident occurrence and traffic congestion are high. Using London as a case study, in this paper we consider the advantages and limitations of data-driven methods for ambulance routing and navigation. Our long-term aim is to enable considerable improvements to their operational efficiency through the automated generation of more effective response strategies and tactics. A key ingredient of our approach is to use a large historical dataset of incidents and ambulance location traces to model route selection and arrival times. Working on the London road network graph modified to reflect the differences between emergency and civilian vehicle traffic, we develop a methodology for the precise estimation of expected ambulance speed at the individual road segment level. We demonstrate how a model that exploits this information achieves best predictive performance by implicitly capturing route-specific persistent patterns in changing traffic conditions. We then present a predictive method that achieves a high route similarity score while minimising journey duration error. This is achieved through the combination of a technique that correctly predicts routes selected by the current navigation system of the London Ambulance Service and our best performing speed estimation model. This hybrid approach outperforms alternative mobility models.


title: Bounds and Prospects for Stable Multiply Charged Particles at the LHC

abstract: Colored and colorless particles that are stable on collider scales and carry exotic electric charges, so-called multiply-charged heavy stable particles (MCHSPs), exist in extensions of the Standard Model, and can include the top partner(s) in solutions of the hierarchy problem. To obtain bounds on color-triplets and color-singlets of charges up to |Q|=8, we recast searches for signatures of two production channels: the "open" channel - where the particles are pair-produced above threshold, and are detectable in dedicated LHC searches for stable multiply charged leptons, and the "closed" channel - where a particle-antiparticle pair is produced as a bound state, detectable in searches for a diphoton resonance. We recast the open lepton searches by incorporating the relevant strong-interaction effects for color-triplets. In both open and closed production, we provide a careful assessment of photon-induced processes using the accurate LUXqed PDF, resulting in substantially weaker bounds than previously claimed in the literature for the colorless case. Our bounds for colored MCHSPs are shown for the first time, as the LHC experiments have not searched for them directly. Generally, we obtain nearly charge-independent lower mass limits of around 970 GeV (color-triplet scalar), 1200 GeV (color-triplet fermion), and 880 - 900 GeV (color-singlet fermion) from open production, and strongly charge-dependent limits from closed production. In all cases there is a crossover between dominance by open and closed searches at some charge. We provide prospective bounds for $\sqrt{s} = 13$ TeV LHC searches at integrated luminosities of 39.5 fb${}^{-1}$, 100 fb${}^{-1}$, and 300 fb${}^{-1}$. Moreover, we show that a joint observation in the open and the closed channels allows to determine the mass, spin, color, and electric charge of the particle.


title: A hybrid machine-learning algorithm for designing quantum experiments

abstract: We introduce a hybrid machine-learning algorithm for designing quantum optics experiments that produce specific quantum states. Our algorithm successfully found experimental schemes to produce all 5 states we asked it to, including Schr\"odinger cat states and cubic phase states, all to a fidelity of over $96\%$. Here we specifically focus on designing realistic experiments, and hence all of the algorithm's designs only contain experimental elements that are available with current technology. The core of our algorithm is a genetic algorithm that searches for optimal arrangements of the experimental elements, but to speed up the initial search we incorporate a neural network that classifies quantum states. The latter is of independent interest, as it quickly learned to accurately classify quantum states given their photon-number distributions.


title: Quantum Fine-Tuning in Stringy Quintessence Models

abstract: We investigate the extent to which quintessence models for dark energy are fine-tuned in the context of recent swampland conjectures. In particular, the issue is whether there is a double fine-tuning in which both $V$ and $|\nabla V|$ are fine-tuned, or whether there is only a single fine-tuning due to the relation $|\nabla V|\sim V/M_{pl}$ arising naturally. We find that indeed this relation arises naturally in simple string compactifications for some scalars, such as the dilaton and volume modulus, when treated classically. However, we find that quantum effects can spoil this natural relation, unless the scalar is conformally coupled to the matter sector. Furthermore, it is well known that such conformal couplings are generically ruled out by fifth force tests. To avoid these fifth forces, an interesting proposal is to assume the scalar (quintessence) only couples to the hidden sector. However, we then find quantum corrections to $V$ from visible sector Standard Model particles generically spoil the relation. A possible way out of all these problems is to have the scalar conformally coupled to a dark sector that is an exact copy of the Standard Model. This ensures the relation $|\nabla V|\sim V/M_{pl}$ is maintained naturally even when matter particles run in the loop. However, we find that quantum corrections from quintessons or gravitons in the loop spoil the relation if the effective theory has a cutoff greater than $\sim 0.1$ GeV.


title: Gauss-Manin Lie algebra of mirror elliptic K3 surfaces

abstract: We study mirror symmetry of families of elliptic K3 surfaces with elliptic fibers of type $E_6,~E_7$ and $E_8$. We consider a moduli space $\mathsf{T}$ of the mirror K3 surfaces enhanced with the choice of differential forms. We show that coordinates on $\mathsf{T}$ are given by the ring of quasi modular forms in two variables, with modular groups adapted to the fiber type. We furthermore introduce an algebraic group $\mathsf{G}$ which acts on $\mathsf{T}$ from the right and construct its Lie algebra $\mathrm{Lie}(\mathsf{G})$. We prove that the extended Lie algebra generated by $\mathrm{Lie}(\mathsf{G})$ together with modular vector fields on $\mathsf{T}$ is isomorphic to $\mathrm{sl}_2(\mathbb{C})\oplus\mathrm{sl}_2(\mathbb{C})$.


title: Removal of Parameter Adjustment of Frangi Filters in Case of Coronary Angiograms

abstract: Frangi Filters are one of the widely used filters for enhancing vessels in medical images. Since they were first proposed, the threshold of the vesselness function of Frangi Filters is to be arranged for each individual application. These thresholds are changed manually for individual fluoroscope, for enhancing coronary angiogram images. Hence it is felt, there is a need of mitigating the tuning procedure of threshold values for every fluoroscope. The current papers approach has been devised in order to treat the coronary angiogram images uniformly, irrespective of the fluoroscopes through which they were obtained and the patient demographics for further stenosis detection. This problem to the best of our knowledge has not been addressed yet. In the approach, before feeding the image to Frangi Filters, non uniform illumination of the input image is removed using homomorphic filters and the image is enhanced using Non Subsampled Contourlet Transform (NSCT). The experiment was conducted on the data that has been accumulated from various hospitals in India and the results obtained verifies dependency removal of parameters without compromising the results obtained by Frangi filters.


title: Towards a priori uncertainty quantification in coarse-grained molecular dynamics: Generalized multipole potentials

abstract: In computational materials science, coarse-graining approaches often lack a priori uncertainty quantification (UQ) tools that estimate the accuracy of a reduced-order model before it is calibrated or deployed. This is especially the case in coarse-grained (CG) molecular dynamics (MD), where "bottom-up" methods need to run expensive atomistic simulations as part of the calibration process. As a result, scientists have been slow to adopt CG techniques in many settings because they do not know in advance whether the cost of developing the CG model is justified. To address this problem, we present an analytical method of coarse-graining rigid-body systems that yields corresponding intermolecular potentials with controllable levels of accuracy relative to their atomistic counterparts. Critically, this analysis: (i) provides a mathematical foundation for assessing the quality of a CG force field without running simulations, and (ii) provides a tool for understanding how atomistic systems can be viewed as appropriate limits of reduced-order models. Simulated results confirm the validity of this approach at the trajectory level and point to issues that must be addressed in coarse-graining fully non-rigid systems.


title: METCC: METric learning for Confounder Control Making distance matter in high dimensional biological analysis

abstract: High-dimensional data acquired from biological experiments such as next generation sequencing are subject to a number of confounding effects. These effects include both technical effects, such as variation across batches from instrument noise or sample processing, or institution-specific differences in sample acquisition and physical handling, as well as biological effects arising from true but irrelevant differences in the biology of each sample, such as age biases in diseases. Prior work has used linear methods to adjust for such batch effects. Here, we apply contrastive metric learning by a non-linear triplet network to optimize the ability to distinguish biologically distinct sample classes in the presence of irrelevant technical and biological variation. Using whole-genome cell-free DNA data from 817 patients, we demonstrate that our approach, METric learning for Confounder Control (METCC), is able to match or exceed the classification performance achieved using a best-in-class linear method (HCP) or no normalization. Critically, results from METCC appear less confounded by irrelevant technical variables like institution and batch than those from other methods even without access to high quality metadata information required by many existing techniques; offering hope for improved generalization.


title: ALMA High Angular Resolution Polarization Study; An Extremely Young Class 0 Source, OMC-3/MMS 6

abstract: Using the $\approx$15km ALMA long baselines, we imaged the Stokes $I$ emission and linearly polarized intensity ($PI$) in the 1.1-mm continuum band of a very young intermediate-mass protostellar source, MMS 6, in the Orion Molecular Cloud-3. The achieved angular resolution, $0''.02{\times}0''.03$ ($\approx$10 AU), shows for the first time a wealth of data on the dust emission polarization in the central 200 AU of a protostar. The $PI$ peak is offset to the south-west (SW) by $\approx$20 AU with respect to the Stokes $I$ peak. Its polarization degree is 11 % with its $E$-vector orientation of P.A.${\approx}135^{\circ}$. A partial ring-like structure with a radius of $\approx$80 AU is detected in $PI$ but not in the Stokes $I$. NW (north-west) and SE (south-east) parts of the ring are bright with a high polarization degree of $\gtrsim$10 %, and their $E$-vector orientations are roughly orthogonal to those observed near the center. We also detected arm-like polarized structures, extending to 1000 AU scale to the north, with the $E$-vectors aligned along the minor axis of the structures. We explored possible origins of the polarized emission comparing with magnetohydrodynamical (MHD) simulations of the toroidal wrapping of the magnetic field. The simulations are consistent with the $PI$ emission in the ring-like and the extended arm-like structures observed with ALMA. However, the current simulations do not completely reproduce observed polarization characteristics in the central 50 AU. Although the self-scattering model can explain the polarization pattern and positional offset between the Stokes $I$ and $PI$, this model is not able to reproduce the observed high degree of polarization.


title: Deep-RBF Networks Revisited: Robust Classification with Rejection

abstract: One of the main drawbacks of deep neural networks, like many other classifiers, is their vulnerability to adversarial attacks. An important reason for their vulnerability is assigning high confidence to regions with few or even no feature points. By feature points, we mean a nonlinear transformation of the input space extracting a meaningful representation of the input data. On the other hand, deep-RBF networks assign high confidence only to the regions containing enough feature points, but they have been discounted due to the widely-held belief that they have the vanishing gradient problem. In this paper, we revisit the deep-RBF networks by first giving a general formulation for them, and then proposing a family of cost functions thereof inspired by metric learning. In the proposed deep-RBF learning algorithm, the vanishing gradient problem does not occur. We make these networks robust to adversarial attack by adding the reject option to their output layer. Through several experiments on the MNIST dataset, we demonstrate that our proposed method not only achieves significant classification accuracy but is also very resistant to various adversarial attacks.


title: Numerical study of dynamical friction with thermal effects -- I. Comparison to linear theory

abstract: We investigate by means of numerical simulations with nested Cartesian meshes the force exerted on a massive and luminous perturber moving at constant speed through a homogeneous and opaque gas, taking into account thermal diffusion in the gas and the radiative feedback from the perturber. The force arising from the release of energy into the ambient medium by the perturber, or heating force, is directed along the direction of motion and induces an acceleration of the perturber. Its value is compared to analytic estimates in the low and high Mach number limits, and found to match those accurately. In addition, the drag exerted on a non-luminous perturber significantly departs from the adiabatic expression when thermal diffusion is taken into account. In the limit of a vanishing velocity, this drag tends to a finite value which we determine using linear perturbation theory and corroborate with numerical simulations. The drag on a non-luminous perturber in a non-adiabatic gas therefore behaves like dry or solid friction. We work out the luminosity threshold to get a net acceleration of the perturber and find it to be generally much smaller than the luminosity of accreting low-mass planetary embryos embedded in a gaseous protoplanetary disc at a few astronomical units. We also present in some detail our implementation of nested meshes, which runs in parallel over several \emph{Graphics Processing Units} (GPUs).


title: A stable added-mass partitioned (AMP) algorithm for elastic solids and incompressible flow: model problem analysis

abstract: A stable added-mass partitioned (AMP) algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and compressible elastic-solids. The AMP scheme remains stable and second-order accurate even when added-mass and added-damping effects are large. The fluid is updated with an implicit-explicit (IMEX) fractional-step scheme whereby the velocity is advanced in one step, treating the viscous terms implicitly, and the pressure is computed in a second step. The AMP interface conditions for the fluid arise from the outgoing characteristic variables in the solid and are partitioned into a Robin (mixed) interface condition for the pressure, and interface conditions for the velocity. The latter conditions include an impedance-weighted average between fluid and solid velocities using a fluid impedance of a special form. A similar impedance-weighted average is used to define interface values for the solid. The fluid impedance is defined using material and discretization parameters and follows from a careful analysis of the discretization of the governing equations and coupling conditions near the interface. A normal mode analysis is performed to show that the AMP scheme is stable for a few carefully-selected model problems. Two extensions of the analysis in Banks et al. are considered, including a first-order accurate discretization of a viscous model problem and a second-order accurate discretization of an inviscid model problem. The AMP algorithm is shown to be stable for any ratio of solid and fluid densities, including when added-mass effects are large. The algorithm is verified for accuracy and stability for a set of new exact benchmark solutions where finite interface deformations are permitted. The AMP scheme is found to be stable and second-order accurate even for very difficult cases of very light solids.


title: A class of weighted Hardy inequalities and applications to evolution problems

abstract: \begin{abstract} We state the following weighted Hardy inequality \begin{equation*} c_{o, \mu}\int_{{\R}^N}\frac{\varphi^2 }{|x|^2}\, d\mu\le \int_{{\R}^N} |\nabla\varphi|^2 \, d\mu + K \int_{\R^N}\varphi^2 \, d\mu \quad \forall\, \varphi \in H_\mu^1 %\qquad c\le c_\mu, \end{equation*} in the context of the study of the Kolmogorov operators \begin{equation*} Lu=\Delta u+\frac{\nabla \mu}{\mu}\cdot\nabla u \end{equation*} perturbed by inverse square potentials and of the related evolution problems. The function $\mu$ in the drift term is a probability density on $\R^N$. We prove the optimality of the constant $c_{o, \mu}$ and state existence and nonexistence results following the Cabr\'e-Martel's approach \cite{CabreMartel} extended to Kolmogorov operators. \end{abstract}


title: The Soft-Excess in Mrk 509: Warm Corona or Relativistic Reflection?

abstract: We present the analysis of the first NuSTAR observations ($\sim 220$ ks), simultaneous with the last SUZAKU observations ($\sim 50$ ks), of the active galactic nucleus of the bright Seyfert 1 galaxy Mrk 509. The time-averaged spectrum in the $1-79$ keV X-ray band is dominated by a power-law continuum ($\Gamma\sim 1.8-1.9$), a strong soft excess around 1 keV, and signatures of X-ray reflection in the form of Fe K emission ($\sim 6.4$ keV), an Fe K absorption edge ($\sim 7.1$ keV), and a Compton hump due to electron scattering ($\sim 20-30$ keV). We show that these data can be described by two very different prescriptions for the soft excess: a warm ($kT\sim 0.5-1$ keV) and optically thick ($\tau\sim10-20$) Comptonizing corona, or a relativistically blurred ionized reflection spectrum from the inner regions of the accretion disk. While these two scenarios cannot be distinguished based on their fit statistics, we argue that the parameters required by the warm corona model are physically incompatible with the conditions of standard coronae. Detailed photoionization calculations show that even in the most favorable conditions, the warm corona should produce strong absorption in the observed spectrum. On the other hand, while the relativistic reflection model provides a satisfactory description of the data, it also requires extreme parameters, such as maximum black hole spin, a very low and compact hot corona, and a very high density for the inner accretion disk. Deeper observations of this source are thus necessary to confirm the presence of relativistic reflection, and to further understand the nature of its soft excess.


title: Counting independent sets in graphs with bounded bipartite pathwidth

abstract: We show that a simple Markov chain, the Glauber dynamics, can efficiently sample independent sets almost uniformly at random in polynomial time for graphs in a certain class. The class is determined by boundedness of a new graph parameter called bipartite pathwidth. This result, which we prove for the more general hardcore distribution with fugacity $\lambda$, can be viewed as a strong generalisation of Jerrum and Sinclair's work on approximately counting matchings, that is, independent sets in line graphs. The class of graphs with bounded bipartite pathwidth includes claw-free graphs, which generalise line graphs. We consider two further generalisations of claw-free graphs and prove that these classes have bounded bipartite pathwidth. We also show how to extend all our results to polynomially-bounded vertex weights.