Galactic Disk Research Articles

Study of the Spectral Energy Distribution Model of the Protoplanetary Disk Geometric around the Brown Dwarf CFHT-BD-Tau 4

The accretion circumstellar disk of young stars and the Brown dwarf plays an essential role in the formation and evaluation of the planet. Our main work in this paper is to investigate the geometrical shape model for the protoplanetary disk around one of the Brown Dwarfs. The photometric measurements for the brown dwarf CFHT-BD-Tau 4 were extracted from the Vizier archive. We used a numerical simulation to build a model of the spectral energy distribution of our target CFHT-BD-Tau 4. The spectral energy distribution model was fitted with observational data for the brown dwarf CFHT-BD-Tau 4. A transitional disk has been assumed around CFHT-BD-Tau 4. We obtained physical properties of the two disks and the size of the gap between them by fitting the SED. The gap in the protoplanetary disk proves that a planet formation process occurred around the Brown dwarf.

Generating X-Ray Transit Profiles with Batman

We present an adaptation of the exoplanet transit model code batman, in order to permit the generation of X-ray transits. Our underlying extended coronal model assumes an isothermal plasma that is radially symmetric. While this ignores the effect of bright, active regions, observations of transits in X-rays will require averaging across multiple epochs of data for the foreseeable future, significantly reducing the importance of more complex modeling. Our publicly available code successfully generates the predicted W-shaped transit profile in X-rays due to the optically thin nature of the emission, which concentrates the expected observational emission around the limb of the photospheric stellar disk. We provide some examples based on the best known X-ray transit target, HD 189733 b, and examine the effect of varying the planet size, coronal temperature, and impact parameter on the resulting transit profile. We also derived scaling relationships for how the overall transit detectability is affected by changing these parameters. Over most of the parameter space, we find that the detectability scales linearly with the cross-sectional area of the planet in X-rays. The relationship with increasing coronal temperature is less fixed, but averages out to a power law with slope −1/4 except when the impact parameter is high. Indeed, varying impact parameter has little effect on detectability at all until it approaches unity.

The Value-added Catalog of OB Stars in LAMOST DR7

In this work, we update the catalog of OB stars based on the Large Sky Area Multi-Object Fiber Spectroscopic Telescope data release 7 and modify the OB stars’ selection criteria in spectral line indices’ space. The new catalog includes 37,778 spectra of 27,643 OB stars, of which 3827 OB stars are newly identified. The spectral subclasses of 27,643 OB stars are obtained using the automatic classification code MKCLASS. We find that the modified OB star selection criteria can better improve the completeness of late B-type stars by analyzing their spectral classification results given by MKCLASS. We also identify 3006 Be-type stars or candidates by examining the Balmer lines in their spectra and find that the frequency of our Be-type stars (10.9%) is consistent with previous results. The spatial distribution of OB stars indicates that they are mainly located in the Galactic disk. This new catalog of OB stars will provide valuable data for studying the structure and evolution of the Milky Way.

Double red giant branch and red clump features of Galactic disc stellar populations with Gaia gspspec

The bimodality of the Milky Way disc, in the form of a thick short disc and a thinner more radially extended one, encrypts the complex internal evolution of our Galaxy and its interaction with the environment. To disentangle the different competing physical processes at play in Galactic evolution, a detailed chrono-chemical-kinematical and dynamical characterisation of the disc bimodality is necessary, including high number statistics. Here, we make use of an extremely precise sub-sample of the Gaia DR3 catalogue of stellar chemo-physical parameters. The selected database is composed of 408 800 stars with a median uncertainty of 10 K, 0.03, and 0.01 dex in and respectively. The stellar parameter precision allows us to break the age-metallicity degeneracy of disc stars. For the first time, the disc bimodality in the Kiel diagram of giant stars is observed, getting rid of interstellar absortion issues. This bimodality produces double red giant branch sequences and red clump features for mono-metallicity populations. A comparison with BaSTI isochrones allows us to demonstrate that an age gap is needed to explain the evolutionary sequence separation, in agreement with previous age-metallicity relations obtained using sub-giant stars. A bimodal distribution in the stellar mass- plane is observed at constant metallicity. Finally, a selection of stars with 0.03 dex shows that the most metal-rich population in the Milky Way disc presents an important proportion of stars with ages in the range of 5-13 Gyr, in agreement with previous literature findings. This old, extremely metal-rich population is possibly a mix of migrated stars from the internal Galactic regions, and old disc stars formed before the last major merger of the Milky Way. The Gaia Kiel diagrams of disc mono-abundance stellar populations reveal a complex, non-linear age-metallicity relation crafted by internal and external processes of Galactic evolution. Their detailed analysis opens new opportunities to reconstruct the puzzle of the Milky Way disc bimodality.

North Polar Spur: Gaseous plume(s) from star-forming regions ~3-5 kpc from the Galactic Center?

We argue that the North Polar Spur (NPS) and many less prominent structures are formed by gaseous metal-rich plumes associated with star-forming regions (SFRs). The SFRs located at the tangent to the 3-5 kpc rings might be particularly relevant to the NPS. A multi-temperature mixture of gaseous components and cosmic rays rises above the Galactic disk under the action of their initial momentum and buoyancy. Eventually, the plume velocity becomes equal to that of the ambient gas, which rotates with different angular speeds than the stars in the disk. As a result, the plumes acquire characteristic bent shapes. An ad hoc model of plumes' trajectories shows an interesting resemblance to the morphology of structures seen in the radio continuum and X-rays.

NEURORETINITIS ASSOCIATED WITH COVID-19 INFECTION AND VARICELLA-ZOSTER VIRUS REACTIVATION.

To report a case of neuroretinitis as the presenting manifestation of COVID-19 infection alongside VZV reactivation, in the absence of other systemic symptoms. A case report. A previously healthy 15-year-old male presented with a one-week history of acute unilateral visual loss and clinical findings consistent with neuroretinitis. Diagnostic work-up revealed a positive nasopharyngeal swab PCR for COVID-19 and serology showed the presence of IgM and IgG antibodies for VZV. Following a comprehensive treatment regimen involving intravenous acyclovir, high-dose intravenous steroids and oral doxycycline, rapid and remarkable improvements were observed. The macular star and optic disc swelling regressed, and visual acuity improved from 20/200 to 20/20. This case offers valuable insights into the neuroretinitis associated with COVID-19 infection and VZV reactivation.

Exoplanets across galactic stellar populations with PLATO. Estimating exoplanet yields around FGK stars for the thin disk, thick disk, and stellar halo

The vast majority of exoplanet discoveries to date have occurred around stars in the solar neighbourhood, with chemical compositions comparable to that of the Sun. However, models suggest that planetary systems in different Galactic environments, with varying dynamical histories and chemical abundances, may exhibit distinct characteristics, which can help improve our understanding of planet formation processes. This study aims to assess the potential of the upcoming PLATO mission to investigate exoplanet populations around stars in diverse Galactic environments, specifically focusing on the Milky Way thin disk, thick disk, and stellar halo. We aim to quantify PLATO's ability to detect planets in each environment and determine how these observations could constrain planet formation models. Beginning with the all-sky PLATO Input Catalogue, we kinematically classified the 2.4 million FGK stars into their respective Galactic components. For the sub-sample of stars in the long-observation LOPS2 and LOPN1 PLATO fields, we estimated planet occurrence rates using the New Generation Planet Population Synthesis (NGPPS) dataset. Combining these estimates with a PLATO detection efficiency model, we predicted the expected planet yields for each Galactic environment during a nominal 2+2 year mission. Based on our analysis, PLATO is likely to detect at least 400 exoplanets around the alpha -enriched thick disk stars. The majority of those planets are expected to be super-Earths and sub-Neptunes with radii between 2 and 10 $R_ and orbital periods between 2 and 50 days, which is ideal for studying the link between the radius valley and stellar chemistry. For the metal-poor halo, PLATO is likely to detect between 1 and 80 planets with periods between 10 and 50 days, depending on the potential existence of a metallicity threshold for planet formation. The PLATO fields contain more than 3,400 potential target stars with Fe/H < -0.6, which will help improve our understanding of planets around metal-poor stars. We identified a specific target list of 47 (kinematically classified) halo stars in the high-priority, high-signal-to-noise PLATO P1 sample, offering prime opportunities in the search for planets in metal-poor environments. PLATO's unique capabilities and large field of view position it as a valuable tool for studying planet formation across the diverse Galactic environments of the Milky Way. By probing exoplanet populations around stars with a varying chemical composition, PLATO will provide helpful insights into the connection between stellar chemistry and planet formation.

