Showing 201–250 of 391 results for author: Fischer, J

HST Scattered Light Imaging and Modeling of the Edge-on Protoplanetary Disk ESO-H$α$ 569

Authors: Schuyler G. Wolff, Marshall D. Perrin, Karl Stapelfeldt, Gaspard Duchene, Francois Menard, Deborah Padgett, Christophe Pinte, Laurent Pueyo, William J. Fischer

Abstract: We present new HST ACS observations and detailed models for a recently discovered edge-on protoplanetary disk around ESO H$α$ 569 (a low-mass T Tauri star in the Cha I star forming region). Using radiative transfer models we probe the distribution of the grains and overall shape of the disk (inclination, scale height, dust mass, flaring exponent and surface/volume density exponent) by model fittin… ▽ More We present new HST ACS observations and detailed models for a recently discovered edge-on protoplanetary disk around ESO H$α$ 569 (a low-mass T Tauri star in the Cha I star forming region). Using radiative transfer models we probe the distribution of the grains and overall shape of the disk (inclination, scale height, dust mass, flaring exponent and surface/volume density exponent) by model fitting to multiwavelength (F606W and F814W) HST observations together with a literature compiled spectral energy distribution. A new tool set was developed for finding optimal fits of MCFOST radiative transfer models using the MCMC code emcee to efficiently explore the high dimensional parameter space. It is able to self-consistently and simultaneously fit a wide variety of observables in order to place constraints on the physical properties of a given disk, while also rigorously assessing the uncertainties in those derived properties. We confirm that ESO H$α$ 569 is an optically thick nearly edge-on protoplanetary disk. The shape of the disk is well described by a flared disk model with an exponentially tapered outer edge, consistent with models previously advocated on theoretical grounds and supported by millimeter interferometry. The scattered light images and spectral energy distribution are best fit by an unusually high total disk mass (gas+dust assuming a ratio of 100:1) with a disk-to-star mass ratio of 0.16. △ Less

Submitted 9 November, 2017; originally announced November 2017.

Comments: Accepted for publication in ApJ

Stellar mass functions and implications for a variable IMF

Authors: M. Bernardi, R. K. Sheth, J. -L. Fischer, A. Meert, K. -H. Chae, H. Dominguez-Sanchez, M. Huertas-Company, F. Shankar, V. Vikram

Abstract: Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with… ▽ More Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $σ_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $σ$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm α_{JAM}}$, from $M_{*}$ and $σ_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,σ_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(σ_R/σ_e) = (R/R_e)^{-γ(n)}$. The resulting stellar mass functions, based on $M_*^{\rm α_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm α_{JAM}}/M_* - σ_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $φ(M_*^{\rm α_{JAM}})$ is an estimate of the stellar mass function in which $σ_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $σ_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form. △ Less

Submitted 15 December, 2017; v1 submitted 19 October, 2017; originally announced October 2017.

Comments: 17 pages, 19 figures, 4 tables. Accepted for publication by MNRAS. Tables 1, C1 and C2 are provided as ancillary files

Feedback and feeding in the context of galaxy evolution with SPICA: direct characterization of molecular outflows and inflows

Authors: E. González-Alfonso, L. Armus, F. J. Carrera, V. Charmandaris, A. Efstathiou, E. Egami, J. A. Fernández-Ontiveros, J. Fischer, G. L. Granato, C. Gruppioni, E. Hatziminaoglou, M. Imanishi, N. Isobe, H. Kaneda, D. Koziel-Wierzbowska, M. A. Malkan, J. Martin-Pintado, S. Mateos, H. Matsuhara, G. Miniutti, T. Nakagawa, F. Pozzi, F. Rico-Villas, G. Rodighiero, P. Roelfsema , et al. (6 additional authors not shown)

Abstract: A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas.… ▽ More A far-infrared observatory such as the {\it SPace Infrared telescope for Cosmology and Astrophysics} ({\it SPICA}), with its unprecedented spectroscopic sensitivity, would unveil the role of feedback in galaxy evolution during the last $\sim10$ Gyr of the Universe ($z=1.5-2$), through the use of far- and mid-infrared molecular and ionic fine structure lines that trace outflowing and infalling gas. Outflowing gas is identified in the far-infrared through P-Cygni line shapes and absorption blueshifted wings in molecular lines with high dipolar moments, and through emission line wings of fine-structure lines of ionized gas. We quantify the detectability of galaxy-scale massive molecular and ionized outflows as a function of redshift in AGN-dominated, starburst-dominated, and main-sequence galaxies, explore the detectability of metal-rich inflows in the local Universe, and describe the most significant synergies with other current and future observatories that will measure feedback in galaxies via complementary tracers at other wavelengths. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: This paper belongs to the SPICA Special Issue on PASA. Accepted for publication in PASA

Galaxy evolution studies with the SPace IR telescope for Cosmology and Astrophysics (SPICA): the power of IR spectroscopy

Authors: L. Spinoglio, A. Alonso-Herrero, L. Armus, M. Baes, J. Bernard-Salas, S. Bianchi, M. Bocchio, A. Bolatto, C. M. Bradford, J. Braine, F. J. Carrera, L. Ciesla, D. L. Clements, H. Dannerbauer, Y. Doi, A. Efstathiou, E. Egami, J. A. Fernandez-Ontiveros, A. Ferrara, J. Fischer, A. Franceschini, S. Gallerani, M. Giard, E. Gonzalez-Alfonso, C. Gruppioni , et al. (43 additional authors not shown)

Abstract: IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Larg… ▽ More IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue

SPICA and the Chemical Evolution of Galaxies: The Rise of Metals and Dust

Authors: J. A. Fernández-Ontiveros, L. Armus, M. Baes, J. Bernard-Salas, A. D. Bolatto, J. Braine, L. Ciesla, I. De Looze, E. Egami, J. Fischer, M. Giard, E. González-Alfonso, G. L. Granato, C. Gruppioni, M. Imanishi, D. Ishihara, H. Kaneda, S. Madden, M. Malkan, H. Matsuhara, M. Matsuura, T. Nagao, F. Najarro, T. Nakagawa, T. Onaka , et al. (11 additional authors not shown)

Abstract: The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryog… ▽ More The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: In press. This paper belongs to the "SPICA Special Issue" on PASA

Dismantling DivSufSort

Authors: Johannes Fischer, Florian Kurpicz

Abstract: We give the first concise description of the fastest known suffix sorting algorithm in main memory, the DivSufSort by Yuta Mori. We then present an extension that also computes the LCP-array, which is competitive with the fastest known LCP-array construction algorithm. We give the first concise description of the fastest known suffix sorting algorithm in main memory, the DivSufSort by Yuta Mori. We then present an extension that also computes the LCP-array, which is competitive with the fastest known LCP-array construction algorithm. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: Presented at the Prague Stringology Conference 2017

Spin Hall magnetoresistance in antiferromagnet/heavy-metal heterostructures

Authors: Johanna Fischer, Olena Gomonay, Richard Schlitz, Kathrin Ganzhorn, Nynke Vlietstra, Matthias Althammer, Hans Huebl, Matthias Opel, Rudolf Gross, Sebastian T. B. Goennenwein, Stephan Geprägs

Abstract: We investigate the spin Hall magnetoresistance in thin film bilayer heterostructures of the heavy metal Pt and the antiferromagnetic insulator NiO. While rotating an external magnetic field in the easy plane of NiO, we record the longitudinal and the transverse resistivity of the Pt layer and observe an amplitude modulation consistent with the spin Hall magnetoresistance. In comparison to Pt on co… ▽ More We investigate the spin Hall magnetoresistance in thin film bilayer heterostructures of the heavy metal Pt and the antiferromagnetic insulator NiO. While rotating an external magnetic field in the easy plane of NiO, we record the longitudinal and the transverse resistivity of the Pt layer and observe an amplitude modulation consistent with the spin Hall magnetoresistance. In comparison to Pt on collinear ferrimagnets, the modulation is phase shifted by 90° and its amplitude strongly increases with the magnitude of the magnetic field. We explain the observed magnetic field-dependence of the spin Hall magnetoresistance in a comprehensive model taking into account magnetic field induced modifications of the domain structure in antiferromagnets. With this generic model we are further able to estimate the strength of the magnetoelastic coupling in antiferromagnets. Our detailed study shows that the spin Hall magnetoresistance is a versatile tool to investigate the magnetic spin structure as well as magnetoelastic effects, even in antiferromagnetic multidomain materials. △ Less

Submitted 15 January, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Journal ref: Phys. Rev. B 97, 014417 (2018)

Review: Far-Infrared Instrumentation and Technology Development for the Next Decade

