Surface Density Profile Research Articles

Where Did the Outskirts Go? Outer Stellar Halos as a Sensitive Probe of Supernova Feedback

A recent comparison by Merritt of simulated and observed Milky Way–mass galaxies has identified a significant tension between the outskirts (r > 20 kpc) of the stellar halos in simulated and observed galaxies. Using observations from the Dragonfly telescope and simulated galaxies from the Illustris-TNG100 project, Merritt found that the outskirts of stellar halos in simulated galaxies have surface densities 1–2 dex higher than those of observed galaxies. In this paper, we compare two suites of 15 simulated Milky Way–like galaxies, each drawn from the same initial conditions, simulated with the same hydrodynamical code, but with two different models for feedback from supernovae. We find that the McMaster Unbiased Galaxy Simulations (MUGS), which use an older “delayed-cooling” model for feedback, also produce too much stellar mass in the outskirts of the halo, with median surface densities well above observational constraints. The MUGS2 simulations, which instead use a new, physically motivated “superbubble” model for stellar feedback, have 1–2 dex lower outer stellar halo masses and surface densities. The MUGS2 simulations generally match both the median surface density profile as well as the scatter in stellar halo surface density profiles seen in observed stellar halos. We conclude that there is no “missing outskirts” problem in cosmological simulations, provided that supernova feedback is modeled in a way that allows it to efficiently regulate star formation in the low-mass progenitor environments of stellar halo outskirts.

The contribution ofin situandex situstar formation in early-type galaxies: MaNGA versus IllustrisTNG

ABSTRACTWe compare stellar mass surface density, metallicity, age, and line-of-sight velocity dispersion profiles in massive ($M_*\ge 10^{10.5}\, \mathrm{M_\odot }$) present-day early-type galaxies (ETGs) from the MaNGA survey with simulated galaxies from the TNG100 simulation of the IllustrisTNG suite. We find an excellent agreement between the stellar mass surface density profiles of MaNGA and TNG100 ETGs, both in shape and normalization. Moreover, TNG100 reproduces the shapes of the profiles of stellar metallicity and age, as well as the normalization of velocity dispersion distributions of MaNGA ETGs. We generally also find good agreement when comparing the stellar profiles of central and satellite galaxies between MaNGA and TNG100. An exception is the velocity dispersion profiles of very massive ($M_*\gtrsim 10^{11.5}\, \mathrm{M_\odot }$) central galaxies, which, on average, are significantly higher in TNG100 than in MaNGA ($\approx 50\, \mathrm{km\, s^{-1}}$). We study the radial profiles of in situ and ex situ stars in TNG100 and discuss the extent to which each population contributes to the observed MaNGA profiles. Our analysis lends significant support to the idea that high-mass ($M_*\gtrsim 10^{11}\, \mathrm{M_\odot }$) ETGs in the present-day Universe are the result of a merger-driven evolution marked by major mergers that tend to homogenize the stellar populations of the progenitors in the merger remnant.

Debiasing the Minimum-mass Extrasolar Nebula: On the Diversity of Solid Disk Profiles

A foundational idea in the theory of in situ planet formation is the “minimum-mass extrasolar nebula” (MMEN), a surface density profile (Σ) of disk solids that is necessary to form the planets in their present locations. While most previous studies have fit a single power law to all exoplanets in an observed ensemble, it is unclear whether most exoplanetary systems form from a universal disk template. We use an advanced statistical model for the underlying architectures of multiplanet systems to reconstruct the MMEN. The simulated physical and Kepler-observed catalogs allow us to directly assess the role of detection biases, and in particular the effect of nontransiting or otherwise undetected planets, in altering the inferred MMEN. We find that fitting a power law of the form to each multiplanet system results in a broad distribution of disk profiles; g cm−2 and encompass the 16th–84th percentiles of the marginal distributions in an underlying population, where is the normalization at a 0 = 0.3 au. Around half of the inner planet-forming disks have minimum solid masses of ≳ 40M ⊕ within 1 au. While transit observations do not tend to bias the median β, they can lead to both significantly over- and underestimated and thus broaden the inferred distribution of disk masses. Nevertheless, detection biases cannot account for the full variance in the observed disk profiles; there is no universal MMEN if all planets formed in situ. The great diversity of solid disk profiles suggests that a substantial fraction (≳23%) of planetary systems experienced a history of migration.

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.

Hierarchical structure of YSO clusters in the W40 and Serpens South region: group extraction and comparison with fractal clusters