Scaling relations between mass components and rotational velocity in disk galaxies

Abstract In this study, we re-evaluate the baryonic Tully–Fisher relation (BTFR) by analyzing the correlations between maximum rotational velocity and various mass components, including stellar mass, atomic hydrogen (H i) mass, baryonic mass, and dark matter mass in a sample of 141 disk galaxies from the SPARC (Spitzer Photometry and Accurate Rotation Curves) database, augmented by recent data on stellar and dark matter masses. We apply multiple statistical methods, including Monte Carlo orthogonal distance regression (MCODR), Monte Carlo least-squares (MCLS), and traditional least-squares (LS), to assess the impact of different fitting techniques on the derived scaling relations between the mass components (stellar, H i, and dark matter) and maximum rotational velocities of these galaxies. We find that the selection of statistical methods significantly influences the derived slopes and intercepts the relation between maximum rotational velocity and mass components. The MCODR method that accounts for errors in both variables consistently produces steeper slopes, suggesting a stronger correlation between stellar mass and rotational dynamics compared to other methods. In contrast, the MCLS method tends to yield flatter slopes, highlighting the sensitivity of this approach to outliers. Despite the variations in slope and intercept across different methods, the fundamental relation between baryonic mass and rotational velocity remains robust. We have also compared dark matter mass derived from different halo models [NFW (Navarro–Frenk–White) versus combined NFW + Dekel–Zhao profiles] and noted that the slope from the NFW profile is slightly steeper than that from the combined profile, highlighting the sensitivity of scaling relations to the selection of halo model. Overall, this study reinforces the robustness of the BTFR across different mass components in disk galaxies while emphasizing the critical role of statistical methods and dark matter profiles in analyzing galactic dynamics.

Disc distortion revisited

Abstract We revisit the dynamics of razor-thin, stone-cold, self-gravitating discs. By recasting the equations into standard cylindrical coordinates, the linearised vertical dynamics of an exponential disc can be followed for several gigayears on a laptop in a few minutes. An initially warped disc rapidly evolves into a flat inner region and an outward-propagating spiral corrugation wave that rapidly winds up and would quickly thicken a disc with non-zero radial velocity dispersion. The Sgr dwarf galaxy generates a similar warp in the Galactic disc as it passes through pericentre, and the warp generated by the dwarf’s last pericentre ∼35 Myr ago is remarkably similar to the warp traced by the Galaxy’s HI disc. The resemblance to the observed warp is fleeting but its timing is perfect. For the adopted parameters the amplitude of the model warp is a factor 3 too small, but there are several reasons for this being so. The marked flaring of our Galaxy’s low-α disc just outside the solar circle can be explained as a legacy of earlier pericentres.

Radial Wave in the Galactic Disk: New Clues to Discriminate Different Perturbations

