To help understand our universe better, researchers and scientists currently run extreme-scale cosmology simulations on leadership supercomputers. However, such simulations can generate large amounts of scientific data, which often result in expensive costs in data associated with data movement and storage. Lossy compression techniques have become attractive because they significantly reduce data size and can maintain high data fidelity for post-analysis. In this paper, we propose to use GPU-based lossy compression for extreme-scale cosmological simulations. Our contributions are threefold: (1) we implement multiple GPU-based lossy compressors to our opensource compression benchmark and analysis framework named Foresight; (2) we use Foresight to comprehensively evaluate the practicality of using GPU-based lossy compression on two real-world extreme-scale cosmology simulations, namely HACC and Nyx, based on a series of assessment metrics; and (3) we develop a general optimization guideline on how to determine the best-fit configurations for different lossy compressors and cosmological simulations. Experiments show that GPU-based lossy compression can provide necessary accuracy on post-analysis for cosmological simulations and high compression ratio of 5~15x on the tested datasets, as well as much higher compression and decompression throughput than CPU-based compressors.

We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged]

We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^{13}Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The inwards deceleration of the filament and the prevalence of blueshifted components in the Lya spectra support a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also recombinations still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside a massive halo at high redshift.

The third-generation ground-based gravitational-wave (GW) detector, Cosmic Explorer (CE), is scheduled to start its observation in the 2030s. In this paper, we make a forecast for cosmological parameter estimation with gravitational-wave standard siren observation from the CE. We use the simulated GW standard siren data of CE to constrain the $\Lambda$CDM, $w$CDM and CPL models. We combine the simulated GW data with the current cosmological electromagnetic observations including the latest cosmic microwave background anisotropies data from Planck, the optical baryon acoustic oscillation measurements, and the type Ia supernovae observation (Pantheon compilation) to do the analysis. We find that the future standard siren observation from CE will improve the cosmological parameter estimation to a great extent, since the future GW standard siren data can well break the degeneracies generated by the optical observations between various cosmological parameters. We also find that the CE's constraining capability on the cosmological parameters is slightly better than that of the same-type GW detector, the Einstein Telescope. In addition, the synergy between the GW standard siren observation from CE and the 21 cm emission observation from SKA is also discussed.

We have obtained spectroscopic confirmation with Hubble Space Telescope WFC3/G141 of a first sizeable sample of nine quiescent galaxies at 2.4<z<3.3. Their average near-UV/optical rest-frame spectrum is characterized by low attenuation (Av$\sim$0.6 mag) and a strong Balmer break, larger than the 4000 A break, corresponding to a fairly young age of $\sim$300 Myr. This formally classifies a substantial fraction of classically selected quiescent galaxies at z $\sim$ 3 as post-starbursts, marking their convergence to the quenching epoch. The rapid spectral evolution with respect to z $\sim$ 1.5 quiescent galaxies is not matched by an increase of residual star-formation, as judged from the weak detection of [O II]3727A emission, pointing to a flattening of the steep increase in gas fractions previously seen from z $\sim$ 0 to 1.8. However, radio 3GHz stacked emission implies either much stronger dust-obscured star formation or substantial further evolution in radio-mode AGN activity with respect to z $\sim$ 1.5.

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

We report Atacama Large Millimetre Array (ALMA) observations of four high-redshift dusty star-forming galaxy candidates selected from far-Infrared (FIR)/submm observations in the COSMOS field. We securely detect all galaxies in the continuum and spectroscopically confirm them at z=3.62--5.85 using ALMA 3mm line scans, detecting multiple CO and/or [CI] transitions. This includes the most distant dusty galaxy currently known in the COSMOS field, ID85001929 at z=5.847. These redshifts are lower than we had expected as these galaxies have substantially colder dust temperatures (i.e., their SEDs peak at longer rest frame wavelengths) than most literature sources at z>4. The observed cold dust temperatures are best understood as evidence for optically thick dust continuum in the FIR, rather than the result of low star formation efficiency with rapid metal enrichment. We provide direct evidence that, given their cold spectral energy distributions, CMB plays a significant role biasing their observed Rayleigh-Jeans (RJ) slopes to unlikely steep values and, possibly, reducing their CO fluxes by a factor of two. We recover standard RJ slopes when the CMB contribution is taken into account. High resolution ALMA imaging shows compact morphology and evidence for mergers. This work reveals a population of cold dusty star-forming galaxies that were under-represented in current surveys, and are even colder than typical Main Sequence galaxies at the same redshift. High FIR dust optical depth might be a widespread feature of compact starbursts at any redshift.

The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. We produce spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assess our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise levels. Then, we analyze such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ~0.1 Gyr for stellar populations older than ~1.5 Gyr, pushing up to an age of ~1 Gyr for stellar populations older than ~5 Gyr. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.

While the formation of stellar clumps in distant galaxies is usually attributed to gravitational violent disk instabilities, we show here that major mergers also represent a competitive mechanism to form bright clumps. Using ~0.1'' resolution ACS F814W images in the entire COSMOS field, we measure the fraction of clumpy emission in 109 main sequence (MS) and 79 Herschel-detected starbursts (off-MS) galaxies at 0.5 < z < 0.9, representative of normal versus merger induced star-forming activity, respectively. We additionally identify merger samples from visual inspection and from Gini-M20 morphological parameters. Regardless of the merger criteria adopted, the clumpiness distribution of merging systems is different from that of normal isolated disks at > 99.5 % confidence level, with the former reaching higher clumpiness values, up to 20 % of the total galaxy emission. We confirm the merger induced clumpiness enhancement with novel hydrodynamical simulations of colliding galaxies with gas fractions typical of z ~ 0.7. Multi-wavelength images of three starbursts in the CANDELS field support the young nature of clumps, which are likely merger products rather than older pre-existing structures. Finally, for a subset of 19 starbursts with existing near-IR rest frame spectroscopy, we find that the clumpiness is mildly anti-correlated with the merger phase, decreasing towards final coalescence. Our result can explain recent ALMA detections of clumps in hyperluminous high-z starbursts, while normal objects are smooth. This work raises a question on the role of mergers on the origin of clumps in high redshift galaxies in general.

