Showing 1–50 of 55 results for author: Alarcón, F

Robust Melanoma Thickness Prediction via Deep Transfer Learning enhanced by XAI Techniques

Authors: Miguel Nogales, Begoña Acha, Fernando Alarcón, José Pereyra, Carmen Serrano

Abstract: This study focuses on analyzing dermoscopy images to determine the depth of melanomas, which is a critical factor in diagnosing and treating skin cancer. The Breslow depth, measured from the top of the granular layer to the deepest point of tumor invasion, serves as a crucial parameter for staging melanoma and guiding treatment decisions. This research aims to improve the prediction of the depth o… ▽ More This study focuses on analyzing dermoscopy images to determine the depth of melanomas, which is a critical factor in diagnosing and treating skin cancer. The Breslow depth, measured from the top of the granular layer to the deepest point of tumor invasion, serves as a crucial parameter for staging melanoma and guiding treatment decisions. This research aims to improve the prediction of the depth of melanoma through the use of machine learning models, specifically deep learning, while also providing an analysis of the possible existance of graduation in the images characteristics which correlates with the depth of the melanomas. Various datasets, including ISIC and private collections, were used, comprising a total of 1162 images. The datasets were combined and balanced to ensure robust model training. The study utilized pre-trained Convolutional Neural Networks (CNNs). Results indicated that the models achieved significant improvements over previous methods. Additionally, the study conducted a correlation analysis between model's predictions and actual melanoma thickness, revealing a moderate correlation that improves with higher thickness values. Explainability methods such as feature visualization through Principal Component Analysis (PCA) demonstrated the capability of deep features to distinguish between different depths of melanoma, providing insight into the data distribution and model behavior. In summary, this research presents a dual contribution: enhancing the state-of-the-art classification results through advanced training techniques and offering a detailed analysis of the data and model behavior to better understand the relationship between dermoscopy images and melanoma thickness. △ Less

Submitted 19 June, 2024; originally announced June 2024.

PDRs4All IX. Sulfur elemental abundance in the Orion Bar

Authors: Asunción Fuente, Evelyne Roueff, Franck Le Petit, Jacques Le Bourlot, Emeric Bron, Mark G. Wolfire, James F. Babb, Pei-Gen Yan, Takashi Onaka, John H. Black, Ilane Schroetter, Dries Van De Putte, Ameek Sidhu, Amélie Canin, Boris Trahin, Felipe Alarcón, Ryan Chown, Olga Kannavou, Olivier Berné, Emilie Habart, Els Peeters, Javier R. Goicoechea, Marion Zannese, Raphael Meshaka, Yoko Okada , et al. (9 additional authors not shown)

Abstract: One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the Ja… ▽ More One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the James Webb Space Telescope provides an unprecedented opportunity to estimate the sulfur abundance through the observation of the [S I] 25.249 $μ$m line. We used the [S III] 18.7 $μ$m, [S IV] 10.5 $μ$m, and [S l] 25.249 $μ$m lines to estimate the amount of sulfur in the ionized and molecular gas along the Orion Bar. For the theoretical part, we used an upgraded version of the Meudon photodissociation region (PDR) code to model the observations. New inelastic collision rates of neutral atomic sulfur with ortho- and para- molecular hydrogen were calculated to predict the line intensities. The [S III] 18.7 $μ$m and [S IV] 10.5 $μ$m lines are detected over the imaged region with a shallow increase (by a factor of 4) toward the HII region. We estimate a moderate sulfur depletion, by a factor of $\sim$2, in the ionized gas. The corrugated interface between the molecular and atomic phases gives rise to several edge-on dissociation fronts we refer to as DF1, DF2, and DF3. The [S l] 25.249 $μ$m line is only detected toward DF2 and DF3, the dissociation fronts located farthest from the HII region. The detailed modeling of DF3 using the Meudon PDR code shows that the emission of the [S l] 25.249 $μ$m line is coming from warm ($>$ 40 K) molecular gas located at A$_{\rm V}$ $\sim$ 1$-$5 mag from the ionization front. Moreover, the intensity of the [S l] 25.249 $μ$m line is only accounted for if we assume the presence of undepleted sulfur. △ Less

Submitted 4 June, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 16 pages, 6 figures. Accepted for publication in Astronomy and Astrophysics

Thermal Structure and Millimeter Emission of Protoplanetary Disk with embedded protoplanets from radiative transfer modeling

Authors: Felipe Alarcón, Edwin Bergin

Abstract: The discovery of protoplanets and circumplanetary disks provides a unique opportunity to characterize planet formation through observations. Massive protoplanets shape the physical and chemical structure of their host circumstellar disk by accretion, localized emission, and disk depletion. In this work, we study the thermal changes induced within the disk by protoplanet accretion and synthetic pre… ▽ More The discovery of protoplanets and circumplanetary disks provides a unique opportunity to characterize planet formation through observations. Massive protoplanets shape the physical and chemical structure of their host circumstellar disk by accretion, localized emission, and disk depletion. In this work, we study the thermal changes induced within the disk by protoplanet accretion and synthetic predictions through hydrodynamical simulations with post-processed radiative transfer with an emphasis on radio millimeter emission. We explored distinct growth conditions and varied both planetary accretion rates and the local dust-to-gas mass ratios for a protoplanet at 1200 K. The radiative transfer models show that beyond the effect of disk gaps, in most cases, the CPD and the planet's emission locally increase the disk temperature. Moreover, depending on the local dust-to-gas depletion and accretion rate, the CPD presence may have detectable signatures in millimeter emission. It also has the power to generate azimuthal asymmetries important for continuum subtraction. Thus, if other means of detection of protoplanets are proven, the lack of corresponding evidence at other wavelengths can set limits on their growth timescales through a combined analysis of the local dust-to-gas ratio and the accretion rate. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 17 pages, 10 figures. Accepted for publication in ApJ

Extinction values towards embedded planets in protoplanetary disks estimated from hydrodynamic simulations

Authors: Felipe Alarcón, Edwin Bergin, Gabriele Cugno

