Search

Context. Recent observations found close-in planets with significant atmospheres of hydrogen and helium in great abundance. These are the so-called super-Earths and mini-Neptunes. Their atmospheric composition suggests that they formed early during the gas-rich phase of the circumstellar disk and were able to avoid becoming hot Jupiters. As a possible explanation, recent studies explored the recycling hypothesis and showed that atmosphere-disk recycling is able to fully compensate for radiative cooling and thereby halt Kelvin-Helmholtz contraction to prevent runaway gas accretion. Aims. To understand the parameters that determine the efficiency of atmospheric recycling, we extend our earlier studies by exploring the effects of the core mass, the effect of circumstellar gas on sub-Keplerian orbits (headwind), and the optical depth of the surrounding gas on the recycling timescale. Additionally, we analyze their effects on the size and mass of the forming atmosphere. Methods. We used three-dimensional (3D) radiation-hydrodynamic simulations to model a local shearing box centered on the planet. Our planet is located at a separation of ap= 0.1 au from its solar-type host star, and we scanned the core mass range from 1 to 10 MEarth. In order to measure and track the recycling of the atmosphere, we employed tracer particles as well as tracer fluids after thermodynamic equilibrium was reached. Results. For the explored parameter space, all simulations eventually reach an equilibrium where heating due to hydrodynamic recycling fully compensates radiative cooling. In this equilibrium, the atmosphere-to-core mass ratio stays well below 10%, preventing the atmosphere from becoming self-gravitating and entering runaway gas accretion. Higher core masses cause the atmosphere to become turbulent, which further enhances recycling. Compared to the core mass, the effect of the headwind on the recycling timescale is negligible. The opacity has no significant effect on the recycling timescale, which demonstrates that the Kelvin-Helmholtz contraction timescale and the atmosphere-disk recycling timescale are independent of each other. Conclusions. Even for our highest core mass of 10 MEarth, atmosphere-disk recycling is efficient enough to fully compensate for radiative cooling and prevent the atmosphere from becoming self-gravitating. Hence, in-situ formation of hot Jupiters is very unlikely, and migration of gas giants is a key process required to explain their existence. Our findings imply that atmosphere-disk recycling is the most natural explanation for the prevalence of close-in super-Earths and mini-Neptunes. [ABSTRACT FROM AUTHOR]

Context. Massive stars shape their surrounding medium through the force of their stellar winds, which collide with the circumstellar medium. Because the characteristics of these stellar winds vary over the course of the evolution of the star, the circumstellar matter becomes a reflection of the stellar evolution and can be used to determine the characteristics of the progenitor star. In particular, whenever a fast wind phase follows a slow wind phase, the fast wind sweeps up its predecessor in a shell, which is observed as a circumstellar nebula. Aims. We make 2D and 3D numerical simulations of fast stellar winds sweeping up their slow predecessors to investigate whether numerical models of these shells have to be 3D, or whether 2D models are sufficient to reproduce the shells correctly. Methods. We use the MPI-AMRVAC code, using hydrodynamics with optically thin radiative losses included, to make numerical models of circumstellar shells around massive stars in 2D and 3D and compare the results. We focus on those situations where a fast Wolf-Rayet star wind sweeps up the slower wind emitted by its predecessor, being either a red supergiant or a luminous blue variable. Results. As the fast Wolf-Rayet wind expands, it creates a dense shell of swept up material that expands outward, driven by the high pressure of the shocked Wolf-Rayet wind. These shells are subject to a fair variety of hydrodynamic-radiative instabilities. If the Wolf-Rayet wind is expanding into the wind of a luminous blue variable phase, the instabilities will tend to form a fairly small-scale, regular filamentary lattice with thin filaments connecting knotty features. If the Wolf-Rayet wind is sweeping up a red supergiant wind, the instabilities will form larger interconnected structures with less regularity. The numerical resolution must be high enough to resolve the compressed, swept-up shell and the evolving instabilities, which otherwise may not even form. Conclusions. Our results show that 3D models, when translated to observed morphologies, give realistic results that can be compared directly to observations. The 3D structure of the nebula will help to distinguish different progenitor scenarios. [ABSTRACT FROM AUTHOR]

The interstellar medium (ISM) is subject, on one hand, to heating and cooling processes that tend to segregate it into distinct phases due to thermal instability (TI), and on the other, to turbulence-driving mechanisms that tend to produce strong nonlinear fluctuations in all the thermodynamic variables. In this regime, large-scale turbulent compressions in the stable warm neutral medium (WNM) dominate the clump-formation process rather than the linear development of TI. Cold clumps formed by this mechanism are still often bounded by sharp density and temperature discontinuities, which however are not contact discontinuities as in the classical 2-phase model, but rather "phase transition fronts", across which there is net mass and momentum flux from the WNM into the clumps. The clumps grow mainly by accretion through their boundaries, are in both thermal and ram pressure balance with their surroundings, and are internally turbulent as well, thus also having significant density fluctuations inside. The temperature and density of the cold and warm gas around the phase transition fronts fluctuate with time and location due to fluctuations in the turbulent pressure. Moreover, shock-compressed diffuse unstable gas can remain in the unstable regime for up to a few Myr before it undergoes a phase transition to the cold phase, and is furthermore constantly replenished by the turbulence. These processes populate the classically forbidden density and temperature ranges. Since gas at all temperatures appears to be present in bi- or tri-stable turbulence, we conclude that the word "phase" applies only locally, surrounding phase transition sites in the gas. Globally, the word "phase" must relax its meaning to simply denote a certain temperature or density range. [ABSTRACT FROM AUTHOR]

Partially ionized plasmas are found in many different astrophysical environments. The study of partially ionized plasmas is of great interest for solar physics because some layers of the solar atmosphere (photosphere and chromosphere) as well as solar structures, such as spicules and prominences, are made of these kinds of plasmas. To our knowledge, despite it being known that the adiabatic coefficient, γ, or the first adiabatic exponent, Γ1, depend on the ionization degree, this fact has been disregarded in all the studies related to magnetohydrodynamic waves in solar partially ionized plasmas. However, in other astrophysical areas, the dependence of γ or Γ1 on the plasma ionization degree has been taken into account. Therefore, our aim here is to study how, in a plasma with prominence physical properties, the joint action of the temperature, density, and ionization degree modifies the numerical values of the first adiabatic exponent Γ1 which affects the adiabatic sound speed and the period of slow waves. In our computations, we have used two different approaches; first of all, we assume local thermodynamic equilibrium (LTE) and, later, we consider a non-local thermodynamic equilibrium (non-LTE) model. When comparing the results in the LTE and non-LTE cases, the numerical values of Γ1 are clearly different for both and they are probably strongly dependent on the assumed model which determines how the ionization degree evolves with temperature. Finally, the effect of the ionization degree dependence of Γ1 on the period of slow waves has been determined showing that it can be of great importance for seismological studies of partially ionized solar structures. [ABSTRACT FROM AUTHOR]

