Showing 1–46 of 46 results for author: White, C

SkyMapper Southern Survey: Data Release 4

Authors: Christopher A. Onken, Christian Wolf, Michael S. Bessell, Seo-Won Chang, Lance C. Luvaul, John L. Tonry, Marc C. White, Gary S. Da Costa

Abstract: We present the fourth data release (DR4) of the SkyMapper Southern Survey (SMSS), the last major step in our hemispheric survey with six optical filters: u, v, g, r, i, z. SMSS DR4 covers 26,000 sq.deg from over 400,000 images acquired by the 1.3m SkyMapper telescope between 2014-03 and 2021-09. The 6-band sky coverage extends from the South Celestial Pole to Dec = +16deg, with some images reachin… ▽ More We present the fourth data release (DR4) of the SkyMapper Southern Survey (SMSS), the last major step in our hemispheric survey with six optical filters: u, v, g, r, i, z. SMSS DR4 covers 26,000 sq.deg from over 400,000 images acquired by the 1.3m SkyMapper telescope between 2014-03 and 2021-09. The 6-band sky coverage extends from the South Celestial Pole to Dec = +16deg, with some images reaching Dec ~ +28deg. In contrast to previous DRs, we include all good-quality images from the facility taken during that time span, not only those explicitly taken for the public Survey. From the image dataset, we produce a catalogue of nearly 13 billion detections made from ~700 million unique astrophysical objects. The typical 10sigma depths for each field range between 18.5 and 20.5 mag, depending on the filter, but certain sky regions include longer exposures that reach as deep as 22 mag in some filters. As with previous SMSS catalogues, we have cross-matched with a host of other imaging and spectroscopic datasets to facilitate additional science outcomes. SMSS DR4 is now available to the worldwide astronomical community. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 32 pages. Submitted to the Publications of the Astronomical Society of Australia

TREASUREHUNT: Transients and Variability Discovered with HST in the JWST North Ecliptic Pole Time Domain Field

Authors: Rosalia O'Brien, Rolf A. Jansen, Norman A. Grogin, Seth H. Cohen, Brent M. Smith, Ross M. Silver, W. P. Maksym III, Rogier A. Windhorst, Timothy Carleton, Anton M. Koekemoer, Nimish P. Hathi, Christopher N. A. Willmer, Brenda L. Frye, M. Alpaslan, M. L. N. Ashby, T. A. Ashcraft, S. Bonoli, W. Brisken, N. Cappelluti, F. Civano, C. J. Conselice, V. S. Dhillon, S. P. Driver, K. J. Duncan, R. Dupke , et al. (34 additional authors not shown)

Abstract: The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with… ▽ More The JWST North Ecliptic Pole (NEP) Time Domain Field (TDF) is a $>$14 arcmin diameter field optimized for multi-wavelength time-domain science with JWST. It has been observed across the electromagnetic spectrum both from the ground and from space, including with the Hubble Space Telescope (HST). As part of HST observations over 3 cycles (the "TREASUREHUNT" program), deep images were obtained with ACS/WFC in F435W and F606W that cover almost the entire JWST NEP TDF. Many of the individual pointings of these programs partially overlap, allowing an initial assessment of the potential of this field for time-domain science with HST and JWST. The cumulative area of overlapping pointings is ~88 arcmin$^2$, with time intervals between individual epochs that range between 1 day and 4$+$ years. To a depth of $m_{AB}$ $\simeq$ 29.5 mag (F606W), we present the discovery of 12 transients and 190 variable candidates. For the variable candidates, we demonstrate that Gaussian statistics are applicable, and estimate that ~80 are false positives. The majority of the transients will be supernovae, although at least two are likely quasars. Most variable candidates are AGN, where we find 0.42% of the general $z$ $<$ 6 field galaxy population to vary at the $~3σ$ level. Based on a 5-year timeframe, this translates into a random supernova areal density of up to ~0.07 transients per arcmin$^2$ (~245 deg$^{-2}$) per epoch, and a variable AGN areal density of ~1.25 variables per arcmin$^2$ (~4500 deg$^{-2}$) to these depths. △ Less

Submitted 2 May, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 32 pages, 11 figures, 5 tables, 1 Appendix

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Wind-Fed GRMHD Simulations of Sagittarius A*: Tilt and Alignment of Jets and Accretion Discs, Electron Thermodynamics, and Multi-Scale Modeling of the Rotation Measure

Authors: Sean M. Ressler, Christopher J. White, Eliot Quataert

Abstract: Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D, wind-fed MHD and GRMHD simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including nonzero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulatio… ▽ More Wind-fed models offer a unique way to form predictive models of the accretion flow surrounding Sagittarius A*. We present 3D, wind-fed MHD and GRMHD simulations spanning the entire dynamic range of accretion from parsec scales to the event horizon. We expand on previous work by including nonzero black hole spin and dynamically evolved electron thermodynamics. Initial conditions for these simulations are generated from simulations of the observed Wolf-Rayet stellar winds in the Galactic Centre. The resulting flow tends to be highly magnetized ($β\approx 2$) with an $\sim$ $r^{-1}$ density profile independent of the strength of magnetic fields in the winds. Our simulations reach the MAD state for some, but not all cases. In tilted flows, SANE jets tend to align with the angular momentum of the gas at large scales, even if that direction is perpendicular to the black hole spin axis. Conversely, MAD jets tend to align with the black hole spin axis. The gas angular momentum shows similar behavior: SANE flows tend to only partially align while MAD flows tend to fully align. With a limited number of dynamical free parameters, our models can produce accretion rates, 230 GHz flux, and unresolved linear polarization fractions roughly consistent with observations for several choices of electron heating fraction. Absent another source of large-scale magnetic field, winds with a higher degree of magnetization (e.g., where the magnetic pressure is 1/100 of the ram pressure in the winds) may be required to get a sufficiently large RM with consistent sign. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: Accepted by MNRAS. Animations for several figures in the paper are available at https://www.youtube.com/playlist?list=PL3pLmTeUPcqSd4jVBnRubYQpa-Dma25ir

The Gravitational-Wave Signature of Core-Collapse Supernovae

Authors: David Vartanyan, Adam Burrows, Tianshu Wang, Matthew S. B. Coleman, Christopher J. White