Abstract: The upcoming new coronographs with deeper contrast limits, together with planned and current high-contrast imaging campaigns will push the detectability limit of protoplanets. These planet-hunting campaigns present a new opportunity to characterize protoplanets and their surrounding environments. However, there are clear uncertainties as to what are the extinction levels at different regions of pr… ▽ More The upcoming new coronographs with deeper contrast limits, together with planned and current high-contrast imaging campaigns will push the detectability limit of protoplanets. These planet-hunting campaigns present a new opportunity to characterize protoplanets and their surrounding environments. However, there are clear uncertainties as to what are the extinction levels at different regions of protoplanetary disks, which will impede our ability to characterize young planets. A correct understanding of the expected extinction together with multiple photometric observations will lead to constraints on the extinction levels, dust growth, disk evolution and protoplanetary accretion rates. In this work, we used hydrodynamic simulations and protoplanetary disk observational constraints obtained from both dust and gas emission to explore the expected extinction maps for continuum filters associated with strong hydrogen lines as tracers of accretion and key broadband photometric filters. We provide a scaling relationship for the extinction as a function of planetary separation and disk mass for three different gas giant masses. We also report values for a subset of disks of interest targetted by multiple imaging campaigns. The described values will be useful for the optimal design of future planet-hunting surveys and for giving context to non-detections in protoplanetary disks and the observed fluxes of point sources along with the birth conditions of protoplanets. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ. 18 pages, 16 figures, 3 tables

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

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

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

Submitted 3 April, 2024; originally announced April 2024.

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

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

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

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

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

Submitted 29 February, 2024; originally announced March 2024.

Journal ref: Science, 383, 6686, 2024

Formation of the Methyl Cation by Photochemistry in a Protoplanetary Disk

Authors: Olivier Berné, Marie-Aline Martin-Drumel, Ilane Schroetter, Javier R. Goicoechea, Ugo Jacovella, Bérenger Gans, Emmanuel Dartois, Laurent Coudert, Edwin Bergin, Felipe Alarcon, Jan Cami, Evelyne Roueff, John H. Black, Oskar Asvany, Emilie Habart, Els Peeters, Amelie Canin, Boris Trahin, Christine Joblin, Stephan Schlemmer, Sven Thorwirth, Jose Cernicharo, Maryvonne Gerin, Alexander Tielens, Marion Zannese , et al. (31 additional authors not shown)

Abstract: Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase orga… ▽ More Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation. △ Less

Submitted 6 January, 2024; originally announced January 2024.

Comments: Published in Nature

Journal ref: Nature 621, 56-59 (2023)

PDRs4All. V. Modelling the dust evolution across the illuminated edge of the Orion Bar

Authors: M. Elyajouri, N. Ysard, A. Abergel, E. Habart, L. Verstraete, A. Jones, M. Juvela, T. Schirmer, R. Meshaka, E. Dartois, J. Lebourlot, G. Rouille, T. Onaka, E. Peeters, O. Berne, F. Alarcon, J. Bernard-Salas, M. Buragohain, J. Cami, A. Canin, R. Chown, K. Demyk, K. Gordon, O. Kannavou, M. Kirsanova , et al. (9 additional authors not shown)

Abstract: We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advant… ▽ More We study the emission of dust grains within the Orion Bar - a well-known, highly far-UV (FUV)-irradiated PDR. The Orion Bar because of its edge-on geometry provides an exceptional benchmark for characterizing dust evolution and the associated driving processes under varying physical conditions. Our goal is to constrain the local properties of dust by comparing its emission to models. Taking advantage of the recent JWST PDRs4All data, we follow the dust emission as traced by JWST NIRCam (at 3.35 and 4.8 micron) and MIRI (at 7.7, 11.3, 15.0, and 25.5 micron), along with NIRSpec and MRS spectroscopic observations. First, we constrain the minimum size and hydrogen content of carbon nano-grains from a comparison between the observed dust emission spectra and the predictions of the THEMIS dust model coupled to the numerical code DustEM. Using this dust model, we then perform 3D radiative transfer simulations of dust emission with the SOC code and compare to data obtained along well chosen profiles across the Orion Bar. The JWST data allows us, for the first time, to spatially resolve the steep variation of dust emission at the illuminated edge of the Orion Bar PDR. By considering a dust model with carbonaceous nano-grains and submicronic coated silicate grains, we derive unprecedented constraints on the properties of across the Orion Bar. To explain the observed emission profiles with our simulations, we find that the nano-grains must be strongly depleted with an abundance (relative to the gas) 15 times less than in the diffuse ISM. The NIRSpec and MRS spectroscopic observations reveal variations in the hydrogenation of the carbon nano-grains. The lowest hydrogenation levels are found in the vicinity of the illuminating stars suggesting photo-processing while more hydrogenated nano-grains are found in the cold and dense molecular region, potentially indicative of larger grains. △ Less

Submitted 2 January, 2024; originally announced January 2024.

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

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

OH as a probe of the warm water cycle in planet-forming disks

Authors: Marion Zannese, Benoît Tabone, Emilie Habart, Javier R. Goicoechea, Alexandre Zanchet, Ewine F. van Dishoeck, Marc C. van Hemert, John H. Black, Alexander G. G. M. Tielens, A. Veselinova, P. G. Jambrina, M. Menendez, E. Verdasco, F. J. Aoiz, L. Gonzalez-Sanchez, Boris Trahin, Emmanuel Dartois, Olivier Berné, Els Peeters, Jinhua He, Ameek Sidhu, Ryan Chown, Ilane Schroetter, Dries Van De Putte, Amélie Canin , et al. (30 additional authors not shown)

Abstract: Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The hi… ▽ More Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The highly excited OH infrared rotational lines are the tell-tale signs of H2O destruction by FUV. The OH infrared ro-vibrational lines are attributed to chemical excitation via the key reaction O+H=OH+H which seeds the formation of water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per month and replenished. These results show that under warm and irradiated conditions water is destroyed and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce the HDO/H2O ratio in the warm regions of planet-forming disks. △ Less

Submitted 22 December, 2023; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: Version submitted to Nature Astronomy

Viscous heating as the dominant heat source inside the water snowline of V883 Ori

Authors: Felipe Alarcón, Simón Casassus, Wladimir lyra, Sebastián Pérez, Lucas Cieza

