Showing 1–50 of 222 results for author: Wolf, S

Feasibility study on retrieving exoplanetary cloud cover distributions using polarimetry

Authors: S. Winning, M. Lietzow-Sinjen, S. Wolf

Abstract: Context. As a new growing field, exocartography aims to map the surface features of exoplanets that are beyond the resolution of traditional observing techniques. While photometric approaches have been discussed extensively, polarimetry has received less attention despite its promising prospects. Aims. We demonstrate that the limb polarization of an exoplanetary atmosphere offers valuable insigh… ▽ More Context. As a new growing field, exocartography aims to map the surface features of exoplanets that are beyond the resolution of traditional observing techniques. While photometric approaches have been discussed extensively, polarimetry has received less attention despite its promising prospects. Aims. We demonstrate that the limb polarization of an exoplanetary atmosphere offers valuable insights into its cloud cover distribution. Specifically, we determine an upper limit for the polarimetric precision, which is required to extract information about the latitudinal cloud cover of temperate Jovian planets for scenarios of observations with and without host stars. Methods. To compute the scattered stellar radiation of an exoplanetary atmosphere and to study the polarization at various planetary phase angles, we used the three-dimensional Monte Carlo radiative transfer code POLARIS. Results. When the planetary signal can be measured separately from the stellar radiation, information about the latitudinal cloud cover for polar cap models is accessible at polarimetric sensitivities of $0.1$ %. In contrast, a precision of about $10^{-3}$ ppm is required when the stellar flux is included to gain this information. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 8 pages, 5 figures

Results from the CsI Calorimeter onboard the 2023 ComPair Balloon Flight

Authors: Daniel Shy, Richard S. Woolf, Clio Sleator, Bernard Phlips, J. Eric Grove, Eric A. Wulf, Mary Johnson-Rambert, Mitch Davis, Emily Kong, Thomas Caligiure, A. Wilder Crosier, Aleksey Bolotnikov, Nicholas Cannady, Gabriella A. Carini, Regina Caputo, Jack Fried, Priyarshini Ghosh, Sean Griffin, Elizabeth Hays, Sven Herrmann, Carolyn Kierans, Nicholas Kirschner, Iker Liceaga-Indart, Zachary Metzler, Julie McEnery , et al. (11 additional authors not shown)

Abstract: The ComPair gamma-ray telescope is a technology demonstrator for a future gamma-ray telescope called the All-sky Medium Energy Gamma-ray Observatory (AMEGO). The instrument is composed of four subsystems, a double-sided silicon strip detector, a virtual Frisch grid CdZnTe calorimeter, a CsI:Tl based calorimeter, and an anti-coincidence detector (ACD). The CsI calorimeter's goal is to measure the p… ▽ More The ComPair gamma-ray telescope is a technology demonstrator for a future gamma-ray telescope called the All-sky Medium Energy Gamma-ray Observatory (AMEGO). The instrument is composed of four subsystems, a double-sided silicon strip detector, a virtual Frisch grid CdZnTe calorimeter, a CsI:Tl based calorimeter, and an anti-coincidence detector (ACD). The CsI calorimeter's goal is to measure the position and energy deposited from high-energy events. To demonstrate the technological readiness, the calorimeter has flown onboard a NASA scientific balloon as part of the GRAPE-ComPair mission and accumulated around 3 hours of float time at an altitude of 40 km. During the flight, the CsI calorimeter observed background radiation, Regener-Pfotzer Maximum, and several gamma-ray activation lines originating from aluminum. △ Less

Submitted 29 May, 2024; v1 submitted 10 May, 2024; originally announced May 2024.

Feasibility of detecting shadows in disks induced by infall

Authors: A. Krieger, M. Kuffmeier, S. Reissl, C. P. Dullemond, C. Ginski, S. Wolf

Abstract: Observations performed with high-resolution imaging techniques revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner with respect to an outer disk. The cause of misalignment, however, is still debated. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which inv… ▽ More Observations performed with high-resolution imaging techniques revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner with respect to an outer disk. The cause of misalignment, however, is still debated. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which involves the late infall of material onto a protostellar system. In particular, we use previously performed hydrodynamical simulations of such events, and generate flux maps in the visible, near-infrared, submillimeter, and millimeter wavelength range using Monte Carlo radiative transfer. Based on that, we derive synthetic observations of these systems performed with the instruments SPHERE/VLT and ALMA, which we use as a basis for our subsequent analysis. We find that near-infrared observations with SPHERE are particularly well suited for detecting shadows via direct imaging alongside other features such as gaps, arcs, and streamers. On the contrary, performing a shadow detection based on reconstructed ALMA observations is very challenging due to the high sensitivity that is required for this task. Thus, in cases that allow for a detection, sophisticated analyses may be needed, for instance by the utilization of carefully constructed azimuthal profiles, aiding the search for potentially shallow shadows. Lastly, we conclude that late infall-induced disk misalignment offers a plausible explanation for the emergence of shadows that are observed in various systems. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A, 16 pages, 16 figures, 1 table

Journal ref: A&A 686, A111 (2024)

Limitations of the modified blackbody fit method for determining molecular cloud properties

Authors: N. Zielinski, S. Wolf

Abstract: Achieving a comprehensive understanding of the star and planet formation process is one of the fundamental tasks of astrophysics, requiring detailed knowledge of the physical conditions during the different phases of this process. During the earliest stages, i.e., concerning physical processes in molecular clouds and filaments, the column density N(H2), dust temperature T and dust emissivity index… ▽ More Achieving a comprehensive understanding of the star and planet formation process is one of the fundamental tasks of astrophysics, requiring detailed knowledge of the physical conditions during the different phases of this process. During the earliest stages, i.e., concerning physical processes in molecular clouds and filaments, the column density N(H2), dust temperature T and dust emissivity index \b{eta} of these objects can be derived by adopting a modified blackbody fit of the far-infrared to (sub-)millimeter spectral energy distributions. However, this often applied method is based on various assumptions. In addition, the observational basis and required, but only assumed cloud properties, such as a limited wavelength-coverage of the spectral energy distribution and dust properties, respectively, may differ between different studies. We review the basic limitations of this method and evaluate their impact on the derived physical properties of the objects of interest, i.e., molecular clouds and filaments. We find that the highest uncertainty when applying this method is introduced by the often poorly constrained dust properties. Therefore, we propose to first derive the optical depth and subsequently the column density with the help of a suitable dust model as the optical depth can be obtained with high accuracy, especially at longer wavelengths. The method provides reliable results up to the high densities and corresponding optical depths observed in molecular clouds. Considering typically used observational data, i.e., measurements obtained with far-infrared instruments like Herschel/PACS, JCMT/SCUBA-2 and SOFIA/HAWC+, data at four wavelengths are sufficient to obtain accurate results. Furthermore, we find that the dust emissivity index \b{eta} derived with this method is not suitable as an indicator of dust grain size. △ Less

Submitted 28 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomische Nachrichten, 19 pages, 12 figures

Mid-infrared evidence for iron-rich dust in the multi-ringed inner disk of HD 144432

Authors: J. Varga, L. B. F. M. Waters, M. Hogerheijde, R. van Boekel, A. Matter, B. Lopez, K. Perraut, L. Chen, D. Nadella, S. Wolf, C. Dominik, Á. Kóspál, P. Ábrahám, J. -C. Augereau, P. Boley, G. Bourdarot, A. Caratti o Garatti, F. Cruz-Sáenz de Miera, W. C. Danchi, V. Gámez Rosas, Th. Henning, K. -H. Hofmann, M. Houllé, J. W. Isbell, W. Jaffe , et al. (18 additional authors not shown)

Abstract: Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interfe… ▽ More Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interferometric data taken by the Very Large Telescope Interferometer (VLTI), combining PIONIER, GRAVITY, and MATISSE observations. Methods. We introduced a new physical disk model, TGMdust, to image the interferometric data, and to fit the disk structure and dust composition. We also performed equilibrium condensation calculations with GGchem. Results. Our best-fit model has three disk zones with ring-like structures at 0.15, 1.3, and 4.1 au. Assuming that the dark regions in the disk at ~0.9 au and at ~3 au are gaps opened by planets, we estimate the masses of the putative gap-opening planets to be around a Jupiter mass. We find evidence for an optically thin emission ($τ<0.4$) from the inner two disk zones ($r<4$ au) at $λ>3\ μ$m. Our silicate compositional fits confirm radial mineralogy gradients. To identify the dust component responsible for the infrared continuum emission, we explore two cases for the dust composition, one with a silicate+iron mixture and the other with a silicate+carbon one. We find that the iron-rich model provides a better fit to the spectral energy distribution. Conclusions. We propose that in the warm inner regions ($r<5$ au) of typical planet-forming disks, most if not all carbon is in the gas phase, while iron and iron sulfide grains are major constituents of the solid mixture along with forsterite and enstatite. Our analysis demonstrates the need for detailed studies of the dust in inner disks with new mid-infrared instruments such as MATISSE and JWST/MIRI. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: 29 pages, 24 figures

Journal ref: A&A, 681, A47 (2024)

Impact of discontinuous grain size distributions on the spectral energy distribution of debris disks

Authors: Minjae Kim, Sebastian Wolf

Abstract: The collisional evolution of debris disks is expected to result in a characteristic wavy pattern of the grain size distributions, i.e., an under/overabundance of particles of specific sizes. This perturbed grain size distribution potentially leaves characteristic patterns in the spectral energy distribution (SED) of the disk system. We aim to quantify and understand the specific influence of disco… ▽ More The collisional evolution of debris disks is expected to result in a characteristic wavy pattern of the grain size distributions, i.e., an under/overabundance of particles of specific sizes. This perturbed grain size distribution potentially leaves characteristic patterns in the spectral energy distribution (SED) of the disk system. We aim to quantify and understand the specific influence of discontinuous particle size distributions on the appearance of debris disks. For this purpose, we consider dust emission models based on two different grain size distributions, i.e., once with a single and once with a broken power law. We compare the spectral index $α$ ($F_ν\,\propto ν^{\rm{\,α}}$) in the case of a continuous grain size distribution with that of a discontinuous grain size distribution. We perform this comparison for central stars with different spectral types and two different disk structures (e.g., slim and broad debris dust rings). Within the considered parameter space, we find a characteristic difference between the spectral slopes of the SED in the different scenarios. More specifically, the overabundance of small grains leads to a steeper slope in the far-infrared/submillimeter regime, while the spectral index in the mm regime is hardly affected. On the other hand, the underabundance of medium-sized grains results in a slight steepening of the far-infrared slope of SED, while its primary effect is on the mm slope of SED, causing it to become shallower. We also find that the impact of an overabundance of small dust particles is more pronounced than that of an underabundance of medium-sized dust particles. Furthermore, we find that the difference between the spectral indices for the two different grain size distributions is largest for debris disks around brighter central stars and broader disks. △ Less