Abstract: We calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova simulations spanning a range of massive stars. Most of the simulations are carried out to times late enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inexorably rises as the… ▽ More We calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova simulations spanning a range of massive stars. Most of the simulations are carried out to times late enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inexorably rises as the proto-neutron star (PNS) core shrinks. We demonstrate that the GW emission is excited mostly by accretion plumes onto the PNS that energize modal oscillations and also high-frequency (``haze") emission correlated with the phase of violent accretion. The duration of the major phase of emission varies with exploding progenitor and there is a strong correlation between the total GW energy radiated and the compactness of the progenitor. Moreover, the total GW emissions vary by as much as three orders of magnitude from star to star. For black-hole formation, the GW signal tapers off slowly and does not manifest the haze seen for the exploding models. For such failed models, we also witness the emergence of a spiral shock motion that modulates the GW emission at a frequency near $\sim$100 Hertz that slowly increases as the stalled shock sinks. We find significant angular anisotropy of both the high- and low-frequency (memory) GW emissions, though the latter have very little power. △ Less

Submitted 30 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: accepted to PRD

An Extension of the Athena++ Code Framework for Radiation-Magnetohydrodynamics in General Relativity Using a Finite-Solid-Angle Discretization

Authors: Christopher J. White, Patrick D. Mullen, Yan-Fei Jiang, Shane W. Davis, James M. Stone, Viktoriya Morozova, Lizhong Zhang

Abstract: We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical… ▽ More We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems. △ Less

Submitted 28 March, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: To be published in ApJS

Report number: LA-UR-23-21167

Millimeter Observational Signatures of Flares in Magnetically Arrested Black Hole Accretion Models

Authors: He Jia, Bart Ripperda, Eliot Quataert, Christopher J. White, Koushik Chatterjee, Alexander Philippov, Matthew Liska

Abstract: In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expe… ▽ More In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating $γ\gtrsim 100$ ($T_e \gtrsim 10^{11}$ K) electrons. For models in which the electrons are heated to $T_e \sim 10^{11}$ K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multi-wavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization. △ Less

Submitted 23 September, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: 11+7 pages, 9+7 figures, 1 table, accepted by MNRAS

Magnetic Flux Transport in Radiatively Inefficient Accretion Flows and the Pathway towards a Magnetically Arrested Disk

Authors: Prasun Dhang, Xue-Ning Bai, Christopher J. White

Abstract: Large-scale magnetic fields play a vital role in determining the angular momentum transport and generating jets/outflows in the accreting systems, yet their origin remains poorly understood. We focus on radiatively inefficient accretion flows (RIAFs) around the black holes (BHs), and conduct three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations using the Athena++ code. We… ▽ More Large-scale magnetic fields play a vital role in determining the angular momentum transport and generating jets/outflows in the accreting systems, yet their origin remains poorly understood. We focus on radiatively inefficient accretion flows (RIAFs) around the black holes (BHs), and conduct three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations using the Athena++ code. We first re-confirm that the MRI dynamo in the RIAF alone does not spontaneously form a magnetically arrested disk (MAD), conducive for strong jet formation. We next investigate the other possibility, where the large-scale magnetic fields are advected inward from external sources (e.g. the companion star in X-ray binaries, magnetized ambient medium in AGNs). Although the actual configuration of the external fields could be complex and uncertain, they are likely to be closed. As a first study, we treat them as closed field loops of different sizes, shapes and field strengths. Unlike earlier studies of flux transport, where magnetic flux is injected in the initial laminar flow, we injected the magnetic field loops in the quasi-stationary turbulent RIAF in inflow equilibrium and followed their evolution. We found that a substantial fraction ($\sim15\%-40\%$) of the flux injected at the large radii reaches the BH with a weak dependence on the loop parameters except when the loops are injected at high latitudes, away from the mid-plane. Relatively high efficiency of flux transport observed in our study hints that a MAD might be formed relatively easily close to the BH, provided that a source of the large-scale field exists at the larger radii. △ Less

Submitted 5 February, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 23 pages, 16 figures, accepted in ApJ

Remarkable correspondence of Sagittarius A* submillimeter variability with a stellar-wind-fed accretion flow model

Authors: Lena Murchikova, Christopher J. White, Sean M. Ressler

Abstract: We compare the 230 GHz near-horizon emission from Sagittarius A* to simulations representing three classes of accretion flows. Using the structure function to capture the variability statistics of the light curve, we find a noticeable discrepancy between the observations and models based on torus-fed accretion disks, whether those disks bring in a small or large amount of net magnetic flux. On the… ▽ More We compare the 230 GHz near-horizon emission from Sagittarius A* to simulations representing three classes of accretion flows. Using the structure function to capture the variability statistics of the light curve, we find a noticeable discrepancy between the observations and models based on torus-fed accretion disks, whether those disks bring in a small or large amount of net magnetic flux. On the other hand, the simulations that are fed more realistically by stellar winds match the observed structure function very well. We describe the differences between models, arguing that feeding by stellar winds may be a critical component in constructing theoretical models for accretion in the Galactic Center. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 13 pages, 4 figures, 1 appendix

Blacklight: A General-Relativistic Ray-Tracing and Analysis Tool

Authors: Christopher J. White

Abstract: We describe the Blacklight code, intended for post-processing general-relativistic magnetohydrodynamic simulation data. Beyond polarized ray tracing of synchrotron radiation, it can produce a number of outputs that aid in analyzing data sets, such as maps of auxiliary quantities and false-color renderings. Additional features include support for adaptive mesh refinement input, slow-light calculati… ▽ More We describe the Blacklight code, intended for post-processing general-relativistic magnetohydrodynamic simulation data. Beyond polarized ray tracing of synchrotron radiation, it can produce a number of outputs that aid in analyzing data sets, such as maps of auxiliary quantities and false-color renderings. Additional features include support for adaptive mesh refinement input, slow-light calculations, and adaptive ray tracing. The code is written with ease of use, readability, and transparency as primary objectives, while it still achieves high performance. Blacklight is publicly available and released into the public domain. △ Less

Submitted 14 June, 2022; v1 submitted 29 March, 2022; originally announced March 2022.

Comments: Code paper; to be published in ApJS; code available at https://github.com/c-white/blacklight

Observational Signatures of Black Hole Accretion: Rotating vs. Spherical Flows with Tilted Magnetic Fields

Authors: He Jia, Christopher J. White, Eliot Quataert, Sean M. Ressler

Abstract: We study the observational signatures of magnetically arrested black hole accretion with non-rotating inflow onto a rotating black hole; we consider a range of angles between the black hole spin and the initial magnetic field orientation. We compare the results of our General Relativistic Magneto-Hydrodynamic simulations to more commonly used rotating initial conditions and to the Event Horizon Te… ▽ More We study the observational signatures of magnetically arrested black hole accretion with non-rotating inflow onto a rotating black hole; we consider a range of angles between the black hole spin and the initial magnetic field orientation. We compare the results of our General Relativistic Magneto-Hydrodynamic simulations to more commonly used rotating initial conditions and to the Event Horizon Telescope (EHT) observations of M87. We find that the mm intensity images, polarization images, and synchrotron emission spectra are very similar among the different simulations when post-processed with the same electron temperature model; observational differences due to different electron temperature models are significantly larger than those due to the different realizations of magnetically arrested accretion. The orientation of the mm synchrotron polarization is particularly insensitive to the initial magnetic field orientation, the electron temperature model, and the rotation of the inflowing plasma. The largest difference among the simulations with different initial rotation and magnetic tilt is in the strength and stability of the jet; spherical inflow leads to kink-unstable jets. We discuss the implications of our results for current and future EHT observations and for theoretical models of event-horizon-scale black hole accretion. △ Less

Submitted 30 May, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 12 pages, 11 figures, accepted by MNRAS

On the Origin of Pulsar and Magnetar Magnetic Fields

Authors: Christopher J. White, Adam Burrows, Matthew S. B. Coleman, David Vartanyan

Abstract: In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue tha… ▽ More In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective proto-neutron stars at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 solar masses to 25 solar masses, including one with initial rotation, have sufficiently vigorous proto-neutron-star convection to generate dipole fields of order ~10^15 gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars. △ Less

Submitted 19 December, 2021; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJ; very minor changes from previous version

The Physical Effects of Progenitor Rotation: Comparing Two Long-Duration 3D Core-Collapse Supernova Simulations

Authors: Matthew S. B. Coleman, Adam Burrows, Christopher J. White

Abstract: We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despit… ▽ More We analyse and determine the effects of modest progenitor rotation in the context of core-collapse supernovae by comparing two separate long-duration three-dimensional simulations of 9 M$_{\odot}$ progenitors, one rotating with an initial spin period of $\sim$60 seconds and the other non-rotating. We determine that both models explode early, though the rotating model explodes a bit earlier. Despite this difference, the asymptotic explosion energies ($\sim$10$^{50}$ ergs) and residual neutron star baryon masses ($\sim$1.3 M$_{\odot}$) are similar. We find that the proto-neutron star (PNS) core can deleptonize and cool significantly more quickly. Soon into the evolution of the rotating model, we witness more vigorous and extended PNS core convection that early in its evolution envelopes the entire inner sphere, not just a shell. Moreover, we see a corresponding excursion in both the $ν_e$ luminosity and gravitational-wave strain that may be diagnostic of this observed dramatic phenomenon. In addition, after bounce the innermost region of the rotating model seems to execute meridional circulation. The rotationally-induced growth of the convective PNS region may facilitate the growth of core B-fields by the dynamo mechanism by facilitating the achievement of the critical Rossby number condition for substantial growth of a dipole field, obviating the need for rapid rotation rates to create dipole fields of significance. The next step is to explore the progenitor-mass and spin dependencies across the progenitor continuum of the supernova explosion, dynamics, and evolution of PNS convection and its potential role in the generation of magnetar and pulsar magnetic fields. △ Less

Submitted 3 May, 2022; v1 submitted 29 October, 2021; originally announced November 2021.

Comments: Withdrawn pending further calculations

Normal, Dust-Obscured Galaxies in the Epoch of Reionization

Authors: Y. Fudamoto, P. A. Oesch, S. Schouws, M. Stefanon, R. Smit, R. J. Bouwens, R. A. A. Bowler, R. Endsley, V. Gonzalez, H. Inami, I. Labbe, D. Stark, M. Aravena, L. Barrufet, E. da Cunha, P. Dayal, A. Ferrara, L. Graziani, J. Hodge, A. Hutter, Y. Li, I. De Looze, T. Nanayakkara, A. Pallottini, D. Riechers , et al. (4 additional authors not shown)

Abstract: Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extrem… ▽ More Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 16 pages, 7 figures, 4 tables, published in Nature

The Effects of Tilt on the Time Variability of Millimeter and Infrared Emission from Sagittarius A*

Authors: Christopher J. White, Eliot Quataert

Abstract: Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24 degrees) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A… ▽ More Using a combination of general-relativistic magnetohydrodynamics simulations and ray tracing of synchrotron emission, we study the effect of modest (24 degrees) misalignment between the black hole spin and plasma angular momentum, focusing on the variability of total flux, image centroids, and image sizes. We consider both millimeter and infrared (IR) observables motivated by Sagittarius A* (Sgr A*), though our results apply more generally to optically thin flows. For most quantities, tilted accretion is more variable, primarily due to a significantly hotter and denser "coronal" region well off the disk midplane. We find (1) a 150% increase in millimeter light curve variability when adding tilt to the flow; (2) the tilted image centroid in the millimeter shifts on a scale of 3.7 microarcseconds over 28 hours (5000 gravitational times) for some electron temperature models; (3) tilted disk image diameters in the millimeter can be 10% larger (52 versus 47 microarcseconds) than those of aligned disks at certain viewing angles; (4) the tilted models produce significant IR flux, similar to that seen in Sgr A*, with comparable or even greater variability than observed; (5) for some electron models, the tilted IR centroid moves by more than 50 microarcseconds over several hours, in a similar fashion to the centroid motion detected by the GRAVITY interferometer. △ Less

