Showing 1–46 of 46 results for author: Vavrek, R

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Differential fuzz testing to detect tampering in sensor systems and its application to arms control authentication

Authors: Jayson R Vavrek, Luozhong Zhou, Joshua Boverhof, Elisa R Heymann, Barton P Miller, Sean Peisert

Abstract: In future nuclear arms control treaties, it will be necessary to authenticate the hardware and software components of verification measurement systems, i.e., to ensure these systems are functioning as intended and have not been tampered with by malicious actors. While methods such as source code hashing and static analysis can help verify the integrity of software components, they may not be capab… ▽ More In future nuclear arms control treaties, it will be necessary to authenticate the hardware and software components of verification measurement systems, i.e., to ensure these systems are functioning as intended and have not been tampered with by malicious actors. While methods such as source code hashing and static analysis can help verify the integrity of software components, they may not be capable of detecting tampering with environment variables, external libraries, or the firmware and hardware of radiation measurement systems. In this article, we introduce the concept of physical differential fuzz testing as a challenge-response-style tamper indicator that can holistically and simultaneously test all the above components in a cyber-physical system. In essence, we randomly sample (or "fuzz") the untampered system's parameter space, including both normal and off-normal parameter values, and consider the time series of outputs as the baseline signature of the system. Re-running the same input sequence on a untampered system will produce an output sequence consistent with this baseline, while running the same input sequence on a tampered system will produce a modified output sequence and raise an alarm. We then apply this concept to authenticating the radiation measurement equipment in nuclear weapon verification systems and conduct demonstration fuzz testing measurements with a sodium iodide (NaI) gamma ray spectrometer. Because there is Poisson noise in the measured output spectra, we also use a mechanism for comparing inherently noisy or stochastic fuzzing sequences. We show that physical differential fuzz testing can detect two types of tamper attempts, and conclude that it is a promising framework for authenticating future cyber-physical systems in nuclear arms control, safeguards, and beyond. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 15 pages, 7 figures, 1 table

Mapping the Minimum Detectable Activities of Gamma-Ray Sources in a 3-D Scene

Authors: Mark S. Bandstra, Daniel Hellfeld, Jaewon Lee, Brian J. Quiter, Marco Salathe, Jayson R. Vavrek, Tenzing H. Y. Joshi

Abstract: The ability to formulate maps of minimum detectable activities (MDAs) that describe the sensitivity of an ad hoc measurement that used one or more freely moving radiation detector systems would be significantly beneficial for the conduct and understanding of many radiological search activities. In a real-time scenario with a free-moving detector system, an MDA map can provide useful feedback to th… ▽ More The ability to formulate maps of minimum detectable activities (MDAs) that describe the sensitivity of an ad hoc measurement that used one or more freely moving radiation detector systems would be significantly beneficial for the conduct and understanding of many radiological search activities. In a real-time scenario with a free-moving detector system, an MDA map can provide useful feedback to the operator about which areas have not been searched as thoroughly as others, thereby allowing the operator to prioritize future actions. Similarly, such a calculation could be used to inform subsequent navigation decisions of autonomous platforms. Here we describe a near real-time MDA mapping approach that can be applied when searching for point sources using detected events in a spectral region of interest while assuming a constant, unknown background rate. We show the application of this MDA mapping method to a real scenario, a survey of the interior of a small building using a handheld detector system. Repeated measurements with no sources and with Cs-137 sources of different strengths yield results consistent with the estimated thresholds and MDA values; namely, that for background-only measurements no sources are seen above threshold anywhere in the scene, while when sources are present they are detected above the thresholds calculated for their locations. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 12 pages, 9 figures, submitted to IEEE Transactions on Nuclear Science

Design and deployment of radiological point-source arrays for the emulation of continuous distributed sources

Authors: Jayson R. Vavrek, C. Corey Hines, Mark S. Bandstra, Daniel Hellfeld, Maddison A. Heine, Zachariah M. Heiden, Nick R. Mann, Brian J. Quiter, Tenzing H. Y. Joshi

Abstract: We demonstrate a method for using arrays of point sources that emulate -- when measured from a standoff of at least several meters -- distributed gamma-ray sources, and present results using this method from outdoor aerial measurements of several planar arrays each comprising up to $100$ ${\sim}7$ mCi Cu-64 sealed sources. The method relies on the Poisson deviance to statistically test whether the… ▽ More We demonstrate a method for using arrays of point sources that emulate -- when measured from a standoff of at least several meters -- distributed gamma-ray sources, and present results using this method from outdoor aerial measurements of several planar arrays each comprising up to $100$ ${\sim}7$ mCi Cu-64 sealed sources. The method relies on the Poisson deviance to statistically test whether the array source ''looks like'' its continuous analogue to a particular gamma-ray detector given the counts recorded as the detector moves about 3D space. We use this deviance metric to design eight different mock distributed sources, ranging in complexity from a $36\times36$ m uniform square grid of sources to a configuration where regions of higher and zero activity are superimposed on a uniform baseline. We then detail the design, manufacture, and testing of the ${\sim} 7$ mCi Cu-64 sealed sources at the Washington State University research reactor, and their deployment during the aerial measurement campaign. We show the results of two such measurements, in which approximate source shapes and qualitative source intensities can be seen. Operationally, we find that the point-source array technique provides high source placement accuracy and ease of quantifying the true source configuration, scalability to source dimensions of ${\lesssim}100$ m, ease of reconfiguration and removal, and relatively low dose to personnel. Finally, we consider potential improvements and generalizations of the point-source array technique for future measurement campaigns. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 15 pages, 12 figures; submitted to IEEE Transactions on Nuclear Science

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

Authors: A. S. Borlaff, P. Gómez-Alvarez, B. Altieri, P. M. Marcum, R. Vavrek, R. Laureijs, R. Kohley, F. Buitrago, J. C. Cuillandre, P. A. Duc, L. M. Gaspar Venancio, A. Amara, S. Andreon, N. Auricchio, R. Azzollini, C. Baccigalupi, A. Balaguera-Antolínez, M. Baldi, S. Bardelli, R. Bender, A. Biviano, C. Bodendorf, D. Bonino, E. Bozzo, E. Branchini , et al. (158 additional authors not shown)

Abstract: While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightn… ▽ More While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 23 pages, 13 figures, Euclid Consortium Key Project, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 657, A92 (2022)

Euclid preparation: I. The Euclid Wide Survey

Authors: R. Scaramella, J. Amiaux, Y. Mellier, C. Burigana, C. S. Carvalho, J. -C. Cuillandre, A. Da Silva, A. Derosa, J. Dinis, E. Maiorano, M. Maris, I. Tereno, R. Laureijs, T. Boenke, G. Buenadicha, X. Dupac, L. M. Gaspar Venancio, P. Gómez-Álvarez, J. Hoar, J. Lorenzo Alvarez, G. D. Racca, G. Saavedra-Criado, J. Schwartz, R. Vavrek, M. Schirmer , et al. (216 additional authors not shown)

Abstract: Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced s… ▽ More Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 43 pages, 51 figures, submitted to A&A

Journal ref: A&A 662, A112 (2022)

Free-moving Quantitative Gamma-ray Imaging

Authors: Daniel Hellfeld, Mark S. Bandstra, Jayson R. Vavrek, Donald L. Gunter, Joseph C. Curtis, Marco Salathe, Ryan Pavlovsky, Victor Negut, Paul J. Barton, Joshua W. Cates, Brian J. Quiter, Reynold J. Cooper, Kai Vetter, Tenzing H. Y. Joshi

Abstract: The ability to map and estimate the activity of radiological source distributions in unknown three-dimensional environments has applications in the prevention and response to radiological accidents or threats as well as the enforcement and verification of international nuclear non-proliferation agreements. Such a capability requires well-characterized detector response functions, accurate time-dep… ▽ More The ability to map and estimate the activity of radiological source distributions in unknown three-dimensional environments has applications in the prevention and response to radiological accidents or threats as well as the enforcement and verification of international nuclear non-proliferation agreements. Such a capability requires well-characterized detector response functions, accurate time-dependent detector position and orientation data, a digitized representation of the surrounding 3D environment, and appropriate image reconstruction and uncertainty quantification methods. We have previously demonstrated 3D mapping of gamma-ray emitters with free-moving detector systems on a relative intensity scale using a technique called Scene Data Fusion (SDF). Here we characterize the detector response of a multi-element gamma-ray imaging system using experimentally benchmarked Monte Carlo simulations and perform 3D mapping on an absolute intensity scale. We present experimental reconstruction results from hand-carried and airborne measurements with point-like and distributed sources in known configurations, demonstrating quantitative SDF in complex 3D environments. △ Less