Decoding the key dynamical processes that shape the Galactic disk structure is crucial for reconstructing the Milky Way’s evolution history. The second Gaia data release unveils a novel wave pattern in the L Z −〈V R 〉 space, but its formation mechanism remains elusive due to the intricate nature of involved perturbations and the challenges in disentangling their effects. Utilizing the latest Gaia DR3 data, we find that the L Z −〈V R 〉 wave systematically shifts toward lower L Z for dynamically hotter stars with larger J Z values. The amplitude of this phase shift between stars of different dynamical hotness (ΔL Z ) peaks at around 2100 km s−1 kpc. To differentiate the role of different perturbations, we perform three sets of test particle simulations, wherein a satellite galaxy, transient spiral arms, and a bar plus the transient spiral arms act as the sole perturber, respectively. Under the satellite impact, the phase shift amplitude ΔL Z decreases toward higher L Z , which we interpret through a toy model of radial phase mixing. While neither the transient spiral arms nor the bar generates an azimuthally universal phase shift variation pattern, combining the bar and spirals generates a characteristic ΔL Z peak at the 2:1 outer Lindblad resonance (OLR) of the bar, qualitatively resembling the observed feature. Therefore, the L Z −〈V R 〉 wave is more likely of internal origin. Furthermore, linking the ΔL Z peak to the 2:1 OLR offers a novel approach to constraining the pattern speed of the Galactic bar, supporting the long/slow bar model.

Gaia/GSP-spec spectroscopic properties of γ Doradus pulsators

Context. The third Data Release of the ESA Gaia mission has provided a large sample of new gravity-mode pulsators, among which more than 11 600 are γ Dor stars. Aims. The goal of the present work is to present the spectroscopic parameters of these γ Dor pulsators estimated by the GSP-Spec module that analysed millions of Gaia spectra. Such a parametrisation could help confirm their γ Dor nature and provide their chemo-physical properties. Methods. The Galactic positions, kinematics, and orbital properties of these new Gaia pulsators were examined in order to define a sub-sample belonging to the Milky Way thin disc, in which these young stars should preferentially be found. The stellar luminosities, radii, and astrometric surface gravities were estimated without adopting any priors from uncertain stellar evolution models. These parameters, combined with the GSP-Spec effective temperatures, spectroscopic gravities, and metallicities were then validated by comparison with recent literature studies. Results. Most stars are found to belong to the Galactic thin disc, as expected. It is also found that the derived luminosities, radii, and astrometric surface gravities are high quality and have values typical of genuine γ Dor pulsators. Moreover, we show that Teff and [M/H] of pulsators with high enough S/N spectra or slow to moderate rotation rates are robust. This allowed to define a sub-sample of genuine slow-rotating Gaiaγ Dor pulsators. Their Teff were found to be between ∼6500 and ∼7800 K, log(g) is around 4.2, and the luminosities and stellar radii peak at ∼5 L⊙ and ∼1.7 R⊙, The median metallicity is close to the Solar value, although γ Dor with higher and lower metallicities by about ±0.5 dex were also identified. The [α/Fe] content is fully consistent with the chemical properties of the Galactic disc. Conclusions.Gaia/DR3 spectroscopic properties of γ Dor stars therefore confirm the nature of these pulsators and allow to chemo-physically parametrise a new large sample of such stars. Moreover, future Gaia data releases should drastically increase the number of γ Dor stars with parameters spectroscopically derived with good precision.

The Stellar Bar–Dark Matter Halo Connection in the TNG50 Simulations

Stellar bars in disk galaxies grow as stars in near-circular orbits lose angular momentum to their environments, including their dark matter (DM) halo, and transform into elongated bar orbits. This angular momentum exchange during galaxy evolution hints at a connection between bar properties and the DM halo spin λ, the dimensionless form of DM angular momentum. We investigate the connection between halo spin λ and galaxy properties in the presence/absence of stellar bars, using the cosmological magnetohydrodynamic TNG50 simulations at multiple redshifts (0 < z r < 1). We determine the bar strength (or bar amplitude, A 2/A 0), using Fourier decomposition of the face-on stellar density distribution. We determine the halo spin for barred and unbarred galaxies (0 < A 2/A 0 < 0.7) in the center of the DM halo, close to the galaxy’s stellar disk. At z r = 0, there is an anticorrelation between halo spin and bar strength. Strongly barred galaxies (A 2/A 0 > 0.4) reside in DM halos with low spin and low specific angular momentum at their centers. In contrast, unbarred/weakly barred galaxies (A 2/A 0 < 0.2) exist in halos with higher central spin and higher specific angular momentum. The anticorrelation is due to the barred galaxies’ higher DM mass and lower angular momentum than the unbarred galaxies at z r = 0, as a result of galaxy evolution. At high redshifts (z r = 1), all galaxies have higher halo spin compared to those at lower redshifts (z r = 0), with a weak anticorrelation for galaxies having A 2/A 0 > 0.2. The formation of DM bars in strongly barred systems highlights how angular momentum transfer to the halo can influence its central spin.

Disk Turbulence and Star Formation Regulation in High-z Main-sequence Analog Galaxies

The gas-phase velocity dispersions in disk galaxies, which trace turbulence in the interstellar medium, are observed to increase with lookback time. However, the mechanisms that set this rise in turbulence are observationally poorly constrained. To address this, we combine kiloparsec-scale Atacama Large Millimeter/submillimeter Array observations of CO(3−2) and CO(4−3) with Hubble Space Telescope observations of Hα to characterize the molecular gas and star formation properties of seven local analogs of main-sequence galaxies at z ∼ 1–2, drawn from the DYNAMO sample. Investigating the “molecular gas main sequence” on kiloparsec scales, we find that galaxies in our sample are more gas-rich than local star-forming galaxies at all disk positions. We measure beam-smearing-corrected molecular gas velocity dispersions and relate them to the molecular gas and star formation rate surface densities. Despite being relatively nearby (z ∼ 0.1), DYNAMO galaxies exhibit high velocity dispersions and gas and star formation rate surface densities throughout their disks, when compared to local star-forming samples. Comparing these measurements to predictions from star formation theory, we find very good agreements with the latest feedback-regulated star formation models. However, we find that theories that combine dissipation of gravitational energy from radial gas transport with feedback overestimate the observed molecular gas velocity dispersions.

PHANGS-MeerKAT and MHONGOOSE HI observations of nearby spiral galaxies: Physical drivers of the molecular gas fraction, Rmol

The molecular-to-atomic gas ratio is crucial to our understanding of the evolution of the interstellar medium (ISM) in galaxies. We investigated the balance between the atomic (ΣHI) and molecular gas (ΣH2) surface densities in eight nearby star-forming galaxies using new high-quality observations from MeerKAT and ALMA (for H I and CO, respectively). We defined the molecular gas ratio as Rmol = ΣH2/ΣHI and measured how Rmol depends on local conditions in the galaxy disks using multiwavelength observations. We find that, depending on the galaxy, H I is detected at &gt; 3σ out to 20 − 120 kpc in galactocentric radius (rgal). The typical radius at which ΣHI reaches 1 M⊙ pc−2 is rH I ≈ 22 kpc, which corresponds to 1 − 3 times the optical radius (r25). We note that, Rmol correlates best with the dynamical equilibrium pressure, PDE, among potential drivers studied, with a median correlation coefficient of ⟨ρ⟩ = 0.89. Correlations between Rmol and the star formation rate surface density, total gas surface density, stellar surface density, metallicity, and ΣSFR/PDE (a proxy for the combined effect of the UV radiation field and number density) are present but somewhat weaker. Our results also show a direct correlation between PDE and ΣSFR, supporting self-regulation models. Quantitatively, we measured similar scalings as previous works, and attribute the modest differences that we do find to the effect of varying resolution and sensitivity. At rgal ≳ 0.4r25, atomic gas dominates over molecular gas among our studied galaxies, and at the balance of these two gas phases (Rmol = 1), we find that the baryon mass is dominated by stars, with Σ* &gt; 5 Σgas. Our study constitutes an important step in the statistical investigation of how local galaxy properties (stellar mass, star formation rate, or morphology) impact the conversion from atomic to molecular gas in nearby galaxies.

