Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

Bias due to imperfect shear calibration is the biggest obstacle when constraints on cosmological parameters are to be extracted from large area weak lensing surveys such as Pan-STARRS-3pi, DES or future satellite missions like Euclid. We demonstrate that bias present in existing shear measurement pipelines (e.g. KSB) can be almost entirely removed by means of neural networks. In this way, bias correction can depend on the properties of the individual galaxy instead on being a single global value. We present a procedure to train neural networks for shear estimation and apply this to subsets of simulated GREAT08 RealNoise data. We also show that circularization of the PSF before measuring the shear reduces the scatter related to the PSF anisotropy correction and thus leads to improved measurements, particularly on low and medium signal-to-noise data. Our results are competitive with the best performers in the GREAT08 competition, especially for the medium and higher signal-to-noise sets. Expressed in terms of the quality parameter defined by GREAT08 we achieve a Q = 40, 140 and 1300 without and 50, 200 and 1300 with circularization for low, medium and high signal-to-noise data sets, respectively.

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.

The embedded phase is a crucial period in the development of a young star. Mid-IR observations, now possible with JWST with unprecedented sensitivity, spectral resolution and sharpness are key for probing many physical and chemical processes on sub-arcsecond scales. JOYS addresses a wide variety of questions, from protostellar accretion and the nature of primeval jets, winds and outflows, to the chemistry of gas and ice, and the characteristics of embedded disks. We introduce the program and show representative results. MIRI-MRS data of 17 low-mass and 6 high-mass protostars show a wide variety of features. Atomic line maps differ among refractory (e.g., Fe), semi-refractory (e.g., S) and volatile elements (e.g., Ne), linked to their different levels of depletion and local (shock) conditions. Nested, stratified jet structures consisting of an inner ionized core seen in [Fe II] with an outer H2 layer are commonly seen. Wide-angle winds are found in low-J H2 lines. [S I] follows the jet in the youngest protostars, but is concentrated on source when more evolved. [Ne II] reveals a mix of jet shock and photoionized emission. H I lines measure accretion, but are also associated with jets. Molecular emission (CO2, C2H2, HCN, H2O, ..) is cool compared with disks, and likely associated with hot cores. Deep ice absorption features reveal not just the major ice components but also ions (as part of salts) and complex organic molecules, with comparable abundances from low- to high-mass sources. A second detection of HDO ice in a solar-mass source is presented with HDO/H2O ~ 0.4%, providing a link with disks and comets. A deep search for solid O2 suggests it is not a significant oxygen reservoir. Only few embedded Class I disks show the same forest of water lines as Class II disks do, perhaps due to significant dust extinction of the upper layers [abridged].

We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four $w$CDM parameters ($\Omega_{\rm m}$, $\sigma_8$, $h$, $w_0$). We also determine the joint covariance matrix with over 1000 realisations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalising over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3% precision on $S_8$ for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

The obliquity between the stellar spin axis and the planetary orbit, detected via the Rossiter-McLaughlin (RM) effect, is a tracer of the formation history of planetary systems. While obliquity measurements have been extensively applied to hot Jupiters and short-period planets, they remain rare for cold and long-period planets due to observational challenges, particularly their long transit durations. We report the detection of the RM effect for the 19-hour-long transit of HIP 41378 f, a temperate giant planet on a 542-day orbit, observed through a worldwide spectroscopic campaign. We measure a slight projected obliquity of 21 $\pm$ 8 degrees and a significant 3D spin-orbit angle of 52 $\pm$ 6 degrees, based on the measurement of the stellar rotation period. HIP 41378 f is part of a 5-transiting planetary system with planets close to mean motion resonances. The observed misalignment likely reflects a primordial tilt of the stellar spin axis relative to the protoplanetary disk, rather than dynamical interactions. HIP 41378 f is the first non-eccentric long-period (P>100 days) planet observed with the RM effect, opening new constraints on planetary formation theories. This observation should motivate the exploration of planetary obliquities across a longer range of orbital distances through international collaboration.

The giant radio relic in CIZA J2242.8+5301 is likely evidence of a Mpc sized shock in a massive merging galaxy cluster. However, the exact shock properties are still not clearly determined. In particular, the Mach number derived from the integrated radio spectrum exceeds the Mach number derived from the X-ray temperature jump by a factor of two. We present here a numerical study, aiming for a model that is consistent with the majority of observations of this galaxy cluster. We first show that in the northern shock upstream X-ray temperature and radio data are consistent with each other. We then derive progenitor masses for the system using standard density profiles, X-ray properties and the assumption of hydrostatic equilibrium. We find a class of models that is roughly consistent with weak lensing data, radio data and some of the X-ray data. Assuming a cool-core versus non-cool-core merger, we find a fiducial model with a total mass of $1.6 \times 10^{15}\,M_\odot$, a mass ratio of 1.76 and a Mach number that is consistent with estimates from the radio spectrum. We are not able to match X-ray derived Mach numbers, because even low mass models over-predict the X-ray derived shock speeds. We argue that deep X-ray observations of CIZA J2242.8+5301 will be able to test our model and potentially reconcile X-ray and radio derived Mach numbers in relics.

The bathtub model of the star forming interstellar medium is based on the powerful constraint that mass has to be conserved when gas flows through its various thermal and density phases, ending up eventually in a young star or being blown away by stellar feedback. It predicts that the star formation rate of a molecular cloud is not determined by the cloud's mass or its internal collapse timescale, but rather by the accretion rate of new gas. For the most simple case of a constant accretion flow an equilibrium state is reached quickly where the star formation rate equals the accretion rate and where the dense gas mass is constant and independent of time. The mass of the young star cluster, on the other hand, increases linearly with time. The stellar mass fraction therefore represents a sensitive clock to measure the age of the star-forming region. The bathtub model predicts that the efficiency of star formation is small, of order 1%, even in the dense filamentary phases of molecular clouds. It provides a simple explanation for the dense gas fraction of order 10% in molecular clouds and for the large gas depletion timescales of star-forming galaxies of order $5 \times 10^8 - 10^9$ yrs.

We present a comprehensive analysis of a planetary microlensing event OGLE-2015-BLG-1609. The planetary anomaly was detected by two survey telescopes, OGLE and MOA. Each of these surveys collected enough data over the planetary anomaly to allow for an unambiguous planet detection. Such survey detections of planetary anomalies are needed to build a robust sample of planets that could improve studies on the microlensing planetary occurrence rate by reducing biases and statistical uncertainties. In this work, we examined different methods for modeling microlensing events using individual datasets, particularly we incorporated a Galactic model prior to better constrain poorly defined microlensing parallax. Ultimately, we fitted a comprehensive model to all available data, identifying three potential typologies, with two showing comparably high Bayesian evidence. Our analysis indicates that the host of the planet is a brown dwarf with a probability of 34%, or a low-mass stellar object (M-dwarf) with the probability of 66%.

The bathtub model of the star forming interstellar medium is based on the powerful constraint that mass has to be conserved when gas flows through its various thermal and density phases, ending up eventually in a young star or being blown away by stellar feedback. It predicts that the star formation rate of a molecular cloud is not determined by the cloud's mass or its internal collapse timescale, but rather by the accretion rate of new gas. For the most simple case of a constant accretion flow an equilibrium state is reached quickly where the star formation rate equals the accretion rate and where the dense gas mass is constant and independent of time. The mass of the young star cluster, on the other hand, increases linearly with time. The stellar mass fraction therefore represents a sensitive clock to measure the age of the star-forming region. The bathtub model predicts that the efficiency of star formation is small, of order 1%, even in the dense filamentary phases of molecular clouds. It provides a simple explanation for the dense gas fraction of order 10% in molecular clouds and for the large gas depletion timescales of star-forming galaxies of order $5 \times 10^8 - 10^9$ yrs.

Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies.

The unique characteristics of the brightest group galaxies (BGGs) link the evolutionary continuum between galaxies like the Milky Way and more massive BCGs in dense clusters. This study investigates the stellar properties of BGGs over cosmic time (z = 0.08-1.30), extending our previous work (Gozaliasl et al. 2016, 2018; Paper I and Paper II). We analyze data of 246 BGGs from our X-ray galaxy group catalog in the COSMOS field, examining stellar age, mass, star formation rate (SFR), specific SFR (sSFR), and halo mass. Comparisons are made with Millennium and Magneticum simulations. We explore the variation of stellar properties with the projected offset from the X-ray peak or host halo center. Using a mock galaxy catalog, we evaluated the accuracy of SED-derived stellar ages, finding a mean absolute error of about one Gyr. Observed BGG age distributions show a bias towards younger ages compared to semi-analytical models and the Magneticum simulation. Our analysis of stellar age versus mass reveals trends with a positive slope, suggesting complex evolutionary pathways across redshifts. We observe a negative correlation between stellar age and SFR across all redshift ranges. Using a cosmic-time-dependent main sequence framework, we identify star-forming BGGs, finding that about 20% of BGGs in the local universe exhibit star-forming characteristics, increasing to 50% at $z=1.0$. Our findings support an inside-out formation scenario for BGGs, where older stellar populations are near the X-ray peak and younger populations at larger offsets indicate ongoing star formation. The distribution of stellar ages for lower-mass BGGs ($10^{10-11} M_\odot$) deviates from constant ages predicted by models, highlighting current models' limitations in capturing galaxies' complex star formation histories.

We present a forward-modelling simulation framework designed to model the data products from the Dark Energy Survey (DES). This forward-model process can be thought of as a transfer function -- a mapping from cosmological and astronomical signals to the final data products used by the scientists. Using output from the cosmological simulations (the Blind Cosmology Challenge), we generate simulated images (the Ultra Fast Image Simulator, Berge et al. 2013) and catalogs representative of the DES data. In this work we simulate the 244 sq. deg coadd images and catalogs in 5 bands for the DES Science Verification (SV) data. The simulation output is compared with the corresponding data to show that major characteristics of the images and catalogs can be captured. We also point out several directions of future improvements. Two practical examples, star/galaxy classification and proximity effects on object detection, are then used to demonstrate how one can use the simulations to address systematics issues in data analysis. With clear understanding of the simplifications in our model, we show that one can use the simulations side-by-side with data products to interpret the measurements. This forward modelling approach is generally applicable for other upcoming and future surveys. It provides a powerful tool for systematics studies which is sufficiently realistic and highly controllable.

Gravitational instabilities play an important role in structure formation of gas-rich high-redshift disc galaxies. In this paper, we revisit the axisymmetric perturbation theory and the resulting growth of structure by taking the realistic thickness of the disc into account. In the unstable regime, which corresponds for thick discs to a Toomre parameter below the critical value Q_0,crit = 0.696, we find a fastest growing perturbation wavelength that is always a factor 1.93 times larger than in the classical razor-thin disc approximation. This result is independent of the adopted disc scaleheight and by this independent of temperature and surface density. In order to test the analytical theory, we compare it with a high-resolution hydrodynamical simulation of an isothermal gravitationally unstable gas disc with the typical vertical sech^2 density profile and study its break up into rings that subsequently fragment into dense clumps. In the first phase, rings form, that organize themselves discretely, with distances corresponding to the local fastest growing perturbation wavelength. We find that the disc scaleheight has to be resolved initially with five or more grid cells in order to guarantee proper growth of the ring structures, which follow the analytical prediction. These rings later on contract to a thin and dense line, while at the same time accreting more gas from the inter-ring region. It is these dense, circular filaments, that subsequently fragment into a large number of clumps. Contrary to what is typically assumed, the clump sizes are therefore not directly determined by the fastest growing wavelength.