Authors: Duncan Farrah, Kimberly Ennico Smith, David Ardila, Charles M. Bradford, Michael Dipirro, Carl Ferkinhoff, Jason Glenn, Paul Goldsmith, David Leisawitz, Thomas Nikola, Naseem Rangwala, Stephen A. Rinehart, Johannes Staguhn, Michael Zemcov, Jonas Zmuidzinas, James Bartlett, Sean Carey, William J. Fischer, Julia Kamenetzky, Jeyhan Kartaltepe, Mark Lacy, Dariusz C. Lis, Lisa Locke, Enrique Lopez-Rodriguez, Meredith MacGregor , et al. (11 additional authors not shown)

Abstract: Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the f… ▽ More Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century. △ Less

Submitted 8 January, 2019; v1 submitted 7 September, 2017; originally announced September 2017.

Comments: Invited review article, submitted to the Journal of Astronomical Telescopes, Instruments, and Systems. Updated in light of referee reports

Journal ref: J. Astron. Telesc. Instrum. Syst. 5(2), 020901 (2019)

ISM conditions in z~0.2 Lyman-Break Analogs

Authors: A. Contursi, A. J. Baker, S. Berta, B. Magnelli, D. Lutz, J. Fischer, A. Verma, M. Nielbock, J. Grácia Carpio, S. Veilleux, E. Sturm, R. Davies, R. Genzel, S. Hailey-Dunsheath, R. Herrera-Camus, A. Janssen, A. Poglitsch, A. Sternberg, L. J. Tacconi

Abstract: We present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar mass… ▽ More We present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar masses between $10^{9.5}$ and $10^{11}$ $M_{\odot}$, which are currently not easily detectable in all these lines even with ALMA and NOEMA. We perform PDR modeling and apply different IR diagnostics to derive the main physical parameters of the FIR emitting gas and dust and we compare the derived ISM properties to those of galaxies on and above the MS at different redshifts. We find that the ISM properties of LBAs are quite extreme (low gas temperature, high density and thermal pressure) with respect to those found in local normal spirals and more active local galaxies. LBAs have no [CII] deficit despite having the high specific star formation rates (sSFRs) typical of starbursts. Although LBAs lie above the local MS, we show that their ISM properties are more similar to those of high-redshift MS galaxies than of local galaxies above the main sequence. This data set represents an important reference for planning future ALMA [CII] observations of relatively low-mass MS galaxies at the epoch of the peak of the cosmic star formation. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 19 pages, 12 Figures,8 Tables, Accepted for publication in A&A

Practical Evaluation of Lempel-Ziv-78 and Lempel-Ziv-Welch Tries

Authors: Johannes Fischer, Dominik Köppl

Abstract: We present the first thorough practical study of the Lempel-Ziv-78 and the Lempel-Ziv-Welch computation based on trie data structures. With a careful selection of trie representations we can beat well-tuned popular trie data structures like Judy, m-Bonsai or Cedar. We present the first thorough practical study of the Lempel-Ziv-78 and the Lempel-Ziv-Welch computation based on trie data structures. With a careful selection of trie representations we can beat well-tuned popular trie data structures like Judy, m-Bonsai or Cedar. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Quasar Feedback in the Ultraluminous Infrared Galaxy F11119+3257: Connecting the Accretion Disk Wind with the Large-Scale Molecular Outflow

Authors: S. Veilleux, A. Bolatto, F. Tombesi, M. Melendez, E. Sturm, E. Gonzalez-Alfonso, J. Fischer, D. S. N. Rupke

Abstract: In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback mode… ▽ More In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the energetics were estimated from the optically thick OH 119 um transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of broad wings in the CO(1-0) profile derived from ALMA observations. The broad CO(1-0) line emission appears to be spatially extended on a scale of at least ~7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80-200) R_7^{-1} M_sun/yr, (1.5-3.0) R_7^{-1} L_AGN/c, and (0.15-0.40)% R_7^{-1} L_AGN are inferred from these data, assuming a CO-to-H_2 conversion factor appropriate for a ULIRG (R_7 is the radius of the outflow normalized to 7 kpc and L_AGN is the AGN luminosity). These rates are time-averaged over a flow time scale of 7x10^6 yrs. They are similar to the OH-based rates time-averaged over a flow time scale of 4x10^5 yrs, but about a factor 4 smaller than the local ("instantaneous"; <10^5 yrs) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also re-examined. △ Less

Submitted 1 June, 2017; originally announced June 2017.

Comments: 15 pages, 6 figures, 4 tables, accepted for publication in ApJ

The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

Authors: William J. Fischer, S. Thomas Megeath, Elise Furlan, Babar Ali, Amelia M. Stutz, John J. Tobin, Mayra Osorio, Thomas Stanke, P. Manoj, Charles A. Poteet, Joseph J. Booker, Lee Hartmann, Thomas L. Wilson, Philip C. Myers, Dan M. Watson

Abstract: The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of… ▽ More The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: Accepted by ApJ, 16 pages, 14 figures

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

Weak-strong uniqueness of solutions to entropy-dissipating reaction-diffusion equations

Authors: Julian Fischer

Abstract: We establish a weak-strong uniqueness principle for solutions to entropy-dissipating reaction-diffusion equations: As long as a strong solution to the reaction-diffusion equation exists, any weak solution and even any renormalized solution must coincide with this strong solution. Our assumptions on the reaction rates are just the entropy condition and local Lipschitz continuity; in particular, we… ▽ More We establish a weak-strong uniqueness principle for solutions to entropy-dissipating reaction-diffusion equations: As long as a strong solution to the reaction-diffusion equation exists, any weak solution and even any renormalized solution must coincide with this strong solution. Our assumptions on the reaction rates are just the entropy condition and local Lipschitz continuity; in particular, we do not impose any growth restrictions on the reaction rates. Therefore, our result applies to any single reversible reaction with mass-action kinetics as well as to systems of reversible reactions with mass-action kinetics satisfying the detailed balance condition. Renormalized solutions are known to exist globally in time for reaction-diffusion equations with entropy-dissipating reaction rates; in contrast, the global-in-time existence of weak solutions is in general still an open problem - even for smooth data - , thereby motivating the study of renormalized solutions. The key ingredient of our result is a careful adjustment of the usual relative entropy functional, whose evolution cannot be controlled properly for weak solutions or renormalized solutions. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: 32 pages

Comparing PyMorph and SDSS photometry. II. The differences are more than semantics and are not dominated by intracluster light

Authors: M. Bernardi, J. -L. Fischer, R. K. Sheth, A. Meert, M. Huertas-Company, F. Shankar, V. Vikram

Abstract: The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMor… ▽ More The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects. △ Less

Submitted 20 March, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: 15 pages, 17 figures, accepted for publication in MNRAS. The PyMorph luminosities and stellar masses are available at https://www.physics.upenn.edu/~ameert/SDSS_PhotDec/

Comparing PyMorph and SDSS photometry. I. Background sky and model fitting effects

Authors: J. -L. Fischer, M. Bernardi, A. Meert

Abstract: A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other tw… ▽ More A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other two are less trivial and tend to be larger; they are due to differences in how the background sky is estimated and what model is fit to the surface brightness profile. We show that PyMorph sky estimates are fainter than those of the SDSS DR7 or DR9 pipelines, but are in excellent agreement with the estimates of Blanton et al. (2011). Using the SDSS sky biases luminosities by more than a few tenths of a magnitude for objects with half-light radii > 7 arcseconds. In the SDSS main galaxy sample these are typically luminous galaxies, so they are not necessarily nearby. This bias becomes worse when allowing the model more freedom to fit the surface brightness profile. When PyMorph sky values are used, then two component Sersic-Exponential fits to E+S0s return more light than single component deVaucouleurs fits (up to ~0.2 mag), but less light than single Sersic fits (0.1 mag). Finally, we show that PyMorph fits of Meert et al. (2015) to DR7 data remain valid for DR9 images. Our findings show that, especially at large luminosities, these PyMorph estimates should be preferred to the SDSS pipeline values. △ Less

Submitted 27 February, 2017; originally announced February 2017.

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

Simple, Fast and Lightweight Parallel Wavelet Tree Construction

Authors: Johannes Fischer, Florian Kurpicz, Marvin Löbel

