Showing 1–50 of 54 results for author: Hartigan, P

PROJECT-J: JWST observations of HH46~IRS and its outflow. Overview and first results

Authors: B. Nisini, M. G. Navarro, T. Giannini, S. Antoniucci, P. J. Kavanagh, P. Hartigan, F. Bacciotti, A. Caratti o Garatti, A. Noriega Crespo, E. van Dishoek, E. Whelan, H. G. Arce, S. Cabrit, D. Coffey, D. Fedele, J. Eisloeffel, M. E. Palumbo, L. Podio, T. P. Ray, M. Schultze, R. G. Urso, J. M. Alcala', M. A. Bautista, C. Codella, T. G. Greene , et al. (1 additional authors not shown)

Abstract: We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We… ▽ More We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We detect, for the first time, the red-shifted jet within ~ 90 au from the source. Dozens of shock-excited forbidden lines are observed, including highly ionized species such as [Ne III] 15.5 micron, suggesting that the gas is excited by high velocity (> 80 km/s) shocks in a relatively high density medium. Images of H2 lines at different excitations outline a complex molecular flow, where a bright cavity, molecular shells, and a jet-driven bow-shock interact with and are shaped by the ambient conditions. Additional NIRCam 2 micron images resolve the HH46 IRS ~ 110 au binary system and suggest that the large asymmetries observed between the jet and the H2 wide angle emission could be due to two separate outflows being driven by the two sources. The spectra of the unresolved binary show deep ice bands and plenty of gaseous lines in absorption, likely originating in a cold envelope or disk. In conclusion, JWST has unraveled for the first time the origin of the HH46 IRS complex outflow demonstrating its capability to investigate embedded regions around young stars, which remain elusive even at near-IR wavelengths. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 28 pages, 15 figures, Accepted for publication on The Astrophysical Journal (9 April 2024)

PDRs4All VIII: Mid-IR emission line inventory of the Orion Bar

Authors: Dries Van De Putte, Raphael Meshaka, Boris Trahin, Emilie Habart, Els Peeters, Olivier Berné, Felipe Alarcón, Amélie Canin, Ryan Chown, Ilane Schroetter, Ameek Sidhu, Christiaan Boersma, Emeric Bron, Emmanuel Dartois, Javier R. Goicoechea, Karl D. Gordon, Takashi Onaka, Alexander G. G. M. Tielens, Laurent Verstraete, Mark G. Wolfire, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Jan Cami, Sara Cuadrado , et al. (113 additional authors not shown)

Abstract: Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observat… ▽ More Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Journal ref: A&A 687, A86 (2024)

A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk

Authors: Olivier Berné, Emilie Habart, Els Peeters, Ilane Schroetter, Amélie Canin, Ameek Sidhu, Ryan Chown, Emeric Bron, Thomas J. Haworth, Pamela Klaassen, Boris Trahin, Dries Van De Putte, Felipe Alarcón, Marion Zannese, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Jan Cami, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea , et al. (121 additional authors not shown)

Abstract: Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of… ▽ More Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Journal ref: Science, 383, 6686, 2024

High-quality Extragalactic Legacy-field Monitoring (HELM) with DECam

Authors: Ming-Yang Zhuang, Qian Yang, Yue Shen, Monika Adamow, Douglas N. Friedel, R. A. Gruendl, Xin Liu, Paul Martini, Timothy M. C. Abbott, Scott F. Anderson, Roberto J. Assef, Franz E. Bauer, Rich Bielby, W. N. Brandt, Colin J. Burke, Jorge Casares, Yu-Ching Chen, Gisella De Rosa, Alex Drlica-Wagner, Tom Dwelly, Alice Eltvedt, Gloria Fonseca Alvarez, Jianyang Fu, Cesar Fuentes, Melissa L. Graham , et al. (23 additional authors not shown)

Abstract: High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the… ▽ More High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the $(u)gri(z)$ bands, over a collective sky area of $\sim 38$ deg${\rm ^2}$. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first $\sim 3.5$ years of light curves during 2019A--2022A. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 11 pages, 8 figures, 4 tables, submitted to ApJS. Median source catalogs and light curves of individual objects are publicly available at https://ariel.astro.illinois.edu/helm/

Cooling and Instabilities in Colliding Radiative Flows with Toroidal Magnetic Fields

Authors: R. N. Markwick, A. Frank, E. G. Blackman, J. Carroll-Nellenback, S. V. Lebedev, D. R. Russell, J. W. D. Halliday, L. G. Suttle, P. M. Hartigan

Abstract: We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calcul… ▽ More We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calculation. We systematically vary plasma beta ($β_m$) in the flows, the strength of the cooling ($Λ_0$) and the exponent ($α$) of temperature-dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a "spine" is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depends on initial flow parameters and provide a simple analytic model that captures the basic features of the flow. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 14 pages, 16 figures. Submitted to MNRAS

PDRs4All VI: Probing the Photochemical Evolution of PAHs in the Orion Bar Using Machine Learning Techniques

Authors: S. Pasquini, E. Peeters, B. Schefter, B. Khan, A. Sidhu, R. Chown, J. Cami, A. Tielens, F. Alarcon, A. Canin, I. Schroetter, B. Trahin, D. Van De Putte, C. Boersma, E. Dartois, T. Onaka, A. Candian, P. Hartigan, T. S. -Y. Lai, G. Rouille, D. A. Sales, Y. Zhang, E. Habart, O. Berne

Abstract: [Abridged] JWST observations of the Orion Bar have shown the incredible richness of PAH bands and their variation on small scales. We aim to probe the photochemical evolution of PAHs across the key zones of the photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning. We use NIRSpec and MIRI IFU data from the JWST ERS Program PDRs4All. We lever bisecting k-means clu… ▽ More [Abridged] JWST observations of the Orion Bar have shown the incredible richness of PAH bands and their variation on small scales. We aim to probe the photochemical evolution of PAHs across the key zones of the photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning. We use NIRSpec and MIRI IFU data from the JWST ERS Program PDRs4All. We lever bisecting k-means clustering to generate detailed spatial maps of the spectral variability in several wavelength regions. We discuss the variations in the cluster profiles and connect them to the local physical conditions. We interpret these variations with respect to the key zones: the HII region, the atomic PDR zone, and the three dissociation fronts. The PAH emission exhibits spectral variation that depends strongly on spatial position in the PDR. We find the 8.6um band to behave differently than all other bands which vary systematically with one another. We find uniform variation in the 3.4-3.6um bands and 3.4/3.3 intensity ratio. We attribute the carrier of the 3.4-3.6um bands to a single side group attached to very similarly sized PAHs. Cluster profiles reveal a transition between characteristic profiles classes of the 11.2um feature from the atomic to the molecular PDR zone. We find the carriers of each of the profile classes to be independent, and reason the latter to be PAH clusters existing solely deep in the molecular PDR. Clustering also reveals a connection between the 11.2 and 6.2um bands; and that clusters generated from variation in the 10.9-11.63um region can be used to recover those in the 5.95-6.6um region. Clustering is a powerful tool for characterizing PAH variability on both spatial and spectral scales. For individual bands as well as global spectral behaviours, we find UV-processing to be the most important driver of the evolution of PAHs and their spectral signatures in the Orion Bar. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 28 pages, 20 figures, submitted to A&A

Journal ref: A&A 685, A77 (2024)

PDRs4All III: JWST's NIR spectroscopic view of the Orion Bar

Authors: Els Peeters, Emilie Habart, Olivier Berne, Ameek Sidhu, Ryan Chown, Dries Van De Putte, Boris Trahin, Ilane Schroetter, Amelie Canin, Felipe Alarcon, Bethany Schefter, Baria Khan, Sofia Pasquini, Alexander G. G. M. Tielens, Mark G. Wolfire, Emmanuel Dartois, Javier R. Goicoechea, Alexandros Maragkoudakis, Takashi Onaka, Marc W. Pound, Silvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma , et al. (113 additional authors not shown)

Abstract: (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion… ▽ More (Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 52 pages, 30 figures, submitted to A&A

Journal ref: A&A 685, A74 (2024)

PDRs4All IV. An embarrassment of riches: Aromatic infrared bands in the Orion Bar

Authors: Ryan Chown, Ameek Sidhu, Els Peeters, Alexander G. G. M. Tielens, Jan Cami, Olivier Berné, Emilie Habart, Felipe Alarcón, Amélie Canin, Ilane Schroetter, Boris Trahin, Dries Van De Putte, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Sara Cuadrado, Emmanuel Dartois, Daniel Dicken, Meriem El-Yajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (114 additional authors not shown)

Abstract: (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory o… ▽ More (Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $μ$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $μ$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $μ$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'. △ Less

Submitted 5 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: 25 pages, 10 figures, to appear in A&A

Journal ref: A&A 685, A75 (2024)

PDRs4All II: JWST's NIR and MIR imaging view of the Orion Nebula

Authors: Emilie Habart, Els Peeters, Olivier Berné, Boris Trahin, Amélie Canin, Ryan Chown, Ameek Sidhu, Dries Van De Putte, Felipe Alarcón, Ilane Schroetter, Emmanuel Dartois, Sílvia Vicente, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Christiaan Boersma, Emeric Bron, Jan Cami, Sara Cuadrado, Daniel Dicken, Meriem Elyajouri, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Lina Issa , et al. (117 additional authors not shown)

Abstract: The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation fron… ▽ More The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate. △ Less

Submitted 2 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Journal ref: A&A 685, A73 (2024)

Young stellar objects, accretion disks, and their variability with Rubin Observatory LSST

Authors: R. Bonito, L. Venuti, S. Ustamujic, P. Yoachim, R. A. Street, L. Prisinzano, P. Hartigan, M. G. Guarcello, K. G. Stassun, T. Giannini, E. D. Feigelson, A. Caratti o Garatti, S. Orlando, W. I. Clarkson, P. McGehee, E. C. Bellm, J. E. Gizis

Abstract: Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, ca… ▽ More Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multi-band filters option available in Rubin LSST, in particular the $u,g,r,i$ filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g. field crowdedness, massive neighbors, metallicity), by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide-Fast-Deep observing mode that will be used to survey the sky over the full duration of the main survey ($\approx$10 years). The combination of these two modes will be vital to investigate the connection between the inner disk dynamics and longer-term eruptive variability behaviors, such as those observed on EXor objects. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 11 pages, 4 figures, 1 table; accepted for publication in The Astrophysical Journal Supplement Series

Morphology of Shocked Lateral Outflows in Colliding Hydrodynamic Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, E. G. Blackman, P. M. Hartigan, S. V. Lebedev, D. R. Russel, J. W. D. Halliday, L. G. Suttle

Abstract: Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding f… ▽ More Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding flows %radius in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows (SLOs) being deflected from the horizontal when the flows are of differing cross-section. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Journal ref: Physics of Plasmas, Volume 29, Issue 10, id.102901, 14 pp (October 2022)

Length-scales and Dynamics of Carina's Western Wall

Authors: Turlough Downes, Patrick Hartigan, Andrea Isella

Abstract: We present a variety of analyses of the turbulent dynamics of the boundary of a photo-dissociation region (PDR) in the Carina Nebula using high resolution ALMA observations. Using Principal Component Analysis we suggest that the turbulence in this molecular cloud is driven at large scales. Analysis of the centroid velocity structure functions indicate that the turbulence is dominated by shocks rat… ▽ More We present a variety of analyses of the turbulent dynamics of the boundary of a photo-dissociation region (PDR) in the Carina Nebula using high resolution ALMA observations. Using Principal Component Analysis we suggest that the turbulence in this molecular cloud is driven at large scales. Analysis of the centroid velocity structure functions indicate that the turbulence is dominated by shocks rather than local (in k-space) transport of energy. We further find that length-scales in the range 0.02 - 0.03 pc are important in the dynamics of this cloud and this finding is supported by analysis of the dominant emission structure length-scale. These length-scales are well resolved by the observational data and we conclude that the apparent importance of this range of scales is physical in origin. Given that it is also well within the range strongly influenced by ambipolar diffusion, we conclude that it is not primarily a product of turbulence alone, but is more likely to be a result of the interplay between gravity and turbulence. Finally, through comparison of these results with previous observations of H2 emission from the Western Wall we demonstrate that observations of a PDR can be used to probe the internal structure of the undisturbed portion of a molecular cloud. △ Less

Submitted 5 December, 2022; v1 submitted 1 December, 2022; originally announced December 2022.

Comments: 13 pages, 9 figures. Accepted to MNRAS

ALMA Datacubes and Continuum Maps of the Irradiated Western Wall in Carina

Authors: Patrick Hartigan, Maxwell Hummel, Andrea Isella, Turlough Downes

Abstract: We present ALMA observations of the continuum and line emission of $^{12}$CO, $^{13}$CO, C$^{18}$O, and [C I] for a portion of the G287.38-0.62 (Car 1-E) region in the Carina star-forming complex. The new data record how a molecular cloud responds on subarcsecond scales when subjected to a powerful radiation front, and provide insights into the overall process of star formation within regions that… ▽ More We present ALMA observations of the continuum and line emission of $^{12}$CO, $^{13}$CO, C$^{18}$O, and [C I] for a portion of the G287.38-0.62 (Car 1-E) region in the Carina star-forming complex. The new data record how a molecular cloud responds on subarcsecond scales when subjected to a powerful radiation front, and provide insights into the overall process of star formation within regions that contain the most massive young stars. The maps show several molecular clouds superpose upon the line of sight, including a portion of the Western Wall, a highly-irradiated cloud situated near the young star cluster Trumpler 14. In agreement with theory, there is a clear progression from fluoresced H$_2$, to [C I], to C$^{18}$O with distance into the PDR front. Emission from optically thick $^{12}$CO extends across the region, while $^{13}$CO, [C I] and especially C$^{18}$O are more optically thin, and concentrate into clumps and filaments closer to the PDR interface. Within the Western Wall cloud itself we identify 254 distinct core-sized clumps in our datacube of C$^{18}$O. The mass distribution of these objects is similar to that of the stellar IMF. Aside from a large-scale velocity gradient, the clump radial velocities lack any spatial coherence size. There is no direct evidence for triggering of star formation in the Western Wall in that its C$^{18}$O clumps and continuum cores appear starless, with no pillars present. However, the densest portion of the cloud lies closest to the PDR, and the C$^{18}$O emission is flattened along the radiation front. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: accepted for publication in the Astronomical Journal. 16 figures and 2 animations. Animations available at https://sparky.rice.edu/~hartigan/pub.html

Rubin Observatory LSST Transients and Variable Stars Roadmap

Authors: Kelly M. Hambleton, Federica B. Bianco, Rachel Street, Keaton Bell, David Buckley, Melissa Graham, Nina Hernitschek, Michael B. Lund, Elena Mason, Joshua Pepper, Andrej Prsa, Markus Rabus, Claudia M. Raiteri, Robert Szabo, Paula Szkody, Igor Andreoni, Simone Antoniucci, Barbara Balmaverde, Eric Bellm, Rosaria Bonito, Giuseppe Bono, Maria Teresa Botticella, Enzo Brocato, Katja Bucar Bricman, Enrico Cappellaro , et al. (57 additional authors not shown)

Abstract: The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the T… ▽ More The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 202 pages (in book format) 34 figures plus chapter heading figures (13)

SpinSpotter: An Automated Algorithm for Identifying Stellar Rotation Periods With Autocorrelation Analysis

Authors: Rae J. Holcomb, Paul Robertson, Patrick Hartigan, Ryan J. Oelkers, Caleb Robinson

Abstract: Spinspotter is a robust and automated algorithm designed to extract stellar rotation periods from large photometric datasets with minimal supervision. Our approach uses the autocorrelation function (ACF) to identify stellar rotation periods up to one-third the observational baseline of the data. Our algorithm also provides a suite of diagnostics that describe the features in the ACF, which allows… ▽ More Spinspotter is a robust and automated algorithm designed to extract stellar rotation periods from large photometric datasets with minimal supervision. Our approach uses the autocorrelation function (ACF) to identify stellar rotation periods up to one-third the observational baseline of the data. Our algorithm also provides a suite of diagnostics that describe the features in the ACF, which allows the user to fine-tune the tolerance with which to accept a period detection. We apply it to approximately 130,000 main-sequence stars observed by the Transiting Exoplanet Survey Satellite (TESS) at 2-minute cadence during Sectors 1-26, and identify rotation periods for 13,504 stars ranging from 0.4 to 14 days. We demonstrate good agreement between our sample and known values from the literature and note key differences between our population of rotators and those previously identified in the Kepler field, most notably a large population of fast-rotating M dwarfs. Our sample of rotating stars provides a data set with coverage of nearly the entire sky that can be used as a basis for future gyrochronological studies, and, when combined with proper motions and distances from Gaia, to search for regions with high densities of young stars, thus identifying areas of recent star formation and undiscovered moving group members. Our algorithm is publicly available for download and use on GitHub. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 14 pages, 5 figures, Accepted for publication in The Astrophysical Journal

PDRs4All: A JWST Early Release Science Program on radiative feedback from massive stars

Authors: Olivier Berné, Émilie Habart, Els Peeters, Alain Abergel, Edwin A. Bergin, Jeronimo Bernard-Salas, Emeric Bron, Jan Cami, Stéphanie Cazaux, Emmanuel Dartois, Asunción Fuente, Javier R. Goicoechea, Karl D. Gordon, Yoko Okada, Takashi Onaka, Massimo Robberto, Markus Röllig, Alexander G. G. M. Tielens, Silvia Vicente, Mark G. Wolfire, Felipe Alarcon, C. Boersma, Ameélie Canin, Ryan Chown, Daniel Dicken , et al. (112 additional authors not shown)

Abstract: Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation… ▽ More Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: Submitted to PASP

Cooling and Instabilities in Colliding Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, B. Liu, E. G. Blackman, S. V. Lebedev, P. M. Hartigan

Abstract: Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature… ▽ More Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period. △ Less

Submitted 7 September, 2021; originally announced September 2021.

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

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 2, pp.2266-2278 (Dec 2021)

Probing jets from young embedded sources: clues from HST near-IR [Fe II] images

Authors: Jessica Erkal, Brunella Nisini, Deirdre Coffey, Francesca Bacciotti, Patrick Hartigan, Simone Antoniucci, Teresa Giannini, Jochen Eislöffel, Carlo Felice Manara

Abstract: We present near-infrared [Fe II] images of four Class 0/I jets (HH 1/2, HH 34, HH 111, HH 46/47) observed with the Hubble Space Telescope Wide Field Camera 3. The unprecedented angular resolution allows us to measure proper motions, jet widths and trajectories, and extinction along the jets. In all cases, we detect the counter-jet which was barely visible or invisible at shorter wavelengths. We me… ▽ More We present near-infrared [Fe II] images of four Class 0/I jets (HH 1/2, HH 34, HH 111, HH 46/47) observed with the Hubble Space Telescope Wide Field Camera 3. The unprecedented angular resolution allows us to measure proper motions, jet widths and trajectories, and extinction along the jets. In all cases, we detect the counter-jet which was barely visible or invisible at shorter wavelengths. We measure tangential velocities of a few hundred km/s, consistent with previous HST measurements over 10 years ago. We measure the jet width as close as a few tens of au from the star, revealing high collimations of about 2 degrees for HH 1, HH 34, HH 111 and about 8 degrees for HH 46, all of which are preserved up to large distances. For HH 34, we find evidence of a larger initial opening angle of about 7 degrees. Measurement of knot positions reveals deviations in trajectory of both the jet and counter-jet of all sources. Analysis of asymmetries in the inner knot positions for HH 111 suggests the presence of a low mass stellar companion at separation 20-30 au. Finally, we find extinction values of 15-20 mag near the source which gradually decreases moving downstream along the jet. These observations have allowed us to study the counter-jet at unprecedented high angular resolution, and will be a valuable reference for planning future JWST mid-infrared observations which will peer even closer into the jet engine. △ Less

Submitted 25 June, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in ApJ - 26 pages, 22 figures, 9 tables

Emission-Line Datacubes of the HH 32 Stellar Jet

Authors: Patrick Hartigan, Lynne Hillenbrand, Matuesz Matuszewski, Arlindo Chan Borges, James D. Neill, D. Christopher Martin, Patrick Morrissey, Anna M. Moore

Abstract: We analyze datacubes of over 60 emission lines in the HH 32 stellar jet acquired with the Keck Cosmic Web Imager (KCWI). The data cover the less-explored blue portion of the spectrum between 3586A and 6351A , and have both high spectral (R $\sim$ 10000) and spatial ($\lesssim$ 1" ) resolution. The study includes all three major ionization states of oxygen, three Balmer lines, multiple lines of Fe… ▽ More We analyze datacubes of over 60 emission lines in the HH 32 stellar jet acquired with the Keck Cosmic Web Imager (KCWI). The data cover the less-explored blue portion of the spectrum between 3586A and 6351A , and have both high spectral (R $\sim$ 10000) and spatial ($\lesssim$ 1" ) resolution. The study includes all three major ionization states of oxygen, three Balmer lines, multiple lines of Fe II and Fe III, as well as the first datacubes ever acquired for important unblended diagnostic lines such as He II $λ$4686, Ca I $λ$3933, and Mg I] $λ$4571. The data cubes generally sort according to excitation, and have a relatively continuous progression from the highest-excitation ions (He II, O III) through the intermediate-excitation ions (O I and H I) to the lowest-excitation ions (Ca II and Mg I). Merging the KCWI cubes with HST images leads to several new insights about the flow, including evidence for bow shocks, partial bow shocks, spur shocks, Mach disks, jet deflection shocks, a wiggling jet, and potential shock precursors. The most surprising result is that one of the velocity components of Fe II in the Mach disk suddenly increases in flux relative to other lines by a factor of two, implying that the Mach disk vaporizes dust in the jet. Hence, jets must accelerate or entrain dust to speeds of over 300 km/s without destroying the grains. △ Less

Submitted 20 August, 2020; v1 submitted 8 August, 2020; originally announced August 2020.

Comments: 16 figures, 1 movie. Accepted in the Astronomical Journal. Movie available at http://sparky.rice.edu/~hartigan/pub.html

A JWST Preview: Adaptive-Optics Images of H$_2$, Br-$γ$, and K-continuum in Carina's Western Wall

Authors: Patrick Hartigan, Turlough Downes, Andrea Isella

Abstract: We present the first wide-field near-infrared adaptive optics images of Carina's Western Wall (G287.38-0.62), one of the brightest and most well-defined irradiated interfaces known in a region of massive star formation. The new narrowband H$_2$ 2.12$μ$m, Br-$γ$ and K-continuum images trace the photoevaporative flow from the cloud and identify locations where UV-radiation from the surrounding massi… ▽ More We present the first wide-field near-infrared adaptive optics images of Carina's Western Wall (G287.38-0.62), one of the brightest and most well-defined irradiated interfaces known in a region of massive star formation. The new narrowband H$_2$ 2.12$μ$m, Br-$γ$ and K-continuum images trace the photoevaporative flow from the cloud and identify locations where UV-radiation from the surrounding massive stars excites molecular hydrogen to fluoresce. With a field of view of $\sim$ 1.5' $\times$ 2.9' and spatial resolution between 60 $-$ 110 mas, the new images show a spectacular level of detail over a large area, and presage what JWST should achieve. The Wall is convex in shape, with a large triangular-shaped extension near its apex. The interface near the apex consists of 3 $-$ 4 regularly-spaced ridges with projected spacings of $\sim$ 2000 AU, suggestive of a large-scale dynamically-important magnetic field. The northern edge of the Wall breaks into several swept-back fragments of width $\sim$ 1800 AU that resemble Kelvin-Helmholtz instabilities, and the southern part of the Wall also shows complex morphologies including a sinusoidal-like variation with a half-wavelength of 2500 AU. Though the dissociation front must increase the density along the surface of the Wall, it does not resolve into pillars that point back to the ionization sources, as could occur if the front triggered new stars to form. We discovered that MHO 1630, an H$_2$ outflow with no clear driving source in the northern portion of the Wall, consists of a series of bow shocks arrayed in a line. △ Less

Submitted 29 September, 2020; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: Accepted in Astrophysical Journal Letters

Proper Motions and Shock Wave Dynamics in the HH 7-11 Stellar Jet

Authors: P. Hartigan, R. Holcomb, A. Frank

Abstract: We have used the Hubble Space Telescope to acquire new broad-band and narrow-band images of the optical line emission and red continuum associated with the HH 7-11 stellar jet in the NGC 1333 star formation region. Combining the new narrow-band images of H$α$, [O~I] $λ$6300 and [S II] $λ$6716 allows us to measure electron densities and excitations at each point in the outflow with the spatial reso… ▽ More We have used the Hubble Space Telescope to acquire new broad-band and narrow-band images of the optical line emission and red continuum associated with the HH 7-11 stellar jet in the NGC 1333 star formation region. Combining the new narrow-band images of H$α$, [O~I] $λ$6300 and [S II] $λ$6716 allows us to measure electron densities and excitations at each point in the outflow with the spatial resolution of HST, while the I-band image traces out the boundary of the cavity evacuated by the outflow. Comparing these images with those taken $\sim$ 20 years ago yields high precision proper motions for all the HH objects in the outflow. HH 11 is a bullet-like clump, and emerges from the exciting source SVS 13A towards the Earth at 24 degrees to line of sight. In contrast, HH 8 and HH 10 consist of two rings of shocked gas that show no bulk proper motions even though the emitting gas is blueshifted. The HH 8 rings are expanding with time. These shocks mark places where ambient material located along the path of the jet redirects the outflow. HH 7 consists of multiple shells, and emits strongly in H$_2$ in what appears to be a terminal bow shock for the outflow, implying that the jet has yet to fully break out of its nascent cloud core. The jet largely fragments into clumps by the time it reaches HH 7. As in the case of HH 110, deflection from ambient material plays a key role in producing observable shock waves in the HH 7-11 outflow. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: 12 figures, plus two animations; Accepted for publication in the Astrophysical Journal

Young Stars and their Variability with LSST

Authors: Rosaria Bonito, Patrick Hartigan, Laura Venuti, Mario Guarcello, Loredana Prisinzano, Costanza Argiroffi, Sergio Messina, Christopher Johns-Krull, Eric Feigelson, John Stauffer, Teresa Giannini, Simone Antoniucci, Salvo Sciortino, Giusi Micela, Ignazio Pillitteri, Davide Fedele, Linda Podio, Francesco Damiani, Peregrine McGehee, Rachel Street, John Gizis, Germano Sacco, Laura Magrini, Ettore Flaccomio, Salvatore Orlando , et al. (13 additional authors not shown)

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

Submitted 7 December, 2018; originally announced December 2018.

Comments: Contribution in response of the Call for White Papers on LSST Cadence Optimization

Photoevaporation of Molecular Clouds in Regions of Massive Star Formation as Revealed Through H$_2$ and Br$γ$ Emission

Authors: S. G. Carlsten, P. M. Hartigan

Abstract: We examine new and pre-existing wide-field, continuum-corrected, narrowband images in H$_2$ 1-0 S(1) and Br$γ$ of three regions of massive star formation: IC 1396, Cygnus OB2, and Carina. These regions contain a variety of globules, pillars, and sheets, so we can quantify how the spatial profiles of emission lines behave in photodissociation regions (PDRs) that differ in their radiation fields and… ▽ More We examine new and pre-existing wide-field, continuum-corrected, narrowband images in H$_2$ 1-0 S(1) and Br$γ$ of three regions of massive star formation: IC 1396, Cygnus OB2, and Carina. These regions contain a variety of globules, pillars, and sheets, so we can quantify how the spatial profiles of emission lines behave in photodissociation regions (PDRs) that differ in their radiation fields and geometries. We have measured 450 spatial profiles of H$_2$ and Br$γ$ along interfaces between HII regions and PDRs. Br$γ$ traces photoevaporative flows from the PDRs, and this emission declines more rapidly with distance as the radius of curvature of the interface decreases, in agreement with models. As noted previously, H$_2$ emission peaks deeper into the cloud relative to Br$γ$, where the molecular gas absorbs far-UV radiation from nearby O-stars. Although PDRs in IC 1396, Cygnus OB2, and Carina experience orders of magnitude different levels of ionizing flux and have markedly differing geometries, all the PDRs have spatial offsets between Br$γ$ and H$_2$ on the order of $10^{17}$cm. There is a weak negative correlation between the offset size and the intensity of ionizing radiation and a positive correlation with the radius of curvature of the cloud. We can reproduce both the size of the offsets and the dependencies of the offsets on these other variables with simple photoevaporative flow models. Both Br$γ$ and H$_2$ 1-0 S(1) will undoubtedly be targeted in future JWST observations of PDRs, so this work can serve as a guide to interpreting these images. △ Less

Submitted 26 October, 2018; originally announced October 2018.

Comments: 19 pages, 17 figures. Accepted to ApJ

A New Look at T Tauri Star Forbidden Lines: MHD Driven Winds from the Inner Disk

Authors: Min Fang, Ilaria Pascucci, Suzan Edwards, Uma Gorti, Andrea Banzatti, Mario Flock, Patrick Hartigan, Gregory J. Herczeg, Andrea K. Dupree

Abstract: Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($Δv\sim$6\kms) analysis of [S II] $λ$4068, [O I] $λ$5577, and [O I] $λ$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity c… ▽ More Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($Δv\sim$6\kms) analysis of [S II] $λ$4068, [O I] $λ$5577, and [O I] $λ$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity component (HVC) associated with jets, and a low-velocity narrow (LVC-NC) and broad (LVC-BC) components. We confirm previous findings that many LVCs are blueshifted by more than 1.5 kms$^{-1}$ thus most likely trace a slow disk wind. We further show that the profiles of individual components are similar in the three lines. We find that most LVC-BC and NC line ratios are explained by thermally excited gas with temperatures between 5,000$-$10,000 K and electron densities $\sim10^{7}-10^{8}$ cm$^{-3}$. The HVC ratios are better reproduced by shock models with a pre-shock H number density of $\sim10^{6}-10^{7}$ cm$^{-3}$. Using these physical properties, we estimate $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ for the LVC and $\dot{M}_{\rm jet}/\dot{M}_{\rm acc}$ for the HVC. In agreement with previous work, the mass carried out in jets is modest compared to the accretion rate. With the likely assumption that the NC wind height is larger than the BC, the LVC-BC $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ is found to be higher than the LVC-NC. These results suggest that most of the mass loss occurs close to the central star, within a few au, through an MHD driven wind. Depending on the wind height, MHD winds might play a major role in the evolution of the disk mass. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 45 pages, 23 figures, and 7 tables, accepted by ApJ

Simulating Radiative Magnetohydrodynamical Flows with AstroBEAR: Implementation and Applications of Non-equilibrium Cooling

Authors: E. C. Hansen, P. Hartigan, A. Frank. A. Wright, J. C. Raymond

Abstract: Radiative cooling plays a crucial role in the dynamics of many astrophysical flows, and is particularly important in the dense shocked gas within Herbig-Haro (HH) objects and stellar jets. Simulating cooling processes accurately is necessary to compare numerical simulations with existing and planned observations of HH objects, such as those from the Hubble Space Telescope and the James Webb Space… ▽ More Radiative cooling plays a crucial role in the dynamics of many astrophysical flows, and is particularly important in the dense shocked gas within Herbig-Haro (HH) objects and stellar jets. Simulating cooling processes accurately is necessary to compare numerical simulations with existing and planned observations of HH objects, such as those from the Hubble Space Telescope and the James Webb Space Telescope. In this paper we discuss a new, non-equilibrium cooling scheme we have implemented into the 3-D magnetohydrodynamic (MHD) code AstroBEAR. The new cooling function includes ionization, recombination, and excitation of all the important atomic species that cool below 10000 K. We tested the routine by comparing its predictions with those from the well-tested 1-D Cox-Raymond shock code (Raymond 1979). The results show thatAstroBEAR accurately tracks the ionization fraction, temperature, and other MHD variables for all low-velocity (.90 km/s) magnetized radiative shock waves. The new routine allows us to predict synthetic emission maps in all the bright forbidden and permitted lines observed in stellar jets, including Hα, [NII], [OI], and [SII]. We present an example as to how these synthetic maps facilitate a direct comparison with narrowband images of HH objects. △ Less

Submitted 15 November, 2018; v1 submitted 6 September, 2018; originally announced September 2018.

Comments: 8 figures

Science-Driven Optimization of the LSST Observing Strategy

Authors: LSST Science Collaboration, Phil Marshall, Timo Anguita, Federica B. Bianco, Eric C. Bellm, Niel Brandt, Will Clarkson, Andy Connolly, Eric Gawiser, Zeljko Ivezic, Lynne Jones, Michelle Lochner, Michael B. Lund, Ashish Mahabal, David Nidever, Knut Olsen, Stephen Ridgway, Jason Rhodes, Ohad Shemmer, David Trilling, Kathy Vivas, Lucianne Walkowicz, Beth Willman, Peter Yoachim, Scott Anderson , et al. (80 additional authors not shown)

Abstract: The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white… ▽ More The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: 312 pages, 90 figures. Browse the current version at https://github.com/LSSTScienceCollaborations/ObservingStrategy, new contributions welcome!

Finding the UV-Visible Path Forward: Proceedings of the Community Workshop to Plan the Future of UV/Visible Space Astrophysics

Authors: Paul A. Scowen, Todd Tripp, Matt Beasley, David Ardila, B-G Andersson, Jesús Maíz Apellániz, Martin Barstow, Luciana Bianchi, Daniela Calzetti, Mark Clampin, Christopher J. Evans, Kevin France, Miriam García García, Ana Gomez de Castro, Walt Harris, Patrick Hartigan, J. Christopher Howk, John Hutchings, Juan Larruquert, Charles F. Lillie, Gary Matthews, Stephan McCandliss, Ron Polidan, Mario R. Perez, Marc Rafelski , et al. (8 additional authors not shown)

Abstract: We present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, a… ▽ More We present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, and the range of possible alternative launch approaches that could also enable some of the science. The recommendations to the Cosmic Origins Program Analysis Group from the workshop attendees on possible future development directions are outlined. △ Less

Submitted 29 November, 2016; originally announced November 2016.

Comments: Proceedings from Workshop held in June 2015 at NASA GSFC on the Future of UV Astronomy from Space

A Candidate Young Massive Planet in Orbit around the Classical T Tauri Star CI Tau

Authors: Christopher M. Johns-Krull, Jacob N. McLane, L. Prato, Christopher J. Crockett, Daniel T. Jaffe, Patrick M. Hartigan, Charles A. Beichman, Naved I. Mahmud, Wei Chen, B. A. Skiff, P. Wilson Cauley, Joshua A. Jones, G. N. Mace

Abstract: The ~2 Myr old classical T Tauri star CI Tau shows periodic variability in its radial velocity (RV) variations measured at infrared (IR) and optical wavelengths. We find that these observations are consistent with a massive planet in a ~9-day period orbit. These results are based on 71 IR RV measurements of this system obtained over 5 years, and on 26 optical RV measurements obtained over 9 years.… ▽ More The ~2 Myr old classical T Tauri star CI Tau shows periodic variability in its radial velocity (RV) variations measured at infrared (IR) and optical wavelengths. We find that these observations are consistent with a massive planet in a ~9-day period orbit. These results are based on 71 IR RV measurements of this system obtained over 5 years, and on 26 optical RV measurements obtained over 9 years. CI Tau was also observed photometrically in the optical on 34 nights over ~one month in 2012. The optical RV data alone are inadequate to identify an orbital period, likely the result of star spot and activity induced noise for this relatively small dataset. The infrared RV measurements reveal significant periodicity at ~9 days. In addition, the full set of optical and IR RV measurements taken together phase coherently and with equal amplitudes to the ~9 day period. Periodic radial velocity signals can in principle be produced by cool spots, hot spots, and reflection of the stellar spectrum off the inner disk, in addition to resulting from a planetary companion. We have considered each of these and find the planet hypothesis most consistent with the data. The radial velocity amplitude yields an Msin(i) of ~8.1 M_Jup; in conjunction with a 1.3 mm continuum emission measurement of the circumstellar disk inclination from the literature, we find a planet mass of ~11.3 M_Jup, assuming alignment of the planetary orbit with the disk. △ Less

Submitted 25 May, 2016; originally announced May 2016.

Comments: 61 pages, 13 figures, accepted for publication in The Astrophysical Journal

When Shock Waves Collide

Authors: P. Hartigan, J. Foster, A. Frank, E. Hansen, K. Yirak, A. S. Liao, P. Graham, B. Wilde, B. Blue, D. Martinez, P. Rosen, D. Farley, R. Paguio

Abstract: Supersonic outflows from objects as varied as stellar jets, massive stars and novae often exhibit multiple shock waves that overlap one another. When the intersection angle between two shock waves exceeds a critical value, the system reconfigures its geometry to create a normal shock known as a Mach stem where the shocks meet. Mach stems are important for interpreting emission-line images of shock… ▽ More Supersonic outflows from objects as varied as stellar jets, massive stars and novae often exhibit multiple shock waves that overlap one another. When the intersection angle between two shock waves exceeds a critical value, the system reconfigures its geometry to create a normal shock known as a Mach stem where the shocks meet. Mach stems are important for interpreting emission-line images of shocked gas because a normal shock produces higher postshock temperatures and therefore a higher-excitation spectrum than an oblique one does. In this paper we summarize the results of a series of numerical simulations and laboratory experiments designed to quantify how Mach stems behave in supersonic plasmas that are the norm in astrophysical flows. The experiments test analytical predictions for critical angles where Mach stems should form, and quantify how Mach stems grow and decay as intersection angles between the incident shock and a surface change. While small Mach stems are destroyed by surface irregularities and subcritical angles, larger ones persist in these situations, and can regrow if the intersection angle changes to become more favorable. The experimental and numerical results show that although Mach stems occur only over a limited range of intersection angles and size scales, within these ranges they are relatively robust, and hence are a viable explanation for variable bright knots observed in HST images at the intersections of some bow shocks in stellar jets. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 2016, accepted by the Astrophysical Journal

The Shock Dynamics of Heterogeneous YSO Jets: 3-D Simulations Meet Multi-Epoch Observations

Authors: E. C. Hansen, A. Frank, P. Hartigan, S. V. Lebedev

Abstract: High resolution observations of Young Stellar Object (YSO) jets show them to be composed of many small-scale knots or clumps. In this paper we report results of 3-D numerical simulations designed to study how such clumps interact and create morphologies and kinematic patterns seen in emission line observations. Our simulations focus on clump scale dynamics by imposing velocity differences between… ▽ More High resolution observations of Young Stellar Object (YSO) jets show them to be composed of many small-scale knots or clumps. In this paper we report results of 3-D numerical simulations designed to study how such clumps interact and create morphologies and kinematic patterns seen in emission line observations. Our simulations focus on clump scale dynamics by imposing velocity differences between spherical, over-dense regions which then lead to the formation of bow shocks as faster clumps overtake slower material. We show that much of the spatial structure apparent in emission line images of jets arises from the dynamics and interactions of these bow shocks. Our simulations show a variety of time-dependent features, including bright knots associated with Mach stems where the shocks intersect, a "frothy" emission structure that arises from the presence of the Non-linear Thin Shell Instability (NTSI) along the surfaces of the bow shocks, and the merging and fragmentation of clumps. Our simulations use a new non-equilibrium cooling method to produce synthetic emission maps in H$α$ and [S II]. These are directly compared with multi-epoch Hubble Space Telescope (HST) observations of Herbig-Haro (HH) jets. We find excellent agreement between features seen in the simulations and the observations in terms of both proper motion and morphologies. Thus we conclude that YSO jets may be dominated by heterogeneous structures and that interactions between these structures and the shocks they produce can account for many details of YSO jet evolution. △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 18 figures

A New Diagnostic of Magnetic Field Strengths in Radiatively-Cooled Shocks

Authors: P. Hartigan, A. Wright

Abstract: We show that it is possible to measure Alfvénic Mach numbers, defined as the shock velocity in the flow divided by the Alfvén velocity, for low-velocity (V$_{shock}$ $\lesssim$ 100 km$\,$s$^{-1}$) radiative shocks. The method combines observations of bright forbidden lines with a measure of the size of the cooling zone, the latter typically obtained from spatial separation between the Balmer emiss… ▽ More We show that it is possible to measure Alfvénic Mach numbers, defined as the shock velocity in the flow divided by the Alfvén velocity, for low-velocity (V$_{shock}$ $\lesssim$ 100 km$\,$s$^{-1}$) radiative shocks. The method combines observations of bright forbidden lines with a measure of the size of the cooling zone, the latter typically obtained from spatial separation between the Balmer emission lines and the forbidden lines. Because magnetic fields become compressed as gas in the postshock region cools, even relatively weak preshock magnetic fields can be detected with this method. We derive analytical formulae that explain how the spatial separations relate to emission-line ratios, and compute a large grid of radiatively-cooled shock models to develop diagnostic diagrams that can be used to derive Alfvénic Mach numbers in flows. Applying the method to existing data for a bright knot in the HH 111 jet, we obtain a relatively low Alfvénic Mach number of $\sim$ 2, indicative of a magnetized jet that has super-magnetosonic velocity perturbations within it. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: Accepted by ApJ

A Survey of Irradiated Pillars, Globules, and Jets in the Carina Nebul

Authors: P. Hartigan, M. Reiter, N. Smith, J. Bally

Abstract: We present wide-field, deep narrowband H$_2$, Br$γ$, H$α$, [S II], [O III], and broadband I and K-band images of the Carina star formation region. The new images provide a large-scale overview of all the H$_2$ and Br$γ$ emission present in over a square degree centered on this signature star forming complex. By comparing these images with archival HST and Spitzer images we observe how intense UV r… ▽ More We present wide-field, deep narrowband H$_2$, Br$γ$, H$α$, [S II], [O III], and broadband I and K-band images of the Carina star formation region. The new images provide a large-scale overview of all the H$_2$ and Br$γ$ emission present in over a square degree centered on this signature star forming complex. By comparing these images with archival HST and Spitzer images we observe how intense UV radiation from O and B stars affects star formation in molecular clouds. We use the images to locate new candidate outflows and identify the principal shock waves and irradiated interfaces within dozens of distinct areas of star-forming activity. Shocked molecular gas in jets traces the parts of the flow that are most shielded from the intense UV radiation. Combining the H$_2$ and optical images gives a more complete view of the jets, which are sometimes only visible in H$_2$. The Carina region hosts several compact young clusters, and the gas within these clusters is affected by radiation from both the cluster stars and the massive stars nearby. The Carina Nebula is ideal for studying the physics of young H II regions and PDR's, as it contains multiple examples of walls and irradiated pillars at various stages of development. Some of the pillars have detached from their host molecular clouds to form proplyds. Fluorescent H$_2$ outlines the interfaces between the ionized and molecular gas, and after removing continuum, we detect spatial offsets between the Br$γ$ and H$_2$ emission along the irradiated interfaces. These spatial offsets can be used to test current models of PDRs once synthetic maps of these lines become available. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: Accepted in the Astronomical Journal

Numerical Simulations of Mach Stem Formation via Intersecting Bow Shocks

Authors: Edward C. Hansen, Adam Frank, Patrick Hartigan

Abstract: Hubble Space Telescope observations show bright knots of H$α$ emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter H$α$ emission since the incoming… ▽ More Hubble Space Telescope observations show bright knots of H$α$ emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter H$α$ emission since the incoming flow to the Mach stem is parallel to the shock normal. In this paper we report first results of a study using 2-D numerical simulations designed to explore Mach stem formation at the intersection of bow shocks formed by hypersonic "bullets" or "clumps". Our 2-D simulations show how the bow shock shapes and intersection angles change as the adiabatic index $γ$ changes. We show that the formation or lack of a Mach stem in our simulations is consistent with the steady-state Mach stem formation theory. Our ultimate goal, which is part of an ongoing research effort, is to characterize the physical and observational consequences of bow shock intersections including the formation of Mach stems. △ Less

Submitted 19 December, 2014; originally announced December 2014.

MHD Effects on Pulsed YSO Jets. I. 2.5-D Simulations

Authors: Edward C. Hansen, Adam Frank, Patrick Hartigan

Abstract: In this paper we explore the dynamics of radiative axisymmetric MHD jets at high resolution using AMR methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end we have implemented microphysics enabling us to produce synthetic maps of H$α$ and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity w… ▽ More In this paper we explore the dynamics of radiative axisymmetric MHD jets at high resolution using AMR methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end we have implemented microphysics enabling us to produce synthetic maps of H$α$ and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity which leads to a complicated structure of internal shocks and rarefactions as has been seen in previous simulations. The high resolution of our simulations allows us to explore in great detail the effect of pinch forces (due to the jet's toroidal magnetic field) within the "working surfaces" where pulses interact. We map the strong H$α$ emission marking shock fronts and the strong [S II] emission inside cooling regions behind shocks as observed with high-resolution images of jets. We find that pinch forces in the stronger field cases produce additional emission regions along the axis as compared with purely hydrodynamic runs. These simulations are a first step to understanding the full 3-D emission properties of radiative MHD jets. △ Less

Submitted 19 December, 2014; originally announced December 2014.

Journal ref: 2015, ApJ, 800, 41

On Sunspot and Starspot Lifetimes

Authors: Stephen J. Bradshaw, Patrick Hartigan

Abstract: We consider the lifetimes of spots on the Sun and other stars from the standpoint of magnetic diffusion. While normal magnetic diffusivity predicts lifetimes of sunspots that are too large by at least two orders of magnitude, turbulent magnetic diffusivity accounts for both the functional form of the solar empirical spot-lifetime relation and for the observed sunspot lifetimes, provided that the r… ▽ More We consider the lifetimes of spots on the Sun and other stars from the standpoint of magnetic diffusion. While normal magnetic diffusivity predicts lifetimes of sunspots that are too large by at least two orders of magnitude, turbulent magnetic diffusivity accounts for both the functional form of the solar empirical spot-lifetime relation and for the observed sunspot lifetimes, provided that the relevant diffusion length is the supergranule size. Applying this relation to other stars, the value of turbulent diffusivity depends almost entirely on supergranule size, with very weak dependence on other variables such as magnetic field strength and density. Overall, the best observational data for other stars is consistent with the extension of the solar relation provided that stellar supergranule sizes for some stars are significantly larger than they are on the Sun. △ Less

Submitted 15 September, 2014; originally announced September 2014.

Jets and Outflows From Star to Cloud: Observations Confront Theory

Authors: A. Frank, T. P. Ray, S. Cabrit, P. Hartigan, H. G. Arce, F. Bacciotti, J. Bally, M. Benisty, J. Eislöffel, M. Güdel, S. Lebedev, B. Nisini, A. Raga

Abstract: In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observat… ▽ More In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows. △ Less

Submitted 14 February, 2014; originally announced February 2014.

Comments: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning

Irradiated Interfaces in the Ara OB1, Carina, Eagle Nebula, and Cyg OB2 Massive Star Formation Regions

Authors: P. Hartigan, J. Palmer, L. I. Cleeves

Abstract: Regions of massive star formation offer some of the best and most easily-observed examples of radiation hydrodynamics. Boundaries where fully-ionized H II regions transition to neutral/molecular photodissociation regions (PDRs) are of particular interest because marked temperature and density contrasts across the boundaries lead to evaporative flows and fluid dynamical instabilities that can evolv… ▽ More Regions of massive star formation offer some of the best and most easily-observed examples of radiation hydrodynamics. Boundaries where fully-ionized H II regions transition to neutral/molecular photodissociation regions (PDRs) are of particular interest because marked temperature and density contrasts across the boundaries lead to evaporative flows and fluid dynamical instabilities that can evolve into spectacular pillar-like structures. When detached from their parent clouds, pillars become ionized globules that often harbor one or more young stars. H2 molecules at the interface between a PDR and an H II region absorb ultraviolet light from massive stars, and the resulting fluoresced infrared emission lines are an ideal way to trace this boundary independent of obscuring dust. This paper presents H2 images of four regions of massive star formation that illustrate different types of PDR boundaries. The Ara OB1 star formation region contains a striking long wall that has several wavy structures which are present in H2, but the emission is not particularly bright because the ambient UV fluxes are relatively low. In contrast, the Carina star formation region shows strong H2 fluorescence both along curved walls and at the edges of spectacular pillars that in some cases have become detached from their parent clouds. The less-spectacular but more well-known Eagle Nebula has two regions that have strong fluorescence in addition to its pillars. While somewhat older than the other regions, Cyg OB2 has the highest number of massive stars of the regions surveyed and contains many isolated, fluoresced globules that have head-tail morphologies which point towards the sources of ionizing radiation. These images provide a collection of potential astrophysical analogs that may relate to ablated interfaces observed in laser experiments of radiation hydrodynamics. △ Less

Submitted 14 April, 2013; originally announced April 2013.

Journal ref: 2012, HEDP 8, 313

A Search for Giant Planet Companions to T Tauri Stars

Authors: Christopher J. Crockett, Naved I. Mahmud, L. Prato, Christopher M. Johns-Krull, Daniel T. Jaffe, Patrick M. Hartigan, Charles A. Beichman

Abstract: We present results from an ongoing multiwavelength radial velocity (RV) survey of the Taurus-Auriga star forming region as part of our effort to identify pre--main sequence giant planet hosts. These 1-3 Myr old T Tauri stars present significant challenges to traditional RV surveys. The presence of strong magnetic fields gives rise to large, cool star spots. These spots introduce significant RV jit… ▽ More We present results from an ongoing multiwavelength radial velocity (RV) survey of the Taurus-Auriga star forming region as part of our effort to identify pre--main sequence giant planet hosts. These 1-3 Myr old T Tauri stars present significant challenges to traditional RV surveys. The presence of strong magnetic fields gives rise to large, cool star spots. These spots introduce significant RV jitter which can mimic the velocity modulation from a planet-mass companion. To distinguish between spot-induced and planet-induced RV modulation, we conduct observations at ~6700 Angstroms and ~2.3 microns and measure the wavelength dependence (if any) in the RV amplitude. CSHELL observations of the known exoplanet host Gl 86 demonstrate our ability to detect not only hot Jupiters in the near infrared but also secular trends from more distant companions. Observations of nine very young stars reveal a typical reduction in RV amplitude at the longer wavelengths by a factor of ~2-3. While we can not confirm the presence of planets in this sample, three targets show different periodicities in the two wavelength regions. This suggests different physical mechanisms underlying the optical and K band variability. △ Less

Submitted 13 November, 2012; v1 submitted 6 November, 2012; originally announced November 2012.

Comments: 18 pages, 20 figures. Accepted to ApJ: Nov 5, 2012

Comparing Poynting flux dominated magnetic tower jets with kinetic-energy dominated jets

Authors: Martín Huarte-Espinosa, Adam Frank, Eric G. Blackman, Andrea Ciardi, Patrick M. Hartigan, Sergey Lebedev, Jeremy P. Chittenden

Abstract: Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental s… ▽ More Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars. △ Less

Submitted 7 September, 2012; originally announced September 2012.

Comments: 11 pages, 4 figures, accepted for publication in the High Energy Density Physics Journal corresponding to the proceedings of the 9th International Conference on High Energy Density Laboratory Astrophysics, May 4, 2012, Tallahassee Florida

On the structure and stability of magnetic tower jets

Authors: Martín Huarte-Espinosa, Adam Frank, Eric G. Blackman, Andrea Ciardi, Patrick M. Hartigan, Sergey Lebedev, Jeremy P. Chittenden

Abstract: Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances fr… ▽ More Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances from the engine distinguishes two classes of models. In magneto-centrifugal launch (MCL) models, magnetic fields dominate only at scales $\lesssim 100$ engine radii, after which the jets become hydrodynamically dominated (HD). By contrast, in Poynting flux dominated (PFD) magnetic tower models, the field dominates even out to much larger scales. To compare the large distance propagation differences of these two paradigms, we perform 3-D ideal MHD AMR simulations of both HD and PFD stellar jets formed via the same energy flux. We also compare how thermal energy losses and rotation of the jet base affects the stability in these jets. For the conditions described, we show that PFD and HD exhibit observationally distinguishable features: PFD jets are lighter, slower, and less stable than HD jets. Unlike HD jets, PFD jets develop current-driven instabilities that are exacerbated as cooling and rotation increase, resulting in jets that are clumpier than those in the HD limit. Our PFD jet simulations also resemble the magnetic towers that have been recently created in laboratory astrophysical jet experiments. △ Less

Submitted 7 September, 2012; v1 submitted 3 April, 2012; originally announced April 2012.

Comments: 16 pages, 11 figures, published in ApJ: ApJ, 757, 66

Starspot-induced optical and infrared radial velocity variability in T Tauri star Hubble 4

Authors: Naved I. Mahmud, Christopher J. Crockett, Christopher M. Johns-Krull, L. Prato, Patrick M. Hartigan, Daniel T. Jaffe, Charles A. Beichman

Abstract: We report optical (6150 Ang) and K-band (2.3 micron) radial velocities obtained over two years for the pre-main sequence weak-lined T Tauri star Hubble I 4. We detect periodic and near-sinusoidal radial velocity variations at both wavelengths, with a semi-amplitude of 1395\pm94 m/s in the optical and 365\pm80 m/s in the infrared. The lower velocity amplitude at the longer wavelength, combined with… ▽ More We report optical (6150 Ang) and K-band (2.3 micron) radial velocities obtained over two years for the pre-main sequence weak-lined T Tauri star Hubble I 4. We detect periodic and near-sinusoidal radial velocity variations at both wavelengths, with a semi-amplitude of 1395\pm94 m/s in the optical and 365\pm80 m/s in the infrared. The lower velocity amplitude at the longer wavelength, combined with bisector analysis and spot modeling, indicates that there are large, cool spots on the stellar surface that are causing the radial velocity modulation. The radial velocities maintain phase coherence over hundreds of days suggesting that the starspots are long-lived. This is one of the first active stars where the spot-induced velocity modulation has been resolved in the infrared. △ Less

Submitted 11 May, 2011; originally announced May 2011.

Comments: Accepted for publication in The Astrophysical Journal

Fluid Dynamics of Stellar Jets in Real Time: Third Epoch HST Images of HH 1, HH 34, and HH 47

Authors: P. Hartigan, A. Frank, J. M. Foster, B. H. Wilde, M. Douglas, P. A. Rosen, R. F. Coker, B. E. Blue, J. F. Hansen

Abstract: We present new, third-epoch HST H-alpha and [S II] images of three HH jets (HH 1 and 2, HH 34, and HH 47) and compare these images with those from the previous epochs. The high-spatial resolution, coupled with a time-series whose cadence is of order both the hydrodynamical and radiative cooling timescales of the flow allows us to follow the hydrodynamical/magnetohydrodynamical evolution of an astr… ▽ More We present new, third-epoch HST H-alpha and [S II] images of three HH jets (HH 1 and 2, HH 34, and HH 47) and compare these images with those from the previous epochs. The high-spatial resolution, coupled with a time-series whose cadence is of order both the hydrodynamical and radiative cooling timescales of the flow allows us to follow the hydrodynamical/magnetohydrodynamical evolution of an astrophysical plasma system in which ionization and radiative cooling play significant roles. Cooling zones behind the shocks are resolved, so it is possible to identify which way material flows through a given shock wave. The images show that heterogeneity is paramount in these jets, with clumps dominating the morphologies of both bow shocks and their Mach disks. This clumpiness exists on scales smaller than the jet widths and determines the behavior of many of the features in the jets. Evidence also exists for considerable shear as jets interact with their surrounding molecular clouds, and in several cases we observe shock waves as they form and fade where material emerges from the source and as it proceeds along the beam of the jet. Fine-structure within two extended bow shocks may result from Mach stems that form at the intersection points of oblique shocks within these clumpy objects. Taken together, these observations represent the most significant foray thus far into the time domain for stellar jets, and comprise one of the richest data sets in existence for comparing the behavior of a complex astrophysical plasma flows with numerical simulations and laboratory experiments. △ Less

Submitted 21 April, 2011; originally announced April 2011.

Comments: 17 figures. High-resolution PDF with better figure quality available at http://sparky.rice.edu/~hartigan/pub.html

Laboratory Experiments, Numerical Simulations, and Astronomical Observations of Deflected Supersonic Jets: Application to HH 110

Authors: P. Hartigan, J. M. Foster, B. H. Wilde, R. F. Coker, P. A. Rosen, J. F. Hansen, B. E. Blue, R. J. R. Williams, R. Carver, A. Frank

Abstract: Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory je… ▽ More Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory jets can be studied for a variety of initial conditions, arbitrary viewing angles, and different times, attributes especially helpful for interpreting astronomical images where the viewing angle and initial conditions are fixed and the time domain is limited. Experiments are also a powerful way to test numerical fluid codes in a parameter range where the codes must perform well. In this paper we combine images from a series of laboratory experiments of deflected supersonic jets with numerical simulations and new spectral observations of an astrophysical example, the young stellar jet HH 110. The experiments provide key insights into how deflected jets evolve in 3-D, particularly within working surfaces where multiple subsonic shells and filaments form, and along the interface where shocked jet material penetrates into and destroys the obstacle along its path. The experiments also underscore the importance of the viewing angle in determining what an observer will see. The simulations match the experiments so well that we can use the simulated velocity maps to compare the dynamics in the experiment with those implied by the astronomical spectra. The experiments support a model where the observed shock structures in HH 110 form as a result of a pulsed driving source rather than from weak shocks that may arise in the supersonic shear layer between the Mach disk and bow shock of the jet's working surface. △ Less

Submitted 1 October, 2009; originally announced October 2009.

Comments: Full resolution figures available at http://sparky.rice.edu/~hartigan/pub.html To appear in ApJ

Journal ref: Astrophys.J.705:1073-1094,2009

Adaptive Optics Echelle Spectroscopy of [Fe II] 1.644 um in the RW Aur Jet: A Narrow Slice Down the Axis of the Flow

Authors: P. Hartigan, L. Hillenbrand

Abstract: We present new adaptive optics echelle spectra of the near-infrared [Fe II] lines in the redshifted and blueshifted jets from the T Tauri star RW Aur. The spectra have an unprecedented combination of high spatial and spectral resolution that makes it possible to trace the dynamics of the flow to a projected distance of only 10 AU from the source. As noted by previous studies, the redshifted flow… ▽ More We present new adaptive optics echelle spectra of the near-infrared [Fe II] lines in the redshifted and blueshifted jets from the T Tauri star RW Aur. The spectra have an unprecedented combination of high spatial and spectral resolution that makes it possible to trace the dynamics of the flow to a projected distance of only 10 AU from the source. As noted by previous studies, the redshifted flow is much slower than its fainter blueshifted counterpart. Our observations clearly show that both the radial velocities and the emission line widths are larger closer to the source on both sides of the jet. The line widths are 20% - 30% of the jet velocity on both sides of the flow, significantly larger than would be produced by a divergent constant velocity flow. The observed line widths could arise from a layered velocity structure in the jet or from magnetic waves. A bright knot in the redshifted jet has no concomitant increase in line width, implying that it is not heated by a bow shock. Alternate heating mechanisms include planar shocks, ambipolar diffusion and magnetic reconnection. △ Less

Submitted 8 September, 2009; originally announced September 2009.

Journal ref: Astrophys.J.705:1388-1394,2009

A Young Planet Search in Visible and IR Light: DN Tau, V836 Tau, and V827 Tau

Authors: L. Prato, M. Huerta, C. M. Johns-Krull, N. Mahmud, D. T. Jaffe, P. Hartigan

Abstract: In searches for low-mass companions to late-type stars, correlation between radial velocity variations and line bisector slope changes indicates contamination by large starspots. Two young stars demonstrate that this test is not sufficient to rule out starspots as a cause of radial velocity variations. As part of our survey for substellar companions to T Tauri stars, we identified the ~2 Myr old… ▽ More In searches for low-mass companions to late-type stars, correlation between radial velocity variations and line bisector slope changes indicates contamination by large starspots. Two young stars demonstrate that this test is not sufficient to rule out starspots as a cause of radial velocity variations. As part of our survey for substellar companions to T Tauri stars, we identified the ~2 Myr old planet host candidates DN Tau and V836 Tau. In both cases, visible light radial velocity modulation appears periodic and is uncorrelated with line bisector span variations, suggesting close companions of several M_Jup in these systems. However, high-resolution, infrared spectroscopy shows that starspots cause the radial velocity variations. We also report unambiguous results for V827 Tau, identified as a spotted star on the basis of both visible light and infrared spectroscopy. Our results suggest that infrared follow up observations are critical for determining the source of radial velocity modulation in young, spotted stars. △ Less

Submitted 21 September, 2008; originally announced September 2008.

Comments: Accepted for publication in the Astrophysical Journal Letters

Star Spot Induced Radial Velocity Variability in LkCa 19

Authors: Marcos Huerta, Christopher M. Johns-Krull, L. Prato, Patrick Hartigan, D. T. Jaffe

Abstract: We describe a new radial velocity survey of T Tauri stars and present the first results. Our search is motivated by an interest in detecting massive young planets, as well as investigating the origin of the brown dwarf desert. As part of this survey, we discovered large-amplitude, periodic, radial velocity variations in the spectrum of the weak line T Tauri star LkCa 19. Using line bisector anal… ▽ More We describe a new radial velocity survey of T Tauri stars and present the first results. Our search is motivated by an interest in detecting massive young planets, as well as investigating the origin of the brown dwarf desert. As part of this survey, we discovered large-amplitude, periodic, radial velocity variations in the spectrum of the weak line T Tauri star LkCa 19. Using line bisector analysis and a new simulation of the effect of star spots on the photometric and radial velocity variability of T Tauri stars, we show that our measured radial velocities for LkCa19 are fully consistent with variations caused by the presence of large star spots on this rapidly rotating young star. These results illustrate the level of activity-induced radial velocity noise associated with at least some very young stars. This activity-induced noise will set lower limits on the mass of a companion detectable around LkCa 19, and similarly active young stars. △ Less

Submitted 15 November, 2007; originally announced November 2007.

Comments: ApJ accepted, 27 pages, 12 figures, aastex

Magnetic Fields in Stellar Jets

Authors: P. Hartigan, A. Frank, P. Varniere, E. Blackman

Abstract: Although several lines of evidence suggest that jets from young stars are driven magnetically from accretion disks, existing observations of field strengths in the bow shocks of these flows imply that magnetic fields play only a minor role in the dynamics at these locations. To investigate this apparent discrepancy we performed numerical simulations of expanding magnetized jets with stochastical… ▽ More Although several lines of evidence suggest that jets from young stars are driven magnetically from accretion disks, existing observations of field strengths in the bow shocks of these flows imply that magnetic fields play only a minor role in the dynamics at these locations. To investigate this apparent discrepancy we performed numerical simulations of expanding magnetized jets with stochastically variable input velocities with the AstroBEAR MHD code. Because the magnetic field B is proportional to the density n within compression and rarefaction regions, the magnetic signal speed drops in rarefactions and increases in the compressed areas of velocity-variable flows. In contrast, B ~ n^0.5 for a steady-state conical flow with a toroidal field, so the Alfven speed in that case is constant along the entire jet. The simulations show that the combined effects of shocks, rarefactions, and divergent flow cause magnetic fields to scale with density as an intermediate power 1 > p > 0.5. Because p > 0.5, the Alfven speed in rarefactions decreases on average as the jet propagates away from the star. This behavior is extremely important to the flow dynamics because it means that a typical Alfven velocity in the jet close to the star is significantly larger than it is in the rarefactions ahead of bow shocks at larger distances, the one place where the field is a measurable quantity. We find that the observed values of weak fields at large distances are consistent with strong fields required to drive the observed mass loss close to the star. For a typical stellar jet the crossover point inside which velocity perturbations of 30 - 40 km/s no longer produce shocks is ~ 300 AU from the source. △ Less

Submitted 6 February, 2007; originally announced February 2007.

Journal ref: Astrophys.J.661:910-918,2007

Collimation, Proper Motions, and Physical Conditions in the HH 30 Jet From HST Slitless Spectroscopy

Authors: P. Hartigan, J. Morse

Abstract: We present Space Telescope Imaging Spectrograph (STIS) spectral images of the HH~30 stellar jet taken through a wide slit over two epochs. The jet is unresolved spectrally, so the observations produce emission-line images for each line in the spectrum. This rich dataset shows how physical conditions in the jet vary with distance and time, produces precise proper motions of knots within the jet,… ▽ More We present Space Telescope Imaging Spectrograph (STIS) spectral images of the HH~30 stellar jet taken through a wide slit over two epochs. The jet is unresolved spectrally, so the observations produce emission-line images for each line in the spectrum. This rich dataset shows how physical conditions in the jet vary with distance and time, produces precise proper motions of knots within the jet, resolves the jet width close to the star, and gives a spectrum of the reflected light from the disk over a large wavelength range at several positions. We introduce a new method for analyzing a set of line ratios based on minimizing a quadratic form between models and data. The method generates images of the density, temperature and ionization fraction computed using all the possible line ratios appropriately weighted. In HH 30, the density declines with distance from the source in a manner consistent with an expanding flow, and is larger by a factor of two along the axis of the jet than it is at the periphery. Ionization in the jet ranges from ~ 5% to 40%, and high ionization/excitation knots form at about 100 AU from the star and propagate outward with the flow. These high-excitation knots are not accompanied by corresponding increases in the density, so if formed by velocity variations the knots must have a strong internal magnetic pressure to smooth out density increases while lengthening recombination times. △ Less

Submitted 19 January, 2007; originally announced January 2007.

Journal ref: Astrophys.J.660:426-440,2007

Infrared [Fe II] emission from P Cygni's nebula: Atomic data, mass, kinematics, and the 1600 A.D. outburst

Authors: Nathan Smith, Patrick Hartigan

Abstract: We present moderate and high-dispersion 1-2.5 micron spectra of the 10 arcsec nebula around P Cygni, dominated by bright emission lines of [Fe II]. Observed [Fe II] line ratios disagree with theoretical transition rates in the literature, so we use the spectrum of P Cygni's nebula to constrain the atomic data for low-lying levels of [Fe II]. Of particular interest is the ratio [Fe II] 12567/1643… ▽ More We present moderate and high-dispersion 1-2.5 micron spectra of the 10 arcsec nebula around P Cygni, dominated by bright emission lines of [Fe II]. Observed [Fe II] line ratios disagree with theoretical transition rates in the literature, so we use the spectrum of P Cygni's nebula to constrain the atomic data for low-lying levels of [Fe II]. Of particular interest is the ratio [Fe II] 12567/16435, often used as a reddening indicator, for which we empirically derive an intrinsic value of 1.49, which is 10--40 per cent higher than previous estimates. High-dispersion spectra of [Fe II] 16435 constrain the geometry, detailed structure, and kinematics of P Cygni's nebula, which is the major product of P Cygni's outburst in 1600 A.D. We use the [N II]/[N I] line ratio to conclude that the nebula is mostly ionized with a total mass of ~0.1 Msun; more than the mass lost by the stellar wind since the eruption. For this mass, we would expect a larger infrared excess than observed. We propose that the dust which obscured the star after the outburst has since been largely destroyed, releasing Fe into the gas phase to produce the bright [Fe II] emission. The kinetic energy of this shell is ~2e46 ergs, far less than the kinetic energy released during the giant eruption of eta Car in the 1840s. The mass and kinetic energy in the nebulae of eta Car and P Cygni give insight to the range of parameters expected for extragalactic eta Car-like eruptions. △ Less

Submitted 30 October, 2005; originally announced October 2005.

Comments: 23 pgs. Accepted by ApJ

Journal ref: Astrophys.J.638:1045-1055,2006

Proper Motions of the HH 47 Jet Observed with the Hubble Space Telescope

Authors: P. Hartigan, S. Heathcote, J. Morse, B. Reipurth, J. Bally

Abstract: We present a proper motion study of the shock waves within the classic stellar jet HH 47 based on Hubble Space Telescope H-alpha and [S II] images of the region taken over two epochs. Individual knots within the jet and in the bow shock/Mach disk working surface of HH 47A move significantly in the five years that separate the images, and the excellent spatial resolution of HST makes it possible… ▽ More We present a proper motion study of the shock waves within the classic stellar jet HH 47 based on Hubble Space Telescope H-alpha and [S II] images of the region taken over two epochs. Individual knots within the jet and in the bow shock/Mach disk working surface of HH 47A move significantly in the five years that separate the images, and the excellent spatial resolution of HST makes it possible to measure the proper motions with enough precision to easily observe differential motions throughout the flow. The bright portion of the jet emerges at 37.5 +/- 2.5 degrees from the plane of the sky with an average velocity of 300 km/s. Dynamical ages of the shock waves in the jet range from a few decades for knots recently ejected by the source to ~ 1300 years for the faint extended bow shock HH 47D. The jet curves, but motions of knots in the jet are directed radially away from the exciting source, and velocity variability in the flow drives the shock waves that heat the jet internally. The jet orientation angle varies with time by about 15 degrees, and currently points to the northwestern portion of a cavity outlined by a reflection nebula, where a quasi-stationary shock deflects the jet. The major working surface HH 47A is more complex than a simple bow shock/Mach disk, and contains numerous clumps that move relative to one another with velocities of ~ +/- 40 km/s. Small clumps or instabilities affect the Mach disk, and dense clumps may move all the way through the working surface to cause the bumpy morphology seen at the bow shock. A localized area between the bow shock and Mach disk varies significantly between the two sets of images. △ Less

Submitted 21 July, 2005; originally announced July 2005.

Comments: full resolution figures available at http://sparky.rice.edu/~hartigan/pub.html; Accepted for publication in the Astronomical Journal

Journal ref: Astron.J.130:2197-2205,2005