Showing 1–50 of 153 results for author: Pérez, L

A Dust-Trapping Ring in the Planet-Hosting Disk of Elias 2-24

Authors: Adolfo S. Carvalho, Laura M. Perez, Anibal Sierra, Maria Jesus Mellado, Lynne A. Hillenbrand, Sean Andrews, Myriam Benisty, Tilman Birnstiel, John M. Carpenter, Viviana V. Guzman, Jane Huang, Andrea Isella, Nicolas Kurtovic, Luca Ricci, David J. Wilner

Abstract: Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid gr… ▽ More Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid growth and possibly planetesimal production via the streaming instability. We present evidence of dust trapping in the bright ring of the planet-hosting disk Elias 2-24, from the analysis of 1.3 mm and 3 mm ALMA observations at high spatial resolution (0.029 arcsec, 4.0 au). We leverage the high spatial resolution to demonstrate that larger grains are more efficiently trapped and place constraints on the local turbulence ($8 \times 10^{-4} < α_\mathrm{turb} < 0.03$) and the gas-to-dust ratio ($Σ_g / Σ_d < 30$) in the ring. Using a scattering-included marginal probability analysis we measure a total dust disk mass of $M_\mathrm{dust} = 13.8^{+0.7}_{-0.5} \times 10^{-4} \ M_\odot$. We also show that at the orbital radius of the proposed perturber, the gap is cleared of material down to a flux contrast of 10$^{-3}$ of the peak flux in the disk. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, accepted to ApJ

LtU-ILI: An All-in-One Framework for Implicit Inference in Astrophysics and Cosmology

Authors: Matthew Ho, Deaglan J. Bartlett, Nicolas Chartier, Carolina Cuesta-Lazaro, Simon Ding, Axel Lapel, Pablo Lemos, Christopher C. Lovell, T. Lucas Makinen, Chirag Modi, Viraj Pandya, Shivam Pandey, Lucia A. Perez, Benjamin Wandelt, Greg L. Bryan

Abstract: This paper presents the Learning the Universe Implicit Likelihood Inference (LtU-ILI) pipeline, a codebase for rapid, user-friendly, and cutting-edge machine learning (ML) inference in astrophysics and cosmology. The pipeline includes software for implementing various neural architectures, training schemata, priors, and density estimators in a manner easily adaptable to any research workflow. It i… ▽ More This paper presents the Learning the Universe Implicit Likelihood Inference (LtU-ILI) pipeline, a codebase for rapid, user-friendly, and cutting-edge machine learning (ML) inference in astrophysics and cosmology. The pipeline includes software for implementing various neural architectures, training schemata, priors, and density estimators in a manner easily adaptable to any research workflow. It includes comprehensive validation metrics to assess posterior estimate coverage, enhancing the reliability of inferred results. Additionally, the pipeline is easily parallelizable and is designed for efficient exploration of modeling hyperparameters. To demonstrate its capabilities, we present real applications across a range of astrophysics and cosmology problems, such as: estimating galaxy cluster masses from X-ray photometry; inferring cosmology from matter power spectra and halo point clouds; characterizing progenitors in gravitational wave signals; capturing physical dust parameters from galaxy colors and luminosities; and establishing properties of semi-analytic models of galaxy formation. We also include exhaustive benchmarking and comparisons of all implemented methods as well as discussions about the challenges and pitfalls of ML inference in astronomical sciences. All code and examples are made publicly available at https://github.com/maho3/ltu-ili. △ Less

Submitted 2 July, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: 22 pages, 10 figures, accepted in the Open Journal of Astrophysics. Code available at https://github.com/maho3/ltu-ili

Journal ref: 2024 OJA, Vol. 7

High-resolution ALMA observations of compact discs in the wide-binary system Sz 65 and Sz 66

Authors: J. M. Miley, J. Carpenter, R. Booth, J. Jennings, T. J. Haworth, M. Vioque, S. Andrews, D. Wilner, M. Benisty, J. Huang, L. Perez, V. Guzman, L. Ricci, A. Isella

Abstract: Substructures in disc density are ubiquitous in the bright extended discs that are observed with high resolution. These substructures are intimately linked to the physical mechanisms driving planet formation and disc evolution. Surveys of star-forming regions find that most discs are in fact compact, less luminous, and do not exhibit these same substructures. It remains unclear whether compact dis… ▽ More Substructures in disc density are ubiquitous in the bright extended discs that are observed with high resolution. These substructures are intimately linked to the physical mechanisms driving planet formation and disc evolution. Surveys of star-forming regions find that most discs are in fact compact, less luminous, and do not exhibit these same substructures. It remains unclear whether compact discs also have similar substructures or if they are featureless. This suggests that different planet formation and disc evolution mechanisms operate in these discs. We investigated evidence of substructure within two compact discs around the stars Sz 65 and Sz 66 using high angular resolution observations with ALMA at 1.3 mm. The two stars form a wide-binary system with 6.36 arcsec separation. The continuum observations achieve a synthesised beam major axis of 0.026 arcsec, equivalent to about 4.0 au, enabling a search for substructure on these spatial scales and a characterisation of the gas and dust disc sizes with high precision. We analysed the data in the image plane through an analysis of reconstructed images, as well as in the uv plane by modelling the visibilities and by an analysis of the 12CO emission line. Comparisons were made with high-resolution observations of compact discs and radially extended discs. We find evidence of substructure in the dust distribution of Sz 65, namely a shallow gap centred at approximately 20 au, with an emission ring exterior to it. Ninety percent of the measured continuum flux is found within 27 au, and the distance for 12CO is 142 au. The observations show that Sz 66 is very compact: 90 per cent of the continuum flux is contained within 16 au, and 48 au for the gas. While the overall prevalence and diversity of substructure in compact discs relative to larger discs is yet to be determined, we find evidence that substructures can exist in compact discs. △ Less

Submitted 13 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Journal ref: A&A, 682, A55 (2024)

Magnetized Strange Stars and Signals of Gravitational Waves

Authors: S. López Pérez, D. Manreza Paret, A. Pérez Martínez

Abstract: We study the emission of gravitational waves from spheroidal magnetized strange stars for both an isolated slowly rotating star and a binary system. In the first case, we compute the quadrupole moment and the amplitude of gravitational waves that may be emitted. For the binary system, the tidal deformability is obtained by solving simultaneously the system of spheroidal structure equations and the… ▽ More We study the emission of gravitational waves from spheroidal magnetized strange stars for both an isolated slowly rotating star and a binary system. In the first case, we compute the quadrupole moment and the amplitude of gravitational waves that may be emitted. For the binary system, the tidal deformability is obtained by solving simultaneously the system of spheroidal structure equations and the Love number equation. These results are compared with the data inferred from the GW170817 event which is also used to calculate the mass and tidal deformability of the companion star in the binary system. Our model supports binary systems formed by magnetized strange stars describing reasonable signals of gravitational waves contrasted with other models of binary systems composed of magnetized hadronic stars and non-magnetized quark stars. △ Less

Submitted 31 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland , et al. (229 additional authors not shown)

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

Field-level simulation-based inference with galaxy catalogs: the impact of systematic effects

Authors: Natalí S. M. de Santi, Francisco Villaescusa-Navarro, L. Raul Abramo, Helen Shao, Lucia A. Perez, Tiago Castro, Yueying Ni, Christopher C. Lovell, Elena Hernandez-Martinez, Federico Marinacci, David N. Spergel, Klaus Dolag, Lars Hernquist, Mark Vogelsberger

Abstract: It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocit… ▽ More It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocities of galaxies that are robust to uncertainties in astrophysics and subgrid models. However, observations are affected by many effects, including 1) masking, 2) uncertainties in peculiar velocities and radial distances, and 3) different galaxy selections. Moreover, observations only allow us to measure redshift, intertwining galaxies' radial positions and velocities. In this paper we train and test our models on galaxy catalogs, created from thousands of state-of-the-art hydrodynamic simulations run with different codes from the CAMELS project, that incorporate these observational effects. We find that, although the presence of these effects degrades the precision and accuracy of the models, and increases the fraction of catalogs where the model breaks down, the fraction of galaxy catalogs where the model performs well is over 90 %, demonstrating the potential of these models to constrain cosmological parameters even when applied to real data. △ Less

Submitted 9 May, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 39 pages, 25 figures. For the reference in the abstract (de Santi et al. 2023) see arXiv:2302.14101

A magnetically driven disc wind in the inner disc of PDS 70

Authors: Justyn Campbell-White, Carlo F. Manara, Myriam Benisty, Antonella Natta, Rik A. B. Claes, Antonio Frasca, Jaehan Bae, Stefano Facchini, Andrea Isella, Laura Pérez, Paola Pinilla, Aurora Sicilia-Aguilar, Richard Teague