Abstract: The wavelet tree (Grossi et al. [SODA, 2003]) and wavelet matrix (Claude et al. [Inf. Syst., 47:15--32, 2015]) are compact indices for texts over an alphabet $[0,σ)$ that support rank, select and access queries in $O(\lg σ)$ time. We first present new practical sequential and parallel algorithms for wavelet tree construction. Their unifying characteristics is that they construct the wavelet tree b… ▽ More The wavelet tree (Grossi et al. [SODA, 2003]) and wavelet matrix (Claude et al. [Inf. Syst., 47:15--32, 2015]) are compact indices for texts over an alphabet $[0,σ)$ that support rank, select and access queries in $O(\lg σ)$ time. We first present new practical sequential and parallel algorithms for wavelet tree construction. Their unifying characteristics is that they construct the wavelet tree bottomup}, i.e., they compute the last level first. We also show that this bottom-up construction can easily be adapted to wavelet matrices. In practice, our best sequential algorithm is up to twice as fast as the currently fastest sequential wavelet tree construction algorithm (Shun [DCC, 2015]), simultaneously saving a factor of 2 in space. This scales up to 32 cores, where we are about equally fast as the currently fastest parallel wavelet tree construction algorithm (Labeit et al. [DCC, 2016]), but still use only about 75 % of the space. An additional theoretical result shows how to adapt any wavelet tree construction algorithm to the wavelet matrix in the same (asymptotic) time, using only little extra space. △ Less

Submitted 10 November, 2017; v1 submitted 24 February, 2017; originally announced February 2017.

Compression with the tudocomp Framework

Authors: Patrick Dinklage, Johannes Fischer, Dominik Köppl, Marvin Löbel, Kunihiko Sadakane

Abstract: We present a framework facilitating the implementation and comparison of text compression algorithms. We evaluate its features by a case study on two novel compression algorithms based on the Lempel-Ziv compression schemes that perform well on highly repetitive texts. We present a framework facilitating the implementation and comparison of text compression algorithms. We evaluate its features by a case study on two novel compression algorithms based on the Lempel-Ziv compression schemes that perform well on highly repetitive texts. △ Less

Submitted 22 April, 2021; v1 submitted 24 February, 2017; originally announced February 2017.

Journal ref: Proc. SEA, LIPIcs 75, pages 13:1-13:22, 2017

Molecular outflows in local ULIRGs: energetics from multi-transition OH analysis

Authors: E. González-Alfonso, J. Fischer, H. W. W. Spoon, K. P. Stewart, M. L. N. Ashby, S. Veilleux, H. A. Smith, E. Sturm, D. Farrah, N. Falstad, M. Meléndez, J. Graciá-Carpio, A. W. Janssen, V. Lebouteiller

Abstract: We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to mo… ▽ More We report on the energetics of molecular outflows in 14 local Ultraluminous Infrared Galaxies (ULIRGs) that show unambiguous outflow signatures (P-Cygni profiles or high-velocity absorption wings) in the far-infrared lines of OH measured with the Herschel/PACS spectrometer. Detection of both ground-state (at 119 and 79 um) and one or more radiatively-excited (at 65 and 84 um) lines allows us to model the nuclear gas (<~300 pc) as well as the more extended components using spherically symmetric radiative transfer models. The highest molecular outflow velocities are found in buried sources, in which slower but massive expansion of the nuclear gas is also observed. With the exception of a few outliers, the outflows have momentum fluxes of (2-5)xL_IR/c and mechanical luminosities of (0.1-0.3)% of L_IR. The moderate momentum boosts in these sources (<~3) suggest that the outflows are mostly momentum-driven by the combined effects of AGN and nuclear starbursts, as a result of radiation pressure, winds, and supernovae remnants. In some sources (~20%), however, powerful (10^{10.5-11} Lsun) AGN feedback and (partially) energy-conserving phases are required, with momentum boosts in the range 3-20. These outflows appear to be stochastic strong-AGN feedback events that occur throughout the merging process. In a few sources, the outflow activity in the innermost regions has subsided in the last ~1 Myr. While OH traces the molecular outflows at sub-kpc scales, comparison of the masses traced by OH with those previously inferred from tracers of more extended outflowing gas suggests that most mass is loaded (with loading factors of Mdot/SFR=1-10) from the central galactic cores (a few x 100 pc). Outflow depletion timescales are <10^8 yr, shorter than the gas consumption timescales by factors of 1.1-15, and are anti-correlated with the AGN luminosity. △ Less

Submitted 11 February, 2017; v1 submitted 24 December, 2016; originally announced December 2016.

Journal ref: The Astrophysical Journal, 2017, 836, 11

Was 49b: An Overmassive AGN in a Merging Dwarf Galaxy?

Authors: Nathan J. Secrest, Henrique R. Schmitt, Laura Blecha, Barry Rothberg, Jacqueline Fischer

Abstract: We present a combined morphological and X-ray analysis of Was 49, an isolated, dual AGN system notable for the presence of a dominant AGN Was 49b in the disk of the primary galaxy Was 49a, at a projected radial distance of 8 kpc from the nucleus. Using X-ray data from Chandra, NuSTAR, and Swift, we find that this AGN has a bolometric luminosity of L_bol ~ 2 x 10^45 erg/s, with a black hole mass of… ▽ More We present a combined morphological and X-ray analysis of Was 49, an isolated, dual AGN system notable for the presence of a dominant AGN Was 49b in the disk of the primary galaxy Was 49a, at a projected radial distance of 8 kpc from the nucleus. Using X-ray data from Chandra, NuSTAR, and Swift, we find that this AGN has a bolometric luminosity of L_bol ~ 2 x 10^45 erg/s, with a black hole mass of M_BH=1.3^{+2.9}_{-0.9} x 10^8 M_Sol. Despite its large mass, our analysis of optical data from the Discovery Channel Telescope shows that the supermassive black hole is hosted by a stellar counterpart with a mass of only 5.6^{+4.9}_{-2.6} x 10^9 M_Sol, making the SMBH potentially larger than expected from SMBH-galaxy scaling relations, and the stellar counterpart exhibits a morphology that is consistent with dwarf elliptical galaxies. Our analysis of the system in the r and K bands indicates that Was 49 is a minor merger, with a mass ratio of Was 49a to Was 49b between 1:7 and 1:15. This is in contrast with findings that the most luminous merger-triggered AGNs are found in major mergers, and that minor mergers predominantly enhance AGN activity in the primary galaxy. △ Less

Submitted 9 December, 2016; originally announced December 2016.

Comments: 11 pages, 6 figures, accepted for publication in ApJ

Inflowing gas onto a compact obscured nucleus in Arp 299A: Herschel spectroscopic studies of H2O and OH

Authors: N. Falstad, E. González-Alfonso, S. Aalto, J. Fischer

Abstract: Aims. We probe the physical conditions in the core of Arp 299A and try to put constraints to the nature of its nuclear power source. Methods. We used Herschel Space Observatory far-infrared and submillimeter observations of H2O and OH rotational lines in Arp 299A to create a multi-component model of the galaxy. In doing this, we employed a spherically symmetric radiative transfer code. Results. Ni… ▽ More Aims. We probe the physical conditions in the core of Arp 299A and try to put constraints to the nature of its nuclear power source. Methods. We used Herschel Space Observatory far-infrared and submillimeter observations of H2O and OH rotational lines in Arp 299A to create a multi-component model of the galaxy. In doing this, we employed a spherically symmetric radiative transfer code. Results. Nine H2O lines in absorption and eight in emission as well as four OH doublets in absorption and one in emission, are detected in Arp 299A. No lines of the 18O isotopologues, which have been seen in compact obscured nuclei of other galaxies, are detected. The absorption in the ground state OH doublet at 119 μm is found redshifted by ~175 km/s compared to other OH and H2O lines, suggesting a low excitation inflow. We find that at least two components are required in order to account for the excited molecular line spectrum. The inner component has a radius of 20-25 pc, a very high infrared surface brightness (> 3e13 Lsun/kpc^2), warm dust (Td > 90 K), and a large H2 column density (NH2 > 1e24 cm^-2). The outer component is larger (50-100 pc) with slightly cooler dust (70-90 K). In addition, a much more extended inflowing component is required to also account for the OH doublet at 119 μm. Conclusions. The Compton-thick nature of the core makes it difficult to determine the nature of the buried power source, but the high surface brightness indicates that it is either an active galactic nucleus and/or a dense nuclear starburst. The high OH/H2O ratio in the nucleus indicates that ion-neutral chemistry induced by X-rays or cosmic-rays is important. Finally we find a lower limit to the 16O/18O ratio of 400 in the nuclear region, possibly indicating that the nuclear starburst is in an early evolutionary stage, or that it is fed through a molecular inflow of, at most, solar metallicity. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 14 pages, 13 figures, Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 597, A105 (2017)

A mechanical model for phase-separation in debris flow

Authors: Shiva P. Pudasaini, Jan-Thomas Fischer