Submitted 19 October, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: 20 pages, 5 figures, 4 supplementary figures, published in Scientific Reports

Journal ref: Scientific Reports, 11, 20515 (2021)

Improved Gamma-Ray Point Source Quantification in Three Dimensions by Modeling Attenuation in the Scene

Authors: M. S. Bandstra, D. Hellfeld, J. R. Vavrek, B. J. Quiter, K. Meehan, P. J. Barton, J. W. Cates, A. Moran, V. Negut, R. Pavlovsky, T. H. Y. Joshi

Abstract: Using a series of detector measurements taken at different locations to localize a source of radiation is a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2D and even 3D,… ▽ More Using a series of detector measurements taken at different locations to localize a source of radiation is a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2D and even 3D, but few have studied the effect of intervening material on the problem. In this work we examine gamma-ray data taken from a freely moving system and develop voxelized 3-D models of the scene using data from the onboard LiDAR. Ray casting is used to compute the distance each gamma ray travels through the scene material, which is then used to calculate attenuation assuming a single attenuation coefficient for solids within the geometry. Parameter estimation using maximum likelihood is performed to simultaneously find the attenuation coefficient, source activity, and source position that best match the data. Using a simulation, we validate the ability of this method to reconstruct the true location and activity of a source, along with the true attenuation coefficient of the structure it is inside, and then we apply the method to measured data with sources and find good agreement. △ Less

Submitted 4 October, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 10 pages, 8 figures, accepted for publication in IEEE Transactions on Nuclear Science

Reconstructing the Position and Intensity of Multiple Gamma-Ray Point Sources with a Sparse Parametric Algorithm

Authors: Jayson R. Vavrek, Daniel Hellfeld, Mark S. Bandstra, Victor Negut, Kathryn Meehan, William J. Vanderlip, Joshua W. Cates, Ryan Pavlovsky, Brian J. Quiter, Reynold J. Cooper, Tenzing H. Y. Joshi

Abstract: We present an experimental demonstration of Additive Point Source Localization (APSL), a sparse parametric imaging algorithm that reconstructs the 3D positions and activities of multiple gamma-ray point sources. Using a handheld gamma-ray detector array and up to four $8$ $μ$Ci $^{137}$Cs gamma-ray sources, we performed both source-search and source-separation experiments in an indoor laboratory e… ▽ More We present an experimental demonstration of Additive Point Source Localization (APSL), a sparse parametric imaging algorithm that reconstructs the 3D positions and activities of multiple gamma-ray point sources. Using a handheld gamma-ray detector array and up to four $8$ $μ$Ci $^{137}$Cs gamma-ray sources, we performed both source-search and source-separation experiments in an indoor laboratory environment. In the majority of the source-search measurements, APSL reconstructed the correct number of sources with position accuracies of ${\sim}20$ cm and activity accuracies (unsigned) of ${\sim}20\%$, given measurement times of two to three minutes and distances of closest approach (to any source) of ${\sim}20$ cm. In source-separation measurements where the detector could be moved freely about the environment, APSL was able to resolve two sources separated by $75$ cm or more given only ${\sim}60$ s of measurement time. In these source-separation measurements, APSL produced larger total activity errors of ${\sim}40\%$, but obtained source separation distances accurate to within $15$ cm. We also compare our APSL results against traditional Maximum Likelihood-Expectation Maximization (ML-EM) reconstructions, and demonstrate improved image accuracy and interpretability using APSL over ML-EM. These results indicate that APSL is capable of accurately reconstructing gamma-ray source positions and activities using measurements from existing detector hardware. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 10 pages, 9 figures. Accepted for publication at IEEE Transactions on Nuclear Science

Origins Space Telescope Mission Concept Study Report

Authors: M. Meixner, A. Cooray, D. Leisawitz, J. Staguhn, L. Armus, C. Battersby, J. Bauer, E. Bergin, C. M. Bradford, K. Ennico-Smith, J. Fortney, T. Kataria, G. Melnick, S. Milam, D. Narayanan, D. Padgett, K. Pontoppidan, A. Pope, T. Roellig, K. Sandstrom, K. Stevenson, K. Su, J. Vieira, E. Wright, J. Zmuidzinas , et al. (44 additional authors not shown)

Abstract: The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and lar… ▽ More The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report. △ Less

Submitted 23 December, 2019; v1 submitted 12 December, 2019; originally announced December 2019.

Comments: 376 pages

Simultaneous LSST and Euclid observations - advantages for Solar System Objects

Authors: C. Snodgrass, B. Carry, J. Berthier, S. Eggl, M. Mommert, J. -M. Petit, F. Spoto, M. Granvik, R. Laureijs, B. Altieri, R. Vavrek, L. Conversi, A. Nucita, M. Popescu, G. Verdoes Kleijn, M. Kidger, G. H. Jones, D. Oszkiewicz, M. Juric, L. Jones

Abstract: The ESA Euclid mission is a space telescope that will survey ~15,000 square degrees of the sky, primarily to study the distant universe (constraining cosmological parameters through the lensing of galaxies). It is also expected to observe ~150,000 Solar System Objects (SSOs), primarily in poorly understood high inclination populations, as it will mostly avoid +/-15 degrees from the ecliptic plane.… ▽ More The ESA Euclid mission is a space telescope that will survey ~15,000 square degrees of the sky, primarily to study the distant universe (constraining cosmological parameters through the lensing of galaxies). It is also expected to observe ~150,000 Solar System Objects (SSOs), primarily in poorly understood high inclination populations, as it will mostly avoid +/-15 degrees from the ecliptic plane. With a launch date of 2022 and a 6 year survey, Euclid and LSST will operate at the same time, and have complementary capabilities. We propose a LSST mini-survey to coordinate quasi-simultaneous observations between these two powerful observatories, when possible, with the primary aim of greatly improving the orbits of SSOs discovered by these facilities. As Euclid will operate from a halo orbit around the Sun-Earth L2 Lagrangian point, there will be significant parallax between observations from Earth and Euclid (0.01 AU). This means that simultaneous observations will give an independent distance measurement to SSOs, giving additional constraints on orbits compared to single Euclid visits. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: white paper submitted as part of the LSST survey strategy call

Validation of Geant4's G4NRF module against nuclear resonance fluorescence data from $^{238}$U and $^{27}$Al

Authors: Jayson R Vavrek, Brian S Henderson, Areg Danagoulian

Abstract: G4NRF is a simulation module for modeling nuclear resonance fluorescence (NRF) interactions in the Geant4 framework. In this work, we validate G4NRF against both absolute and relative measurements of three NRF interactions near 2.2 MeV in $^{238}$U and $^{27}$Al using the transmission NRF data from the experiments described in arXiv:1712.02904. Agreement between the absolute NRF count rates observ… ▽ More G4NRF is a simulation module for modeling nuclear resonance fluorescence (NRF) interactions in the Geant4 framework. In this work, we validate G4NRF against both absolute and relative measurements of three NRF interactions near 2.2 MeV in $^{238}$U and $^{27}$Al using the transmission NRF data from the experiments described in arXiv:1712.02904. Agreement between the absolute NRF count rates observed in the data and predicted by extensive Geant4+G4NRF modeling validate the combined Geant4+G4NRF to within $15$--$20\%$ in the $^{238}$U NRF transitions and $8\%$ in $^{27}$Al, for an average $13\%$ discrepancy across the entire study. The difference between simulation and experiment in relative NRF rates, as expressed as ratios of count rates in various NRF lines, is found at the level of ${\lesssim}4\%$, and is statistically identical to zero. Inverting the analysis, approximate values of the absolute level widths and branching ratios for $^{238}$U and $^{27}$Al are also obtained. △ Less

Submitted 5 July, 2019; v1 submitted 6 July, 2018; originally announced July 2018.

Comments: 12 pages, 4 figures, 4 tables; revisions after peer review comments, chiefly making the paper more concise and the reporting of results more clear

Journal ref: NIM B 459 188-195 (2019)

High-accuracy Geant4 simulation and semi-analytical modeling of nuclear resonance fluorescence

Authors: Jayson R Vavrek, Brian S Henderson, Areg Danagoulian