Abstract: PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant r… ▽ More PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O I] 6300AA line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O I] line luminosity is high with respect to the accretion luminosity when compared to a large sample of discs with cavities in nearby star-forming regions. The FWHM and blue-shifted peak of the [O I] line suggest an emission in a region very close to the star, favouring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He I 10830AA line, which is unusual for low-accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner disc wind could also have had a significant contribution to clearing the inner-disc. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Constraints on the Epoch of Reionization with Roman Space Telescope and the Void Probability Function of Lyman-Alpha Emitters

Authors: Lucia A. Perez, Sangeeta Malhotra, James E. Rhoads, Isak G. B. Wold

Abstract: We use large simulations of Lyman-Alpha Emitters with different fractions of ionized intergalactic medium to quantify the clustering of Ly$α$ emitters as measured by the Void Probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace R… ▽ More We use large simulations of Lyman-Alpha Emitters with different fractions of ionized intergalactic medium to quantify the clustering of Ly$α$ emitters as measured by the Void Probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace Roman Space Telescope. Since Roman will be able to carry out blind spectroscopic surveys of Ly$α$ emitters continuously between $7<z<12$ to sensitivities of at least $10^{-17}$ erg sec$^{-1}$ over a wide field of view, it can measure the epoch of reionization as well as the pace of ionization of the intergalactic medium (IGM). We compare deep Roman surveys covering roughly 1, 4, and 16 deg$^2$, and quantify what constraints on reionization the VPF may find for these surveys. A survey of 1 deg$^2$ would distinguish between very late reionization and early reionization to 3$σ$ near $z=7.7$ with the VPF. The VPF of a 4 deg$^2$ survey can distinguish between slow vs.\ fast, and early vs.\ late, reionization at $> 3-4σ$ at several redshifts between $7<z<9$. However, a survey of 13-16 deg$^2$ would allow the VPF to give several robust constraints ($>5-8σ$) across the epoch of reionization, and would yield a detailed history of the reionization of the IGM and its effect on Lyman-$α$ Emitter clustering. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 19 total pages, 7 total figures, appendix with 3 additional figures

Kinematic and thermal signatures of the directly imaged protoplanet candidate around Elias 2-24

Authors: Christophe Pinte, Iain Hammond, Daniel J. Price, Valentin Christiaens, Sean M. Andrews, Gaël Chauvin, Laura M. Pérez, Sebastián Jorquera, Himanshi Garg, Brodie J. Norfolk, Josh Calcino, Mickaël Bonnefoy

Abstract: We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require loc… ▽ More We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models show this could either be a $\gtrsim 5$M$_\mathrm{Jup}$, or a lower mass ( $\gtrsim 2$M$_\mathrm{Jup}$) but accreting proto-planet. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

A kinematically detected planet candidate in a transition disk

Authors: Jochen Stadler, Myriam Benisty, Andrés F. Izquierdo, Stefano Facchini, Richard Teague, Nicolas Kurtovic, Paola Pinilla, Jaehan Bae, Megan Ansdell, Ryan Loomis, Satoshi Mayama, Laura M. Pérez, Leonardo Testi

Abstract: Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps an… ▽ More Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps and moment maps of the line of sight velocity and peak intensity. We fitted a Keplerian model of the channel-by-channel emission to study line profile differences and produced deprojected radial profiles for all velocity components. The $^{12}$CO emission shows a cavity inwards of $\sim$56 au and within the dust continuum ring at 81 au. The azimuthal brightness variations in the $^{12}$CO line and dust continuum are broadly aligned with the shadows detected in scattered-light observations. We find a strong localized non-Keplerian feature toward the west within the continuum ring (at $R=41\pm 10$ au and $PA=280\pm 2 ^\circ$). A tightly wound spiral is also detected which extends over 300$^\circ$ in azimuth, possibly connected to the localized non-Keplerian feature. Finally, a bending of the iso-velocity contours within the gas cavity indicates a highly perturbed inner region, possibly related to the presence of a misaligned inner disk. While broadly aligned with the scattered-light shadows, the localized non-Keplerian feature cannot be solely due to changes in temperature. Instead, we interpret the kinematical feature as tracing a massive companion located at the edge of the dust continuum ring. We speculate that the spiral is caused by buoyancy resonances driven by planet-disk interactions. However, this potential planet at $\sim$41 au cannot explain the gas-depleted cavity, the low accretion rate, and the misaligned inner disk, which suggests the presence of another companion closer in. △ Less

Submitted 11 January, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 13 pages, 14 figures. To be published in Astronomy & Astrophysics Letters

Journal ref: A&A 670, L1 (2023)

Mapping Protoplanetary Disk Vertical Structure with CO Isotopologue Line Emission

Authors: Charles J. Law, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Myriam Benisty, Stefano Facchini, Kevin Flaherty, Andrea Isella, Sheng Jin, Jun Hashimoto, Jane Huang, Ryan A. Loomis, Feng Long, Carlos E. Muñoz-Romero, Teresa Paneque-Carreño, Laura M. Pérez, Chunhua Qi, Kamber R. Schwarz, Jochen Stadler, Takashi Tsukagoshi, David J. Wilner, Gerrit van der Plas