Submitted 17 December, 2021; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJ; figures and numbers slightly updated from previous version with additional data, but conclusions are unchanged

Numerical simulation of hot accretion flows (IV): effects of black hole spin and magnetic field strength on the wind and the comparison between wind and jet properties

Authors: Hai Yang, Feng Yuan, Ye-fei Yuan, Christopher J. White

Abstract: This is the fourth paper of our series of works studying winds from hot accretion flows around black holes. In the first two papers, we have shown the existence of strong winds in hot accretion flows using hydrodynamical and magnetohydrodynamical (MHD) simulations. In the third paper, by using three dimensional general relativity MHD numerical simulation data of hot accretion flows and adopting a… ▽ More This is the fourth paper of our series of works studying winds from hot accretion flows around black holes. In the first two papers, we have shown the existence of strong winds in hot accretion flows using hydrodynamical and magnetohydrodynamical (MHD) simulations. In the third paper, by using three dimensional general relativity MHD numerical simulation data of hot accretion flows and adopting a "virtual particle trajectory" data analysis approach, we have calculated the properties of wind, such as its mass flux and velocity. However, that paper focuses only on a non-spinning black hole and SANE (standard and normal accretion). In the present paper, we extend the third paper by including cases of a rapidly rotating black hole and MAD (magnetically arrested disk). We focus on investigating the effect of spin and magnetic field on the properties of wind and jet. It is found that a larger spin and stronger magnetic field usually enhance the wind and jet. The formulae describing the mass flux, poloidal velocity, and fluxes of momentum, kinetic energy, and total energy of wind and jet are presented. One interesting finding, among others, is that even in the case of very rapidly spinning black hole where the jet is supposed to be the strongest, the momentum flux of jet is smaller than that of wind, while the total energy flux of jet is larger than that of wind by at most a factor of 10. This result suggests that wind potentially plays a more important role than jet at least for some problems in active galactic nuclei feedback. △ Less

Submitted 7 May, 2021; v1 submitted 5 February, 2021; originally announced February 2021.

Comments: 21 pages, 16 figures, accepted for publication in ApJ

Magnetically Modified Spherical Accretion in GRMHD: Reconnection-Driven Convection and Jet Propagation

Authors: Sean M. Ressler, Eliot Quataert, Christopher J. White, Omer Blaes

Abstract: We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulenc… ▽ More We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: Accepted by MNRAS. For animations see [https://smressle.bitbucket.io/animations.html]

The Effects of Resolution on Black Hole Accretion Simulations of Jets

Authors: Christopher J. White, Fiona Chrystal

Abstract: We perform three general-relativistic magnetohydrodynamic simulations of black hole accretion designed to test how sensitive results are to grid resolution in the jet region. The cases differ only in numerics, modelling the same physical scenario of a radiatively inefficient, geometrically thick, magnetically arrested flow onto a rapidly spinning black hole. Properties inferred with the coarsest g… ▽ More We perform three general-relativistic magnetohydrodynamic simulations of black hole accretion designed to test how sensitive results are to grid resolution in the jet region. The cases differ only in numerics, modelling the same physical scenario of a radiatively inefficient, geometrically thick, magnetically arrested flow onto a rapidly spinning black hole. Properties inferred with the coarsest grid generally agree with those found with higher resolutions, including total jet power and its decomposition into different forms, velocity structure, nonaxisymmetric structure, and the appearance of resolved millimetre images. Some measures of variability and magnetization are sensitive to resolution. We conclude that most results obtained by limiting resolution near the jet for computational expediency should still be reliable, at least insofar as they would not be improved with a finer grid. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: Accepted by MNRAS

Search For Electron-Antineutrinos Associated With Gravitational-Wave Events GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817 at Daya Bay

Authors: F. P. An, A. B. Balantekin, H. R. Band, M. Bishai, S. Blyth, G. F. Cao, J. Cao, J. F. Chang, Y. Chang, H. S. Chen, S. M. Chen, Y. Chen, Y. X. Chen, J. Cheng, Z. K. Cheng, J. J. Cherwinka, M. C. Chu, J. P. Cummings, O. Dalager, F. S. Deng, Y. Y. Ding, M. V. Diwan, T. Dohnal, J. Dove, M. Dvorak , et al. (161 additional authors not shown)

Abstract: Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW1… ▽ More Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows. △ Less

Submitted 14 September, 2020; v1 submitted 27 June, 2020; originally announced June 2020.

Comments: 16 pages, 12 figures, 9 tables

Looking for the underlying cause of black hole X-ray variability in GRMHD simulations

Authors: D. A. Bollimpalli, R. Mahmoud, C. Done, P. C. Fragile, W. Kluźniak, R. Narayan, C. J. White

Abstract: Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrin… ▽ More Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags. △ Less

Submitted 27 July, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: 22 pages, 20 figures, 4 tables

Journal ref: MNRAS 496 (2020), 3808-3828

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

Authors: Sean M. Ressler, Christopher J. White, Eliot Quataert, James M. Stone

Abstract: We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained pro… ▽ More We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation. △ Less

Submitted 12 June, 2020; v1 submitted 29 May, 2020; originally announced June 2020.

Comments: Accepted by ApJL. 14 pages, 8 figures. Animations available at https://smressle.bitbucket.io/animations.html

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

Authors: James M. Stone, Kengo Tomida, Christopher J. White, Kyle G. Felker

Abstract: The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic exe… ▽ More The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: 50 pages, 41 figures, accepted for publication in American Astronomical Society journals

The Lyman Continuum Escape Fraction of Galaxies and AGN in the GOODS Fields

Authors: Brent M. Smith, Rogier A. Windhorst, Seth H. Cohen, Anton M. Koekemoer, Rolf A. Jansen, Cameron White, Sanchayeeta Borthakur, Nimish Hathi, Linhua Jiang, Michael Rutkowski, Russell E. Ryan Jr., Akio K. Inoue, Robert W. O'Connell, John W. MacKenty, Christopher J. Conselice, Joseph I. Silk

Abstract: We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AG… ▽ More We present our analysis of the LyC emission and escape fraction of 111 spectroscopically verified galaxies with and without AGN from $2.26<z<4.3$. We extended our ERS sample from Smith et al. (2018; arXiv:1602.01555) with 64 galaxies in the GOODS North and South fields using WFC3/UVIS F225W, F275W, and F336W mosaics we independently drizzled using the HDUV, CANDELS, and UVUDF data. Among the 17 AGN from the 111 galaxies, one provided a LyC detection in F275W at $m_{AB}=23.19$ mag (S/N $\simeq$ 133) and $GALEX$ NUV at $m_{AB}=23.77$ mag (S/N $\simeq$ 13). We simultaneously fit $SDSS$ and $Chandra$ spectra of this AGN to an accretion disk and Comptonization model and find $f_{esc}$ values of $f_{esc}^{F275W}\simeq 28^{+20}_{-4}$% and $f_{esc}^{NUV}\simeq 30^{+22}_{-5}$%. For the remaining 110 galaxies, we stack image cutouts that capture their LyC emission using the F225W, F275W, and F336W data of the GOODS and ERS samples, and both combined, as well as subsamples of galaxies with and without AGN, and $all$ galaxies. We find the stack of 17 AGN dominate the LyC production from $\langle z\rangle\simeq 2.3-4.3$ by a factor of $\sim$10 compared to all 94 galaxies without AGN. While the IGM of the early universe may have been reionized mostly by massive stars, there is evidence that a significant portion of the ionizing energy came from AGN. △ Less

Submitted 10 April, 2020; v1 submitted 9 April, 2020; originally announced April 2020.

Comments: 36 pages, 17 figures, 5 tables. Accepted for publication in The Astrophysical Journal. v1 abstract latex errors corrected, minor changes to table 5, orcid ID corrected for one author

The Structure of Radiatively Inefficient Black Hole Accretion Flows

Authors: Christopher J. White, Eliot Quataert, Charles F. Gammie

Abstract: We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/… ▽ More We run three long-timescale general-relativistic magnetohydrodynamic simulations of radiatively inefficient accretion flows onto non-rotating black holes. Our aim is to achieve steady-state behavior out to large radii and understand the resulting flow structure. A simulation with adiabatic index Gamma = 4/3 and small initial alternating poloidal magnetic field loops is run to a time of 440,000 GM/c^3, reaching inflow equilibrium inside a radius of 370 GM/c^2. Variations with larger alternating field loops and with Gamma = 5/3 are run to 220,000 GM/c^3, attaining equilibrium out to 170 GM/c^2 and 440 GM/c^2. There is no universal self-similar behavior obtained at radii in inflow equilibrium: the Gamma = 5/3 simulation shows a radial density profile with power law index ranging from -1 in the inner regions to -1/2 in the outer regions, while the others have a power-law slope ranging from -1/2 to close to -2. Both simulations with small field loops reach a state with polar inflow of matter, while the more ordered initial field has polar outflows. However, unbound outflows remove only a factor of order unity of the inflowing material over a factor of ~300 in radius. Our results suggest that the dynamics of radiatively inefficient accretion flows are sensitive to how the flow is fed from larger radii, and may differ appreciably in different astrophysical systems. Millimeter images appropriate for Sgr A* are qualitatively (but not quantitatively) similar in all simulations, with a prominent asymmetric image due to Doppler boosting. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 22 pages, 21 figures, accepted by ApJ

The Effects of Tilt on the Images of Black Hole Accretion Flows

Authors: Christopher J. White, Jason Dexter, Omer Blaes, Eliot Quataert

Abstract: We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed… ▽ More We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets. △ Less

Submitted 30 March, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 15 pages, 15 figures, 1 table; accepted by ApJ

Development and Application of Numerical Techniques for General-Relativistic Magnetohydrodynamics Simulations of Black Hole Accretion

Authors: Christopher J. White

Abstract: We describe the implementation of sophisticated numerical techniques for general-relativistic magnetohydrodynamics simulations in the Athena++ code framework. Improvements over many existing codes include the use of advanced Riemann solvers and of staggered-mesh constrained transport. Combined with considerations for computational performance and parallel scalability, these allow us to investigate… ▽ More We describe the implementation of sophisticated numerical techniques for general-relativistic magnetohydrodynamics simulations in the Athena++ code framework. Improvements over many existing codes include the use of advanced Riemann solvers and of staggered-mesh constrained transport. Combined with considerations for computational performance and parallel scalability, these allow us to investigate black hole accretion flows with unprecedented accuracy. The capability of the code is demonstrated by exploring magnetically arrested disks. △ Less

Submitted 23 June, 2019; originally announced June 2019.

Comments: PhD thesis from 2016. 99 pages, 46 figures. A more succinct but less complete presentation of methods can be found in arXiv:1511.00943. More recent and improved application results can be found in arXiv:1903.01509. It is hoped this document can be a pedagogical resource for anyone wishing to understand the details involved in GRMHD codes. Caveat lector: there may well be typos in some formulas

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Authors: Evan B. Bauer, Christopher J. White, Lars Bildsten

Abstract: Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motiv… ▽ More Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motivated by the recent discovery of fast-moving objects that occupy unusual locations on the HR diagram, we explore the impact of the thermonuclear supernovae on the donors in this specific double detonation scenario. We use MESA to model the binary up to the moment of detonation, then 3D Athena++ to model the hydrodynamic interaction of the supernova ejecta with the donor star, calculating the amount of mass that is stripped and the entropy deposited in the deep stellar interior by the strong shock that traverses it. We show that these donor remnants are ejected with velocities primarily set by their orbital speeds: $700-900\ \rm km\ s^{-1}$. We model the long-term thermal evolution of remnants by introducing the shock entropy into MESA models. In response to this entropy change, donor remnants expand and brighten for timescales ranging from $10^6-10^8$ years, giving ample time for these runaway stars to be observed in their inflated state before they leave the galaxy. Even after surface layers are stripped, some donors retain enough mass to resume core helium burning and further delay fading for more than $10^8$ years. △ Less

Submitted 6 October, 2019; v1 submitted 21 June, 2019; originally announced June 2019.

Comments: 12 pages, 8 figures, 3 tables, accepted for publication in ApJ. Movies for figures 1 and 2 available in source

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

A Resolution Study of Magnetically Arrested Disks

Authors: Christopher J. White, James M. Stone, Eliot Quataert

Abstract: We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree ac… ▽ More We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford-Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: 19 pages, 18 figures, to be published in ApJ

Journal ref: ApJ 874 (2019) 168

Tilted Disks around Black Holes: A Numerical Parameter Survey for Spin and Inclination Angle

Authors: Christopher J. White, Eliot Quataert, Omer Blaes

Abstract: We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quick… ▽ More We conduct a systematic study of the properties of tilted accretion flows around spinning black holes, covering a range of tilt angles and black hole spins, using the general-relativistic magnetohydrodynamics code Athena++. The same initial magnetized torus is evolved around black holes with spins ranging from 0 to 0.9, with inclinations ranging from 0 degrees to 24 degrees. The tilted disks quickly reach a warped and twisted shape that rigidly precesses about the black hole spin axis with deformations in shape large enough to hinder the application of linear bending wave theory. Magnetized polar outflows form, oriented along the disk rotation axes. At sufficiently high inclinations a pair of standing shocks develops in the disks. These shocks dramatically affect the flow at small radii, driving angular momentum transport. At high spins they redirect material more effectively than they heat it, reducing the dissipation rate relative to the mass accretion rate and lowering the radiative efficiency of the flow. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 20 pages, 14 figures, to be published in ApJ

Journal ref: ApJ 878 (2019) 51

SkyMapper Southern Survey: First Data Release (DR1)

Authors: Christian Wolf, Christopher A. Onken, Lance C. Luvaul, Brian P. Schmidt, Michael S. Bessell, Seo-Won Chang, Gary S. Da Costa, Dougal Mackey, Tony Martin-Jones, Simon J. Murphy, Tim Preston, Richard A. Scalzo, Li Shao, Jon Smillie, Patrick Tisserand, Marc C. White, Fang Yuan

Abstract: We present the first data release (DR1) of the SkyMapper Southern Survey, a hemispheric survey carried out with the SkyMapper Telescope at Siding Spring Observatory in Australia. Here, we present the survey strategy, data processing, catalogue construction and database schema. The DR1 dataset includes over 66,000 images from the Shallow Survey component, covering an area of 17,200 deg$^2$ in all s… ▽ More We present the first data release (DR1) of the SkyMapper Southern Survey, a hemispheric survey carried out with the SkyMapper Telescope at Siding Spring Observatory in Australia. Here, we present the survey strategy, data processing, catalogue construction and database schema. The DR1 dataset includes over 66,000 images from the Shallow Survey component, covering an area of 17,200 deg$^2$ in all six SkyMapper passbands $uvgriz$, while the full area covered by any passband exceeds 20,000 deg$^2$. The catalogues contain over 285 million unique astrophysical objects, complete to roughly 18 mag in all bands. We compare our $griz$ point-source photometry with PanSTARRS1 DR1 and note an RMS scatter of 2%. The internal reproducibility of SkyMapper photometry is on the order of 1%. Astrometric precision is better than 0.2 arcsec based on comparison with Gaia DR1. We describe the end-user database, through which data are presented to the world community, and provide some illustrative science queries. △ Less

Submitted 23 January, 2018; originally announced January 2018.

Comments: 31 pages, 19 figures, 10 tables, PASA, accepted

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields: Dependence on Jet Radius

Authors: Ken-Ichi Nishikawa, Yosuke Mizuno, Jose L. Gomez, Ioana Dutan, Athina Meli, Charley White, Jacek Niemiec, Oleh Kobzar, Martin Pohl, Asaf Pe'er, Jacob Trier Frederiksen, AAke Nordlund, Helene Sol, Philip E. Hardee, Dieter H. Hartmann

Abstract: In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with… ▽ More In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Comments: 7 pages, 4 figures, submitted to special issue "Polarised Emission from Astrophysical Jets", June 12-16, 2017, Ierapetra, Greece. arXiv admin note: substantial text overlap with arXiv:1609.09363

Journal ref: Galaxies, 5, 58, 2017

Stability Aspects of Wormholes in $R^2$ Gravity

Authors: James B. Dent, Damien A. Easson, Thomas W. Kephart, Sara C. White

Abstract: We study radial perturbations of a wormhole in $R^2$ gravity to determine regions of stability. We also investigate massive and massless particle orbits and tidal forces in this space-time for a radially infalling observer. We study radial perturbations of a wormhole in $R^2$ gravity to determine regions of stability. We also investigate massive and massless particle orbits and tidal forces in this space-time for a radially infalling observer. △ Less

Submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 9 figures

EGG: hatching a mock Universe from empirical prescriptions

Authors: C. Schreiber, D. Elbaz, M. Pannella, E. Merlin, M. Castellano, A. Fontana, N. Bourne, K. Boutsia, F. Cullen, J. Dunlop, H. C. Ferguson, M. J. Michalowski, K. Okumura, P. Santini, X. W. Shu, T. Wang, C. White

Abstract: This paper introduces EGG, the Empirical Galaxy Generator, a tool designed within the ASTRODEEP collaboration to generate mock galaxy catalogs for deep fields with realistic fluxes and simple morphologies. The simulation procedure is based exclusively on empirical prescriptions -- rather than first principles -- to provide the most accurate match with observations at 0<z<7. In particular, we consi… ▽ More This paper introduces EGG, the Empirical Galaxy Generator, a tool designed within the ASTRODEEP collaboration to generate mock galaxy catalogs for deep fields with realistic fluxes and simple morphologies. The simulation procedure is based exclusively on empirical prescriptions -- rather than first principles -- to provide the most accurate match with observations at 0<z<7. In particular, we consider that galaxies can be either quiescent or star-forming, and use their stellar mass (M*) and redshift (z) as the fundamental properties from which all the other observables can be statistically derived. Drawing z and M* from the observed galaxy stellar mass functions, we associate a star formation rate (SFR) to each galaxy from the tight SFR-M* main sequence, while dust attenuation, optical colors and morphologies (including bulge-to-total ratios, sizes and aspect ratios) are obtained from empirical relations that we establish from the high quality Hubble and Herschel observations available in the CANDELS fields. Random scatter is introduced in each step to reproduce the observed distributions of each parameter. Based on these observables, a panchromatic spectral energy distribution (SED) is selected for each galaxy and synthetic photometry is produced by integrating the redshifted SED in common broad-band filters. Finally, the mock galaxies are placed on the sky at random positions with a fixed angular two-point correlation function to implement basic clustering. The resulting flux catalogs reproduce accurately the observed number counts in all broad bands from the ultraviolet up to the sub-millimeter, and can be directly fed to image simulators such as Skymaker. The images can then be used to test source extraction softwares and image-based techniques such as stacking. EGG is open-source, and is made available to the community together with a set of pre-generated catalogs and images. △ Less

Submitted 10 January, 2017; v1 submitted 16 June, 2016; originally announced June 2016.

Comments: 24 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 602, A96 (2017)

An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport

Authors: Christopher J. White, James M. Stone, Charles F. Gammie

Abstract: We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the diverge… ▽ More We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the divergence-free constraint of the magnetic field. Our code is designed to work with arbitrary stationary spacetimes in one, two, or three dimensions, and we demonstrate its reliability in a number of tests. We also report on its promising performance and scalability. △ Less

Submitted 3 July, 2016; v1 submitted 3 November, 2015; originally announced November 2015.

Comments: 30 pages, 17 figures, 1 table; this version accepted to ApJS; changes from previous version include additional tests and more thorough literature comparisons

Turbulent mixing layers in supersonic protostellar outflows, with application to DG Tauri

Authors: Marc C. White, Geoffrey V. Bicknell, Ralph S. Sutherland, Raquel Salmeron, Peter J. McGregor

Abstract: Turbulent entrainment processes may play an important role in the outflows from young stellar objects at all stages of their evolution. In particular, lateral entrainment of ambient material by high-velocity, well-collimated protostellar jets may be the cause of the multiple emission-line velocity components observed in the microjet-scale outflows driven by classical T Tauri stars. Intermediate-ve… ▽ More Turbulent entrainment processes may play an important role in the outflows from young stellar objects at all stages of their evolution. In particular, lateral entrainment of ambient material by high-velocity, well-collimated protostellar jets may be the cause of the multiple emission-line velocity components observed in the microjet-scale outflows driven by classical T Tauri stars. Intermediate-velocity outflow components may be emitted by a turbulent, shock- excited mixing layer along the boundaries of the jet. We present a formalism for describing such a mixing layer based on Reynolds decomposition of quantities measuring fundamental properties of the gas. In this model, the molecular wind from large disc radii provides a continual supply of material for entrainment. We calculate the total stress profile in the mixing layer, which allows us to estimate the dissipation of turbulent energy, and hence the luminosity of the layer. We utilize MAPPINGS IV shock models to determine the fraction of total emission that occurs in [Fe II] 1.644 μm line emission in order to facilitate comparison to previous observations of the young stellar object DG Tauri. Our model accurately estimates the luminosity and changes in mass outflow rate of the intermediate-velocity component of the DG Tau approaching outflow. Therefore, we propose that this component represents a turbulent mixing layer surrounding the well-collimated jet in this object. Finally, we compare and contrast our model to previous work in the field. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 18 pages, 13 figures, accepted for publication in MNRAS

The mass evolution of the first galaxies: stellar mass functions and star formation rates at $4 < z < 7$ in the CANDELS GOODS-South field

Authors: Kenneth Duncan, Christopher J. Conselice, Alice Mortlock, William G. Hartley, Yicheng Guo, Henry C. Ferguson, Romeel Davé, Yu Lu, Jamie Ownsworth, Matthew L. N. Ashby, Avishai Dekel, Mark Dickinson, Sandra M. Faber, Mauro Giavalisco, Norman A. Grogin, Dale Kocevski, Anton M. Koekemoer, Rachel S. Somerville, Catherine E. White

Abstract: We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributi… ▽ More We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($α\approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} ρ_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($β$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift. △ Less

Submitted 11 August, 2014; originally announced August 2014.

Comments: 28 pages, 23 figures, 2 appendices. Accepted for publication in MNRAS, August 7 2014

A parametric study of possible solutions to the high-redshift overproduction of stars in modeled dwarf galaxies

Authors: Catherine E. White, Rachel S. Somerville, Henry C. Ferguson

Abstract: Both numerical hydrodynamic and semi-analytic cosmological models of galaxy formation struggle to match observed star formation histories of galaxies in low mass halos (M$_{\rm{H}} \lesssim 10^{11} M_\odot$), predicting more star formation at high redshift and less star formation at low redshift than observed. The fundamental problem is that galaxies' gas accretion and star formation rates are too… ▽ More Both numerical hydrodynamic and semi-analytic cosmological models of galaxy formation struggle to match observed star formation histories of galaxies in low mass halos (M$_{\rm{H}} \lesssim 10^{11} M_\odot$), predicting more star formation at high redshift and less star formation at low redshift than observed. The fundamental problem is that galaxies' gas accretion and star formation rates are too closely coupled in the models: the accretion rate largely drives the star formation rate. Observations point to gas accretion rates that outpace star formation at high redshift, resulting in a buildup of gas and a delay in star formation until lower redshifts. We present three empirical adjustments of standard recipes in a semi-analytic model motivated by three physical scenarios that could cause this decoupling: 1) the mass-loading factors of outflows driven by stellar feedback may have a steeper dependence on halo mass at earlier times, 2) the efficiency of star formation may be lower in low mass halos at high redshift, and 3) gas may not be able to accrete efficiently onto the disk in low mass halos at high redshift. These new recipes, once tuned, better reproduce the evolution of $f_\star \equiv M_\star/M_{\rm{H}}$ as a function of halo mass as derived from abundance matching over redshifts $z=0$ to $3$, though they have different effects on cold gas fractions, star formation rates, and metallicities. Changes to gas accretion and stellar-driven winds are promising, while direct modification of the star formation timescale requires drastic measures that are not physically well-motivated. △ Less

Submitted 31 December, 2014; v1 submitted 7 July, 2014; originally announced July 2014.

Comments: 27 pages, 29 figures

Observation of muon intensity variations by season with the MINOS Near Detector

Authors: P. Adamson, I. Anghel, A. Aurisano, G. Barr, M. Bishai, A. Blake, G. J. Bock, D. Bogert, S. V. Cao, C. M. Castromonte, S. Childress, J. A. B. Coelho, L. Corwin, D. Cronin-Hennessy, J. K. de Jong, A. V. Devan, N. E. Devenish, M. V. Diwan, C. O. Escobar, J. J. Evans, E. Falk, G. J. Feldman, T. H. Fields, M. V. Frohne, H. R. Gallagher , et al. (87 additional authors not shown)

Abstract: A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $α_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428… ▽ More A sample of 1.53$\times$10$^{9}$ cosmic-ray-induced single muon events has been recorded at 225 meters-water-equivalent using the MINOS Near Detector. The underground muon rate is observed to be highly correlated with the effective atmospheric temperature. The coefficient $α_{T}$, relating the change in the muon rate to the change in the vertical effective temperature, is determined to be 0.428$\pm$0.003(stat.)$\pm$0.059(syst.). An alternative description is provided by the weighted effective temperature, introduced to account for the differences in the temperature profile and muon flux as a function of zenith angle. Using the latter estimation of temperature, the coefficient is determined to be 0.352$\pm$0.003(stat.)$\pm$0.046(syst.). △ Less

