Astrophysics

The Baryonic Tully-Fisher Relation (BTFR) links baryonic mass of rotationally supported galaxies to their flat disk velocities. A popular form of the BTFR linked to MOND is based on an empirically determined characteristic acceleration, a0 that serves as the constant of proportionality. In this work, we propose an alternative, parametric form of the BTFR employing individual galactic properties from the SPARC galaxy database. Based on this data we find that a precise mass-velocity correlation is possible taking into account galactic disk radius and two dynamical related properties; mass discrepancy and disk surface density. We find no need to invoke a characteristic acceleration constant although its ansatz can be extracted and compared to several recent analyses also arguing against the MOND interpretation of a0. This improved BTFR finding has ramifications for the Radial Acceleration Relation (RAR) as well. Rather than a universal relation that describes the dynamics of all rotationally-supported galaxies, we find the RAR consists of a statistically distributed family of curves, reflecting the unique properties attributed to individual galaxies.

We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass the footprints of observations associated with the five CANDELS fields. In order to allow field-to-field variance to be explored, we have created eight realizations of each field. In this paper, we present comparisons with observable global galaxy properties, including counts in observed frame bands, luminosity functions, color-magnitude distributions and color-color distributions. We additionally present comparisons with physical galaxy parameters derived from SED fitting for the CANDELS observations, such as stellar masses and star formation rates. We find relatively good agreement between the model predictions and CANDELS observations for luminosity and stellar mass functions. We find poorer agreement for colors and star formation rate distributions. All of the mock lightcones as well as curated "theory friendly" versions of the observational CANDELS catalogs are made available through a web-based data hub.

We have studied the properties of molecular clouds in the second quadrant of the Milky Way Mid-plane from l = 104 . ??75 to l = 119 . ??75 and b =??5 . ??25 to b = 5 . ??25 using the 12 CO, 13 CO, and C 18 O J=1?? emission line data from the Milky Way Imaging Scroll Painting project (MWISP). We have identified 857 and 300 clouds in the 12 CO and 13 CO spectral cubes, respectively, using the DENDROGRAM + SCIMES algorithms. The distances of the molecular clouds are estimated and the physical properties like masses, sizes, and surface densities of the clouds are tabulated. The molecular clouds in the Perseus arm are about 30 ??50 times more massive and 4 ??6 times larger than the clouds in the Local arm. This result, however, is likely biased by distance selection effects. The surface densities of the clouds are enhanced in the Perseus arm with an average value of ??100 M ??pc ?? . We selected the 40 most extended ( > 0.35 arcdeg 2 ) molecular clouds from the 12 CO catalog to build the H 2 column density probability distribution function (N-PDF). About 78\% of the N-PDFs of the selected molecular clouds are well fitted with log-normal functions with only small deviations at high-densities which correspond to star-forming regions with scales of ??1-5 pc in the Local arm and ??5-10 pc in the Perseus arm. About 18\% of the selected molecular clouds have power-law N-PDFs at high-densities. In these molecular clouds, the majority of the regions fitted with the power-law correspond to molecular clumps of sizes of ??1 pc or filaments of widths of ??1 pc.

We use the ~370 square degrees data from the MWISP CO survey to study the vertical distribution of the molecular clouds (MCs) toward the tangent points in the region of l=[16,52]deg and |b|<5.1deg. The molecular disk consists of two components with the layer thickness (FWHM) of ~85pc and ~280pc, respectively. In the inner Galaxy, the molecular mass in the thin disk is dominant, while the molecular mass traced by the discrete MCs with weak CO emission in the thick disk is probably <10% of the whole molecular disk. For the CO gas in the thick disk, we identified 1055 high-z MCs that are >100pc from the Galactic plane. However, only a few samples (i.e., 32 MCs) are located in the |z|>360pc region. Typically, the discrete MCs of the thick disk population have a median peak temperature of 2.1 K, a median velocity dispersion of 0.8km/s, and a median effective radius of 2.5pc. The median surface density of these MCs is 6.8 Msun/pc^2, indicating very faint CO emission for these high-z MCs. The cloud-cloud velocity dispersion is 4.9+-1.3 km/s and a linear variation with a slope of -0.4 km/s/kpc is obtained in the region of R_GC=2.2-6.4kpc. Assuming that these clouds are supported by their turbulent motions against the gravitational pull of the disk, a model of rho0(R) = 1.28exp(-R/3.2kpc) Msun/pc^3 can be used to describe the distribution of the total mass density in the Galactic midplane.

I simulate the evolution of a stellar system hosting two stellar populations whose initial set up is defined according to the two main scenarios proposed for the origin of multiple populations in Galactic globular clusters: (i) formation of a second generation from a cooling flow of pristine+polluted gas and (ii) accretion of polluted gas onto the proto-stellar disks of a fraction of low-mass stars. For this purpose, Monte Carlo simulations containing from 10 5 up to 3??10 6 particles have been run including the effect of stellar evolution, binary interactions, external tidal field and a detailed modelling of the proto-stellar disk structure. The early accretion of gas onto proto-stellar disks is unable to produce discrete populations and to alter the chemical composition of a significant ( >10% ) fraction of stars unless a disk lifetime larger ( t disk ??0 Myr ) than that predicted by models is assumed. Moreover, in this scenario the mixing timescale of the two populations is too short to reproduce the observed segregation of the chemically enriched population. On the other hand, simulations run within the cooling flow scenario can evolve after a Hubble time into stellar systems with a first-to-second population mass ratio similar to that observed in globular clusters, provided that an initial filling-factor r h / r J >0.15 is adopted. However, in the weak tidal field regime a radial segregation of the second population stronger than what observed in Milky Way globular clusters at large Galactocentric distances is predicted. This discrepancy disappears in simulations following eccentric orbits in a realistic axisymmetric potential.

We present a new study of archival ALMA observations of the CO(2-1) line emission of the host galaxy of quasar RX J1131 at redshift z =0.654, lensed by a foreground galaxy. A simple lens model is shown to well reproduce the optical images obtained by the Hubble Space Telescope. Clear evidence for rotation of the gas contained in the galaxy is obtained and a simple rotating disc model is shown to give an excellent overall description of the morpho-kinematics of the source. The possible presence of a companion galaxy suggested by some previous authors is not confirmed. Detailed comparison between model and observations gives evidence for a more complex dynamics than implied by the model. Doppler velocity dispersion within the beam size in the image plane is found to account for the observed line width.

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively.

We present a new suite of mock galaxy catalogs mimicking the low-redshift Universe, based on an updated halo occupation distribution (HOD) model and a scaling relation between optical properties and the neutral hydrogen (HI) content of galaxies. Our algorithm is constrained by observations of the luminosity function and luminosity- and colour-dependent clustering of SDSS galaxies, as well as the HI mass function and HI-dependent clustering of massive HI-selected galaxies in the ALFALFA survey. Mock central and satellite galaxies with realistic values of r -band luminosity, g?�r and u?�r colour, stellar mass and HI mass are populated in an N -body simulation, inheriting a number of properties of the density and tidal environment of their host halos. The host halo of each central galaxy is also `baryonified' with realistic spatial distributions of stars as well as hot and cold gas, along with the corresponding rotation curve. Our default HOD assumes that galaxy properties are a function of group halo mass alone, and can optionally include effects such as galactic conformity and colour-dependent galaxy assembly bias. The mocks predict the relation between the stellar mass and HI mass of massive HI galaxies, as well as the 2-point cross-correlation function of spatially co-located optical and HI-selected samples. They enable novel null tests for galaxy assembly bias, provide predictions for the HI velocity width function, and clarify the origin and universality of the radial acceleration relation in the ? CDM framework.

Light scattering at near-infrared wavelengths has been used to study the optical properties of the interstellar dust grains, but these studies are limited by the assumptions on the strength of the radiation field. On the other hand, thermal dust emission can be used to constrain the properties of the radiation field, although this is hampered by uncertainty about the dust emissivity. We test if current dust models allow us to model a molecular cloud simultaneously in the near infrared (NIR) and far infrared (FIR) wavelengths and compare the results with observations. Our aim is to place constraints on the properties of the dust grains and the strength of the radiation field. We present computations of dust emission and scattered light of a quiescent molecular cloud LDN1512. We construct radiative transfer models for LDN1512 that include an anisotropic radiation field and a three-dimensional cloud model. We are able to reproduce the observed FIR observations, with a radiation field derived from the DIRBE observations, with all of the tested dust models. However, with the same density distribution and the assumed radiation field, the models fail to reproduce the observed NIR scattering in all cases except for models that take into account dust evolution via coagulation and mantle formation. We find that the column densities derived from our radiative transfer modelling can differ by a factor of up to two, compared to the column densities derived from the observations with modified blackbody fits. The discrepancy in the column densities is likely caused because of temperature difference between a modified blackbody fit and the real spectra. We show that the observed dust emission can be reproduced with several different assumptions about the properties of the dust grains. However, in order to reproduce the observed scattered surface brightness dust evolution must be taken into account.

The mechanism of thermal driving for launching mass outflows is interconnected with classical thermal instability (TI). In a recent paper, we demonstrated that as a result of this interconnectedness, radial wind solutions of X-ray heated flows are prone to becoming clumpy. In this paper, we first show that the Bernoulli function determines whether or not the entropy mode can grow due to TI in dynamical flows. Based on this finding, we identify a critical `unbound' radius beyond which TI should accompany thermal driving. Our numerical disk wind simulations support this result and reveal that clumpiness is a consequence of buoyancy disrupting the stratified structure of steady state solutions. Namely, instead of a smooth transition layer separating the highly ionized disk wind from the cold phase atmosphere below, hot bubbles formed from TI rise up and fragment the atmosphere. These bubbles first appear within large scale vortices that form below the transition layer, and they result in the episodic production of distinctive cold phase structures referred to as irradiated atmospheric fragments (IAFs). Upon interacting with the wind, IAFs advect outward and develop extended crests. The subsequent disintegration of the IAFs takes place within a turbulent wake that reaches high elevations above the disk. We show that this dynamics has the following observational implications: dips in the absorption measure distribution are no longer expected within TI zones and there can be a less sudden desaturation of X-ray absorption lines such as \OVIII as well as multiple absorption troughs in \FeXXVK.