In this paper we address the question why cluster radio relics that are connected to shock acceleration, so-called radio gischt, have preferentially been found in the outskirts of galaxy clusters. By identifying merger shock waves in cosmological grid simulations, we explore several prescriptions for relating the energy dissipated in shocks to the energy emitted in the radio band. None of the investigated models produce detectable radio relics within 100-200 kpc from the cluster centre. All models cause > 50 per cent of the detectable relic emission at projected distances > 800 kpc. Central radio relics caused by shocks that propagate along the line-of-sight are rare events for simple geometrical reasons, and they have a low surface brightness making them elusive for current instruments. Our simulations show that the radial distribution of observed relics can be explained by the radial trend of dissipated kinetic energy in shocks, that increases with distance from the cluster centre up until half of the virial radius.

We present a distance determination to the Small Magellanic Cloud (SMC) based on an analysis of four detached, long period, late type eclipsing binaries discovered by the OGLE Survey. The components of the binaries show negligible intrinsic variability. A consistent set of stellar parameters was derived with low statistical and systematic uncertainty. The absolute dimensions of the stars are calculated with a precision of better than 3%. The surface brightness - infrared color relation was used to derive the distance to each binary. The four systems clump around a distance modulus of (m - M)=18.99 with a dispersion of only 0.05 mag. Combining these results with the distance published by Graczyk et al. for the eclipsing binary OGLE SMC113.3 4007 we obtain a mean distance modulus to the SMC of 18.965 +/- 0.025 (stat.) +/- 0.048 (syst.) mag. This corresponds to a distance of 62.1 +/- 1.9 kpc, where the error includes both uncertainties. Taking into account other recent published determinations of the SMC distance we calculated the distance modulus difference between the SMC and the LMC equal to 0.458 +/- 0.068 mag. Finally we advocate mu_{SMC}=18.95 +/- 0.07 as a new "canonical" value of the distance modulus to this galaxy.

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

We present a nuclear transient event, PS1-13cbe, that was first discovered in the Pan-STARRS1 survey in 2013. The outburst occurred in the nucleus of the galaxy SDSS J222153.87+003054.2 at $z = 0.12355$, which was classified as a Seyfert 2 in a pre-outburst archival Sloan Digital Sky Survey (SDSS) spectrum. PS1-13cbe showed the appearance of strong broad H$\alpha$ and H$\beta$ emission lines and a non-stellar continuum in a Magellan spectrum taken 57 days after the peak of the outburst that resembled the characteristics of a Seyfert 1. These broad lines were not present in the SDSS spectrum taken a decade earlier and faded away within two years, as observed in several late-time MDM spectra. We argue that the dramatic appearance and disappearance of the broad lines and factor of $\sim 8$ increase in the optical continuum is most likely caused by variability in the pre-existing accretion disk than a tidal disruption event, supernova, or variable obscuration. The timescale for the turn-on of the optical emission of $\sim 70$ days observed in this transient is among the shortest observed in a "changing look" active galactic nucleus.

Over the past 10 Gyr, star-forming galaxies have changed dramatically, from clumpy and gas rich, to rather quiescent stellar-dominated disks with specific star formation rates lower by factors of a few tens. We present a general theoretical model for how this transition occurs, and what physical processes drive it, making use of 1D axisymmetric thin disk simulations with an improved version of the Gravitational Instability-Dominated Galaxy Evolution Tool (GIDGET) code. We show that at every radius galaxies tend to be in a slowly evolving equilibrium state wherein new accretion is balanced by star formation, galactic winds, and radial transport of gas through the disk by gravitational instability (GI) -driven torques. The gas surface density profile is determined by which of these terms are in balance at a given radius - direct accretion is balanced by star formation and galactic winds near galactic centers, and by transport at larger radii. We predict that galaxies undergo a smooth transition from a violent disk instability phase to secular evolution. This model provides a natural explanation for the high velocity dispersions and large clumps in z~2 galaxies, the growth and subsequent quenching of bulges, and features of the neutral gas profiles of local spiral galaxies.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.