Abstract: FU Orionis-type objects(FUors) are embedded protostars that undergo episodes of high accretion, potentially indicating a widespread but poorly understood phase in the formation of low-mass stars. Gaining a better understanding of the influence exerted by these outbursts on the evolution of the surrounding protoplanetary disc may hold significant implications for the process of planet formation and… ▽ More FU Orionis-type objects(FUors) are embedded protostars that undergo episodes of high accretion, potentially indicating a widespread but poorly understood phase in the formation of low-mass stars. Gaining a better understanding of the influence exerted by these outbursts on the evolution of the surrounding protoplanetary disc may hold significant implications for the process of planet formation and the evolution of disc chemistry. The heating due to outbursts of high accretion in FUors pushes the snowlines of key volatiles farther out in the disc, so they become easier to observe and study. Among the known FUors, V883 Ori is of particular interest. V883 Ori was the first FUor to show indirect evidence of a resolvable snowline beyond 40 au. By introducing a radial-dependent model of this source including viscous heating, we show that active heating is needed to reproduce the steep thermal profile of dust in the inner disc of V883 Ori. Our disc modeling combines the effect of stellar irradiation and the influence on the disc shape caused by the outburst of accretion. The accuracy of our model is tested by comparing synthetic ALMA images with continuum observations of V883 Ori, showing that the model successfully reproduces the 1.3 mm emission of V883 Ori at high spatial resolution. Our final predictions underline the importance of viscous heating as a predominant heat source for this type of object, changing the physical conditions (shape and temperature) of the disc, and influencing its evolution. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted for publication in MNRAS

C I Traces the Disk Atmosphere in the IM Lup Protoplanetary Disk

Authors: Charles J. Law, Felipe Alarcón, L. Ilsedore Cleeves, Karin I. Öberg, Teresa Paneque-Carreño

Abstract: The central star and its energetic radiation fields play a vital role in setting the vertical and radial chemical structure of planet-forming disks. We present observations that, for the first time, clearly reveal the UV-irradiated surface of a protoplanetary disk. Specifically, we spatially resolve the atomic-to-molecular (C I-to-CO) transition in the IM Lup disk with ALMA archival observations o… ▽ More The central star and its energetic radiation fields play a vital role in setting the vertical and radial chemical structure of planet-forming disks. We present observations that, for the first time, clearly reveal the UV-irradiated surface of a protoplanetary disk. Specifically, we spatially resolve the atomic-to-molecular (C I-to-CO) transition in the IM Lup disk with ALMA archival observations of [C I] $^3$P$_1$-$^3$P$_0$. We derive a C I emitting height of z/r $\gtrsim$ 0.5 with emission detected out to a radius of ${\approx}$600 au. Compared to other systems with C I heights inferred from unresolved observations or models, the C I layer in the IM Lup disk is at scale heights almost double that of other disks, confirming its highly flared nature. C I arises from a narrow, optically-thin layer that is substantially more elevated than that of $^{12}$CO (z/r $\approx$ 0.3-0.4), which allows us to directly constrain the physical gas conditions across the C I-to-CO transition zone. We also compute a radially-resolved C I column density profile and find a disk-averaged C I column density of 2$\times10^{16}$ cm$^{-2}$, which is ${\approx}$3-20$\times$ lower than that of other disks with spatially-resolved C I detections. We do not find evidence for vertical substructures or spatially-localized deviations in C I due, e.g., to either an embedded giant planet or a photoevaporative wind that have been proposed in the IM Lup disk, but emphasize that deeper observations are required for robust constraints. △ Less

Submitted 13 December, 2023; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 12 pages, 6 figures, accepted for publication in ApJL; references added

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

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

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

Submitted 2 November, 2023; originally announced November 2023.

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

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

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

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

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

Submitted 12 October, 2023; originally announced October 2023.

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

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

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

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

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

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

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

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

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

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

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

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

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

Rheology of Pseudomonas fluorescens biofilms: from experiments to predictive DPD mesoscopic modelling

Authors: Jose Mart. n-Roca, Valentino Bianco, Francisco Alarcon, Ajay K. Monnappa, Paolo Natale, Francisco Monroy, Belen Orgaz, Ivan L. pez-Montero, Chantal Valeriani

Abstract: Bacterial biofilms mechanically behave as viscoelastic media consisting of micron-sized bacteria crosslinked to a selfproduced network of extracellular polymeric substances (EPS) embedded in water. Structural principles for numerical modelling aim at describing mesoscopic viscoelasticity without loosing detail on the underlying interactions existing in wide regimes of deformation under hydrodynami… ▽ More Bacterial biofilms mechanically behave as viscoelastic media consisting of micron-sized bacteria crosslinked to a selfproduced network of extracellular polymeric substances (EPS) embedded in water. Structural principles for numerical modelling aim at describing mesoscopic viscoelasticity without loosing detail on the underlying interactions existing in wide regimes of deformation under hydrodynamic stress. Here we approach the computational challenge to model bacterial biofilms for predictive mechanics in silico under variable stress conditions. Up-to-date models are not entirely satisfactory due to the plethora of parameters required to make them functioning under the effects of stress. As guided by the structural depiction gained in a previous work with Pseudomonas fluorescens (Jara et al. Front. Microbiol. (2021)), we propose a mechanical modeling by means of Dissipative Particle Dynamics (DPD), which captures the essentials of the topological and compositional interactions between bacteria particles and crosslinked EPS-embedding under imposed shear. The P. fluorescens biofilms have been modeled under mechanical stress mimicking shear stresses as undergone in vitro. The predictive capacity for mechanical features in DPD-simulated biofilms has been investigated by varying the externally imposed field of shear strain at variable amplitude and frequency. The parametric map of essential biofilm ingredients has been explored by making the rheological responses to emerge among conservative mesoscopic interactions and frictional dissipation in the underlying microscale. The proposed coarse grained DPD simulation qualitatively catches the rheology of the P. fluorescens biofilm over several decades of dynamic scaling. △ Less

Submitted 31 July, 2023; originally announced July 2023.

Journal ref: J. Chem. Phys. 158, 074902 (2023)

Molecules with ALMA at Planet-forming Scales (MAPS). Complex Kinematics in the AS 209 Disk Induced by a Forming Planet and Disk Winds

Authors: Maria Galloway-Sprietsma, Jaehan Bae, Richard Teague, Myriam Benisty, Stefano Facchini, Yuri Aikawa, Felipe Alarcón, Sean M. Andrews, Edwin Bergin, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzmán, Jane Huang, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Karin I. Öberg, Catherine Walsh, David J. Wilner

Abstract: We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a cohere… ▽ More We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets. △ Less

Submitted 12 May, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: This paper has been accepted for publication in the Astrophysical Journal (ApJ)

A localized kinematic structure detected in atomic carbon emission spatially coincident with a proposed protoplanet in the HD 163296 disk