In this work, we use the simulated gravitational wave (GW) standard siren data from the future observation of the Einstein Telescope (ET) to constrain various dark energy cosmological models, including the $\Lambda$CDM, $w$CDM, CPL, $\alpha$DE, GCG, and NGCG models. We also use the current mainstream cosmological electromagnetic observations, i.e., the cosmic microwave background anisotropies data, the baryon acoustic oscillations data, and the type Ia supernovae data, to constrain these models. We find that the GW standard siren data could tremendously improve the constraints on the cosmological parameters for all these dark energy models. For all the cases, the GW standard siren data can be used to break the parameter degeneracies generated by the current cosmological electromagnetic observational data. Therefore, it is expected that the future GW standard siren observation from the ET would play a crucial role in the cosmological parameter estimation in the future. The conclusion of this work is quite solid because it is based on the analysis for various dark energy models.

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4<z<6.0 in the COSMOS field with robust spectroscopic redshifts. We achieve high photo-z precision ($\sigma_{\Delta z/(1+z_{spec})}\lesssim0.04$), high accuracy (i.e., minimal offset biases; median$(\Delta z/(1+z_{spec}))\lesssim0.02$), and low catastrophic failure rates ($\eta\simeq4\%$) over all redshifts. Interestingly, we find that a weak 2175A absorption feature in the attenuation curve models is required to remove a subtle systematic photo-z offset ($z_{phot}-z_{spec}\simeq-0.03$) that occurs when this feature is not included. As expected, the accuracy of derived physical properties in MAGPHYS+photo-z decreases strongly as redshift uncertainty increases. The all-in-one treatment of uncertainties afforded with this code is beneficial for accurately interpreting physical properties of galaxies in large photometric datasets. Finally, we emphasize that MAGPHYS+photo-z is not intended to replace existing photo-z codes, but rather offer flexibility to robustly interpret physical properties when spectroscopic redshifts are unavailable. The MAGPHYS+photo-z code is publicly available online.

We derive the geometric albedo of a near-Earth asteroid, (4179) Toutatis, to investigate its surface physical conditions. The asteroid has been studied rigorously not only via ground-based photometric, spectrometric, polarimetric, and radar observations but also via \textit{in situ} observation by the Chinese Chang'e-2 space probe; however, its geometric albedo is not well understood. We conducted V-band photometric observations when the asteroid was at opposition in April 2018 using the three telescopes in the southern hemisphere that compose the Korea Microlensing Telescope Network (KMTNet). The observed time-variable cross section was corrected using the radar shape model. We find that Toutatis has a geometric albedo $p_\mathrm{V} = 0.185^{+0.045}_{-0.039} $, which is typical of S-type asteroids. We compare the geometric albedo with archival polarimetric data and further find that the polarimetric slope--albedo law provides a reliable estimate for the albedo of this S-type asteroid. The thermal infrared observation also produced similar results if the size of the asteroid is updated to match the results from Chang'e-2. We conjecture that the surface of Toutatis is covered with grains smaller than that of the near-Sun asteroids including (1566) Icarus and (3200) Phaethon.

We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk surface density profile follows a power-law of the form $\rho_r \propto r^{-4}$. For the disk density profiles that can lead to a small reduced $\chi^2$ of goodness-of-fit, we find that there is a clear linear correlation between the disk mass and the power-law index $\gamma$ in the equation of disk density profile. This suggests that the $^{12}$CO disk of LkCa 15 is optically thick and we can fit its $^{12}$CO radial intensity profile using either a lower disk mass with a smaller $\gamma$ or a higher disk mass with a bigger $\gamma$. By comparing the $^{12}$CO channel maps of the best-fit model with disk models with higher or lower masses, we find that a disk mass of $\sim$ 0.1 $M_{\odot}$ can best reproduce the observed morphology of the $^{12}$CO channel maps. The dust continuum map at 0.87 mm of the LkCa 15 disk shows an inner cavity of the similar size of the best-fit gas model, but its out edge is at $\sim$ 200 AU, much smaller than the fitted gas disk. Such a discrepancy between the outer edges of the gas and dust disks is consistent with dust drifting and trapping models.

ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$.

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr\`o et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase.

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs.

We present the first results of an ALMA survey of the lower fine structure line of atomic carbon [C I]$(^3P_1\,-\,^{3}P_0)$ in far infrared-selected galaxies on the main sequence at $z\sim1.2$ in the COSMOS field. We compare our sample with a comprehensive compilation of data available in the literature for local and high-redshift starbursting systems and quasars. We show that the [C I]($^3P_1$$\rightarrow$$^3P_0$) luminosity correlates on global scales with the infrared luminosity $L_{\rm IR}$ similarly to low-$J$ CO transitions. We report a systematic variation of $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L_{\rm IR}$ as a function of the galaxy type, with the ratio being larger for main-sequence galaxies than for starbursts and sub-millimeter galaxies at fixed $L_{\rm IR}$. The $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$/$L'_{\rm CO(2-1)}$ and $M_{\rm{[C I]}}$/$M_{\rm dust}$ mass ratios are similar for main-sequence galaxies and for local and high-redshift starbursts within a 0.2 dex intrinsic scatter, suggesting that [C I] is a good tracer of molecular gas mass as CO and dust. We derive a fraction of $f_{\rm{[C\,I]}} = M_{\rm{[C\,I]}} / M_{\rm{C}}\sim3-13$% of the total carbon mass in the atomic neutral phase. Moreover, we estimate the neutral atomic carbon abundance, the fundamental ingredient to calibrate [C I] as a gas tracer, by comparing $L'_{\rm [C\,I]^3P_1\,-\, ^3P_0}$ and available gas masses from CO lines and dust emission. We find lower [C I] abundances in main-sequence galaxies than in starbursting systems and sub-millimeter galaxies, as a consequence of the canonical $\alpha_{\rm CO}$ and gas-to-dust conversion factors. This argues against the application to different galaxy populations of a universal standard [C I] abundance derived from highly biased samples.

We present an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed.

We present a "super-deblended" far-infrared to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al. 2018, with key adaptations. We obtain point spread function (PSF) fitting photometry at fixed prior positions including 88,008 galaxies detected in either VLA 1.4~GHz, 3~GHz and/or MIPS 24~$\mu$m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed ``active' removal of non relevant priors at FIR/(sub)mm bands using spectral energy distribution (SED) fitting and redshift information. In order to cope with the shallower COSMOS data we subtract from the maps the flux of faint non-fitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer, Herschel, SCUBA2, AzTEC, MAMBO and NSF's Karl G. Jansky Very Large Array at 3~GHz and 1.4~GHz) displays well behaved quasi-Gaussian uncertainties, calibrated from Monte Carlo simulations and tailored to observables (crowding, residual maps). Comparison to ALMA photometry for hundreds of sources provide a remarkable validation of the technique. We detect 11,220 galaxies over the 100--1200~$\mu$m range, extending to $z_{\rm phot}\sim7$. We conservatively selected a sample of 85 $z>4$ high redshift candidates, significantly detected in the FIR/(sub)mm, often with secure radio and/or Spitzer/IRAC counterparts. This provides a chance to investigate the first generation of vigorous starburst galaxies (SFRs$\sim1000\mathrm{M}_\odot$~yr$^{-1}$). The photometric and value added catalogs are publicly released.

We obtained optical/near-IR rest-frame Magellan FIRE spectra (including Pa$\beta$ and Pa$\gamma$) of 25 starburst galaxies at 0.5<z<0.9, with average star formation rates (SFR) x7 above the Main Sequence (MS). We find that Paschen-to-Balmer line ratios saturate around a constant value corresponding to $A_{\rm V}\sim$2-3 mag, while line to IR luminosity ratios suggest a large range of more extreme obscurations and appear to be uncorrelated to the former. This behavior is not consistent with standard attenuation laws derived for local and distant galaxies, while being remarkably consistent with observations of starburst cores in which young stars and dust are homogeneously mixed. This model implies $A_{\rm V}=$2-30 mag attenuation to the center of starburst cores, with a median of ~9 mag (a factor of 4000). X-ray hardness ratios for 6 AGNs in our sample and column densities derived from observed dust masses and radio sizes independently confirm this level of attenuation. In these conditions observed optical/near-IR emission comes from surface regions, while inner starburst cores are invisible. We thus attribute the high [NII]/H$\alpha$ ratios to widespread shocks from accretion, turbulence and dynamic disturbances rather than to AGNs. The large range of optical depths demonstrates that substantial diversity is present within the starburst population, possibly connected to different merger phases or progenitor properties. The majority of our targets are, in fact, morphologically classified as mergers. We argue that the extreme obscuration provides in itself smoking gun evidence of their merger origin, and a powerful tool for identifying mergers at even higher redshifts.

Recent Atacama Large Millimeter and Submillimeter Array (ALMA) observations of the protoplanetary disk around the Herbig Ae star HD 163296 revealed three depleted dust gaps at 60, 100 and 160 au in the 1.3 mm continuum as well as CO depletion in the middle and outer dust gaps. However, no CO depletion was found in the inner dust gap. To examine the planet--disk interaction model, we present results of two-dimensional two fluid (gas + dust) hydrodynamic simulations coupled with three-dimensional radiative transfer simulations. In order to fit the high gas-to-dust ratio of the first gap, we find the Shakura--Sunyaev viscosity parameter $\alpha$ must be very small ($\lesssim 10^{-4}$) in the inner disk. On the other hand, a relatively large $\alpha$ ($\sim 7.5\times 10^{-3}$) is required to reproduce the dust surface density in the outer disk. We interpret the variation of $\alpha$ as an indicator of the transition from an inner dead zone to the outer magnetorotational instability (MRI) active zone. Within $\sim 100$ au, the HD 163296 disk's ionization level is low, and non-ideal magnetohydrodynamic (MHD) effects could suppress the MRI, so the disk can be largely laminar. The disk's ionization level gradually increases toward larger radii, and the outermost disk ($r > 300$ au) becomes turbulent due to MRI. Under this condition, we find that the observed dust continuum and CO gas line emissions can be reasonably fit by three half-Jovian-mass planets (0.46, 0.46 and 0.58 $M_\textrm{J}$) at 59, 105 and 160 au, respectively.

We present a new technique to measure multi-wavelength "Super-deblended" photometry from highly confused images, which we apply to Herschel and ground-based far-infrared (FIR) and (sub-)millimeter (mm) data in the northern field of the Great Observatories Origins Deep Survey (GOODS). There are two key novelties. First, starting with a large database of deep Spitzer 24{\mu}m and VLA 20cm detections that are used to define prior positions for fitting the FIR/submm data, we perform an active selection of useful priors independently at each frequency band, moving from less to more confused bands. Exploiting knowledge of redshift and all available photometry, we identify hopelessly faint priors that we remove from the fitting pool. This approach significantly reduces blending degeneracies and allows reliable photometry to be obtained for galaxies in FIR+mm bands. Second, we obtain well-behaved, nearly Gaussian flux density uncertainties, individually tailored to all fitted priors in each band. This is done by exploiting extensive simulations that allow us to calibrate the conversion of formal fitting uncertainties to realistic uncertainties depending on quantities directly measurable. We achieve deeper detection limits with high fidelity measurements and uncertainties at FIR+mm bands. As an illustration of the utility of these measurements, we identify 70 galaxies with z>3 and reliable FIR+mm detections. We present new constraints on the cosmic star formation rate density at 3<z<6, finding a significant contribution from z>3 dusty galaxies that are missed by optical-to-near-infrared color selection. Photometric measurements for 3306 priors, including over 1000 FIR+mm detections are released publicly with our catalog.

Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time.

We use an end-to-end model of planet formation, thermodynamic evolution, and atmospheric escape to investigate how the statistical imprints of evaporation depend on the bulk composition of planetary cores (rocky vs. icy). We find that the population-wide imprints like the location of the "evaporation valley" in the distance-radius plane and the corresponding bimodal radius distribution clearly differ depending on the bulk composition of the cores. Comparison with the observed position of the valley (Fulton et al. 2017) suggests that close-in low-mass Kepler planets have a predominately Earth-like rocky composition. Combined with the excess of period ratios outside of MMR, this suggests that low-mass Kepler planets formed inside of the water iceline, but still undergoing orbital migration. The core radius becomes visible for planets losing all primordial H/He. For planets in this "triangle of evaporation" in the distance-radius plane, the degeneracy in compositions is reduced. In the observed diagram, we identify a trend to more volatile-rich compositions with increasing radius (R/R_Earth<1.6 rocky; 1.6-3.0 ices and/or H/He; >3: H/He). The mass-density diagram contains important information about formation and evolution. Its characteristic broken V-shape reveals the transitions from solid planets to low-mass core-dominated planets with H/He and finally to gas-dominated giants. Evaporation causes density and orbital distance to be anti-correlated for low-mass planets, in contrast to giants, where closer-in planets are less dense, likely due to inflation. The temporal evolution of the statistical properties reported here will be of interest for the PLATO 2.0 mission which will observe the temporal dimension.

We present the discovery of a very faint stellar system, SMASH 1, that is potentially a satellite of the Large Magellanic Cloud. Found within the Survey of the MAgellanic Stellar History (SMASH), SMASH 1 is a compact ($r_h = 9.1^{+5.9}_{-3.4}$ pc) and very low luminosity (M_V = -1.0 +/- 0.9, $L_V=10^{2.3 +/- 0.4}$ Lsun) stellar system that is revealed by its sparsely populated main sequence and a handful of red-giant-branch candidate member stars. The photometric properties of these stars are compatible with a metal-poor ([Fe/H]=-2.2) and old (13 Gyr) isochrone located at a distance modulus of ~18.8, i.e. a distance of ~57 kpc. Situated at 11.3$^\circ$ from the LMC in projection, its 3-dimensional distance from the Cloud is ~13 kpc, consistent with a connection to the LMC, whose tidal radius is at least 16 kpc. Although the nature of SMASH 1 remains uncertain, its compactness favors it being a stellar cluster and hence dark-matter free. If this is the case, its dynamical tidal radius is only <19 pc at this distance from the LMC, and smaller than the system's extent on the sky. Its low luminosity and apparent high ellipticity ($\epsilon=0.62^{+0.17}_{-0.21}$) with its major axis pointing toward the LMC may well be the tell-tale sign of its imminent tidal demise.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We show that the most distant X-ray detected cluster known to date, ClJ1001 at z=2.506, hosts a strong overdensity of radio sources. Six of them are individually detected (within 10") in deep 0.75" resolution VLA 3GHz imaging, with S(3GHz)>8uJy. Of the six, AGN likely affects the radio emission in two galaxies while star formation is the dominant source powering the remaining four. We searched for cluster candidates over the full COSMOS 2-square degree field using radio-detected 3GHz sources and looking for peaks in Sigma5 density maps. ClJ1001 is the strongest overdensity by far with >10sigma, with a simple z_phot>1.5 preselection. A cruder photometric rejection of z<1 radio foregrounds leaves ClJ1001 as the second strongest overdensity, while even using all radio sources ClJ1001 remains among the four strongest projected overdensities. We conclude that there are great prospects for future, deep and wide-area radio surveys to discover large samples of the first generation of forming galaxy clusters. In these remarkable structures widespread star formation and AGN activity of massive galaxy cluster members, residing within the inner cluster core, will ultimately lead to radio continuum as one of the most effective means for their identification, with detection rates expected in the ballpark of 0.1-1 per square degree at z>2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models.

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event.

In this work, we extensively investigate the formation of near 4:2:1 mean motion resonances (MMRs) configuration by performing two sets of N-body simulations. We model the eccentricity damping, gas drag, type I and type II planetary migration of planetesimals, planetary embryos and giant planets in the first sets. For the simulations of giant planets with type II migration, the massive terrestrial planets, with a mass up to several Earth masses, are likely produced in the systems. We further show that by shepherding and/or scattering mechanisms through Jovian planet's type II migration, the terrestrial planets and giant planets in the systems can be evolved into a near 4:2:1 MMRs. Moreover, the models are applicable to the formation of Kepler-238 and 302 systems. In the second set, we study the 4:2:1 MMRs formation in the terrestrial planetary systems, where the planets undergo type I migration and eccentricity damping. By considering type I migration, $\sim$ 17.1\% of the simulations indicate that terrestrial planets are evolved into 4:2:1 MMRs. However, this probability should depend on the initial conditions of planets. Hence, we conclude that both type I and type II migration can play a crucial role in close-in terrestrial planet formation.

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages <1 Gyr is at least 3 times higher in the outer regions of the pBzK galaxies than in their cores. This results in a mild age gradient of ~<0.4 Gyr over ~6 kpc and suggests either the occurrence of widespread rejuvenation episodes or that inside-out quenching played a role in the passivization of this galaxy population. We also report on low-level star formation rates derived from the [OII]3727A emission line, with SFR_OII = 3.7-4.5 Msun/yr. This estimate is confirmed by an independent measurement on a separate sample of similarly-selected quiescent galaxies in the COSMOS field, using stacked far-infrared data (SFR_FIR = 2-4 Msun/yr). This second, photometric sample also displays significant excess at 1.4 GHz, suggestive of the presence of radio-mode AGN activity.

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after $10^4$ orbits, or $2.5$ Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

We report the discovery of one RR Lyrae star in the ultra--faint satellite galaxy Hydra II based on time series photometry in the g, r and i bands obtained with the Dark Energy Camera at Cerro Tololo Interamerican Observatory, Chile. The RR Lyrae star has a mean magnitude of $i = 21.30\pm 0.04$ which translates to a heliocentric distance of $151\pm 8$ kpc for Hydra II; this value is $\sim 13\%$ larger than the estimate from the discovery paper based on the average magnitude of several blue horizontal branch star candidates. The new distance implies a slightly larger half-light radius of $76^{+12}_{-10}$ pc and a brighter absolute magnitude of $M_V = -5.1 \pm 0.3$, which keeps this object within the realm of the dwarf galaxies. The pulsational properties of the RR Lyrae star ($P=0.645$ d, $\Delta g = 0.68$ mag) suggest Hydra II may be a member of the intermediate Oosterhoff or Oosterhoff II group. A comparison with other RR Lyrae stars in ultra--faint systems indicates similar pulsational properties among them, which are different to those found among halo field stars and those in the largest of the Milky Way satellites. We also report the discovery of 31 additional short period variables in the field of view (RR Lyrae, SX Phe, eclipsing binaries, and a likely anomalous cepheid). However, given their magnitudes and large angular separation from Hydra II, they must be field stars not related to Hydra II.

We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices of disk parameters, assuming an already settled dust distribution and no planet migration. We can obtain quite adequate fits to the observed dust emission using three planets with masses 0.35, 0.17, and 0.26 $M_{Jup}$ at 13.1, 33.0, and 68.6 AU, respectively. Implications for the planet formation as well as the limitations of this scenario are discussed.

We present the simplest nuclear energy density functional (NEDF) to date, determined by only 4 significant phenomenological parameters, yet capable of fitting measured nuclear masses with better accuracy than the Bethe-Weizs\"acker mass formula, while also describing density structures (charge radii, neutron skins etc.) and time-dependent phenomena (induced fission, giant resonances, low energy nuclear collisions, etc.). The 4 significant parameters are necessary to describe bulk nuclear properties (binding energies and charge radii); an additional 2 to 3 parameters have little influence on the bulk nuclear properties, but allow independent control of the density dependence of the symmetry energy and isovector excitations, in particular the Thomas-Reiche-Kuhn sum rule. This Hohenberg-Kohn-style of density functional theory successfully realizes Weizs\"acker's ideas and provides a computationally tractable model for a variety of static nuclear properties and dynamics, from finite nuclei to neutron stars, where it will also provide a new insight into the physics of the r-process, nucleosynthesis, and neutron star crust structure. This new NEDF clearly separates the bulk geometric properties - volume, surface, symmetry, and Coulomb energies which amount to 8MeV per nucleon or up to 2000MeV per nucleus for heavy nuclei - from finer details related to shell effects, pairing, isospin breaking, etc. which contribute at most a few MeV for the entire nucleus. Thus it provides a systematic framework for organizing various contributions to the NEDF. Measured and calculated physical observables - symmetry and saturation properties, the neutron matter equation of state, and the frequency of giant dipole resonances - lead directly to new terms not considered in current NEDF parameterizations.

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

We present a proper-motion study on models of the dwarf spheroidal galaxy Sculptor, based on the predicted proper-motion accuracy of Gaia measurements. Gaia will measure proper motions of several hundreds of stars for a Sculptor-like system. Even with an uncertainty on the proper motion of order 1.5 times the size of an individual proper-motion value of ~10 mas/century, we find that it is possible to recover Sculptor's systemic proper motion at its distance of 79 kpc.

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.

Despite the increase in observational data on exoplanets, the processes that lead to the formation of planets are still not well understood. But thanks to the high number of known exoplanets, it is now possible to look at them as a population that puts statistical constraints on theoretical models. A method that uses these constraints is planetary population synthesis. Its key element is a global model of planet formation and evolution that directly predicts observable planetary properties based on properties of the natal protoplanetary disk. To do so, global models build on many specialized models that address one specific physical process. We thoroughly review the physics of the sub-models included in global formation models. The sub-models can be classified as models describing the protoplanetary disk (gas and solids), the (proto)planet (solid core, gaseous envelope, and atmosphere), and finally the interactions (migration and N-body interaction). We compare the approaches in different global models and identify physical processes that require improved descriptions in future. We then address important results of population synthesis like the planetary mass function or the mass-radius relation. In these results, the global effects of physical mechanisms occurring during planet formation and evolution become apparent, and specialized models describing them can be put to the observational test. Due to their nature as meta models, global models depend on the development of the field of planet formation theory as a whole. Because there are important uncertainties in this theory, it is likely that global models will in future undergo significant modifications. Despite this, they can already now yield many testable predictions. With future global models addressing the geophysical characteristics, it should eventually become possible to make predictions about the habitability of planets.

We have analyzed the distributions in the color-magnitude diagram (CMD) of a large sample of face-on galaxies to minimize the effect of dust extinctions on galaxy color. About 300 thousand galaxies with $log(a/b) < $ 0.2 and redshift $z < 0.2$ are selected from the SDSS DR7 catalog. Two methods are employed to investigate the distributions of galaxies in the CMD including 1-D Gaussian fitting to the distributions in individual magnitude bins and 2-D Gaussian mixture model (GMM) fitting to galaxies as a whole. We find that in the 1-D fitting only two Gaussians are not enough to fit galaxies with the excess present between the blue cloud and the red sequence. The fitting to this excess defines the centre of the green-valley in the local universe to be $(u-r)_{0.1} = -0.121M_{r,0.1}-0.061$. The fraction of blue cloud and red sequence galaxies turns over around $M_{r,0.1} \sim -20.1$ mag, corresponding to stellar mass of $3\times10^{10}M_\odot$. For the 2-D GMM fitting, a total of four Gaussians are required, one for the blue cloud, one for the red sequence and the additional two for the green valley. The fact that two Gaussians are needed to describe the distributions of galaxies in the green valley is consistent with some models that argue for two different evolutionary paths from the blue cloud to the red sequence.

We present a new measurement of the $\alpha$-spectroscopic factor ($S_\alpha$) and the asymptotic normalization coefficient (ANC) for the 6.356 MeV 1/2$^+$ subthreshold state of $^{17}$O through the $^{13}$C($^{11}$B, $^{7}$Li)$^{17}$O transfer reaction and we determine the $\alpha$-width of this state. This is believed to have a strong effect on the rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction, the main neutron source for {\it slow} neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. Based on the new width we derive the astrophysical S-factor and the stellar rate of the $^{13}$C($\alpha$, $n$)$^{16}$O reaction. At a temperature of 100 MK our rate is roughly two times larger than that by \citet{cau88} and two times smaller than that recommended by the NACRE compilation. We use the new rate and different rates available in the literature as input in simulations of AGB stars to study their influence on the abundances of selected $s$-process elements and isotopic ratios. There are no changes in the final results using the different rates for the $^{13}$C($\alpha$, $n$)$^{16}$O reaction when the $^{13}$C burns completely in radiative conditions. When the $^{13}$C burns in convective conditions, as in stars of initial mass lower than $\sim$2 $M_\sun$ and in post-AGB stars, some changes are to be expected, e.g., of up to 25% for Pb in our models. These variations will have to be carefully analyzed when more accurate stellar mixing models and more precise observational constraints are available.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

We present a revised analysis of a speculated stellar bridge between the Milky Way dwarf galaxies Leo IV and Leo V. Using data acquired with Subaru/Suprime-Cam over a 1deg x 4deg field encompassing the two satellites and the region in between, we confirm our previous detection of a stellar overdensity between Leo IV and Leo V (de Jong et al. 2010). The larger area coverage and improved depth of our current dataset allow for an improved analysis of the stellar overdensity that had previously appeared to bridge the two galaxies. A main-sequence turn-off feature visible in the stacked colour-magnitude diagram of the contiguously observed Subaru fields reveals an extended stellar structure at a distance of approximately 20 kpc. Its angular proximity to the Virgo overdensity, as well as a good correspondence in distance and metallicity, suggests that the smaller structure we detect may be associated with the much larger Virgo stellar overdensity.

We present the first metallicity distribution functions of the old field populations in the Magellanic Clouds. Our metallicities are based on the Fourier decomposition of Type ab RR Lyrae light curves from the Optical Gravitational Lensing Experiment (OGLE-III). On the metallicity scale of Zinn & West; we find a mean metallicity of [Fe/H] = -1.50 +/- 0.24 dex based on 16776 RR Lyrae stars in the Large Magellanic Cloud (LMC). For the Small Magellanic Cloud (SMC) we obtain -1.70 +/- 0.27 dex based on 1831 RR Lyrae stars. These uncertainties represent the intrinsic spread in the population rather than the standard deviation of the mean. Our results are in good agreement with the few existing spectroscopic metallicity determinations for LMC RR Lyrae stars from the literature. For both the LMC and the SMC the metallicity spread exceeds 1 dex in [Fe/H]. The distribution of metallicities in both Clouds is very uniform, and no significant metallicity gradient is detectable. We also do not find any pronounced populations of extremely metal-poor RR Lyrae candidates with metallicities well below -2 dex, although we need to caution that the photometric method used may overestimate the metallicities of metal-deficient stars. Moreover, because of stellar evolutionary effects one does not expect to observe many RR Lyrae stars among very metal-poor horizontal branch stars. We suggest that the Magellanic Clouds experienced fairly rapid and efficient early enrichment involving pre-enriched gas as well as possibly gas infall, while metal loss through outflows does not seem to have played a significant role. Moreover we suggest that the differences in the metallicities of the old population of LMC and SMC make an origin from a single, common progenitor unlikely, unless the separation happened very early on.

Stellar streams are key players in many aspects of Milky Way studies and, in particular, studying their orbital dynamics is crucial for furthering our understanding of the Milky Way's gravitational potential. Although this is not a trivial task when faced with incomplete dynamical phase-space information, transverse motions of streams can nevertheless be comprehended by harnessing the information contained within their radial velocity gradients. Such methods are not only applicable to stellar streams, but also to HI streams residing in the Milky Way halo. Here, I present the results of two studies that use radial velocity gradients to determine the system's orbit: for Hercules, one of the 'ultra-faint' dwarf galaxies exhibiting a large ellipticity and located at a distance of 140kpc, showing that it may in fact be a stellar stream, and for a string of high-velocity HI clouds belonging to the GCN complex, indicating its likelihood for being a gaseous stream at a distance of approximately 20kpc.

In this work, we extensively investigate the terrestrial planetary formation for the inclined planetary systems (considering the OGLE-2006-BLG-109L system as prototype) in the late stage. In the simulations, we show that the occurrence of terrestrial planets is quite common, in the final assembly stage. Moreover, we find that 40% of the runs finally occupy one planet in the habitable zone (HZ). On the other hand, the numerical results also indicate that the inner region of the planetesimal disk, ranging from $\sim 0.1$ to 0.3 AU, plays an important role in building up terrestrial planets. By examining all simulations, we note that the survivals are located either between 0.1$\sim$1.0 AU or beyond 7 AU, or at the 1:1 mean motion resonance of OGLE-2006-BLG-109Lb at $\sim$2.20 AU. The outcomes suggest that it may exist moderate possibility for the inclined systems to harbor terrestrial planets, even planets in the HZs.

The large number of exoplanets found to orbit their host stars in very close orbits have significantly advanced our understanding of the planetary formation process. It is now widely accepted that such short-period planets cannot have formed {\em in situ}, but rather must have migrated to their current orbits from a formation location much farther from their host star. In the late stages of planetary formation, once the gas in the proto-planetary disk has dissipated and migration has halted, gas-giants orbiting in the inner disk regions will excite planetesimals and planetary embryos, resulting in an increased rate of orbital crossings and large impacts. We present the results of dynamical simulations for planetesimal evolution in this later stage of planet formation. We find that a mechanism is revealed by which the collision-merger of planetary embryos can kick terrestrial planets directly into orbits extremely close to their parent stars.

We investigate the possibility that Hercules, a recently discovered Milky Way (MW) satellite, is a stellar stream in the process of formation. This hypothesis is motivated by Hercules' highly elongated shape as well as the measurement of a tentative radial velocity gradient along its body. The application of simple analytical techniques on radial velocity data of its member stars provides tight constraints on the tangential velocity of the system (v_t = -16^{+6}_{-22} km/s, relative to the Galactic Standard of Rest). Combining this with its large receding velocity (145 km/s) and distance (138 kpc) yields an orbit that would have taken Hercules to within 6^{+9}_{-2} kpc of the Galactic centre approximately 0.6 Gyr ago. This very small perigalacticon can naturally explain the violent tidal destruction of the dwarf galaxy in the MW's gravitational potential, inducing its transformation into a stellar stream.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

We show that a string of HI clouds that form part of the high-velocity cloud complex known as GCN is a probable gaseous stream extending over more than 50 deg in the Galactic halo. The radial velocity gradient along the stream is used to deduce transverse velocities as a function of distance, enabling a family of orbits to be computed. We find that a direction of motion towards the Galactic disk coupled with a mid-stream distance of ~20 kpc provides a good match to the observed sky positions and radial velocities of the HI clouds comprising the stream. With an estimated mass of 10^5 Msun, its progenitor is likely to be a dwarf galaxy. However, no stellar counterpart has been found amongst the currently known Galactic dwarf spheroidal galaxies or stellar streams and the exact origin of the stream is therefore currently unknown.

We investigate the possibility that the recently discovered Hercules Milky Way satellite is in fact a stellar stream in formation, thereby explaining its very elongated shape with an axis ratio of 3 to 1. Under the assumption that Hercules is a stellar stream and that its stars are flowing along the orbit of its progenitor, we find an orbit that would have recently brought the system close enough to the Milky Way to induce its disruption and transformation from a bound dwarf galaxy into a stellar stream. The application of simple analytical techniques to the tentative radial velocity gradient observed in the satellite provides tight constraints on the tangential velocity of the system (v_t = -16^{+6}_{-22} km/s in the Galactic Standard of Rest). Combined with its large receding velocity, the determined tangential velocity yields an orbit with a small pericentric distance (R_peri = 6^{+9}_{-2} kpc). Tidal disruption is therefore a valid scenario for explaining the extreme shape of Hercules. The increase in the mean flattening of dwarf galaxies as one considers fainter systems could therefore be the impact of a few of these satellites not being bound stellar systems dominated by dark matter but, in fact, stellar streams in formation, shedding their stars in the Milky Way's stellar halo.

We have obtained deep photometry in two 1x1 degree fields covering the close pair of dwarf spheroidal galaxies (dSph) Leo IV and Leo V and part of the area in between. We find that both systems are significantly larger than indicated by previous measurements based on shallower data and also significantly elongated. With half-light radii of r_h=4'.6 +- 0'.8 (206 +- 36 pc) and r_h=2'.6 +- 0'.6 (133 +- 31 pc), respectively, they are now well within the physical size bracket of typical Milky Way dSph satellites. Their ellipticities of epsilon ~0.5 are shared by many faint (M_V>-8) Milky Way dSphs. The large spatial extent of our survey allows us to search for extra-tidal features with unprecedented sensitivity. The spatial distribution of candidate red giant branch and horizontal branch stars is found to be non-uniform at the ~3 sigma level. This substructure is aligned along the direction connecting the two systems, indicative of a possible `bridge' of extra-tidal material. Fitting the stellar distribution with a linear Gaussian model yields a significance of 4 sigma for this overdensity, a most likely FWHM of ~16 arcmin and a central surface brightness of ~32 mag arcsec^{-2}. We investigate different scenarios to explain the close proximity of Leo IV and Leo V and the possible tidal bridge between them. Orbit calculations demonstrate that they are unlikely to be remnants of a single disrupted progenitor, while a comparison with cosmological simulations shows that a chance collision between unrelated subhalos is negligibly small. Leo IV and Leo V could, however, be a bound `tumbling pair' if their combined mass exceeds 8 +- 4 x 10^9 M_sun. The scenario of an internally interacting pair appears to be the most viable explanation for this close celestial companionship. (abridged)