Submitted 23 December, 2023; originally announced December 2023.

Comments: 14 pages, 10 figures

Improving Monte Carlo radiative transfer in the regime of high optical depths: The minimum scattering order

Authors: Anton Krieger, Sebastian Wolf

Abstract: Radiative transfer (RT) simulations are a powerful tool that enables the calculation of synthetic images of a wide range of astrophysical objects. These simulations are often based on the Monte Carlo (MC) method, as it provides the needed versatility that allows the consideration of the diverse and often complex conditions found in those objects. However, this method faces fundamental problems in… ▽ More Radiative transfer (RT) simulations are a powerful tool that enables the calculation of synthetic images of a wide range of astrophysical objects. These simulations are often based on the Monte Carlo (MC) method, as it provides the needed versatility that allows the consideration of the diverse and often complex conditions found in those objects. However, this method faces fundamental problems in the regime of high optical depths which may result in noisy images and underestimated flux values. In this study, we propose an advanced MCRT method, i.e., an enforced minimum scattering order that is aimed at providing a minimum quality of determined flux estimates. For that purpose, we extended our investigations of the scattering order problem and derived an analytic expression for the minimum number of interactions that depends on the albedo and optical depth of the system, which needs to be considered to achieve a certain coverage of the scattering order distribution. The method is based on the utilization of this estimated minimum scattering order and enforces the consideration of a sufficient number of interactions during a simulation. Moreover, we identified two notably distinct cases that shape the kind of complexity that arises in MCRT simulations: the albedo-dominated and the optical depth-dominated case. Based on that, we analyzed implications regarding the best usage of a stretching method as a means to alleviate the scattering order problem. We find that its most suitable application requires taking into account the albedo and the optical depth. Then, we argue that the derived minimum scattering order can be used to assess the performance of a stretching method with regard to the scattering orders its usage promotes. Finally, we stress the need for developing advanced pathfinding techniques to fully solve the problem of MCRT simulations in the regime of high optical depths. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A, 10 pages, 5 figures

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Improving Monte Carlo radiative transfer simulations: A shift of framework

Authors: Anton Krieger, Sebastian Wolf

Abstract: Monte Carlo radiative transfer (MCRT) simulations are a powerful tool for determining the appearance of astrophysical objects, analyzing the prevalent physical conditions within them, and inferring their properties on the basis of real observations. Consequently, a broad variety of codes has been implemented and optimized with the goal of solving this task efficiently. To that end, two distinct fr… ▽ More Monte Carlo radiative transfer (MCRT) simulations are a powerful tool for determining the appearance of astrophysical objects, analyzing the prevalent physical conditions within them, and inferring their properties on the basis of real observations. Consequently, a broad variety of codes has been implemented and optimized with the goal of solving this task efficiently. To that end, two distinct frameworks have emerged, namely, the extinction and the scattering framework, which form the basis of the path determination procedures of those codes. These procedures affect the step length of simulated photon packages and are used for determining flux estimates. Despite the fact that these simulations play an important role at present and thus require significant computational resources, little attention has been paid to the benefits and the drawbacks of both frameworks so far. In this study, we investigate their differences and assess their performance with regard to the quality of thereby obtained flux estimates, with a particular focus on the required computational demand. To that end, we use a testbed composed of an infinite plane-parallel slab, illuminated from one side, and we determine transmitted intensity using MCRT simulations for both frameworks. We find that there are vast differences between the frameworks with regard to their convergence speed. The scattering framework outperforms the extinction framework across all considered optical depths and albedos when solving this task, particularly in the regime of high optical depths. Its implementation can therefore greatly benefit all modern MCRT codes as it has the potential to significantly reduce required computation times. Thus, we highly recommend its consideration for various tasks that require MCRT simulations. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A, 10 pages, 8 figures

Journal ref: A&A 680, A67 (2023)

Shaping the CO snowline in protoplanetary disks

Authors: S. Gavino, J. Kobus, A. Dutrey, S. Guilloteau, S. Wolf, J. K. Jørgensen, R. Sharma

Abstract: Characterizing the dust thermal structure in protoplanetary disks is a fundamental task as the dust surface temperature can affect both the planetary formation and the chemical evolution. Since the temperature is dependent on many parameters, including the grain size, properly modeling the grain temperature structure can be challenging. Many chemistry disk models usually employ a sophisticated sin… ▽ More Characterizing the dust thermal structure in protoplanetary disks is a fundamental task as the dust surface temperature can affect both the planetary formation and the chemical evolution. Since the temperature is dependent on many parameters, including the grain size, properly modeling the grain temperature structure can be challenging. Many chemistry disk models usually employ a sophisticated single dust structure designed to reproduce the effect of a realistic population presumably composed of a large diversity of sizes. This generally represents a good approximation in most cases. Nonetheless, this dilutes the effects of the complex radiative interactions between the different grain populations on the resulting dust temperature, and thus the chemistry. We seek to show that the radiative interactions between dust grains of different sizes can induce a non-trivial dust temperature structure that cannot be reproduced by a single dust population and that can significantly affect the chemical outcome. The disk thermal structures are computed using the Monte-Carlo radiative transfer code RADMC-3D. The thermal structures are post-processed using the gas-grain code NAUTILUS to calculate the evolution of the chemical abundance. We find that simultaneously using at least two independent dust grain populations in disk models produces a complex temperature structure due to the starlight intercepted by the upper layers of the disk. In particular, we find that micron-sized dust grains are warmer than larger grains and can even show a radial temperature bump in some conditions. This dust temperature spread between the grains populations results in the segregation of the CO snowline and the presence of an unexpected CO gas hole along the midplane. We compare the results with observed close to edge-on class I/II disks. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 21 pages, 19 figures

Journal ref: A&A 680, A59 (2023)

Protoplanetary and debris disks in the $η$ Chamaeleontis Association: A sub-millimeter survey obtained with APEX/LABOCA

Authors: V. Roccatagliata, A. Sicilia-Aguilar, M. Kim, J. Campbell-White, M. Fang, S. J. Murphy, S. Wolf, W. A. Lawson, Th. Henning, J. Bouwman

Abstract: Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5-10,Myrs, where most of the disk should have already dissipated forming planets, are of particular interest. We present the first complete study of both protoplanetary and debris disks in a young region, using the $η$ Chamaeleontis association as a test bench to s… ▽ More Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5-10,Myrs, where most of the disk should have already dissipated forming planets, are of particular interest. We present the first complete study of both protoplanetary and debris disks in a young region, using the $η$ Chamaeleontis association as a test bench to study the cold disk content. We obtained sub-millimeter data for the entire core population down to late M-type stars, plus a few halo members. We performed a continuum sub-millimeter survey with APEX/LABOCA of all the core populations of $η$ Cha association. Disk properties have been derived by modeling protoplanetary and debris disks using RADMC 2D and DMS, respectively. We find that protoplanetary disks in $η$ Cha typically have holes with radii of the order of 0.01 to 0.03 AU, while ring-like emission from the debris disks is located between 20 au and 650 au from the central star. The parallaxes and Gaia eDR3 photometry, in combination with the PARSEC and COLIBRI isochrones, enable us to confirm an age of $η$ Cha between 7 and 9 Myrs. In general, the disk mass seems insufficient to support accretion over a long time, even for the lowest mass accretors, a clear difference compared with other regions and also a sign that the mass budget is further underestimated. We do not find a correlation between the stellar masses, accretion rates, and disk masses, although this could be due to sample issues. We confirm that the presence of inner holes is not enough to stop accretion unless accompanied by dramatic changes to the total disk mass content. Comparing $η$ Cha with other regions at different ages, we find that the physical processes responsible for debris disks (e.g., dust growth, dust trapping) efficiently act in less than 5 Myrs. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 20 pages, 9 figures, 7 tables

Impact of hot exozodiacal dust on the polarimetric analysis of close-in exoplanets

Authors: Kevin Ollmann, Sebastian Wolf, Moritz Lietzow, Thomas A. Stuber

Abstract: Hot exozodiacal dust (HEZD) found around main-sequence stars through interferometric observations in the photometric bands H to L is located close to the dust sublimation radius, potentially at orbital radii comparable to those of close-in exoplanets. Consequently, HEZD has a potential influence on the analysis of the scattered-light polarization of close-in exoplanets and vice versa. We analyze t… ▽ More Hot exozodiacal dust (HEZD) found around main-sequence stars through interferometric observations in the photometric bands H to L is located close to the dust sublimation radius, potentially at orbital radii comparable to those of close-in exoplanets. Consequently, HEZD has a potential influence on the analysis of the scattered-light polarization of close-in exoplanets and vice versa. We analyze the impact of HEZD on the polarimetric characterization of close-in exoplanets. This study is motivated in particular by the recently proven feasibility of exoplanet polarimetry. Applying the 3D Monte Carlo radiative transfer code POLARIS in an extended and optimized version for radiative transfer in exoplanetary atmospheres and an analytical tool for modeling the HEZD, we simulated and compared the polarization characteristics of the wavelength-dependent scattered-light polarization of HEZD and close-in exoplanets. The varied parameters are the planetary phase angle ($0^\circ-180^\circ$), the dust grain radius ($0.02\ μ$m $- \ 10\ μ$m), the HEZD mass ($10^{-10}$$\rm{M}_{\oplus}$ $-\ 10^{-8}$$\rm{M}_{\oplus}$), the orbital inclination ($0^\circ-90^\circ$), the composition of the planetary atmosphere (Mie and Rayleigh scattering atmosphere), the orbital radius of the HEZD ($0.02$ au $-\ 0.4$ au), and the planetary orbital radius ($0.01$ au $-\ 0.05$ au). The dust grain radius has the strongest influence on the polarimetric analysis due to its significant impact on the wavelength-dependent polarization characteristics and the total order of magnitude of the scattered-light polarization. In certain scenarios, the scattered-light polarization of the HEZD even exceeds that of the close-in exoplanet. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy \& Astrophysics. 14 pages, 13 figures

Journal ref: A&A 677, A187 (2023)

The potential of VLTI observations for the study of circumstellar disk variability

Authors: A. Bensberg, J. Kobus, S. Wolf