Authors: Felipe Alarcón, Edwin Bergin, Richard Teague

Abstract: Over the last five years, studies of the kinematics in protoplanetary disks have led to the discovery of new protoplanet candidates and several structures linked to possible planet-disk interactions. We detect a localized kinematic bipolar structure in the HD 163296 disk present inside the deepest dust gap at 48 au from atomic carbon line emission. HD 163296's stellar jet and molecular winds have… ▽ More Over the last five years, studies of the kinematics in protoplanetary disks have led to the discovery of new protoplanet candidates and several structures linked to possible planet-disk interactions. We detect a localized kinematic bipolar structure in the HD 163296 disk present inside the deepest dust gap at 48 au from atomic carbon line emission. HD 163296's stellar jet and molecular winds have been described in detail in the literature; however, the kinematic anomaly in C I emission is not associated with either of them. Further, the velocity of the kinematic structure points indicates a component fast enough to differentiate it from the Keplerian profile of the disk; and its atomic nature hints at a localized UV source strong enough to dissociate CO and launch a C I outflow, or a strong polar flow from the upper layers of the disk. By discarding the stellar jet and previously observed molecular winds, we explore different sources for this kinematic feature in C I emission that could be associated with a protoplanet inflow/outflow, or disk winds. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in ApJL, 11 pages, 6 figures

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/

Rheology of Pseudomonas fluorescens biofilms: from experiments to DPD mesoscopic modelling

Authors: Jose Martin Roca, Valentino Bianco, Francisco Alarcon, Ajay K. Monnappa, Paolo Natale, Francisco Monroy, Belen Orgaz, Ivan Lopez-Montero, Chantal Valeriani

Abstract: The presence of bacterial biofilms in clinical and industrial settings is a major issue worldwide. A biofilm is a viscoelastic matrix, composed of bacteria producing a network of Extracellular Polymeric Substances (EPS) to which bacteria crosslink. Modelling of complex biofilms is relevant to provide accurate descriptions and predictions that include parameters as hydrodynamics, dynamics of the ba… ▽ More The presence of bacterial biofilms in clinical and industrial settings is a major issue worldwide. A biofilm is a viscoelastic matrix, composed of bacteria producing a network of Extracellular Polymeric Substances (EPS) to which bacteria crosslink. Modelling of complex biofilms is relevant to provide accurate descriptions and predictions that include parameters as hydrodynamics, dynamics of the bacterial population and solute mass transport. However, up-to-date numerical modelling, even at a coarse-grained level, is not satisfactorily. In this work, we present a numerical coarse-grain model of a bacterial biofilm, consisting of bacteria immersed in an aqueous matrix of a polymer network, that allows to study rheological properties of a biofilm. We study its viscoelastic modulus, varying topology and composition (such as the number of crosslinks between EPS polymers, the number of bacteria and the amount of solvent), and compare the numerical results with experimental rheological data of Pseudomonas fluorescens biofilms grown under static and shaking conditions, as previously described by Jara et al, Frontiers in Microbiology (2021). △ Less

Submitted 22 April, 2022; originally announced April 2022.

Intrinsic structure perspective for MIPS interfaces in two dimensional systems of Active Brownian Particles

Authors: Enrique Chacon, Francisco Alarcon, Jorge Ramirez, Pedro Tarazona, Chantal Valeriani

Abstract: Suspensions of Active Brownian Particles (ABP) undergo motility induced phase separation (MIPS) over a wide range of mean density and activity strength [1], even in the absence of an explicit attraction. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while the stable fluctuations of the interfacial line would be interpreted as a po… ▽ More Suspensions of Active Brownian Particles (ABP) undergo motility induced phase separation (MIPS) over a wide range of mean density and activity strength [1], even in the absence of an explicit attraction. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while the stable fluctuations of the interfacial line would be interpreted as a positive capillary surface tension [2], while in equilibrium liquid surfaces these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide a new insight in the origin of the MIPS from the local rectification of the random active force on the particles near the interface. As it had been pointed, that effect acts as an external potential [3] that produces a pressure gradient across the interface, so that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; but our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We test that with the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABP with purely repulsive (WCA-LJ model truncated at its minimum) and when the interaction includes a range of the LJ attractive force. △ Less

Submitted 4 February, 2022; originally announced February 2022.

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

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

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

Submitted 13 January, 2022; originally announced January 2022.

Comments: Submitted to PASP

Tests of General Relativity with GWTC-3

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, P. F. de Alarcón, S. Albanesi, R. A. Alfaidi, A. Allocca , et al. (1657 additional authors not shown)

Abstract: The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of th… ▽ More The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of those detectors. We restrict our analysis to the 15 confident signals that have false alarm rates $\leq 10^{-3}\, {\rm yr}^{-1}$. In addition to signals consistent with binary black hole (BH) mergers, the new events include GW200115_042309, a signal consistent with a neutron star--BH merger. We find the residual power, after subtracting the best fit waveform from the data for each event, to be consistent with the detector noise. Additionally, we find all the post-Newtonian deformation coefficients to be consistent with the predictions from GR, with an improvement by a factor of ~2 in the -1PN parameter. We also find that the spin-induced quadrupole moments of the binary BH constituents are consistent with those of Kerr BHs in GR. We find no evidence for dispersion of GWs, non-GR modes of polarization, or post-merger echoes in the events that were analyzed. We update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.27 \times 10^{-23} \mathrm{eV}/c^2$. The final mass and final spin as inferred from the pre-merger and post-merger parts of the waveform are consistent with each other. The studies of the properties of the remnant BHs, including deviations of the quasi-normal mode frequencies and damping times, show consistency with the predictions of GR. In addition to considering signals individually, we also combine results from the catalog of GW signals to calculate more precise population constraints. We find no evidence in support of physics beyond GR. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Report number: LIGO-P2100275

GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, S. Akcay, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin , et al. (1637 additional authors not shown)

Abstract: The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There ar… ▽ More The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There are 35 compact binary coalescence candidates identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} > 0.5$. Of these, 18 were previously reported as low-latency public alerts, and 17 are reported here for the first time. Based upon estimates for the component masses, our O3b candidates with $p_\mathrm{astro} > 0.5$ are consistent with gravitational-wave signals from binary black holes or neutron star-black hole binaries, and we identify none from binary neutron stars. However, from the gravitational-wave data alone, we are not able to measure matter effects that distinguish whether the binary components are neutron stars or black holes. The range of inferred component masses is similar to that found with previous catalogs, but the O3b candidates include the first confident observations of neutron star-black hole binaries. Including the 35 candidates from O3b in addition to those from GWTC-2.1, GWTC-3 contains 90 candidates found by our analysis with $p_\mathrm{astro} > 0.5$ across the first three observing runs. These observations of compact binary coalescences present an unprecedented view of the properties of black holes and neutron stars. △ Less