We present a simple result in which the distance gradient along a stream can be used to derive the transverse velocity (i.e. proper motion) along it, if the line-of-sight velocity is also known. We show its application to a mock orbit to illustrate its validity and usage. For less extended objects, such as globular clusters and satellite galaxies being tidally disrupted, the same result can be applied in its small-angle approximation. The procedure does not rely on energy or angular momentum conservation and hence does not require a Galactic model in order to deduce the local velocity vector of the stream.

A useful crude approximation for Abelian functions is developed and applied to orbits. The bound orbits in the power-law potentials A*r^{-alpha} take the simple form (l/r)^k = 1 + e cos(m*phi), where k = 2 - alpha > 0 and 'l' and 'e' are generalisations of the semi-latus-rectum and the eccentricity. 'm' is given as a function of 'eccentricity'. For nearly circular orbits 'm' is sqrt{k}, while the above orbit becomes exact at the energy of escape where 'e' is one and 'm' is 'k'. Orbits in the logarithmic potential that gives rise to a constant circular velocity are derived via the limit of small alpha. For such orbits, r^2 vibrates almost harmonically whatever the 'eccentricity'. Unbound orbits in power-law potentials are given in an appendix. The transformation of orbits in one potential to give orbits in a different potential is used to determine orbits in potentials that are positive powers of r. These transformations are extended to form a group which associates orbits in sets of six potentials, e.g. there are corresponding orbits in the potentials proportional to r, r^{-2/3}, r^{-3}, r^{-6}, r^{4/3} and r^{-4}. A degeneracy reduces this to three, which are r^{-1}, r^2 and r^{-4} for the Keplerian case. A generalisation of this group includes the isochrone with the Kepler set.