Abstract: Nuclear resonance fluorescence (NRF) is a photonuclear interaction that enables highly isotope-specific measurements in both pure and applied physics scenarios. High-accuracy design and analysis of NRF measurements in complex geometries is aided by Monte Carlo simulations of photon physics and transport, motivating Jordan and Warren (2007) to develop the G4NRF codebase for NRF simulation in Geant4… ▽ More Nuclear resonance fluorescence (NRF) is a photonuclear interaction that enables highly isotope-specific measurements in both pure and applied physics scenarios. High-accuracy design and analysis of NRF measurements in complex geometries is aided by Monte Carlo simulations of photon physics and transport, motivating Jordan and Warren (2007) to develop the G4NRF codebase for NRF simulation in Geant4. In this work, we enhance the physics accuracy of the G4NRF code and perform improved benchmarking simulations. The NRF cross section calculation in G4NRF, previously a Gaussian approximation, has been replaced with a full numerical integration for improved accuracy in thick-target scenarios. A high-accuracy semi-analytical model of expected NRF count rates in a typical NRF measurement is then constructed and compared against G4NRF simulations for both simple homogeneous and more complex heterogeneous geometries. Agreement between rates predicted by the semi-analytical model and G4NRF simulation is found at a level of ${\sim}1\%$ in simple test cases and ${\sim}3\%$ in more realistic scenarios, improving upon the ${\sim}20\%$ level of the initial benchmarking study and establishing a highly-accurate NRF framework for Geant4. △ Less

Submitted 30 July, 2018; v1 submitted 3 July, 2018; originally announced July 2018.

Comments: 16 pages, 6 figures, revised for peer review

Journal ref: NIM B 433 34-42 (2018)

Experimental demonstration of an isotope-sensitive warhead verification technique using nuclear resonance fluorescence

Authors: Jayson R Vavrek, Brian S Henderson, Areg Danagoulian

