Radial Surface Density Profiles Research Articles

Disk Assembly of the Milky Way Suggested from the Time-resolved Chemical Abundance

Both simulations and observations suggest that the disk assembly of galaxies is governed by the interplay between coplanar gas inflow, ex-planar gas outflow, and in situ star formation on the disk, known as the leaky accretion disk. This scenario predicts a strong connection between radial distributions of star formation and chemical abundances. The Milky Way, being the sole Galaxy where we can reliably measure star formation histories and the corresponding temporally resolved chemical abundances with individual stars, provides a unique opportunity to scrutinize this scenario. Based on the recent large spectroscopic and photometric surveys of Milky Way stars, we obtain the radial profiles of magnesium abundance ([Mg/H]) and star formation rate surface density at different lookback times. We find the radial profiles of [Mg/H] can be well-reproduced using the leaky accretion disk model with only two free parameters for stars formed within 4 Gyr, as well as the flattening at large radii of metallicity profiles traced by H ii regions and Cepheids. Furthermore, the constraint effective yield of the Milky Way and nearby galaxies shows broad consistency with the theoretical predictions from the stellar chemical evolution model with a mass-loading factor of 0–2. These results support that the recent assembly of the Milky Way adheres to the leaky accretion disk scenario, bridging the disk formation of our home Galaxy to the big picture of disk formation in the Universe.

A statistical and multiwavelength photometric analysis of a young embedded open star cluster: IC 1590

ABSTRACT We present a statistical and multiwavelength photometric studies of young open cluster IC 1590. We identified 91 cluster members using Gaia DR3 astrometry data using ensemble-based unsupervised machine learning techniques. From Gaia EDR3 data, we estimate the best-fitting parameters for IC 1590 using the Automated Stellar Cluster Analysis package (asteca) yielding the distance d ∼ 2.87 ± 0.02 kpc, age ∼ 3.54 ± 0.05 Myr, metallicity z ∼ 0.0212 ± 0.003, binarity value of ∼ 0.558, and extinction Av ∼ 1.252 ± 0.4 mag for an Rv value of ∼ 3.322 ± 0.23. We estimate the initial mass function slope of the cluster to be α = 1.081 ± 0.112 for single stars and α = 1.490 ± 0.051 for a binary fraction of ∼ 0.558 in the mass range 1 M⊙ ≤ m (M⊙) ≤ 100 M⊙. The G-band luminosity function slope is estimated to be ∼ 0.33 ± 0.09. We use (J − H) versus (H − Ks) colour–colour diagram to identify young stellar objects (YSOs). We found that all the identified YSOs have ages ≤ 2 Myr and masses ∼ 0.35 – 5.5 M⊙. We also fit the radial surface density profile. Using the galpy, we performed orbit analysis of the cluster. The extinction map for the cluster region has been generated using the PNICER technique, and it is almost similar to the dust structure obtained from the 500 μm dust continuum emissions map of Herschel SPIRE. We finally at the end discussed the star formation scenario in the cluster region.

The effect of sculpting planets on the steepness of debris-disc inner edges

ABSTRACT Debris discs are our best means to probe the outer regions of planetary systems. Many studies assume that planets lie at the inner edges of debris discs, akin to Neptune and the Kuiper Belt, and use the disc morphologies to constrain those otherwise-undetectable planets. However, this produces a degeneracy in planet mass and semimajor axis. We investigate the effect of a sculpting planet on the radial surface-density profile at the disc inner edge, and show that this degeneracy can be broken by considering the steepness of the edge profile. Like previous studies, we show that a planet on a circular orbit ejects unstable debris and excites surviving material through mean-motion resonances. For a non-migrating, circular-orbit planet, in the case where collisions are negligible, the steepness of the disc inner edge depends on the planet-to-star mass ratio and the initial-disc excitation level. We provide a simple analytic model to infer planet properties from the steepness of ALMA-resolved disc edges. We also perform a collisional analysis, showing that a purely planet-sculpted disc would be distinguishable from a purely collisional disc and that, whilst collisions flatten planet-sculpted edges, they are unlikely to fully erase a planet’s signature. Finally, we apply our results to ALMA-resolved debris discs and show that, whilst many inner edges are too steep to be explained by collisions alone, they are too flat to arise through completed sculpting by non-migrating, circular-orbit planets. We discuss implications of this for the architectures, histories, and dynamics in the outer regions of planetary systems.

The Global Structure of Molecular Clouds. I. Trends with Mass and Star Formation Rate

We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby clouds and find solutions that fit each of the observed radial surface density profiles remarkably well. Our most salient findings are that clouds with higher central volume densities are more compact and also have lower total mass. These same lower-mass clouds tend to have shorter gas depletion times, regardless of whether we consider their total mass or dense mass. Our analyses lead us to conclude that cylindrical clouds can be characterized by a universal structure that sets the timescale on which they form stars.

Periodic X-ray sources in the massive globular cluster 47 Tucanae: Evidence for dynamically formed cataclysmic variables

ABSTRACT We present a systematic study of periodic X-ray sources in the massive globular cluster 47 Tuc, utilizing deep archival Chandra observations that resolve the cluster core and recently available eROSITA observations that cover the cluster outskirt. By applying the Gregory-Loredo algorithm, we detect 20 periodic signals among 18 X-ray sources, ranging between 205–95731 s. Fourteen periods are newly discovered in the X-ray band. We classify these periodic sources into four quiescent low-mass X-ray binaries, 1 ms pulsar, two coronally-active binaries, and eleven cataclysmic variables (CVs) based on their X-ray temporal and spectral properties, as well as multiband information. Despite a small sample subject to potential selection bias against faint and non-magnetic CVs, the 11 CVs together define an orbital period distribution significantly different from that of the CVs previously found in the solar neighbourhood and the Galactic bulge. In particular, there exists in 47 Tuc an apparent paucity of short-period CVs below the period gap, which might be attributed to a high occupation fraction of non-magnetic CVs. Also characteristic of the 47 Tuc CVs are an overabundance of long-period CVs with a subgiant donor, a substantial fraction of CVs within the period gap, and a steep radial surface density profile. These are best understood as a group of CVs having recently formed via dynamical interactions in the dense cluster core. Despite sufficient sensitivity of the X-ray data, only one periodic source is found between one-third of the half-light radius and the tidal radius, the nature of which is unclear.

Radial profiles of surface density in debris discs

ABSTRACT Resolved observations of debris discs can be used to derive radial profiles of azimuthally averaged surface density (ASD), which carries important information about the disc structure even in the presence of non-axisymmetric features and has improved signal-to-noise characteristics relative to the two-dimensional surface brightness measurements. We develop a (semi-)analytical formalism allowing one to relate ASD to the underlying semimajor axis and eccentricity distributions of the debris particles in a straightforward manner. This approach does not involve the distribution of particle apsidal angles, thus simplifying calculations. It is a much faster, more flexible and effective way of calculating ASD than the Monte Carlo sampling of orbital parameters of debris particles. We present explicit analytical results based on this technique for a number of particle eccentricity distributions, including two cases of particular practical importance: a prescribed radial profile of eccentricity, and the Rayleigh distribution of eccentricities. We then show how our framework can be applied to observations of debris discs and rings for retrieving either the semimajor axis distribution or (in some cases) the eccentricity distribution of debris, thus providing direct information about the architecture and dynamical processes operating in debris discs. Our approach also provides a fast and efficient way of forward modelling observations. Applications of this technique to other astrophysical systems, e.g. the nuclear stellar disc in M31 or tenuous planetary rings, are also discussed.

Towards a fully consistent Milky Way disk model