Submitted 23 October, 2023; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 88 pages (10 pages author list, 31 pages main text, 1 page acknowledgements, 24 pages appendices, 22 pages bibliography), 17 figures, 16 tables. Update to match version to be published in Physical Review X. Data products available from https://gwosc.org/GWTC-3/

Report number: LIGO-P2000318

Journal ref: Phys. Rev. X; 13(4):041039; 2023

Constraints on the cosmic expansion history from GWTC-3

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

Abstract: We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog.… ▽ More We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and $H(z)$. The source mass distribution displays a peak around $34\, {\rm M_\odot}$, followed by a drop-off. Assuming this mass scale does not evolve with redshift results in a $H(z)$ measurement, yielding $H_0=68^{+12}_{-7} {\rm km\,s^{-1}\,Mpc^{-1}}$ ($68\%$ credible interval) when combined with the $H_0$ measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the $H_0$ estimate from GWTC-1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of $H_0=68^{+8}_{-6} {\rm km\,s^{-1}\,Mpc^{-1}}$ with the galaxy catalog method, an improvement of 42% with respect to our GWTC-1 result and 20% with respect to recent $H_0$ studies using GWTC-2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about $H_0$) is the well-localized event GW190814. △ Less

Submitted 19 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Main paper: 30 pages, 15 figure, 7 tables

Report number: LIGO-P2100185-v6

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) VIII: CO Gap in AS 209--Gas Depletion or Chemical Processing?

Authors: Felipe Alarcón, Arthur Bosman, Edwin Bergin, Ke Zhang, Richard Teague, Jaehan Bae, Yuri Aikawa, Sean M. Andrews, Alice Booth, Jenny Calahan, Gianni Cataldi, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, François Ménard, Karin Öberg, Kamber R. Schwarz, Merel L. R. Van't Hoff, Catherine Walsh, David J. Wilner

Abstract: Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a… ▽ More Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C$_2$H column density traces changes in the C/O ratio rather than the H$_2$ gas surface density. However, the presence of a massive planet (> 0.2 M$_{Jup}$) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H$_2$ surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series. △ Less

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

Comments: 21 pages, 18 figures

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. XX. The Massive Disk Around GM Aurigae

Authors: Kamber R. Schwarz, Jenny K. Calahan, Ke Zhang, Felipe Alarcón, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, Enrique Macías, Melissa McClure, François Ménard, Karin I. Öberg, Richard Teague , et al. (3 additional authors not shown)

Abstract: Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk tempe… ▽ More Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (24-73 au at 159 pc) and the spatially unresolved HD J=1-0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 solar masses, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures

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). VII. Sub-stellar O/H and C/H and super-stellar C/O in planet feeding gas

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

Abstract: The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, wi… ▽ More The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, with variable C/O ratios and small-grain abundances, from the DALI code with CO and C2H column densities derived from the high-resolution observations of the disks of AS 209, HD 163296, and MWC 480. Elevated C/O ratios (~2.0), even within the CO ice line, are necessary to match the inferred C2H column densities, over most of the pebble disk. Combined with constraints on the CO abundances in these systems, this implies that both the O/H and C/H ratios in the gas are substellar by a factor of 4-10, with the O/H depleted by a factor of 20-50, resulting in the high C/O ratios. This necessitates that even within the CO ice line, most of the volatile carbon and oxygen is still trapped on grains in the midplane. Planets accreting gas in the gaps of the AS 209, HD 163296, and MWC 480 disks will thus acquire very little carbon and oxygen after reaching the pebble isolation mass. In the absence of atmosphere-enriching events, these planets would thus have a strongly substellar O/H and C/H and superstellar C/O atmospheric composition. △ Less

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

Comments: 19 pages, 8 figures This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Updates references for other MAPS papers

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) III: Characteristics of Radial Chemical Substructures

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

Abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We i… ▽ More The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 May, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 62 pages, 31 figures, accepted for publication in ApJS, MAPS cross-references updated, corrected Figure 21, updated gas disk sizes (Table 2, Figures 15-16) from associated Erratum

Molecules with ALMA at Planet-forming Scales (MAPS) XVII: Determining the 2D Thermal Structure of the HD 163296 Disk

Authors: Jenny K. Calahan, Edwin A. Bergin, Ke Zhang, Kamber R. Schwarz, Karin I. Oberg, Viviana V. Guzman, Catherine Walsh, Yuri Aikawa, Felipe Alarcon, Sean M. Andrews, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Feng Long, Ryan A. Loomis, Francois Menard, Hideko Nomura, Chunhua Qi , et al. (4 additional authors not shown)

Abstract: Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spa… ▽ More Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spatially resolved ALMA observations (~0.12 arcsec (13 au) - 0.25 arcsec (26 au)) of CO J = 2-1, 13CO J = 1-0, 2-1, C18O J = 1-0, 2-1, and C17O J = 1-0, the HD J = 1-0 flux upper limit, the spectral energy distribution (SED), and continuum morphology. The final model incorporates both a radial depletion of CO motivated by a time scale shorter than typical CO gas-phase chemistry (0.01 Myr) and an enhanced temperature near the surface layer of the the inner disk (z/r <= 0.21). This model agrees with the majority of the empirically derived temperatures and observed emitting surfaces derived from the J = 2-1 CO observations. We find an upper limit for the disk mass of 0.35 Msun, using the upper limit of the HD J = 1-0 and J = 2-1 flux. With our final thermal structure, we explore the impact that gaps have on the temperature structure constrained by observations of the resolved gaps. Adding a large gap in the gas and small dust additionally increases gas temperature in the gap by only 5-10%. 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: 15 pages + 11 pages of appendix, accepted to ApJS, part of MAPS collaboration

Study of the Parent-of-origin effect in monogenic diseases with variable age of onset. Application on ATTRv

Authors: Flora Alarcon, Violaine Planté-Bordeneuve, Gregory Nuel