Abstract: Understanding the physics of phase-separation between solid and fluid phases as a mixture mass moves down slope is a long-standing challenge. Here, we propose an extension of the two phase mass flow model (Pudasaini, 2012) by including a new mechanism, called separation-flux, that leads to strong phase-separation in avalanche and debris flows while balancing the enhanced solid flux with the reduce… ▽ More Understanding the physics of phase-separation between solid and fluid phases as a mixture mass moves down slope is a long-standing challenge. Here, we propose an extension of the two phase mass flow model (Pudasaini, 2012) by including a new mechanism, called separation-flux, that leads to strong phase-separation in avalanche and debris flows while balancing the enhanced solid flux with the reduced fluid flux. The separation flux mechanism is capable of describing the dynamically evolving phase-separation and levee formation in a multi-phase, geometrically three-dimensional debris flow. These are often observed phenomena in natural debris flows and industrial processes that involve the transportation of particulate solid-fluid mixture material. The novel separation-flux model includes several dominant physical and mechanical aspects such as pressure gradients, volume fractions of solid and fluid phases and their gradients, shear-rates, flow depth, material friction, viscosity, material densities, topographic constraints, grain size, etc. Due to the inherent separation mechanism, as the mass moves down slope, more and more solid particles are transported to the front and the sides, resulting in solid-rich and mechanically strong frontal surge head, and lateral levees followed by a weaker tail largely consisting of viscous fluid. The primary frontal solid-rich surge head followed by secondary fluid-rich surges is the consequence of phase-separation. Such typical and dominant phase-separation phenomena are revealed for two-phase debris flow simulations. Finally, changes in flow composition, that are explicitly considered by the new modelling approach, result in significant changes of impact pressure estimates. These are highly important in hazard assessment and mitigation planning and highlight the application potential of the new approach. △ Less

Submitted 12 October, 2016; originally announced October 2016.

Journal ref: International Journal of Multiphase Flow, 129, 103292, 2020

Engineering a Distributed Full-Text Index

Authors: Johannes Fischer, Florian Kurpicz, Peter Sanders

Abstract: We present a distributed full-text index for big data applications in a distributed environment. Our index can answer different types of pattern matching queries (existential, counting and enumeration). We perform experiments on inputs up to 100 GiB using up to 512 processors, and compare our index with the distributed suffix array by Arroyuelo et al. [Parall. Comput. 40(9): 471--495, 2014]. The r… ▽ More We present a distributed full-text index for big data applications in a distributed environment. Our index can answer different types of pattern matching queries (existential, counting and enumeration). We perform experiments on inputs up to 100 GiB using up to 512 processors, and compare our index with the distributed suffix array by Arroyuelo et al. [Parall. Comput. 40(9): 471--495, 2014]. The result is that our index answers counting queries up to 5.5 times faster than the distributed suffix array, while using about the same space. We also provide a succinct variant of our index that uses only one third of the memory compared with our non-succinct variant, at the expense of only 20% slower query times. △ Less

Submitted 5 December, 2016; v1 submitted 11 October, 2016; originally announced October 2016.

A mechanical erosion model for two-phase mass flows

Authors: Shiva P. Pudasaini, Jan-Thomas Fischer

Abstract: Erosion, entrainment and deposition are complex and dominant, but yet poorly understood, mechanical processes in geophysical mass flows. Here, we propose a novel, process-based, two-phase, erosion-deposition model capable of adequately describing these complex phenomena commonly observed in landslides, avalanches, debris flows and bedload transport. The model is based on the jump in the momentum f… ▽ More Erosion, entrainment and deposition are complex and dominant, but yet poorly understood, mechanical processes in geophysical mass flows. Here, we propose a novel, process-based, two-phase, erosion-deposition model capable of adequately describing these complex phenomena commonly observed in landslides, avalanches, debris flows and bedload transport. The model is based on the jump in the momentum flux including changes of material and flow properties along the flow-bed interface and enhances an existing general two-phase mass flow model (Pudasaini, 2012). A two-phase variably saturated erodible basal morphology is introduced and allows for the evolution of erosion-deposition-depths, incorporating the inherent physical process including momentum and rheological changes of the flowing mixture. By rigorous derivation, we show that appropriate incorporation of the mass and momentum productions or losses in conservative model formulation is essential for the physically correct and mathematically consistent description of erosion-entrainment-deposition processes. We show that mechanically deposition is the reversed process of erosion. We derive mechanically consistent closures for coefficients emerging in the erosion rate models. We prove that effectively reduced friction in erosion is equivalent to the momentum production. With this, we solve the long standing dilemma of mass mobility, and show that erosion enhances the mass flow mobility. The novel enhanced real two-phase model reveals some major aspects of the mechanics associated with erosion, entrainment and deposition. The model appropriately captures the emergence and propagation of complex frontal surge dynamics associated with the frontal ambient-drag with erosion. △ Less

Submitted 6 October, 2016; originally announced October 2016.

Journal ref: International Journal of Multiphase Flow, 132, 103416, 2020

Effect of adding nanometre-sized heterogeneities on the structural dynamics and the excess wing of a molecular glass former

Authors: S. Gupta, J. K. H. Fischer, P. Lunkenheimer, A. Loidl, E. Novak, N. Jalarvo, M. Ohl

Abstract: We present the relaxation dynamics of glass-forming glycerol mixed with 1.1 nm sized polyhedral oligomeric silsesquioxane (POSS) molecules using dielectric spectroscopy (DS) and two different neutron scattering (NS) techniques. Both, the reorientational dynamics as measured by DS and the density fluctuations detected by NS reveal a broadening of the alpha relaxation when POSS molecules are added.… ▽ More We present the relaxation dynamics of glass-forming glycerol mixed with 1.1 nm sized polyhedral oligomeric silsesquioxane (POSS) molecules using dielectric spectroscopy (DS) and two different neutron scattering (NS) techniques. Both, the reorientational dynamics as measured by DS and the density fluctuations detected by NS reveal a broadening of the alpha relaxation when POSS molecules are added. Moreover, we find a significant slowing down of the alpha-relaxation time. These effects are in accord with the heterogeneity scenario considered for the dynamics of glasses and supercooled liquids. The addition of POSS also affects the excess wing in glycerol arising from a secondary relaxation process, which seems to exhibit a dramatic increase in relative strength compared to the alpha-relaxation. △ Less

Submitted 22 September, 2016; originally announced September 2016.

Comments: 32 pages, 7 figures, accepted for publication in the journal Scientific Reports

Journal ref: Scientific Reports 6, 35034 (2016)

The Mid-Infrared Evolution of the FU Orionis Disk

Authors: Joel D. Green, Olivia C. Jones, Luke D. Keller, Charles A. Poteet, Yao-Lun Yang, William J. Fischer, Neal J. Evans II, Benjamin A. Sargent, Luisa M. Rebull

Abstract: We present new SOFIA-FORCAST observations obtained in Feburary 2016 of the archetypal outbursting low mass young stellar object FU Orionis, and compare the continuum, solid state, and gas properties with mid-IR data obtained at the same wavelengths in 2004 with Spitzer-IRS. In this study, we conduct the first mid-IR spectroscopic comparison of an FUor over a long time period. Over a 12 year period… ▽ More We present new SOFIA-FORCAST observations obtained in Feburary 2016 of the archetypal outbursting low mass young stellar object FU Orionis, and compare the continuum, solid state, and gas properties with mid-IR data obtained at the same wavelengths in 2004 with Spitzer-IRS. In this study, we conduct the first mid-IR spectroscopic comparison of an FUor over a long time period. Over a 12 year period, UBVR monitoring indicates that FU Orionis has continued its steady decrease in overall brightness by ~ 14%. We find that this decrease in luminosity occurs only at wavelengths < 20 microns. In particular, the continuum short ward of the silicate emission complex at 10 microns exhibits a ~ 12% (~ 3 sigma) drop in flux density, but no apparent change in slope; both the Spitzer and SOFIA spectra are consistent with a 7200 K blackbody. Additionally, the detection of water absorption is consistent with the Spitzer spectrum. The silicate emission feature at 10 microns continues to be consistent with unprocessed grains, unchanged over 12 years. We conclude that either the accretion rate in FU Orionis has decreased by ~ 12-14% over this time baseline, or that the inner disk has cooled, but the accretion disk remains in a superheated state outside of the innermost region. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 12 pages, 7 figures

Stochastic homogenization of linear elliptic equations: Higher-order error estimates in weak norms via second-order correctors

Authors: Peter Bella, Benjamin Fehrman, Julian Fischer, Felix Otto