Abstract: Future nuclear arms reduction efforts will require technologies to verify that warheads slated for dismantlement are authentic without revealing any sensitive weapons design information to international inspectors. Despite several decades of research, no technology has met these requirements simultaneously. Recent work by Kemp et al. [Kemp RS, Danagoulian A, Macdonald RR, Vavrek JR (2016) Proc Nat… ▽ More Future nuclear arms reduction efforts will require technologies to verify that warheads slated for dismantlement are authentic without revealing any sensitive weapons design information to international inspectors. Despite several decades of research, no technology has met these requirements simultaneously. Recent work by Kemp et al. [Kemp RS, Danagoulian A, Macdonald RR, Vavrek JR (2016) Proc Natl Acad Sci USA 113:8618--8623] has produced a novel physical cryptographic verification protocol that approaches this treaty verification problem by exploiting the isotope-specific nature of nuclear resonance fluorescence (NRF) measurements to verify the authenticity of a warhead. To protect sensitive information, the NRF signal from the warhead is convolved with that of an encryption foil that contains key warhead isotopes in amounts unknown to the inspector. The convolved spectrum from a candidate warhead is statistically compared against that from an authenticated template warhead to determine whether the candidate itself is authentic. Here we report on recent proof-of-concept warhead verification experiments conducted at the Massachusetts Institute of Technology. Using high-purity germanium (HPGe) detectors, we measured NRF spectra from the interrogation of proxy 'genuine' and 'hoax' objects by a 2.52 MeV endpoint bremsstrahlung beam. The observed differences in NRF intensities near 2.2 MeV indicate that the physical cryptographic protocol can distinguish between proxy genuine and hoax objects with high confidence in realistic measurement times. △ Less

Submitted 9 April, 2018; v1 submitted 7 December, 2017; originally announced December 2017.

Comments: 38 pages, 19 figures; revised for peer review and copy editing; addition to SI for realistic scenario projections; minor length reduction for journal requirements

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 115 (17) 4363-4368 (2018)

Star Formation Under the Outflow: The Discovery of a Non-Thermal Jet from OMC-2 FIR 3 and its Relationship to the Deeply Embedded FIR 4 Protostar

Authors: Mayra Osorio, Ana K. Diaz-Rodriguez, Guillem Anglada, S. Thomas Megeath, Luis F. Rodriguez, John J. Tobin, Amelia M. Stutz, Elise Furlan, William J. Fischer, P. Manoj, Jose F. Gomez, Beatriz Gonzalez-Garcia, Thomas Stanke, Dan M. Watson, Laurent Loinard, Roland Vavrek, Carlos Carrasco-Gonzalez

Abstract: We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Clas… ▽ More We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects. △ Less

Submitted 22 March, 2017; originally announced March 2017.

Comments: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

The Euclid mission design

Authors: Giuseppe D Racca, Rene Laureijs, Luca Stagnaro, Jean Christophe Salvignol, Jose Lorenzo Alvarez, Gonzalo Saavedra Criado, Luis Gaspar Venancio, Alex Short, Paolo Strada, Tobias Boenke, Cyril Colombo, Adriano Calvi, Elena Maiorano, Osvaldo Piersanti, Sylvain Prezelus, Pierluigi Rosato, Jacques Pinel, Hans Rozemeijer, Valentina Lesna, Paolo Musi, Marco Sias, Alberto Anselmi, Vincent Cazaubiel, Ludovic Vaillon, Yannick Mellier , et al. (17 additional authors not shown)

Abstract: Euclid is a space-based optical/near-infrared survey mission of the European Space Agency (ESA) to investigate the nature of dark energy, dark matter and gravity by observing the geometry of the Universe and on the formation of structures over cosmological timescales. Euclid will use two probes of the signature of dark matter and energy: Weak gravitational Lensing, which requires the measurement o… ▽ More Euclid is a space-based optical/near-infrared survey mission of the European Space Agency (ESA) to investigate the nature of dark energy, dark matter and gravity by observing the geometry of the Universe and on the formation of structures over cosmological timescales. Euclid will use two probes of the signature of dark matter and energy: Weak gravitational Lensing, which requires the measurement of the shape and photometric redshifts of distant galaxies, and Galaxy Clustering, based on the measurement of the 3-dimensional distribution of galaxies through their spectroscopic redshifts. The mission is scheduled for launch in 2020 and is designed for 6 years of nominal survey operations. The Euclid Spacecraft is composed of a Service Module and a Payload Module. The Service Module comprises all the conventional spacecraft subsystems, the instruments warm electronics units, the sun shield and the solar arrays. In particular the Service Module provides the extremely challenging pointing accuracy required by the scientific objectives. The Payload Module consists of a 1.2 m three-mirror Korsch type telescope and of two instruments, the visible imager and the near-infrared spectro-photometer, both covering a large common field-of-view enabling to survey more than 35% of the entire sky. All sensor data are downlinked using K-band transmission and processed by a dedicated ground segment for science data processing. The Euclid data and catalogues will be made available to the public at the ESA Science Data Centre. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: 23 pages, 19 figures, Presented at the SPIE Astronomical Telescopes and Instrumentation conference in Edinburgh, Scotland, United Kingdom, 6 June 1 July 2016

Herschel/PACS far-IR spectral imaging of a jet from an intermediate mass protostar in the OMC-2 region

Authors: B. Gonzalez-Garcia, P. Manoj, D. M. Watson, R. Vavrek, S. T. Megeath, A. M. Stutz, M. Osorio, F. Wyrowski, W. J. Fischer, J. J. Tobin, M. Sanchez-Portal, A. K. Diaz Rodriguez, T. L. Wilson

Abstract: We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region… ▽ More We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region (PDR) and a narrow, but intense jet connecting the two protostars. The jet seen in [O I] emission is spatially aligned with the Spitzer/IRAC 4.5 micron jet and the CO (6-5) molecular outflow centered on FIR 3. The mass loss rate derived from the total [O I] 63 micron line luminosity of the jet is 7.7 x 10^-6 M_sun/yr, more than an order of magnitude higher than that measured for typical low mass class 0 protostars. The implied accretion luminosity is significantly higher than the observed bolometric luminosity of FIR 4, indicating that the [O I] jet is unlikely to be associated with FIR 4. We argue that the peak line emission seen toward FIR 4 originates in the terminal shock produced by the jet driven by FIR 3. The higher mass-loss rate that we find for FIR 3 is consistent with the idea that intermediate mass protostars drive more powerful jets than their low-mass counterparts. Our results also call into question the nature of FIR 4. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: Accepted for publication in A & A

The Far-Infrared Surveyor Mission Study: Paper I, the Genesis

Authors: M. Meixner, A. Cooray, R. Carter, M. DiPirro, A. Flores, D. Leisawitz, L. Armus, C. Battersby, E. Bergin, C. M. Bradford, K. Ennico, G. J. Melnick, S. Milam, D. Narayanan, K. Pontoppidan, A. Pope, T. Roellig, K. Sandstrom, K. Y. L. Su, J. Vieira, E. Wright, J. Zmuidzinas, S. Alato, S. Carey, M. Gerin , et al. (5 additional authors not shown)

Abstract: This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to ma… ▽ More This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality. △ Less

Submitted 12 August, 2016; originally announced August 2016.

Comments: 8 pages, SPIE proceedings of the Space Telescopes and Instrumentation 2016: Optical, Infrared, and Millimeter Wave conference

Discovery of molecular gas around HD 131835 in an APEX molecular line survey of bright debris disks

Authors: A. Moór, Th. Henning, A. Juhász, P. Ábrahám, Z. Balog, Á. Kóspál, I. Pascucci, Gy. M. Szabó, R. Vavrek, M. Curé, T. Csengeri, C. Grady, R. Güsten, Cs. Kiss

Abstract: Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other indi… ▽ More Debris disks are considered to be gas-poor, but recent observations revealed molecular or atomic gas in several 10-40 Myr old systems. We used the APEX and IRAM 30m radiotelescopes to search for CO gas in 20 bright debris disks. In one case, around the 16 Myr old A-type star HD 131835, we discovered a new gas-bearing debris disk, where the CO 3-2 transition was successfully detected. No other individual system exhibited a measurable CO signal. Our Herschel Space Observatory far-infrared images of HD 131835 marginally resolved the disk both at 70 and 100$μ$m, with a characteristic radius of ~170 au. While in stellar properties HD 131835 resembles $β$ Pic, its dust disk properties are similar to those of the most massive young debris disks. With the detection of gas in HD 131835 the number of known debris disks with CO content has increased to four, all of them encircling young ($\leq$40 Myr) A-type stars. Based on statistics within 125 pc, we suggest that the presence of detectable amount of gas in the most massive debris disks around young A-type stars is a common phenomenon. Our current data cannot conclude on the origin of gas in HD 131835. If the gas is secondary, arising from the disruption of planetesimals, then HD 131835 is a comparably young and in terms of its disk more massive analogue of the $β$ Pic system. However, it is also possible that this system similarly to HD 21997 possesses a hybrid disk, where the gas material is predominantly primordial, while the dust grains are mostly derived from planetesimals. △ Less

Submitted 13 October, 2015; originally announced October 2015.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures, 5 tables

Probing the Isotropy in the Sky Distribution of the Gamma-Ray Bursts

Authors: A. Meszaros, Z. Bagoly, L. G. Balazs, I. Horvath, R. Vavrek

Abstract: The statistical tests - done by the authors - are surveyed, which verify the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at BATSE Catalog. The tests use the counts-in-cells method, an analysis of spherical harmonics, a test based on the two-point correlation function and a method based on multiscale methods. The tests suggest that the inter… ▽ More The statistical tests - done by the authors - are surveyed, which verify the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at BATSE Catalog. The tests use the counts-in-cells method, an analysis of spherical harmonics, a test based on the two-point correlation function and a method based on multiscale methods. The tests suggest that the intermediate subclass of gamma-ray bursts are distributed anisotropically. △ Less

Submitted 6 May, 2015; originally announced May 2015.

Comments: Gamma-Ray Bursts in the Afterglow Era: Proceedings of the International Workshop Held in Rome, Italy, 17-20 October 2000, ESO ASTROPHYSICS SYMPOSIA. ISBN 3-540-42771-6. Edited by E. Costa, F. Frontera, and J. Hjorth. Springer-Verlag, 2001, p. 47

Herschel Far-Infrared Spectral-mapping of Orion BN/KL Outflows: Spatial distribution of excited CO, H2O, OH, O and C+ in shocked gas

Authors: Javier R. Goicoechea, Luis Chavarria, Jose Cernicharo, David A. Neufeld, Roland Vavrek, Edwin A. Bergin, Sara Cuadrado, Pierre Encrenaz, Mireya Etxaluze, Gary J. Melnick, Edward Polehampton

Abstract: We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, mos… ▽ More We present ~2'x2' spectral-maps of Orion BN/KL outflows taken with Herschel at ~12'' resolution. For the first time in the far-IR domain, we spatially resolve the emission associated with the bright H2 shocked regions "Peak 1" and "Peak 2" from that of the Hot Core and ambient cloud. We analyze the ~54-310um spectra taken with the PACS and SPIRE spectrometers. More than 100 lines are detected, most of them rotationally excited lines of 12CO (up to J=48-47), H2O, OH, 13CO, and HCN. Peaks 1/2 are characterized by a very high L(CO)/L(FIR)~5x10^{-3} ratio and a plethora of far-IR H2O emission lines. The high-J CO and OH lines are a factor ~2 brighter toward Peak 1 whereas several excited H2O lines are ~50% brighter toward Peak 2. A simplified non-LTE model allowed us to constrain the dominant gas temperature components. Most of the CO column density arises from Tk~200-500 K gas that we associate with low-velocity shocks that fail to sputter grain ice mantles and show a maximum gas-phase H2O/CO~10^{-2} abundance ratio. In addition, the very excited CO (J>35) and H2O lines reveal a hotter gas component (Tk~2500 K) from faster (v_S>25 km/s) shocks that are able to sputter the frozen-out H2O and lead to high H2O/CO>~1 abundance ratios. The H2O and OH luminosities cannot be reproduced by shock models that assume high (undepleted) abundances of atomic oxygen in the preshock gas and/or neglect the presence of UV radiation in the postshock gas. Although massive outflows are a common feature in other massive star-forming cores, Orion BN/KL seems more peculiar because of its higher molecular luminosities and strong outflows caused by a recent explosive event. △ Less

Submitted 11 November, 2014; originally announced November 2014.

Comments: 21 pages, 11 figures. Accepted for publication in ApJ (english not edited)

Herschel's view of the large-scale structure in the Chamaeleon dark clouds

Authors: C. Alves de Oliveira, N. Schneider, B. Merín, T. Prusti, Á. Ribas, N. L. J. Cox, R. Vavrek, V. Könyves, D. Arzoumanian, E. Puga, G. L. Pilbratt, Á. Kóspál, Ph. André, P. Didelon, A. Men'shchikov, R. Royer, C. Waelkens, S. Bontemps, E. Winston, L. Spezzi

Abstract: The Chamaeleon molecular cloud complex is one of the nearest star-forming sites encompassing three molecular clouds with a different star-formation history, from quiescent (Cha III) to actively forming stars (Cha II), and reaching the end of star-formation (Cha I). To charactize its large-scale structure, we derived column density and temperature maps using PACS and SPIRE observations from the Her… ▽ More The Chamaeleon molecular cloud complex is one of the nearest star-forming sites encompassing three molecular clouds with a different star-formation history, from quiescent (Cha III) to actively forming stars (Cha II), and reaching the end of star-formation (Cha I). To charactize its large-scale structure, we derived column density and temperature maps using PACS and SPIRE observations from the Herschel Gould Belt Survey, and applied several tools, such as filament tracing, power-spectra, Δ-variance, and probability distribution functions of column density (PDFs), to derive physical properties. The column density maps reveal a different morphological appearance for the three clouds, with a ridge-like structure for Cha I, a clump-dominated regime for Cha II, and an intricate filamentary network for Cha III. The filament width is measured to be around 0.12\pm0.04 pc in the three clouds, and the filaments found to be gravitationally unstable in Cha I and II, but mostly subcritical in Cha III. Faint filaments (striations) are prominent in Cha I showing a preferred alignment with the large-scale magnetic field. The PDFs of all regions show a lognormal distribution at low column densities. For higher densities, the PDF of Cha I shows a turnover indicative of an extended higher density component, culminating with a power-law tail. Cha II shows a power-law tail with a slope characteristic of gravity. The PDF of Cha III can be best fit by a single lognormal. The turbulence properties of the three regions are found to be similar, pointing towards a scenario where the clouds are impacted by large-scale processes. The magnetic field could possibly play an important role for the star-formation efficiency in the Chamaeleon clouds if proven that it can effectively channel material on Cha I, and possibly Cha II, but probably less efficiently on the quiescent Cha III cloud. △ Less

Submitted 25 April, 2014; originally announced April 2014.

Comments: Accepted for publication in A&A

Journal ref: A&A 568, A98 (2014)

Characterizing the structure of diffuse emission in Hi-GAL maps

Authors: D. Elia, F. Strafella, N. Schneider, R. Paladini, R. Vavrek, Y. Maruccia, S. Molinari, S. Pezzuto, A. Noriega-Crespo, K. L. J. Rygl, A. Di Giorgio, A. Traficante, E. Schisano, L. Calzoletti, M. Pestalozzi, J. S. Liu, P. Natoli, M. Huang, P. Martin, Y. Fukui, T. Hayakawa

Abstract: We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this… ▽ More We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey carried out at five photometric bands from 70 to 500 μm constitutes an unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2°x2° Hi-GAL tiles located in the third Galactic quadrant (217° < l < 225°, -2° < b < 0°). The low level of confusion along the line of sight testified by CO observations makes this region an ideal case. We find very different values of the power spectrum slope from tile to tile but also from wavelength to wavelength (2 < β< 3), with similarities between fields attributable to components located at the same distance. Thanks to the comparison with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed. △ Less

Submitted 8 April, 2014; originally announced April 2014.

Comments: Accepted by ApJ

First results with the boloSource() algorithm: Photometry of faint standard stars observed by Herschel/PACS

Authors: Gábor Marton, Roland Vavrek, Csaba Kiss, Thomas G. Müller

Abstract: The boloSource() algorithm is a tool to separate the signal of compact sources from that of the diffuse background in the timeline of far-infrared measurements performed by the PACS camera of the Herschel Space Observatory. An important characteristic and quality indicator of this method is that how well it can reproduce the flux of faint standard stars which have reliable flux estimates. For this… ▽ More The boloSource() algorithm is a tool to separate the signal of compact sources from that of the diffuse background in the timeline of far-infrared measurements performed by the PACS camera of the Herschel Space Observatory. An important characteristic and quality indicator of this method is that how well it can reproduce the flux of faint standard stars which have reliable flux estimates. For this propose we selected a few calibrator targets and constructed light curves by extracting point source flux for each repetition of the measurements independently using standard aperture photometry methods. These were compared with the light curves obtained using the boloSource() method on the same dataset. The results indicate that boloSource() provides a similar level of photometric accuracy and reproducibility as the usual flux extraction and photometry methods. This new technique will be developed further and also tested against other methods in more complex fields with the goal to make it usable for large-scale studies in the future. △ Less

Submitted 20 October, 2013; originally announced October 2013.

Comments: 10 pages, 7 figures, 4 tables

A resolved debris disk around the candidate planet-hosting star HD95086

Authors: A. Moór, P. Ábrahám, Á. Kóspál, Gy. M. Szabó, D. Apai, Z. Balog, T. Csengeri, C. Grady, Th. Henning, A. Juhász, Cs. Kiss, I. Pascucci, J. Szulágyi, R. Vavrek

Abstract: Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, ß Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, h… ▽ More Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, ß Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, here we present new optical, infrared and millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory reveal a characteristic disk size of ~6.0x5.4 arcsec (540x490 AU) and disk inclination of ~25 degree. Assuming the same inclination for the planet candidate's orbit, its re-projected radial distance from the star is 62 AU, very close to the blackbody radius of the outer cold dust ring. The structure of the planetary system at HD95086 resembles the one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) between the two dusty regions. Modelling implies that the candidate planet can dynamically excite the motion of planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks co-exist. △ Less

Submitted 6 September, 2013; originally announced September 2013.

Comments: Accepted for publication in ApJ Letters, 6 pages, 2 figures, 3 tables (emulateapj style)

The Herschel Gould Belt Survey in Chamaeleon II - Properties of cold dust in disks around young stellar objects

Authors: L. Spezzi, N. L. J. Cox, T. Prusti, B. Merin, A. Ribas, C. Alves de Oliveira, E. Winston, A. Kospal, P. Royer, R. Vavrek, Ph. Andre, G. L. Pilbratt, L. Testi, E. Bressert, L. Ricci, A. Menshchikov, V. Konyves

Abstract: We report on the Herschel Gould Belt survey (HGBS) of ChaII, focusing on the detection of Class I to III young stellar objects (YSOs). We aim at characterizing the circumstellar material around these YSOs and understanding which disk parameters are most likely constrained by the new HGBS data. We recovered 29 out of the 63 known YSOs in ChaII with a detection in at least one of the PACS/SPIRE pass… ▽ More We report on the Herschel Gould Belt survey (HGBS) of ChaII, focusing on the detection of Class I to III young stellar objects (YSOs). We aim at characterizing the circumstellar material around these YSOs and understanding which disk parameters are most likely constrained by the new HGBS data. We recovered 29 out of the 63 known YSOs in ChaII with a detection in at least one of the PACS/SPIRE pass-bands: 3 Class I YSOs (i.e.,100%), 1 Flat source (i.e., 50%), 21 Class II objects (i.e., 55%), 3 Class III objects (i.e, 16%) and the unclassified source IRAS 12522-7640. We explore PACS/SPIRE colors of this sample and present modeling of their SEDs using the RADMC-2D radiative transfer code. We find that YSO colors are confined in specific regions of PACS/SPIRE color-color diagrams. These color ranges are expected to be only marginally contaminated by extragalactic sources and field stars and, hence, provide a useful YSO selection tool when applied altogether. We were able to model the SED of 26 out of the 29 detected YSOs. We discuss the degeneracy/limitations of our SED fitting results and adopt the Bayesian method to estimate the probability of different values for the derived disk parameters. The Cha II YSOs present typical disk inner radii around 0.1 AU. The lower limit to Rc is typically around 50 AU. The lower limits to Mdisk are proportional to the stellar masses with a typical 0.3% ratio, i.e., in the range estimated in the literature for young Class II stars and brown dwarfs across a broad range of stellar masses. The estimated flaring angles, although very uncertain, point towards rather flat disks (1+phi less than 1.2), as found for low-mass M-type YSO samples in other star forming regions. Thus, our results support the idea that disk properties show a dependence on stellar properties. △ Less

Submitted 17 May, 2013; originally announced May 2013.

Comments: 17 pages, 9 figures

Identification of transitional disks in Chamaeleon with Herschel

Authors: Á. Ribas, B. Merín, H. Bouy, C. Alves de Oliveira, D. R. Ardila, E. Puga, Á. Kóspál, L. Spezzi, N. L. J. Cox, T. Prusti, G. L. Pilbratt, Ph. André, L. Matrà, R. Vavrek

Abstract: Transitional disks are circumstellar disks with inner holes that in some cases are produced by planets and/or substellar companions in these systems. For this reason, these disks are extremely important for the study of planetary system formation. The Herschel Space Observatory provides an unique opportunity for studying the outer regions of protoplanetary disks. In this work we update previous kn… ▽ More Transitional disks are circumstellar disks with inner holes that in some cases are produced by planets and/or substellar companions in these systems. For this reason, these disks are extremely important for the study of planetary system formation. The Herschel Space Observatory provides an unique opportunity for studying the outer regions of protoplanetary disks. In this work we update previous knowledge on the transitional disks in the Chamaeleon I and II regions with data from the Herschel Gould Belt Survey. We propose a new method for transitional disk classification based on the WISE 12 micron-PACS 70 micron color, together with inspection of the Herschel images. We applied this method to the population of Class II sources in the Chamaeleon region and studied the spectral energy distributions of the transitional disks in the sample. We also built the median spectral energy distribution of Class II objects in these regions for comparison with transitional disks. The proposed method allows a clear separation of the known transitional disks from the Class II sources. We find 6 transitional disks, all previously known, and identify 5 objects previously thought to be transitional as possibly non-transitional. We find higher fluxes at the PACS wavelengths in the sample of transitional disks than those of Class II objects. We show the Herschel 70 micron band to be an efficient tool for transitional disk identification. The sensitivity and spatial resolution of Herschel reveals a significant contamination level among the previously identified transitional disk candidates for the two regions, which calls for a revision of previous samples of transitional disks in other regions. The systematic excess found at the PACS bands could be a result of the mechanism that produces the transitional phase, or an indication of different evolutionary paths for transitional disks and Class II sources. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: Accepted for publication in A&A: 11 March 2013 11 pages, 15 figures

Herschel/PACS Spectroscopic Survey of Protostars in Orion: The Origin of Far-Infrared CO Emission

Authors: P. Manoj, D. M. Watson, D. A. Neufeld, S. T. Megeath, R. Vavrek, Vincent Yu, R. Visser, E. A. Bergin, W. J. Fischer, J. J. Tobin, A. M. Stutz, B. Ali, T. L. Wilson, J. Di Francesco, M. Osorio, S. Maret, C. A. Poteet

Abstract: We present far-IR (57-196 mu) spectra of 21 protostars in the Orion molecular clouds, obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space observatory, as part of the Herschel Orion Protostar Survey (HOPS) program. We analyzed the CO emission lines (J_up = 14-46) in the PACS spectra, extracted within a projected distance of <= 2000 AU centered on the prot… ▽ More We present far-IR (57-196 mu) spectra of 21 protostars in the Orion molecular clouds, obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space observatory, as part of the Herschel Orion Protostar Survey (HOPS) program. We analyzed the CO emission lines (J_up = 14-46) in the PACS spectra, extracted within a projected distance of <= 2000 AU centered on the protostar. The total luminosity of the CO lines observed with PACS (L(CO)) is found to increase with increasing L_bol. The CO rotational temperature implied by the line ratios increases with J, and at least 3-4 rotational temperature components are required to fit the observed rotational diagram. The rotational temperature components are remarkably invariant between protostars and show no dependence on L_bol, T_bol or envelope density, implying that if the emitting gas is in LTE, the CO emission must arise in multiple temperature components that remain independent of L_bol over two orders of magnitudes. The observed CO emission can also be modeled as arising from a single temperature gas component or from a medium with a power-law temperature distribution; both of these require sub-thermally excited molecular gas at low densities (n(H_2) <= 10^6 cm^-3) and high temperatures (T >= 2000 K). Our results suggest that the contribution from PDRs along the envelope cavity walls is unlikely to be the dominant component of the CO emission observed with PACS. Instead, the "universality" of the rotational temperatures and the observed correlation between L(CO) and L_bol can most easily be explained if the observed CO emission originates in shock-heated, hot (T >= 2000 K), sub-thermally excited (n(H_2) <= 10^6 cm^-3) molecular gas. Post-shock gas at these densities is more likely to be found within the outflow cavities along the molecular outflow or along the cavity walls at radii >= several 100-1000 AU. △ Less

Submitted 9 November, 2012; originally announced November 2012.

Comments: 30 pages, 20 figures, 4 tables, accepted for publication in ApJ

A new Herschel view of the young star T54: not a transitional disk?

Authors: L. Matrà, B. Merín, C. Alves de Oliveira, N. Huélamo, A. Kóspál, N. L. J. Cox, Á. Ribas, E. Puga, R. Vavrek, P. Royer, T. Prusti, G. L. Pilbratt, P. André

Abstract: Context: Observations of transitional disks give us an understanding of the formation of planets and planetary systems such as our own. But care must be taken in the identification of such sources: the higher spatial resolution of the Herschel Space Observatory provides a new view on the origin of the far-infrared and sub-millimeter excesses observed. Aims: We review the nature of previously kno… ▽ More Context: Observations of transitional disks give us an understanding of the formation of planets and planetary systems such as our own. But care must be taken in the identification of such sources: the higher spatial resolution of the Herschel Space Observatory provides a new view on the origin of the far-infrared and sub-millimeter excesses observed. Aims: We review the nature of previously known transitional disks in the Chamaeleon I star-forming region with Herschel data. Methods: We analyze Herschel PACS and SPIRE images of the young star T54 together with ancillary images. We also analyze its spectral energy distribution and indications from optical and mid-infrared spectroscopy. Results: We detect extended emission in the PACS 70 μm image ~6" off source at a position angle of 196° from T54. The emission detected at longer wavelength (PACS 100, 160, SPIRE 250 and 350 μm) is also offset from the position of the star. This suggests that the excess observed in the far-infrared part of the SED is not fully associated with T54. Conclusions: Herschel images show that the far-infrared excess seen in T54 is not due to a transitional disk but to extended emission south-west of the source. The object still shows point-like and now downscaled excess at mid-infrared wavelengths, but its origin cannot be constrained without higher spatial resolution data. However, different indications point towards an evolved disk or extended unresolved emission close to the source. △ Less

Submitted 29 October, 2012; originally announced October 2012.

Comments: 5 pages, 2 figures, Accepted for publication in Astronomy & Astrophysics on 25 October 2012

The dusty environment of HD 97300 as seen by Herschel and Spitzer

Authors: Á. Kóspál, T. Prusti, N. L. J. Cox, G. L. Pilbratt, Ph. André, C. Alves de Oliveira, E. Winston, B. Merín, A. Ribas, P. Royer, R. Vavrek, C. Waelkens

Abstract: Aims. We analyze the surroundings of HD 97300, one of two intermediate-mass stars in the Chamaeleon I star-forming region. The star is known to be surrounded by a conspicuous ring of polycyclic aromatic hydrocarbons (PAHs). Methods. We present infrared images taken with Herschel and Spitzer using 11 different broad-band filters between 3.6 um and 500 um. We compare the morphology of the emission… ▽ More Aims. We analyze the surroundings of HD 97300, one of two intermediate-mass stars in the Chamaeleon I star-forming region. The star is known to be surrounded by a conspicuous ring of polycyclic aromatic hydrocarbons (PAHs). Methods. We present infrared images taken with Herschel and Spitzer using 11 different broad-band filters between 3.6 um and 500 um. We compare the morphology of the emission using cuts along different position angles. We construct spectral energy distributions, which we compare to different dust models, and calculate dust temperatures. We also derive opacity maps and analyze the density structure of the environment of HD 97300. Results. We find that HD 97300 has no infrared excess at or below 24 um, confirming its zero-age main-sequence nature. The morphology of the ring is very similar between 3.6 um and 24 um. The emission at these wavelengths is dominated by either PAH features or PAH continuum. At longer wavelengths, only the northwestern part of the ring is visible. A fit to the 100-500 um observations suggests that the emission is due to relatively warm (~26 K) dust. The temperature gradually decreases with increasing distance from the ring. We find a general decrease in the density from north to south, and an approximate 10% density increase in the northeastern part of the ring. Conclusions. Our results are consistent with the theory that the ring around HD 97300 is essentially a bubble blown into the surrounding interstellar matter and heated by the star. △ Less

Submitted 24 March, 2012; originally announced March 2012.

Comments: 6 pages, 3 figures, accepted for publication in A&A

The Herschel view of star formation in the Rosette molecular cloud under the influence of NGC 2244

Authors: N. Schneider, F. Motte, S. Bontemps, M. Hennemann, J. DiFrancesco, Ph. Andr, A. Zavagno, T. Csengeri, A. Men'shchikov, A. Abergel, J. -P. Baluteau, J. -Ph. Bernard, P. Cox, P. Didelon, A. -M diGiorgio, R. Gastaud, M. Griffin, P. Hargrave, T. Hill, M. Huang, J. Kirk, V. Konyves, S. Leeks, J. Z. Li, A. Marston , et al. (21 additional authors not shown)

Abstract: The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investig… ▽ More The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result. △ Less

Submitted 21 May, 2010; originally announced May 2010.

Comments: Letter in press for Astronomy & Astrophysics Herschel Special Issue

Small-scale structure in the Rosette molecular cloud revealed by Herschel

Authors: J. Di Francesco, S. Sadavoy, F. Motte, N. Schneider, M. Hennemann, S. Bontemps, T. Csengeri, Z. Balog, A. Zavagno, Ph. Andre, P. Saraceno, M. Griffin, A. Men'shchikov, A. Abergel, J. -P. Baluteau, J. -Ph. Bernard, P. Cox, L. Deharveng, P. Didelon, A. -M. di Giorgio, P. Hargrave, M. Huang, J. Kirk, S. Leeks, J. Z. Li , et al. (18 additional authors not shown)

Abstract: We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 cl… ▽ More We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups. △ Less

Submitted 20 May, 2010; originally announced May 2010.

Comments: 6 pages, to be published in Astronomy & Astrophysics

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

Authors: S. Molinari, B. Swinyard, J. Bally, M. Barlow, J. P. Bernard, P. Martin, T. Moore, A. Noriega-Crespo, R. Plume, L. Testi, A. Zavagno, A. Abergel, B. Ali, L. Anderson, P. André, J. P. Baluteau, C. Battersby, M. T. Beltrán, M. Benedettini, N. Billot, J. Blommaert, S. Bontemps, F. Boulanger, J. Brand, C. Brunt , et al. (99 additional authors not shown)

Abstract: We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structure… ▽ More We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows. △ Less

Submitted 18 May, 2010; originally announced May 2010.

Comments: A&A, accepted

The Herschel first look at protostars in the Aquila Rift

Authors: S. Bontemps, Ph. Andre, V. Konyves, A. Men'shchikov, N. Schneider, A. Maury, N. Peretto, D. Arzoumanian, M. Attard, F. Motte, V. Minier, P. Didelon, P. Saraceno, A. Abergel, J. -P. Baluteau, J. -Ph. Bernard, L. Cambresy, P. Cox, J. Di Francesco, A. M. Di Giorgo, M. Griffin, P. Hargrave, M. Huang, J. Kirk, J. Li , et al. (18 additional authors not shown)

Abstract: As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the col… ▽ More As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: 6 pages, 6 figures, publication in A&A

From filamentary clouds to prestellar cores to the stellar IMF: Initial highlights from the Herschel Gould Belt survey

Authors: Ph. André, A. Men'shchikov, S. Bontemps, V. Könyves, F. Motte, N. Schneider, P. Didelon, V. Minier, P. Saraceno, D. Ward-Thompson, J. Di Francesco, G. White, S. Molinari, L. Testi, A. Abergel, M. Griffin, Th. Henning, P. Royer, B. Merín, R. Vavrek, M. Attard, D. Arzoumanian, C. D. Wilson, P. Ade, H. Aussel , et al. (32 additional authors not shown)

Abstract: We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-… ▽ More We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability. △ Less

Submitted 25 May, 2010; v1 submitted 14 May, 2010; originally announced May 2010.

Comments: 7 pages, 4 figures. Accepted for publication in Astronomy and Astrophysics, Herschel Special Issue

The Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory

Authors: A. Poglitsch, C. Waelkens, N. Geis, H. Feuchtgruber, B. Vandenbussche, L. Rodriguez, O. Krause, E. Renotte, C. van Hoof, P. Saraceno, J. Cepa, F. Kerschbaum, P. Agnese, B. Ali, B. Altieri, P. Andreani, J. -L. Augueres, Z. Balog, L. Barl, O. H. Bauer, N. Belbachir, M. Benedettini, N. Billot, O. Boulade, H. Bischof , et al. (58 additional authors not shown)

Abstract: The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210… ▽ More The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210μ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85μ m or 85-125μ\m and 125-210μ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions. △ Less

Submitted 20 May, 2010; v1 submitted 10 May, 2010; originally announced May 2010.

New Statistical Results on the Angular Distribution of Gamma-Ray Bursts

Authors: L. G. Balazs, I. Horvath, R. Vavrek, Z. Bagoly, A. Meszaros

Abstract: We presented the results of several statistical tests of the randomness in the angular sky-distribution of gamma-ray bursts in BATSE Catalog. Thirteen different tests were presented based on Voronoi tesselation, Minimal spanning tree and Multifractal spectrum for five classes (short1, short2, intermediate, long1, long2) of gamma-ray bursts, separately. The long1 and long2 classes are distributed… ▽ More We presented the results of several statistical tests of the randomness in the angular sky-distribution of gamma-ray bursts in BATSE Catalog. Thirteen different tests were presented based on Voronoi tesselation, Minimal spanning tree and Multifractal spectrum for five classes (short1, short2, intermediate, long1, long2) of gamma-ray bursts, separately. The long1 and long2 classes are distributed randomly. The intermediate subclass, in accordance with the earlier results of the authors, is distributed non-randomly. Concerning the short subclass earlier statistical tests also suggested some departure from the random distribution, but not on a high enough confidence level. The new tests presented in this article suggest also non-randomness here. △ Less

Submitted 27 February, 2009; originally announced February 2009.

Comments: in GAMMA-RAY BURSTS 2007: Proceedings of the Santa Fe Conference

Testing the Randomness in the Sky-Distribution of Gamma-Ray Bursts

Authors: R. Vavrek, L. G. Balázs, A. Mészáros, I. Horváth, Z. Bagoly

Abstract: We studied the complete randomness of the angular distribution of gamma-ray bursts (GRBs) detected by BATSE. Since GRBs seem to be a mixture of objects of different physical nature we divided the BATSE sample into 5 subsamples (short1, short2, intermediate, long1, long2) based on their durations and peak fluxes and studied the angular distributions separately. We used three methods, Voronoi tess… ▽ More We studied the complete randomness of the angular distribution of gamma-ray bursts (GRBs) detected by BATSE. Since GRBs seem to be a mixture of objects of different physical nature we divided the BATSE sample into 5 subsamples (short1, short2, intermediate, long1, long2) based on their durations and peak fluxes and studied the angular distributions separately. We used three methods, Voronoi tesselation, minimal spanning tree and multifractal spectra to search for non-randomness in the subsamples. To investigate the eventual non-randomness in the subsamples we defined 13 test-variables (9 from the Voronoi tesselation, 3 from the minimal spanning tree and one from the multifractal spectrum). Assuming that the point patterns obtained from the BATSE subsamples are fully random we made Monte Carlo simulations taking into account the BATSE's sky-exposure function. The MC simulations enabled us to test the null hypothesis i.e. that the angular distributions are fully random. We tested the randomness by binomial test and introducing squared Euclidean distances in the parameter space of the test-variables. We concluded that the short1, short2 groups deviate significantly (99.90%, 99.98%) from the fully randomness in the distribution of the squared Euclidean distances but it is not the case at the long samples. At the intermediate group the squared Euclidean distances also give significant deviation (98.51%). △ Less

Submitted 24 June, 2008; originally announced June 2008.

Comments: 17 pages, 4 figures, accepted for publication in MNRAS

Redshifts of the Long Gamma-Ray Bursts

Authors: Z. Bagoly, I. Csabai, A. Meszaros, P. Meszaros, I. Horvath, L. G. Balazs, R. Vavrek

Abstract: The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma… ▽ More The low energy spectra of some gamma-ray bursts' show excess components beside the power-law dependence. The consequences of such a feature allows to estimate the gamma photometric redshift of the long gamma-ray bursts in the BATSE Catalog. There is good correlation between the measured optical and the estimated gamma photometric redshifts. The estimated redshift values for the long bright gamma-ray bursts are up to z=4, while for the the faint long bursts - which should be up to z=20 - the redshifts cannot be determined unambiguously with this method. The redshift distribution of all the gamma-ray bursts with known optical redshift agrees quite well with the BATSE based gamma photometric redshift distribution. △ Less

Submitted 6 April, 2007; originally announced April 2007.

Comments: published in Baltic Astronomy

Journal ref: Baltic Astron. 13 (2004) 227

HST/NICMOS observations of a proto-brown dwarf candidate

Authors: D. Apai, L. V. Toth, Th. Henning, R. Vavrek, Z. Kovacs, D. Lemke

Abstract: We present deep HST/NICMOS observations peering through the outflow cavity of the protostellar candidate IRAS 04381+2540 in the Taurus Molecular Cloud-1. A young stellar object as central source, a jet and a very faint and close (0.6") companion are identified. The primary and the companion have similar colours, consistent with strong reddening. We argue that the companion is neither a shock-exc… ▽ More We present deep HST/NICMOS observations peering through the outflow cavity of the protostellar candidate IRAS 04381+2540 in the Taurus Molecular Cloud-1. A young stellar object as central source, a jet and a very faint and close (0.6") companion are identified. The primary and the companion have similar colours, consistent with strong reddening. We argue that the companion is neither a shock-excited knot nor a background star. The colour/magnitude information predicts a substellar upper mass limit for the companion, but the final confirmation will require spectroscopic information. Because of its geometry, young age and its rare low-mass companion, this system is likely to provide a unique insight into the formation of brown dwarfs. △ Less

Submitted 23 February, 2005; originally announced February 2005.

Comments: Astronomy & Astrophysics Letters, in press; 4 pages, 2 figures

Gamma photometric redshifts for long gamma-ray bursts

Authors: Z. Bagoly, I. Csabai, A. Mészáros, P. Mészáros, I. Horváth, L. G. Balázs, R. Vavrek

Abstract: It is known that the soft tail of the gamma-ray bursts' spectra show excesses from the exact power-law dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the BATSE Catalog. A verification of these redshifts is obtained for the… ▽ More It is known that the soft tail of the gamma-ray bursts' spectra show excesses from the exact power-law dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the BATSE Catalog. A verification of these redshifts is obtained for the 8 GRB which have both BATSE DISCSC data and measured optical spectroscopic redshifts. There is good correlation between the measured and esti redshifts, and the average error is $Δz \approx 0.33$. The method is similar to the photometric redshift estimation of galaxies in the optical range, hence it can be called as "gamma photometric redshift estimation". The estimated redshifts for the long bright gamma-ray bursts are up to $z \simeq 4$. For the the faint long bursts - which should be up to $z \simeq 20$ - the redshifts cannot be determined unambiguously with this method. △ Less

Submitted 25 November, 2002; originally announced November 2002.

Comments: accepted in A&A, 7 pages incl. 7 figures

Journal ref: Astron.Astrophys. 398 (2003) 919-926

A very young star forming region detected by the ISOPHOT Serendipity Survey

Authors: O. Krause, D. Lemke, L. V. Toth, U. Klaas, M. Haas, M. Stickel, R. Vavrek

Abstract: We present a multi-wavelength study of the star forming region ISOSS J 20298+3559, which was identified by a cross-correlation of cold compact sources from the 170 micron ISOPHOT Serendipity Survey (ISOSS) database coinciding with objects detected by the MSX, 2MASS and IRAS infrared surveys. ISOSS J 20298+3559 is associated with a massive dark cloud complex (M ~ 760 M$_{\odot}$) and located in t… ▽ More We present a multi-wavelength study of the star forming region ISOSS J 20298+3559, which was identified by a cross-correlation of cold compact sources from the 170 micron ISOPHOT Serendipity Survey (ISOSS) database coinciding with objects detected by the MSX, 2MASS and IRAS infrared surveys. ISOSS J 20298+3559 is associated with a massive dark cloud complex (M ~ 760 M$_{\odot}$) and located in the Cygnus X giant molecular cloud. We derive a distance of 1800 pc on the basis of optical extinction data. The low average dust temperature (T ~ 16 K) and large mass (M ~ 120 M$_{\odot}$) of the dense inner part of the cloud, which has not been dispersed, indicates a recent begin of star formation. The youth of the region is supported by the early evolutionary stage of several pre- and protostellar objects discovered across the regio n: I) Two candidate Class 0 objects with masses of 8 and 3.5 M$_{\odot}$, II) a gravitationally bound, cold (T ~ 12 K) and dense (n(H$_{2}$) \~ 2 x 10$^{5}$ cm$^{-3}$) cloud core with a mass of 50 M$_{\odot}$ and III) a Herbig B2 star with a mass of 6.5 M$_{\odot}$ and a bolometric luminosity of 2200 L$_{\odot}$, showing evidence for ongoing accretion and a stellar age of less than 40000 years. The dereddened SED of the Herbig star is well reproduced by an accretion disc + star model. The externally heated cold cloud core is a good candidate for a massive pre-protostellar object. The star formation efficiency in the central cloud region is about 14 %. △ Less