Abstract: In genetic diseases with variable age of onset, an accurate estimation of the survival function for the mutation carriers and also modifying factors effects estimations are important for the management of asymptomatic gene carriers across life. Among the modifying factors, the gender of the parent transmitting the mutation (i.e. the parent-of-origin effect) has been shown to have a significant eff… ▽ More In genetic diseases with variable age of onset, an accurate estimation of the survival function for the mutation carriers and also modifying factors effects estimations are important for the management of asymptomatic gene carriers across life. Among the modifying factors, the gender of the parent transmitting the mutation (i.e. the parent-of-origin effect) has been shown to have a significant effect on survival curve estimation on transthyretin familial amyloid polyneuropathy (ATTRv) families. However, as most genotypes are unknown, the parent-of-origin must be calculated through a probability estimated from the pedigree. We propose in this article to extend the method providing mutation carrier survival estimates in order to estimate the parent-of-origin effect. The method is both validated on simulated data and applied to familly samples with ATTRv. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 14 pages, 4 figures,

Self-assembly of active amphiphilic Janus particles

Authors: S. A. Mallory, F. Alarcon, A. Cacciuto, C. Valeriani

Abstract: In this article, we study the phenomenology of a two dimensional dilute suspension of active amphiphilic Janus particles. We analyze how the morphology of the aggregates emerging from their self-assembly depends on the strength and the direction of the active forces. We systematically explore and contrast the phenomenologies resulting from particles with a range of attractive patch coverages. Fina… ▽ More In this article, we study the phenomenology of a two dimensional dilute suspension of active amphiphilic Janus particles. We analyze how the morphology of the aggregates emerging from their self-assembly depends on the strength and the direction of the active forces. We systematically explore and contrast the phenomenologies resulting from particles with a range of attractive patch coverages. Finally, we illustrate how the geometry of the colloids and the directionality of their interactions can be used to control the physical properties of the assembled active aggregates and suggest possible strategies to exploit self-propulsion as a tunable driving force for self-assembly. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Journal ref: New journal of Physics, 19 125014 (2017)

Destruction of refractory carbon grains drives the final stage of proto-planetary disk chemistry

Authors: Arthur D. Bosman, Felipe Alarcon, Ke Zhang, Edwin A. Bergin

Abstract: Here we aim to explore the origin of the strong C2H lines to reimagine the chemistry of protoplanetary disks. There are a few key aspects that drive our analysis. First, C2H is detected in young and old systems, hinting at a long-lived chemistry. Second, as a radical, C2H is rapidly destroyed, within <1000 yr. These two statements hint that the chemistry responsible for C2H emission must be predom… ▽ More Here we aim to explore the origin of the strong C2H lines to reimagine the chemistry of protoplanetary disks. There are a few key aspects that drive our analysis. First, C2H is detected in young and old systems, hinting at a long-lived chemistry. Second, as a radical, C2H is rapidly destroyed, within <1000 yr. These two statements hint that the chemistry responsible for C2H emission must be predominantly in the gas-phase and must be in equilibrium. Combining new and published chemical models we find that elevating the total volatile (gas and ice) C/O ratio is the only natural way to create a long lived, high C2H abundance. Most of the \ce{C2H} resides in gas with a Fuv/n-gas ~ 10^-7 G0 cm^3. To elevate the volatile C/O ratio, additional carbon has to be released into the gas to enable an equilibrium chemistry under oxygen-poor conditions. Photo-ablation of carbon-rich grains seems the most straightforward way to elevate the C/O ratio above 1.5, powering a long-lived equilibrium cycle. The regions at which the conditions are optimal for the presence of high C/O ratio and elevated C2H abundances in the gas disk set by the Fuv/n-gas ~ 10^-7 G0 cm^3 condition lie just outside the pebble disk as well as possibly in disk gaps. This process can thus also explain the (hints of) structure seen in C2H observations. △ Less

Submitted 29 January, 2021; originally announced January 2021.

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

Characterization of MIPS in a suspension of repulsive Active Brownian Particles through dynamical features

Authors: Jose Martin, Raul Martinez, Lachlan C. Alexander, Angel Luis Diez, Dirk G. A. L. Aarts, Francisco Alarcon, Jorge Ramirez, Chantal Valeriani

Abstract: The two-dimensional Active Brownian Particles system is meant to be composed of hard disks, that show excluded volume interactions, usually simulated via molecular dynamics using pure repulsive potentials. We show that the softness of the chosen potential plays a role in the result of the simulation, focusing on the case of the emergence of Motility Induced Phase Separation. In a pure hard-sphere… ▽ More The two-dimensional Active Brownian Particles system is meant to be composed of hard disks, that show excluded volume interactions, usually simulated via molecular dynamics using pure repulsive potentials. We show that the softness of the chosen potential plays a role in the result of the simulation, focusing on the case of the emergence of Motility Induced Phase Separation. In a pure hard-sphere system with no traslational diffusion,the phase diagram should be completely determined by their density and Péclet number. However, we have found two additional effects that affect the phase diagram in the ABP model we simulate: the relative strength of the traslational diffusion compared to the propulsion term and the overlapping of the particles. As we show, the second effect can be strongly mitigated if we use, instead of the standard Weeks-Chandler- Andersen potential, a harder one, the pseudo-hard spheres potential. Moreover, in determining the boundary of our phase space, we have tried different approaches to detect MIPS and concluded that observing dynamical features, via the non-Gaussian parameter, is more efficient than observing structural ones, via the local density distribution function. We also demonstrate that the Vogel-Fulcher equation successfully reproduces the decay of the diffusion as a function of density, except for very high density cases. Thus, the ABP system behaves similarly to a fragile glass in this regard. △ Less

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

Journal ref: Journal of Chemical Physics, 154, 164901 (2021)

Chemical Evolution in a Protoplanetary Disk within Planet Carved Gaps and Dust Rings

Authors: Felipe Alarcon, Richard Teague, Ke Zhang, Edwin Bergin, Marcelo Barraza-Alfaro