Submitted 26 June, 2014; originally announced June 2014.

Comments: 9 pages 10 figures

Report number: FERMILAB-PUB-14-209

Journal ref: Phys. Rev. D 90, 012010 (2014)

Slow-Speed Supernovae from the Palomar Transient Factory: Two Channels

Authors: Christopher J. White, Mansi M. Kasliwal, Peter E. Nugent, Avishay Gal-Yam, D. Andrew Howell, Mark Sullivan, Ariel Goobar, Anthony L. Piro, Joshua S. Bloom, Shrinivas R. Kulkarni, Russ R. Laher, Frank Masci, Eran O. Ofek, Jason Surace, Sagi Ben-Ami, Yi Cao, S. Bradley Cenko, Isobel M. Hook, Jakob Jönsson, Thomas Matheson, Assaf Sternberg, Robert M. Quimby, Ofer Yaron

Abstract: Since the discovery of the unusual prototype SN 2002cx, the eponymous class of low-velocity, hydrogen-poor supernovae has grown to include at most another two dozen members identified from several heterogeneous surveys, in some cases ambiguously. Here we present the results of a systematic study of 1077 hydrogen-poor supernovae discovered by the Palomar Transient Factory, leading to nine new membe… ▽ More Since the discovery of the unusual prototype SN 2002cx, the eponymous class of low-velocity, hydrogen-poor supernovae has grown to include at most another two dozen members identified from several heterogeneous surveys, in some cases ambiguously. Here we present the results of a systematic study of 1077 hydrogen-poor supernovae discovered by the Palomar Transient Factory, leading to nine new members of this peculiar class. Moreover we find there are two distinct subclasses based on their spectroscopic, photometric, and host galaxy properties: The "SN 2002cx-like" supernovae tend to be in later-type or more irregular hosts, have more varied and generally dimmer luminosities, have longer rise times, and lack a Ti II trough when compared to the "SN 2002es-like" supernovae. None of our objects show helium, and we counter a previous claim of two such events. We also find that these transients comprise 5.6+17-3.7% (90% confidence) of all SNe Ia, lower compared to earlier estimates. Combining our objects with the literature sample, we propose that these subclasses have two distinct physical origins. △ Less