Context. The semi-analytic Just-Jahreiß (JJ) model of the Galactic disk is a flexible tool for stellar population synthesis with a fine age resolution of 25 Myr. The model has recently been calibrated in the solar neighbourhood against the Gaia DR2 stars. We have identified two star-formation bursts within the last ∼4 Gyr of the local star-formation rate (SFR) history. Aims. In this work we present a generalised version of the JJ model that incorporates our findings for the solar neighbourhood and is applicable to a wide range of galactocentric distances, 4 kpc ≲R ≲ 14 kpc. Methods. The JJ model includes the four flattened and two spheroidal components of the Milky Way, describing it as an axisymmetric system. The thin and thick disks, as well as atomic and molecular gas layers, are assumed to have exponential radial surface density profiles. Spherical stellar halo and dark matter in the form of a cored isothermal sphere are also added to the model. The overall thickness of the thin disk is assumed to be constant at all radii, though model realisations with a flaring disk can also be tested. The adopted radial variation in the thin-disk SFR reflects the inside-out disk growth scenario. Motivated by our findings for the solar neighbourhood, we allow a smooth power-law SFR continuum to be modified by an arbitrary number of Gaussian peaks. Additionally, the vertical kinematics of the stellar populations associated with these episodes of star-formation excess is allowed to differ from the kinematics prescribed by the age-velocity dispersion relation for the thin-disk populations of the same age. Results. We present a public code of the JJ model complemented by the three sets of isochrones generated by the stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code (PARSEC), the Modules and Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST), and a Bag of Stellar Tracks and Isochrones (BaSTI). Assuming a plausible set of parameters, we take the first step towards calibrating the JJ model at non-solar radii. Using metallicity distributions of the red clump giants from the Apache Point Observatory Galactic Evolution Experiment (APOGEE), we constrain the radial variation of the JJ-model age-metallicity relation and propose a new analytic form for the age-metallicity relation function. Conclusions. The generalised JJ model is a publicly available tool for studying different stellar populations across the Milky Way disk. With its fine age resolution and flexibility, it can be particularly useful for reconstructing the thin-disk SFR, as a variety of different SFR shapes can be constructed within its framework.

The globular clusters and star formation history of the isolated, quiescent ultra-diffuse galaxy DGSAT I

ABSTRACT We investigate the isolated, quiescent ultra-diffuse galaxy (UDG) DGSAT I and its globular cluster (GC) system using two orbits of Hubble Space Telescope Advanced Camera for Surveys imaging in the F606W and F814W filters. This is the first study of GCs around a UDG in a low-density environment. DGSAT I was previously found to host an irregular blue low surface brightness clump, which we confirm as very likely belonging to the galaxy rather than being a chance projection, and represents a recent episode of star formation (∼500 Myr) that challenges some UDG formation scenarios. We select GC candidates based on colours and magnitudes, and construct a self-consistent model of the GC radial surface density profile along with the background. We find a half-number radius of RGC = 2.7 ± 0.1 kpc (more compact than the diffuse starlight) and a total of 12 ± 2 GCs. The total mass fraction in GCs is relatively high, supporting an overmassive dark matter halo as also implied by the high velocity dispersion previously measured. The GCs extend to higher luminosities than expected, and have colours that are unusually similar to their host galaxy colour, with a very narrow spread, all of which suggest an early, intense burst of cluster formation. The nature and origin of this galaxy remain puzzling, but the most likely scenario is a ‘failed galaxy’ that formed relatively few stars for its halo mass, and could be related to cluster UDGs whose size and quiescence pre-date their infall.

Determining the full satellite population of a Milky Way-mass halo in a highly resolved cosmological hydrodynamic simulation

ABSTRACT We investigate the formation of the satellite galaxy population of a Milky Way-mass halo in a very highly resolved magnetohydrodynamic cosmological zoom-in simulation (baryonic mass resolution mb = 800 $\rm M_{\odot }$). We show that the properties of the central star-forming galaxy, such as the radial stellar surface density profile and star formation history, are (i) robust to stochastic variations associated with the so-called Butterfly Effect and (ii) well converged over 3.5 orders of magnitude in mass resolution. We find that there are approximately five times as many satellite galaxies at this high resolution compared to a standard ($m_b\sim 10^{4-5}\, \rm M_{\odot }$) resolution simulation of the same system. This is primarily because two-thirds of the high-resolution satellites do not form at standard resolution. A smaller fraction (one-sixth) of the satellites present at high-resolution form and disrupt at standard resolution; these objects are preferentially low-mass satellites on intermediate- to low-eccentricity orbits with impact parameters ≲30 kpc. As a result, the radial distribution of satellites becomes substantially more centrally concentrated at higher resolution, in better agreement with recent observations of satellites around Milky Way-mass haloes. Finally, we show that our galaxy formation model successfully forms ultra-faint galaxies and reproduces the stellar velocity dispersion, half-light radii, and V-band luminosities of observed Milky Way and Local Group dwarf galaxies across six orders of magnitude in luminosity (103–$10^{9}\, \rm L_{\odot }$).

Mapping the Galactic Disk with the LAMOST and Gaia Red Clump Sample. VII. The Stellar Disk Structure Revealed by the Mono-abundance Populations