Abstract: Recent surveys of protoplanetary disks show that substructure in dust thermal continuum emission maps is common in protoplanetary disks. These substructures, most prominently rings and gaps, shape and change the chemical and physical conditions of the disk, along with the dust size distributions. In this work, we use a thermochemical code to focus on the chemical evolution that is occurring within… ▽ More Recent surveys of protoplanetary disks show that substructure in dust thermal continuum emission maps is common in protoplanetary disks. These substructures, most prominently rings and gaps, shape and change the chemical and physical conditions of the disk, along with the dust size distributions. In this work, we use a thermochemical code to focus on the chemical evolution that is occurring within the gas-depleted gap and the dust-rich ring often observed behind it. The composition of these spatial locations are of great import, as the gas and ice-coated grains will end up being part of the atmospheres of gas giants and/or the seeds of rocky planets. Our models show that the dust temperature at the midplane of the gap increases, enough to produce local sublimation of key volatiles and pushing the molecular layer closer to the midplane, while it decreases in the dust-rich ring, causing a higher volatile deposition onto the dust grain surfaces. Further, the ring itself presents a freeze-out trap for volatiles in local flows powered by forming planets, becoming a site of localized volatile enhancement. Within the gas depleted gap, the line emission depends on several different parameters, such as: the depth of the gap in surface density, the location of the dust substructure, and the abundance of common gas tracers, such as CO. In order to break this uncertainty between abundance and surface density, other methods such as disk kinematics, become necessary to constrain the disk structure and its chemical evolution. △ Less

Submitted 15 October, 2020; originally announced October 2020.

Comments: 21 pages,16 figures

Trapping flocking particles with asymmetric obstacles

Authors: Raul Martinez, Francisco Alarcon, Juan Luis Aragones, Chantal Valeriani

Abstract: Asymmetric obstacles can be exploited to direct the motion and induce sorting of run-and-tumbling particles. In this work, we show that flocking particles which follow the Vicsek model aligning rules experience a collective trapping in the presence of a wall of funnels made of chevrons, concentrating at the opposite side of a wall of funnels than run-and-tumbling particles. Flocking particles can… ▽ More Asymmetric obstacles can be exploited to direct the motion and induce sorting of run-and-tumbling particles. In this work, we show that flocking particles which follow the Vicsek model aligning rules experience a collective trapping in the presence of a wall of funnels made of chevrons, concentrating at the opposite side of a wall of funnels than run-and-tumbling particles. Flocking particles can be completely trapped or exhibit a dynamical trapping behaviour; these two regimes open the door to the design of a system with two perpendicular flows of active particles. This systematic study broaden our understanding about the emergence of collective motion of microorganisms in confined environments and direct the design of new microfluidics devices able to controlthese collective behaviours. △ Less

Submitted 29 November, 2019; originally announced November 2019.

Collective behavior of Vicsek particles without and with obstacles

Authors: Raul Martinez, Francisco Alarcon, Diego Rogel Rodriguez, Juan Luis Aragones, Chantal Valeriani

Abstract: In our work we have studied a two-dimensional suspension of finite-size Vicsek hard-disks, whose time evolution follows an event-driven dynamics between subsequent time steps. Having compared its collective behaviour with the one expected for a system of scalar Vicsek point-like particles, we have analysed the effect of considering two possible bouncing rules between the disks: a Vicsek-like rule… ▽ More In our work we have studied a two-dimensional suspension of finite-size Vicsek hard-disks, whose time evolution follows an event-driven dynamics between subsequent time steps. Having compared its collective behaviour with the one expected for a system of scalar Vicsek point-like particles, we have analysed the effect of considering two possible bouncing rules between the disks: a Vicsek-like rule and a pseudo-elastic one, focusing on the order-disorder transition. Next, we have added to the two-dimensional suspension of hard-disk Vicsek particles disk-like passive obstacles of two types: either fixed in space or moving according to the same event-driven dynamics. We have performed a detailed analysis of the particles' collective behaviour observed for both fixed and moving obstacles. In the fixed obstacles case, we have observed formation of clusters at low noise, in agreement with previous studies. When using moving passive obstacles, we found that that order of active particles is better destroyed as the drag of obstacles increases. In the no drag limit an interesting result was found: introduction of low drag passive particles can lead in some cases to a more ordered state of active flocking particles than what they show in bulk. △ Less

Submitted 23 November, 2019; originally announced November 2019.

Long Baseline Observations of the HD100546 Protoplanetary Disk with ALMA

Authors: Sebastián Pérez, Simon Casassus, Antonio Hales, Sebastián Marino, Anthony Cheetham, Alice Zurlo, Lucas Cieza, Ruobing Dong, Felipe Alarcón, Pablo Benítez-Llambay, Ed Fomalont, Henning Avenhaus

Abstract: Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or… ▽ More Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening. △ Less

Submitted 16 January, 2020; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 10 pages, 6 figures. Accepted for publication in ApJ Letters

Orientational order and morphology of clusters of self-assembled Janus swimmers

Authors: Francisco Alarcon, Eloy Navarro-Argemí, Chantal Valeriani, Ignacio Pagonabarraga

Abstract: Due to the combined effect of anisotropic interactions and activity, Janus swimmers are capable to self-assemble in a wide variety of structures, many more than their equilibrium counterpart. This might lead to the development of novel active materials capable of performing tasks without any central control. Their potential application in designing such materials endows trying to understand the fu… ▽ More Due to the combined effect of anisotropic interactions and activity, Janus swimmers are capable to self-assemble in a wide variety of structures, many more than their equilibrium counterpart. This might lead to the development of novel active materials capable of performing tasks without any central control. Their potential application in designing such materials endows trying to understand the fundamental mechanism in which these swimmers self-assemble. In the present work, we study a quasi-two dimensional semi-dilute suspensions of two classes of amphiphilic spherical swimmers whose direction of motion can be tuned: either swimmers propelling in the direction of the hydrophobic patch (WP), or swimmers propelling in the opposite direction (towards the hydrophilic side) (AP). In both systems we have systematically tuned swimmers' hydrophobic strength and signature, and observed that the anisotropic interactions, characterized by the angular attractive potential and its interaction range, in competition with the active stress, pointing towards or against the attractive patch gives rise to a rich aggregation phenomenology. △ Less

Submitted 21 May, 2019; originally announced May 2019.

Comments: Submitted in Physical Review E

Dust traps in the protoplanetary disc MWC 758: two vortices produced by two giant planets?

Authors: Clément Baruteau, Marcelo Barraza, Sebastián Pérez, Simon Casassus, Ruobing Dong, Wladimir Lyra, Sebastián Marino, Valentin Christiaens, Zhaohuan Zhu, Andrés Carmona, Florian Debras, Felipe Alarcon

Abstract: Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1… ▽ More Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario. △ Less