Submitted 28 May, 2014; originally announced May 2014.

Comments: 49 pages, 36 figures, submitted to ApJ

Multi-epoch Sub-arcsecond [Fe II] Spectroimaging of the DG Tau Outflows with NIFS. II. On the Nature of the Bipolar Outflow Asymmetry

Authors: Marc C White, Geoffrey V Bicknell, Peter J McGregor, Raquel Salmeron

Abstract: The origin of bipolar outflow asymmetry in young stellar objects (YSOs) remains poorly understood. It may be due to an intrinsically asymmetric outflow launch mechanism, or it may be caused by the effects of the ambient medium surrounding the YSO. Answering this question is an important step in understanding outflow launching. We have investigated the bipolar outflows driven by the T Tauri star DG… ▽ More The origin of bipolar outflow asymmetry in young stellar objects (YSOs) remains poorly understood. It may be due to an intrinsically asymmetric outflow launch mechanism, or it may be caused by the effects of the ambient medium surrounding the YSO. Answering this question is an important step in understanding outflow launching. We have investigated the bipolar outflows driven by the T Tauri star DG Tauri on scales of hundreds of AU, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The approaching outflow consists of a well-collimated jet, nested within a lower-velocity disc wind. The receding outflow is composed of a single-component bubble-like structure. We analyse the kinemat- ics of the receding outflow using kinetic models, and determine that it is a quasi-stationary bubble with an expanding internal velocity field. We propose that this bubble forms because the receding counterjet from DG Tau is obstructed by a clumpy ambient medium above the circumstellar disc surface, based on similarities between this structure and those found in the modeling of active galactic nuclei outflows. We find evidence of interaction between the obscured counterjet and clumpy ambient material, which we attribute to the large molecular envelope around the DG Tau system. An analytical model of a momentum-driven bubble is shown to be consistent with our interpretation. We conclude that the bipolar outflow from DG Tau is intrinsically symmetric, and the observed asymmetries are due to environmental effects. This mechanism can potentially be used to explain the observed bipolar asymmetries in other YSO outflows. △ Less