Magnetohydrodynamics (MHD)-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC). [ABSTRACT FROM AUTHOR]

Context. Massive stars form in magnetized and turbulent environments and are often located in stellar clusters. The accretion and outflows mechanisms associated with forming massive stars and the origin of the stellar multiplicity of their system are poorly understood. Aims. We study the effect of magnetic fields and turbulence on the accretion mechanism of massive protostars and their multiplicity. We also focus on disk formation as a prerequisite for outflow launching. Methods. We present a series of four radiation-magnetohydrodynamical simulations of the collapse of a massive magnetized, turbulent core of 100 M⊙ with the adaptive-mesh-refinement code RAMSES, including a hybrid radiative transfer method for stellar irradiation and ambipolar diffusion. We varied the Mach and Alfvénic Mach numbers to probe sub- and super-Alfvénic turbulence and sub- and supersonic turbulence regimes. Results. Sub-Alfvénic turbulence leads to single stellar systems, and super-Alfvénic turbulence leads to binary formation from disk fragmentation following the collision of spiral arms, with mass ratios of 1.1–1.6 and a separation of several hundred AU that increases with initial turbulent support and with time. In these runs, infalling gas reaches the individual disks through a transient circumbinary structure. Magnetically regulated, thermally dominated (plasma beta β > 1) Keplerian disks form in all runs, with sizes 100–200 AU and masses 1–8 M⊙. The disks around primary and secondary sink particles have similar properties. We obtain mass accretion rates of ~10−4M⊙ yr−1 onto the protostars and observe higher accretion rates onto the secondary stars than onto their primary star companion. The primary disk orientation is found to be set by the initial angular momentum carried by turbulence rather than by magnetic fields. Even without turbulence, axisymmetry and north–south symmetry with respect to the disk plane are broken by the interchange instability and thermally dominated streamers, respectively. Conclusions. Small (≲300 AU) massive protostellar disks such as those that are frequently observed today can so far only be reproduced in the presence of (moderate) magnetic fields with ambipolar diffusion, even in a turbulent medium. The interplay between magnetic fields and turbulence sets the multiplicity of stellar clusters. A plasma beta β > 1 is a good indicator for distinguishing streamers and individual disks from their surroundings. [ABSTRACT FROM AUTHOR]

We report the discovery of the massive hot Jupiter NGTS-13b by the Next Generation Transit Survey (NGTS). The V = 12.7 host star is likely in the subgiant evolutionary phase with logg* = 4.04 ± 0.05, Teff = 5819 ± 73 K, M* = 1.30 −0.18+0.11 $^{+0.11}_{-0.18}$ − 0.18 + 0.11 M⊙, and R* = 1.79 ± 0.06 R⊙. The NGTS detected a transiting planet with a period of P = 4.12 days around the star, which was later validated with the Transiting Exoplanet Survey Satellite (TESS; TIC 454069765). We confirm the planet using radial velocities from the CORALIE spectrograph. Using NGTS and TESS full-frame image photometry combined with CORALIE radial velocities, we determine NGTS-13b to have a radius of RP = 1.142 ± 0.046 RJup, a mass of MP = 4.84 ± 0.44 MJup, and an eccentricity of e = 0.086 ± 0.034. Previous studies have suggested that ~4 MJup may be the border separating two formation scenarios (e.g., core accretion and disk instability) and that massive giant planets share similar formation mechanisms as lower-mass brown dwarfs. NGTS-13b is just above 4 MJup, making it an important addition to the statistical sample needed to understand the differences between various classes of substellar companions. The high metallicity of NGTS-13, [Fe/H] = 0.25 ± 0.17, does not support previous suggestions that massive giants are found preferentially around lower metallicity host stars, but NGTS-13b does support findings that more massive and evolved hosts may have a higher occurrence of close-in massive planets than lower-mass unevolved stars. [ABSTRACT FROM AUTHOR]

The interplanetary and magnetospheric phenomena time-coincident with intense geomagnetically induced current (GIC) > 10 A and > 30 A events during 21 years (1999 through 2019) at the Mantsala, Finland (57.9° magnetic latitude) gas pipeline have been studied. Although forward shocks and substorms are predominant causes of intense GICs, some newly discovered geoeffective interplanetary features are: solar wind plasma parcel (PP) impingements, possible interplanetary magnetic field (IMF) northward (Bn) and southward (Bs) turnings, and reverse shocks. The PPs are possibly the loop and filament portions of coronal mass ejections (CMEs). From a study of > 30 A GIC events, it is found that supersubstorm (SSS: SML < --2500 nT) and intense substorm (--2500 nT < SML < --2000 nT) auroral electrojet intensifications are the most frequent (76%) cause of all of these GIC events. These events occur most often (76%) in superstorm (SYM-H < --250 nT) main phases, but they can occur in other storm phases and lesser intensity storms as well. After substorms, PPs were the most frequent causes of Mantsala GIC > 30 A events. Forward shocks were the third most frequent cause of the > 30 A events. Shock-related GICs were observed to occur at all local times. The two "Halloween" superstorms of 29-30 and 30-31 October 2003 produced by far the greatest number of GICs in the interval of study (9 > 30 A GICs and 168 > 10 A GICs). In the first Halloween superstorm, a shock-triggered SSS (SML < --3548 nT) caused 33, 57, 51 and 52 A GICs. The 57 A GIC was the most intense event of the superstorm and of this study. It is possible that this SSS is a new form of substorm. Equally intense magnetic storms were also studied but their related GICs were far less numerous and less intense. [ABSTRACT FROM AUTHOR]

The development of numerical models and tools which have operational space weather potential is an increasingly important area of research. This study presents recent Canadian efforts toward the development of a numerical framework for Sun-to-Earth simulations of solar wind disturbances. This modular three-dimensional (3D) simulation framework is based on a semi-empirical data-driven approach to describe the solar corona and an MHD-based description of the heliosphere. In the present configuration, the semi-empirical component uses the potential field source surface (PFSS) and Schatten current sheet (SCS) models to derive the coronal magnetic field based on observed magnetogram data. Using empirical relations, solar wind properties are associated with this coronal magnetic field. Together with a coronal mass ejection (CME) model, this provides inner boundary conditions for a global MHD model which is used to describe interplanetary propagation of the solar wind and CMEs. The proposed MHD numerical approach makes use of advanced numerical techniques. The 3D MHD code employs a finite-volume discretization procedure with limited piecewise linear reconstruction to solve the governing partial-differential equations. The equations are solved on a body-fitted hexahedral multi-block cubed-sphere mesh and an efficient iterative Newton method is used for time-invariant simulations and an explicit time-marching scheme is applied for unsteady cases. Additionally, an efficient anisotropic block-based refinement technique provides significant reductions in the size of the computational mesh by locally refining the grid in selected directions as dictated by the flow physics. The capabilities of the framework for accurately capturing solar wind structures and forecasting solar wind properties at Earth are demonstrated. Furthermore, a comparison with previously reported results and future space weather forecasting challenges are discussed. [ABSTRACT FROM AUTHOR]

The reconstruction of plasma parameters in the interplanetary medium is very important to understand the interplanetary propagation of solar eruptions and for Space Weather application purposes. Because only a few spacecraft are measuring in situ these parameters, reconstructions are currently performed by running complex numerical Magneto-hydrodynamic (MHD) simulations starting from remote sensing observations of the Sun. Current models apply full 3D MHD simulations of the corona or extrapolations of photospheric magnetic fields combined with semi-empirical relationships to derive the plasma parameters on a sphere centered on the Sun (inner boundary). The plasma is then propagated in the interplanetary medium up to the Earth's orbit and beyond. Nevertheless, this approach requires significant theoretical and computational efforts, and the results are only in partial agreement with the in situ observations. In this paper we describe a new approach to this problem called RIMAP -- Reverse In situ data and MHD APproach. The plasma parameters in the inner boundary at 0.1 AU are derived directly from the in situ measurements acquired at 1 AU, by applying a back reconstruction technique to remap them into the inner heliosphere. This remapping is done by using the Weber and Davies solar wind theoretical model to reconstruct the wind flowlines. The plasma is then re-propagated outward from 0.1 AU by running a MHD numerical simulation based on the PLUTO code. The interplanetary spiral reconstructions obtained with RIMAP are not only in a much better agreement with the in situ observations, but are also including many more small-scale longitudinal features in the plasma parameters that are not reproduced with the approaches developed so far. [ABSTRACT FROM AUTHOR]

Context. Dense molecular filaments are ubiquituous in the interstellar medium, yet their internal physical conditions and the role of gravity, turbulence, the magnetic field, radiation, and the ambient cloud during their evolution remain debated. Aims. We study the kinematics and physical conditions in the Musca filament, the ambient cloud, and the Chamaeleon-Musca complex to constrain the physics of filament formation. Methods. We produced CO(2–1) isotopologue maps with the APEX telescope that cut through the Musca filament. We further study a NANTEN2 12CO(1–0) map of the full Musca cloud, H I emission of the Chamaeleon-Musca complex, a Planck polarisation map, line radiative transfer models, Gaia data, and synthetic observations from filament formation simulations. Results. The Musca cloud, with a size of ~3–6 pc, contains multiple velocity components. Radiative transfer modelling of the CO emission indicates that the Musca filament consists of a cold (~10 K), dense (nH2 ∼ 104 cm−3) crest, which is best described with a cylindrical geometry. Connected to the crest, a separate gas component at T ~ 15 K and nH2 ∼ 103 cm−3 is found, the so-called strands. The velocity-coherent filament crest has an organised transverse velocity gradient that is linked to the kinematics of the nearby ambient cloud. This velocity gradient has an angle ≥30° with respect to the local magnetic field orientation derived from Planck, and the magnitude of the velocity gradient is similar to the transonic linewidth of the filament crest. Studying the large scale kinematics, we find coherence of the asymmetric kinematics from the 50 pc H I cloud down to the Musca filament. We also report a strong [C18O]/[13CO] abundance drop by an order of magnitude from the filament crest to the strands over a distance <0.2 pc in a weak ambient far-ultraviolet (FUV) field. Conclusions. The dense Musca filament crest is a long-lived (several crossing times), dynamic structure that can form stars in the near future because of continuous mass accretion replenishing the filament. This mass accretion on the filament appears to be triggered by a H I cloud–cloud collision, which bends the magnetic field around dense filaments. This bending of the magnetic field is then responsible for the observed asymmetric accretion scenario of the Musca filament, which is, for instance, seen as a V-shape in the position–velocity (PV) diagram. [ABSTRACT FROM AUTHOR]

Context. The density of the ambient medium around a supermassive black hole (SMBH) and the way it varies with distance plays an important role in our understanding of the inflow-outflow mechanisms in the Galactic centre (GC). This dependence is often fitted by spherical power-law profiles based on observations in the X-ray, infrared (IR), submillimetre (submm), and radio domains. Aims. Nevertheless, the density profile is poorly constrained at the intermediate scales of 1000 Schwarzschild radii (Rs). Here we independently constrain the spherical density profile using the stellar bow shock of the star S2 which orbits the SMBH at the GC with the pericentre distance of 14.4 mas (∼1500 Rs). Methods. Assuming an elliptical orbit, we apply celestial mechanics and the theory of bow shocks that are at ram pressure equilibrium. We analyse the measured IR flux density and magnitudes of S2 in the L′-band (3.8 micron) obtained over seven epochs in the years between 2004–2018. We put an upper limit on the emission from S2's associated putative bow shock and constrain the density profile of the ambient medium. Results. We detect no significant change in S2 flux density until the recent periapse in May 2018. The intrinsic flux variability of S2 is at the level of 2–3%. Based on the dust-extinction model, the upper limit on the number density at the S2 periapse is ∼1.87  ×  109 cm−3, which yields a density slope of at most 3.20. Using the synchrotron bow-shock emission, we obtain the ambient density of ≲1.01  ×  105 cm−3 and a slope of ≲1.47. These values are consistent with a wide variety of media from hot accretion flows to potentially colder and denser media comparable in properties to broad-line-region clouds. However, a standard thin disc can be excluded at the distance of S2's pericentre. Conclusions. With the current photometry sensitivity of 0.01 mag, we are not able to make stringent constraints on the density of the ambient medium in the GC using S2-star observations. We can distinguish between hot accretion flows and thin, cold discs, where the latter can be excluded at the scale of the S2 periapse. Future observations of stars in the S cluster using instruments such as Mid-IR Extremely Large Telescope Imager and Spectrograph at Extremely Large Telescope with the photometric sensitivity of as much as 10−3 mag will allow the GC medium to be probed at intermediate scales at densities as low as ∼700 cm−3 in case of non-thermal bow-shock emission. The new instrumentation, in combination with discoveries of stars with smaller pericentre distances, will help to independently constrain the density profile around Sagittarius A* (Sgr A*). [ABSTRACT FROM AUTHOR]

Context. The mechanism(s) behind coronal heating still elude(s) direct observation and modelling of viable theoretical processes and the subsequent effect on coronal structures is one of the key tools available to assess possible heating mechanisms. Wave heating via the phase mixing of magnetohydrodynamic (MHD) transverse waves has been proposed as a possible way to convert magnetic energy into thermal energy, but MHD models increasingly suggest this is not an efficient enough mechanism. Aims. We modelled heating by phase mixing transverse MHD waves in various configurations in order to investigate whether certain circumstances can enhance the heating sufficiently to sustain the million degree solar corona and to assess the impact of the propagation and phase mixing of transverse MHD waves on the structure of the boundary shell of coronal loops. Methods. We used 3D MHD simulations of a pre-existing density enhancement in a magnetised medium and a boundary driver to trigger the propagation of transverse waves with the same power spectrum as measured by the Coronal Multi-Channel Polarimeter. We consider different density structures, boundary conditions at the non-drive footpoint, characteristics of the driver, and different forms of magnetic resistivity. Results. We find that different initial density structures significantly affect the evolution of the boundary shell and that some driver configurations can enhance the heating generated from the dissipation of the MHD waves. In particular, drivers coherent on a larger spatial scale and higher dissipation coefficients can generate significant heating, although it is still insufficient to balance the radiative losses in this setup. Conclusions. We conclude that while phase mixing of transverse MHD waves is unlikely to sustain the thermal structure of the corona, there are configurations that allow for an enhanced efficiency of this mechanism. We provide possible signatures to identify the presence of such configurations, such as the location of where the heating is deposited along the coronal loop. [ABSTRACT FROM AUTHOR]

Context. Uranus is the only planet in the Solar System whose rotation axis and orbital plane are nearly parallel to each other. Uranus is also the planet with the largest angle between the rotation axis and the direction of its magnetic dipole (roughly 59°). Consequently, the shape and structure of its magnetospheric tail is very different to those of all other planets in whichever season one may consider. The only in situ measurements were obtained in January 1986 during a flyby of the Voyager II spacecraft. At that date, Uranus was near solstice time, but unfortunately the data collected by the spacecraft were much too sparse to allow for a clear view of the structure and dynamics of its extended magnetospheric tail. Later numerical simulations revealed that the magnetic tail of Uranus at solstice time is helically shaped with a characteristic pitch of the order of 1000 planetary radii. Aims. We aim to propose a magnetohydrodynamic model for the magnetic tail of Uranus at solstice time. Methods. We constructed our model based on a symmetrised version of the Uranian system by assuming an exact alignment of the solar wind and the planetary rotation axis and an angle of 90° between the planetary magnetic dipole and the rotation axis. We do also postulate that the impinging solar wind is steady and unmagnetised, which implies that the magnetosphere is quasi-steady in the rotating planetary frame and that there is no magnetic reconnection at the magnetopause. Results. One of the main conclusions is that all magnetic field lines forming the extended magnetic tail follow the same qualitative evolution from the time of their emergence through the planet's surface and the time of their late evolution after having been stretched and twisted several times downstream of the planet. In the planetary frame, these field lines move on magnetic surfaces that wind up to form a tornado-shaped vortex with two foot points on the planet (one in each magnetic hemisphere). The centre of the vortex (the eye of the tornado) is a simple double helix with a helical pitch (along the symmetry axis z) λ = τ[vz+Bz/(μ0ρ)1/2] $\lambda\,{=}\,\tau[v_z+B_z/(\mu_0\rho)^{1/2}]$ λ   =   τ [ v z + B z / (μ 0 ρ) 1 / 2 ] , where τ is the rotation period of the planet, μ0 the permeability of vacuum, ρ the mass density, vz the fluid velocity, and Bz the magnetic field where all quantities have to be evaluated locally at the centre of the vortex. In summary, in the planetary frame, the motion of a typical magnetic field of the extended Uranian magnetic tail is a vortical motion, which asymptotically converges towards the single double helix, regardless of the line's emergence point on the planetary surface. [ABSTRACT FROM AUTHOR]

Context. The amount of data collected by spectrometers from radio to ultraviolet (UV) wavelengths opens a new era where the statistical and chemical information contained in the observations can be used concomitantly to investigate the thermodynamical state and the evolution of the interstellar medium (ISM). Aims. In this paper, we study the statistical properties of the HI-to-H2 transition observed in absorption in the local diffuse and multiphase ISM. Our goal is to identify the physical processes that control the probability of occurrence of any line of sight and the origins of the variations of the integrated molecular fraction from one line of sight to another. Methods. The turbulent diffuse ISM is modeled using the RAMSES code, which includes detailed treatments of the magnetohydrodynamics, the thermal evolution of the gas, and the chemistry of H2. The impacts of the UV radiation field, the mean density, the turbulent forcing, the integral scale, the magnetic field, and the gravity on the molecular content of the gas are explored through a parametric study that covers a wide range of physical conditions. The statistics of the HI-to-H2 transition are interpreted through analytical prescriptions and compared with the observations using a modified and robust version of the Kolmogorov-Smirnov test. Results. The analysis of the observed background sources shows that the lengths of the lines of sight follow a flat distribution in logarithmic scale from ~100 pc to ~3 kpc. Without taking into account any variation of the parameters along a line of sight or from one line of sight to another, the results of one simulation, convolved with the distribution of distances of the observational sample, are able to simultaneously explain the position, the width, the dispersion, and most of the statistical properties of the HI-to-H2 transition observed in the local ISM. The tightest agreement is obtained for a neutral diffuse gas modeled over ~200 pc, with a mean density n̅H̅ = 1−2 $\overline{n_{\textrm{H}}} =1{-}2$ n H ¯ = 1 − 2 cm−3, illuminated by the standard interstellar UV radiation field, and stirred up by a large-scale compressive turbulent forcing. Within this configuration, the 2D probability histogram of the column densities of H and H2, poetically called the kingfisher diagram, is remarkably stable and is almost unaltered by gravity, the strength of the turbulent forcing, the resolution of the simulation, or the strength of the magnetic field Bx, as long as Bx < 4 μG. The weak effect of the resolution and our analytical prescription suggest that the column densities of HI are likely built up in large-scale warm neutral medium and cold neutral medium (CNM) structures correlated in density over ~20 pc and ~10 pc, respectively, while those of H2 are built up in CNM structures between ~3 and ~10 pc. Conclusions. Combining the chemical and statistical information contained in the observations of HI and H2 sheds new light on the study of the diffuse matter. Applying this new tool to several atomic and molecular species is a promising perspective to understanding the effects of turbulence, magnetic field, thermal instability, and gravity on the formation and evolution of molecular clouds. [ABSTRACT FROM AUTHOR]

Context. Classical T Tauri stars (CTTs) magnetically interact with their surrounding disks, a process that is thought to regulate their rotational evolution. Aims. We compute torques acting on the stellar surface of CTTs that arise from different accreting (accretion funnels) and ejecting (stellar winds and magnetospheric ejections) flow components. Furthermore, we compare the magnetic braking due to stellar winds in two different systems: isolated (i.e., weak-line T Tauri and main-sequence) and accreting (i.e., classical T Tauri) stars. Methods. We use 2.5D magnetohydrodynamic, time-dependent, axisymmetric simulations that were computed with the PLUTO code. For both systems, the stellar wind is thermally driven. In the star-disk-interaction (SDI) simulations, the accretion disk is Keplerian, viscous, and resistive, and is modeled with an alpha prescription. Two series of simulations are presented, one for each system (i.e., isolated and accreting stars). Results. In classical T Tauri systems, the presence of magnetospheric ejections confines the stellar-wind expansion, resulting in an hourglass-shaped geometry of the outflow, and the formation of the accretion columns modifies the amount of open magnetic flux exploited by the stellar wind. These effects have a strong impact on the stellar-wind properties, and we show that the stellar-wind braking is more efficient in the SDI systems than in the isolated ones. We further derive torque scalings over a wide range of magnetic field strengths for each flow component in an SDI system (i.e., magnetospheric accretion and ejections, and stellar winds), which directly applies a torque on the stellar surface. Conclusions. In all the performed SDI simulations, the stellar wind extracts less than 2% of the mass accretion rate and the disk is truncated by up to 66% of the corotation radius. All simulations show a net spin-up torque. We conclude that in order to achieve a stellar-spin equilibrium, we need either more massive stellar winds or disks that are truncated closer to the corotation radius, which increases the torque efficiency of the magnetospheric ejections. [ABSTRACT FROM AUTHOR]

Context. The chemical composition of protoplanetary disks is expected to impact the composition of the forming planets. Characterizing the diversity of chemical composition in disks and the physicochemical factors that lead to this diversity is consequently of high interest. Aims. The aim of this study is to investigate the chemical evolution from the prestellar phase to the formation of the disk, and to determine the impact that the chemical composition of the cold and dense core has on the final composition of the disk. Methods. We performed 3D nonideal magneto-hydrodynamic (MHD) simulations of a dense core collapse using the adaptive-mesh-refinement RAMSES code. For each particle ending in the young rotationally supported disk, we ran chemical simulations with the three-phase gas-grain chemistry code Nautilus. Two different sets of initial abundances, which are characteristic of cold cores, were considered. The final distributions of the abundances of common species were compared to each other, as well as with the initial abundances of the cold core. Results. We find that the spatial distributions of molecules reflect their sensitivity to the temperature distribution. The main carriers of the chemical elements in the disk are usually the same as the ones in the cold core, except for the S-bearing species, where HS is replaced by H2S3, and the P-bearing species, where atomic P leads to the formation of PO, PN, HCP, and CP. However, the abundances of less abundant species change over time. This is especially the case for "large" complex organic molecules (COMs) such as CH3CHO, CH3NH2, CH3OCH3, and HCOOCH3 which see their abundances significantly increase during the collapse. These COMs often present similar abundances in the disk despite significantly different abundances in the cold core. In contrast, the abundances of many radicals decrease with time. A significant number of species still show the same abundances in the cold core and the disk, which indicates efficient formation of these molecules in the cold core. This includes H2O, H2CO, HNCO, and "small" COMs such as CH3OH, CH3CN, and NH2CHO. We computed the MHD resistivities within the disk for the full gas–grain chemical evolution and find results in qualitative agreement with the literature assuming simpler chemical networks. Conclusions. In conclusion, the chemical content of prestellar cores is expected to affect the chemical composition of disks. The impact is more or less important depending on the type of species. Users of stand-alone chemical models of disks should pay special attention to the initial abundances they choose. [ABSTRACT FROM AUTHOR]

The blind H I survey Arecibo Galaxy Environment Survey (AGES) detected several unresolved sources in the Virgo cluster, which do not have optical counterparts in the Sloan Digital Sky Survey. The origin of these dark clouds is unknown. They might be crucial objects since they could be the so-called dark galaxies, that is, the dark matter halos without stellar content that are expected from cosmological simulations. In order to reveal the nature of the dark clouds, we took a deep optical image of one them, AGESVC1 282, with the newly-commissioned 1.4 m Milanković Telescope. After observing it for 10.4 h in the L-filter, the image reached a surface-brightness limit of about 29.1 mag arcsec−2 in V. No optical counterpart was detected. We placed an upper limit on the V-band luminosity of the object of 1.1 × 107 L⊙, giving a stellar mass below 1.4 × 107 M⊙ and a H I-to-stellar mass ratio above 3.1. By inspecting archival H I observations of the surrounding region, we found that none of the standard explanations for optically dark H I clouds fits the available constraints on this object. [ABSTRACT FROM AUTHOR]

Context. AB Aur is a Herbig Ae star that hosts a prototypical transition disk. The disk shows a plethora of features connected with planet formation mechanisms, such as spiral arms, dust cavities, and dust traps. Understanding the physical and chemical characteristics of these features is crucial to advancing our knowledge of the planet formation processes. Aims. We aim to characterize the gaseous disk around the Herbig Ae star AB Aur. A complete spectroscopic study was performed using NOEMA to determine the physical and chemical conditions with high spatial resolution. Methods. We present new NOrthern Extended Millimeter Array (NOEMA) interferometric observations of the continuum and 12CO, 13CO, C18O, H2CO, and SO lines obtained at high resolution. We used the integrated intensity maps and stacked spectra to derive reliable estimates of the disk temperature. By combining our 13CO and C18O observations, we computed the gas-to-dust ratio along the disk. We also derived column density maps for the different species and used them to compute abundance maps. The results of our observations were compared with a set of Nautilus astrochemical models to obtain insight into the disk properties. Results. We detected continuum emission in a ring that extends from 0.6′′ to ~2.0′′, peaking at 0.97′′ and with a strong azimuthal asymmetry. The molecules observed show different spatial distributions, and the peaks of the distributions are not correlated with the binding energy. Using H2CO and SO lines, we derived a mean disk temperature of 39 K. We derived a gas-to-dust ratio that ranges from 10 to 40 along the disk. Abundance with respect to 13CO for SO (~2 × 10−4) is almost one order of magnitude greater than the value derived for H2CO (1.6 × 10−5). The comparison with Nautilus models favors a disk with a low gas-to-dust ratio (40) and prominent sulfur depletion. Conclusions. From a very complete spectroscopic study of the prototypical disk around AB Aur, we derived, for the first time, the gas temperature and the gas-to-dust ratio along the disk, providing information that is essential to constraining hydrodynamical simulations. Moreover, we explored the gas chemistry and, in particular, the sulfur depletion. The derived sulfur depletion is dependent on the assumed C/O ratio. Our data are better explained with C/O ~ 0.7 and S/H = 8 × 10−8. [ABSTRACT FROM AUTHOR]

Aims. Investigating the process of matter accretion onto forming stars through scaled experiments in the laboratory is important in order to better understand star and planetary system formation and evolution. Such experiments can indeed complement observations by providing access to the processes with spatial and temporal resolution. A previous investigation revealed the existence of a two-component stream: a hot shell surrounding a cooler inner stream. The shell was formed by matter laterally ejected upon impact and refocused by the local magnetic field. That laboratory investigation was limited to normal incidence impacts. However, in young stellar objects, the complex structure of magnetic fields causes variability of the incident angles of the accretion columns. This led us to undertake an investigation, using laboratory plasmas, of the consequence of having a slanted accretion impacting a young star. Methods. Here, we used high power laser interactions and strong magnetic field generation in the laboratory, complemented by numerical simulations, to study the asymmetry induced upon accretion structures when columns of matter impact the surface of young stars with an oblique angle. Results. Compared to the scenario where matter accretes perpendicularly to the star surface, we observe a strongly asymmetric plasma structure, strong lateral ejecta of matter, poor confinement of the accreted material, and reduced heating compared to the normal incidence case. Thus, slanted accretion is a configuration that seems to be capable of inducing perturbations of the chromosphere and hence possibly influencing the level of activity of the corona. [ABSTRACT FROM AUTHOR]

High-mass clusters at low redshifts have been intensively studied at various wavelengths. However, while more distant objects at lower masses constitute the bulk population of future surveys, their physical state remain poorly explored to date. In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC 102, a relatively low-mass system (M500 ∼ 2 × 1014 M⊙) at z = 0.97 detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC 102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC 102, obtaining relatively tight constraints up to about ∼r500, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC 102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z ∼ 1, especially with low signal-to-noise ratio data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data. This is promising for the study of high redshift clusters from the combination of eROSITA and high resolution SZ instruments and will complement the new generation of optical surveys from facilities such as LSST and Euclid. [ABSTRACT FROM AUTHOR]

The interaction between the solar wind and the Martian induced magnetosphere can lead to the formation of various regions with different plasma and magnetic-field characteristics. In this paper, these structures are investigated based on the plasma and magnetic-field measurements from Mars Atmosphere and Volatile EvolutioN (MAVEN). We find that the structures upstream of Mars are similar to those around Earth: both have a bow shock, magnetosheath, magnetopause, and a magnetosphere or induced magnetosphere. The inner part of Martian magnetosheath is called a plasma depletion region (PDR), similar to the plasma depletion layer upstream of the Earth's magnetopause, in which the magnetosheath magnetic fields are piled up and the magnetosheath plasmas (including ions and electrons) are partially depleted. Several cases of PDRs are examined in detail. The hotter plasmas in PDRs are squeezed out along the enhanced magnetic field, resulting in the decrease of the plasma beta, the plasma density, and the ion temperature. The boundary between the magnetosheath and the induced magnetosphere is called the magnetopause, which can be identified as a magnetohydrodynamic discontinuity, either tangential discontinuity (TD) or rotational discontinuity, where the magnetic field changes its orientation. Tangential discontinuities with an insignificant normal component (BN ≈ 0) of the magnetic field are the focus of this study. This discontinuity separates the magnetosheath H+ ions from the heavy ions (e.g. O+, O 2+ $_{2}^{+}$ 2 + ) in the induced magnetosphere. Inside a TD, ions from both sides are mixed. There are 3332 boundary crossings by MAVEN in 2015, 1075 cases of which are identified as the TD (including the potential TD). Tangential discontinuities at Mars are at higher locations in the southern hemisphere and have an average thickness of ~200 km, mostly ranging from 50 to 400 km. The sample of TD is a decreasing function of θ (θ is the magnetic field rotation angle on the two sides of the TD). The PDRs in front of TDs are thicker in the northern hemisphere. From the sub-solar point to the Mars tail, PDR thickness increases and the proton number density and temperature decrease. [ABSTRACT FROM AUTHOR]

Context. Due to the effects that they can have on the atmospheres of exoplanets, stellar winds have recently received significant attention in the literature. Alfvén-wave-driven 3D magnetohydrodynamic models, which are increasingly used to predict stellar wind properties, contain unconstrained parameters and rely on low-resolution stellar magnetograms. Aims. In this paper, we explore the effects of the input Alfvén wave energy flux and the surface magnetogram on the wind properties predicted by the Alfvén Wave Solar Model (AWSoM) model for both the solar and stellar winds. Methods. We lowered the resolution of two solar magnetograms during solar cycle maximum and minimum using spherical harmonic decomposition. The Alfvén wave energy was altered based on non-thermal velocities determined from a far ultraviolet spectrum of the solar twin 18 Sco. Additionally, low-resolution magnetograms of three solar analogues, 18 Sco, HD 76151, and HN Peg, were obtained using Zeeman Doppler imaging and used as a proxy for the solar magnetogram. Finally, the simulated wind properties were compared to Advanced Composition Explorer (ACE) observations. Results. AWSoM simulations using well constrained input parameters taken from solar observations can reproduce the observed solar wind mass loss and angular momentum loss rates. The simulated wind velocity, proton density, and ram pressure differ from ACE observations by a factor of approximately two. The resolution of the magnetogram has a small impact on the wind properties and only during cycle maximum. However, variation in Alfvén wave energy influences the wind properties irrespective of the solar cycle activity level. Furthermore, solar wind simulations carried out using the low-resolution magnetogram of the three stars instead of the solar magnetogram could lead to an order of a magnitude difference in the simulated solar wind properties. Conclusions. The choice in Alfvén energy has a stronger influence on the wind output compared to the magnetogram resolution. The influence could be even stronger for stars whose input boundary conditions are not as well constrained as those of the Sun. Unsurprisingly, replacing the solar magnetogram with a stellar magnetogram could lead to completely inaccurate solar wind properties, and should be avoided in solar and stellar wind simulations. Further observational and theoretical work is needed to fully understand the complexity of solar and stellar winds. [ABSTRACT FROM AUTHOR]

The mass function of clumps observed in molecular clouds raises interesting theoretical issues, especially in its relation to the stellar initial mass function (IMF). We propose a statistical model of the mass function of prestellar cores (CMF), formed in self-gravitating isothermal clouds at a given stage of their evolution. The latter is characterized by the mass-density probability distribution function (-PDF), which is a power-law with slope q. The different molecular clouds are divided into ensembles according to the PDF slope and each ensemble is represented by a single spherical cloud. The cores are considered as elements of self-similar structure typical for fractal clouds and are modeled by spherical objects populating each cloud shell. Our model assumes relations between size, mass, and density of the statistical cores. Out of these, a core mass-density relationship p / mx is derived where x = 1=(1 + q). We find that q determines the existence or nonexistence of a threshold density for core collapse. The derived general CMF is a power law of slope -1 while the CMF of gravitationally unstable cores has a slope (-1 + x=2), comparable with the slopes of the high-mass part of the stellar IMF and of observational CMFs. [ABSTRACT FROM AUTHOR]

To understand the role that active galactic nuclei (AGN) feedback plays in galaxy evolution, we need in-depth studies of the multiphase structure and energetics of galaxy-wide outflows. In this work, we present new, deep (~50 h) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately ±1200 km s-1 and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located ~6 kpc northwest from the galaxy moving away at ~900 km s-1, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, and warm and cold molecular gas phases, is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for ~40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central ~1.5 kpc region of only ~3 Myr, a factor of ~2 shorter than the depletion time by star formation activity. [ABSTRACT FROM AUTHOR]

Aims. We present a fluid model that has been developed for multicomponent two-temperature magnetized plasmas in chemical nonequilibrium for the partially to fully ionized collisional regimes. We focus on transport phenomena with the aim of representing the atmosphere of the Sun. Methods. This study is based on an asymptotic fluid model for multicomponent plasmas derived from kinetic theory, yielding a rigorous description of the dissipative effects. The governing equations and consistent transport properties are obtained using a multiscale Chapman-Enskog perturbative solution to the Boltzmann equation based on a dimensional analysis. The mass disparity between free electrons and heavy particles is accounted for, as well as the influence of the electromagnetic field. We couple this model to the Maxwell equations for the electromagnetic field and derive the generalized Ohm's law for multicomponent plasmas. The model inherits a well-identified mathematical structure leading to an extended range of validity for the Sun's atmospheric conditions. We compute consistent transport properties by means of a spectral Galerkin method using the Laguerre-Sonine polynomial approximation. Two non-vanishing polynomial terms are used when deriving the transport systems for electrons, whereas only one term is retained for heavy particles. Results. In a simplified framework where the plasma is fully ionized, we compare the transport properties for the lower solar atmosphere to conventional expressions for magnetized plasmas attributed to Braginskii, showing a good agreement between both results. For more general partially ionized conditions, representative of the lower solar atmosphere, we compute the muticomponent transport properties corresponding to the species diffusion velocities, heavy-particle and electron heat fluxes, and viscous stress tensor of the model for a helium-hydrogen mixture in local thermodynamic equilibrium. The model is assessed for the 3D radiative magnetohydrodynamic simulation of a pore at the Sun photosphere. The resistive term is found to dominate mainly the dynamics of the electric field at the pore location. The battery term for heavy particles appears to be higher at the pore location and at some intergranulation boundaries. [ABSTRACT FROM AUTHOR]

Context. Cloud-scale surveys of molecular gas reveal the link between giant molecular cloud properties and star formation across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the normal star formation process, while galaxy centers tend to harbor denser gas that exhibits more extreme star formation. At high resolution, however, molecular clouds with exceptional gas properties and star formation activity may also be observed in normal disk environments. In this paper we study the brightest cloud traced in CO(2–1) emission in the disk of nearby spiral galaxy NGC 628. Aims. We characterize the properties of the molecular and ionized gas that is spatially coincident with an extremely bright H II region in the context of the NGC 628 galactic environment. We investigate how feedback and large-scale processes influence the properties of the molecular gas in this region. Methods. High-resolution ALMA observations of CO(2–1) and CO(1−0) emission were used to characterize the mass and dynamical state of the "headlight" molecular cloud. The characteristics of this cloud are compared to the typical properties of molecular clouds in NGC 628. A simple large velocity gradient (LVG) analysis incorporating additional ALMA observations of 13CO(1−0), HCO+(1−0), and HCN(1−0) emission was used to constrain the beam-diluted density and temperature of the molecular gas. We analyzed the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the H II region. Results. The unusually bright headlight cloud is massive (1 − 2 × 107 M⊙), with a beam-diluted density of nH2 = 5 × 104 cm−3 based on LVG modeling. It has a low virial parameter, suggesting that the CO emission associated with this cloud may be overluminous due to heating by the H II region. A young (2 − 4 Myr) stellar population with mass 3 × 105 M⊙ is associated. Conclusions. We argue that the headlight cloud is currently being destroyed by feedback from young massive stars. Due to the large mass of the cloud, this phase of the its evolution is long enough for the impact of feedback on the excitation of the gas to be observed. The high mass of the headlight cloud may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk and receives a sustained inflow of gas that can promote the mass growth of the cloud. [ABSTRACT FROM AUTHOR]