Abstract Using a sample of 96,201 primary red clump stars selected from the LAMOST and Gaia surveys, we investigate the stellar structure of the Galactic disk. The sample stars show two separated sequences of high-[α/Fe] and low-[α/Fe] in the [α/Fe]–[Fe/H] plane. We divide the sample stars into five mono-abundance populations (MAPs) with different ranges of [α/Fe] and [Fe/H], named as the high-[α/Fe], high-[α/Fe] and high-[Fe/H], low-[Fe/H], solar, high-[Fe/H] MAPs, respectively. We present the stellar number density distributions in the R–Z plane, and the scale heights and scale lengths of the individual MAPs by fitting their vertical and radial density profiles. The vertical profiles, the variation trend of scale height with the Galactocentric radius, indicate that there is a clear disk flare in the outer disk both for the low-[α/Fe] and the high-[α/Fe] MAPs. While the radial surface-density profiles show a peak radius of 7 kpc and 8 kpc for the high-[α/Fe] and low-[α/Fe] MAPs, respectively. We also investigate the correlation between the mean rotation velocity and metallicity of the individual MAPs, and find that the mean rotation velocities are well separated and show different trends between the high-[α/Fe] and the low-[α/Fe] MAPs. Finally, we discuss the character of the high-[α/Fe] and high-[Fe/H] MAP and find that it is more similar to the high-[α/Fe] MAP either in the radial and vertical density profiles or in the rotation velocity.

VINTERGATAN – I. The origins of chemically, kinematically, and structurally distinct discs in a simulated Milky Way-mass galaxy

ABSTRACT Spectroscopic surveys of the Milky Way’s stars have revealed spatial, chemical, and kinematical structures that encode its history. In this work, we study their origins using a cosmological zoom simulation, VINTERGATAN, of a Milky Way-mass disc galaxy. We find that in connection to the last major merger at z ∼ 1.5, cosmological accretion leads to the rapid formation of an outer, metal-poor, low-[α/Fe] gas disc around the inner, metal-rich galaxy containing the old high-[α/Fe] stars. This event leads to a bimodality in [α/Fe] over a range of [Fe/H]. A detailed analysis of how the galaxy evolves since z ∼ 1 is presented. We demonstrate the way in which inside-out growth shapes the radial surface density and metallicity profile and how radial migration preferentially relocates stars from the inner disc to the outer disc. Secular disc heating is found to give rise to increasing velocity dispersions and scale heights with stellar age, which together with disc flaring explains several trends observed in the Milky Way, including shallower radial [Fe/H] profiles above the mid-plane. We show how the galaxy formation scenario imprints non-trivial mappings between structural associations (i.e. thick and thin discs), velocity dispersions, α-enhancements, and ages of stars; e.g. the most metal-poor stars in the low-[α/Fe] sequence are found to have a scale height comparable to old high-[α/Fe] stars. Finally, we illustrate how at low spatial resolution, comparable to the thickness of the galaxy, the proposed pathway to distinct sequences in [α/Fe]–[Fe/H] cannot be captured.

New Determination of Fundamental Properties of Palomar 5 Using Deep DESI Imaging Data

Abstract The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of , tidal radius of , and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H] = −1.798 ± 0.014, reddening of E(B − V) = 0.0552 ± 0.0005, and distance modulus of . The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

A thick-disk galaxy model and simulations of equal-mass galaxy pair collisions

We implement a numerical model reported in the literature to simulate the evolution of a galaxy composed of four matter components, such as: a dark-matter halo; a rotating disk of stars; a spherical bulge of stars and a ring of molecular gas. We show that the evolution of this galaxy model is stable at least for 10 Gyr (Gyr = 109 years). We characterize the resulting configuration of this galaxy model by figures of the circular velocity and angular momentum distribution; the tangential and radial components of the velocity; the peak density evolution and the radial density profile. Additionally, we calculate several models of equal-mass galaxy binary collisions, such as: (i) frontal and (ii) oblique (with an impact parameter), (iii) two models with initial conditions taken from a 2-body orbit and (iv) a very close passage. To allow comparison with the galaxy model, we characterize the dynamics of the collision models in an analogous way. Finally, we determine the de Vaucouleurs fitting curves of the radial density profile, on a radial scale of 0–100 kpc, for all the collision models irrespective of the pre-collision trajectory. To study the radial mass density and radial surface density profiles at a smaller radial scale, 0–20 kpc, we use a four-parameter fitting curve.