ABSTRACT Young stellar clusters are believed to inherit the spatial distribution like hierarchical structures of their natal molecular cloud during their formation. However, the change of the structures between the cloud and the young clusters is not well constrained observationally. We select the W40–Serpens South region (∼7 × 9 pc2) of the Aquila Rift as a testbed and investigate hierarchical properties of spatial distribution of young stellar objects (YSOs) in this region. We develop a minimum spanning tree (MST) based method to group stars into several levels by successively cutting down edges longer than an algorithmically determined critical value. A total of 832 YSOs are divided into 5 levels with 23 groups. For describing the hierarchical properties in a controlled way, we construct a set of synthetic source distributions at various fractal dimensions, and apply the same method to explore their group characters. By comparing the Q parameter and the surface density profiles of the observed and the synthetic data, we find that the YSO observation matches spatial patterns from multifractal dimensions. In the periphery region where the molecular clouds are more diffuse, the YSO structure is close to a fractal dimension of 2.0. While in the core regions, the fractal dimensions are close to 1.6 and 1.4 for the W40 and the Serpens South regions, respectively. Therefore, the YSOs may inherit the fractal pattern of the dense part of the molecular clouds, but such pattern dissipates slowly in several Myr.

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.

XGASS: the connection between angular momentum, mass, and atomic gas fraction in nearby galaxies

ABSTRACT We use a sample of 559 disc galaxies extracted from the eXtended GALEX Arecibo SDSS Survey to study the connection between baryonic angular momentum, mass, and atomic gas fraction in the local Universe. Baryonic angular momenta are determined by combining H i and H2 integrated profiles with 2D stellar mass surface density profiles. In line with previous work, we confirm that specific angular momentum and atomic gas fraction are tightly correlated, but we find a larger scatter than previously observed. This is most likely due to the wider range of galaxy properties covered by our sample. We compare our findings with the predictions of the analytical stability model developed by Obreschkow et al. and find that, while the model provides a very good first-order approximation for the connection between baryonic angular momentum, mass, and gas fraction, it does not fully match our data. Specifically, we find that at fixed baryonic mass, the dependence of specific angular momentum on gas fraction is significantly weaker, and at fixed gas fraction, the slope of the angular momentum versus mass relation is shallower than what was predicted by the model. The reasons behind this tension remain unclear, but we speculate that multiple factors may simultaneously play a role, all related to the fact that the model is not able to encapsulate the full diversity of galaxy properties in our sample.

Discovery of Line Pressure Broadening and Direct Constraint on Gas Surface Density in a Protoplanetary Disk

The gas surface density profile of protoplanetary disks is one of the most fundamental physical properties to understanding planet formation. However, it is challenging to determine the surface density profile observationally, because the H2 emission cannot be observed in low-temperature regions. We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) archival data of the 12CO J = 3 − 2 line toward the protoplanetary disk around TW Hya and discovered extremely broad line wings due to the pressure broadening. In conjunction with a previously reported optically thin CO isotopologue line, the pressure broadened line wings enabled us to directly determine the midplane gas density for the first time. The gas surface density at ∼5 au from the central star reaches ∼103 g cm−2, which suggests that the inner region of the disk has enough mass to form a Jupiter-mass planet. Additionally, the gas surface density drops at the inner cavity by ∼2 orders of magnitude compared to outside the cavity. We also found a low CO abundance of ∼10−6 with respect to H2, even inside the CO snow line, which suggests conversion of CO to less volatile species. Combining our results with previous studies, the gas surface density jumps at r ∼ 20 au, suggesting that the inner region (3 < r < 20 au) might be the magnetorotational instability dead zone. This study sheds light on the direct gas surface density constraint without assuming the CO/H2 ratio using ALMA.

Metamorphism of snow on Arctic sea ice during the melt season: impact on spectral albedo and radiative fluxes through snow