Aims. We study the properties of the coronal magnetic field according to the current sheet source surface (CSSS) model in 1976–2017 for all physically reasonable values of the three model parameters (cusp surface radius Rcs, source surface radius Rss, and current parameter a), and compare the CSSS field with the potential field source surface (PFSS) model field. Methods. We used the synoptic maps of the photospheric magnetic field from the Wilcox Solar Observatory (WSO), National Solar Observatory/Kitt Peak (NSO/KP), and the NSO Synoptic Optical Long-term Investigations of the Sun Vector Spectromagnetograph (SOLIS/VSM) in order to calculate the coronal magnetic field according to the CSSS and PFSS models. We calculated the coronal field strength, its latitudinal variation and neutral line location, as well as its polarity match with the heliospheric magnetic field. Results. The CSSS model can correct the erroneous latitudinal variation of the PFSS model if the source surface is sufficiently far out with respect to the cusp surface (Rss ≥ 3 ⋅ Rcs). The topology of the neutral line only slightly depends on source surface radius or current parameter, but excludes very low values of the cusp surface (Rcs ≤ 1.5). A comparison of the polarities gives an optimum cusp surface radius that varies in time between 2 and 5; a stronger current yields a larger optimum Rcs. Interestingly, the optimum polarity match percentages and optimum radii vary very similarly in the two models over the four solar cycles we studied. Conclusions. The CSSS model can produce a stronger total coronal flux than the PFSS model and correct its latitudinal variation. However, the topology of the CSSS model is rather independent of horizontal currents and remains very similar to that of the PFSS model. Therefore, the CSSS model cannot improve the match of field polarities between corona and heliosphere. [ABSTRACT FROM AUTHOR]

Multi-phase filamentary structures around brightest cluster galaxies (BCG) are likely a key step of AGN-feedback. We observed molecular gas in three cool cluster cores, namely Centaurus, Abell S1101, and RXJ1539.5, and gathered ALMA (Atacama Large Millimeter/submillimeter Array) and MUSE (Multi Unit Spectroscopic Explorer) data for 12 other clusters. Those observations show clumpy, massive, and long (3−25 kpc) molecular filaments, preferentially located around the radio bubbles inflated by the AGN. Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the Hα/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 and 6 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100 and 400 km s−1, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin (and as yet) unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets, or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the radial extent of the filaments, rfil, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short-cooling-time region, where tcool/tff <  20 (9 of 13 sources). The range of tcool/tff of 8−23 at rfil, is likely due to (i) a more complex gravitational potential affecting the free-fall time tff (sloshing, mergers, etc.) and (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time tcool. For some of the sources, rfil lies where the ratio of the cooling time to the eddy-turnover time, tcool/teddy, is approximately unity. [ABSTRACT FROM AUTHOR]

Context. Worldwide very long baseline radio interferometry (VLBI) arrays are expected to obtain horizon-scale images of supermassive black hole candidates and of relativistic jets in several nearby active galactic nuclei. This, together with the expected detection of electromagnetic counterparts of gravitational-wave signals, motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. Aims. The Black Hole Accretion Code (BHAC) is a publicliy available code intended to aid with the modeling of such sources by performing general relativistic magnetohydrodynamical simulations in arbitrary stationary spacetimes. New additions to the code are required in order to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. Methods. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, which are essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field. We also present a general class of prolongation operators for face-allocated variables compatible with them. Results. After presenting several standard tests for the new implementation, we show that the choice of the divergence-control method can produce qualitative differences in the simulation results for scientifically relevant accretion problems. We demonstrate the ability of AMR to decrease the computational costs of black hole accretion simulations while sufficiently resolving turbulence arising from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up for which the new AMR implementation is particularly advantageous. Conclusions. The CT methods and AMR strategies discussed here are currently being used in the simulations performed with BHAC for the generation of theoretical models for the Event Horizon Telescope collaboration. [ABSTRACT FROM AUTHOR]

We present interferometric observations with the Atacama Large Millimeter Array (ALMA) of the free–free continuum and recombination line emission at 1 and 3 mm of the Red Square Nebula surrounding the B[e]-type star MWC 922. The distance to the source, which is unknown, is usually taken to be d = 1.7–3 kpc. The unprecedented angular resolution (up to ~0.′′02) and exquisite sensitivity of these data reveal for the first time the structure and kinematics of the nascent compact ionized region at its center. We imaged the line emission of H30α and H39α, previously detected with single-dish observations, and of H51ɛ, H55γ, and H63δ, detected for the first time in this work. The line emission is seen over a full velocity range of ~180 km s−1 arising in a region of diameter <0.′′14 (less than a few hundred au) in the maser line H30α, which is the most intense transition reported here. We resolve the spatio-kinematic structure of a nearly edge-on disk rotating around a central mass of ~10 M⊙ (d = 1.7 kpc) or ~18 M⊙ (d = 3 kpc), assuming Keplerian rotation. Our data also reveal a fast (~100 km s−1) bipolar ejection (possibly a jet) orthogonal to the disk. In addition, a slow (<15 km s−1) wind may be emanating from the disk. Both, the slow and the fast winds are found to be rotating in a similar manner to the ionized layers of the disk. This represents the first empirical proof of rotation in a bipolar wind expanding at high velocity (~100 km s−1). The launching radius of the fast wind is found to be <30–51 au (i.e., smaller than the inner rim of the ionized disk probed by our observations). We believe that the fast wind is actively being launched, probably by a disk-mediated mechanism in a (possibly accretion) disk around a possible compact companion. We have modeled our observations using the radiative transfer code MORELI. This enables us to describe with unparalleled detail the physical conditions and kinematics in the inner layers of MWC 922, which has revealed itself as an ideal laboratory for studying the interplay of disk rotation and jet-launching. Although the nature of MWC 922 remains unclear, we believe it could be a ~15 M⊙ post-main sequence star in a mass-exchanging binary system. If this is the case, a more realistic value of the distance may be d ~ 3 kpc. [ABSTRACT FROM AUTHOR]