Resonant Friction on Disks in Galactic Nuclei

We argue that resonant friction has a dramatic effect on a disk whose rotation direction is misaligned with that of its host nuclear star cluster. The disk’s gravity causes gravitational perturbation of the cluster that in turn exerts a strong torque back onto the disk. We argue that this torque may be responsible for the observed disruption of the clockwise disk of young stars in the Galactic center, and show in numerical experiments that it produces the observed features in the distribution of the stars’ angular momenta. More generally, we speculate that the rotation of nuclear star clusters has a stabilizing effect on the orientation of transient massive accretion disks around the supermassive black holes residing in their centers, and thus on the directions and magnitudes of the black hole spins.

Propagation of Gamma-Ray Burst Relativistic Jets in Active Galactic Nucleus Disks and Its Implication for Gamma-Ray Burst Detection

The accretion disks of supermassive black holes (SMBHs) harboring in active galactic nuclei (AGN) are considered to be an ideal site for producing different types of gamma-ray bursts (GRBs). The detectability of these GRB phenomena hidden in AGN disks is highly dependent on the dynamical evolution of the GRB relativistic jets. By investigating the reverse- and forward-shock dynamics due to the interaction between the jets and AGN disk material, we find that the relativistic jets can successfully break out from the disks only for a sufficiently high luminosity and a long enough duration. In comparison, relatively normal GRB jets are inclined to be choked in the disks unless the GRBs occur near an SMBH with relatively low mass (e.g., ∼106 M ⊙). For the choked jets, unlike normal GRB prompt and afterglow emission, we can only expect to detect emission from the forward shock when the shock is very close to the edge of the disks, i.e., at the shock breakout emission and subsequent cooling of the shock.

Bursty Star Formation in Dwarfs is Sensitive to Numerical Choices in Supernova Feedback Models

Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed, and the burstiness of the star formation history, are highly sensitive to the numerical choices adopted in these subgrid models. Using the SMUGGLE stellar feedback model, we carry out idealized simulations of an M vir ∼ 1010 M ⊙ dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of ∼3. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies, and a retention of cuspy dark matter density profiles. An improved understanding of the nonlinear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment.

EDIG-CHANGES

Context. The kinematic information of the extraplanar diffuse ionized gas (eDIG) around galaxies provides clues to the origin of the gas. Aims. The eDIG-CHANGES project studies the physical and kinematic properties of the eDIG around the CHANG-ES sample of nearby edge-on disk galaxies. Methods. We use a novel multi-slit narrow-band spectroscopy technique to obtain the spatial distribution of the spectral properties of the ionized gas around NGC 891, which is often regarded as an analogue of the Milky Way. We developed specific data reduction procedures for the multi-slit narrow-band spectroscopy data taken with the MDM 2.4 m telescope. The data presented in this paper cover the Hα and [N II]λλ6548, 6583 Å emission lines. Results. The eDIG traced by the Hα and [N II] lines shows an obvious asymmetric morphology, being brighter in the northeastern part of the galactic disk and extending a few kiloparsecs above and below the disk. Global variations in the [N II]/Hα line ratio suggest additional heating mechanisms for the eDIG at large heights beyond photoionization. We also construct position-velocity (PV) diagrams of the eDIG based on our optical multi-slit spectroscopy data and compare them to similar PV diagrams constructed with the H I data. The dynamics of the two gas phases are generally consistent with each other. Modelling the rotation curves at different heights from the galactic mid-plane suggests a vertical negative gradient in turnover radius and maximum rotation velocity, with magnitudes of approximately 3 kpc kpc−1 and 22 − 25 km s−1 kpc−1, respectively. Conclusions. Our measured vertical gradients of the rotation curve parameters suggest significant differential rotation of the ionized gas in the halo, often referred to as the lagging eDIG. Systematic study of the lagging eDIG, using the multi-slit narrow-band spectroscopy technique developed in our eDIG-CHANGES project, will help us to better understand the dynamics of the ionized gas in the halo.