Abstract. The energy budget of Arctic sea ice is strongly affected by the snow cover. Intensive sampling of snow properties was conducted near Qikiqtarjuaq in Baffin Bay on typical landfast sea ice during two melt seasons in 2015 and 2016. The sampling included stratigraphy, vertical profiles of snow specific surface area (SSA), density and irradiance, and spectral albedo (300–1100 nm). Both years featured four main phases: (I) dry snow cover, (II) surface melting, (III) ripe snowpack, and (IV) melt pond formation. Each phase was characterized by distinctive physical and optical properties. A high SSA value of 49.3 m2 kg−1 was measured during phase I on surface wind slabs together with a corresponding broadband albedo (300–3000 nm) of 0.87. Phase II was marked by alternating episodes of surface melting, which dramatically decreased the SSA below 3 m2 kg−1, and episodes of snowfall re-establishing pre-melt conditions. Albedo was highly time-variable, with minimum broadband values around 0.70. In phase III, continued melting led to a fully ripe snowpack composed of clustered rounded grains. Albedo began to decrease in the visible as snow thickness decreased but remained steady at longer wavelengths. Moreover, significant spatial variability appeared for the first time following snow depth heterogeneity. Spectral albedo was simulated by radiative transfer using measured SSA and density vertical profiles and estimated impurity contents based on limited measurements. Simulations were most of the time within 1 % of measurements in the visible and within 2 % in the infrared. Simulations allowed the calculations of albedo and of the spectral flux at the snow–ice interface. These showed that photosynthetically active radiation fluxes at the bottom of the snowpack durably exceeded 5 W m−2 (∼9.2 µmol m−2 s−1) only when the snowpack thickness started to decrease at the end of phase II.

Molecular Dynamics Study on Properties of Hydration Layers above Polymer Antifouling Membranes.

Zwitterionic polymers as crucial antifouling materials exhibit excellent antifouling performance due to their strong hydration ability. The structure–property relationship at the molecular level still remains to be elucidated. In this work, the surface hydration ability of three antifouling polymer membranes grafting on polysiloxane membranes Poly(sulfobetaine methacrylate) (T4-SB), poly(3-(methacryloyloxy)propane-1-sulfonate) (T4-SP), and poly(2-(dimethylamino)ethyl methacrylate) (T4-DM) was investigated. An orderly packed, and tightly bound surface hydration layer above T4-SP and T4-SB antifouling membranes was found by means of analyzing the dipole orientation distribution, diffusion coefficient, and average residence time. To further understand the surface hydration ability of three antifouling membranes, the surface structure, density profile, roughness, and area percentage of hydrophilic surface combining electrostatic potential, RDFs, SDFs, and noncovalent interactions of three polymers’ monomers were studied. It was concluded that the broadest distribution of electrostatic potential on the surface and the nature of anionic SO3- groups led to the following antifouling order of T4-SB > T4-SP > T4-DM. We hope that this work will gain some insight for the rational design and optimization of ecofriendly antifouling materials.

Spreading pressure bumps in gas-dust discs can stall planet migration via planet-vortex interactions

ABSTRACT We investigate the gravitational interaction between low- to intermediate-mass planets ($M_p \in [0.06-210]\, \mathrm{ M}_{\oplus }$) and two previously formed pressure bumps in a gas-dust protoplanetary disc. We explore how the disc structure changes due to planet-induced perturbations and also how the appearance of vortices affects planet migration. We use multifluid 2D hydrodynamical simulations and the dust is treated in the pressureless-fluid approximation, assuming a single grain size of $5\, \mu {\mathrm{m}}$. The initial surface density profiles containing two bumps are motivated by recent observations of the protoplanetary disc HD163296. When planets are allowed to migrate, either a single planet from the outer pressure maximum or two planets from each pressure maximum, the initial pressure bumps quickly spread and merge into a single bump which is radially wide and has a very low amplitude. The redistribution of the disc material is accompanied by the Rossby Wave Instability and an appearance of mini-vortices that merge in a short period of time to form a large vortex. The large vortex induces perturbations with a spiral wave pattern that propagate away from the vortex as density waves. We found that these vortex-induced spiral waves strongly interact with the spiral waves generated by the planet and we called this mechanism the ‘Faraway Interaction’. It facilitates much slower and/or stagnant migration of the planets and it excites their orbital eccentricities in some cases. Our study provides a new explanation for how rocky planets can come to have a slow migration in protoplanetary discs where vortex formation occurs.

Forming young and hypervelocity stars in the Galactic Centre via tidal disruption of a molecular cloud