We present a geometrodynamical method for determining distances to orbital streams of HI gas in the Galaxy. The method makes use of our offset from the Galactic centre and assumes that the gas comprising the stream nearly follows a planar orbit about the Galactic centre. We apply this technique to the Magellanic Stream and determine the distances to all points along it; a consistency check shows that the angular momentum is approximately constant. Applying this technique to the Large Magellanic Cloud itself gives an independent distance which agrees within its accuracy of around 10%. Relaxing the demand for exact conservation of energy and angular momentum at all points along the stream allows for an increase in orbital period between the lagging end and the front end led by the Magellanic Clouds. Similar methods are applicable to other long streams of high-velocity clouds, provided they also nearly follow planar orbits; these would allow otherwise unknown distances to be determined.

We consider the possibility that the Galactic neutral hydrogen stream Complex A and the stellar Orphan stream are related, and use this hypothesis to determine possible distances to Complex A and the Orphan stream, and line-of-sight velocities for the latter. The method presented uses our current knowledge of the projected positions of the streams, as well as line-of-sight velocities for Complex A, and we show that a solution exists in which the two streams share the same orbit. If Complex A and the Orphan stream are on this orbit, our calculations suggest the Orphan stream to be at an average distance of 20 kpc, with heliocentric radial velocities of approximately -110 km/s. Complex A would be ahead of the Orphan stream in the same wrap of the orbit, with an average distance of 10 kpc, which is consistent with the distance constraints determined through interstellar absorption line techniques.

It is now established that several of the Local Group dwarf Spheroidal galaxies (dSphs) have large mass-to-light ratios. We consider the possibility that the dark matter in the halos of dSphs is composed of massive black holes with masses in the range 10^5 to 10^7 solar masses. We use direct N-body simulations to determine the long term evolution of a two-component dSph composed of a pressure-supported stellar population in a black hole halo. The black holes are initially distributed according to a Navarro, Frenk & White profile. For black hole masses between 10^5 and 10^6 solar masses, the dark matter halo evolves towards a shallower inner profile in less than a Hubble time. This suggests that black holes in this mass range might provide an explanation for the origin of the dark matter cores inferred from observations of Low Surface Brightness galaxy rotation curves. We compare the simulated evolution of the stellar population with observed data for the Draco dSph. The dependence of the heating rate on the black hole mass is determined, and an upper limit of 10^5 solar masses is placed on the individual black holes comprising the dark matter halo of Draco, if its present stellar distribution is representative of the initial one. We also present a simple scaling argument which demonstrates that the dynamical heating of an initially compact, though not self-gravitating, stellar distribution might produce a remnant distribution similar in extent to Draco after 10 Gyr, even for black hole masses somewhat in excess of 10^5 solar masses.