Abstract: We are concerned with the homogenization of second-order linear elliptic equations with random coefficient fields. For symmetric coefficient fields with only short-range correlations, quantified through a logarithmic Sobolev inequality for the ensemble, we prove that when measured in weak spatial norms, the solution to the homogenized equation provides a higher-order approximation of the solution… ▽ More We are concerned with the homogenization of second-order linear elliptic equations with random coefficient fields. For symmetric coefficient fields with only short-range correlations, quantified through a logarithmic Sobolev inequality for the ensemble, we prove that when measured in weak spatial norms, the solution to the homogenized equation provides a higher-order approximation of the solution to the equation with oscillating coefficients. In the case of nonsymmetric coefficient fields, we provide a higher-order approximation (in weak spatial norms) of the solution to the equation with oscillating coefficients in terms of solutions to constant-coefficient equations. In both settings, we also provide optimal error estimates for the two-scale expansion truncated at second order. Our results rely on novel estimates on the second-order homogenization corrector, which we establish via sensitivity estimates for the second-order corrector and a large-scale $L^p$ theory for elliptic equations with random coefficients. Our results also cover the case of elliptic systems. △ Less

Submitted 4 November, 2016; v1 submitted 6 September, 2016; originally announced September 2016.

Comments: 35 pages

About Learning in Recurrent Bistable Gradient Networks

Authors: J. Fischer, S. Lackner

Abstract: Recurrent Bistable Gradient Networks are attractor based neural networks characterized by bistable dynamics of each single neuron. Coupled together using linear interaction determined by the interconnection weights, these networks do not suffer from spurious states or very limited capacity anymore. Vladimir Chinarov and Michael Menzinger, who invented these networks, trained them using Hebb's lear… ▽ More Recurrent Bistable Gradient Networks are attractor based neural networks characterized by bistable dynamics of each single neuron. Coupled together using linear interaction determined by the interconnection weights, these networks do not suffer from spurious states or very limited capacity anymore. Vladimir Chinarov and Michael Menzinger, who invented these networks, trained them using Hebb's learning rule. We show, that this way of computing the weights leads to unwanted behaviour and limitations of the networks capabilities. Furthermore we evince, that using the first order of Hintons Contrastive Divergence algorithm leads to a quite promising recurrent neural network. These findings are tested by learning images of the MNIST database for handwritten numbers. △ Less

Submitted 29 August, 2016; originally announced August 2016.

Comments: 3 pages, 4 figures

Reconstructing Neural Parameters and Synapses of arbitrary interconnected Neurons from their Simulated Spiking Activity

Authors: J. Fischer, P. Manoonpong, S. Lackner

Abstract: To understand the behavior of a neural circuit it is a presupposition that we have a model of the dynamical system describing this circuit. This model is determined by several parameters, including not only the synaptic weights, but also the parameters of each neuron. Existing works mainly concentrate on either the synaptic weights or the neural parameters. In this paper we present an algorithm to… ▽ More To understand the behavior of a neural circuit it is a presupposition that we have a model of the dynamical system describing this circuit. This model is determined by several parameters, including not only the synaptic weights, but also the parameters of each neuron. Existing works mainly concentrate on either the synaptic weights or the neural parameters. In this paper we present an algorithm to reconstruct all parameters including the synaptic weights of a spiking neuron model. The model based on works of Eugene M. Izhikevich (Izhikevich 2007) consists of two differential equations and covers different types of cortical neurons. It combines the dynamical properties of Hodgkin-Huxley-type dynamics with a high computational efficiency. The presented algorithm uses the recordings of the corresponding membrane potentials of the model for the reconstruction and consists of two main components. The first component is a rank based Genetic Algorithm (GA) which is used to find the neural parameters of the model. The second one is a Least Mean Squares approach which computes the synaptic weights of all interconnected neurons by minimizing the squared error between the calculated and the measured membrane potentials for each time step. In preparation for the reconstruction of the neural parameters and of the synaptic weights from real measured membrane potentials, promising results based on simulated data generated with a randomly parametrized Izhikevich model are presented. The reconstruction does not only converge to a global minimum of neural parameters, but also approximates the synaptic weights with high precision. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: 6 pages, 7 figures

The evolution of far-infrared CO emission from protostars

Authors: P. Manoj, J. D. Green, S. T. Megeath, N. J. Evans II, A. M. Stutz, J. J. Tobin, D. M. Watson, W. J. Fischer, E. Furlan, T. Henning

Abstract: We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar… ▽ More We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar properties. We find a strong and tight correlation between far-IR CO luminosity ($L^{\rm fir}_{\rm CO}$) and the bolometric luminosity ($L_{\rm{bol}}$) of the protostars with $L^{\rm fir}_{\rm CO}$ $\propto$ $L_{\rm{bol}}^{0.7}$. We, however, do not find a strong correlation between $L^{\rm fir}_{\rm CO}$ and protostellar evolutionary indicators, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$. FIR CO emission from protostars traces the currently shocked gas by jets/outflows, and $L^{\rm fir}_{\rm CO}$ is proportional to the instantaneous mass loss rate, $\dot{M}_{\rm{out}}$. The correlation between $L^{\rm fir}_{\rm CO}$ and $L_{\rm{bol}}$ is indicative of instantaneous $\dot{M}_{\rm{out}}$ tracking instantaneous $\dot{M}_{\rm{acc}}$. The lack of correlation between $L^{\rm fir}_{\rm CO}$ and evolutionary indicators $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ suggests that $\dot{M}_{\rm{out}}$ and, therefore, $\dot{M}_{\rm{acc}}$ do not show any clear evolutionary trend. These results are consistent with mass accretion/ejection in protostars being episodic. Taken together with the previous finding that the time-averaged mass ejection/accretion rate declines during the protostellar phase (e.g., Bontemps et al. 1996), our results suggest that the instantaneous accretion/ejection rate of protostars is highly time variable and episodic, but the amplitude and/or frequency of this variability decreases with time such that the time averaged accretion/ejection rate declines with system age. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: Accepted for publication in ApJ

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

On the Benefit of Merging Suffix Array Intervals for Parallel Pattern Matching

Authors: Johannes Fischer, Dominik Köppl, Florian Kurpicz

Abstract: We present parallel algorithms for exact and approximate pattern matching with suffix arrays, using a CREW-PRAM with $p$ processors. Given a static text of length $n$, we first show how to compute the suffix array interval of a given pattern of length $m$ in $O(\frac{m}{p}+ \lg p + \lg\lg p\cdot\lg\lg n)$ time for $p \le m$. For approximate pattern matching with $k$ differences or mismatches, we s… ▽ More We present parallel algorithms for exact and approximate pattern matching with suffix arrays, using a CREW-PRAM with $p$ processors. Given a static text of length $n$, we first show how to compute the suffix array interval of a given pattern of length $m$ in $O(\frac{m}{p}+ \lg p + \lg\lg p\cdot\lg\lg n)$ time for $p \le m$. For approximate pattern matching with $k$ differences or mismatches, we show how to compute all occurrences of a given pattern in $O(\frac{m^kσ^k}{p}\max\left(k,\lg\lg n\right)\!+\!(1+\frac{m}{p}) \lg p\cdot \lg\lg n + \text{occ})$ time, where $σ$ is the size of the alphabet and $p \le σ^k m^k$. The workhorse of our algorithms is a data structure for merging suffix array intervals quickly: Given the suffix array intervals for two patterns $P$ and $P'$, we present a data structure for computing the interval of $PP'$ in $O(\lg\lg n)$ sequential time, or in $O(1+\lg_p\lg n)$ parallel time. All our data structures are of size $O(n)$ bits (in addition to the suffix array). △ Less

Submitted 8 June, 2016; originally announced June 2016.

A WISE Census of Young Stellar Objects in Canis Major

Authors: William J. Fischer, Deborah L. Padgett, Karl L. Stapelfeldt, Marta Sewilo