Abstract: High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we… ▽ More High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the 2D (r, z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high angular resolution (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=3-2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=3-2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that $^{13}$CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 31 pages, 15 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/7430257

Probing Patchy Reionization with the Void Probability Function of Lyman-$α$ Emitters

Authors: Lucia A. Perez, Sangeeta Malhotra, James E. Rhoads, Peter Laursen, Isak G. B. Wold

Abstract: We probe what constraints for the global ionized hydrogen fraction the Void Probability Function (VPF) clustering can give for the Lyman-Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyman-alpha emission, and measuring the drop in the luminosity function of Lyman-$α$ emitters (LAEs) has been used to constrain the i… ▽ More We probe what constraints for the global ionized hydrogen fraction the Void Probability Function (VPF) clustering can give for the Lyman-Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyman-alpha emission, and measuring the drop in the luminosity function of Lyman-$α$ emitters (LAEs) has been used to constrain the ionization fraction in narrowband surveys. However, the clustering of LAEs is independent from the luminosity function's inherent evolution, and can offer additional constraints for reionization under different models. The VPF measures how likely a given circle is to be empty. It is a volume-averaged clustering statistic that traces the behavior of higher order correlations, and its simplicity offers helpful frameworks for planning surveys. Using the \citet{Jensen2014} simulations of LAEs within various amount of ionized intergalactic medium, we predict the behavior of the VPF in one (301x150.5x30 Mpc$^3$), four (5.44$\times 10^6$ Mpc$^3$), or eight (1.1$\times 10^7$ Mpc$^3$) fields of LAGER imaging. We examine the VPF at 5 and 13 arcminutes, corresponding to the minimum scale implied by the LAE density and the separation of the 2D VPF from random, and the maximum scale from the 8-field 15.5 deg$^2$ LAGER area. We find that even a single DECam field of LAGER (2-3 deg$^2$) could discriminate between mostly neutral vs. ionized. Additionally, we find four fields allows the distinction between 30, 50, and 95 percent ionized; and that eight fields could even distinguish between 30, 50, 73, and 95 percent ionized. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: 14 pages, 5 figures

Vertically extended and asymmetric CN emission in the Elias 2-27 protoplanetary disk

Authors: T. Paneque-Carreño, A. Miotello, E. F. van Dishoeck, L. M. Pérez, S. Facchini, A. Izquierdo, L. Tychoniec, L. Testi

Abstract: Elias 2-27 is a young star that hosts an extended, bright and inclined disk of dust and gas. The inclination and extreme flaring of the disk make Elias 2-27 an ideal target to study the vertical distribution of molecules, particularly CN. We directly trace the emission of CN in Elias 2-27 and compare it to previously published CO isotopologue data. CN $N = 3-2$ emission is analyzed in two differen… ▽ More Elias 2-27 is a young star that hosts an extended, bright and inclined disk of dust and gas. The inclination and extreme flaring of the disk make Elias 2-27 an ideal target to study the vertical distribution of molecules, particularly CN. We directly trace the emission of CN in Elias 2-27 and compare it to previously published CO isotopologue data. CN $N = 3-2$ emission is analyzed in two different transitions $J = 7/2 - 5/2$ and $J = 5/2 - 3/2$, for which we detect two hyperfine group transitions. The vertical location of CN emission is traced directly from the channel maps, following geometrical methods that have been previously used to analyze the CO emission of Elias 2-27. Analytical models are used to parametrize the vertical profile of each molecule and study the extent of each tracer, additionally we compute radial profiles of column density and optical depth. We show that the vertical location of CN and CO isotopologues in Elias 2-27 is layered and consistent with predictions from thermochemical models. A north/south asymmetry in the radial extent of CN is detected and we find that the CN emission is mostly optically thin and constrained vertically to a thin slab at $z/r \sim$0.5. A column density of 10$^{14}$\,cm$^{-2}$ is measured in the inner disk which for the north side decreases to 10$^{12}$\,cm$^{-2}$ and for the south side to 10$^{13}$\,cm$^{-2}$ in the outer regions. In Elias 2-27, CN traces a vertically elevated region above the midplane, very similar to that traced by $^{12}$CO. The inferred CN properties are consistent with thermo-chemical disk models, in which CN formation is initiated by the reaction of N with UV-pumped H$_2$. The observed north/south asymmetry may be caused by either ongoing infall or by a warped inner disk. This study highlights the importance of tracing the vertical location of various molecules to constrain the disk physical conditions. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A, 19 pages, 14 figures

Journal ref: A&A 666, A168 (2022)

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

Authors: Jaehan Bae, Richard Teague, Sean M. Andrews, Myriam Benisty, Stefano Facchini, Maria Galloway-Sprietsma, Ryan A. Loomis, Yuri Aikawa, Felipe Alarcon, Edwin Bergin, Jennifer B. Bergner, Alice S. Booth, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzman, Jane Huang, John D. Ilee, Nicolas T. Kurtovic, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Francois Menard, Karin I. Oberg , et al. (7 additional authors not shown)

Abstract: We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is… ▽ More We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted for publication in the ApJ Letters (July 7, 2022), 19 pages, 13 figures, interactive figures (Figure 7, 8, 9) are available at http://jaehanbae.com/as209/

Distributions of gas and small and large grains in the LkH$α\,330$ disk trace a young planetary system

Authors: P. Pinilla, M. Benisty, N. T. Kurtovic, J. Bae, R. Dong, Z. Zhu, S. Andrews, J. Carpenter, C. Ginski, J. Huang, A. Isella, L. Pérez, L. Ricci, G. Rosotti, M. Villenave, D. Wilner

Abstract: [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''… ▽ More [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$α\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction. △ Less

Submitted 25 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A, minor changes after language edition

Journal ref: A&A 665, A128 (2022)

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

Abstract: We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of… ▽ More We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

The morphology of CSCha circumbinary disk suggesting the existence of a Saturn-mass planet

Authors: N. T. Kurtovic, P. Pinilla, Anna B. T. Penzlin, M. Benisty, L. Pérez, C. Ginski, A. Isella, W. Kley, F. Menard, S. Pérez, A. Bayo

Abstract: Planets have been detected in circumbinary orbits in several different systems, despite the additional challenges faced during their formation in such an environment. We investigate the possibility of planetary formation in the spectroscopic binary CS Cha by analyzing its circumbinary disk. The system was studied with high angular resolution ALMA observations at 0.87mm. Visibilities modeling and K… ▽ More Planets have been detected in circumbinary orbits in several different systems, despite the additional challenges faced during their formation in such an environment. We investigate the possibility of planetary formation in the spectroscopic binary CS Cha by analyzing its circumbinary disk. The system was studied with high angular resolution ALMA observations at 0.87mm. Visibilities modeling and Keplerian fitting are used to constrain the physical properties of CS Cha, and the observations were compared to hydrodynamic simulations. Our observations are able to resolve the disk cavity in the dust continuum emission and the 12CO J:3-2 transition. We find the dust continuum disk to be azimuthally axisymmetric (less than 9% of intensity variation along the ring) and of low eccentricity (of 0.039 at the peak brightness of the ring). Under certain conditions, low eccentricities can be achieved in simulated disks without the need of a planet, however, the combination of low eccentricity and axisymmetry is consistent with the presence of a Saturn-like planet orbiting near the edge of the cavity. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 18 pages, 13 figures, accepted in A&A. Watch a 3min summary by the first author in https://youtu.be/6j5ROSj2N-U

Journal ref: A&A 664, A151 (2022)

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

Authors: Charles J. Law, Sage Crystian, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Kevin Flaherty, Viviana V. Guzmán, Jane Huang, John D. Ilee, Joel H. Kastner, Ryan A. Loomis, Feng Long, Laura M. Pérez, Sebastián Pérez, Chunhua Qi, Giovanni P. Rosotti, Dary Ruíz-Rodríguez, Takashi Tsukagoshi, David J. Wilner

Abstract: High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure o… ▽ More High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

Constraining cosmology with machine learning and galaxy clustering: the CAMELS-SAM suite

Authors: Lucia A. Perez, Shy Genel, Francisco Villaescusa-Navarro, Rachel S. Somerville, Austen Gabrielpillai, Daniel Anglés-Alcázar, Benjamin D. Wandelt, L. Y. Aaron Yung

Abstract: As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophy… ▽ More As the next generation of large galaxy surveys come online, it is becoming increasingly important to develop and understand the machine learning tools that analyze big astronomical data. Neural networks are powerful and capable of probing deep patterns in data, but must be trained carefully on large and representative data sets. We developed and generated a new `hump' of the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project: CAMELS-SAM, encompassing one thousand dark-matter only simulations of (100 $h^{-1}$ cMpc)$^3$ with different cosmological parameters ($Ω_m$ and $σ_8$) and run through the Santa Cruz semi-analytic model for galaxy formation over a broad range of astrophysical parameters. As a proof-of-concept for the power of this vast suite of simulated galaxies in a large volume and broad parameter space, we probe the power of simple clustering summary statistics to marginalize over astrophysics and constrain cosmology using neural networks. We use the two-point correlation function, count-in-cells, and the Void Probability Function, and probe non-linear and linear scales across $0.68<$ R $<27\ h^{-1}$ cMpc. Our cosmological constraints cluster around 3-8$\%$ error on $Ω_{\text{M}}$ and $σ_8$, and we explore the effect of various galaxy selections, galaxy sampling, and choice of clustering statistics on these constraints. We additionally explore how these clustering statistics constrain and inform key stellar and galactic feedback parameters in the Santa Cruz SAM. CAMELS-SAM has been publicly released alongside the rest of CAMELS, and offers great potential to many applications of machine learning in astrophysics: https://camels-sam.readthedocs.io. △ Less

Submitted 22 May, 2023; v1 submitted 5 April, 2022; originally announced April 2022.

Comments: 28 pages main text, additional 12 in appendices, 23 figures, 9 tables

Detection of Near-Infrared Water Ice at the Surface of the (pre)Transitional Disk of AB Aur: Informing Icy Grain Abundance, Composition, and Size

Authors: S. K. Betti, K. Follette, S. Jorquera, G. Duchêne, J. Mazoyer, M. Bonnefoy, G. Chauvin, L. M. Pérez, A. Boccaletti, C. Pinte, A. J. Weinberger, C. Grady, L. M. Close, D. Defrère, E. C. Downey, P. M. Hinz, F. Ménard, G. Schneider, A. J. Skemer, A. Vaz

Abstract: We present near-infrared Large Binocular Telescope Interferometer LMIRCam imagery of the disk around the Herbig Ae/Be star AB Aurigae. A comparison of surface brightness at Ks (2.16 $μ$m), H2O narrowband (3.08 $μ$m), and L' (3.7 $μ$m) allows us to probe the presence of icy grains in this (pre)transitional disk environment. By applying Reference Differential Imaging PSF subtraction, we detect the d… ▽ More We present near-infrared Large Binocular Telescope Interferometer LMIRCam imagery of the disk around the Herbig Ae/Be star AB Aurigae. A comparison of surface brightness at Ks (2.16 $μ$m), H2O narrowband (3.08 $μ$m), and L' (3.7 $μ$m) allows us to probe the presence of icy grains in this (pre)transitional disk environment. By applying Reference Differential Imaging PSF subtraction, we detect the disk at high signal to noise in all three bands. We find strong morphological differences between bands, including asymmetries consistent with observed spiral arms within 100 AU in L'. An apparent deficit of scattered light at 3.08 $μ$m relative to bracketing wavelengths (Ks and L') is evocative of ice absorption at the disk surface layer. However, the $Δ$(Ks-H2O) color is consistent with grains with little to no ice (0-5% by mass). The $Δ$(H2O-L') color, conversely, suggests grains with a much higher ice mass fraction (~0.68), and the two colors cannot be reconciled under a single grain population model. Additionally, we find the extremely red $Δ$(Ks-L') disk color cannot be reproduced under conventional scattered light modeling with any combination of grain parameters or reasonable local extinction values. We hypothesize that the scattering surfaces at the three wavelengths are not co-located, and optical depth effects result in each wavelength probing the grain population at different disk surface depths. The morphological similarity between Ks and H2O suggests their scattering surfaces are near one another, lending credence to the $Δ$(Ks-H2O) disk color constraint of < 5% ice mass fraction for the outermost scattering disk layer. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: Accepted for publication in AJ

LBT search for companions and sub-structures in the (pre)transitional disk of AB Aurigae

Authors: Sebastián Jorquera, Mickaël Bonnefoy, Sarah Betti, Gaël Chauvin, Esther Buenzli, Laura M. Pérez, Katherine B. Follette, Philip M. Hinz, Anthony Boccaletti, Vanessa Bailey, Beth Biller, Denis Defrère, Josh Eisner, Thomas Henning, Hubert Klahr, Jarron Leisenring, Johan Olofsson, Joshua E. Schlieder, Andrew J. Skemer, Michael F. Skrutskie, Roy Van Boekel

Abstract: Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and ne… ▽ More Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and near-infrared (SPHERE; 1.5-2.5$μ$m) wavelengths which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep $K_s$ (2.16$μ$m) and L' (3.7$μ$m) band images of AB Aurigae obtained with LMIRCam on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate ($ρ= 0.681", PA = 7.6^{\circ}$) had been previously claimed. The nature of a second innermost planet candidate ($ρ= 0.16'', PA = 203.9^{\circ}$) can not be investigated by the new data. We are able to derive 5$σ$ detection limits in both magnitude and mass for the system, going from 14 \Mjup at 0.3'' (49 au) down to 3-4 \Mjup at 0.6'' (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (< 0.5'') resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.5-0.7") showing a clear south-east/north-west asymmetry. This structure, observed for the first time at L'-band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk. △ Less

Submitted 10 February, 2022; v1 submitted 21 January, 2022; originally announced January 2022.

Comments: 13 pages, 5 figures, accepted on ApJ

An XMM-Newton catalogue of BL Lacs

Authors: I. de la Calle Pérez, N. Álvarez Crespo, E. Racero, A. Rouco

Abstract: We present an XMM-Newton catalogue of BL Lac X-ray, optical, and UV properties based on cross-correlation with the 1151 BL Lacs listed in the fifth edition of the Roma-BZCAT. We searched for the X-ray counterparts to these objects in the field of view of all pointed observations in the XMM-Newton archive over nearly 20 years of mission. The cross-correlation yields a total of 310 XMM-Newton fields… ▽ More We present an XMM-Newton catalogue of BL Lac X-ray, optical, and UV properties based on cross-correlation with the 1151 BL Lacs listed in the fifth edition of the Roma-BZCAT. We searched for the X-ray counterparts to these objects in the field of view of all pointed observations in the XMM-Newton archive over nearly 20 years of mission. The cross-correlation yields a total of 310 XMM-Newton fields which correspond to 103 different BL Lacs. We homogeneously analysed data from the three EPIC cameras (X-ray) and OM(optical/UV) using the XMM-Newton SAS software, and produced images, light curves, and spectral products for BL Lacs detected in any of the three EPIC cameras. We tested two different phenomenological models, log parabola and power law, with different variations of the absorbing column density and extracted their parameters. We derived time-variability information from the light curves following well-established statistical methods and quantified variability through statistical indicators. OM magnitudes and fluxes were computed wherever possible. We see that the log parabola model is preferred over the power law model for sources showing higher fluxes, which might indicate that curvature is intrinsic to BL Lacs and is only seen when the flux is high. We present the results of our analysis as a catalogue of X-ray spectral properties of the sample in the 0.2-10 keV energy band as well as in the optical/UV band. We complete the catalogue with multi-wavelength information at radio and gamma-ray energies. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 17 pages, 10 figures, 3 Tables. Tables B.1-B.4 and C.1 are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/655/A48

Journal ref: A&A 655, A48 (2021)

The CAMELS project: public data release

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch , et al. (22 additional authors not shown)

Abstract: The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present… ▽ More The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$α$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at https://www.camel-simulations.org

Probing inner and outer disk misalignments in transition disks

Authors: A. J. Bohn, M. Benisty, K. Perraut, N. van der Marel, L. Wölfer, E. F. van Dishoeck, S. Facchini, C. F. Manara, R. Teague, L. Francis, J-P. Berger, R. Garcia-Lopez, C. Ginski, T. Henning, M. Kenworthy, S. Kraus, F. Ménard, A. Mérand, L. M. Pérez

Abstract: For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and… ▽ More For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and $^{13}$CO molecular line data from ALMA to derive the orientation of the outer disk regions. We fit simple models to the GRAVITY data to derive the inner disks inclination and position angles. The outer disk geometries were derived from Keplerian fits to the ALMA velocity maps and compared to the inner disk constraints. We also predicted the locations of shadows for significantly misaligned systems. Our analysis reveals six disks to exhibit significant misalignments between their inner and outer disks. The predicted shadow positions agree well with the scattered light images of HD100453 and HD142527, and we find supporting evidence for a shadow in the disk around CQ Tau. In the other three targets for which we infer significantly misaligned disks, V1247 Ori, V1366 Ori, and RY Lup, we do not see any evident sign of shadows in the scattered light images. The scattered light shadows observed in DoAr44, HD135344B, and HD139614 are consistent with our observations, yet the underlying morphology is likely too complex to be described by our models and the accuracy achieved by our observations. Whereas we can derive precise constraints on the potential shadow positions for well-resolved inner disks around HAeBe stars, the statistical uncertainties for the marginally resolved inner disks around the TTS of our sample make it difficult to extract conclusive constraints for the presence of shadows in these systems. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: Accepted for publication in Astronomy \& Astrophysics

Journal ref: A&A 658, A183 (2022)

New spectroscopic confirmations of Lyman-$α$ emitters at z $\sim$ 7 from the LAGER survey

Authors: Santosh Harish, Isak G. B. Wold, Sangeeta Malhotra, James Rhoads, Weida Hu, Junxian Wang, Zhen-ya Zheng, L. Felipe Barrientos, Jorge González-López, Lucia A. Perez, Ali Ahmad Khostovan, Leopoldo Infante, Chunyan Jiang, Cristóbal Moya-Sierralta, John Pharo, Francisco Valdes, Huan Yang

Abstract: We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sen… ▽ More We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sensitive probes of cosmic reionization and narrow-band imaging is a robust and effective method for selecting a large number of LAEs. In this work, we present results from the spectroscopic follow-up of LAE candidates in two LAGER fields, COSMOS and WIDE-12, using observations from Keck/LRIS. We report the successful detection of Ly$α$ emission in 15 candidates (11 in COSMOS and 4 in WIDE-12 fields). Three of these in COSMOS have matching confirmations from a previous LAGER spectroscopic follow-up and are part of the overdense region, LAGER-$z7$OD1. Additionally, two candidates that were not detected in the LRIS observations have prior spectroscopic confirmations from Magellan. Including these, we obtain a spectroscopic confirmation success rate of $\sim$$80$% for LAGER LAE candidates. Apart from Ly$α$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$σ$ upper limit for the ratio of NV/Ly$α$, $f_{NV}/f_{Lyα} \lesssim 0.27$, which implies that ionizing emission from these sources is mostly dominated by star formation. Including confirmations from this work, a total of 33 LAE sources from LAGER are now spectroscopically confirmed. LAGER has more than doubled the sample of spectroscopically confirmed LAE sources at $z\sim7$. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 8 pages, 4 figures; Submitted to ApJ

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

Authors: Francisco Villaescusa-Navarro, Shy Genel, Daniel Angles-Alcazar, Leander Thiele, Romeel Dave, Desika Narayanan, Andrina Nicola, Yin Li, Pablo Villanueva-Domingo, Benjamin Wandelt, David N. Spergel, Rachel S. Somerville, Jose Manuel Zorrilla Matilla, Faizan G. Mohammad, Sultan Hassan, Helen Shao, Digvijay Wadekar, Michael Eickenberg, Kaze W. K. Wong, Gabriella Contardo, Yongseok Jo, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Lucia A. Perez , et al. (3 additional authors not shown)

Abstract: We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light year… ▽ More We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: 17 pages, 1 figure. Third paper of a series of four. Hundreds of thousands of labeled 2D maps and 3D grids from thousands of simulated universes publicly available at https://camels-multifield-dataset.readthedocs.io

Molecules with ALMA at Planet-forming Scales (MAPS) XIV: Revealing disk substructures in multi-wavelength continuum emission

Authors: Anibal Sierra, Laura M. Pérez, Ke Zhang, Charles J. Law, Viviana V. Guzmán, Chunhua Qi, Arthur D. Bosman, Karin I. Öberg, Sean M. Andrews, Feng Long, Richard Teague, Alice S. Booth, Catherine Walsh, David J. Wilner, François Ménard, Gianni Cataldi, Ian Czekala, Jaehan Bae, Jane Huang, Jennifer B. Bergner, John D. Ilee, Myriam Benisty, Romane Le Gal, Ryan A. Loomis, Takashi Tsukagoshi , et al. (3 additional authors not shown)

Abstract: Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distributi… ▽ More Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distribution of solid mass and maximum grain size in each disk, including or excluding dust scattering. In the nonscattering model, the dust surface density and maximum grain size profiles decrease from the inner disks to the outer disks, with local maxima at the bright ring locations, as expected from dust trapping models. The inferred maximum grain sizes from the inner to outer disks decrease from ~1 cm to 1 mm. For IM Lup, HD 163296, and MWC 480 in the scattering model, two solutions are compatible with their observed inner disk emission: one solution corresponding to a maximum grain size of a few millimeters (similar to the nonscattering model), and the other corresponding to a few hundred micrometer sizes. Based on the estimated Toomre parameter, only IM Lup -- which shows a prominent spiral morphology in millimeter dust -- is found to be gravitationally unstable. The estimated maximum Stokes number in all the disks lies between 0.01 and 0.3, and the estimated turbulence parameters in the rings of AS 209 and HD 163296 are close to the threshold where dust growth is limited by turbulent fragmentation. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 20 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Molecules with ALMA at Planet-forming Scales (MAPS) VI: Distribution of the small organics HCN, C2H, and H2CO

Authors: Viviana V. Guzmán, Jennifer B. Bergner, Charles J. Law, Karin I. Oberg, Catherine Walsh, Gianni Cataldi, Yuri Aikawa, Edwin A. Bergin, Ian Czekala, Jane Huang, Sean M. Andrews, Ryan A. Loomis, Ke Zhang, Romane Le Gal, Felipe Alarcón, John D. Ilee, Richard Teague, L. Ilsedore Cleeves, David J. Wilner, Feng Long, Kamber R. Schwarz, Arthur D. Bosman, Laura M. Pérez, François Ménard, Yao Liu

Abstract: Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find… ▽ More Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find that the HCN emission arises in a temperate (20-30 K) layer in the disk, while C2H is present in relatively warmer (20-60 K) layers. In the case of HD 163296, we find a decrease in column density for HCN and C2H inside one of the dust gaps near 83 au, where a planet has been proposed to be located. We derive H2CO column density profiles assuming temperatures between 20 and 50 K, and find slightly higher column densities in the colder disks around T Tauri stars than around Herbig Ae stars. The H2CO column densities rise near the location of the CO snowline and/or millimeter dust edge, suggesting an efficient release of H2CO ices in the outer disk. Finally, we find that the inner 50 au of these disks are rich in organic species, with abundances relative to water that are similar to cometary values. Comets could therefore deliver water and key organics to future planets in these disks, similar to what might have happened here on Earth. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 23 pages, 10 figures, 4 tables. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS) I: Program Overview and Highlights

Authors: Karin I. Oberg, Viviana V. Guzman, Catherine Walsh, Yuri Aikawa, Edwin A. Bergin, Charles J. Law, Ryan A. Loomis, Felipe Alarcon, Sean M. Andrews, Jaehan Bae, Jennifer B. Bergner, Yann Boehler, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Kenji Furuya, Jane Huang, John D. Ilee, Nicolas T. Kurtovic, Romane Le Gal, Yao Liu, Feng Long , et al. (13 additional authors not shown)

Abstract: Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission fro… ▽ More Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks - around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 - in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical sub-structures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Accepted for publication in the ApJS MAPS Special Issue. v2 has updated MAPS references and a correction to Fig. 3

Molecules with ALMA at Planet-forming Scales (MAPS) V: CO gas distributions

Authors: Ke Zhang, Alice S. Booth, Charles J. Law, Arthur D. Bosman, Kamber R. Schwarz, Edwin A. Bergin, Karin I. Öberg, Sean M. Andrews, Viviana V. Guzmán, Catherine Walsh, Chunhua Qi, Merel L. R. van 't Hoff, Feng Long, David J. Wilner, Jane Huang, Ian Czekala, John D. Ilee, Gianni Cataldi, Jennifer B. Bergner, Yuri Aikawa, Richard Teague, Jaehan Bae, Ryan A. Loomis, Jenny K. Calahan, Felipe Alarcón , et al. (10 additional authors not shown)

Abstract: Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling f… ▽ More Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 23 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. 36 pages, 21 figures, accepted for publication in ApJS

Molecules with ALMA at Planet-forming Scales (MAPS). XV. Tracing protoplanetary disk structure within 20 au

Authors: Arthur D. Bosman, Edwin A. Bergin, Ryan A. Loomis, Sean M. Andrews, Merel L. R. van 't Hoff, Richard Teague, Karin I. Öberg, Viviana V. Guzmán, Catherine Walsh, Yuri Aikawa, Felipe Alarcón, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Hideko Nomura , et al. (3 additional authors not shown)

Abstract: Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large p… ▽ More Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large program to construct radial surface brightness profiles with a ~3 au effective resolution for the CO isotopologue J=2-1 lines using the line velocity profile. IM Lup reveals a central depression in 13CO and C18O that is ascribed to a pileup of ~500 $M_\oplus$ of dust in the inner 20 au, leading to a gas-to-dust ratio of around <10. This pileup is consistent with efficient drift of grains ($\gtrsim$ 100 $M_\oplus$ Myr$^{-1}$) and a local gas-to-dust ratio that suggests that the streaming instability could be active. The CO isotopologue emission in the GM Aur disk is consistent with a small (~15 au), strongly depleted gas cavity within the ~40 au dust cavity. The radial surface brightness profiles for both the AS 209 and HD 163296 disks show a local minimum and maximum in the C18O emission at the location of a known dust ring (~14 au) and gap (~10 au), respectively. This indicates that the dust ring has a low gas-to-dust ratio ($>$ 10) and that the dust gap is gas-rich enough to have optically thick C18O. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 22 pages, 14 figures, accepted by ApJS. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Update with correct references to other MAPS papers

Molecules with ALMA at Planet-forming Scales (MAPS XVIII): Kinematic Substructures in the Disks of HD 163296 and MWC 480

Authors: Richard Teague, Jaehan Bae, Yuri Aikawa, Sean M. Andrews, Edwin A. Bergin, Jennifer B. Bergner, Yann Boehler, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, Ian Czekala, Viviana V. Guzmán, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Feng Long, Ryan A. Loomis, François Ménard, Karin I. Öberg, Laura M. Pérez, Kamber R. Schwarz, Anibal Sierra, Catherine Walsh, David J. Wilner , et al. (2 additional authors not shown)

Abstract: We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Kepleria… ▽ More We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, $v_φ$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 5 \%$ relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 20 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 32 pages, 21 figures, accepted for publication in ApJS, MAPS cross-references updated

Molecules with ALMA at Planet-forming Scales (MAPS) IV: Emission Surfaces and Vertical Distribution of Molecules

Authors: Charles J. Law, Richard Teague, Ryan A. Loomis, Jaehan Bae, Karin I. Öberg, Ian Czekala, Sean M. Andrews, Yuri Aikawa, Felipe Alarcón, Edwin A. Bergin, Jennifer B. Bergner, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Kenji Furuya, Viviana V. Guzmán, Jane Huang, John D. Ilee, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Hideko Nomura , et al. (10 additional authors not shown)

Abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a funct… ▽ More The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a function of radius ($r$). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C$_2$H. We find that $^{12}$CO emission traces the most elevated regions with $z/r > 0.3$, while emission from the less abundant $^{13}$CO and C$^{18}$O probes deeper into the disk at altitudes of $z/r \lesssim 0.2$. C$_2$H and HCN have lower opacities and SNRs, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with $z/r \lesssim 0.1$, i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermo-chemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 20 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 31 pages, 20 figures, accepted for publication in ApJS, MAPS cross-references updated

Molecules with ALMA at Planet-forming Scales (MAPS) II: CLEAN Strategies for Synthesizing Images of Molecular Line Emission in Protoplanetary Disks

Authors: Ian Czekala, Ryan A. Loomis, Richard Teague, Alice S. Booth, Jane Huang, Gianni Cataldi, John D. Ilee, Charles J. Law, Catherine Walsh, Arthur D. Bosman, Viviana V. Guzmán, Romane Le Gal, Karin I. Öberg, Yoshihide Yamato, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, L. Ilsedore Cleeves, Nicolas T. Kurtovic, François Ménard, Hideko Nomura, Laura M. Pérez, Chunhua Qi , et al. (5 additional authors not shown)

Abstract: The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our mu… ▽ More The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our multi-stage workflow -- built around the CASA tclean image deconvolution procedure -- that we used to generate the core data product of the MAPS LP: the position-position-velocity image cubes for each spectral line. Owing to the expansive nature of the survey, we encountered a range of imaging challenges; some are familiar to the sub-mm protoplanetary disk community, like the benefits of using an accurate CLEAN mask, and others less well-known, like the incorrect default flux scaling of the CLEAN residual map first described in Jorsater & van Moorsel 1995 (the "JvM effect"). We distill lessons learned into recommended workflows for synthesizing image cubes of molecular emission. In particular, we describe how to produce image cubes with accurate fluxes via the "JvM correction," a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low S/N emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of signal to noise ratio when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 24 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: ApJS accepted. Part of the MAPS ALMA large program series: http://www.alma-maps.info/. Updated bibliography with MAPS LP arXiv references

A Circumplanetary Disk Around PDS70c

Authors: Myriam Benisty, Jaehan Bae, Stefano Facchini, Miriam Keppler, Richard Teague, Andrea Isella, Nicolas T. Kurtovic, Laura M. Perez, Anibal Sierra, Sean M. Andrews, John Carpenter, Ian Czekala, Carsten Dominik, Thomas Henning, Francois Menard, Paola Pinilla, Alice Zurlo

Abstract: PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$α$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter… ▽ More PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$α$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 855$μ$m at high angular resolution ($\sim$20mas, 2.3au) that aim to resolve the circumplanetary disks and constrain their dust masses. Our observations confirm the presence of a compact source of emission co-located with PDS70c, spatially separated from the circumstellar disk and less extended than $\sim$1.2au in radius, a value close to the expected truncation radius of the cicumplanetary disk at a third of the Hill radius. The emission around PDS70c has a peak intensity of $\sim$86$\pm$16 $μ\mathrm{Jy}~\mathrm{beam}^{-1}$ which corresponds to a dust mass of $\sim$0.031M$_{\oplus}$ or $\sim$0.007M$_{\oplus}$, assuming that it is only constituted of 1 $μ$m or 1 mm sized grains, respectively. We also detect extended, low surface brightness continuum emission within the cavity near PDS70b. We observe an optically thin inner disk within 18au of the star with an emission that could result from small micron-sized grains transported from the outer disk through the orbits of b and c. In addition, we find that the outer disk resolves into a narrow and bright ring with a faint inner shoulder. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: ApJ Letters, in press; 19 pages, 9 figures

Journal ref: 2021, ApJ, 916L, 2B

LAGER Ly$α$ Luminosity Function at $z\sim7$, Implications for Reionization

Authors: Isak G. B. Wold, Sangeeta Malhotra, James Rhoads, Junxian Wang, Weida Hu, Lucia A. Perez, Zhen-Ya Zheng, Ali Ahmad Khostovan, Alistair R. Walker, L. Felipe Barrientos, Jorge González-López, Santosh Harish, Leopoldo Infante, Chunyan Jiang, John Pharo, Cristóbal Moya-Sierralta, Franz E. Bauer, Gaspar Galaz, Francisco Valdes, Huan Yang

Abstract: We present a new measurement of the Ly$α$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest… ▽ More We present a new measurement of the Ly$α$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest $z\sim 7$ survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify $N=95$ new $z=6.9$ Ly$α$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$α$ luminosity functions. The best-fit Schechter luminosity function parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Ly$α$ luminosity density declines at the same rate as the UV continuum LF from $z=5.7$ to $z=6.9$. This is consistent with an intergalactic medium that was fully ionized as early as redshift $z\sim 7$, or with a volume-averaged neutral hydrogen fraction of $x_{HI} < 0.33$ at $1σ$. △ Less

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

Comments: 19 pages, 11 figures; accepted by ApJ

Limits on Millimeter Continuum Emission from Circumplanetary Material in the DSHARP Disks

Authors: Sean M. Andrews, William Elder, Shangjia Zhang, Jane Huang, Myriam Benisty, Nicolás T. Kurtovic, David J. Wilner, Zhaohuan Zhu, John M. Carpenter, Laura M. Pérez, Richard Teague, Andrea Isella, Luca Ricci

Abstract: We present a detailed analysis for a subset of the high resolution (~35 mas, or 5 au) ALMA observations from the Disk Substructures at High Angular Resolution Project (DSHARP) to search for faint 1.3 mm continuum emission associated with dusty circumplanetary material located within the narrow annuli of depleted emission (gaps) in circumstellar disks. This search used the Jennings et al. (2020)… ▽ More We present a detailed analysis for a subset of the high resolution (~35 mas, or 5 au) ALMA observations from the Disk Substructures at High Angular Resolution Project (DSHARP) to search for faint 1.3 mm continuum emission associated with dusty circumplanetary material located within the narrow annuli of depleted emission (gaps) in circumstellar disks. This search used the Jennings et al. (2020) $\tt{frank}$ modeling methodology to mitigate contamination from the local disk emission, and then deployed a suite of injection-recovery experiments to statistically characterize point-like circumplanetary disks in residual images. While there are a few putative candidates in this sample, they have only marginal local signal-to-noise ratios and would require deeper measurements to confirm. Associating a 50% recovery fraction with an upper limit, we find these data are sensitive to circumplanetary disks with flux densities $\gtrsim 50-70$ $μ$Jy in most cases. There are a few examples where those limits are inflated ($\gtrsim 110$ $μ$Jy) due to lingering non-axisymmetric structures in their host circumstellar disks, most notably for a newly identified faint spiral in the HD 143006 disk. For standard assumptions, this analysis suggests that these data should be sensitive to circumplanetary disks with dust masses $\gtrsim 0.001-0.2$ M$_\oplus$. While those bounds are comparable to some theoretical expectations for young giant planets, we discuss how plausible system properties (e.g., relatively low host planet masses or the efficient radial drift of solids) could require much deeper observations to achieve robust detections. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: ApJ, in press; 30 pages, 15 figures (one is a 4-part figure set)

The Architecture of the V892 Tau System: the Binary and its Circumbinary Disk

Authors: Feng Long, Sean M. Andrews, Justin Vega, David J. Wilner, Claire J. Chandler, Enrico Ragusa, Richard Teague, Laura M. Pérez, Nuria Calvet, John M. Carpenter, Thomas Henning, Woojin Kwon, Hendrik Linz, Luca Ricci

Abstract: We present high resolution millimeter continuum and CO line observations for the circumbinary disk around V892 Tau to constrain the stellar and disk properties. The total mass of the two near-equal-mass A stars is estimated to be $6.0\pm0.2\,M_{\odot}$ based on our models of the Keplerian-dominated gas disk rotation. The detection of strong ionized gas emission associated with the two stars at 8 m… ▽ More We present high resolution millimeter continuum and CO line observations for the circumbinary disk around V892 Tau to constrain the stellar and disk properties. The total mass of the two near-equal-mass A stars is estimated to be $6.0\pm0.2\,M_{\odot}$ based on our models of the Keplerian-dominated gas disk rotation. The detection of strong ionized gas emission associated with the two stars at 8 mm, when combined with previous astrometric measurements in the near-infrared, provides an updated view of the binary orbit with $a=7.1\pm0.1$ au, $e=0.27\pm0.1$, and $P=7.7\pm0.2$ yr, which is about half of a previously reported orbital period. The binary orbital plane is proposed to be near coplanar to the circumbinary disk plane (with a mutual inclination of only $Δ=8\pm4.2$ deg; another solution with $Δ=113$ deg is less likely given the short re-alignment timescale). An asymmetric dust disk ring peaking at a radius of 0.''2 is detected at 1.3 mm and its fainter counterparts are also detected at the longer 8 and 9.8 mm. The CO gas disk, though dominated by Keplerian rotation, presents a mild inner and outer disk misalignment, such that the inner disk to the SW and outer disk to the NE appear brighter than their counterparts at the opposite disk sides. The radial extension of the disk, its asymmetric dust ring, and the presence of a disk warp could all be explained by the interaction between the eccentric binary and the circumbinary disk, which we assume were formed with non-zero mutual inclination. Some tentatively detected gas spirals in the outer disk are likely produced by interactions with the low mass tertiary component located 4'' to the northeast. Our analyses demonstrate the promising usage of V892 Tau as an excellent benchmark system to study the details of binary--disk interactions. △ Less

Submitted 6 May, 2021; originally announced May 2021.

Comments: Accepted for publication in ApJ

A dynamical measurement of the disk mass in Elias 2-27

Authors: Benedetta Veronesi, Teresa Paneque-Carreno, Giuseppe Lodato, Leonardo Testi, Laura M. Pérez, Giuseppe Bertin, Cassandra Hall

Abstract: Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is parti… ▽ More Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is particularly interesting if we consider that gravitational instabilities could be a fundamental first step for the formation of planetesimals and planets. In this Letter, we model the rotation curve obtained by CO data of Elias 2-27 with a theoretical rotation curve including both the disk self-gravity and the star contribution to the gravitational potential. We compare this model with a purely Keplerian one and with a simple power-law function. We find that (especially for the $^{13}$CO isotopologue) the rotation curve is better described by considering not only the star, but also the disk self-gravity. We are thus able to obtain for the first time a dynamical estimate of the disk mass of $0.08\pm0.04\,M_{\odot}$ and the star mass of $0.46\pm0.03\,M_{\odot}$ (in the more general case), the latter being comparable with previous estimates. From these values, we derive that the disk is 17$\%$ of the star mass, meaning that it could be prone to gravitational instabilities. This result would strongly support the hypothesis that the two spiral arms are generated by gravitational instabilities. △ Less

Submitted 7 May, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 14 pages, 9 figures, 2 tables, accepted for publication on ApJL

Spiral Arms and a Massive Dust Disk with non-Keplerian Kinematics: Possible Evidence for Gravitational Instability in the Disk of Elias 2-27

Authors: T. Paneque-Carreño, L. M. Perez, M. Benisty, C. Hall, B. Veronesi, G. Lodato, A. Sierra, J. M. Carpenter, S. M. Andrews, Jaehan Bae, Th. Henning, W. Kwon, H. Linz, L. Loinard, C. Pinte, L. Ricci, M. Tazzari, L. Testi, D. Wilner

Abstract: To determine the origin of the spiral structure observed in the dust continuum emission of Elias 2-27 we analyze multi-wavelength continuum ALMA data with a resolution of $\sim$0.2 arcsec ($\sim$23au) at 0.89, 1.3 and 3.3mm. We also study the kinematics of the disk with $^{13}$CO and C$^{18}$O ALMA observations in the $J=$3-2 transition. The spiral arm morphology is recovered at all wavelengths in… ▽ More To determine the origin of the spiral structure observed in the dust continuum emission of Elias 2-27 we analyze multi-wavelength continuum ALMA data with a resolution of $\sim$0.2 arcsec ($\sim$23au) at 0.89, 1.3 and 3.3mm. We also study the kinematics of the disk with $^{13}$CO and C$^{18}$O ALMA observations in the $J=$3-2 transition. The spiral arm morphology is recovered at all wavelengths in the dust continuum observations, where we measure contrast and spectral index variations along the spiral arms and detect subtle dust-trapping signatures. We determine that the emission from the midplane is cold and interpret the optical depth results as signatures of a higher disk mass than previous constraints. From the gas data, we search for deviations from Keplerian motion and trace the morphology of the emitting surfaces and the velocity profiles. We find an azimuthally varying emission layer height in the system, large-scale emission surrounding the disk, and strong perturbations in the channel maps, co-located with the spirals. Additionally, we develop multigrain dust and gas SPH simulations of a gravitationally unstable disk and compare them to the observations. Given the large scale emission and highly perturbed gas structure, together with the comparison of continuum observations to theoretical predictions, we propose infall-triggered gravitational instabilities as origin for the observed spiral structure. △ Less

Submitted 24 May, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: Accepted for publication in ApJ, 38 pages, 23 Figures

Correlations between H$α$ Equivalent Width and Galaxy Properties at $z = 0.47$: Physical or Selection-driven?

Authors: Ali Ahmad Khostovan, Sangeeta Malhotra, James E. Rhoads, Santosh Harish, Chunyan Jiang, Junxian Wang, Isak Wold, Zhen-Ya Zheng, L. Felipe Barrientos, Alicia Coughlin, Weida Hu, Leopoldo Infante, Lucia A. Perez, John Pharo, Francisco Valdes, Alistair R. Walker

Abstract: The H$α$ equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic star-formation activity. Previous assessments show that EW strongly anti-correlates with stellar mass as $M^{-0.25}$ similar to the sSFR -- stellar mass relation. However, such a correlation may be driven/formed by selection effects. In this study, we investigate how H$α$ EWs c… ▽ More The H$α$ equivalent width (EW) is an observational proxy for specific star formation rate (sSFR) and a tracer of episodic star-formation activity. Previous assessments show that EW strongly anti-correlates with stellar mass as $M^{-0.25}$ similar to the sSFR -- stellar mass relation. However, such a correlation may be driven/formed by selection effects. In this study, we investigate how H$α$ EWs correlate with galaxy properties and how selection biases could alter such correlations using a $z = 0.47$ narrowband-selected sample of 1572 H$α$ emitters from the Ly$α$ Galaxies in the Epoch of Reionization (LAGER) survey. The sample covers 3 deg$^2$ of COSMOS and $1.1\times10^5$ cMpc$^3$. We assume an intrinsic EW distribution to form mock samples of H$α$ emitters (HAEs) and propagate the selection criteria to match observations, giving us control on how selection biases can affect the underlying results. We find EW intrinsically correlates with stellar mass as $W_0 \propto M^{-0.16\pm0.03}$ and decreases by a factor of $\sim 3$ from $10^{7}$ to $10^{10}$ M$_\odot$. We find low-mass HAEs to be $\sim 320$ times more likely to have rest-frame EW$ > 200$Åcompared to high-mass HAEs. Combining the intrinsic EW -- stellar mass correlation with an observed SMF correctly reproduces the observed H$α$ LF, while not correcting for selection effects underestimates the number of bright HAEs. This suggests that the intrinsic EW -- stellar mass correlation is physically significant and reproduces three statistical distributions of galaxy populations (LF, SMF, EW distribution). At lower masses, we find there are more high-EW outliers compared to high masses, even after taking into account selection effects. Our results suggest that high sSFR outliers indicative of bursty SF activity are intrinsically more prevalent in low-mass HAEs and not a byproduct of selection effects. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 21 Pages, 13 Figures, 2 Tables, accepted to MNRAS

A Lyman-α protocluster at redshift 6.9

Authors: Weida Hu, Junxian Wang, Leopoldo Infante, James E. Rhoads, Zhen-Ya Zheng, Huan Yang, Sangeeta Malhotra, L. Felipe Barrientos, Chunyan Jiang, Jorge González-López, Gonzalo Prieto, Lucia A. Perez, Pascale Hibon, Gaspar Galaz, Alicia Coughlin, Santosh Harish, Xu Kong, Wenyong Kang, Ali Ahmad Khostovan, John Pharo, Francisco Valdes, Isak Wold, Alistair R. Walker, XianZhong Zheng

Abstract: Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million ye… ▽ More Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million years old and could be experiencing rapid evolution of the neutral hydrogen fraction in the intergalactic medium. The protocluster is identified by an overdensity of 6 times the average galaxy density, and with 21 narrowband selected Lyman-$α$ galaxies, among which 16 have been spectroscopically confirmed. At redshifts similar to or above this record, smaller protogroups with fewer members have been reported. LAGER-z7OD1 shows an elongated shape and consists of two subprotoclusters, which would have merged into one massive cluster with a present-day mass of $3.7 \times 10^{15}$ solar masses. The total volume of the ionized bubbles generated by its member galaxies is found to be comparable to the volume of the protocluster itself, indicating that we are witnessing the merging of the individual bubbles and that the intergalactic medium within the protocluster is almost fully ionized. LAGER-z7OD1 thus provides a unique natural laboratory to investigate the reionization process. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: 9 pages, 5 figures, and 1 table. Accepted by Nature Astronomy

A search for companions via direct imaging in the DSHARP planet-forming disks

Authors: Sebastián Jorquera, Laura M. Pérez, Gaël Chauvin, Myriam Benisty, Zhaohuan Zhu, Andrea Isella, Jane Huang, Luca Ricci, Sean M. Andrews, Shangjia Zhang, John Carpenter, Nicolás T. Kurtovic, Tilman Birnstiel

Abstract: The "Disk Substructures at High Angular Resolution Project" (DSHARP) has revealed an abundance and ubiquity of rings and gaps over a large sample of young planet-forming disks, which are hypothesised to be induced by the presence of forming planets. In this context, we present the first attempt to directly image these young companions for 10 of the DSHARP disks, by using NaCo/VLT high contrast obs… ▽ More The "Disk Substructures at High Angular Resolution Project" (DSHARP) has revealed an abundance and ubiquity of rings and gaps over a large sample of young planet-forming disks, which are hypothesised to be induced by the presence of forming planets. In this context, we present the first attempt to directly image these young companions for 10 of the DSHARP disks, by using NaCo/VLT high contrast observations in L'-band instrument and angular differential imaging techniques. We report the detection of a point-like source candidate at 1.1" (174.9 au) for RU Lup, and at 0.42" (55 AU) for Elias 24. In the case of RU Lup, the proper motion of the candidate is consistent with a stationary background contaminant, based on the astrometry derived from our observations and available archival data. For Elias 24 the point-like source candidate is located in one of the disk gaps at 55 AU. Assuming it is a planetary companion, our analysis suggest a mass ranging from $0.5 M_J$ up to $5 M_J$, depending on the presence of a circumplanetary disk and its contribution to the luminosity of the system. However, no clear confirmation is obtained at this stage, and follow-up observations are mandatory to verify if the proposed source is physical, comoving with the stellar host, and associated with a young massive planet sculpting the gap observed at 55\,AU. For all the remaining systems, the lack of detections suggests the presence of planetary companions with masses lower than $5M_J$, based on our derived mass detection limits. This is consistent with predictions of both hydrodynamical simulations and kinematical signatures on the disk, and allows us to set upper limits on the presence of massive planets in these young disks. △ Less

Submitted 31 March, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 19 pages, 10 figures; Published by AJ

Journal ref: The Astronomical Journal (2021), Volume 161, Number 3

The 2017 May 20$^{\rm th}$ stellar occultation by the elongated centaur (95626) 2002 GZ$_{32}$

Authors: P. Santos-Sanz, J. L. Ortiz, B. Sicardy, G. Benedetti-Rossi, N. Morales, E. Fernández-Valenzuela, R. Duffard, R. Iglesias-Marzoa, J. L. Lamadrid, N. Maícas, L. Pérez, K. Gazeas, J. C. Guirado, V. Peris, F. J. Ballesteros, F. Organero, L. Ana-Hernández, F. Fonseca, A. Alvarez-Candal, Y. Jiménez-Teja, M. Vara-Lubiano, F. Braga-Ribas, J. I. B. Camargo, J. Desmars, M. Assafin , et al. (34 additional authors not shown)

Abstract: We predicted a stellar occultation of the bright star Gaia DR1 4332852996360346368 (UCAC4 385-75921) (m$_{\rm V}$= 14.0 mag) by the centaur 2002 GZ$_{32}$ for 2017 May 20$^{\rm th}$. Our latest shadow path prediction was favourable to a large region in Europe. Observations were arranged in a broad region inside the nominal shadow path. Series of images were obtained with 29 telescopes throughout E… ▽ More We predicted a stellar occultation of the bright star Gaia DR1 4332852996360346368 (UCAC4 385-75921) (m$_{\rm V}$= 14.0 mag) by the centaur 2002 GZ$_{32}$ for 2017 May 20$^{\rm th}$. Our latest shadow path prediction was favourable to a large region in Europe. Observations were arranged in a broad region inside the nominal shadow path. Series of images were obtained with 29 telescopes throughout Europe and from six of them (five in Spain and one in Greece) we detected the occultation. This is the fourth centaur, besides Chariklo, Chiron and Bienor, for which a multi-chord stellar occultation is reported. By means of an elliptical fit to the occultation chords we obtained the limb of 2002 GZ$_{32}$ during the occultation, resulting in an ellipse with axes of 305 $\pm$ 17 km $\times$ 146 $\pm$ 8 km. From this limb, thanks to a rotational light curve obtained shortly after the occultation, we derived the geometric albedo of 2002 GZ$_{32}$ ($p_{\rm V}$ = 0.043 $\pm$ 0.007) and a 3-D ellipsoidal shape with axes 366 km $\times$ 306 km $\times$ 120 km. This shape is not fully consistent with a homogeneous body in hydrostatic equilibrium for the known rotation period of 2002 GZ$_{32}$. The size (albedo) obtained from the occultation is respectively smaller (greater) than that derived from the radiometric technique but compatible within error bars. No rings or debris around 2002 GZ$_{32}$ were detected from the occultation, but narrow and thin rings cannot be discarded. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: Accepted for publication in MNRAS (8-Dec.-2020), 15 pages, 9 figures

The Ophiuchus DIsc Survey Employing ALMA (ODISEA)-III: the evolution of substructures in massive discs at 3-5 au resolution

Authors: Lucas A. Cieza, Camilo González-Ruilova, Antonio S. Hales, Paola Pinilla, Dary Ruíz-Rodríguez, Alice Zurlo, Simón Casassus, Sebastián Pérez, Hector Cánovas, Carla Arce-Tord, Mario Flock, Nicolas Kurtovic, Sebastian Marino, Pedro H. Nogueira, Laura Perez, Daniel J. Price, David A. Principe, Jonathan P. Williams

Abstract: We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in… ▽ More We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in more evolved sources with Class II SEDs. Our sample includes 5 objects with previously known large dust cavities (r $>$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques. △ Less

Submitted 5 December, 2020; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

Void Probability Function of Simulated Surveys of high-redshift Lyman-Alpha Emitters

Authors: Lucia A. Perez, Sangeeta Malhotra, James E. Rhoads, Vithal Tilvi

Abstract: We calculate the void probability function (VPF) in simulations of Lyman-$α$ emitters (LAEs) across a wide redshift range ($z=3.1,\ 4.5,\ 5.7,\ 6.6$). The VPF measures the zero-point correlation function (i.e. places devoid of galaxies) and naturally connects to higher order correlation functions while being computationally simple to calculate. We explore the Poissonian and systematic errors on th… ▽ More We calculate the void probability function (VPF) in simulations of Lyman-$α$ emitters (LAEs) across a wide redshift range ($z=3.1,\ 4.5,\ 5.7,\ 6.6$). The VPF measures the zero-point correlation function (i.e. places devoid of galaxies) and naturally connects to higher order correlation functions while being computationally simple to calculate. We explore the Poissonian and systematic errors on the VPF, specify its accuracy as a function of average source density and the volume probed, and provide the appropriate size scales to measure the VPF. At small radii the accuracy of the VPF is limited by galaxy density, while at large radii the VPF is limited by the number of independent volumes probed. We also offer guidelines for understanding and quantifying the error in the VPF. We approximate the error in the VPF by using independent sub-volumes of the catalogs, after finding that jackknife statistics underestimate the uncertainty. We use the VPF to probe the strength of higher order correlation functions by measuring and examining the hierarchical scaling between the correlation functions using count-in-cells. The negative binomial model (NBM) has been shown to best describe the scaling between the two point correlation function and VPF for low-redshift galaxy observations. We further test the fit of the NBM by directly deriving the volume averaged two-point correlation function from the VPF and vice versa. We find the NBM best describes the $z=3.1, 4.5, 5.7$ simulated LAEs, with a 1$σ$ deviation from the model in the $z=6.6$ catalog. This suggests that LAEs show higher order clustering terms similar to those of normal low redshift galaxies. △ Less

Submitted 6 November, 2020; originally announced November 2020.

Comments: 22 pages, 6 figures

Observables of spheroidal magnetized Strange Stars

Authors: Diana Alvear Terrero, Samantha López Pérez, Daryel Manreza Paret, Aurora Pérez Martínez, Gretel Quintero Angulo

Abstract: We study stable spheroidal configurations of magnetized Strange Stars using an axially symmetric metric in spherical coordinates that uses a gamma parameter to link the anisotropy in the Equation of State due to the magnetic field with the deformation of the star. The stars are composed by magnetized Strange Quark Matter described within the framework of the MIT-Bag model. Their masses, radii, ecc… ▽ More We study stable spheroidal configurations of magnetized Strange Stars using an axially symmetric metric in spherical coordinates that uses a gamma parameter to link the anisotropy in the Equation of State due to the magnetic field with the deformation of the star. The stars are composed by magnetized Strange Quark Matter described within the framework of the MIT-Bag model. Their masses, radii, eccentricity, redshift and mass quadrupole moment are computed. Results are compared with spherical Strange Stars solutions obtained with TOV equations and observational data of Strange Stars candidates. In the spheroidal model the observables depend directly on the deformation of the stars, and even though it is small, the observables strongly deviate from the corresponding spherical configurations. Thus, the highest values of the mass quadrupole moment correspond to the intermediate mass regime. These differences might allow to discriminate between models with/without magnetic field when compared with observations. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 20 pages, 16 figures

Journal ref: Phys. Rev. C 103, 045807 (2021)

Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup

Authors: Jane Huang, Sean M. Andrews, Karin I. Öberg, Megan Ansdell, Myriam Benisty, John M. Carpenter, Andrea Isella, Laura M. Pérez, Luca Ricci, Jonathan P. Williams, David J. Wilner, Zhaohuan Zhu

Abstract: While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolu… ▽ More While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3''$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like'' structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup's peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup's extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems. △ Less

Submitted 31 July, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 32 pages, 13 figures, accepted by ApJ. Typos corrected. Data are available at https://zenodo.org/record/3928862

Annular substructures in the transition disks around LkCa 15 and J1610

Authors: S. Facchini, M. Benisty, J. Bae, R. Loomis, L. Perez, M. Ansdell, S. Mayama, P. Pinilla, R. Teague, A. Isella, A. Mann

Abstract: We present high resolution millimeter continuum ALMA observations of the disks around the T Tauri stars LkCa 15 and J1610. These disks host dust-depleted inner regions, possibly carved by massive planets, and are of prime interest to study the imprints of planet-disk interactions. While at moderate angular resolution they appear as a broad ring surrounding a cavity, the continuum emission resolves… ▽ More We present high resolution millimeter continuum ALMA observations of the disks around the T Tauri stars LkCa 15 and J1610. These disks host dust-depleted inner regions, possibly carved by massive planets, and are of prime interest to study the imprints of planet-disk interactions. While at moderate angular resolution they appear as a broad ring surrounding a cavity, the continuum emission resolves into multiple rings at a resolution of ~60$\times$40 mas (~7.5 au for LkCa 15, ~6 au for J1610) and ~$7\,μ$Jy beam$^{-1}$ rms at 1.3 mm. In addition to a broad extended component, LkCa 15 and J1610 host 3 and 2 narrow rings, respectively, with two bright rings in LkCa 15 being radially resolved. The rings look marginally optically thick, with peak optical depths of ~0.5 (neglecting scattering), in agreement with high angular resolution observations of full disks. We perform hydrodynamical simulations with an embedded, sub-Jovian-mass planet and show that the observed multi-ringed substructure can be qualitatively explained as the outcome of the planet-disk interaction. We note however that the choice of the disk cooling timescale alone can significantly impact the resulting gas and dust distributions around the planet, leading to different numbers of rings and gaps and different spacings between them. We propose that the massive outer disk regions of transition disks are favorable places for planetesimals and possibly second generation planet formation of objects with a lower mass than the planets carving the inner cavity (typically few $M_{\rm Jup}$), and that the annular substructures observed in LkCa 15 and J1610 may be indicative of planetary core formation within dust-rich pressure traps. Current observations are compatible with other mechanisms being at the origin of the observed substructures, in particular with narrow rings generated at the edge of the CO and N$_2$ snowlines. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 17 pages, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 639, A121 (2020)