Submitted 23 April, 2014; originally announced April 2014.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Multi-epoch Sub-arcsecond [Fe II] Spectroimaging of the DG Tau Outflows with NIFS. I. First data epoch

Authors: Marc C. White, Peter J. McGregor, Geoffrey V. Bicknell, Raquel Salmeron, Tracy L. Beck

Abstract: Investigating the outflows emanating from young stellar objects (YSOs) on sub-arcsecond scales provides important clues to the nature of the underlying accretion-ejection process occurring near the central protostar. We have investigated the structures and kinematics of the outflows driven by the YSO DG Tauri, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The blueshif… ▽ More Investigating the outflows emanating from young stellar objects (YSOs) on sub-arcsecond scales provides important clues to the nature of the underlying accretion-ejection process occurring near the central protostar. We have investigated the structures and kinematics of the outflows driven by the YSO DG Tauri, using the Near-infrared Integral Field Spectrograph (NIFS) on Gemini North. The blueshifted outflow shows two distinct components in [Fe II] 1.644 micron emission, which are separated using multi-component line fitting. A stationary recollimation shock is observed, in agreement with previous X-ray and FUV observations. The presence of this shock indicates that the innermost streamlines of the high-velocity component are launched at a very small radius, 0.01-0.15 AU, from the central star. The jet accelerates and expands downstream of the recollimation shock; the 'acceleration' is likely a sign of velocity variations in the jet. No evidence of rotation is found, and we compare this non-detection to previous counter-claims. Moving jet knots, likely the result of the jet velocity variations, are observed. One of these knots moves more slowly than previously observed knots, and the knot ejection interval appears to be non-periodic. An intermediate-velocity component surrounds this central jet, and is interpreted as the result of a turbulent mixing layer along the jet boundaries. Such lateral entrainment requires the presence of a magnetic field of strength a few mG or less at hundreds of AU above the disc surface, which is argued to be a reasonable proposition. In H2 1-0 S(1) 2.1218 micron emission, a wide-angle, intermediate-velocity blueshifted outflow is observed. Both outflows are consistent with being launched by a magnetocentrifugal disc wind, although an X-wind origin for the high-velocity jet cannot be ruled out. △ Less

Submitted 2 April, 2014; originally announced April 2014.

Comments: 28 pages, 19 figures, accepted for publication in MNRAS

Stability of metal-rich very massive stars

Authors: J. Goodman, Christopher J. White

Abstract: We revisit the stability of very massive nonrotating main-sequence stars at solar metallicity, with the goal of understanding whether radial pulsations set a physical upper limit to stellar mass. Models of up to 938 solar masses are constructed with the Mesa code, and their linear stability in the fundamental mode, assumed to be the most dangerous, is analysed with a fully nonadiabatic method. Mod… ▽ More We revisit the stability of very massive nonrotating main-sequence stars at solar metallicity, with the goal of understanding whether radial pulsations set a physical upper limit to stellar mass. Models of up to 938 solar masses are constructed with the Mesa code, and their linear stability in the fundamental mode, assumed to be the most dangerous, is analysed with a fully nonadiabatic method. Models above 100 MSun have extended tenuous atmospheres ("shelves") that affect the stability of the fundamental. Even when positive, this growth rate is small, in agreement with previous results. We argue that small growth rates lead to saturation at small amplitudes that are not dangerous to the star. A mechanism for saturation is demonstrated involving nonlinear parametric coupling to short-wavelength g modes and the damping of the latter by radiative diffusion. The shelves are subject to much more rapidly growing strange modes. This also agrees with previous results but is extended here to higher masses. The strange modes probably saturate via shocks rather than mode coupling but have very small amplitudes in the core, where almost all of the stellar mass resides. Although our stellar models are hydrostatic, the structure of their outer parts suggests that optically thick winds, driven by some combination of radiation pressure, transsonic convection, and strange modes, are more likely than pulsation in the fundamental mode to limit the main-sequence lifetime. △ Less

Submitted 17 November, 2015; v1 submitted 27 March, 2014; originally announced March 2014.

Comments: 14 pages, 8 figures, 1 appendix; this version to be published in MNRAS

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

Authors: Yu Lu, Risa H. Wechsler, Rachel S. Somerville, Darren Croton, Lauren Porter, Joel Primack, Peter S. Behroozi, Henry C. Ferguson, David C. Koo, Yicheng Guo, Mohammadtaher Safarzadeh, Kristian Finlator, Marco Castellano, Catherine E. White, Veronica Sommariva, Chris Moody

Abstract: We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian… ▽ More We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED △ Less

Submitted 11 December, 2013; originally announced December 2013.

Comments: 29 pages, 14 figures, 1 table, submitted to ApJ

The MICA Experiment: Astrophysics in Virtual Worlds

Authors: S. G. Djorgovski, Piet Hut, Rob Knop, Giuseppe Longo, Steve McMillan, Enrico Vesperini, Ciro Donalek, Matthew Graham, Asish Mahabal, Franz Sauer, Charles White, Crista Lopes

Abstract: We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. W… ▽ More We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted. △ Less

Submitted 28 January, 2013; originally announced January 2013.

Comments: 10 pages, 8 figures; invited paper for the refereed proc. of the SLACTIONS 2012 conference

Journal ref: Ref. proc. of the SLACTIONS 2012 International Research Conference on Virtual Worlds,ceds. L. Morgado, Y. Sivan, A.M. Maia, et al., pp. 49-58 (2012). Vila Real, PT: UTAD, ISBN 978-989-704-102-0, avail. at http://hdl.handle.net/10348/2305

DBI Inflation using a One-Parameter Family of Throat Geometries

Authors: Florian Gmeiner, Chris D. White

Abstract: We demonstrate the possibility of examining cosmological signatures in the DBI inflation setup using the BGMPZ solution, a one-parameter family of geometries for the warped throat which interpolate between the Maldacena-Nunez and Klebanov-Strassler solutions. The warp factor is determined numerically and subsequently used to calculate cosmological observables including the scalar and tensor spec… ▽ More We demonstrate the possibility of examining cosmological signatures in the DBI inflation setup using the BGMPZ solution, a one-parameter family of geometries for the warped throat which interpolate between the Maldacena-Nunez and Klebanov-Strassler solutions. The warp factor is determined numerically and subsequently used to calculate cosmological observables including the scalar and tensor spectral indices, for a sample point in the parameter space. As one moves away from the KS solution for the throat the warp factor is qualitatively different, which leads to a significant change for the observables, but also generically increases the non-Gaussianity of the models. We argue that the different models can potentially be differentiated by current and future experiments. △ Less

Submitted 6 April, 2008; v1 submitted 10 October, 2007; originally announced October 2007.

Comments: 17 pages, 10 figures; v2: section 4 expanded, references added; v3: typos fixed

Report number: NIKHEF/2007-022

Journal ref: JCAP 0802:012,2008