Submitted 15 September, 2002; originally announced September 2002.

Comments: 14 pages, 14 figures, 7 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: Astron.Astrophys. 398 (2003) 1007-1020

A Remarkable Angular Distribution of the Intermediate Subclass of Gamma-Ray Bursts

Authors: A. Meszaros, Z. Bagoly, I. Horvath, L. G. Balazs, R. Vavrek

Abstract: In the article a test is developed, which allows to test the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at the Current BATSE Catalog. The method is a modified version of the well-known counts-in-cells test, and fully eliminates the non-uniform sky-exposure function of BATSE instrument. Applying this method to the case of all gamma-ray bu… ▽ More In the article a test is developed, which allows to test the null-hypothesis of the intrinsic randomness in the angular distribution of gamma-ray bursts collected at the Current BATSE Catalog. The method is a modified version of the well-known counts-in-cells test, and fully eliminates the non-uniform sky-exposure function of BATSE instrument. Applying this method to the case of all gamma-ray bursts no intrinsic non-randomness is found. The test also did not find intrinsic non-randomnesses for the short and long gamma-ray bursts, respectively. On the other hand, using the method to the new intermediate subclass of gamma-ray bursts, the null-hypothesis of the intrinsic randomness for 181 intermediate gamma-ray bursts is rejected on the 96.4% confidence level. Taking 92 dimmer bursts from this subclass itself, we obtain the surprising result: This "dim" subclass of the intermediate subclass has an intrinsic non-randomness on the 99.3% confidence level. On the other hand, the 89 "bright" bursts show no intrinsic non-randomness. △ Less