Abstract: Context. A characteristic feature of young stellar objects is their variability, which is caused by a variety of different physical processes. High-resolution interferometric observations in the near- and mid-infrared wavelength ranges spanning multiple epochs allow the detailed study of these processes. Aims. We aim at investigating the expected variations of the interferometric observables con… ▽ More Context. A characteristic feature of young stellar objects is their variability, which is caused by a variety of different physical processes. High-resolution interferometric observations in the near- and mid-infrared wavelength ranges spanning multiple epochs allow the detailed study of these processes. Aims. We aim at investigating the expected variations of the interferometric observables connected to changes in the measured photometric fluxes of a typical variable accreting central young stellar object with a circumstellar disk. Methods. We calculated visibilities and closure phases as well as the photometric flux of brightness distributions obtained using 3D Monte Carlo radiative transfer simulations for a model of a circumstellar disk with an accreting central star. Results. Changes in the accretion luminosity of the central object, that is, an accreting pre-main-sequence star, can lead to significant variations in the visibility and closure phase of the star-disk system measured with instruments at the Very Large Telescope Interferometer (VLTI) that can be related to changes in the photometric flux. Taking into account additional effects due to baseline variation, interferometric observations can provide valuable contributions to the understanding of the underlying processes. Additionally, we provide the web application VLTI B-VAR that allows the impact of the hour angle on the visibility and closure phase for customized intensity maps to be estimated. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 14 pages, 15 figures

Journal ref: A&A 677, A126 (2023)

How much large dust could be present in hot exozodiacal dust systems?

Authors: T. A. Stuber, F. Kirchschlager, T. D. Pearce, S. Ertel, A. V. Krivov, S. Wolf

Abstract: An infrared excess over the stellar photospheric emission of main-sequence stars has been found in interferometric surveys, commonly attributed to the presence of hot exozodiacal dust (HEZD). While submicrometer-sized grains in close vicinity to their host star have been inferred to be responsible for the found near-infrared excesses, the presence and amount of larger grains as part of the dust di… ▽ More An infrared excess over the stellar photospheric emission of main-sequence stars has been found in interferometric surveys, commonly attributed to the presence of hot exozodiacal dust (HEZD). While submicrometer-sized grains in close vicinity to their host star have been inferred to be responsible for the found near-infrared excesses, the presence and amount of larger grains as part of the dust distributions are weakly constrained. We quantify how many larger grains (above-micrometer-sized) could be present in addition to submicrometer-sized grains, while being consistent with observational constraints. This is important in order to distinguish between various scenarios for the origin of HEZD and to better estimate its observational appearance when observed with future instruments. We extended a model suitable to reproduce current observations of HEZD to investigate a bimodal size distribution. By deriving the characteristics of dust distributions whose observables are consistent with observational limits from interferometric measurements in the $K$ and $N$ bands we constrained the radii of sub- and above-micrometer-sized grains as well as their mass, number, and flux density ratios. In the most extreme cases of some of the investigated systems, large grains $\gtrsim 10\,μ$m might dominate the mass budget of HEZD while contributing up to 25$\,$% of the total flux density originating from the dust at a wavelength of 2.13$\,μ$m and up to 50$\,$% at a wavelength of 4.1$\,μ$m; at a wavelength of 11.1$\,μ$m their emission might clearly dominate over the emission of small grains. While it is not possible to detect such hot-dust distributions using ALMA, the ngVLA might allow us to detect HEZD at millimeter wavelengths. Large dust grains might have a more important impact on the observational appearance of HEZD than previously assumed, especially at longer wavelengths. △ Less

Submitted 20 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 18 pages, 7 figures

Journal ref: A&A 678, A121 (2023)

Development of Dual-Gain SiPM Boards for Extending the Energy Dynamic Range

Authors: Daniel Shy, Richard S. Woolf, Eric A. Wulf, Clio C. Sleator, Mary Johnson-Rambert, W. Neil Johnson, J. Eric Grove, Bernard F. Phlips

Abstract: Astronomical observations with gamma rays in the range of several hundred keV to hundreds of MeV currently represent the least explored energy range. To address this so-called MeV gap, we designed and built a prototype CsI:Tl calorimeter instrument using a commercial off-the-shelf (COTS) SiPMs and front-ends which may serve as a subsystem for a larger gamma-ray mission concept. During development,… ▽ More Astronomical observations with gamma rays in the range of several hundred keV to hundreds of MeV currently represent the least explored energy range. To address this so-called MeV gap, we designed and built a prototype CsI:Tl calorimeter instrument using a commercial off-the-shelf (COTS) SiPMs and front-ends which may serve as a subsystem for a larger gamma-ray mission concept. During development, we observed significant non-linearity in the energy response. Additionally, using the COTS readout, the calorimeter could not cover the four orders of magnitude in energy range required for the telescope. We, therefore, developed dual-gain silicon photomultiplier (SiPM) boards that make use of two SiPM species that are read out separately to increase the dynamic energy range of the readout. In this work, we investigate the SiPM's response with regards to active area ($3\times3 \ \mathrm{mm}^2$ and $1 \times 1 \ \mathrm{mm}^2$) and various microcell sizes ($10$, $20$, and $35 \ μ\mathrm{m}$). We read out $3\times3\times6 \ \mathrm{cm}^3$ CsI:Tl chunks using dual-gain SiPMs that utilize $35 \ μ\mathrm{m}$ microcells for both SiPM species and demonstrate the concept when tested with high-energy gamma-ray and proton beams. We also studied the response of $17 \times 17 \times 100 \ \mathrm{mm}^3$ CsI bars to high-energy protons. With the COTS readout, we demonstrate a sensitivity to $60 \ \mathrm{MeV}$ protons with the two SiPM species overlapping at a range of around $2.5-30 \ \mathrm{MeV}$. This development aims to demonstrate the concept for future scintillator-based high-energy calorimeters with applications in gamma-ray astrophysics. △ Less

Submitted 31 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Development of a CsI Calorimeter for the Compton-Pair (ComPair) Balloon-Borne Gamma-Ray Telescope

Authors: Daniel Shy, Richard S. Woolf, Clio C. Sleator, Eric A. Wulf, Mary Johnson-Rambert, Emily Kong, J. Mitch Davis, Thomas J. Caligiure, J. Eric Grove, Bernard F. Phlips

Abstract: There is a growing interest in astrophysics to fill in the observational gamma-ray MeV gap. We, therefore, developed a CsI:Tl calorimeter prototype as a subsystem to a balloon-based Compton and Pair-production telescope known as ComPair. ComPair is a technology demonstrator for a gamma-ray telescope in the MeV range that is comprised of 4 subsystems: the double-sided silicon detector, virtual Fris… ▽ More There is a growing interest in astrophysics to fill in the observational gamma-ray MeV gap. We, therefore, developed a CsI:Tl calorimeter prototype as a subsystem to a balloon-based Compton and Pair-production telescope known as ComPair. ComPair is a technology demonstrator for a gamma-ray telescope in the MeV range that is comprised of 4 subsystems: the double-sided silicon detector, virtual Frisch grid CdZnTe, CsI calorimeter, and a plastic-based anti-coincidence detector. The prototype CsI calorimeter is composed of thirty CsI logs, each with a geometry of $1.67 \times 1.67 \times 10 \ \mathrm{cm^3}$. The logs are arranged in a hodoscopic fashion with 6 in a row that alternate directions in each layer. Each log has a resolution of around $8 \%$ full-width-at-half-maximum (FWHM) at $662 \ \mathrm{keV}$ with a dynamic energy range of around $250\ \mathrm{keV}-30 \ \mathrm{MeV}$. A $2\times2$ array of SensL J-series SiPMs read out each end of the log to estimate the depth of interaction and energy deposition with signals read out with an IDEAS ROSSPAD. We also utilize an Arduino to synchronize with the other ComPair subsystems that comprise the full telescope. This work presents the development and performance of the calorimeter, its testing in thermal and vacuum conditions, and results from irradiation by $2-25 \ \mathrm{MeV}$ monoenergetic gamma-ray beams. The CsI calorimeter will fly onboard ComPair as a balloon experiment in the summer of 2023. △ Less

Submitted 31 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Scattered polarized radiation of extrasolar circumplanetary rings

Authors: M. Lietzow, S. Wolf

Abstract: We have investigated the impact of circumplanetary rings consisting of spherical micrometer-sized particles on the net scattered light polarization of extrasolar gas giants. Using the three-dimensional Monte Carlo radiative transfer code POLARIS, we studied the impact of the macroscopic parameters that define the ring, such as its radius and inclination, and the chemical composition of the ring pa… ▽ More We have investigated the impact of circumplanetary rings consisting of spherical micrometer-sized particles on the net scattered light polarization of extrasolar gas giants. Using the three-dimensional Monte Carlo radiative transfer code POLARIS, we studied the impact of the macroscopic parameters that define the ring, such as its radius and inclination, and the chemical composition of the ring particles on the net scattered polarization. For the spherical ring particles, we applied the Mie scattering theory. We studied the flux and polarization of the scattered stellar radiation as a function of planetary phase angle and wavelength from the optical to the near-infrared. For the chosen grain size distribution, the dust particles in the ring show strong forward scattering at the considered wavelengths. Thus, the reflected flux of the planet dominates the total reflected and polarized flux at small phase angles. However, the scattered and polarized flux of the ring increase at large phase angles and exceeds the total reflected planetary flux. For large rings that contain silicate particles, the total reflected flux is dominated by the radiation scattered by the dust in the ring at all phase angles. As a result, the orientation of polarization is parallel to the scattering plane at small phase angles. In contrast, for a ring that contains water ice particles, the orientation of polarization is parallel to the scattering plane at large phase angles. Depending on the ring inclination and orientation, the total reflected and polarized flux show a specific distribution as well. Large particles show a strong polarization at large phase angles compared to smaller particles. For a Jupiter-like atmosphere that contains methane and aerosols, methane absorption features are missing in the spectrum of a ringed planet. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 12 pages, 12 figures

Journal ref: A&A 671, A113 (2023)

Sub-mm/mm optical properties of real protoplanetary matter derived from Rosetta/MIRO observations of comet 67P

Authors: Johanna Bürger, Thilo Glißmann, Anthony Lethuillier, Dorothea Bischoff, Bastian Gundlach, Harald Mutschke, Sonja Höfer, Sebastian Wolf, Jürgen Blum