Submitted 15 March, 2019; originally announced March 2019.

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

Planet formation: The case for large efforts on the computational side

Authors: Wladimir Lyra, Thomas Haworth, Bertram Bitsch, Simon Casassus, Nicolás Cuello, Thayne Currie, Andras Gáspár, Hannah Jang-Condell, Hubert Klahr, Nathan Leigh, Giuseppe Lodato, Mordecai-Mark Mac Low, Sarah Maddison, George Mamatsashvili, Colin McNally, Andrea Isella, Sebastián Pérez, Luca Ricci, Debanjan Sengupta, Dimitris Stamatellos, Judit Szulágyi, Richard Teague, Neal Turner, Orkan Umurhan, Jacob White , et al. (32 additional authors not shown)

Abstract: Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been a… ▽ More Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 decadal survey

Fast quasi-conformal regional flattening of the left atrium

Authors: Marta Nuñez-Garcia, Gabriel Bernardino, Francisco Alarcón, Gala Caixal, Lluís Mont, Oscar Camara, Constantine Butakoff

Abstract: Two-dimensional representation of 3D anatomical structures is a simple and intuitive way for analysing patient information across populations and image modalities. It also allows convenient visualizations that can be included in clinical reports for a fast overview of the whole structure. While cardiac ventricles, especially the left ventricle, have an established standard representation (e.g. bul… ▽ More Two-dimensional representation of 3D anatomical structures is a simple and intuitive way for analysing patient information across populations and image modalities. It also allows convenient visualizations that can be included in clinical reports for a fast overview of the whole structure. While cardiac ventricles, especially the left ventricle, have an established standard representation (e.g. bull's eye plot), the 2D depiction of the left atrium (LA) is challenging due to its sub-structural complexity including the pulmonary veins (PV) and the left atrial appendage (LAA). Quasi-conformal flattening techniques, successfully applied to cardiac ventricles, require additional constraints in the case of the LA to place the PV and LAA in the same geometrical 2D location for different cases. Some registration-based methods have been proposed but 3D (or 2D) surface registration is time-consuming and prone to errors if the geometries are very different. We propose a novel atrial flattening methodology where a quasi-conformal 2D map of the LA is obtained quickly and without errors related to registration. In our approach, the LA is divided into 5 regions which are then mapped to their analogue two-dimensional regions. A dataset of 67 human left atria from magnetic resonance images (MRI) was studied to derive a population-based 2D LA template representing the averaged relative locations of the PVs and LAA. The clinical application of the proposed methodology is illustrated on different use cases including the integration of MRI and electroanatomical data. △ Less

Submitted 16 January, 2020; v1 submitted 16 November, 2018; originally announced November 2018.

Causal mediation analysis in presence of multiple mediators uncausally related

Authors: Allan Jerolon, Laura Baglietto, Etienne Birmele, Vittorio Perduca, Flora Alarcon

Abstract: Mediation analysis aims at disentangling the effects of a treatment on an outcome through alternative causal mechanisms and has become a popular practice in biomedical and social science applications. The causal framework based on counterfactuals is currently the standard approach to mediation, with important methodological advances introduced in the literature in the last decade, especially for s… ▽ More Mediation analysis aims at disentangling the effects of a treatment on an outcome through alternative causal mechanisms and has become a popular practice in biomedical and social science applications. The causal framework based on counterfactuals is currently the standard approach to mediation, with important methodological advances introduced in the literature in the last decade, especially for simple mediation, that is with one mediator at the time. Among a variety of alternative approaches, K. Imai et al. showed theoretical results and developed an R package to deal with simple mediation as well as with multiple mediation involving multiple mediators conditionally independent given the treatment and baseline covariates. This approach does not allow to consider the often encountered situation in which an unobserved common cause induces a spurious correlation between the mediators. In this context, which we refer to as mediation with uncausally related mediators, we show that, under appropriate hypothesis, the natural direct and joint indirect effects are non-parametrically identifiable. Moreover, we adopt the quasi-Bayesian algorithm developed by Imai et al. and propose a procedure based on the simulation of counterfactual distributions to estimate not only the direct and joint indirect effects but also the indirect effects through individual mediators. We study the properties of the proposed estimators through simulations. As an illustration, we apply our method on a real data set from a large cohort to assess the effect of hormone replacement treatment on breast cancer risk through three mediators, namely dense mammographic area, nondense area and body mass index. △ Less

Submitted 19 August, 2019; v1 submitted 21 September, 2018; originally announced September 2018.

Comments: 47 pages, 7 Tables, 9 Figures

Cm-wavelength observations of MWC758: resolved dust trapping in a vortex

Authors: Simon Casassus, Sebastian Marino, Wladimir Lyra, Clement Baruteau, Matias Vidal, Alwyn Wootten, Sebastian Perez, Felipe Alarcon, Marcelo Barraza, Miguel Carcamo, Ruobing Dong, Anibal Sierra, Zhaohuan Zhu, Luca Ricci, Valentin Christiaens, Lucas Cieza

Abstract: The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates fo… ▽ More The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains. △ Less

Submitted 30 November, 2018; v1 submitted 8 May, 2018; originally announced May 2018.

Comments: accepted to MNRAS

An inner warp in the DoAr 44 T Tauri transition disk

Authors: Simon Casassus, Henning Avenhaus, Sebastian Perez, Victor Navarro, Miguel Carcamo, Sebastian Marino, Lucas Cieza, Sascha P. Quanz, Felipe Alarcon, Alice Zurlo, Axel Osses, Fernando R. Rannou, Pablo E. Romaan, Marcelo Barraza

Abstract: Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring… ▽ More Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD142527 and HD100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr44. With a fairly axially symmetric ring in the submm radio continuum, DoAr44 is likely also a warped system. We constrain the warp geometry by comparing radiative transfer predictions with the DPI data in H band (Q_φ(H)) and with a re-processing of archival 336GHz ALMA observations. The observed DPI shadows have coincident radio counterparts, but the intensity drops are much deeper in Q_φ(H) (~88%), compared to the shallow drops at 336GHz (~24%). Radiative transfer predictions with an inner disk tilt of ~30+-5deg approximately account for the observations. ALMA long-baseline observations should allow the observation of the warped gas kinematics inside the cavity of DoAr44. △ Less

Submitted 6 April, 2018; originally announced April 2018.

Comments: 11 pages, 10 figures, 3 tables, accepted for publication in MNRAS