Submitted 9 May, 2000; v1 submitted 20 March, 2000; originally announced March 2000.

Comments: 14 pages, 2 figures, Accepted for publication in Astrophysical Journal, typo changes

Journal ref: Astrophys.J. 539 (2000) 98

Testing the Intrinsic Randomnesses in the Angular Distributions of Gamma-Ray Bursts

Authors: A. Meszaros, Z. Bagoly, I. Horvath, L. G. Balazs, R. Vavrek

Abstract: The counts-in-cells and the two-point angular correlation function method are used to test the randomnesses in the angular distributions of both the all gamma-ray bursts collected at BATSE Catalog, and also their three subclasses ("short", "intermediate", "long"). The methods elimate the non-zero sky-exposure function of BATSE instrument. Both tests suggest intrinsic non-randomnesses for the int… ▽ More The counts-in-cells and the two-point angular correlation function method are used to test the randomnesses in the angular distributions of both the all gamma-ray bursts collected at BATSE Catalog, and also their three subclasses ("short", "intermediate", "long"). The methods elimate the non-zero sky-exposure function of BATSE instrument. Both tests suggest intrinsic non-randomnesses for the intermediate subclass; for the remaining three cases only the correlation function method. The confidence levels are between 95% and 99.9%. Separating the GRBs into two parts ("dim half" and "bright half", respectively) we obtain the result that the "dim" half shows a non-randomness on the 99.3% confidence level from the counts-in-cells test. △ Less

