Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the \textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $b\approx0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias.

Shells are low surface brightness tidal debris that appear as interleaved caustics with large opening angles, often situated on both sides of the galaxy center. In this paper, we study the incidence and formation processes of shell galaxies in the cosmological gravity+hydrodynamics Illustris simulation. We identify shells at redshift z=0 using stellar surface density maps, and we use stellar history catalogs to trace the birth, trajectory and progenitors of each individual star particle contributing to the tidal feature. Out of a sample of the 220 most massive galaxies in Illustris ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_{\odot}$), $18\%\pm3\%$ of the galaxies exhibit shells. This fraction increases with increasing mass cut: higher mass galaxies are more likely to have stellar shells. Furthermore, the fraction of massive galaxies that exhibit shells decreases with increasing redshift. We find that shell galaxies observed at redshift $z=0$ form preferentially through relatively major mergers ($\gtrsim$1:10 in stellar mass ratio). Progenitors are accreted on low angular momentum orbits, in a preferred time-window between $\sim$4 and 8 Gyrs ago. Our study indicates that, due to dynamical friction, more massive satellites are allowed to probe a wider range of impact parameters at accretion time, while small companions need almost purely radial infall trajectories in order to produce shells. We also find a number of special cases, as a consequence of the additional complexity introduced by the cosmological setting. These include galaxies with multiple shell-forming progenitors, satellite-of-satellites also forming shells, or satellites that fail to produce shells due to multiple major mergers happening in quick succession.

Stellar shells are low surface brightness arcs of overdense stellar regions, extending to large galactocentric distances. In a companion study, we identified 39 shell galaxies in a sample of 220 massive ellipticals ($\mathrm{M}_{\mathrm{200crit}}>6\times10^{12}\,\mathrm{M}_\odot$) from the Illustris cosmological simulation. We used stellar history catalogs to trace the history of each individual star particle inside the shell substructures, and we found that shells in high-mass galaxies form through mergers with massive satellites (stellar mass ratios $\mu_{\mathrm{stars}}\gtrsim1:10$). Using the same sample of shell galaxies, the current study extends the stellar history catalogs in order to investigate the metallicity of stellar shells around massive galaxies. Our results indicate that outer shells are often times more metal-rich than the surrounding stellar material in a galaxy's halo. For a galaxy with two different satellites forming $z=0$ shells, we find a significant difference in the metallicity of the shells produced by each progenitor. We also find that shell galaxies have higher mass-weighted logarithmic metallicities ([Z/H]) at $2$-$4\,\mathrm{R}_{\mathrm{eff}}$ compared to galaxies without shells. Our results indicate that observations comparing the metallicities of stars in tidal features, such as shells, to the average metallicities in the stellar halo can provide information about the assembly histories of galaxies.

Self-gravitating bosonic fields can support stable and localised field configurations. For real fields, these solutions oscillate in time and are known as oscillatons. The density profile is static, and is soliton. Such solitons should be ubiquitous in models of axion dark matter, with the soliton characteristic mass and size depending on some inverse power of the axion mass. Stable configurations of non-relativistic axions are studied numerically using the Schr\"{o}dinger-Poisson system. This method, and the resulting soliton density profiles, are reviewed. Using a scaling symmetry and the uncertainty principle, the core size of the soliton can be related to the central density and axion mass, $m_a$, in a universal way. Solitons have a constant central density due to pressure-support, unlike the cuspy profile of cold dark matter (CDM). One consequence of this fact is that solitons composed of ultra-light axions (ULAs) may resolve the `cusp-core' problem of CDM. In DM halos, thermodynamics will lead to a CDM-like Navarro-Frenk-White profile at large radii, with a central soliton core at small radii. Using Monte-Carlo techniques to explore the possible density profiles of this form, a fit to stellar-kinematical data of dwarf spheroidal galaxies is performed. In order for ULAs to resolve the cusp-core problem (without recourse to baryon feedback or other astrophysical effects) the axion mass must satisfy $m_a<1.1\times 10^{-22}\text{ eV}$ at 95\% C.L. On the other hand, ULAs with $m_a\lesssim 1\times 10^{-22}\text{ eV}$ are in some tension with cosmological structure formation. An axion solution to the cusp-core problem thus makes novel predictions for future measurements of the epoch of reionisation. On the other hand, this can be seen as evidence that structure formation could soon impose a \emph{Catch 22} on axion/scalar field DM, similar to the case of warm DM.

We use kinetic hybrid simulations (kinetic ions - fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.

We use high-resolution cosmological hydrodynamic galaxy formation simulations to gain insights on how galaxies lose their cold gas at low redshift as they migrate from the field to the high density regions of clusters of galaxies. We find that beyond three cluster virial radii, the fraction of gas-rich galaxies is constant, representing the field. Within three cluster-centric radii, the fraction of gas-rich galaxies declines steadily with decreasing radius, reaching $\mathbf{<10%}$ near the cluster center. Our results suggest that galaxies start to feel the impact of the cluster environment on their gas content well beyond the cluster virial radius. We show that almost all gas-rich galaxies at the cluster virial radius are falling in for the first time at nearly radial orbits. Furthermore, we find that almost no galaxy moving outward at the cluster virial radius is gas-rich (with gas to baryon ratio greater than $\mathbf{1%}$). These results suggest that galaxies that fall into clusters lose their cold gas within a single radial round-trip.

We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented.