Showing 1–50 of 129 results for author: Larsson, S

Towards Negotiative Dialogue for the Talkamatic Dialogue Manager

Authors: Staffan Larsson, Alexander Berman, David Hjelm

Abstract: The paper describes a number of dialogue phenomena associated with negotiative dialogue, as implemented in a development version of the Talkamatic Dialogue Manager (TDM). This implementation is an initial step towards full coverage of general features of negotiative dialogue in TDM. The paper describes a number of dialogue phenomena associated with negotiative dialogue, as implemented in a development version of the Talkamatic Dialogue Manager (TDM). This implementation is an initial step towards full coverage of general features of negotiative dialogue in TDM. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Topic and genre in dialogue

Authors: Amandine Decker, Ellen Breitholtz, Christine Howes, Staffan Larsson

Abstract: In this paper we argue that topic plays a fundamental role in conversations, and that the concept is needed in addition to that of genre to define interactions. In particular, the concepts of genre and topic need to be separated and orthogonally defined. This would enable modular, reliable and controllable flexible-domain dialogue systems. In this paper we argue that topic plays a fundamental role in conversations, and that the concept is needed in addition to that of genre to define interactions. In particular, the concepts of genre and topic need to be separated and orthogonally defined. This would enable modular, reliable and controllable flexible-domain dialogue systems. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Journal ref: SEMDIAL 2023, Aug 2023, Maribor, Slovenia. pp.143-145

On the importance of AI research beyond disciplines

Authors: Virginia Dignum, Donal Casey, Teresa Cerratto-Pargman, Frank Dignum, Valentina Fantasia, Bodil Formark, Björn Hammarfelt, Gunnar Holmberg, André Holzapfel, Stefan Larsson, Amanda Lagerkvist, Nicolette Lakemond, Helena Lindgren, Fabian Lorig, Ana Marusic, Lina Rahm, Yulia Razmetaeva, Sverker Sikström, Kıvanç Tatar, Jason Tucker

Abstract: As the impact of AI on various scientific fields is increasing, it is crucial to embrace interdisciplinary knowledge to understand the impact of technology on society. The goal is to foster a research environment beyond disciplines that values diversity and creates, critiques and develops new conceptual and theoretical frameworks. Even though research beyond disciplines is essential for understand… ▽ More As the impact of AI on various scientific fields is increasing, it is crucial to embrace interdisciplinary knowledge to understand the impact of technology on society. The goal is to foster a research environment beyond disciplines that values diversity and creates, critiques and develops new conceptual and theoretical frameworks. Even though research beyond disciplines is essential for understanding complex societal issues and creating positive impact it is notoriously difficult to evaluate and is often not recognized by current academic career progression. The motivation for this paper is to engage in broad discussion across disciplines and identify guiding principles fir AI research beyond disciplines in a structured and inclusive way, revealing new perspectives and contributing to societal and human wellbeing and sustainability. △ Less

Submitted 13 February, 2023; originally announced February 2023.

The Fourth Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope -- Data Release 3

Authors: The Fermi-LAT collaboration, :, Marco Ajello, Luca Baldini, Jean Ballet, Denis Bastieri, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Elisabetta Bissaldi, Raffaella Bonino, Ari Brill, Philippe Bruel, Sara Buson, Regina Caputo, Patrizia Caraveo, Teddy Cheung, Graziano Chiaro, Nicolo Cibrario, Stefano Ciprini, Milena Crnogorcevic, Sara Cutini, Filippo D'Ammando, Salvatore De Gaetano, Niccolo Di Lalla , et al. (79 additional authors not shown)

Abstract: An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all source… ▽ More An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data. △ Less

Submitted 6 October, 2022; v1 submitted 24 September, 2022; originally announced September 2022.

Comments: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal Supplement Series. Fits files are available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/4LACDR3/ or alternatively at https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-h.fits and https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-l.fits

Search for new cosmic-ray acceleration sites within the 4FGL catalog Galactic plane sources

Authors: Fermi-LAT Collaboration, S. Abdollahi, F. Acero, M. Ackermann, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, B. Berenji, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, P. Bruel, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, D. Castro, G. Chiaro, N. Cibrario, S. Ciprini, J. Coronado-Blázquez, M. Crnogorcevic , et al. (95 additional authors not shown)

Abstract: Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Superno… ▽ More Cosmic rays are mostly composed of protons accelerated to relativistic speeds. When those protons encounter interstellar material, they produce neutral pions which in turn decay into gamma rays. This offers a compelling way to identify the acceleration sites of protons. A characteristic hadronic spectrum, with a low-energy break around 200 MeV, was detected in the gamma-ray spectra of four Supernova Remnants (SNRs), IC 443, W44, W49B and W51C, with the Fermi Large Area Telescope. This detection provided direct evidence that cosmic-ray protons are (re-)accelerated in SNRs. Here, we present a comprehensive search for low-energy spectral breaks among 311 4FGL catalog sources located within 5 degrees from the Galactic plane. Using 8 years of data from the Fermi Large Area Telescope between 50 MeV and 1 GeV, we find and present the spectral characteristics of 56 sources with a spectral break confirmed by a thorough study of systematic uncertainty. Our population of sources includes 13 SNRs for which the proton-proton interaction is enhanced by the dense target material; the high-mass gamma-ray binary LS~I +61 303; the colliding wind binary eta Carinae; and the Cygnus star-forming region. This analysis better constrains the origin of the gamma-ray emission and enlarges our view to potential new cosmic-ray acceleration sites. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: Accepted for publication in The Astrophysical Journal

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

Authors: M. Ajello, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, A. Berretta, B. Bhattacharyya, E. Bissaldi, R. D. Blandford, E. Bloom, R. Bonino, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, N. Cibrario, S. Ciprini, C. J. Clark, I. Cognard, J. Coronado-Blázquez , et al. (107 additional authors not shown)

Abstract: After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to… ▽ More After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

An energy-based deep splitting method for the nonlinear filtering problem

Authors: Kasper Bågmark, Adam Andersson, Stig Larsson

Abstract: The purpose of this paper is to explore the use of deep learning for the solution of the nonlinear filtering problem. This is achieved by solving the Zakai equation by a deep splitting method, previously developed for approximate solution of (stochastic) partial differential equations. This is combined with an energy-based model for the approximation of functions by a deep neural network. This res… ▽ More The purpose of this paper is to explore the use of deep learning for the solution of the nonlinear filtering problem. This is achieved by solving the Zakai equation by a deep splitting method, previously developed for approximate solution of (stochastic) partial differential equations. This is combined with an energy-based model for the approximation of functions by a deep neural network. This results in a computationally fast filter that takes observations as input and that does not require re-training when new observations are received. The method is tested on four examples, two linear in one and twenty dimensions and two nonlinear in one dimension. The method shows promising performance when benchmarked against the Kalman filter and the bootstrap particle filter. △ Less

Submitted 21 December, 2022; v1 submitted 31 March, 2022; originally announced March 2022.

Comments: 22 pages, 7 figures

MSC Class: 60G35; 62F15; 62G07; 62M20; 65C30; 65M75; 68T07

Incremental Fermi Large Area Telescope Fourth Source Catalog

Authors: Fermi-LAT collaboration, :, Soheila Abdollahi, Fabio Acero, Luca Baldini, Jean Ballet, Denis Bastieri, Ronaldo Bellazzini, Bijan Berenji, Alessandra Berretta, Elisabetta Bissaldi, Roger D. Blandford, Elliott Bloom, Raffaella Bonino, Ari Brill, Richard J. Britto, Philippe Bruel, Toby H. Burnett, Sara Buson, Rob A. Cameron, Regina Caputo, Patrizia A. Caraveo, Daniel Castro, Sylvain Chaty, Teddy C. Cheung , et al. (116 additional authors not shown)

Abstract: We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral param… ▽ More We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources. △ Less

Submitted 10 May, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: accepted in ApJS; follow-up paper to 1902.10045

Journal ref: ApJS 260, 53 (2022)

Achieving a Data-driven Risk Assessment Methodology for Ethical AI

Authors: Anna Felländer, Jonathan Rebane, Stefan Larsson, Mattias Wiggberg, Fredrik Heintz

Abstract: The AI landscape demands a broad set of legal, ethical, and societal considerations to be accounted for in order to develop ethical AI (eAI) solutions which sustain human values and rights. Currently, a variety of guidelines and a handful of niche tools exist to account for and tackle individual challenges. However, it is also well established that many organizations face practical challenges in n… ▽ More The AI landscape demands a broad set of legal, ethical, and societal considerations to be accounted for in order to develop ethical AI (eAI) solutions which sustain human values and rights. Currently, a variety of guidelines and a handful of niche tools exist to account for and tackle individual challenges. However, it is also well established that many organizations face practical challenges in navigating these considerations from a risk management perspective. Therefore, new methodologies are needed to provide a well-vetted and real-world applicable structure and path through the checks and balances needed for ethically assessing and guiding the development of AI. In this paper we show that a multidisciplinary research approach, spanning cross-sectional viewpoints, is the foundation of a pragmatic definition of ethical and societal risks faced by organizations using AI. Equally important is the findings of cross-structural governance for implementing eAI successfully. Based on evidence acquired from our multidisciplinary research investigation, we propose a novel data-driven risk assessment methodology, entitled DRESS-eAI. In addition, through the evaluation of our methodological implementation, we demonstrate its state-of-the-art relevance as a tool for sustaining human values in the data-driven AI era. △ Less

Submitted 29 November, 2021; originally announced December 2021.

Comments: 29 pages, 5 figures

We went to look for meaning and all we got were these lousy representations: aspects of meaning representation for computational semantics

Authors: Simon Dobnik, Robin Cooper, Adam Ek, Bill Noble, Staffan Larsson, Nikolai Ilinykh, Vladislav Maraev, Vidya Somashekarappa

Abstract: In this paper we examine different meaning representations that are commonly used in different natural language applications today and discuss their limits, both in terms of the aspects of the natural language meaning they are modelling and in terms of the aspects of the application for which they are used. In this paper we examine different meaning representations that are commonly used in different natural language applications today and discuss their limits, both in terms of the aspects of the natural language meaning they are modelling and in terms of the aspects of the application for which they are used. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 14 pages

Fermi Large Area Telescope Performance After 10 Years Of Operation

Authors: The Fermi LAT Collaboration, M. Ajello, W. B. Atwood, M. Axelsson, R. Bagagli, M. Bagni, L. Baldini, D. Bastieri, F. Bellardi, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, J. Bregeon, A. Brez, P. Bruel, R. Buehler, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, M. Ceccanti, S. Chen, C. C. Cheung, S. Ciprini , et al. (104 additional authors not shown)

Abstract: The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase… ▽ More The Large Area Telescope (LAT), the primary instrument for the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy gamma-ray telescope, covering the energy range from 30 MeV to more than 300 GeV. We describe the performance of the instrument at the 10-year milestone. LAT performance remains well within the specifications defined during the planning phase, validating the design choices and supporting the compelling case to extend the duration of the Fermi mission. The details provided here will be useful when designing the next generation of high-energy gamma-ray observatories. △ Less

Submitted 6 September, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 60 pages, 28 figures. Published in ApJS

Journal ref: ApJS 256 (2021) 12

Catalog of Long-Term Transient Sources in the First 10 Years of Fermi-LAT Data

Authors: L. Baldini, J. Ballet, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, A. Berretta, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, P. Bruel, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, S. Chen, G. Chiaro, D. Ciangottini, S. Ciprini, P. Cristarella Orestano, M. Crnogorcevic, S. Cutini, F. D'Ammando, P. de la Torre Luque , et al. (90 additional authors not shown)

Abstract: We present the first Fermi Large Area Telescope (LAT) catalog of long-term $γ$-ray transient sources (1FLT). This comprises sources that were detected on monthly time intervals during the first decade of Fermi-LAT operations. The monthly time scale allows us to identify transient and variable sources that were not yet reported in other Fermi-LAT catalogs. The monthly datasets were analyzed using a… ▽ More We present the first Fermi Large Area Telescope (LAT) catalog of long-term $γ$-ray transient sources (1FLT). This comprises sources that were detected on monthly time intervals during the first decade of Fermi-LAT operations. The monthly time scale allows us to identify transient and variable sources that were not yet reported in other Fermi-LAT catalogs. The monthly datasets were analyzed using a wavelet-based source detection algorithm that provided the candidate new transient sources. The search was limited to the extragalactic regions of the sky to avoid the dominance of the Galactic diffuse emission at low Galactic latitudes. The transient candidates were then analyzed using the standard Fermi-LAT Maximum Likelihood analysis method. All sources detected with a statistical significance above 4$σ$ in at least one monthly bin were listed in the final catalog. The 1FLT catalog contains 142 transient $γ$-ray sources that are not included in the 4FGL-DR2 catalog. Many of these sources (102) have been confidently associated with Active Galactic Nuclei (AGN): 24 are associated with Flat Spectrum Radio Quasars; 1 with a BL Lac object; 70 with Blazars of Uncertain Type; 3 with Radio Galaxies; 1 with a Compact Steep Spectrum radio source; 1 with a Steep Spectrum Radio Quasar; 2 with AGN of other types. The remaining 40 sources have no candidate counterparts at other wavelengths. The median $γ$-ray spectral index of the 1FLT-AGN sources is softer than that reported in the latest Fermi-LAT AGN general catalog. This result is consistent with the hypothesis that detection of the softest $γ$-ray emitters is less efficient when the data are integrated over year-long intervals. △ Less

Submitted 31 May, 2021; originally announced June 2021.

Comments: 41 pages, 17 figures, 7 tables; Accepted by ApJS on 24 May 2021; Contact Authors: I. Mereu, S. Cutini, E. Cavazzuti, G. Tosti

The First Fermi-LAT Solar Flare Catalog

Authors: M. Ajello, L. Baldini, D. Bastieri, R. Bellazzini, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, P. Bruel, S. Buson, R. A. Cameron, R. Caputo, E. Cavazzuti, C. C. Cheung, G. Chiaro, D. Costantin, S. Cutini, F. D'Ammando, F. de Palma, R. Desiante, N. Di Lalla, L. Di Venere, F. Fana Dirirsa, S. J. Fegan, Y. Fukazawa , et al. (60 additional authors not shown)

Abstract: We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive har… ▽ More We present the first Fermi - Large Area Telescope (LAT) solar flare catalog covering the 24 th solar cycle. This catalog contains 45 Fermi -LAT solar flares (FLSFs) with emission in the gamma-ray energy band (30 MeV - 10 GeV) detected with a significance greater than 5 sigma over the years 2010-2018. A subsample containing 37 of these flares exhibit delayed emission beyond the prompt-impulsive hard X-ray phase with 21 flares showing delayed emission lasting more than two hours. No prompt-impulsive emission is detected in four of these flares. We also present in this catalog the observations of GeV emission from 3 flares originating from Active Regions located behind the limb (BTL) of the visible solar disk. We report the light curves, spectra, best proton index and localization (when possible) for all the FLSFs. The gamma-ray spectra is consistent with the decay of pions produced by >300 MeV protons. This work contains the largest sample of high-energy gamma-ray flares ever reported and provides the unique opportunity to perform population studies on the different phases of the flare and thus allowing to open a new window in solar physics. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: Published in ApJS

Journal ref: ApJS 252 13 (2021)

Fermi and Swift Observations of GRB 190114C: Tracing the Evolution of High-Energy Emission from Prompt to Afterglow

Authors: M. Ajello, M. Arimoto, M. Axelsson, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, A. Berretta, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, E. Cavazzuti, S. Chen, G. Chiaro, S. Ciprini, J. Cohen-Tanugi , et al. (125 additional authors not shown)

Abstract: We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transiti… ▽ More We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process. △ Less

Submitted 23 January, 2020; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: 24 pages, 9 figures, Accepted to ApJ

Error estimates of the backward Euler-Maruyama method for multi-valued stochastic differential equations

Authors: Monika Eisenmann, Mihály Kovács, Raphael Kruse, Stig Larsson

Abstract: In this paper, we derive error estimates of the backward Euler-Maruyama method applied to multi-valued stochastic differential equations. An important example of such an equation is a stochastic gradient flow whose associated potential is not continuously differentiable, but assumed to be convex. We show that the backward Euler-Maruyama method is well-defined and convergent of order at least… ▽ More In this paper, we derive error estimates of the backward Euler-Maruyama method applied to multi-valued stochastic differential equations. An important example of such an equation is a stochastic gradient flow whose associated potential is not continuously differentiable, but assumed to be convex. We show that the backward Euler-Maruyama method is well-defined and convergent of order at least $1/4$ with respect to the root-mean-square norm. Our error analysis relies on techniques for deterministic problems developed in [Nochetto, Savaré, and Verdi, Comm.\ Pure Appl.\ Math., 2000]. We verify that our setting applies to an overdamped Langevin equation with a discontinuous gradient and to a spatially semi-discrete approximation of the stochastic $p$-Laplace equation. △ Less

Submitted 27 June, 2019; originally announced June 2019.

MSC Class: 65C30; 60H10; 34A60

Journal ref: BIT Numer. Math. (2021)

A greedy algorithm for optimal heating in powder-bed-based additive manufacturing

Authors: Robert Forslund, Anders Snis, Stig Larsson

Abstract: Powder-bed-based additive manufacturing involves melting of a powder bed using a moving laser or electron beam as a heat source. In this paper, we formulate an optimization scheme that aims to control this type of melting. The goal consists of tracking maximum temperatures on lines that run along the beam path. Time-dependent beam parameters (more specifically, beam power, spot size, and speed) ac… ▽ More Powder-bed-based additive manufacturing involves melting of a powder bed using a moving laser or electron beam as a heat source. In this paper, we formulate an optimization scheme that aims to control this type of melting. The goal consists of tracking maximum temperatures on lines that run along the beam path. Time-dependent beam parameters (more specifically, beam power, spot size, and speed) act as control functions. The scheme is greedy in the sense that it exploits local properties of the melt pool in order to divide a large optimization problem into several small ones. As illustrated by numerical examples, the scheme can resolve heat conduction issues such as concentrated heat accumulation at turning points and non-uniform melt depths. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Journal ref: J. Math. Ind. 11 (2021), Paper No. 14, 23 pp

MAGIC and Fermi-LAT gamma-ray results on unassociated HAWC sources

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak, P. Colin , et al. (318 additional authors not shown)

Abstract: The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the fir… ▽ More The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail. △ Less

Submitted 13 January, 2019; originally announced January 2019.

Comments: 12 pages, 3 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 485, Issue 1, May 2019, Pages 356-366

Unresolved Gamma-Ray Sky through its Angular Power Spectrum

Authors: M. Ackermann, M. Ajello, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, E. Cavazzuti, S. Chen, G. Chiaro, S. Ciprini, D. Costantin, A. Cuoco , et al. (85 additional authors not shown)

Abstract: The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluc… ▽ More The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$σ$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15. △ Less

Submitted 3 May, 2019; v1 submitted 5 December, 2018; originally announced December 2018.

Comments: 8 + 12 (SOM) pages, Contact author: Michela Negro (michela.negro@to.infn.it). (Typos in Tab. VII in previous versions)

Journal ref: Phys. Rev. Lett. 121, 241101 (2018)

VERITAS and Fermi-LAT observations of new HAWC sources

Authors: VERITAS Collaboration, A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, A. J. Chromey, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, M. Hutten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson , et al. (259 additional authors not shown)

Abstract: The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detect… ▽ More The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Comments: Accepted for publication in the ApJ, Corresponding author: Nahee Park (VERITAS Collaboration), John W. Hewitt (Fermi-LAT Collaboration), Ignacio Taboada (HAWC Collaboration), 30 pages, 12 figures

Journal ref: ApJ 866 (2018) no.1, 24

Analytical solution for heat conduction due to a moving Gaussian heat flux with piecewise constant parameters

Authors: Robert Forslund, Anders Snis, Stig Larsson

Abstract: We provide an analytical solution of the heat equation in the half-space subject to a moving Gaussian heat flux with piecewise constant parameters. The solution is of interest in powder bed fusion applications where these parameters can be used to control the conduction of heat due to a scanning beam of concentrated energy. The analytical solution is written in a dimensionless form as a sum of int… ▽ More We provide an analytical solution of the heat equation in the half-space subject to a moving Gaussian heat flux with piecewise constant parameters. The solution is of interest in powder bed fusion applications where these parameters can be used to control the conduction of heat due to a scanning beam of concentrated energy. The analytical solution is written in a dimensionless form as a sum of integrals over (dimensionless) time. For the numerical computation of these integrals we suggest a quadrature scheme that utilizes pre-calculated look-up tables for the required quadrature orders. Such a scheme is efficient because the required quadrature orders are strongly dependent on the parameters in the heat flux. The possibilities of using the obtained computational technique for the control and optimization of powder bed fusion processes are discussed. △ Less

Submitted 28 March, 2018; originally announced March 2018.

MSC Class: 65M80

Journal ref: Appl. Math. Model. 66 (2019), 227-240

Mittag-Leffler Euler integrator for a stochastic fractional order equation with additive noise

Authors: Mihály Kovács, Stig Larsson, Fardin Saedpanah

Abstract: Motivated by fractional derivative models in viscoelasticity, a class of semilinear stochastic Volterra integro-differential equations, and their deterministic counterparts, are considered. A generalized exponential Euler method, named here as the Mittag-Leffler Euler integrator, is used for the temporal discretization, while the spatial discretization is performed by the spectral Galerkin method.… ▽ More Motivated by fractional derivative models in viscoelasticity, a class of semilinear stochastic Volterra integro-differential equations, and their deterministic counterparts, are considered. A generalized exponential Euler method, named here as the Mittag-Leffler Euler integrator, is used for the temporal discretization, while the spatial discretization is performed by the spectral Galerkin method. The temporal rate of strong convergence is found to be (almost) twice compared to when the backward Euler method is used together with a convolution quadrature for time discretization. Numerical experiments that validate the theory are presented. △ Less

Submitted 16 January, 2020; v1 submitted 12 March, 2018; originally announced March 2018.

Comments: 20 pages, 5 figures

Journal ref: SIAM J. Numer. Anal. 58 (2020), 66-85

$J$-factors for self-interacting dark matter in 20 dwarf spheroidal galaxies

Authors: Sebastian Bergström, Riccardo Catena, Andrea Chiappo, Jan Conrad, Björn Eurenius, Magdalena Eriksson, Michael Högberg, Susanna Larsson, Emelie Olsson, Andreas Unger, Rikard Wadman

Abstract: Dwarf spheroidal galaxies are among the most promising targets for indirect dark matter (DM) searches in $γ$-rays. The $γ$-ray flux from DM annihilation in a dwarf spheroidal galaxy is proportional to the $J$-factor of the source. The $J$-factor of a dwarf spheroidal galaxy is the line-of-sight integral of the DM mass density squared times… ▽ More Dwarf spheroidal galaxies are among the most promising targets for indirect dark matter (DM) searches in $γ$-rays. The $γ$-ray flux from DM annihilation in a dwarf spheroidal galaxy is proportional to the $J$-factor of the source. The $J$-factor of a dwarf spheroidal galaxy is the line-of-sight integral of the DM mass density squared times $\langle σ_{\rm ann} v_{\rm rel} \rangle/(σ_{\rm ann} v_{\rm rel})_0$, where $σ_{\rm ann} v_{\rm rel}$ is the DM annihilation cross-section times relative velocity $v_{\rm rel}=|{\bf v}_{\rm rel}|$, angle brackets denote average over ${\bf v}_{\rm rel}$, and $(σ_{\rm ann} v_{\rm rel})_0$ is the $v_{\rm rel}$-independent part of $σ_{\rm ann} v_{\rm rel}$. If $σ_{\rm ann} v_{\rm rel}$ is constant in $v_{\rm rel}$, $J$-factors only depend on the DM space distribution in the source. However, if $σ_{\rm ann} v_{\rm rel}$ varies with $v_{\rm rel}$, as in the presence of DM self-interactions, $J$-factors also depend on the DM velocity distribution, and on the strength and range of the DM self-interaction. Models for self-interacting DM are increasingly important in the study of the small scale clustering of DM, and are compatible with current cosmological observations. Here we derive the $J$-factor of 20 dwarf spheroidal galaxies from stellar kinematic data under the assumption of Yukawa DM self-interactions. $J$-factors are derived through a profile Likelihood approach, assuming either NFW or cored DM profiles. We also compare our results with $J$-factors derived assuming the same velocity for all DM particles in the target galaxy. We find that this common approximation overestimates the $J$-factors by up to one order of magnitude. $J$-factors for a sample of DM particle masses, self-interaction coupling constants and density profiles are provided electronically, ready to be used in other projects. △ Less

Submitted 8 December, 2017; originally announced December 2017.

Comments: 10 pages, 3 figures and 2 tables

Journal ref: Phys. Rev. D 98, 043017 (2018)

On a randomized backward Euler method for nonlinear evolution equations with time-irregular coefficients

Authors: Monika Eisenmann, Mihály Kovács, Raphael Kruse, Stig Larsson

Abstract: In this paper we introduce a randomized version of the backward Euler method, that is applicable to stiff ordinary differential equations and nonlinear evolution equations with time-irregular coefficients. In the finite-dimensional case, we consider Carathéodory type functions satisfying a one-sided Lipschitz condition. After investigating the well-posedness and the stability properties of the ran… ▽ More In this paper we introduce a randomized version of the backward Euler method, that is applicable to stiff ordinary differential equations and nonlinear evolution equations with time-irregular coefficients. In the finite-dimensional case, we consider Carathéodory type functions satisfying a one-sided Lipschitz condition. After investigating the well-posedness and the stability properties of the randomized scheme, we prove the convergence to the exact solution with a rate of $0.5$ in the root-mean-square norm assuming only that the coefficient function is square integrable with respect to the temporal parameter. These results are then extended to the numerical solution of infinite-dimensional evolution equations under monotonicity and Lipschitz conditions. Here we consider a combination of the randomized backward Euler scheme with a Galerkin finite element method. We obtain error estimates that correspond to the regularity of the exact solution. The practicability of the randomized scheme is also illustrated through several numerical experiments. △ Less

Submitted 11 December, 2018; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: 37 pages, 3 figures

MSC Class: 65C05; 65L05; 65L20; 65M12; 65M60

Journal ref: Found. Comput. Math. 19 (2019), 1387-1430

MAGIC observations of variable very-high-energy gamma-ray emission from PKS1510-089 during May 2015 outburst

Authors: Julian Sitarek, Josefa Becerra González, Vandad Fallah Ramazani Elina Lindfors, Giovanna Pedaletti, Fabrizio Tavecchio, Monica Vazquez Acosta, Stefan Larsson, the MAGIC Collaboration, the Fermi-LAT Collaboration, Kiran Baliyan, Navpreet Kaur, Sameer, Svetlana Jorstad, Claudia Raiteri

Abstract: PKS1510-089 is a flat spectrum radio quasar located at a redshift of 0.36. It is one of only a few such sources detected in very-high-energy (VHE, >100 GeV) gamma rays. Though PKS1510-089 is highly variable at GeV energies, until recently no variability has been observed in the VHE band. In 2015 May PKS1510-089 showed a high state in optical and in the GeV range. A VHE gamma-ray flare was detected… ▽ More PKS1510-089 is a flat spectrum radio quasar located at a redshift of 0.36. It is one of only a few such sources detected in very-high-energy (VHE, >100 GeV) gamma rays. Though PKS1510-089 is highly variable at GeV energies, until recently no variability has been observed in the VHE band. In 2015 May PKS1510-089 showed a high state in optical and in the GeV range. A VHE gamma-ray flare was detected with MAGIC at that time, showing the first instance of VHE gamma-ray flux variability on the time scale of days in this source. We will present the MAGIC results from this observation, discuss their temporal and spectral properties in the multi-wavelength context and present modelling of such emission in the external Compton scenario. △ Less

Submitted 21 August, 2017; v1 submitted 14 August, 2017; originally announced August 2017.

Comments: Proceedings of the 35th International Cosmic Ray Conference (ICRC 2017), Bexco, Busan, Korea (arXiv:1708.05153)

Report number: MAGIC-ICRC/2017/04

Femtosecond X-ray Fourier holography imaging of free-flying nanoparticles

Authors: Tais Gorkhover, Anatoli Ulmer, Ken Ferguson, Max Bucher, Filipe Maia, Johan Bielecki, Tomas Ekeberg, Max F. Hantke, Benedikt J. Daurer, Carl Nettelblad, Jakob Andreasson, Anton Barty, Petr Bruza, Sebastian Carron, Dirk Hasse, Jacek Krzywinski, Daniel S. D. Larsson, Andrew Morgan, Kerstin Muehlig, Maria Mueller, Kenta Okamoto, Alberto Pietrini, Daniela Rupp, Mario Sauppe, Gijs van der Schot , et al. (13 additional authors not shown)

Abstract: Ultrafast X-ray imaging provides high resolution information on individual fragile specimens such as aerosols, metastable particles, superfluid quantum systems and live biospecimen, which is inaccessible with conventional imaging techniques. Coherent X-ray diffractive imaging, however, suffers from intrinsic loss of phase, and therefore structure recovery is often complicated and not always unique… ▽ More Ultrafast X-ray imaging provides high resolution information on individual fragile specimens such as aerosols, metastable particles, superfluid quantum systems and live biospecimen, which is inaccessible with conventional imaging techniques. Coherent X-ray diffractive imaging, however, suffers from intrinsic loss of phase, and therefore structure recovery is often complicated and not always uniquely-defined. Here, we introduce the method of in-flight holography, where we use nanoclusters as reference X-ray scatterers in order to encode relative phase information into diffraction patterns of a virus. The resulting hologram contains an unambiguous three-dimensional map of a virus and two nanoclusters with the highest lat- eral resolution so far achieved via single shot X-ray holography. Our approach unlocks the benefits of holography for ultrafast X-ray imaging of nanoscale, non-periodic systems and paves the way to direct observation of complex electron dynamics down to the attosecond time scale. △ Less

Submitted 9 August, 2017; v1 submitted 28 July, 2017; originally announced July 2017.

Fermi-LAT Observations of High-energy Behind-the-limb Solar Flares

Authors: M. Ackermann, A. Allafort, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, S. Ciprini, F. Costanza, S. Cutini, F. D'Ammando, F. de Palma, R. Desiante, S. W. Digel , et al. (64 additional authors not shown)

Abstract: We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (… ▽ More We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (CMEs) and strong solar energetic particle events. We present updated localizations of the >100 MeV photon emission, hard X-ray (HXR)and EUV images, and broadband spectra from 10 keV to 10 GeV, as well as microwave spectra. We also provide a comparison of the BTL flares detected by Fermi-LAT with three on-disk flares and present a study of some of the significant quantities of these flares as an attempt to better understand the acceleration mechanisms at work during these occulted flares. We interpret the HXR emission to be due to electron bremsstrahlung from a coronal thin-target loop top with the accelerated electron spectra steepening at semirelativistic energies. The >100 MeV gamma-rays are best described by a pion-decay model resulting from the interaction of protons (and other ions) in a thick-target photospheric source. The protons are believed to have been accelerated (to energies >10 GeV) in the CME environment and precipitate down to the photosphere from the downstream side of the CME shock and landed on the front side of the Sun, away from the original flare site and the HXR emission. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 13 pages, 14 figures published on ApJ

Journal ref: ApJ 835, 219, 1 February 2017

Search for extended sources in the Galactic Plane using 6 years of Fermi-Large Area Telescope Pass 8 data above 10 GeV

Authors: The Fermi LAT Collaboration, M. Ackermann, M. Ajello, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, R. A. Cameron, M. Caragiulo, P. A. Caraveo, D. Castro, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, C. C. Cheung , et al. (95 additional authors not shown)

Abstract: The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two… ▽ More The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain. △ Less

Submitted 11 April, 2018; v1 submitted 1 February, 2017; originally announced February 2017.

Comments: 33 pages, 22 figures & 3 tables. Published by The Astrophysical Journal. Available on the Fermi Science Support Center (FSSC) together with the 3FHL catalog

Strong convergence of a fully discrete finite element approximation of the stochastic Cahn-Hilliard equation

Authors: Daisuke Furihata, Mihály Kovács, Stig Larsson, Fredrik Lindgren

Abstract: We consider the stochastic Cahn-Hilliard equation driven by additive Gaussian noise in a convex domain with polygonal boundary in dimension $d\le 3$. We discretize the equation using a standard finite element method in space and a fully implicit backward Euler method in time. By proving optimal error estimates on subsets of the probability space with arbitrarily large probability and uniform-in-ti… ▽ More We consider the stochastic Cahn-Hilliard equation driven by additive Gaussian noise in a convex domain with polygonal boundary in dimension $d\le 3$. We discretize the equation using a standard finite element method in space and a fully implicit backward Euler method in time. By proving optimal error estimates on subsets of the probability space with arbitrarily large probability and uniform-in-time moment bounds we show that the numerical solution converges strongly to the solution as the discretization parameters tend to zero. △ Less

Submitted 2 January, 2018; v1 submitted 30 December, 2016; originally announced December 2016.

Comments: 25 pages

MSC Class: 60H15; 60H35; 65C30

Journal ref: SIAM J. Numer. Anal. 56(2) (2018) 708-731

The second catalog of flaring gamma-ray sources from the Fermi All-sky Variability Analysis

Authors: S. Abdollahi, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, S. Buson, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman , et al. (102 additional authors not shown)

Abstract: We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data relea… ▽ More We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6σ$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2. △ Less

Submitted 12 September, 2017; v1 submitted 9 December, 2016; originally announced December 2016.

Comments: Accepted for publication by The Astrophysical Journal Supplement. Online analysis results available at http://fermi.gsfc.nasa.gov/ssc/data/access/lat/FAVA/index.php . Corresponding authors: M. Giomi, R. Buehler, D. Kocevski, and M. Ajello

Gamma-ray blazar spectra with H.E.S.S. II mono analysis: the case of PKS 2155-304 and PG 1553+113

Authors: H. E. S. S. Collaboration, :, H. Abdalla, A. Abramowski, F. Aharonian, F. Ait Benkhali, A. G. Akhperjanian, T. Andersson, E. O. Angüner, M. Arrieta, P. Aubert, M. Backes, A. Balzer, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, P. Bordas, F. Brun , et al. (311 additional authors not shown)

Abstract: The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft sp… ▽ More The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 600, A89 (2017)

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, R. Carosi, A. Carosi, A. Chatterjee, P. Colin , et al. (151 additional authors not shown)

Abstract: Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In… ▽ More Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models. △ Less

Submitted 28 February, 2017; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: 12 pages, 10 figures, accepted for a publication in A&A

Journal ref: A&A 603, A29 (2017)

The F-GAMMA program: Multi-frequency study of Active Galactic Nuclei in the Fermi era. Program description and the first 2.5 years of monitoring

Authors: L. Fuhrmann, E. Angelakis, J. A. Zensus, I. Nestoras, N. Marchili, V. Pavlidou, V. Karamanavis, H. Ungerechts, T. P. Krichbaum, S. Larsson, S. S. Lee, W. Max-Moerbeck, I. Myserlis, T. J. Pearson, A. C. S. Readhead, J. L. Richards, A. Sievers, B. W. Sohn

Abstract: To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally… ▽ More To fully exploit the scientific potential of the Fermi mission, we initiated the F-GAMMA program. Between 2007 and 2015 it was the prime provider of complementary multi-frequency monitoring in the radio regime. We quantify the radio variability of gamma-ray blazars. We investigate its dependence on source class and examine whether the radio variability is related to the gamma-ray loudness. Finally, we assess the validity of a putative correlation between the two bands. The F-GAMMA monitored monthly a sample of about 60 sources at up to twelve radio frequencies between 2.64 and 228.39 GHz. We perform a time series analysis on the first 2.5-year dataset to obtain variability parameters. A maximum likelihood analysis is used to assess the significance of a correlation between radio and gamma-ray fluxes. We present light curves and spectra (coherent within ten days) obtained with the Effelsberg 100-m and IRAM 30-m telescopes. All sources are variable across all frequency bands with amplitudes increasing with frequency up to rest frame frequencies of around 60 - 80 GHz as expected by shock-in-jet models. Compared to FSRQs, BL Lacs show systematically lower variability amplitudes, brightness temperatures and Doppler factors at lower frequencies, while the difference vanishes towards higher ones. The time scales appear similar for the two classes. The distribution of spectral indices appears flatter or more inverted at higher frequencies for BL Lacs. Evolving synchrotron self-absorbed components can naturally account for the observed spectral variability. We find that the Fermi-detected sources show larger variability amplitudes as well as brightness temperatures and Doppler factors, than non-detected ones. Flux densities at 86.2 and 142.3 GHz correlate with 1 GeV fluxes at a significance level better than 3sigma, implying that gamma rays are produced very close to the mm-band emission region. △ Less

Submitted 7 August, 2016; originally announced August 2016.

Comments: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics (18 pages, 9 figures)

Journal ref: A&A 596, A45 (2016)

Searching the Gamma-ray Sky for Counterparts to Gravitational Wave Sources: Fermi GBM and LAT Observations of LVT151012 and GW151226

Authors: J. L. Racusin, E. Burns, A. Goldstein, V. Connaughton, C. A. Wilson-Hodge, P. Jenke, L. Blackburn, M. S. Briggs, J. Broida, J. Camp, N. Christensen, C. M. Hui, T. Littenberg, P. Shawhan, L. Singer, J. Veitch, P. N. Bhat, W. Cleveland, G. Fitzpatrick, M. H. Gibby, A. von Kienlin, S. McBreen, B. Mailyan, C. A. Meegan, W. S. Paciesas , et al. (116 additional authors not shown)

Abstract: We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techn… ▽ More We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: 27 pages, 1 table, 11 figures, Submitted to ApJ

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the dif… ▽ More This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands. △ Less

Submitted 21 July, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: For the main Letter, see arXiv:1602.08492

Report number: LIGO-P1600137-v2

Journal ref: The Astrophysical Journal Supplement Series, 225:8 (15pp), 2016 July

Quasi-optimality of Petrov-Galerkin discretizations of parabolic problems with random coefficients

Authors: Stig Larsson, Christian Mollet, Matteo Molteni

Abstract: We consider a linear parabolic problem with random elliptic operator in the usual Gelfand triple setting. We do not assume uniform bounds on the coercivity and boundedness constants, but allow them to be random variables. The parabolic problem is studied in a weak space-time formulation, where we can derive explicit formulas for the inf-sup constants. Under suitable assumptions we prove existence… ▽ More We consider a linear parabolic problem with random elliptic operator in the usual Gelfand triple setting. We do not assume uniform bounds on the coercivity and boundedness constants, but allow them to be random variables. The parabolic problem is studied in a weak space-time formulation, where we can derive explicit formulas for the inf-sup constants. Under suitable assumptions we prove existence of moments of the solution. We also prove quasi-optimal error estimates for piecewise polynomial Petrov-Galerkin discretizations. △ Less

Submitted 22 April, 2016; originally announced April 2016.

Comments: 20 pages, 6 figures

Discrete Variational Derivative Methods for the EPDiff equation

Authors: Stig Larsson, Takayasu Matsuo, Klas Modin, Matteo Molteni

Abstract: The aim of this paper is the derivation of structure preserving schemes for the solution of the EPDiff equation, with particular emphasis on the two dimensional case. We develop three different schemes based on the Discrete Variational Derivative Method (DVDM) on a rectangular domain discretized with a regular, structured, orthogonal grid. We present numerical experiments to support our claims:… ▽ More The aim of this paper is the derivation of structure preserving schemes for the solution of the EPDiff equation, with particular emphasis on the two dimensional case. We develop three different schemes based on the Discrete Variational Derivative Method (DVDM) on a rectangular domain discretized with a regular, structured, orthogonal grid. We present numerical experiments to support our claims: we investigate the preservation of energy and linear momenta, the reversibility, and the empirical convergence of the schemes. The quality of our schemes is finally tested by simulating the interaction of singular wave fronts. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Comments: 41 pages, 41 figures

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

Authors: M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprini , et al. (90 additional authors not shown)

Abstract: We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interac… ▽ More We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data. △ Less

Submitted 13 April, 2016; v1 submitted 12 April, 2016; originally announced April 2016.

Comments: Published by PRD, 16 pages, 11 figures, corresponding authors: F. Loparco and M. N. Mazziotta

Numerical solution of parabolic problems based on a weak space-time formulation

Authors: Stig Larsson, Matteo Molteni

Abstract: We investigate a weak space-time formulation of the heat equation and its use for the construction of a numerical scheme. The formulation is based on a known weak space-time formulation, with the difference that a pointwise component of the solution, which in other works is usually neglected, is now kept. We investigate the role of such a component by first using it to obtain a pointwise bound on… ▽ More We investigate a weak space-time formulation of the heat equation and its use for the construction of a numerical scheme. The formulation is based on a known weak space-time formulation, with the difference that a pointwise component of the solution, which in other works is usually neglected, is now kept. We investigate the role of such a component by first using it to obtain a pointwise bound on the solution and then deploying it to construct a numerical scheme. The scheme obtained, besides being quasi-optimal in the $L^2$ sense, is also pointwise superconvergent in the temporal nodes. We prove a priori error estimates and we present numerical experiments to empirically support our findings. △ Less

Submitted 10 March, 2016; originally announced March 2016.

Comments: 23 pages, 8 figures

MSC Class: 65M15; 65M60

Localization and broadband follow-up of the gravitational-wave transient GW150914

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, M. R. Abernathy, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, A. Allocca, P. A. Altin, S. B. Anderson, W. G. Anderson, K. Arai , et al. (1522 additional authors not shown)

Abstract: A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared wit… ▽ More A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams. △ Less

Submitted 21 July, 2016; v1 submitted 26 February, 2016; originally announced February 2016.

Comments: For Supplement, see https://arxiv.org/abs/1604.07864

Report number: LIGO-P1500227-v12

Journal ref: The Astrophysical Journal Letters, 826:L13 (8pp), 2016 July 20

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Authors: F. Acero, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, C. Cecchi , et al. (109 additional authors not shown)

Abstract: Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Ga… ▽ More Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center. △ Less

Submitted 23 February, 2016; originally announced February 2016.

Comments: Accepted for publication in ApJS

Journal ref: ApJS 223 26 (2016)

Radio follow-up of the gamma-ray flaring gravitational lens JVAS B0218+357

Authors: Cristiana Spingola, D. Dallacasa, M. Orienti, M. Giroletti, J. P. McKean, C. C. Cheung, T. Hovatta, S. Ciprini, F. D'Ammando, E. Falco, S. Larsson, W. Max-Moerbeck, R. Ojha, A. C. S. Readhead, J. L. Richards, J. Scargle

Abstract: We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the gamma-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density va… ▽ More We present results on multifrequency Very Long Baseline Array (VLBA) monitoring observations of the double-image gravitationally lensed blazar JVAS B0218+357. Multi-epoch observations started less than one month after the gamma-ray flare detected in 2012 by the Large Area Telescope on board Fermi, and spanned a 2-month interval. The radio light curves did not reveal any significant flux density variability, suggesting that no clear correlation between the high energy and low-energy emission is present. This behaviour was confirmed also by the long-term Owens Valley Radio Observatory monitoring data at 15 GHz. The milliarcsecond-scale resolution provided by the VLBA observations allowed us to resolve the two images of the lensed blazar, which have a core-jet structure. No significant morphological variation is found by the analysis of the multi-epoch data, suggesting that the region responsible for the gamma-ray variability is located in the core of the AGN, which is opaque up to the highest observing frequency of 22 GHz. △ Less

Submitted 14 January, 2016; originally announced January 2016.

Comments: 9 pages, 6 figures, 3 tables. Accepted for publication in MNRAS

Very-high-energy gamma-rays from the Universe's middle age: detection of the z=0.940 blazar PKS 1441+25 with MAGIC

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, W. Bednarek, E. Bernardini, B. Biassuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, A. Chatterjee, R. Clavero, P. Colin, E. Colombo , et al. (229 additional authors not shown)

Abstract: The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies… ▽ More The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 σ using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy. △ Less

Submitted 12 January, 2018; v1 submitted 14 December, 2015; originally announced December 2015.

Comments: Corresponding Authors: J. Becerra (josefa.becerra@nasa.gov), M. Nievas Rosillo (mnievas@ucm.es), M. Manganaro (manganaro@iac.es), F. Tavecchio (fabrizio.tavecchio@brera.inaf.it) Published in The Astrophysical Journal Letters (ApJL), 2015ApJ...815L..23A, DOI: 10.1088/2041-8205/815/2/L23

Journal ref: ApJL, 2015, 815, L23

The 1st Fermi Lat Supernova Remnant Catalog

Authors: Fabio Acero, Markus Ackermann, Marco Ajello, Luca Baldini, Jean Ballet, Guido Barbiellini, Denis Bastieri, Ronaldo Bellazzini, E. Bissaldi, Roger Blandford, E. D. Bloom, Raffaella Bonino, Eugenio Bottacini, J. Bregeon, Philippe Bruel, Rolf Buehler, S. Buson, G. A. Caliandro, Rob A. Cameron, R Caputo, Micaela Caragiulo, Patrizia A. Caraveo, Jean Marc Casandjian, Elisabetta Cavazzuti, Claudia Cecchi , et al. (134 additional authors not shown)

Abstract: To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245… ▽ More To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: Resubmitted to ApJS

Journal ref: ApJS 224 8 (2016)

On the discretisation in time of the stochastic Allen-Cahn equation

Authors: Mihály Kovács, Stig Larsson, Fredrik Lindgren

Abstract: We consider the stochastic Allen--Cahn equation perturbed by smooth additive Gaussian noise in a spatial domain with smooth boundary in dimension $d\le 3$, and study the semidiscretisation in time of the equation by an Euler type split-step method. We show that the method converges strongly with a rate $O(Δt^{\frac12}) $. By means of a perturbation argument, we also establish the strong convergenc… ▽ More We consider the stochastic Allen--Cahn equation perturbed by smooth additive Gaussian noise in a spatial domain with smooth boundary in dimension $d\le 3$, and study the semidiscretisation in time of the equation by an Euler type split-step method. We show that the method converges strongly with a rate $O(Δt^{\frac12}) $. By means of a perturbation argument, we also establish the strong convergence of the standard backward Euler scheme with the same rate. △ Less

Submitted 2 January, 2018; v1 submitted 13 October, 2015; originally announced October 2015.

Comments: 32 pages. Several typos fixed, the statement and proof of Lemma 3.2 are updated and so is the proof of Lemma 5.2. The proof of Lemma 5.1 is omitted

MSC Class: 60H15; 60H35; 65C30

Journal ref: Mathematische Nachrichten, Volume 291, Issue 5-6, April 2018, Pages 966-995

Search for extended gamma-ray emission from the Virgo galaxy cluster with Fermi-LAT

Authors: M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti , et al. (96 additional authors not shown)

Abstract: Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for search… ▽ More Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $γ$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3°that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $γ$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle σv \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $γ$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$. △ Less

Submitted 30 September, 2015; originally announced October 2015.

Comments: 15 pages, 11 figures, 4 tables, accepted for publication in ApJ; corresponding authors: T. Jogler, S. Zimmer & A. Pinzke

Journal ref: ApJ 812 159 (2015)

Multiwavelength Evidence for Quasi-periodic Modulation in the Gamma-ray Blazar PG 1553+113

Authors: The Fermi LAT collaboration, M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo , et al. (117 additional authors not shown)

Abstract: We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated osci… ▽ More We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity. △ Less

Submitted 12 October, 2015; v1 submitted 7 September, 2015; originally announced September 2015.

Comments: 8 pages, 5 figures. Accepted to The Astrophysical Journal Letters. Corresponding authors: S. Ciprini (ASDC/INFN), S. Cutini (ASDC/INFN), S. Larsson (Stockholm Univ/KTH), A. Stamerra (INAF/SNS), D. J. Thompson (NASA GSFC)

PSR J1906+0722: An Elusive Gamma-ray Pulsar

Authors: C. J. Clark, H. J. Pletsch, J. Wu, L. Guillemot, M. Ackermann, B. Allen, A. de Angelis, C. Aulbert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, O. Bock, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo , et al. (95 additional authors not shown)

Abstract: We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT s… ▽ More We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: 7 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

Journal ref: ApJ 809 L2 (2015)

Search for Early Gamma-ray Production in Supernovae Located in a Dense Circumstellar Medium with the Fermi LAT

Authors: M. Ackermann, I. Arcavi, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang , et al. (86 additional authors not shown)

Abstract: Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We se… ▽ More Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf). △ Less

Submitted 26 June, 2015; v1 submitted 4 June, 2015; originally announced June 2015.

Comments: Accepted for publication in ApJ. Corresponding author: A. Franckowiak (afrancko@slac.stanford.edu), updated author list and acknowledgements

Journal ref: ApJ, 807, 169 (2015)

Covariance structure of parabolic stochastic partial differential equations with multiplicative Lévy noise

Authors: Kristin Kirchner, Annika Lang, Stig Larsson

Abstract: The characterization of the covariance function of the solution process to a stochastic partial differential equation is considered in the parabolic case with multiplicative Lévy noise of affine type. For the second moment of the mild solution, a well-posed deterministic space-time variational problem posed on projective and injective tensor product spaces is derived, which subsequently leads to a… ▽ More The characterization of the covariance function of the solution process to a stochastic partial differential equation is considered in the parabolic case with multiplicative Lévy noise of affine type. For the second moment of the mild solution, a well-posed deterministic space-time variational problem posed on projective and injective tensor product spaces is derived, which subsequently leads to a deterministic equation for the covariance function. △ Less

Submitted 16 April, 2016; v1 submitted 1 June, 2015; originally announced June 2015.

Comments: 28 pages

MSC Class: 60H15; 35K90; 35R60

Journal ref: J. Differential Equations 262 (2017), 5896-5927

Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise

Authors: Rikard Anton, David Cohen, Stig Larsson, Xiaojie Wang

Abstract: A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size… ▽ More A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size restriction as in the often used Störmer-Verlet-leap-frog scheme. Furthermore, it satisfies an almost trace formula (i.e., a linear drift of the expected value of the energy of the problem). Numerical experiments are presented and confirm the theoretical results. △ Less

Submitted 25 November, 2015; v1 submitted 28 February, 2015; originally announced March 2015.

MSC Class: 65C20; 60H10; 60H15; 60H35; 65C30