Submitted 3 December, 1999; originally announced December 1999.

Comments: 5 pages, 1 figure, contribution in 5th Huntsville Gamma-Ray Burst Symposium Proceedings

On the existence of the intrinsic anisotropies in the angular distributions of gamma-ray bursts

Authors: A. Meszaros, Z. Bagoly, R. Vavrek

Abstract: This articles is concerned primarily with the intrinsic anisotropy in the angular distribution of 2281 gamma-ray bursts (GRBs) collected in Current BATSE Gamma-Ray Burst Catalog until the end of year 1998, and, second, with intrinsic anisotropies of three subclasses ("short", "intermediate", "long") of GRBs. Testing based on spherical harmonics of each class, in equatorial coordinates, is presen… ▽ More This articles is concerned primarily with the intrinsic anisotropy in the angular distribution of 2281 gamma-ray bursts (GRBs) collected in Current BATSE Gamma-Ray Burst Catalog until the end of year 1998, and, second, with intrinsic anisotropies of three subclasses ("short", "intermediate", "long") of GRBs. Testing based on spherical harmonics of each class, in equatorial coordinates, is presented. Because the sky exposure function of BATSE instrument is not dependent on the right ascension $α$, any non-zero spherical harmonic proportional either to $P_n^m(\sin δ) \sin m α$ or to $P_n^m(\sin δ) \cos mα$ with $m \neq 0$ ($δ$ is the declination), immediately indicates an intrinsic non-zero term. It is a somewhat surprising result that the "intermediate" subclass shows an intrinsic anisotropy at the $97\%$ significance level caused by the high non-zero $P_3^1(\sin δ) \sin m α$ harmonic. The remaining two subclasses, and the full sample of GRBs, remain isotropic. △ Less