Abstract: Optical properties are required for the correct understanding and modelling of protoplanetary and debris discs. By assuming that comets are the most pristine bodies in the solar system, our goal is to derive optical constants of real protoplanetary material. We determine the complex index of refraction of the near-surface material of comet 67P/Churyumov-Gerasimenko by fitting the sub-millimetre/mi… ▽ More Optical properties are required for the correct understanding and modelling of protoplanetary and debris discs. By assuming that comets are the most pristine bodies in the solar system, our goal is to derive optical constants of real protoplanetary material. We determine the complex index of refraction of the near-surface material of comet 67P/Churyumov-Gerasimenko by fitting the sub-millimetre/millimetre observations of the thermal emission of the comet's sub-surface made by the Microwave Instrument for the Rosetta Orbiter (MIRO) with synthetic temperatures derived from a thermophysical model and radiative-transfer models. According to the two major formation scenarios of comets, we model the sub-surface layers to consist of pebbles as well as of homogeneously packed dust grains. In the case of a homogeneous dusty surface material, we find a solution for the length-absorption coefficient of $α\approx 0.22~\mathrm{cm^{-1}}$ for a wavelength of 1.594 mm and $α\geq 3.84~\mathrm{cm^{-1}}$ for a wavelength of 0.533 mm and a constant thermal conductivity of $0.006~\mathrm{Wm^{-1}K^{-1}}$. For the pebble scenario, we find for the pebbles and a wavelength of 1.594 mm a complex refractive index of $n = (1.074 - 1.256) + \mathrm{i} \, (2.580 - 7.431)\cdot 10^{-3}$ for pebble radii between 1 mm and 6 mm. Taking into account other constraints, our results point towards a pebble makeup of the cometary sub-surface with pebble radii between 3 mm and 6 mm. The derived real part of the refractive index is used to constrain the composition of the pebbles and their volume filling factor. The optical and physical properties are discussed in the context of protoplanetary and debris disc observations. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

Time-dependent Monte Carlo continuum radiative transfer

Authors: A. Bensberg, S. Wolf

Abstract: Aims. We present an implementation of an algorithm for 3D time-dependent Monte Carlo radiative transfer. It allows one to simulate temperature distributions as well as images and spectral energy distributions of the scattered light and thermal reemission radiation for variable illuminating and heating sources embedded in dust distributions, such as circumstellar disks and dust shells on time scale… ▽ More Aims. We present an implementation of an algorithm for 3D time-dependent Monte Carlo radiative transfer. It allows one to simulate temperature distributions as well as images and spectral energy distributions of the scattered light and thermal reemission radiation for variable illuminating and heating sources embedded in dust distributions, such as circumstellar disks and dust shells on time scales up to weeks. Methods. We extended the publicly available 3D Monte Carlo radiative transfer code POLARIS with efficient methods for the simulation of temperature distributions, scattering, and thermal reemission of dust distributions illuminated by temporally variable radiation sources. The influence of the chosen temporal step width and the number of photon packages per time step as key parameters for a given configuration is shown by simulating the temperature distribution in a spherical envelope around an embedded central star. The effect of the optical depth on the temperature simulation is discussed for the spherical envelope as well as for a model of a circumstellar disk with an embedded star. Finally, we present simulations of an outburst of a star surrounded by a circumstellar disk. Results. The presented algorithm for time-dependent 3D continuum Monte Carlo radiative transfer is a valuable basis for preparatory studies as well as for the analysis of continuum observations of the dusty environment around variable sources, such as accreting young stellar objects. In particular, the combined study of light echos in the optical and near-infrared wavelength range and the corresponding time-dependent thermal reemission observables of variable, for example outbursting sources, becomes possible on all involved spatial scales. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics. 15 pages, 18 figures

Journal ref: A&A 668, A120 (2022)

On the underestimation of dust mass in protoplanetary disks: Effects of disk structure and dust properties

Authors: Yao Liu, Hendrik Linz, Min Fang, Thomas Henning, Sebastian Wolf, Mario Flock, Giovanni P. Rosotti, Hongchi Wang, Dafa Li

Abstract: The total amount of dust grains in protoplanetary disks is one of the key properties that characterize the potential for planet formation. With (sub-)millimeter flux measurements, literature studies usually derive the dust mass using an analytic form under the assumption of optically thin emission, which may lead to substantial underestimation. In this work, we conduct a parameter study with the g… ▽ More The total amount of dust grains in protoplanetary disks is one of the key properties that characterize the potential for planet formation. With (sub-)millimeter flux measurements, literature studies usually derive the dust mass using an analytic form under the assumption of optically thin emission, which may lead to substantial underestimation. In this work, we conduct a parameter study with the goal of investigating the effects of disk structure and dust properties on the underestimation through self-consistent radiative transfer models. Different dust models, scattering modes and approaches for dust settling are considered and compared. The influences of disk substructures, such as rings and crescents, on the mass derivation are investigated as well. The results indicate that the traditional analytic method can underestimate the mass by a factor of a few to hundreds, depending on the optical depth along the line of sight set mainly by the true dust mass, disk size and inclination. As an application, we perform a detailed radiative transfer modeling of the spectral energy distribution of DoAr 33, one of the observed DSHARP disks. When the DSHARP dust opacities are adopted, the most probable dust mass returned from the Bayesian analysis is roughly 7 times higher than the value given by the analytic calculation. Our study demonstrates that estimating disk dust masses from radiative transfer modeling is one solution for alleviating the problem of insufficient mass for planet formation raised in the ALMA era. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 16 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 668, A175 (2022)

Using debris disk observations to infer substellar companions orbiting within or outside a parent planetesimal belt

Authors: T. A. Stuber, T. Löhne, S. Wolf

Abstract: Aims. We analyze whether the effects of secular perturbations, originating from a substellar companion, on the dust dynamics in a debris disk can be investigated with spatially resolved observations. Methods. We numerically simulated the collisional evolution of narrow and eccentric cold planetesimal belts around a star of spectral type A3V that are secularly perturbed by a companion that orbits… ▽ More Aims. We analyze whether the effects of secular perturbations, originating from a substellar companion, on the dust dynamics in a debris disk can be investigated with spatially resolved observations. Methods. We numerically simulated the collisional evolution of narrow and eccentric cold planetesimal belts around a star of spectral type A3V that are secularly perturbed by a companion that orbits either closer to or farther from the star than the belt. Based on the resulting spatial dust distributions, we simulated spatially resolved maps of their surface brightness in the $K$, $N$, and $Q$ bands and at wavelengths of 70$μ$m and 1300$μ$m. Results. Assuming a nearby debris disk seen face-on, we find that the brightness distribution varies significantly with observing wavelength, for example between the $N$ and $Q$ band. This can be explained by the varying relative contribution of the emission of the smallest grains near the blowout limit. The orbits of both the small grains that form the halo and the large grains close to the parent belt precess due to the secular perturbations induced by a companion orbiting inward of the belt. The halo, being composed of older grains, trails the belt. The magnitude of the trailing decreases with increasing perturber mass and hence with increasing strength of the perturbation, a trend we recovered in synthetic maps of surface brightness by fitting ellipses to lines of constant brightness. Systems with an outer perturber do not show a uniform halo precession since the orbits of small grains are strongly altered. We identified features of the brightness distributions suitable for distinguishing between systems with a potentially detectable inner or outer perturber, especially with a combined observation with JWST/MIRI in the $Q$ band tracing small grain emission and with ALMA at mm wavelengths tracing the position of the parent planetesimal belt. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics. 25 pages, 21 figures

Journal ref: A&A 669, A3 (2023)

Hot exozodis: cometary supply without trapping is unlikely to be the mechanism

Authors: Tim D. Pearce, Florian Kirchschlager, Gaël Rouillé, Steve Ertel, Alexander Bensberg, Alexander V. Krivov, Mark Booth, Sebastian Wolf, Jean-Charles Augereau

