Search

We present a study connecting the physical properties of protostellar envelopes to the morphology of the envelope-scale magnetic field. We use the ALMA polarization observations of 55 young prtostars at 0.87 mm on $\sim400-3000$ au scales from the {\em B}-field Orion Protostellar Survey (BOPS) to infer the envelope-scale magnetic field and both dust and gas emission on comparable scales to measure the envelope properties. We find that the protostellar envelopes with compact polarized dust emission tend to have lower envelope masses, than the sources with more extended envelopes. We also find that protostars showing hourglass-field morphologies tend to have lower velocity dispersions in their envelopes, whereas systems with spiral-field morphologies have higher velocity dispersion. Combining with the disk properties taken from the Orion VLA/ALMA Nascent Disk and Multiplicity (VANDAM) survey, we connect envelope properties to fragmentation. Our results suggest that envelope mass may not correlate with fragmentation, whereas turbulence appears to promote fragmentation. On the other hand, we find that fragmentation is suppressed in systems with pinched magnetic fields, suggesting that the magnetic field play a role on providing additional support against gravitational collapse, and the formation of an hourglass-like field may coincide with enhanced magnetic braking that removes angular momentum and hinders the formation of embedded disks. Nevertheless, significant misalignment between magnetic field and outflow axes tends to reduce magnetic braking, leading to the formation of larger disks., Comment: 18 pages, 6 figures, submitted to ApJ

The polarization spectrum, or wavelength dependence of the polarization fraction, of interstellar dust emission provides important insights into the grain alignment mechanism of interstellar dust grains. We investigate the far-infrared polarization spectrum of a realistic simulated high-mass star forming cloud under various models of grain alignment and emission. We find that neither a homogeneous grain alignment model nor a grain alignment model that includes collisional dealignment is able to produce the falling spectrum seen in observations. On the other hand, we find that a grain alignment model with grain alignment efficiency dependent on local temperature is capable of producing a falling spectrum that is in qualitative agreement with observations of OMC-1. For the model most in agreement with OMC-1, we find no correlation between temperature and the slope of the polarization spectrum. However, we do find a positive correlation between column density and the slope of the polarization spectrum. We suggest this latter correlation to be the result of wavelength-dependent polarization by absorption., Comment: 14 pages, 6 figures, 2 tables; Accepted to ApJ

Magnetic fields likely play an important role in the formation of young protostars. Multiscale and multiwavelength dust polarization observations can reveal the inferred magnetic field from scales of the cloud to core to protostar. We present continuum polarization observations of the young protostellar triple system IRAS 16293-2422 at 89 $\mu$m using HAWC+ on SOFIA. The inferred magnetic field is very uniform with an average field angle of 89$^\circ\pm$23$^\circ$ (E of N), which is different from the $\sim$170$^\circ$ field morphology seen at 850 $\mu$m at larger scales (> 2000 au) with JCMT POL-2 and at 1.3 mm on smaller scales (< 300 au) with ALMA. The HAWC+ magnetic field direction is aligned with the known E-W outflow. This alignment difference suggests that the shorter wavelength HAWC+ data is tracing the magnetic field associated with warmer dust likely from the outflow cavity, whereas the longer wavelength data are tracing the bulk magnetic field from cooler dust. Also, we show in this source the dust emission peak is strongly affected by the observing wavelength. The dust continuum peaks closer to source B (northern source) at shorter wavelengths and progressively moves toward the southern A source with increasing wavelength (from 22 $\mu$m to 850 $\mu$m)., Comment: 10 pages, 5 figures. Accepted to ApJ

We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking., Comment: 12 pages, 5 figures (16 pages and 7 figures for the supplementary materials) Accepted for publication on ApJ Letters