Submitted 1 December, 1999; originally announced December 1999.

Comments: 6 pages. Astronomy and Astrophysics, accepted for publication

An intrinsic anisotropy in the angular distribution of gamma-ray bursts

Authors: L. G. Balazs, A. Meszaros, I. Horvath, R. Vavrek

Abstract: The anisotropy of the sky distribution of 2025 gamma-ray bursts (GRBs) collected in Current BATSE catalog is confirmed. It is shown that the quadrupole term being proportional to $\sim \sin 2b \sin l$ is non-zero with a probability 99.9%. The occurrence of this anisotropy term is then supported by the binomial test even with the probability 99.97%. It is also argued that this anisotropy cannot b… ▽ More The anisotropy of the sky distribution of 2025 gamma-ray bursts (GRBs) collected in Current BATSE catalog is confirmed. It is shown that the quadrupole term being proportional to $\sim \sin 2b \sin l$ is non-zero with a probability 99.9%. The occurrence of this anisotropy term is then supported by the binomial test even with the probability 99.97%. It is also argued that this anisotropy cannot be caused exclusively by instrumental effects due to the non-uniform sky exposure of BATSE instrument; there should exist also some intrinsic anisotropy in the angular distribution of GRBs. Separating GRBs into short and long subclasses, it is shown that the 251 short ones are distributed anisotropically, but the 681 long ones seem to be distributed still isotropically. The 2-sample Kolmogorov-Smirnov test shows that they are distributed differently with a 98.7% probability. The character of anisotropy suggests that the cosmological origin of short GRBs further holds, and there is no evidence for their Galactical origin. The work in essence contains the key ideas and results of a recently published paper (\cite{balazs}), to which the new result following from the 2-sample Kolmogorov-Smirnov test is added, too. △ Less

Submitted 22 September, 1999; originally announced September 1999.

Comments: 5 pages, Accepted in A&A

Journal ref: Astron.Astrophys.Suppl.Ser. 138 (1999) 417