ABSTRACT The Milky Way Galaxy hosts a four million solar mass black hole, Sgr A*, that underwent a major accretion episode approximately 3–6 Myr ago. During the episode, hundreds of young massive stars formed in a disc orbiting Sgr A* in the central half parsec. The recent discovery of a hypervelocity star (HVS) S5-HVS1, ejected by Sgr A* five Myr ago with a velocity vector consistent with the disc, suggests that this event also produced binary star disruptions. The initial stellar disc has to be rather eccentric for this to occur. Such eccentric discs can form from the tidal disruptions of molecular clouds. Here, we perform simulations of such disruptions, focusing on gas clouds on rather radial initial orbits. As a result, stars formed in our simulations are on very eccentric orbits ($\bar{e}\sim 0.6$) with a lopsided configuration. For some clouds, counterrotating stars are formed. As in previous work, we find that such discs undergo a secular gravitational instability that leads to a moderate number of particles obtaining eccentricities of 0.99 or greater, sufficient for stellar binary disruption. We also reproduce the mean eccentricity of the young disc in the Galactic Centre, though not the observed surface density profile. We discuss missing physics and observational biases that may explain this discrepancy. We conclude that observed S-stars, HVSs, and disc stars tightly constrain the initial cloud parameters, indicating a cloud mass between a few × 104 and $10^5\, {\rm M}_{\odot }$, and a velocity between ∼40 and 80 km s−1 at 10 pc.

AMICO galaxy clusters in KiDS-DR3: The impact of estimator statistics on the luminosity-mass scaling relation

Context. As modern-day precision cosmology aims for statistical uncertainties of the percent level or lower, it becomes increasingly important to reconsider estimator assumptions at each step of the process, along with their consequences on the statistical variability of the scientific results. Aims. We compare L1 regression statistics to the weighted mean, the canonical L2 method based on Gaussian assumptions, to infer the weak gravitational shear signal from a catalog of background ellipticity measurements around a sample of clusters, which has been a standard step in the processes of many recent analyses. Methods. We use the shape measurements of background sources around 6925 AMICO clusters detected in the KiDS third data release. We investigate the robustness of our results and the dependence of uncertainties on the signal-to-noise ratios of the background source detections. Using a halo model approach, we derive lensing masses from the estimated excess surface density profiles. Results. The highly significant shear signal allows us to study the scaling relation between the r-band cluster luminosity, L200, and the derived lensing mass, M200. We show the results of the scaling relations derived in 13 bins in L200, with a tightly constrained power-law slope of ∼1.24 ± 0.08. We observe a small, but significant, relative bias of a few percent in the recovered excess surface density profiles between the two regression methods, which translates to a 1σ difference in M200. The efficiency of L1 is at least that of the weighted mean and increases with higher signal-to-noise shape measurements. Concluions. Our results indicate the relevance of optimizing the estimator for inferring the gravitational shear from a distribution of background ellipticities. The interpretation of measured relative biases can be gauged by deeper observations, and the increased computation times remain feasible.

NProFit: A Tool for Fitting the Surface Brightness Profiles of Star Clusters with Dynamical Models

The surface brightness profiles (SBPs) of star clusters hold invaluable information on the dynamical state of clusters. The observed SBPs of star clusters, especially that of globular clusters, are in good agreement with the SBPs expected for isothermal spheres containing stars of reduced kinetic energies. However, the SBPs of configurations that satisfy these theoretical criteria cannot be uniquely expressed by analytical formulae, which had hindered the analysis of dynamical state of observed clusters in external galaxies. To counter this shortcoming, it has become a practice to use empirical fitting formulae that best represent the core and halo characteristics of theoretical models. We here present a general purpose code, named nProFit, that allows fitting of the surface brightness profiles of extragalactic star clusters to theoretical star clusters, defined by dynamical models of King and Wilson. In addition, we also incorporated theoretical models that result in power-law surface brightness profiles represented by Elson et al. The code returns the basic size parameters such as core radius, half-light radius and tidal radius, as well as dynamically relevant parameters, such as the volume and surface density profiles, velocity dispersion profile, total mass and the binding energy for a user-fixed mass-to-light ratio. The usefulness of the code in the dynamical study of extragalactic clusters has been already illustrated in Cuevas-Otahola et al. The code, which is python-based at the user end, but makes calls to advanced routines in Pyraf and Fortran, is now available for public use. We provide example scripts and mock clusters in the installation package as guide to users.

Exploiting non-linear scales in galaxy–galaxy lensing and galaxy clustering: A forecast for the dark energy survey