Abstract: Excess near-infrared emission is detected around one fifth of main-sequence stars, but its nature is a mystery. These excesses are interpreted as thermal emission from populations of small, hot dust very close to their stars (`hot exozodis'), but such grains should rapidly sublimate or be blown out of the system. To date, no model has fully explained this phenomenon. One mechanism commonly suggest… ▽ More Excess near-infrared emission is detected around one fifth of main-sequence stars, but its nature is a mystery. These excesses are interpreted as thermal emission from populations of small, hot dust very close to their stars (`hot exozodis'), but such grains should rapidly sublimate or be blown out of the system. To date, no model has fully explained this phenomenon. One mechanism commonly suggested in the literature is cometary supply, where star-grazing comets deposit dust close to the star, replenishing losses from grain sublimation and blowout. However, we show that this mechanism alone is very unlikely to be responsible for hot exozodis. We model the trajectory and size evolution of dust grains released by star-grazing comets, to establish the dust and comet properties required to reproduce hot-exozodi observations. We find that cometary supply alone can only reproduce observations if dust ejecta has an extremely steep size distribution upon release, and the dust-deposition rate is extraordinarily high. These requirements strongly contradict our current understanding of cometary dust and planetary systems. Cometary supply is therefore unlikely to be solely responsible for hot exozodis, so may need to be combined with some dust-trapping mechanism (such as gas or magnetic trapping) if it is to reproduce observations. △ Less

Submitted 22 September, 2022; originally announced September 2022.

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

Locating dust and molecules in the inner circumstellar environment of R~Sculptoris with MATISSE

Authors: Julien Drevon, Florentin Millour, Pierre Cruzalèbes, Claudia Paladini, Josef Hron, A. Meilland, F. Allouche, K. -H. Hofmann, S. Lagarde, B. Lopez, A. Matter, R. Petrov, S. Robbe-Dubois, D. Schertl, M. Wittkowski, G. Zins, P. Ábrahám, P. Antonelli, U. Beckmann, P. Berio, F. Bettonvil, A. Glindemann, U. Graser, M. Heininger, Thomas Henning , et al. (27 additional authors not shown)

Abstract: AGB stars are one of the main sources of dust production in the Galaxy. However, it is not clear what this process looks like and where the dust is condensing in the circumstellar environment. By characterizing the location of the dust and the molecules in the close environment of an AGB star, we aim to achieve a better understanding the history of the dust formation process. We observed the carbo… ▽ More AGB stars are one of the main sources of dust production in the Galaxy. However, it is not clear what this process looks like and where the dust is condensing in the circumstellar environment. By characterizing the location of the dust and the molecules in the close environment of an AGB star, we aim to achieve a better understanding the history of the dust formation process. We observed the carbon star R Scl with the VLTI-MATISSE instrument in L- and N-bands. The high angular resolution of the VLTI observations, combined with a large uv-plane coverage allowed us to use image reconstruction methods. To constrain the dust and molecules' location, we used two different methods: MIRA image reconstruction and the 1D code RHAPSODY. We found evidence of C2H2 and HCN molecules between 1 and 3.4 Rstar which is much closer to the star than the location of the dust (between 3.8 and 17.0 Rstar). We also estimated a mass-loss rate of 1.2+-0.4x10-6 Msun per yr. In the meantime, we confirmed the previously published characteristics of a thin dust shell, composed of amorphous carbon (amC) and silicon carbide (SiC). However, no clear SiC feature has been detected in the MATISSE visibilities. This might be caused by molecular absorption that can affect the shape of the SiC band at 11.3 micron. The appearance of the molecular shells is in good agreement with predictions from dynamical atmosphere models. For the first time, we co-located dust and molecules in the environment of an AGB star. We confirm that the molecules are located closer to the star than the dust. The MIRA images unveil the presence of a clumpy environment in the fuzzy emission region beyond 4.0 Rstar. Furthermore, with the available dynamic range and angular resolution, we did not detect the presence of a binary companion. Additional observations combining MATISSE and SAM-VISIR instrument should enable this detection in future studies. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 19 pages, published in A&A

VapoRock: Thermodynamics of vaporized silicate melts for modeling volcanic outgassing and magma ocean atmospheres

Authors: Aaron S. Wolf, Noah Jäggi, Paolo A. Sossi, Dan J. Bower

Abstract: Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software "VapoRock" combines the MELTS liquid model (Ghiorso et al., 1995) with gas-species properties from multiple thermochemistry tables (e.g., Chase et al., 1998). VapoRock calculates the partial p… ▽ More Silicate vapors play a key role in planetary evolution, especially dominating early stages of rocky planet formation through outgassed magma ocean atmospheres. Our open-source thermodynamic modeling software "VapoRock" combines the MELTS liquid model (Ghiorso et al., 1995) with gas-species properties from multiple thermochemistry tables (e.g., Chase et al., 1998). VapoRock calculates the partial pressures of 34 gaseous species in equilibrium with magmatic liquid in the system Si-Mg-Fe-Al-Ca-Na-K-Ti-Cr-O at desired temperatures and oxygen fugacities (fO2, or partial pressure of O2). Comparison with experiments shows that pressures and melt-oxide activities (which vary over many orders of magnitude) are reproduced to within a factor of ~3, consistent with measurement uncertainties. We also benchmark the model against a wide selection of igneous rock compositions including bulk silicate Earth, predicting elemental vapor abundances that are comparable (Na, Ca, & Al) or more realistic (K, Si, Mg, Fe, & Ti) than those of the closed-source MAGMA code (with maximum deviations by factors of 10-300 for K & Si). Vapor abundances depend critically on the activities of liquid components. The MELTS model underpinning VapoRock was calibrated and extensively tested on natural igneous liquids. In contrast, MAGMA's liquid model assumes ideal mixtures of a limited set of chemically simplified pseudo-species, which only roughly approximates the non-ideal compositional interactions typical of many-component natural silicate melts. Finally, we explore how relative abundances of SiO and SiO2 provide a spectroscopically measurable proxy for oxygen fugacity in devolatilized exoplanetary atmospheres, potentially constraining fO2 in outgassed exoplanetary mantles. △ Less

Submitted 20 December, 2022; v1 submitted 19 August, 2022; originally announced August 2022.

Millimeter Gap Contrast as a Probe for Turbulence Level in Protoplanetary Disks

Authors: Yao Liu, Gesa H. -M. Bertrang, Mario Flock, Giovanni P. Rosotti, Ewine F. van Dishoeck, Yann Boehler, Stefano Facchini, Can Cui, Sebastian Wolf, Min Fang

Abstract: Turbulent motions are believed to regulate angular momentum transport and influence dust evolution in protoplanetary disks. Measuring the strength of turbulence is challenging through gas line observations because of the requirement for high spatial and spectral resolution data, and an exquisite determination of the temperature. In this work, taking the well-known HD 163296 disk as an example, we… ▽ More Turbulent motions are believed to regulate angular momentum transport and influence dust evolution in protoplanetary disks. Measuring the strength of turbulence is challenging through gas line observations because of the requirement for high spatial and spectral resolution data, and an exquisite determination of the temperature. In this work, taking the well-known HD 163296 disk as an example, we investigated the contrast of gaps identified in high angular resolution continuum images as a probe for the level of turbulence. With self-consistent radiative transfer models, we simultaneously analyzed the radial brightness profiles along the disk major and minor axes, and the azimuthal brightness profiles of the B67 and B100 rings. By fitting all the gap contrasts measured from these profiles, we constrained the gas-to-dust scale height ratio $Λ$ to be $3.0_{-0.8}^{+0.3}$, $1.2_{-0.1}^{+0.1}$ and ${\ge}\,6.5$ for the D48, B67 and B100 regions, respectively. The varying gas-to-dust scale height ratios indicate that the degree of dust settling changes with radius. The inferred values for $Λ$ translate into a turbulence level of $α_{\rm turb}\,{<}\,3\times10^{-3}$ in the D48 and B100 regions, which is consistent with previous upper limits set by gas line observations. However, turbulent motions in the B67 ring are strong with $α_{\rm turb}\,{\sim}1.2\,{\times}\,10^{-2}$. Due to the degeneracy between $Λ$ and the depth of dust surface density drops, the turbulence strength in the D86 gap region is not constrained. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 13 pages, 7 figures, accepted for publication in Science China Physics, Mechanics & Astronomy

Journal ref: Sci. China-Phys. Mech. Astron. 65, 129511 (2022)

Synthetic observations using POLARIS: an application to simulations of massive prestellar cores

Authors: Joaquin Zamponi, Andrea Giannetti, Stefano Bovino, Giovanni Sabatini, Dominik R. G. Schleicher, Bastian Körtgen, Stefan Reissl, Sebastian Wolf

Abstract: Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make d… ▽ More Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make detections in massive sources hard to achieve. We present an application of the radiative transfer code POLARIS, with the goal to test the observability of the ortho-H$_2$D$^+$ transition $1_{10}$-$1_{11}$ (372.42 GHz) using simulations of high-mass collapsing cores that include deuteration chemistry. We analyzed an early and a late stage of the collapse of a 60 M$_{\odot}$ core, testing different source distances. For all cases, we generated synthetic single-dish and interferometric observations and studied the differences in both techniques. The column densities we derive are comparable to values reported for similar sources. These estimates depend on the extent over which they are averaged, and sources with compact emission they can be highly affected by beam dilution. Combined ALMA-ACA observations improve in signal-to-noise ratio and lead to better column density estimates as compared to ALMA alone. We confirm the feasibility to study ortho-H$_2$D$^+$ emission up to distances of 7 kpc. We provide a proof-of-concept of our framework for synthetic observations and highlight its importance when comparing numerical simulations with real observations. This work also proves how relevant it is to combine single-dish and interferometric measurements to derive appropriate source column densities. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astrophysics & Space Science. arXiv admin note: substantial text overlap with arXiv:2009.00407

The Twisted Magnetic Field of the Protobinary L483

Authors: Erin G. Cox, Giles Novak, Sarah Sadavoy, Leslie W. Looney, Dennis Lee, Marc Berthoud, Tyler L. Bourke, Simon Coudé, Frankie Encalada, Laura M. Fissel, Rachel Harrison, Martin Houde, Zhi-Yun Li, Philip C. Myers, Kate Pattle, Fabio P. Santos, Ian W. Stephens, Hailin Wang, Sebastian Wolf

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

Submitted 1 June, 2022; originally announced June 2022.

Comments: Accepted for publication in ApJ

Journal ref: 2022 ApJ 932 34

The disk of FU Orionis viewed with MATISSE/VLTI: first interferometric observations in $L$ and $M$ bands

Authors: F. Lykou, P. Ábrahám, L. Chen, J. Varga, Á. Kóspál, A. Matter, M. Siwak, Zs. M. Szabó, Z. Zhu, H. B. Liu, B. Lopez, F. Allouche, J. -C. Augereau, P. Berio, P. Cruzalèbes, C. Dominik, Th. Henning, K. -H. Hofmann, M. Hogerheijde, W. J. Jaffe, E. Kokoulina, S. Lagarde, A. Meilland, F. Millour, E. Pantin , et al. (8 additional authors not shown)

Abstract: The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data i… ▽ More The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data indicate that the dusty passive disk is silicate-rich. Only the innermost region of said dusty disk is found to emit strongly in the $N$ band, and it is resolved at an angular size of $\sim5$ mas, which translates to a diameter of about 2 au. The observations therefore place stringent constraints for the outer radius of the inner accretion disk. Dust radiative transfer simulations with RADMC-3D provide adequate fits to the spectral energy distribution from the optical to the submillimeter and to the interferometric observables when opting for an accretion rate $\dot{M}\sim 2\times 10^{-5}\, M_\odot$ yr$^{-1}$ and assuming $M_*=0.6\, M_\odot$. Most importantly, the hot inner accretion disk's outer radius can be fixed at 0.3 au. The outer radius of the dusty disk is placed at 100 au, based on constraints from scattered-light images in the literature. The dust mass contained in the disk is $2.4\times10^{-4}\, M_\odot$, and for a typical gas-to-dust ratio of 100, the total mass in the disk is approximately 0.02 $M_\odot$. We did not find any evidence for a nearby companion in the current interferometric data, and we tentatively explored the case of disk misalignment. For the latter, our modeling results suggest that the disk orientation is similar to that found in previous imaging studies by ALMA. Should there be an asymmetry in the very compact, inner accretion disk, this might be resolved at even smaller spatial scales ($\leq1$ mas). △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 17 pages, 14 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 663, A86 (2022)

Magnetic fields and outflows in the large Bok globule CB 54

Authors: Kate Pattle, Shih-Ping Lai, Sarah Sadavoy, Simon Coudé, Sebastian Wolf, Ray Furuya, Woojin Kwon, Chang Won Lee, Niko Zielinski

Abstract: We have observed the large Bok globule CB 54 in 850$μ$m polarised light using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We find that the magnetic field in the periphery of the globule shows significant, ordered deviation from the mean field direction in the globule centre. This deviation appears to correspond with the extended but relatively weak $^{12}$CO outflow emanatin… ▽ More We have observed the large Bok globule CB 54 in 850$μ$m polarised light using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We find that the magnetic field in the periphery of the globule shows significant, ordered deviation from the mean field direction in the globule centre. This deviation appears to correspond with the extended but relatively weak $^{12}$CO outflow emanating from the Class 0 sources at the centre of the globule. Energetics analysis suggests that if the outflow is reshaping the magnetic field in the globule's periphery, then we can place an upper limit of $<27\,μ$G on the magnetic field strength in the globule's periphery. Comparison with archival Planck and CARMA measurements shows that the field in the centre of the globule is consistent over several orders of magnitude in size scale, and oriented parallel to the density structure in the region in projection. We thus hypothesise that while non-thermal motions in the region may be sub-Alfvénic, the magnetic field is subdominant to gravity over a wide range of size scales. Our results suggest that even a relatively weak outflow may be able to significantly reshape magnetic fields in star-forming regions on scales $> 0.1$ pc. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 11 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The relative orientation between local magnetic field and Galactic plane in low latitude dark clouds

Authors: Gulafsha B. Choudhury, Himadri S. Das, B. J. Medhi, J. C. Pandey, S. Wolf, T. K. Dhar, A. M. Mazarbhuiya

Abstract: In this work, we study the magnetic field morphology of selected star-forming clouds spread over the galactic latitude ($b$) range, $-10^\circ$ to $10^\circ$. The polarimetric observation of clouds CB24, CB27 and CB188 are conducted to study the magnetic field geometry of those clouds from ARIES, Manora Peak, Nainital, India. These observations are combined with those of 14 further low latitude cl… ▽ More In this work, we study the magnetic field morphology of selected star-forming clouds spread over the galactic latitude ($b$) range, $-10^\circ$ to $10^\circ$. The polarimetric observation of clouds CB24, CB27 and CB188 are conducted to study the magnetic field geometry of those clouds from ARIES, Manora Peak, Nainital, India. These observations are combined with those of 14 further low latitude clouds available in the literature. Analyzing the polarimetric data of 17 clouds, we find that the alignment between the envelope magnetic field ($θ_{B}^{env}$) and Galactic plane ($θ_{GP}$) of the low-latitude clouds varies with their galactic longitudes ($l$). We observe a strong correlation between the longitude (\textit{l}) and the offset ($θ_{off}=|θ_B^{env}-θ_{GP}|$) which shows that $θ_{B}^{env}$ is parallel to the Galactic plane (GP) when the clouds are situated in the region, $115^\circ<l<250^\circ$. However, $θ_{B}^{env}$ has its own local deflection irrespective of the orientation of $θ_{GP}$ when the clouds are at $l<100^\circ$ and $l>250^\circ$. To check the consistency of our results, the stellar polarization data available at Heiles (2000) catalogue are overlaid on DSS image of the clouds having mean polarization vector of field stars. The results are almost consistent with the Heiles data. The effect of turbulence of the cloud is also studied which may play an important role in causing the misalignment phenomenon observed between $θ_{B}^{env}$ and $θ_{GP}$. We have used \textit{Herschel} \textit{SPIRE} 500 $μm$ and \textit{SCUBA} 850 $μm$ dust continuum emission maps in our work to understand the density structure of the clouds. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 22 pages, 7 figures. This is the Accepted Manuscript version of an article accepted for publication in Research in Astronomy and Astrophysics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. This Accepted Manuscript is published under a CC BY licence. The Version of Record is available soon online

Polarimetric investigation of selected cloud compositions in exoplanetary atmospheres

Authors: Moritz Lietzow, Sebastian Wolf

Abstract: We investigated the impact of selected cloud condensates in exoplanetary atmospheres on the polarization of scattered stellar radiation. We considered a selection of 25 cloud condensates that are expected to be present in extrasolar planetary atmospheres. Using the three-dimensional Monte Carlo radiative transfer code POLARIS and assuming Mie scattering theory, we calculated and studied the net po… ▽ More We investigated the impact of selected cloud condensates in exoplanetary atmospheres on the polarization of scattered stellar radiation. We considered a selection of 25 cloud condensates that are expected to be present in extrasolar planetary atmospheres. Using the three-dimensional Monte Carlo radiative transfer code POLARIS and assuming Mie scattering theory, we calculated and studied the net polarization of scattered radiation as a function of planetary phase angle at optical to near-infrared wavelengths. In addition to the well-known characteristics in the state of polarization, such as the rainbow determined by the real part of the refractive index, the behavior of the underlying imaginary part of the refractive index causes an increase or decrease in the degree of polarization and a change of sign in the polarization at a characteristic wavelength. In contrast to Al$_2$O$_3$ and MgFeSiO$_4$, clouds composed of SiO, MnS, Na$_2$S, or ZnS produce a rapidly decreasing degree of polarization with increasing wavelength in the context of an exoplanetary atmosphere. Furthermore, the sign of the polarization changes at a wavelength of about 0.5 $μ$m to 0.6 $μ$m, depending on the specific cloud condensate. The resulting net polarization is mainly positive for cloud compositions with large imaginary parts of the refractive index, such as Fe, FeS, and FeO. In addition, for Fe and FeS clouds, the maximum degree of polarization at long wavelengths is shifted to larger phase angles than for FeO. We found that most of these cloud condensates are distinguishable from each other due to their unique wavelength-dependent complex refractive index. In particular, an increase or decrease of the net polarization as a function of wavelength and a change of sign in the polarization at specific wavelengths are important features for characterizing cloud compositions in exoplanetary atmospheres. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: Accepted for publication in Astronomy & Astrophysics, 12 pages, 6 figures

Journal ref: A&A 663, A55 (2022)

The Future of Gamma-Ray Experiments in the MeV-EeV Range

Authors: Kristi Engel, Jordan Goodman, Petra Huentemeyer, Carolyn Kierans, Tiffany R. Lewis, Michela Negro, Marcos Santander, David A. Williams, Alice Allen, Tsuguo Aramaki, Rafael Alves Batista, Mathieu Benoit, Peter Bloser, Jennifer Bohon, Aleksey E. Bolotnikov, Isabella Brewer, Michael S. Briggs, Chad Brisbois, J. Michael Burgess, Eric Burns, Regina Caputo, Gabriella A. Carini, S. Bradley Cenko, Eric Charles, Stefano Ciprini , et al. (74 additional authors not shown)

Abstract: Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are… ▽ More Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Feasibility of detecting and characterizing embedded low-mass giant planets in gaps in the VIS/NIR wavelength range

Authors: A. Krieger, S. Wolf

Abstract: High-contrast imaging in the visible and near-infrared (VIS/NIR) has revealed the presence of a plethora of substructures in circumstellar disks (CSDs). One of the most commonly observed substructures are concentric gaps that may hint at the presence of embedded forming planets. However, direct detections of them are rare, and thus ambiguity regarding the origin of most gap features remains. The a… ▽ More High-contrast imaging in the visible and near-infrared (VIS/NIR) has revealed the presence of a plethora of substructures in circumstellar disks (CSDs). One of the most commonly observed substructures are concentric gaps that may hint at the presence of embedded forming planets. However, direct detections of them are rare, and thus ambiguity regarding the origin of most gap features remains. The aim of this study is to investigate the capabilities of high-contrast VIS/NIR imaging of directly detecting and characterizing low-mass giant planets in gaps. To this end, a grid of models of protoplanetary disks was generated. The models include a central T Tauri star surrounded by a face-on CSD harboring an accreting planet, which itself is surrounded by a circumplanetary disk (CPD) and carves a gap. We use the Monte Carlo radiative transfer code Mol3D to generate temperature distributions and synthetic observations. Based on these simulations, we measured the impact the planet and its CPD have on contrast curves and quantified the impact of the observing wavelength and of five key parameters on the determined signal strength. Then, we applied a detection criterion on our results and assess the capabilities of SPHERE/VLT of detecting the embedded planets. We find that a part of the investigated parameter space includes detectable planets, and we elaborate on the implications of a non-detection. Furthermore, we analyze the potential loss of valuable information by the use of a too small coronagraphic mask. However, we find this outcome to be very unlikely. Finally, within the VIS/NIR we identify for each of the investigated basic properties of the planets and the disks the most promising observing wavelengths that enable us to distinguish between different underlying parameter values. We find, that the detectability and the characterization often benefit from different observing wavelengths. △ Less

Submitted 7 March, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A, 30 pages, 13 figures, 8 tables

Journal ref: A&A 662, A99 (2022)

Thermal imaging of dust hiding the black hole in the Active Galaxy NGC 1068

Authors: Violeta Gamez Rosas, Jacob W. Isbell, Walter Jaffe, Romain G. Petrov, James H. Leftley, Karl-Heinz Hofmann, Florentin Millour, Leonard Burtscher, Klaus Meisenheimer, Anthony Meilland, Laurens B. F. M. Waters, Bruno Lopez, Stephane Lagarde, Gerd Weigelt, Philippe Berio, Fatme Allouche, Sylvie Robbe-Dubois, Pierre Cruzalebes, Felix Bettonvil, Thomas Henning, Jean-Charles Augereau, Pierre Antonelli, Udo Beckmann, Roy van Boekel, Philippe Bendjoya , et al. (27 additional authors not shown)

Abstract: In the widely accepted 'Unified Model' solution of the classification puzzle of Active Galactic Nuclei, the orientation of a dusty accretion torus around the central black hole dominates their appearance. In 'type-1' systems, the bright nucleus is visible at the centre of a face-on torus. In 'type-2' systems the thick, nearly edge-on torus hides the central engine. Later studies suggested evolutio… ▽ More In the widely accepted 'Unified Model' solution of the classification puzzle of Active Galactic Nuclei, the orientation of a dusty accretion torus around the central black hole dominates their appearance. In 'type-1' systems, the bright nucleus is visible at the centre of a face-on torus. In 'type-2' systems the thick, nearly edge-on torus hides the central engine. Later studies suggested evolutionary effects and added dusty clumps and polar winds but left the basic picture intact. However, recent high-resolution images of the archetypal type-2 galaxy NGC 1068 suggested a more radical revision. They displayed a ring-like emission feature which the authors advocated to be hot dust surrounding the black hole at the radius where the radiation from the central engine evaporates the dust. That ring is too thin and too far tilted from edge-on to hide the central engine, and ad hoc foreground extinction is needed to explain the type-2 classification. These images quickly generated reinterpretations of the type 1-2 dichotomy. Here we present new multi-band mid-infrared images of NGC1068 that detail the dust temperature distribution and reaffirm the original model. Combined with radio data, our maps locate the central engine below the previously reported ring and obscured by a thick, nearly edge-on disk, as predicted by the Unified Model. We also identify emission from polar flows and absorbing dust that is mineralogically distinct from that towards the Milky Way centre. △ Less

Submitted 27 December, 2021; originally announced December 2021.

Comments: In press at Nature

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Authors: A. Chiavassa, K. Kravchenko, M. Montargès, F. Millour, A. Matter, B. Freytag, M. Wittkowski, V. Hocdé, P. Cruzalèbes, F. Allouche, B. Lopez, S. Lagarde, R. G. Petrov, A. Meilland, S. Robbe-Dubois, K. -H. Hofmann, G. Weigelt, P. Berio, P. Bendjoya, F. Bettonvil, A. Domiciano de Souza, M. Heininger, Th. Henning, J. W. Isbell, W. Jaffe , et al. (28 additional authors not shown)

Abstract: Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reco… ▽ More Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: Accepted for publication on Astronomy and Astrophysics

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

Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)

