ALBERT

Description: AU Mic is a young planetary system with a resolved debris disc showing signs of planet formation and two transiting warm Neptunes near mean-motion resonances. Here we analyse three transits of AU Mic b observed with the CHaracterising ExOPlanet Satellite (CHEOPS), supplemented with sector 1 and 27 Transiting Exoplanet Survey Satellite (TESS) photometry, and the All-Sky Automated Survey from the ground. The refined orbital period of AU Mic b is 8.462995 ± 0.000003 d, whereas the stellar rotational period is Prot = 4.8367 ± 0.0006 d. The two periods indicate a 7:4 spin–orbit commensurability at a precision of 0.1%. Therefore, all transits are observed in front of one of the four possible stellar central longitudes. This is strongly supported by the observation that the same complex star-spot pattern is seen in the second and third CHEOPS visits that were separated by four orbits (and seven stellar rotations). Using a bootstrap analysis we find that flares and star spots reduce the accuracy of transit parameters by up to 10% in the planet-to-star radius ratio and the accuracy on transit time by 3–4 min. Nevertheless, occulted stellar spot features independently confirm the presence of transit timing variations (TTVs) with an amplitude of at least 4 min. We find that the outer companion, AU Mic c, may cause the observed TTVs.

Description: The grand-design face-on spiral galaxy M 51 is an excellent laboratory for studying magnetic fields in galaxies. Due to wavelength-dependent Faraday depolarization, linearly polarized synchrotron emission at different radio frequencies yields a picture of the galaxy at different depths: observations in the L-band (1–2 GHz) probe the halo region, while at 4.85 GHz (C-band) and 8.35 GHz (X-band), the linearly polarized emission mostly emerges from the disk region of M 51. We present new observations of M 51 using the Karl G. Jansky Very Large Array at the intermediate frequency range of the S-band (2–4 GHz), where previously no high-resolution broadband polarization observations existed, to shed new light on the transition region between the disk and the halo. We present the S-band radio images of the distributions of the total intensity, polarized intensity, degree of polarization, and rotation measure (RM). The RM distribution in the S-band shows a fluctuating pattern without any apparent large-scale structure. We discuss a model of the depolarization of synchrotron radiation in a multi-layer magneto-ionic medium and compare the model predictions to the multi-frequency polarization data of M 51 between 1–8 GHz. The model makes distinct predictions of a two-layer (disk–halo) and three-layer (far-side halo “disk” near-side halo) system. Since the model predictions strongly differ within the wavelength range of the S-band, the new S-band data are essential for distinguishing between the different systems. A two-layer model of M 51 is preferred. The parameters of the model are adjusted to fit to the data of polarization fractions in a few selected regions. In three spiral arm regions, the turbulent field in the disk dominates with strengths between 18 μG and 24 μG, while the regular field strengths are 8 − 16 μG. In one inter-arm region, the regular field strength of 18 μG exceeds that of the turbulent field of 11 μG. The regular field strengths in the halo are 3 − 5 μG. The observed RMs in the disk-halo transition region are probably dominated by tangled regular fields, as predicted from models of evolving dynamos, and/or vertical fields, as predicted from numerical simulations of Parker instabilities or galactic winds. Both types of magnetic fields have frequent reversals on scales similar to or larger than the beam size (∼550 pc) that contribute to an increase of the RM dispersion and to distortions of any large-scale pattern of the regular field. Our study devises new ways of analyzing and interpreting broadband multi-frequency polarization data that will be applicable to future data from, for example, the Square Kilometre Array.

Description: Context. Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the isotopologue ratio in the circumstellar envelope (CSE) may deviate considerably from the stellar atmospheric value. Aims. In this paper, we aim to examine how the 12CO/13CO and H12CN/H13CN abundance ratios vary radially due to chemical reactions in the outflows of asymptotic giant branch (AGB) stars and the effect of excitation and optical depth on the resulting line intensity ratios. We study both carbon-rich (C-type) and oxygen-rich (O-type) CSEs. Methods. We performed chemical modeling to derive radial abundance distributions of our selected species in the CSEs over a wide range of mass-loss rates (10−8 〈 Ṁ 〈 10−4 M⊙ yr−1). We used these as input in a non-local thermodynamic equilibrium radiative transfer code to derive the line intensities of several ground-state rotational transitions. We also test the influence of stellar parameters, physical conditions in the outflows, the intensity of the interstellar radiation field, and the importance of considering the chemical networks in our model results. Results. We quantified deviations from the atmospheric value for typical outflows. We find that the circumstellar value of 12CO/13CO can deviate from its atmospheric value by up to 25–94% and 6–60% for C- and O-type CSEs, respectively, in radial ranges that depend on the mass-loss rate. We show that variations of the intensity of the interstellar radiation field and the gas kinetic temperature can significantly influence the CO isotopologue abundance ratio in the outer CSEs of both C-type and O-type. On the contrary, the H12CN/H13CN abundance ratio is stable throughout the CSEs for all tested mass-loss rates. The radiative transfer modeling shows that the integrated line intensity ratio I12CO/I13CO of different rotational transitions varies significantly for stars with mass-loss rates in the range from 10−7 to 10−6 M⊙ yr−1 due to combined chemical and excitation effects. In contrast, the excitation conditions for the HCN isotopologues are the same for both isotopologues. Conclusions. We demonstrate the importance of using the isotopologue abundance profiles from detailed chemical models as inputs to radiative transfer models in the interpretation of isotopologue observations. Previous studies of circumstellar CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN molecule as a tracer of C isotope ratios is affected by fewer complicating problems, but we note that the corrections for high optical depths are very large in the case of high-mass-loss-rate C-type CSEs; and in O-type CSEs the H13CN lines may be too weak to detect.

Description: Context. The magnetic field in spiral galaxies is known to have a large-scale spiral structure along the galactic disk and is observed as X-shaped in the halo of some galaxies. While the disk field can be well explained by dynamo action, the three-dimensional structure of the halo field and its physical nature are still unclear. Aims. As first steps towards understanding the halo fields, we want to clarify whether or not the observed X-shaped field is a wide-spread pattern in the halos of spiral galaxies. We also aim to investigate whether these halo fields are simply turbulent fields ordered by compression or shear (anisotropic turbulent fields), or have a large-scale regular structure. Methods. Analysis of the Faraday rotation in the halo is used as a tool to distinguish anisotropic turbulent fields from large-scale magnetic fields. However, this has been challenging until recently because of the faint halo emission in linear polarization. Our sensitive VLA broadband observations in C-band and L-band of 35 spiral galaxies seen edge-on (called CHANG-ES) allowed us to perform rotation measure synthesis (RM synthesis) in their halos and to analyze the results. We further accomplished a stacking of the observed polarization maps of 28 CHANG-ES galaxies in C-band. Results. Though the stacked edge-on galaxies were of different Hubble type, and had differing star formation activity and interaction activity, the stacked image clearly reveals an X-shaped structure of the apparent magnetic field. We detected a large-scale (coherent) halo field in all 16 galaxies that have extended polarized emission in their halos. We detected large-scale field reversals in all of their halos. In six galaxies, these are along lines that are approximately perpendicular to the galactic midplane (vertical RMTL) with about 2 kpc separation. Only in NGC 3044 and possibly in NGC 3448 did we observe vertical giant magnetic ropes (GMR) similar to those detected recently in NGC 4631. Conclusions. The observed X-shaped structure of the halo field seems to be an underlying feature of spiral galaxies. It can be regarded as the two-dimensional projection of the regular magnetic field which we found to have scales of typically 1 kpc or larger observed over several kiloparsecs. The ordered magnetic field extends far out in the halo and beyond. We detected large-scale magnetic field reversals in the halo that may indicate that GMR are more or less tightly wound. With these discoveries, we hope to stimulate model simulations for the halo magnetic field that should also explain the determined asymmetry of the polarized intensity (PI).

Description: Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the ∼120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope to perform near-infrared high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4 yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions.

Description: Context. This is the first publication from the DEATHSTAR project. The overall goal of the project is to reduce the uncertainties of the observational estimates of mass-loss rates from evolved stars on the Asymptotic Giant Branch (AGB). Aim. The aim in this first publication is to constrain the sizes of the 12CO emitting region from the circumstellar envelopes around 42 mostly southern AGB stars, of which 21 are M-type and 21 are C-type, using the Atacama Compact Array (ACA) at the Atacama Large Millimeter/submillimeter Array. The symmetry of the outflows is also investigated. Methods. Line emission from 12CO J = 2→1 and 3→2 from all of the sources were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. A detailed radiative transfer analysis will be presented in a future publication. The major and minor axis of the best-fit Gaussian at the line center velocity of the 12CO J = 2→1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the Gaussian major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. For a subsample of sources, the measured emission distribution is compared to predictions from previous best-fit radiative transfer modeling results. Results. We find that the CO envelope sizes are, in general, larger for C-type than for M-type AGB stars, which is as expected if the CO/H2 ratio is larger in C-type stars. Furthermore, the measurements show a relation between the measured (Gaussian) 12CO J = 2→1 size and circumstellar density that, while in broad agreement with photodissociation calculations, reveals large scatter and some systematic differences between the different stellar types. For lower mass-loss-rate irregular and semi-regular variables of both M- and C-type AGB stars, the 12CO J = 2→1 size appears to be independent of the ratio of the mass-loss rate to outflow velocity, which is a measure of circumstellar density. For the higher mass-loss-rate Mira stars, the 12CO J = 2→1 size clearly increases with circumstellar density, with larger sizes for the higher CO-abundance C-type stars. The M-type stars appear to be consistently smaller than predicted from photodissociation theory. The majority of the sources have CO envelope sizes that are consistent with a spherically symmetric, smooth outflow, at least on larger scales. For about a third of the sources, indications of strong asymmetries are detected. This is consistent with what was found in previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. Conclusions. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to deduce the nature and origin of the complex circumstellar dynamics revealed by our ACA observations.