ABSTRACT The combination of galaxy–galaxy lensing (GGL) and galaxy clustering is a powerful probe of low-redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we use an N-body and halo occupation distribution (HOD) emulator method to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias Alens. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile ΔΣ(rp) and the projected galaxy correlation function wp, gg(rp) in the final (Year 6) DES data set over scales $r_p=0.3\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$. For a $3{{\ \rm per\ cent}}$ prior on Alens we forecast precisions of $1.9{{\ \rm per\ cent}}$, $2.0{{\ \rm per\ cent}}$, and $1.9{{\ \rm per\ cent}}$ on Ωm, σ8, and $S_8 \equiv \sigma _8\Omega _m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as Alens. Adding scales $r_p=0.3\!-\!3.0\, h^{-1} \, \mathrm{Mpc}$ improves the S8 precision by a factor of ∼1.6 relative to a large scale ($3.0\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$) analysis, equivalent to increasing the survey area by a factor of ∼2.6. Sharpening the Alens prior to $1{{\ \rm per\ cent}}$ further improves the S8 precision to $1.1{{\ \rm per\ cent}}$, and it amplifies the gain from including non-linear scales. Our emulator achieves per cent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.

The outermost stellar halo of NGC 5128 (Centaurus A): Radial structure

Context.The extended stellar halos of galaxies contain important clues for investigating their assembly history and evolution.Aims.We investigate the resolved stellar content and the extended halo of NGC 5128 as a function of galactocentric distance, and trace the halo outward to its currently detectable limits.Methods.We usedHubbleSpace Telescope images obtained with the WFPC2, ACS, and WFC3 cameras equipped withF606WandF814Wfilters to resolve individual red giant branch (RGB) stars in 28 independent pointings across the halo of NGC 5128. The stellar halo analysis for 14 of these pointings is presented here for the first time. Star counts from deepVIcolor-magnitude diagrams reaching at least 1.5 mag below the tip of the RGB are used to derive the surface density distribution of the halo. The contamination by Milky Way stars is assessed with a new control field, with models, and by combining optical and near-IR photometry.Results.We present a new calibration of the WFC3F606W + F814Wphotometry to the ground-basedVIphotometric system. The photometry shows that the stellar halo of NGC 5128 is dominated by old RGB stars that are present in all fields. TheV-band surface brightness of fields changes from 23 to 32 mag arcsec−2between the innermost field only 8.3 kpc from the galaxy center to our outermost halo fields, which are located 140 kpc away from the center along the major axis and 92 kpc along the minor axis. Within the inner ∼30 kpc, we also find evidence for a 2 − 3 Gyr old population traced by asymptotic giant branch stars that are brighter than the tip of the RGB. This population contributes only up to 10% in total stellar mass if it is 2 Gyr old, but a larger fraction of 30 − 40% is required if its age is 3 Gyr. The stellar surface density profile is well fit by a classicr1/4curve or a simple power-law form ∼r−3.1over the full radial range, with no obvious break in the slope, but with large field-to-field scatter. The ellipticity measured from integrated-light photometry in the inner parts,e = (b/a) = 0.77, flattens toe = 0.54 ± 0.02 beyond 30 kpc. Considering the flattening of the outer halo, the projection of the elliptical isophote on the semimajor axis for our most distant field reaches nearly 30 effective radii.

Radial profiles of lensed z ∼ 1 galaxies on sub-kiloparsec scales

We study the spatially resolved physical properties of the Cosmic Snake arc in MACS J1206.2–0847 and the arc in Abell 0521 (A521). These are two strongly lensed galaxies at redshifts z = 1.036 and z = 1.044. We used observations of the Hubble Space Telescope (HST) and the Atacama Large Millimeter/submillimeter Array (ALMA). The former gives access to the star formation rate (SFR) and stellar mass (M⋆), and the latter to the H2 molecular gas mass (Mmol). HST and ALMA observations have similar angular resolutions of 0.15″ − 0.2″, which with the help of strong gravitational lensing enable us to reach spatial resolutions down to ∼30 pc and ∼100 pc in these two galaxies, respectively. These resolutions are close to the resolution of observations of nearby galaxies. We study the radial profiles of SFR, M⋆, and Mmol surface densities of these high-redshift galaxies and compare the corresponding exponential scale lengths with those of local galaxies. We find that the scale lengths in the Cosmic Snake are about 0.5 kpc − 1.5 kpc, and they are 3–10 times larger in A521. This is a significant difference knowing that the two galaxies have comparable integrated properties. These high-redshift scale lengths are nevertheless comparable to those of local galaxies, which cover a wide distribution. The particularity of our high-redshift radial profiles is the normalisation of the Mmol surface density profiles (ΣMmol), which are offset by up to a factor of 20 with respect to the profiles of z = 0 counterparts. The SFR surface density profiles are also offset by the same factor as ΣMmol, as expected from the Kennicutt-Schmidt law.