Sensitive wide-field observations of polarized thermal emission from interstellar dust grains will allow astronomers to address key outstanding questions about the life cycle of matter and energy driving the formation of stars and the evolution of galaxies. Stratospheric balloon-borne telescopes can map this polarized emission at far-infrared wavelengths near the peak of the dust thermal spectrum - wavelengths that are inaccessible from the ground. In this paper we address the sensitivity achievable by a Super Pressure Balloon (SPB) polarimetry mission, using as an example the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) Observatory. By launching from Wanaka, New Zealand, BLAST Observatory can obtain a 30-day flight with excellent sky coverage - overcoming limitations of past experiments that suffered from short flight duration and/or launch sites with poor coverage of nearby star-forming regions. This proposed polarimetry mission will map large regions of the sky at sub-arcminute resolution, with simultaneous observations at 175, 250, and 350 $\mu m$, using a total of 8274 microwave kinetic inductance detectors. Here, we describe the scientific motivation for the BLAST Observatory, the proposed implementation, and the forecasting methods used to predict its sensitivity. We also compare our forecasted experiment sensitivity with other facilities., Comment: Published in PASP

The Next Generation Deep Extragalactic Exploratory Public (NGDEEP) survey program was designed specifically to include Near Infrared Slitless Spectroscopic observations (NGDEEP-NISS) to detect multiple emission lines in as many galaxies as possible and across a wide redshift range using the Near Infrared Imager and Slitless Spectrograph (NIRISS). We present early results obtained from the the first set of observations (Epoch 1, 50$\%$ of the allocated orbits) of this program (NGDEEP-NISS1). Using a set of independently developed calibration files designed to deal with a complex combination of overlapping spectra, multiple position angles, and multiple cross filters and grisms, in conjunction with a robust and proven algorithm for quantifying contamination from overlapping dispersed spectra, NGDEEP-NISS1 has achieved a 3$\sigma$ sensitivity limit of 2 $\times$ 10$^{-18}$ erg/s/cm$^2$. We demonstrate the power of deep wide field slitless spectroscopy (WFSS) to characterize the star-formation rates, and metallicity ([OIII]/H$\beta$), and dust content, of galaxies at $1

Understanding the formation of substructures in protoplanetary disks is vital for gaining insights into dust growth and the process of planet formation. Studying these substructures in highly embedded Class 0 objects using the Atacama Large Millimeter/submillimeter Array (ALMA), however, poses significant challenges. Nonetheless, it is imperative to do so to unravel the mechanisms and timing behind the formation of these substructures. In this study, we present high-resolution ALMA data at Bands 6 and 4 of the NGC1333 IRAS4A Class 0 protobinary system. This system consists of two components, A1 and A2, separated by 1.8" and located in the Perseus molecular cloud at $\sim$293 pc distance. To gain a comprehensive understanding of the dust properties and formation of substructures in the early stages, we conducted a multi-wavelength analysis of IRAS4A1. Additionally, we sought to address whether the lack of observed substructures in very young disks, could be attributed to factors such as high degrees of disk flaring and large scale heights. To explore this phenomenon, we employed radiative transfer models using RADMC-3D. Our multi-wavelength analysis of A1 discovered characteristics such as high dust surface density, substantial dust mass within the disk, and elevated dust temperatures. These findings suggest the presence of large dust grains compared to the ones in the interstellar medium (ISM), greater than 100 microns in size within the region. Furthermore, while there's no direct detection of any substructure, our models indicate that some, such as a small gap, must be present. In summary, this result implies that disk substructures may be masked or obscured by a large scale height in combination with a high degree of flaring in Class 0 disks. [Abridged], Comment: Accepted for publication in A&A, 13 pages, 7 figures