Abstract: With the Wide-Field Infrared Survey Explorer (WISE), we searched for young stellar objects (YSOs) in a 100 deg^2 region centered on the lightly studied Canis Major star forming region. Applying stringent magnitude cuts to exclude the majority of extragalactic contaminants, we find 144 Class I candidates and 335 Class II candidates. The sensitivity to Class II candidates is limited by their faintne… ▽ More With the Wide-Field Infrared Survey Explorer (WISE), we searched for young stellar objects (YSOs) in a 100 deg^2 region centered on the lightly studied Canis Major star forming region. Applying stringent magnitude cuts to exclude the majority of extragalactic contaminants, we find 144 Class I candidates and 335 Class II candidates. The sensitivity to Class II candidates is limited by their faintness at the distance to Canis Major (assumed as 1000 pc). More than half the candidates (53%) are found in 16 groups of more than four members, including four groups with more than 25 members each. The ratio of Class II to Class I objects, N_II/N_I, varies from 0.4 to 8.3 in just the largest four groups. We compare our results to those obtainable with combined Two Micron All Sky Survey (2MASS) and post-cryogenic Spitzer Space Telescope data; the latter approach recovers missing Class II sources. Via a comparison to protostars characterized with the Herschel Space Observatory, we propose new WISE color criteria for flat-spectrum and Class 0 protostars, finding 80 and seven of these, respectively. The distribution of YSOs in CMa OB1 is consistent with supernova-induced star formation, although the diverse N_II/N_I ratios are unexpected if this parameter traces age and the YSOs are due to the same supernova. Less massive clouds feature larger N_II/N_I ratios, suggesting that initial conditions play a role in determining this quantity. △ Less

Submitted 6 June, 2016; originally announced June 2016.

Comments: Accepted by ApJ, 25 pages, 28 figures

A $Herschel-SPIRE$ Survey of the Mon R2 Giant Molecular Cloud: Analysis of the Gas Column Density Probability Density Function

Authors: R. Pokhrel, R. Gutermuth, B. Ali, T. Megeath, J. Pipher, P. Myers, W. J. Fischer, T. Henning, S. J. Wolk, L. Allen, J. J. Tobin

Abstract: We present a far-IR survey of the entire Mon R2 GMC with $Herschel-SPIRE$ cross-calibrated with $Planck-HFI$ data. We fit the SEDs of each pixel with a greybody function and an optimal beta value of 1.8. We find that mid-range column densities obtained from far-IR dust emission and near-IR extinction are consistent. For the entire GMC, we find that the column density histogram, or N-PDF, is lognor… ▽ More We present a far-IR survey of the entire Mon R2 GMC with $Herschel-SPIRE$ cross-calibrated with $Planck-HFI$ data. We fit the SEDs of each pixel with a greybody function and an optimal beta value of 1.8. We find that mid-range column densities obtained from far-IR dust emission and near-IR extinction are consistent. For the entire GMC, we find that the column density histogram, or N-PDF, is lognormal below $\sim$10$^{21}$ cm$^{-2}$. Above this value, the distribution takes a power law form with an index of -2.16. We analyze the gas geometry, N-PDF shape, and YSO content of a selection of subregions in the cloud. We find no regions with pure lognormal N-PDFs. The regions with a combination of lognormal and one power law N-PDF have a YSO cluster and a corresponding centrally concentrated gas clump. The regions with a combination of lognormal and two power law N-PDF have significant numbers of typically younger YSOs but no prominent YSO cluster. These regions are composed of an aggregate of closely spaced gas filaments with no concentrated dense gas clump. We find that for our fixed scale regions, the YSO count roughly correlates with the N-PDF power law index. The correlation appears steeper for single power law regions relative to two power law regions with a high column density cut-off, as a greater dense gas mass fraction is achieved in the former. A stronger correlation is found between embedded YSO count and the dense gas mass among our regions. △ Less

Submitted 6 June, 2016; originally announced June 2016.

Comments: 15 pages, 15 figures, 3 tables. Accepted in MNRAS

Magnetic properties and Curie temperatures of disordered Heusler compounds: Co(1+x)Fe(2-x)Si

Authors: Julia Erika Fischer, Julie Karel, Simone Fabbrici, Peter Adler, Siham Ouardi, Gerhard H. Fecher, Franca Albertini, Claudia Felser

Abstract: The local atomic environments and magnetic properties were investigated for a series of Co(1+x)Fe(2-x)Si (0<x<1) Heusler compounds. While the total magnetic moment in these compounds increases with the number of valance electrons, the highest Curie temperature (Tc) in this series was found for Co1.5Fe1.5Si, with a Tc of 1069 K (24 K higher than the well known Co2FeSi). 57Fe Mössbauer spectroscopy… ▽ More The local atomic environments and magnetic properties were investigated for a series of Co(1+x)Fe(2-x)Si (0<x<1) Heusler compounds. While the total magnetic moment in these compounds increases with the number of valance electrons, the highest Curie temperature (Tc) in this series was found for Co1.5Fe1.5Si, with a Tc of 1069 K (24 K higher than the well known Co2FeSi). 57Fe Mössbauer spectroscopy was used to characterize the local atomic order and to estimate the Co and Fe magnetic moments. Consideration of the local magnetic moments and the exchange integrals is necessary to understand the trend in Tc. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Journal ref: Phys. Rev. B 94, 024418 (2016)

Host Galaxy Identification for Supernova Surveys

Authors: Ravi R. Gupta, Steve Kuhlmann, Eve Kovacs, Harold Spinka, Richard Kessler, Daniel A. Goldstein, Camille Liotine, Katarzyna Pomian, Chris B. D'Andrea, Mark Sullivan, Jorge Carretero, Francisco J. Castander, Robert C. Nichol, David A. Finley, John A. Fischer, Ryan J. Foley, Alex G. Kim, Andreas Papadopoulos, Masao Sako, Daniel M. Scolnic, Mathew Smith, Brad E. Tucker, Syed Uddin, Rachel C. Wolf, Fang Yuan , et al. (35 additional authors not shown)

Abstract: Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve p… ▽ More Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey. △ Less

Submitted 24 June, 2016; v1 submitted 20 April, 2016; originally announced April 2016.

Comments: 22 pages, 18 figures; Accepted by AJ, revised to incorporate referee comments

Liouville principles and a large-scale regularity theory for random elliptic operators on the half-space

Authors: Julian Fischer, Claudia Raithel

Abstract: We consider the large-scale regularity of solutions to second-order linear elliptic equations with random coefficient fields. In contrast to previous works on regularity theory for random elliptic operators, our interest is in the regularity at the boundary: We consider problems posed on the half-space with homogeneous Dirichlet boundary conditions and derive an associated $C^{1,α}$-type large-sca… ▽ More We consider the large-scale regularity of solutions to second-order linear elliptic equations with random coefficient fields. In contrast to previous works on regularity theory for random elliptic operators, our interest is in the regularity at the boundary: We consider problems posed on the half-space with homogeneous Dirichlet boundary conditions and derive an associated $C^{1,α}$-type large-scale regularity theory in the form of a corresponding decay estimate for the homogenization-adapted tilt-excess. This regularity theory entails an associated Liouville-type theorem. The results are based on the existence of homogenization correctors adapted to the half-space setting, which we construct - by an entirely deterministic argument - as a modification of the homogenization corrector on the whole space. This adaption procedure is carried out inductively on larger scales, crucially relying on the regularity theory already established on smaller scales. △ Less

Submitted 25 October, 2016; v1 submitted 10 April, 2016; originally announced April 2016.

Comments: 34 pages; minor revision

Real-Time Computation of Parameter Fitting and Image Reconstruction Using Graphical Processing Units

Authors: Uldis Locans, Andreas Adelmann, Andreas Suter, Jannis Fischer, Werner Lustermann, Gunther Dissertori, Qiulin Wang

Abstract: In recent years graphical processing units (GPUs) have become a powerful tool in scientific computing. Their potential to speed up highly parallel applications brings the power of high performance computing to a wider range of users. However, programming these devices and integrating their use in existing applications is still a challenging task. In this paper we examined the potential of GPUs f… ▽ More In recent years graphical processing units (GPUs) have become a powerful tool in scientific computing. Their potential to speed up highly parallel applications brings the power of high performance computing to a wider range of users. However, programming these devices and integrating their use in existing applications is still a challenging task. In this paper we examined the potential of GPUs for two different applications. The first application, created at Paul Scherrer Institut (PSI), is used for parameter fitting during data analysis of muSR (muon spin rotation, relaxation and resonance) experiments. The second application, developed at ETH, is used for PET (Positron Emission Tomography) image reconstruction and analysis. Applications currently in use were examined to identify parts of the algorithms in need of optimization. Efficient GPU kernels were created in order to allow applications to use a GPU, to speed up the previously identified parts. Benchmarking tests were performed in order to measure the achieved speedup During this work, we focused on single GPU systems to show that real time data analysis of these problems can be achieved without the need for large computing clusters. The results show that the currently used application for parameter fitting, which uses OpenMP to parallelize calculations over multiple CPU cores, can be accelerated around 40 times through the use of a GPU. The speedup may vary depending on the size and complexity of the problem. For PET image analysis, the obtained speedups of the GPU version was more than x40 larger compared to a single core CPU implementation. The achieved results show that it is possible to improve the execution time by orders of magnitude. △ Less