Authors: Elena Orlando, Eugenio Bottacini, Alexander Moiseev, Arash Bodaghee, Werner Collmar, Torsten Ensslin, Igor V. Moskalenko, Michela Negro, Stefano Profumo, Matthew G. Baring, Aleksey Bolotnikov, Nicholas Cannady, Gabriella A. Carini, Seth Digel, Isabelle A. Grenier, Alice K. Harding, Dieter Hartmann, Sven Herrmann, Matthew Kerr, Roman Krivonos, Philippe Laurent, Francesco Longo, Aldo Morselli, Makoto Sasaki, Peter Shawhan , et al. (11 additional authors not shown)

Abstract: The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel… ▽ More The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds on and improves past and currently proposed missions at such energies. We outline the motivations for combining a coded mask and a Compton telescope and we define the scientific goals of such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which enable identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena. △ Less

Submitted 22 April, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 16 pages. Contact authors: Elena Orlando, Eugenio Bottacini

The impact of the presence of water ice on the analysis of debris disk observations

Authors: Thomas A. Stuber, Sebastian Wolf

Abstract: The analysis of debris disk observations is often based on the assumption of a dust phase composed of compact spherical grains consisting of astronomical silicate. Instead, observations indicate the existence of water ice in debris disks. We quantify the impact of water ice as a potential grain constituent in debris disks on the disk parameter values estimated from photometric and spatially resolv… ▽ More The analysis of debris disk observations is often based on the assumption of a dust phase composed of compact spherical grains consisting of astronomical silicate. Instead, observations indicate the existence of water ice in debris disks. We quantify the impact of water ice as a potential grain constituent in debris disks on the disk parameter values estimated from photometric and spatially resolved observations in the mid- and far-infrared. We simulated photometric measurements and radial profiles of debris disks containing water ice and analyzed them by applying a disk model purely consisting of astronomical silicate. Subsequently, we quantified the deviations between the derived and the true parameter values. As stars in central positions we discuss a $β$ Pic sibling and main-sequence stars with spectral types ranging from A0 to K5. To simulate observable quantities we employed selected observational scenarios regarding the choice of wavelengths and instrument characteristics. For the $β$ Pic stellar model and ice fractions $\geq 50\ \%$ the derived inner disk radius is biased by ice sublimation toward higher values. However, the derived slope of the radial density profile is mostly unaffected. Along with an increasing ice fraction, the slope of the grain size distribution is overestimated by up to a median factor of $\sim 1.2$ for an ice fraction of $90\ \%$ while the total disk mass is underestimated by a factor of $\sim 0.4$. The reliability of the derived minimum grain size strongly depends on the spectral type of the central star. For an A0-type star the minimum grain size can be underestimated by a factor of $\sim 0.2$, while for solar-like stars it is overestimated by up to a factor of $\sim 4 - 5$. Neglecting radial profile measurements and using solely photometric measurements, the factor of overestimation increases for solar-like stars up to $\sim 7 - 14$. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Astronomy $\&$ Astrophysics. 12 pages, 4 figures

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

VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris

Authors: K. -H. Hofmann, A. Bensberg, D. Schertl, G. Weigelt, S. Wolf, A. Meilland, F. Millour, L. B. F. M. Waters, S. Kraus, K. Ohnaka, B. Lopez, R. G. Petrov, S. Lagarde, Ph. Berio, F. Allouche, S. Robbe-Dubois, W. Jaffe, Th. Henning, C. Paladini, M. Schöller, A. Mérand, A. Glindemann, U. Beckmann, M. Heininger, F. Bettonvil , et al. (36 additional authors not shown)

Abstract: Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the di… ▽ More Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: accepted for publication in A&A

Magnetic field structure of OMC-3 in the far infrared revealed by SOFIA/HAWC+

Authors: Niko Zielinski, Sebastian Wolf

Abstract: We report the SOFIA/HAWC+ band D (154$\,μ$m) and E (214$\,μ$m) polarimetric observations of the filamentary structure OMC-3 that is part of the Orion molecular cloud. The polarization pattern is uniform for both bands and parallel to the filament structure. The polarization degree decreases toward regions with high intensity for both bands, revealing a so called "polarization hole." We identified… ▽ More We report the SOFIA/HAWC+ band D (154$\,μ$m) and E (214$\,μ$m) polarimetric observations of the filamentary structure OMC-3 that is part of the Orion molecular cloud. The polarization pattern is uniform for both bands and parallel to the filament structure. The polarization degree decreases toward regions with high intensity for both bands, revealing a so called "polarization hole." We identified an optical depth effect in which polarized emission and extinction act as counteracting mechanisms as a potential contributor to this phenomenon. Assuming that the detected polarization is caused by the emission of magnetically aligned non-spherical dust grains, the inferred magnetic field is uniform and oriented perpendicular to the filament. The magnetic field strength derived from the polarization patterns at 154$\,μ$m and 214$\,μ$m amounts to 202$\,μ$G and 261$\,μ$G, respectively. The derived magnetic field direction is consistent with that derived from previous polarimetric observations in the far infrared and submillimeter (submm) wavelength range. Investigating the far-infrared polarization spectrum derived from the SOFIA/HAWC+ observations, we do not find a clear correlation between the polarization spectrum and cloud properties, namely, the column density, $N(H_2$), and temperature, $T$. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A, 13 pages, 13 figures

Journal ref: A&A 659, A22 (2022)

Predicting resolved galaxy properties from photometric images using convolutional neural networks

Authors: Tobias Buck, Steffen Wolf

Abstract: Multi-band images of galaxies reveal a huge amount of information about their morphology and structure. However, inferring properties of the underlying stellar populations such as age, metallicity or kinematics from those images is notoriously difficult. Traditionally such information is best extracted from expensive spectroscopic observations. Here we present the… ▽ More Multi-band images of galaxies reveal a huge amount of information about their morphology and structure. However, inferring properties of the underlying stellar populations such as age, metallicity or kinematics from those images is notoriously difficult. Traditionally such information is best extracted from expensive spectroscopic observations. Here we present the $Painting\, IntrinsiC\, Attributes\, onto\, SDSS\, Objects$ (PICASSSO) project and test the information content of photometric multi-band images of galaxies. We train a convolutional neural network on 27,558 galaxy image pairs to establish a connection between broad-band images and the underlying physical stellar and gaseous galaxy property maps. We test our machine learning (ML) algorithm with SDSS $ugriz$ mock images for which uncertainties and systematics are exactly known. We show that multi-band galaxy images contain enough information to reconstruct 2d maps of stellar mass, metallicity, age and gas mass, metallicity as well as star formation rate. We recover the true stellar properties on a pixel by pixel basis with only little scatter, $\lesssim20\%$ compared to $\sim50\%$ statistical uncertainty from traditional mass-to-light-ratio based methods. We further test for any systematics of our algorithm with image resolution, training sample size or wavelength coverage. We find that galaxy morphology alone constrains stellar properties to better than $\sim20\%$ thus highlighting the benefits of including morphology into the parameter estimation. The machine learning approach can predict maps of high resolution, only limited by the resolution of the input bands, thus achieving higher resolution than IFU observations. The network architecture and all code is publicly available on GitHub. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 14 pages main text, 9 figures, 3 pages appendix with 5 figures, submitted to MNRAS

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

Authors: B. Lopez, S. Lagarde, R. G. Petrov, W. Jaffe, P. Antonelli, F. Allouche, P. Berio, A. Matter, A. Meilland, F. Millour, S. Robbe-Dubois, Th. Henning, G. Weigelt, A. Glindemann, T. Agocs, Ch. Bailet, U. Beckmann, F. Bettonvil, R. van Boekel, P. Bourget, Y. Bresson, P. Bristow, P. Cruzalèbes, E. Eldswijk, Y. Fanteï Caujolle , et al. (128 additional authors not shown)

Abstract: Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have develope… ▽ More Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now. △ Less

Submitted 2 March, 2022; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: 24 pages, 26 figures, submitted to Astronomy & Astrophysics

Journal ref: A&A 659, A192 (2022)

Evolution of Mercury's Earliest Atmosphere

Authors: Noah Jäggi, Diana Gamborino, Dan J. Bower, Paolo A. Sossi, Aaron S. Wolf, Apurva V. Oza, Audrey Vorburger, André Galli, Peter Wurz

Abstract: MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, a… ▽ More MESSENGER observations suggest a magma ocean formed on proto-Mercury, during which evaporation of metals and outgassing of C- and H-bearing volatiles produced an early atmosphere. Atmospheric escape subsequently occurred by plasma heating, photoevaporation, Jeans escape, and photoionization. To quantify atmospheric loss, we combine constraints on the lifetime of surficial melt, melt composition, and atmospheric composition. Consideration of two initial Mercury sizes and four magma ocean compositions determine the atmospheric speciation at a given surface temperature. A coupled interior-atmosphere model determines the cooling rate and therefore the lifetime of surficial melt. Combining the melt lifetime and escape flux calculations provide estimates for the total mass loss from early Mercury. Loss rates by Jeans escape are negligible. Plasma heating and photoionization are limited by homopause diffusion rates of $\sim10^{6}$ kg/s. Loss by photoevaporation depends on the timing of Mercury formation and assumed heating efficiency and ranges from $\sim10^{6.6}$ to $\sim10^{9.6}$ kg/s. The material for photoevaporation is sourced from below the homopause and is therefore energy-limited rather than diffusion-limited. The timescale for efficient interior-atmosphere chemical exchange is less than ten thousand years. Therefore, escape processes only account for an equivalent loss of less than 2.3 km of crust ($0.3\%$ of Mercury's mass). Accordingly, $\leq0.02\%$ of the total mass of H$_2$O and Na is lost. Therefore, cumulative loss cannot significantly modify Mercury's bulk mantle composition during the magma ocean stage. Mercury's high core:mantle ratio and volatile-rich surface may instead reflect chemical variations in its building blocks resulting from its solar-proximal accretion environment. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 25 pages, 8 figures, accepted to be published in The Planetary Science Journal

Journal ref: Planet. Sci. J. 2 (2021) 230

The formation of planetary systems with SPICA

Authors: I. Kamp, M. Honda, H. Nomura, M. Audard, D. Fedele, L. B. F. M. Waters, Y. Aikawa, A. Banzatti, J. E. Bowey, M. Bradford, C. Dominik, K. Furuya, E. Habart, D. Ishihara, D. Johnstone, G. Kennedy, M. Kim, Q. Kral, S. P. Lai, B. Larsson, M. McClure, A. Miotello, M. Momose, T. Nakagawa, D. Naylor , et al. (16 additional authors not shown)

Abstract: In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared spac… ▽ More In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated) △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: accepted for publication in PASA