The terrestrial planet formation around M dwarfs: insitu, inward migration, or reversed migration

ABSTRACT Terrestrial planets are commonly observed to orbit M dwarfs with close-in trajectories. In this work, we extensively perform N-body simulations of planetesimal accretion with three models of in situ, inward migration, and reversed migration to explore terrestrial formation in tightly compact systems of M dwarfs. In the simulations, the solid discs are assumed to be 0.01 per cent of the masses of host stars and spread from 0.01 to 0.5 au with the surface density profile scaling with r−k according to the observations. Our results show that the in-situ scenario may produce $7.77^{+3.23}_{-3.77}$ terrestrial planets with an average mass of $1.23^{+4.01}_{-0.93}\,{\rm M}_{\oplus }$ around M dwarfs. The number of planets tends to increase as the disc slope is steeper or with a larger stellar mass. Moreover, we show that $2.55^{+1.45}_{-1.55}$ planets with a mass of $3.76^{+8.77}_{-3.46} \ {\rm M}_{\oplus }$ are formed in the systems via inward migration, while $2.85^{+1.15}_{-0.85}$ planets with $3.01^{+13.77}_{-2.71} \ {\rm M}_{\oplus }$ are yielded under reversed migration. Migration scenarios can also deliver plentiful water from the exterior of the ice line to the interior due to more efficient accretion. The simulation outcomes of the reversed migration model produce the best match with observations, being suggestive of a likely mechanism for planetary formation around M dwarfs.

Constraining the properties of the potential embedded planets in the disk around HD 100546

Context. The protoplanetary disk around the star HD 100546 displays prominent substructures in the form of two concentric rings. Recent observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) have revealed these features with high angular resolution and have resolved the faint outer ring well. This allows us to study the nature of the system further. Aims. Our aim is to constrain some of the properties of potential planets embedded in the disk, assuming that they induce the observed rings and gaps. Methods. We present the self-calibrated 0.9 mm ALMA observations of the dust continuum emission from the circumstellar disk around HD 100546. These observations reveal substructures in the disk that are consistent with two rings, the outer ring being much fainter than the inner one. We reproduced this appearance closely with a numerical model that assumes two embedded planets. We varied planet and disk parameters in the framework of the planet-disk interaction code FARGO3D and used the outputs for the gas and dust distribution to generate synthetic observations with the code RADMC-3D. Results. From this comparison, we find that an inner planet located at r1 = 13 au with a mass M1 = 8 MJup and an outer planet located at r2 = 143 au with a mass M2 = 3 MJup leads to the best agreement between synthetic and ALMA observations (deviation less than 3σ for the normalized radial profiles). To match the very low brightness of the outer structure relative to the inner ring, the initial disk gas surface density profile needs to follow an exponentially tapered power law (self-similar solution), rather than a simple power-law profile.

Mercury as the Relic of Earth and Venus Outward Migration

Abstract In spite of substantial advancements in simulating planet formation, the planet Mercury’s diminutive mass and isolated orbit and the absence of planets with shorter orbital periods in the solar system continue to befuddle numerical accretion models. Recent studies have shown that if massive embryos (or even giant planet cores) formed early in the innermost parts of the Sun’s gaseous disk, they would have migrated outward. This migration may have reshaped the surface density profile of terrestrial planet-forming material and generated conditions favorable to the formation of Mercury-like planets. Here we continue to develop this model with an updated suite of numerical simulations. We favor a scenario where Earth’s and Venus’s progenitor nuclei form closer to the Sun and subsequently sculpt the Mercury-forming region by migrating toward their modern orbits. This rapid formation of ∼0.5 M ⊕ cores at ∼0.1–0.5 au is consistent with modern high-resolution simulations of planetesimal accretion. In successful realizations, Earth and Venus accrete mostly dry, enstatite chondrite–like material as they migrate, thus providing a simple explanation for the masses of all four terrestrial planets, the inferred isotopic differences between Earth and Mars, and Mercury’s isolated orbit. Furthermore, our models predict that Venus’s composition should be similar to the Earth’s and possibly derived from a larger fraction of dry material. Conversely, Mercury analogs in our simulations attain a range of final compositions.