Submitted 22 November, 2016; v1 submitted 8 April, 2016; originally announced April 2016.

The high mass end of the stellar mass function: Dependence on stellar population models and agreement between fits to the light profile

Authors: M. Bernardi, A. Meert, R. K. Sheth, J. -L. Fischer, M. Huertas-Company, C. Maraston, F. Shankar, V. Vikram

Abstract: We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from di… ▽ More We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared. △ Less

Submitted 27 February, 2017; v1 submitted 4 April, 2016; originally announced April 2016.

Comments: 18 pages, 17 figures, accepted for publication in MNRAS. The PyMorph luminosities and stellar masses are available at https://www.physics.upenn.edu/~ameert/SDSS_PhotDec/

Broad [CII] line wings as tracer of molecular and multi-phase outflows in infrared bright galaxies

Authors: A. W. Janssen, N. Christopher, E. Sturm, S. Veilleux, A. Contursi, E. González-Alfonso, J. Fischer, R. Davies, A. Verma, J. Graciá-Carpio, R. Genzel, D. Lutz, A. Sternberg, L. Tacconi, L. Burtscher, A. Poglitsch

Abstract: We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,μ\text{m}$ and [CII] emission at $158\,μ\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample… ▽ More We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,μ\text{m}$ and [CII] emission at $158\,μ\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample have a broad wing component as traced by [CII], and all of these also show OH119 absorption indicative of an outflow (in 1 case an inflow). The other eight cases, where no broad [CII] component was found, are predominantly objects with no OH outflow or a low-velocity ($\leq 100\,\text{km s}^{-1}$) OH outflow. The full width at half maximum (FWHM) of the broad [CII] component shows a trend with the OH119 blue-shifted velocity, although with significant scatter. Moreover, and despite large uncertainties, the outflow masses derived from OH and broad [CII] show a 1:1 relation. The main conclusion is therefore that broad [CII] wings can be used to trace molecular outflows. This may be particularly relevant at high redshift, where the usual tracers of molecular gas (like low-J CO lines) become hard to observe. Additionally, observations of blue-shifted Na I D $λλ5890,5896$ absorption are available for ten of our sources. Outflow velocities of Na I D show a trend with OH velocity and broad [CII] FWHM. These observations suggest that the atomic and molecular gas phases of the outflow are connected. △ Less

Submitted 1 April, 2016; originally announced April 2016.

Comments: accepted for publication in ApJ

Operator Norm Inequalities between Tensor Unfoldings on the Partition Lattice

Authors: Miaoyan Wang, Khanh Dao Duc, Jonathan Fischer, Yun S. Song

Abstract: Interest in higher-order tensors has recently surged in data-intensive fields, with a wide range of applications including image processing, blind source separation, community detection, and feature extraction. A common paradigm in tensor-related algorithms advocates unfolding (or flattening) the tensor into a matrix and applying classical methods developed for matrices. Despite the popularity of… ▽ More Interest in higher-order tensors has recently surged in data-intensive fields, with a wide range of applications including image processing, blind source separation, community detection, and feature extraction. A common paradigm in tensor-related algorithms advocates unfolding (or flattening) the tensor into a matrix and applying classical methods developed for matrices. Despite the popularity of such techniques, how the functional properties of a tensor changes upon unfolding is currently not well understood. In contrast to the body of existing work which has focused almost exclusively on matricizations, we here consider all possible unfoldings of an order-$k$ tensor, which are in one-to-one correspondence with the set of partitions of $\{1,\ldots,k\}$. We derive general inequalities between the $l^p$-norms of arbitrary unfoldings defined on the partition lattice. In particular, we demonstrate how the spectral norm ($p=2$) of a tensor is bounded by that of its unfoldings, and obtain an improved upper bound on the ratio of the Frobenius norm to the spectral norm of an arbitrary tensor. For specially-structured tensors satisfying a generalized definition of orthogonal decomposability, we prove that the spectral norm remains invariant under specific subsets of unfolding operations. △ Less

Submitted 15 January, 2017; v1 submitted 17 March, 2016; originally announced March 2016.

Comments: 17 pages, 1 figure

MSC Class: 15A60; 15A69; 06B99; 05A18

Journal ref: Linear Algebra and its Applications, Vol. 520 (2017) 44-66

An HST Survey for 100-1000 AU Companions around Young Stellar Objects in the Orion Molecular Clouds: Evidence for Environmentally Dependent Multiplicity

Authors: M. Kounkel, S. T. Megeath, C. A. Poteet, W. J. Fischer, L. Hartmann

Abstract: We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$μ$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$μ$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-mai… ▽ More We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$μ$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$μ$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-main sequence stars with disks observed with WFC3, we detect 21 and 28 candidate companions between the projected separations of 100---1000 AU, of which less than 5 and 8, respectively, are chance line of sight coincidences. The resulting companion fraction ($CF$) after the correction for the line of sight contamination is 14.4$^{+1.1}_{-1.3}$% for protostars and 12.5$^{+1.2}_{-0.8}$% for the pre-main sequence stars. These values are similar to those found for main sequence stars, suggesting that there is little variation in the $CF$ with evolution, although several observational biases may mask a decrease in the $CF$ from protostars to the main sequence stars. After segregating the sample into two populations based on the surrounding surface density of YSOs, we find that the $CF$ in the high stellar density regions ($Σ_{YSO} > 45$~pc$^{-2}$) is approximately 50% higher than that found in the low stellar density regions ($Σ_{YSO} < 45$~pc$^{-2}$). We interpret this as evidence for the elevated formation of companions at 100 to 1000 AU in the denser environments of Orion. We discuss possible reasons for this elevated formation. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 20 pages, 7 tables, 22 figures Accepted to ApJ

The Herschel Orion Protostar Survey: Spectral Energy Distributions and Fits Using a Grid of Protostellar Models

Authors: E. Furlan, W. J. Fischer, B. Ali, A. M. Stutz, T. Stanke, J. J. Tobin, S. T. Megeath, M. Osorio, L. Hartmann, N. Calvet, C. A. Poteet, J. Booker, P. Manoj, D. M. Watson, L. Allen

Abstract: We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel,… ▽ More We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $μ$m and sample the peak of the protostellar envelope emission at ~100 $μ$m. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well-fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs. △ Less

Submitted 30 March, 2016; v1 submitted 23 February, 2016; originally announced February 2016.

Comments: Accepted by ApJS; 87 pages; 2 online-only figure sets; corrected typos and other minor fixes

A Dark Energy Camera Search for Missing Supergiants in the LMC After the Advanced LIGO Gravitational Wave Event GW150914

Authors: J. Annis, M. Soares-Santos, E. Berger, D. Brout, H. Chen, R. Chornock, P. S. Cowperthwaite, H. T. Diehl, Z. Doctor, A. Drlica-Wagner, M. R. Drout, B. Farr, D. A. Finley, B. Flaugher, R. J. Foley, J. Frieman, R. A. Gruendl, K. Herner, D. Holz, R. Kessler, H. Lin, J. Marriner, E. Neilsen, A. Rest, M. Sako , et al. (86 additional authors not shown)

Abstract: The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a bur… ▽ More The collapse of the core of a star is expected to produce gravitational radiation. While this process will usually produce a luminous supernova, the optical signatue could be subluminous and a direct collapse to a black hole, with the star just disappearing, is possible. The gravitational wave event GW150914 reported by the LIGO Virgo Collaboration (LVC) on 2015 September 16, was detected by a burst analysis and whose high probability spatial localization included the Large Magellanic Cloud. Shortly after the announcement of the event, we used the Dark Energy Camera to observe 102 deg$^2$ of the localization area, including a 38 deg$^2$ area centered on the LMC. Using a catalog of 152 LMC luminous red supergiants, candidates to undergo a core collapse without a visible supernova, we find that the positions of 144 of these are inside our images, and that all are detected - none have disappeared. There are other classes of candidates: we searched existing catalogs of red supergiants, yellow supergiants, Wolf-Rayet stars, and luminous blue variable stars, recovering all that were inside the imaging area. Based on our observations, we conclude that it is unlikely that GW150914 was caused by the core collapse of a supergiant in the LMC, consistent with the LIGO Collaboration analyses of the gravitational wave form as best described by a binary black hole merger. We discuss how to generalize this search for future very nearby core collapse candidates. △ Less

Submitted 18 February, 2016; v1 submitted 12 February, 2016; originally announced February 2016.

Comments: 6 pages, 3 figures, as submitted

Report number: FERMILAB-PUB-16-032-AE-PPD