Contemporaneous high-angular-resolution imaging of the AGB star W Hya in vibrationally excited H2O lines and visible polarized light with ALMA and VLT/SPHERE-ZIMPOL

Aims. We present contemporaneous high-angular-resolution millimeter imaging and visible polarimetric imaging of the nearby asymptotic giant branch (AGB) star W Hya to better understand the dynamics and dust formation within a few stellar radii. Methods. The star W Hya was observed in two vibrationally excited H2O lines at 268 and 251 GHz with Atacama Large Millimeter/submillimeter Array (ALMA) at a spatial resolution of 16 × 20 mas and at 748 and 820 nm at a resolution of 26 × 27 mas with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE)-Zurich Imaging Polarimeter (ZIMPOL). Results. ALMA’s high spatial resolution allowed us to image strong emission of the vibrationally excited H2O line at 268 GHz (v2 = 2, JKa, Kc = 65, 2–74, 3) over the stellar surface instead of absorption against the continuum, which is expected for thermal excitation. Strong, spotty emission was also detected along and just outside the stellar disk limb at an angular distance of ∼40 mas (∼1.9 R⋆), extending to ∼60 mas (∼2.9 R⋆). Another H2O line (v2 = 2, JKa, Kc = 92, 8–83, 5) at 251 GHz with a similar upper-level energy was tentatively identified, which shows absorption over the stellar surface. This suggests that the emission over the surface seen in the 268 GHz H2O line is suprathermal or even maser emission. The estimated gas temperature and H2O density are consistent with the radiatively pumped masers. The 268 GHz H2O line reveals global infall at up to ∼15 km s−1 within 2–3 R⋆, but outflows at up to ∼8 km s−1 are also present. The polarized intensity maps obtained in the visible reveal clumpy dust clouds forming within ∼40 mas (∼1.9 R⋆) with a particularly prominent cloud in the SW quadrant and a weaker cloud in the east. The 268 GHz H2O emission overlaps very well with the visible polarized intensity maps, which suggests that both the nonthermal and likely maser H2O emission and the dust originate from dense, cool pockets in the inhomogeneous atmosphere within ∼2–3 R⋆.

When LAMOST meets Gaia DR3. Exploring the metallicity of open clusters

Open clusters (OCs) are excellent probes as their age and abundance can be tightly constrained, allowing us to explore the distribution of metallicity and composition across the disk of the Milky Way. By conducting a comprehensive analysis of the metallicity of OCs, we can obtain valuable information about the history of their chemical enrichment. Moreover, by observing stars in different regions of the Milky Way, we can identify significant spatial structures in their chemical composition and abundance. This enables us to understand stellar birth radii through chemical tagging. Nevertheless, it remains challenging to infer the original positions of OCs using current data alone. The aim of this study is to investigate the distribution of metallicity in the solar neighborhood using a large dataset from Gaia DR3 combined with LAMOST spectra. With accurate ages and metallicity measurements, we can determine birth radii for the stars and attempt to understand their migration pattern. We chose a total of 1131 OCs within 3 Kpc of the Sun from the Gaia DR3 and LAMOST DR8 low-resolution spectral database (R=1800). We used an artificial neural network to correct the LAMOST data by incorporating high-resolution spectral data from GALAH DR3 (R=28000). The average metallicity of the OCs was determined based on the reliable Fe/H values for their members. We then examined the distribution of metallicity across different regions within the Galaxy and inferred birth radii of the OCs from their age and metallicity. The correction method presented here can partially eliminate the systematic offset for LAMOST data. We discuss the metallicity trend as a function of Galactocentric distance and the guiding radii. We also compare these observational results with those from chemo-dynamic simulations. Values derived from observational metallicity data are slightly lower than predicted values when the uncertainties are not considered. However, the metallicity gradients are consistent with previous calculations. Finally, we investigated the birthplace of OCs and find hints that the majority of OCs near the Sun have migrated from the outer Galactic disk.