First MATISSE L-band observations of HD 179218. Is the inner 10 au region rich in carbon dust particles?

Authors: E. Kokoulina, A. Matter, B. Lopez, E. Pantin, N. Ysard, G. Weigelt, E. Habart, J. Varga, A. Jones, A. Meilland, E. Dartois, L. Klarmann, J. -C. Augereau, R. van Boekel, M. Hogerheijde, G. Yoffe, L. B. F. M. Waters, C. Dominik, W. Jaffe, F. Millour, Th. Henning, K. -H. Hofmann, D. Schertl, S. Lagarde, R. G. Petrov , et al. (36 additional authors not shown)

Abstract: Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as… ▽ More Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons (PAHs)), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry. △ Less

Submitted 29 July, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Journal ref: A&A 652, A61 (2021)

Impact of Size-dependent Grain Temperature on Gas-Grain Chemistry in Protoplanetary Disks: the case of low mass star disks

Authors: S. Gavino, A. Dutrey, V. Wakelam, S. Guilloteau, J. Kobus, S. Wolf, W. Iqbal, E. Di Folco, E. Chapillon, V. Piétu

Abstract: Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-… ▽ More Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-phase astrochemical code NAUTILUS. Photoprocesses are handled using frequency-dependent cross-sections, and a new method to account for self and mutual shielding. The multi-grain model outputs are compared to those of single-grain size models (0.1 $μ$m), with two different assumptions for their equivalent temperature. We find that the Langmuir-Hinshelwood (LH) mechanism at equilibrium temperature is not efficient to form H$_2$ at 3-4 scale heights ($H$), and adopt a parametric fit to a stochastic method to model H$_2$ formation instead. We find the molecular layer composition (1-3 $H$) to depend on the amount of remaining H atoms. Differences in molecular surface densities between single and multi-grain models are mostly due to what occurs above 1.5 $H$. At 100 au, models with colder grains produce H$_2$O and CH$_4$ ices in the midplane, and warmer ones produce more CO$_2$ ices, both allowing efficient depletion of C and O as soon as CO sticks on grain surfaces. Complex organic molecules (COMs) production is enhanced by the presence of warmer grains in the multi-grain models. Using a single grain model mimicking grain growth and dust settling fails to reproduce the complexity of gas-grain chemistry. Chemical models with a single grain size are sensitive to the adopted grain temperature, and cannot account for all expected effects. A spatial spread of the snowlines is expected to result from the ranges in grain temperature. The amplitude of the effects will depend on the dust disk mass. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 39 pages, 33 figures

Journal ref: A&A 654, A65 (2021)

Mid-infrared circumstellar emission of the long-period Cepheid l Carinae resolved with VLTI/MATISSE

Authors: V. Hocdé, N. Nardetto, A. Matter, E. Lagadec, A. Mérand, P. Cruzalèbes, A. Meilland, F. Millour, B. Lopez, P. Berio, G. Weigelt, R. Petrov, J. W. Isbell, W. Jaffe, P. Kervella, A. Glindemann, M. Schöller, F. Allouche, A. Gallenne, A. Domiciano de Souza, G. Niccolini, E. Kokoulina, J. Varga, S. Lagarde, J. -C. Augereau , et al. (129 additional authors not shown)

Abstract: The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. W… ▽ More The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

Self-scattering on large, porous grains in protoplanetary disks with dust settling

Authors: Robert Brunngräber, Sebastian Wolf

Abstract: Observations of protoplanetary disks (PPDs) at sub-mm wavelengths suggest that polarisation is caused by scattering of thermal radiation. Most of the dust models that are used to explain these observations have major drawbacks: They either use much smaller grains than expected from dust evolution models, or result in lower polarisations than observed. We investigate the effect of dust grain porosi… ▽ More Observations of protoplanetary disks (PPDs) at sub-mm wavelengths suggest that polarisation is caused by scattering of thermal radiation. Most of the dust models that are used to explain these observations have major drawbacks: They either use much smaller grains than expected from dust evolution models, or result in lower polarisations than observed. We investigate the effect of dust grain porosity on the polarisation due to scattering at sub-mm wavelengths arising from grains of up to mm sizes, as they are expected to be present in the midplane of PPDs. Using the effective medium theory, we calculated the optical properties of porous dust and used them to predict the behaviour of the scattering polarisation. Subsequently, radiative transfer simulations for PPDs with porous dust were performed to analyse the additional effect of the optical depth structure, and thus the effect of multiple scattering events and inhomogeneous temperature distributions on the polarisation. We find that porous dust with moderate filling factors of about 10\% increase the polarisation compared to compact grains. For higher grain porosities, that is, grains with a filling factor of 1\% or lower, the extinction opacity decreases, as does the optical depth of a disk with constant mass. Consequently, the unpolarised direct radiation dominates the scattered flux, and the polarisation drops rapidly. Even though the simulated polarisation is higher than in the case of compact grains, it is still below the typical observed values for face-on disks. However, the polarisation can be increased when crucial model assumptions derived from disk and dust evolution theories, for instance, dust settling and mm-sized dust grains, are dropped. In the case of inclined disks, however, our reference model is able to achieve polarisation degrees of about 1\%, and using higher disk masses, even higher than this. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: Accepted by A&A

Journal ref: A&A 648, A87 (2021)

Hunting for Dark Matter and New Physics with GECCO

Authors: Adam Coogan, Alexander Moiseev, Logan Morrison, Stefano Profumo, Matthew G. Baring, Aleksey Bolotnikov, Gabriella A. Carini, Sven C. Herrmann, Francesco Longo, Floyd W. Stecker, Alessandro Armando Vigliano, Richard S. Woolf

Abstract: We outline the science opportunities in the areas of searches for dark matter and new physics offered by a proposed future MeV gamma-ray telescope, the Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO). We point out that such an instrument would play a critical role in opening up a discovery window for particle dark matter with mass in the MeV or sub-MeV range, in disentanglin… ▽ More We outline the science opportunities in the areas of searches for dark matter and new physics offered by a proposed future MeV gamma-ray telescope, the Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO). We point out that such an instrument would play a critical role in opening up a discovery window for particle dark matter with mass in the MeV or sub-MeV range, in disentangling the origin of the mysterious 511 keV line emission in the Galactic Center region, and in potentially discovering Hawking evaporation from light primordial black holes. △ Less

Submitted 2 May, 2023; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: 27 pages, 10 figures, ~matches published version

Journal ref: Phys. Rev. D 107, 023022 (2023)

Three-dimensional continuum radiative transfer of polarized radiation in exoplanetary atmospheres

Authors: Moritz Lietzow, Sebastian Wolf, Robert Brunngräber

Abstract: Polarimetry is about to become a powerful tool for determining the atmospheric properties of exoplanets. To provide the basis for the interpretation of observational results and for predictive studies to guide future observations, sophisticated analysis tools are required. Our goal is to develop a radiative transfer tool that contains all the relevant continuum polarization mechanisms for the comp… ▽ More Polarimetry is about to become a powerful tool for determining the atmospheric properties of exoplanets. To provide the basis for the interpretation of observational results and for predictive studies to guide future observations, sophisticated analysis tools are required. Our goal is to develop a radiative transfer tool that contains all the relevant continuum polarization mechanisms for the comprehensive analysis of the polarized flux resulting from the scattering in the atmosphere of, on the surface of, and in the local planetary environment (e.g., planetary rings, exomoons) of extra-solar planets. Furthermore, our goal is to avoid common simplifications such as locally plane-parallel planetary atmospheres, the missing cross-talk between latitudinal and longitudinal regions, or the assumption of either a point-like star or plane-parallel illumination. As a platform for the newly developed numerical algorithms, we use the 3D Monte Carlo radiative transfer code POLARIS. The code is extended and optimized for the radiative transfer in exoplanetary atmospheres. We investigate the reflected flux and its degree of polarization for different phase angles for a homogeneous cloud-free atmosphere and an inhomogeneous cloudy atmosphere. To take advantage of the 3D radiative transfer and to demonstrate the potential of the code, the impact of an additional circumplanetary ring on the reflected polarized flux is studied. The presence of a circumplanetary ring consisting of small water-ice particles has a noticeable impact on the reflected polarized radiation. In particular, the reflected flux strongly increases at larger phase angles if the planetary orbit is seen edge-on because the considered particles tend to scatter forwards. In contrast, the degree of polarization decreases at these phase angles. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted for publication in A&A, 12 pages, 10 figures

Journal ref: A&A 645, A146 (2021)

Constraining the magnetic field properties of Bok globule B335 using SOFIA/HAWC+

Authors: Niko Zielinski, Sebastian Wolf, Robert Brunngräber

Abstract: Thanks to their well-defined shape and mostly isolated locations, Bok globules are suitable objects for studying the physics of low-mass star formation. To study the magnetic field of the prototypical Bok globule B335, we obtained a spatially resolved polarization map with SOFIA/HAWC+ at a wavelength of 214$\,μ$m. For the first time, these observations reveal that polarization holes in Bok globule… ▽ More Thanks to their well-defined shape and mostly isolated locations, Bok globules are suitable objects for studying the physics of low-mass star formation. To study the magnetic field of the prototypical Bok globule B335, we obtained a spatially resolved polarization map with SOFIA/HAWC+ at a wavelength of 214$\,μ$m. For the first time, these observations reveal that polarization holes in Bok globules, that is, the decrease in polarization degree towards their dense centers, also occur in the far-infrared wavelength regime. The observed polarization pattern is uniform with a mean polarization angle of 48$^\circ\pm $26$^\circ$ and a magnetic field strength of $\sim$ 142$\,μ$G. Moreover, we use complementary polarimetic data for B335 obtained at near-infrared to millimeter wavelengths to analyze and constrain the magnetic field across different scales. By applying the 3D Monte-Carlo radiative transfer code POLARIS (Reissl et al. 2016), we developed a model for the density and magnetic field structure as well as for the dust properties of this globule. We conclude that the column density towards the center of B335 is too low to cause the observed polarization hole in B335 via dichroic absorption (Brauer et al. 2016). Furthermore, we conclude that the effect of self-scattering has no significant impact on the observed polarization. Adopting dust-grain alignment via the radiative torque mechanism, a combination of the interstellar radiation field and the central star as radiation sources is consistent with the decrease in polarization degree at the outer regions of B335 ($\approx\,$10$^4\,$au from the core). However, the model fails to explain the low polarization degree within the inner 5000 au. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: Accepted for publication in A&A, 11 pages, 10 figures

Journal ref: A&A 645, A125 (2021)