A Dark Energy Camera Search for an Optical Counterpart to the First Advanced LIGO Gravitational Wave Event GW150914

Authors: M. Soares-Santos, R. Kessler, E. Berger, J. Annis, D. Brout, E. Buckley-Geer, H. Chen, P. S. Cowperthwaite, H. T. Diehl, Z. Doctor, A. Drlica-Wagner, B. Farr, D. A. Finley, B. Flaugher, R. J. Foley, J. Frieman, R. A. Gruendl, K. Herner, D. Holz, H. Lin, J. Marriner, E. Neilsen, A. Rest, M. Sako, D. Scolnic , et al. (94 additional authors not shown)

Abstract: We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z ban… ▽ More We report initial results of a deep search for an optical counterpart to the gravitational wave event GW150914, the first trigger from the Advanced LIGO gravitational wave detectors. We used the Dark Energy Camera (DECam) to image a 102 deg$^2$ area, corresponding to 38% of the initial trigger high-probability sky region and to 11% of the revised high-probability region. We observed in i and z bands at 4-5, 7, and 24 days after the trigger. The median $5σ$ point-source limiting magnitudes of our search images are i=22.5 and z=21.8 mag. We processed the images through a difference-imaging pipeline using templates from pre-existing Dark Energy Survey data and publicly available DECam data. Due to missing template observations and other losses, our effective search area subtends 40 deg$^{2}$, corresponding to 12% total probability in the initial map and 3% of the final map. In this area, we search for objects that decline significantly between days 4-5 and day 7, and are undetectable by day 24, finding none to typical magnitude limits of i= 21.5,21.1,20.1 for object colors (i-z)=1,0,-1, respectively. Our search demonstrates the feasibility of a dedicated search program with DECam and bodes well for future research in this emerging field. △ Less

Submitted 18 February, 2016; v1 submitted 12 February, 2016; originally announced February 2016.

Comments: 6 pages, 1 figure. Updated references. Submitted to ApJL

Report number: FERMILAB-PUB-16-031-AE-PPD

SDSS-II Supernova Survey: An Analysis of the Largest Sample of Type Ia Supernovae and Correlations with Host-Galaxy Spectral Properties

Authors: Rachel C. Wolf, Chris B. D'Andrea, Ravi R. Gupta, Masao Sako, John A. Fischer, Rick Kessler, Saurabh W. Jha, Marisa C. March, Daniel M. Scolnic, Johanna-Laina Fischer, Heather Campbell, Robert C. Nichol, Matthew D. Olmstead, Michael Richmond, Donald P. Schneider, Mathew Smith

Abstract: Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis u… ▽ More Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis utilizes host-galaxy spectroscopy obtained during the SDSS-I/II spectroscopic survey and from an ancillary program on the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) that obtained spectra for nearly all host galaxies of SDSS-II SN candidates. In addition, we use photometric host-galaxy properties from the SDSS-SNS data release (Sako et al. 2014) such as host stellar mass and star-formation rate. We confirm the well-known relation between HR and host-galaxy mass and find a 3.6σ significance of a non-zero linear slope. We also recover correlations between HR and host-galaxy gas-phase metallicity and specific star-formation rate as they are reported in the literature. With our large dataset, we examine correlations between HR and multiple host-galaxy properties simultaneously and find no evidence of a significant correlation. We also independently analyze our spectroscopically-confirmed and photometrically-classified SNe Ia and comment on the significance of similar combined datasets for future surveys. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: Accepted for publication in ApJ

A Coherent Polariton Laser

Authors: Seonghoon Kim, Bo Zhang, Zhaorong Wang, Julian Fischer, Sebastian Brodbeck, Martin Kamp, Christian Schneider, Sven Höfling, Hui Deng

Abstract: The semiconductor polariton laser promises a new source of coherent light, which, compared to conventional semiconductor photon lasers, has input-energy threshold orders of magnitude lower. However, intensity stability, a defining feature of a coherent state, has remained poor. Intensity noise at many times of the shot-noise of a coherent state has persisted, which has been attributed to multiple… ▽ More The semiconductor polariton laser promises a new source of coherent light, which, compared to conventional semiconductor photon lasers, has input-energy threshold orders of magnitude lower. However, intensity stability, a defining feature of a coherent state, has remained poor. Intensity noise at many times of the shot-noise of a coherent state has persisted, which has been attributed to multiple mechanisms that are difficult to separate in conventional polariton systems. The large intensity noise in turn limited the phase coherence. These limit the capability of the polariton laser as a source of coherence light. Here, we demonstrate a polariton laser with shot-noise limited intensity stability, as expected of a fully coherent state. This is achieved by using an optical cavity with high mode selectivity to enforce single-mode lasing, suppress condensate depletion, and establish gain saturation. The absence of spurious intensity fluctuations moreover enabled measurement of a transition from exponential to Gaussian decay of the phase coherence of the polariton laser. It suggests large self-interaction energies in the polariton condensate, exceeding the laser bandwidth. Such strong interactions are unique to matter-wave laser and important for nonlinear polariton devices. The results will guide future development of polariton lasers and nonlinear polariton devices. △ Less

Submitted 1 February, 2016; v1 submitted 1 February, 2016; originally announced February 2016.

Comments: 23 pages, 6 figures

Journal ref: Phys. Rev. X 6, 011026 (2016)

Evolution of Mass Outflow in Protostars

Authors: Dan M. Watson, Nuria P. Calvet, William J. Fischer, W. J. Forrest, P. Manoj, S. Thomas Megeath, Gary J. Melnick, Joan Najita, David A. Neufeld, Patrick D. Sheehan, Amelia M. Stutz, John J. Tobin

Abstract: We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues… ▽ More We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them. △ Less

Submitted 15 November, 2015; originally announced November 2015.

Comments: Version submitted to ApJ: 36 pages, 3 tables, 8 figures

Large anomalous Hall effect driven by non-vanishing Berry curvature in non-collinear antiferromagnetic Mn3Ge

Authors: Ajaya K. Nayak, Julia Fischer, Yan Sun, Binghai Yan, Julie Karel, Alexander Komarek, Chandra Shekhar, Nitesh Kumar, Walter Schnelle, Juergen Kuebler, Claudia Felser, Stuart S. P. Parkin

Abstract: It is well established that the anomalous Hall effect that a ferromagnet displays scales with its magnetization. Therefore, an antiferromagnet that has no net magnetization should exhibit no anomalous Hall effect. Here we show that the non-collinear triangular antiferromagnet Mn3Ge exhibits a large anomalous Hall effect comparable to that of ferromagnetic metals; the magnitude of the anomalous con… ▽ More It is well established that the anomalous Hall effect that a ferromagnet displays scales with its magnetization. Therefore, an antiferromagnet that has no net magnetization should exhibit no anomalous Hall effect. Here we show that the non-collinear triangular antiferromagnet Mn3Ge exhibits a large anomalous Hall effect comparable to that of ferromagnetic metals; the magnitude of the anomalous conductivity is 500 per ohm per cm at 2 K and 50 per ohm per cm at room temperature. The angular dependence of the anomalous Hall effect measurements confirm that the small residual in-plane magnetic moment has no role in the observed effect. Our theoretical calculations demonstrate that the large anomalous Hall effect in Mn3Ge originates from a non-vanishing Berry curvature that arises from the chiral spin structure, and which also results in a large spin Hall effect, comparable to that of platinum. The present results pave the way to realize room temperature antiferromagnetic spintronics and spin Hall effect based data storage devices. △ Less

Submitted 18 April, 2016; v1 submitted 10 November, 2015; originally announced November 2015.

Comments: Published in Science Advances

Journal ref: Science Advances 2, e1501870 (2016)

The WISE Census of Young Stellar Objects in Canis Major

Authors: William J. Fischer, Deborah L. Padgett, Karl R. Stapelfeldt

Abstract: While searches for young stellar objects (YSOs) with the Spitzer Space Telescope focused on known molecular clouds, photometry from the Wide-field Infrared Survey Explorer (WISE) can be used to extend the search to the entire sky. As a precursor to more expansive searches, we present results for a 100 square degree region centered on the Canis Major clouds. While searches for young stellar objects (YSOs) with the Spitzer Space Telescope focused on known molecular clouds, photometry from the Wide-field Infrared Survey Explorer (WISE) can be used to extend the search to the entire sky. As a precursor to more expansive searches, we present results for a 100 square degree region centered on the Canis Major clouds. △ Less

Submitted 20 October, 2015; originally announced October 2015.

Comments: 2 pages, 1 figure; to appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314 (Cambridge University Press), J.H. Kastner, B. Stelzer, S.A. Metchev, eds