High-resolution, millimeter observations of disks at the protoplanetary stage reveal substructures such as gaps, rings, arcs, spirals, and cavities. While many protoplanetary disks host such substructures, only a few at the younger protostellar stage have shown similar features. We present a detailed search for early disk substructures in ALMA 1.3 and 0.87~mm observations of ten protostellar disks in the Ophiuchus star-forming region. Of this sample, four disks have identified substructure, two appear to be smooth disks, and four are considered ambiguous. The structured disks have wide Gaussian-like rings ($\sigma_R/R_{\mathrm{disk}}\sim0.26$) with low contrasts ($C<0.2$) above a smooth disk profile, in comparison to protoplanetary disks where rings tend to be narrow and have a wide variety of contrasts ($\sigma_R/R_{\mathrm{disk}}\sim0.08$ and $C$ ranges from $0-1$). The four protostellar disks with the identified substructures are among the brightest sources in the Ophiuchus sample, in agreement with trends observed for protoplanetary disks. These observations indicate that substructures in protostellar disks may be common in brighter disks. The presence of substructures at the earliest stages suggests an early start for dust grain growth and, subsequently, planet formation. The evolution of these protostellar substructures is hypothesized in two potential pathways: (1) the rings are the sites of early planet formation, and the later observed protoplanetary disk ring-gap pairs are secondary features, or (2) the rings evolve over the disk lifetime to become those observed at the protoplanetary disk stage., Comment: Accepted by ApJ, 22 pages, 10 figures

We report the discovery of an accreting supermassive black hole at z=8.679, in CEERS_1019, a galaxy previously discovered via a Ly$\alpha$-break by Hubble and with a Ly$\alpha$ redshift from Keck. As part of the Cosmic Evolution Early Release Science (CEERS) survey, we observed this source with JWST/NIRSpec spectroscopy, MIRI and NIRCam imaging, and NIRCam/WFSS slitless spectroscopy. The NIRSpec spectra uncover many emission lines, and the strong [O III] emission line confirms the ground-based Ly$\alpha$ redshift. We detect a significant broad (FWHM~1200 km/s) component in the H$\beta$ emission line, which we conclude originates in the broad-line region of an active galactic nucleus (AGN), as the lack of a broad component in the forbidden lines rejects an outflow origin. This hypothesis is supported by the presence of high-ionization lines, as well as a spatial point-source component embedded within a smoother surface brightness profile. The mass of the black hole is log($M_{BH}/M_{\odot})=6.95{\pm}0.37$, and we estimate that it is accreting at 1.2 ($\pm$0.5) x the Eddington limit. The 1-8 $\mu$m photometric spectral energy distribution (SED) from NIRCam and MIRI shows a continuum dominated by starlight and constrains the host galaxy to be massive (log M/M$_{\odot}$~9.5) and highly star-forming (SFR~30 M$_{\odot}$ yr$^{-1}$). Ratios of the strong emission lines show that the gas in this galaxy is metal-poor (Z/Z$_{\odot}$~0.1), dense (n$_{e}$~10$^{3}$ cm$^{-3}$), and highly ionized (log U~-2.1), consistent with the general galaxy population observed with JWST at high redshifts. We use this presently highest-redshift AGN discovery to place constraints on black hole seeding models and find that a combination of either super-Eddington accretion from stellar seeds or Eddington accretion from massive black hole seeds is required to form this object by the observed epoch., Comment: 26 pages, 14 figures, 6 tables, published in ApJL