Volumetric star formation prescriptions in vertically resolved edge-on galaxies

ABSTRACT We measure the gas disc thicknesses of the edge-on galaxy NGC 4013 and the less edge-on galaxies (NGC 4157 and 5907) using CO (CARMA/OVRO) and/or H i (EVLA) observations. We also estimate the scale heights of stars and/or the star formation rate (SFR) for our sample of five galaxies using Spitzer IR data (3.6 and 24 µm). We derive the average volume densities of the gas and the SFR using the measured scale heights along with radial surface density profiles. Using the volume density that is more physically relevant to the SFR than the surface density, we investigate the existence of a volumetric star formation law (SFL), how the volumetric SFL is different from the surface-density SFL, and how the gas pressure regulates the SFR based on our galaxy sample. We find that the volumetric and surface SFLs in terms of the total gas have significantly different slopes, while the volumetric and surface SFLs in terms of the molecular gas do not show any noticeable difference. The volumetric SFL for the total gas has a flatter power-law slope of 1.26 with a smaller scatter of 0.19 dex compared to the slope (2.05) and the scatter (0.25 dex) of the surface SFL. The molecular gas SFLs have similar slopes of 0.78 (volume density) and 0.77 (surface density) with similar rms scatters. We show that the interstellar gas pressure is strongly correlated with the SFR but find no significant difference between the correlations based on the volume and surface densities.

Dust modeling of the combined ALMA and SPHERE datasets of HD 163296

Context. Multiwavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aim to construct a three-dimensional model of HD 163296 that is capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We want to determine why the spectral energy distribution of HD 163296 is intermediary between the otherwise well-separated group I and group II Herbig stars. Methods. The disk was modeled using the Monte Carlo radiative transfer code MCMax3D. The radial dust surface density profile was modeled after the ALMA observations, while the polarized scattered light observations were used to constrain the inclination of the inner disk component and turbulence and grain growth in the outer disk. Results. While three rings are observed in the disk midplane in millimeter thermal emission at ~80, 124, and 200 AU, only the innermost of these is observed in polarized scattered light, indicating a lack of small dust grains on the surface of the outer disk. We provide two models that are capable of explaining this difference. The first model uses increased settling in the outer disk as a mechanism to bring the small dust grains on the surface of the disk closer to the midplane and into the shadow cast by the first ring. The second model uses depletion of the smallest dust grains in the outer disk as a mechanism for decreasing the optical depth at optical and near-infrared wavelengths. In the region outside the fragmentation-dominated regime, such depletion is expected from state-of-the-art dust evolution models. We studied the effect of creating an artificial inner cavity in our models, and conclude that HD 163296 might be a precursor to typical group I sources.

On the diversity and statistical properties of protostellar discs

Abstract We present results from the first population synthesis study of protostellar discs. We analyse the evolution and properties of a large sample of protostellar discs formed in a radiation hydrodynamical simulation of star cluster formation. Due to the chaotic nature of the star formation process, we find an enormous diversity of young protostellar discs, including misaligned discs, and discs whose orientations vary with time. Star–disc interactions truncate discs and produce multiple systems. Discs may be destroyed in dynamical encounters and/or through ram-pressure stripping, but reform by later gas accretion. We quantify the distributions of disc mass and radii for protostellar ages up to ≈105 yr. For low-mass protostars, disc masses tend to increase with both age and protostellar mass. Disc radii range from of order 10 to a few hundred au, grow in size on time-scales ≲ 104 yr, and are smaller around lower mass protostars. The radial surface density profiles of isolated protostellar discs are flatter than the minimum mass solar nebula model, typically scaling as Σ ∝ r−1. Disc to protostar mass ratios rarely exceed two, with a typical range of Md/M* = 0.1–1 to ages ≲ 104 yr and decreasing thereafter. We quantify the relative orientation angles of circumstellar discs and the orbit of bound pairs of protostars, finding a preference for alignment that strengths with decreasing separation. We also investigate how the orientations of the outer parts of discs differ from the protostellar and inner disc spins for isolated protostars and pairs.

The structure of young embedded protostellar discs

Young protostellar discs provide the initial conditions for planet formation. The properties of these discs may be different from those of late-phase (T Tauri) discs due to continuing infall from the envelope and protostellar variability resulting from irregular gas accretion. We use a set of hydrodynamic simulations to determine the structure of discs forming in collapsing molecular clouds. We examine how radiative feedback from the host protostar affects the disc properties by examining three regimes: without radiative feedback, with continuous radiative feedback and with episodic feedback, similar to FU Ori-type outbursts. We find that the radial surface density and temperature profiles vary significantly as the disc accretes gas from the infalling envelope. These profiles are sensitive to the presence of spiral structure, induced by gravitational instabilities, and the radiative feedback provided by the protostar, especially in the case when the feedback is episodic. We also investigate whether mass estimates from position-velocity (PV) diagrams are accurate for early-phase discs. We find that the protostellar system mass (i.e. the mass of the protostar and its disc) is underestimated by up to 20%, due to the impact of an enhanced radial pressure gradient on the gas. The mass of early-phase discs is a significant fraction of the mass of the protostar, so position-velocity diagrams cannot accurately provide the mass of the protostar alone. The enhanced radial pressure gradient expected in young discs may lead to an increased rate of dust depletion due to gas drag, and therefore to a reduced dust-to-gas ratio.

Red Clump Stars from LAMOST II: the outer disc of the Milky Way

We present stellar density maps of the Galactic outer disc with red clump stars from LAMOST data. These samples are separated into younger (mean age ∼ 2.7 Gyr) and older (mean age ∼ 4.6 Gyr) populations so that they can trace the variation of structures with ages in the range of Galactocentric radius R from 9 to 13.5 kpc. We show that scale heights for both of the two populations increase with R and display radial gradients of 48±6 and 40±4 pc kpc for the older and younger populations, respectively. It is evident that flaring occurs in the thin disc populations with a wide range of ages. Moreover, the intensity of flaring does not seem to be significantly related to the age of the thin disc populations. On the other hand, scale lengths of the radial surface density profiles are 4.7±0.5 kpc for the younger population and 3.4±0.2 kpc for the older population, meaning that the younger disc population is more radially extended than the older one. Although the fraction of the younger population mildly increases from 28% at R ∼ 9 to about 35% at R ∼ 13 kpc, the older population is prominent with a fraction of no less than 65% in the outer disc.

The Spatial Structure of the Galactic outer disk with LAMOST DR3 K giant stars

AbstractThe spatial structure of the Milky Way outer disk is characterized with ~ 70000 LAMOST DR3 K giants after carefully correction of their selection effects. By slicing the data into various Galactocentric radius bins, we are able to fit the vertical stellar density profile with a models composed of two isothermal-sheet disks and an oblate power-law halo. We find that although the thin disk is significantly flared, the radial surface density profile can extend to as far as 19 kpc. Beyond 12 kpc, only one thicker disk, rather than two disk components, are found in the samples. Moreover, the residual of the density profiles after subtracting the best fit models show different oscillation patterns in almost all range of detecting radius.

Radial Surface Density Profiles of Gas and Dust in the Debris Disk around 49 Ceti

Abstract We present ∼0.″4 resolution images of CO(3–2) and associated continuum emission from the gas-bearing debris disk around the nearby A star 49 Ceti, observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We analyze the ALMA visibilities in tandem with the broadband spectral energy distribution to measure the radial surface density profiles of dust and gas emission from the system. The dust surface density decreases with radius between ∼100 and 310 au, with a marginally significant enhancement of surface density at a radius of ∼110 au. The SED requires an inner disk of small grains in addition to the outer disk of larger grains resolved by ALMA. The gas disk exhibits a surface density profile that increases with radius, contrary to most previous spatially resolved observations of circumstellar gas disks. While ∼80% of the CO flux is well described by an axisymmetric power-law disk in Keplerian rotation about the central star, residuals at ∼20% of the peak flux exhibit a departure from axisymmetry suggestive of spiral arms or a warp in the gas disk. The radial extent of the gas disk (∼220 au) is smaller than that of the dust disk (∼300 au), consistent with recent observations of other gas-bearing debris disks. While there are so far only three broad debris disks with well characterized radial dust profiles at millimeter wavelengths, 49 Ceti’s disk shows a markedly different structure from two radially resolved gas-poor debris disks, implying that the physical processes generating and sculpting the gas and dust are fundamentally different.