Description: Context. The asymptotic giant branch (AGB) marks the final evolutionary stage of stars with initial masses between ~0.8 and 8 M⊙. During this phase, stars undergo copious mass loss, which contributes significantly to the enrichment of the interstellar medium. The well-accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the density-enhanced and extended AGB stellar atmospheres. The details of the mass-loss process are not yet well understood, however. For oxygen-rich AGB stars, which are the focus of this study, the dust grains that drive the wind are expected to scatter visible light very efficiently because their sizes are relative large. Aims. We study the distribution of dust in the inner wind of oxygen-rich AGB stars to advance our understanding of the wind-driving process. Methods. We observed light scattered off dust grains that form around three oxygen-rich AGB stars (W Hya, SW Vir, and R Crt) with mass-loss rates between 10−7 and 10−6 M⊙ yr−1 using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL with three filters centred at 0.65, 0.75 and 0.82 μm. We compared the observed morphologies and the spectral dependence of the scattered light between the three sources and determined the radial profile, per image octant, of the dust density distribution around the closest target, W Hya. Results. We find the distribution of dust to be asymmetric for the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be inferred from the ZIMPOL data, but we speculate that it might be the consequence of a circumstellar disc or of the action of strong magnetic fields. The dust grains polarise light more efficiently at 0.65 μm for R Crt and SW Vir and at 0.82 μm for W Hya. This indicates that at the time of the observations, the grains around SW Vir and R Crt had sizes

Description: Context. The CHaracterising ExOPlanet Satellite (CHEOPS) is a mission dedicated to the search for exoplanetary transits through high precision photometry of bright stars already known to host planets. The telescope will provide the unique capability of determining accurate radii for planets whose masses have already been measured from ground-based spectroscopic surveys. This will allow a first-order characterisation of the planets’ internal structure through the determination of the bulk density, providing direct insight into their composition. By identifying transiting exoplanets with high potential for in-depth characterisation, CHEOPS will also provide prime targets for future instruments suited to the spectroscopic characterisation of exoplanetary atmospheres. Aims. The CHEOPS simulator has been developed to perform detailed simulations of the data which is to be received from the CHEOPS satellite. It generates accurately simulated images that can be used to explore design options and to test the on-ground data processing, in particular, the pipeline producing the photometric time series. It is, thus, a critical tool for estimating the photometric performance expected in flight and to guide photometric analysis. It can be used to prepare observations, consolidate the noise budget, and asses the performance of CHEOPS in realistic astrophysical fields that are difficult to reproduce in the laboratory. Methods. The simulator has been implemented as a highly configurable tool called CHEOPSim, with a web-based user interface. Images generated by CHEOPSim take account of many detailed effects, including variations of the incident signal flux and backgrounds, and detailed modelling of the satellite orbit, pointing jitter and telescope optics, as well as the CCD response, noise and readout. Results. The simulator results presented in this paper have been used in the context of validating the data reduction processing chain, in which image time series generated by CHEOPSim were used to generate light curves for simulated planetary transits across real and simulated targets. Independent analysts were successfully able to detect the planets and measure their radii to an accuracy within the science requirements of the mission: for an Earth-sized planet with an orbital period of 50 days orbiting a Sun-like target with magnitude V = 6, the median measured value of the planet to star radius ratio, Rp/Rs, was 0.00923 ± 0.00054(stat) ± 0.00019(syst), compared to a true input value of 0.00916. For a Neptune-sized planet with an orbital period of 13 days orbiting a target with spectral type K5V and magnitude V = 12, the median measured value of Rp/Rs was 0.05038 ± 0.00061(stat) ± 0.00031(syst), compared to a true input value of 0.05.

Description: Context. ALMA high angular resolution observations of the dust and CO emission have already revealed signatures of protoplanets embedded in protoplanetary disks. These detections are around single T Tauri stars, while exoplanet surveys reveal that planets can also form in binary (or multiple) systems, either in circumstellar or circumbinary orbits. Aims. We searched for indirect evidence for planet formation in the multiple system GG Tau A, which harbors the most massive circumbinary disk among T Tauri stars. Methods. We performed CO(2–1) ALMA Cycle 6 observations of GG Tau A at 0.3″ resolution. The images confirm the “hot spot” detected at higher frequencies, but also reveal prominent spiral-like features. We modeled these features using the analytic prescription for the linear perturbation regime induced by low-mass planets. Results. The brightest spiral is well reproduced by a density wave excited by a protoplanet (GG Tau Ac) at the hot-spot location (290 au), just outside the dust ring. The absence of a clear gap (in gas or dust) at the planet location implies that its mass is significantly lower than that of Jupiter, i.e., of about the mass of Neptune or lower. Furthermore, other prominent (trailing) spiral patterns can be represented by adding one (or more) planet(s) at larger orbital radii, with the most obvious candidate located near the 2:1 mean-motion resonance with GG Tau Ac. Conclusions. The (proto-)planet GG Tau Ac appears to externally confine the ring in a stable configuration, explaining its high mass. Our results also suggest that planets similar in mass to Neptune may form in dense circumbinary disks orbiting (wide) binary stars. In the GG Tau case, orbital resonances appear to play an important role in shaping this multiple circumbinary planet system.