We present the Next Generation Deep Extragalactic Exploratory Public (NGDEEP) Survey, a deep slitless spectroscopic and imaging Cycle 1 JWST treasury survey designed to constrain feedback mechanisms in low-mass galaxies across cosmic time. NGDEEP targets the Hubble Ultra Deep Field (HUDF) with NIRISS slitless spectroscopy (f~1.2e-18 erg/s/cm^2, 5sigma) to measure metallicities and star-formation rates (SFRs) for low-mass galaxies through the peak of the cosmic SFR density (0.512, constraining the slope of the faint-end of the rest-ultraviolet luminosity function. NGDEEP overlaps with the deepest HST ACS optical imaging in the sky: F435W in the HUDF (m=29.6), and F814W in HUDF-Par2 (m=30), making this a premier HST+JWST Deep Field. As a treasury survey, NGDEEP data is public immediately, and we will rapidly release data products and catalogs in the spirit of previous deep field initiatives. In this paper we present the NGDEEP survey design, summarize the science goals, and detail plans for the public release of NGDEEP reduced data products., Comment: 21 pages, 9 figures, submitted to AAS Journals

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

We present a comprehensive analysis of the evolution of the morphological and structural properties of a large sample of galaxies at z=3-9 using early JWST CEERS NIRCam observations. Our sample consists of 850 galaxies at z>3 detected in both CANDELS HST imaging and JWST CEERS NIRCam images to enable a comparison of HST and JWST morphologies. Our team conducted a set of visual classifications, with each galaxy in the sample classified by three different individuals. We also measure quantitative morphologies using the publicly available codes across all seven NIRCam filters. Using these measurements, we present the fraction of galaxies of each morphological type as a function of redshift. Overall, we find that galaxies at z>3 have a wide diversity of morphologies. Galaxies with disks make up a total of 60\% of galaxies at z=3 and this fraction drops to ~30% at z=6-9, while galaxies with spheroids make up ~30-40% across the whole redshift range and pure spheroids with no evidence for disks or irregular features make up ~20%. The fraction of galaxies with irregular features is roughly constant at all redshifts (~40-50%), while those that are purely irregular increases from ~12% to ~20% at z>4.5. We note that these are apparent fractions as many selection effects impact the visibility of morphological features at high redshift. The distributions of S\'ersic index, size, and axis ratios show significant differences between the morphological groups. Spheroid Only galaxies have a higher S\'ersic index, smaller size, and higher axis ratio than Disk/Irregular galaxies. Across all redshifts, smaller spheroid and disk galaxies tend to be rounder. Overall, these trends suggest that galaxies with established disks and spheroids exist across the full redshift range of this study and further work with large samples at higher redshift is needed to quantify when these features first formed., Comment: Accepted for publication to ApJL, 24 pages, 14 figures

Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1

We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. As a result of these observations % that the observations provided, the state-of-the-art theoretical understanding of the formation and evolution of young stellar objects is described. In particular, we emphasize the importance of the coupling, decoupling, and re-coupling between weakly ionized gas and the magnetic field on appropriate scales. This highlights the complex and intimate relationship between gravitational collapse and magnetic fields in young protostars., Comment: 32 pages, 15 figures, accepted as a chapter of Protostars and Planets VII

VLA 1623 West is an ambiguous source that has been described as a shocked cloudlet as well as a protostellar disk. We use deep ALMA 1.3 and 0.87 millimeter observations to constrain its shape and structure to determine its origins better. We use a series of geometric models to fit the uv visibilities at both wavelengths with GALARIO. Although the Real visibilities show structures similar to what has been identified as gaps and rings in protoplanetary disks, we find that a modified Flat-Topped Gaussian at high inclination provides the best fit to the observations. This fit agrees well with expectations for an optically thick, highly inclined disk. Nevertheless, we find that the geometric models consistently yield positive residuals at the four corners of the disk at both wavelengths. We interpret these residuals as evidence that the disk is flared in the millimeter dust. We use a simple toy model for an edge-on flared disk and find that the residuals best match a disk with flaring that is mainly restricted to the outer disk at $R \gtrsim 30$ au. Thus, VLA 1623W may represent a young protostellar disk where the large dust grains have not yet had enough time to settle into the mid-plane. This result may have implications for how disk evolution and vertical dust settling impact the initial conditions leading to planet formation., Comment: Accepted by ApJ, 12 pages, 13 figures

Lyman Break Galaxy (LBG) candidates at z>10 are rapidly being identified in JWST/NIRCam observations. Due to the (redshifted) break produced by neutral hydrogen absorption of rest-frame UV photons, these sources are expected to drop out in the bluer filters while being well detected in redder filters. However, here we show that dust-enshrouded star-forming galaxies at lower redshifts (z<7) may also mimic the near-infrared (near-IR) colors of z>10 LBGs, representing potential contaminants in LBG candidate samples. First, we analyze CEERS-DSFG-1, a NIRCam dropout undetected in the F115W and F150W filters but detected at longer wavelengths. Combining the JWST data with (sub)millimeter constraints, including deep NOEMA interferometric observations, we show that this source is a dusty star-forming galaxy (DSFG) at z~5.1. We also present a tentative 2.6sigma SCUBA-2 detection at 850um around a recently identified z~16 LBG candidate in the same field and show that, if the emission is real and associated with this candidate, the available photometry is consistent with a z~5 dusty galaxy with strong nebular emission lines despite its blue near-IR colors. Further observations on this candidate are imperative to mitigate the low confidence of this tentative submillimeter emission and its positional uncertainty. Our analysis shows that robust (sub)millimeter detections of NIRCam dropout galaxies likely imply z=4-6 redshift solutions, where the observed near-IR break would be the result of a strong rest-frame optical Balmer break combined with high dust attenuation and strong nebular line emission, rather than the rest-frame UV Lyman break. This provides evidence that DSFGs may contaminate searches for ultra high-redshift LBG candidates from JWST observations., Comment: Published in The Astrophysical Journal Letters (updated to match the published version)

We report the discovery of a candidate galaxy with a photo-z of z~12 in the first epoch of the JWST Cosmic Evolution Early Release Science (CEERS) Survey. Following conservative selection criteria we identify a source with a robust z_phot = 11.8^+0.3_-0.2 (1-sigma uncertainty) with m_F200W=27.3, and >7-sigma detections in five filters. The source is not detected at lambda < 1.4um in deep imaging from both HST and JWST, and has faint ~3-sigma detections in JWST F150W and HST F160W, which signal a Ly-alpha break near the red edge of both filters, implying z~12. This object (Maisie's Galaxy) exhibits F115W-F200W > 1.9 mag (2-sigma lower limit) with a blue continuum slope, resulting in 99.6% of the photo-z PDF favoring z > 11. All data quality images show no artifacts at the candidate's position, and independent analyses consistently find a strong preference for z > 11. Its colors are inconsistent with Galactic stars, and it is resolved (r_h = 340 +/- 14 pc). Maisie's Galaxy has log M*/Msol ~ 8.5 and is highly star-forming (log sSFR ~ -8.2 yr^-1), with a blue rest-UV color (beta ~ -2.5) indicating little dust though not extremely low metallicity. While the presence of this source is in tension with most predictions, it agrees with empirical extrapolations assuming UV luminosity functions which smoothly decline with increasing redshift. Should followup spectroscopy validate this redshift, our Universe was already aglow with galaxies less than 400 Myr after the Big Bang., Comment: 17 pages, 6 figures, 3 tables, ApJL in press. Summary of changes from original submission: Improvements in astrometry generated a weak detection in F150W that reduces the photo-z to 11.8 but does not increase the likelihood of lower-z solutions. A full discussion of changes from the original version is available at: https://web.corral.tacc.utexas.edu/ceersdata/papers/Maisie_update.pdf

We present H-band (1.65 $\mu$m) and SOFIA HAWC+ 154 $\mu$m polarization observations of the low-mass core L483. Our H-band observations reveal a magnetic field that is overwhelmingly in the E-W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish $^{12}$CO observations. From our 154 $\mu$m data, we infer a $\sim$ 45$^{\circ}$ twist in the magnetic field within the inner 5" (1000 au) of L483. We compare these new observations with published single-dish 350 $\mu$m polarimetry and find that the 10,000 au scale H-band data match the smaller scale 350 $\mu$m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm ALMA data of L483 which reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the $\sim$ 1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario., Comment: Accepted for publication in ApJ

We present VLA observations toward the Class 0 protostar L1157 MMS at 6.8 mm and 9 mm with a resolution of ~0.04" (14 au). We detect two sources within L1157 MMS and interpret these sources as a binary protostar with a separation of ~16 au. The material directly surrounding the binary system within the inner 50 au radius of the system has an estimated mass of 0.11 M_sun, calculated from the observed dust emission. We interpret the observed binary system in the context of previous observations of its flattened envelope structure, low rates of envelope rotation from 5000 to 200 au scales, and an ordered, poloidal magnetic field aligned with the outflow. Thus, L1157 MMS is a prototype system for magnetically-regulated collapse and the presence of a compact binary within L1157 MMS demonstrates that multiple star formation can still occur within envelopes that likely have dynamically important magnetic fields., Comment: 21 pages, 6 figures, 4 tables, accepted to ApJ

Polarized dust continuum emission has been observed with ALMA in an increasing number of deeply embedded protostellar systems. It generally shows a sharp transition going from the protostellar envelope to the disk scale, with the polarization fraction typically dropping from ${\sim} 5\%$ to ${\sim} 1\%$ and the inferred magnetic field orientations becoming more aligned with the major axis of the system. We quantitatively investigate these observational trends using a sample of protostars in the Perseus molecular cloud and compare these features with a non-ideal MHD disk formation simulation. We find that the gas density increases faster than the magnetic field strength in the transition from the envelope to the disk scale, which makes it more difficult to magnetically align the grains on the disk scale. Specifically, to produce the observed ${\sim} 1\%$ polarization at ${\sim} 100\,\mathrm{au}$ scale via grains aligned with the B-field, even relatively small grains of $1\,\mathrm{\mu m}$ in size need to have their magnetic susceptibilities significantly enhanced (by a factor of ${\sim} 20$) over the standard value, potentially through superparamagnetic inclusions. This requirement is more stringent for larger grains, with the enhancement factor increasing linearly with the grain size, reaching ${\sim} 2\times 10^4$ for millimeter-sized grains. Even if the required enhancement can be achieved, the resulting inferred magnetic field orientation in the simulation does not show a preference for the major axis, which is inconsistent with the observed pattern. We thus conclude that the observed trends are best described by the model where the polarization on the envelope scale is dominated by magnetically aligned grains and that on the disk scale by scattering., Comment: 13 pages, 6 figures. To be published in MNRAS

We present a study of the relative orientation between the magnetic field and elongated cloud structures for the $\rho$ Oph A and $\rho$ Oph E regions in L1688 in the Ophiuchus molecular cloud. Combining inferred magnetic field orientation from HAWC+ 154 $\mu$m observations of polarized thermal emission with column density maps created using Herschel submillimeter observations, we find consistent perpendicular relative alignment at scales of $0.02$ pc ($33.6"$ at $d \approx 137$ pc) using the histogram of relative orientations (HRO) technique. This supports the conclusions of previous work using Planck polarimetry and extends the results to higher column densities. Combining this HAWC+ HRO analysis with a new Planck HRO analysis of L1688, the transition from parallel to perpendicular alignment in L1688 is observed to occur at a molecular hydrogen column density of approximately $10^{21.7}$ cm$^{-2}$. This value for the alignment transition column density agrees well with values found for nearby clouds via previous studies using only Planck observations. Using existing turbulent, magnetohydrodynamic simulations of molecular clouds formed by colliding flows as a model for L1688, we conclude that the molecular hydrogen volume density associated with this transition is approximately $\sim10^{4}$ cm$^{-3}$. We discuss the limitations of our analysis, including incomplete sampling of the dense regions in L1688 by HAWC+., Comment: To be published in ApJ

The BLAST Observatory is a proposed superpressure balloon-borne polarimeter designed for a future ultra-long duration balloon campaign from Wanaka, New Zealand. To maximize scientific output while staying within the stringent superpressure weight envelope, BLAST will feature new 1.8m off-axis optical system contained within a lightweight monocoque structure gondola. The payload will incorporate a 300L $^4$He cryogenic receiver which will cool 8,274 microwave kinetic inductance detectors (MKIDs) to 100mK through the use of an adiabatic demagnetization refrigerator (ADR) in combination with a $^3$He sorption refrigerator all backed by a liquid helium pumped pot operating at 2K. The detector readout utilizes a new Xilinx RFSOC-based system which will run the next-generation of the BLAST-TNG KIDPy software. With this instrument we aim to answer outstanding questions about dust dynamics as well as provide community access to the polarized submillimeter sky made possible by high-altitude observing unrestricted by atmospheric transmission. The BLAST Observatory is designed for a minimum 31-day flight of which 70$\%$ will be dedicated to observations for BLAST scientific goals and the remaining 30$\%$ will be open to proposals from the wider astronomical community through a shared-risk proposals program., Comment: Presented at SPIE Ground-based and Airborne Telescopes VIII, December 13-18, 2020

The Next Generation Balloon-borne Large Aperture Submillimeter Telescope (BLAST-TNG) is a submillimeter polarimeter designed to map interstellar dust and galactic foregrounds at 250, 350, and 500 microns during a 24-day Antarctic flight. The BLAST-TNG detector arrays are comprised of 918, 469, and 272 MKID pixels, respectively. The pixels are formed from two orthogonally oriented, crossed, linear-polarization sensitive MKID antennae. The arrays are cooled to sub 300mK temperatures and stabilized via a closed cycle $^3$He sorption fridge in combination with a $^4$He vacuum pot. The detectors are read out through a combination of the second-generation Reconfigurable Open Architecture Computing Hardware (ROACH2) and custom RF electronics designed for BLAST-TNG. The firmware and software designed to readout and characterize these detectors was built from scratch by the BLAST team around these detectors, and has been adapted for use by other MKID instruments such as TolTEC and OLIMPO. We present an overview of these systems as well as in-depth methodology of the ground-based characterization and the measured in-flight performance., Comment: Presented at SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, December 13-18, 2020

The Next Generation Balloon-Borne Large Aperture Submillimeter Telescope (BLAST-TNG) was a unique instrument for characterizing the polarized submillimeter sky at high-angular resolution. BLAST-TNG flew from the Long Duration Balloon Facility in Antarctica in January 2020. Despite the short flight duration, the instrument worked very well and is providing significant information about each subsystem that will be invaluable for future balloon missions. In this contribution, we discuss the performance of telescope and gondola., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation, Ground-based and Airborne Telescopes VIII

Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $\Delta V <5000$ km s$^{-1}$) at $0.5

We analyze several Galerkin approximations of a Gaussian random field $\mathcal{Z}\colon\mathcal{D}\times\Omega\to\mathbb{R}$ indexed by a Euclidean domain $\mathcal{D}\subset\mathbb{R}^d$ whose covariance structure is determined by a negative fractional power $L^{-2\beta}$ of a second-order elliptic differential operator $L:= -\nabla\cdot(A\nabla) + \kappa^2$. Under minimal assumptions on the domain $\mathcal{D}$, the coefficients $A\colon\mathcal{D}\to\mathbb{R}^{d\times d}$, $\kappa\colon\mathcal{D}\to\mathbb{R}$, and the fractional exponent $\beta>0$, we prove convergence in $L_q(\Omega; H^\sigma(\mathcal{D}))$ and in $L_q(\Omega; C^\delta(\overline{\mathcal{D}}))$ at (essentially) optimal rates for (i) spectral Galerkin methods and (ii) finite element approximations. Specifically, our analysis is solely based on $H^{1+\alpha}(\mathcal{D})$-regularity of the differential operator $L$, where $0<\alpha\leq 1$. For this setting, we furthermore provide rigorous estimates for the error in the covariance function of these approximations in $L_{\infty}(\mathcal{D}\times\mathcal{D})$ and in the mixed Sobolev space $H^{\sigma,\sigma}(\mathcal{D}\times\mathcal{D})$, showing convergence which is more than twice as fast compared to the corresponding $L_q(\Omega; H^\sigma(\mathcal{D}))$-rate. For the well-known example of such Gaussian random fields, the original Whittle-Mat\'ern class, where $L=-\Delta + \kappa^2$ and $\kappa \equiv \operatorname{const.}$, we perform several numerical experiments which validate our theoretical results., Comment: 41 pages, 2 figures

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions., Comment: 28 pages, 14 figures, ApJ, accepted

We present high-sensitivity ($\sigma_I\sim0.2-0.5$ mJy, $\sigma_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $\mu$m dust emission from the young multiple system VLA 1623 in $\rho$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W, and L1527 strongly in the polarized emission, at the $\sim1-3\%$ level. We spatially resolve VLA 1623A into two sources, VLA 1623Aa and VLA 1623Ab, separated by $\sim30$ au and located within a cavity of radius $\sim50$ au within the circumbinary Keplerian disk, as well as the edge-on disk of VLA 1623W. The polarization angle of the emission is uniform across each protostellar source and nearly coincides with each disk's minor axis. The offsets between the minor axis and the polarization angle are not uniformly distributed at the $P\lesssim2\times10^{-4}$ level. The circumbinary disk surrounding VLA 1623Aab is azimuthally symmetrically polarized. Each compact source's emission is partially optically thick ($\tau\gtrsim1$) at 872 $\mu$m, complicating interpretations of polarization involving aligned grains. We find evidence against alignment by radiative flux in each source, particularly in the edge-on VLA 1623W and L1527. We detect astrometric offsets between the polarized emission and the total intensity in VLA 1623Aa, VLA 1623Ab, and VLA 1623B, as predicted if self-scattering in the optically thick limit operates. We conclude that self-scattering is likely responsible for disk-scale polarization at 872 $\mu$m in these systems., Comment: 22 pages, 7 figures, accepted to Astrophysical Journal

In this paper we develop a stochastic integration theory for processes with values in a quasi-Banach space. The integrator is a cylindrical Brownian motion. The main results give sufficient conditions for stochastic integrability. They are natural extensions of known results in the Banach space setting. We apply our main results to the stochastic heat equation where the forcing terms are assumed to have Besov regularity in the space variable with integrability exponent $p\in (0,1]$. The latter is natural to consider for its potential application to adaptive wavelet methods for stochastic partial differential equations.

We analyze the dust morphology of 29 transition disks (TDs) observed with ALMA at (sub-) millimeter-emission. We perform the analysis in the visibility plane to characterize the total flux, cavity size, and shape of the ring-like structure. First, we found that the $M_{\rm{dust}}-M_\star$ relation is much flatter for TDs than the observed trends from samples of class II sources in different star forming regions. This relation demonstrates that cavities open in high (dust) mass disks, independent of the stellar mass. The flatness of this relation contradicts the idea that TDs are a more evolved set of disks. Two potential reasons (not mutually exclusive) may explain this flat relation: the emission is optically thick or/and millimeter-sized particles are trapped in a pressure bump. Second, we discuss our results of the cavity size and ring width in the context of different physical processes for cavity formation. Photoevaporation is an unlikely leading mechanism for the origin of the cavity of any of the targets in the sample. Embedded giant planets or dead zones remain as potential explanations. Although both models predict correlations between the cavity size and the ring shape for different stellar and disk properties, we demonstrate that with the current resolution of the observations, it is difficult to obtain these correlations. Future observations with higher angular resolution observations of TDs with ALMA will help to discern between different potential origins of cavities in TDs., Comment: Minor changes after language edition. Accepted for publication in ApJ