Search

Galaxy evolution emerges from the balance between cosmic gas accretion, fueling star formation, and supernova (SN) feedback, regulating the metal enrichment. Hence, the stellar mass ($M_*$) - gas metallicity relation (MZR) is key to understand the physics of galaxies. High-quality JWST data enable accurate measurements of the MZR up to redshift z=10. Our aims are to understand the observed MZR, its connection with the star formation rate (SFR), the role played by SFR stochasticity, and how it is regulated by SN feedback. We compare the MZR from the JADES, CEERS, and UNCOVER surveys, which comprise about 180 galaxies at $z=3-10$ with $10^6

The gas in the interstellar medium (ISM) of galaxies is supersonically turbulent. Measurements of turbulence typically rely on cold gas emission lines for low-z galaxies and warm ionized gas observations for z>0 galaxies. Studies of warm gas kinematics at z>0 conclude that the turbulence strongly evolves as a function of redshift, due to the increasing impact of gas accretion and mergers in the early Universe. However, recent findings suggest potential biases in turbulence measurements derived from ionized gas at high-z, impacting our understanding of turbulence origin, ISM physics and disk formation. We investigate the evolution of turbulence using velocity dispersion ($\sigma$) measurements from cold gas tracers (i.e., CO, [CI], [CII]) derived from a sample of 57 galaxy disks spanning the redshift range z=0-5. This sample consists of main-sequence and starburst galaxies with stellar masses $\gtrsim 10^{10} M_{\odot}$. The comparison with current H$\alpha$ kinematic observations and existing models demonstrates that the velocity dispersion inferred from cold gas tracers differ by a factor of $\approx 3$ from those obtained using emission lines tracing warm gas. We show that stellar feedback is the main driver of turbulence measured from cold gas tracers. This is fundamentally different from the conclusions of studies based on warm gas, which had to consider additional turbulence drivers to explain the high values of $\sigma$. We present a model predicting the redshift evolution of turbulence in galaxy disks, attributing the increase of $\sigma$ with redshift to the higher energy injected by supernovae due to the elevated star-formation rate in high-z galaxies. This supernova-driven model suggests that turbulence is lower in galaxies with lower stellar mass compared to those with higher stellar mass. Additionally, it forecasts the evolution of $\sigma$ in Milky-Way like progenitors., Comment: Accepted for publication in A&A. The abstract has been modified to comply with arXiv's character limit

The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R\sim100$) spectroscopic observations. These are among the highest-$z$ sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kpc, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log($M_*/M_\odot) \sim 7.6-9.3$, star formation rates (SFRs) $\sim1-15$ $M_\odot$ yr$^{-1}$, and gas-phase metallicities 12+log(O/H) $\sim 7.7-8.3$, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-$z$ samples to date ($z\sim3-10$), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the $z\sim3-10$ sources. Moreover, the targets in this work follow the predicted slope of the MZR at $z\sim 6-8$ from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and SN winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at $z\sim6-8$ sets also important constraints on cosmological simulations and chemical evolution models, whose predictions on the cosmic evolution of metallicity gradients differ significantly, but are mostly limited to $z<3$ so far., Comment: 15 pages (+5 appendix); 7 figures (+8 in appendix); in production on A&A

Early galaxy formation, initiated by the dark matter and gas assembly, evolves through frequent mergers and feedback processes into dynamically hot, chaotic structures. In contrast, dynamically cold, smooth rotating disks have been observed in massive evolved galaxies merely 1.4 billion years after the Big Bang, suggesting rapid morphological and dynamical evolution in the early Universe. Probing this evolution mechanism necessitates studies of young galaxies, yet efforts have been hindered by observational limitations in both sensitivity and spatial resolution. Here we report high-resolution observations of a strongly lensed and quintuply imaged, low-luminosity, young galaxy at $z=6.072$ (dubbed the Cosmic Grapes), 930 million years after the Big Bang. Magnified by gravitational lensing, the galaxy is resolved into at least 15 individual star-forming clumps with effective radii of $r_{\rm e}\simeq$ 10--60 parsec (pc), which dominate $\simeq$ 70\% of the galaxy's total flux. The cool gas emission unveils a smooth, underlying rotating disk characterized by a high rotational-to-random motion ratio and a gravitationally unstable state (Toomre $Q \simeq$ 0.2--0.3), with high surface gas densities comparable to local dusty starbursts with $\simeq10^{3-5}$ $M_{\odot}$/pc$^{2}$. These gas properties suggest that the numerous star-forming clumps are formed through disk instabilities with weak feedback effects. The clumpiness of the Cosmic Grapes significantly exceeds that of galaxies at later epochs and the predictions from current simulations for early galaxies. Our findings shed new light on internal galaxy substructures and their relation to the underlying dynamics and feedback mechanisms at play during their early formation phases, potentially explaining the high abundance of bright galaxies observed in the early Universe and the dark matter core-cusp problem., Comment: Submitted. 44 pages, 14 figures, 5 tables. Comments are welcome! See also the companion papers on arXiv. Valentino+2024: arXiv:2402.17845 Gim\'enez-Arteaga+2024: arXiv:2402.17875

Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and dust morphology of a z~3.4 lensed galaxy hosting four clumps detected in the UV continuum. We aimed to observe the [CII] emission of individual clumps that, unlike the UV, is not affected by dust extinction, to probe their nature and cold gas content. We conducted ALMA observations probing scales down to ~300 pc and detected three [CII] clumps. One (dubbed "NE") coincides with the brightest UV clump, while the other two ("SW" and "C") are not detected in the UV continuum. We do not detect the dust continuum. We converted the [CII] luminosity of individual clumps into molecular gas mass and found Mmol~10^8 Msun. By complementing it with the star formation rate (SFR) estimate from the UV continuum, we estimated the gas depletion time (tdep) of clumps and investigated their location in the Schmidt-Kennicutt plane. While the NE clump has a short tdep=0.16 Gyr, comparable with high-redshift starbursts, the SW and C clumps instead have longer tdep>0.65 Gyr and are likely probing the initial phases of star formation. The lack of dust continuum detection is consistent with the blue UV continuum slope estimated for this galaxy (beta~-2.5) and it indicates that dust inhomogeneities do not significantly affect the detection of UV clumps in this target. We pushed the observation of the cold gas content of individual clumps up to z~3.4 and showed that the [C II] line emission is a promising tracer of molecular clouds at high redshift, allowing the detection of clumps with a large range of depletion times., Comment: Accepted for publication in A&A

A sizable fraction of the heavy elements synthesized by stars in galaxies condenses into sub-micron-sized solid-state particles, known as dust grains. Dust produces a wavelength-dependent attenuation, $A_\lambda$, of the galaxy emission, thereby significantly altering its observed properties. Locally, $A_\lambda$ is in general the sum of a power-law and a UV feature ('bump') produced by small, carbon-based grains. However, scant information exists regarding its evolution across cosmic time. Here, leveraging data from 173 galaxies observed by the James Webb Space Telescope in the redshift range z = 2 - 12, we report the most distant detection of the UV bump in a z ~ 7.55 galaxy (when the Universe was only ~ 700 Myr old), and show for the first time that the power-law slope and the bump strength decrease towards high redshifts. We propose that the flat $A_\lambda$ shape at early epochs is produced by large grains newly formed in supernova ejecta, which act as the main dust factories at such early epochs. Importantly, these grains have undergone minimal reprocessing in the interstellar medium due to the limited available cosmic time. This discovery opens new perspectives in the study of cosmic dust origin and evolution., Comment: 16 pages, 5 figures

High-redshift galaxies are expected to be more turbulent than local galaxies because of their smaller size and higher star formation and thus stronger feedback from star formation, frequent mergers events, and gravitational instabilities. However, this scenario has recently been questioned by the observational evidence of a few galaxies at z~4-5 with a gas velocity dispersion similar to what is observed in the local population. Our goal is to determine whether galaxies in the first Gyrs of the Universe have already formed a dynamically cold rotating disk similar to the local counterparts. We studied the gas kinematic of 22 main-sequence star-forming galaxies at z > 5 and determined their dynamical state by estimating the ratio of the rotational velocity and of the gas velocity dispersion. We mined the ALMA archive and exploited the [CII] and [OIII] observations to perform a kinematic analysis of the cold and warm gas of z>5 main-sequence galaxies. The gas kinematics of the high-z galaxies is consistent within the errors with rotating but turbulent disks. We infer a velocity dispersion that is systematically higher by 4 times than the local galaxy population and the z~5 dust-obscured galaxies reported in the literature. The difference between our results and those reported at similar redshift can be ascribed to the systematic difference in the galaxy properties in the two samples: the disks of massive dusty galaxies are dynamically colder than the disks of dust-poor galaxies. The comparison with the theoretical predictions suggests that the main driver of the velocity dispersion in high-z galaxies is the gravitational energy that is released by the transport of mass within the disk. Finally, we stress that future deeper ALMA high-angular resolution observations are crucial to constrain the kinematic properties of high-z galaxies and to distinguish rotating disks from kpc-scale mergers., Comment: 14 pages, 11 figures, 1 tables, accepted for publication in A&A

Spatially-resolved studies of the kinematics of galaxies provide crucial insights into their assembly and evolution, enabling to infer the properties of the dark matter halos, derive the impact of feedback on the ISM, characterize the outflow motions. To date, most of the kinematic studies at z=0.5-3.5 were obtained using emission lines tracing the warm, ionized gas. However, whether these provide an exhaustive or only a partial view of the dynamics of galaxies and of the properties of the ISM is still debated. Complementary insights on the cold gas kinematics are therefore needed. We present ALPAKA, a project aimed at gathering high-resolution observations of CO and [CI] emission lines of star-forming galaxies at z=0.5-3.5 from the ALMA public archive. With 147 hours of total integration time, ALPAKA assembles ~0.25'' observations for 28 star-forming galaxies, the largest sample with spatially-resolved cold gas kinematics as traced by either CO or [CI] at z>0.5. By combining multi-wavelength ancillary data, we derive the stellar masses ($M_{\star}$) and star-formation rates (SFR) for our targets, finding values of $M_{\star}\gtrsim 10^{10}$ M$_{\odot}$ and SFR of 10-3000 M$_{\odot}$/yr. A large fraction of ALPAKA galaxies (19/28) lie in overdense regions (clusters, groups, and protoclusters). We exploit the ALMA data to infer their dynamical state and we find that 19/28 ALPAKA galaxies are rotating disks, 2 are interacting systems, while for the remaining 7 sources the classification is uncertain. The disks have velocity dispersion values that are typically larger in the innermost regions than in the outskirts, with a median value for the entire disk sample of 35$^{+11}_{-9}$ km/s. Despite the bias of our sample towards galaxies hosting very energetic mechanisms, the ALPAKA disks have high ratios of ordered-to-random motion ($V/\sigma$) with a median value of 9$^{+7}_{-2}$., Comment: 36 pages, 23 figures, 5 tables; accepted for publication in A&A. The data and the outputs of the kinematic analysis will be made available at https://alpaka-survey.github.io/index.html

ALMA observations have revealed that many high redshift galaxies are surrounded by extended (10-15 kpc) [CII]-emitting halos which are not predicted by even the most advanced zoom-in simulations. Using a semi-analytical model, in a previous work we suggested that such halos are produced by starburst-driven, catastrophically cooling outflows. Here, we further improve the model and compare its predictions with data from 7 star-forming ($10\lesssim \rm SFR/ M_\odot \rm yr^{-1}<100$) galaxies at z=4-6, observed in the ALPINE survey. We find that (a) detected [CII] halos are a natural by-product of starburst-driven outflows; (b) the outflow mass loading factors are in the range $4\lesssim\eta\lesssim 7$, with higher $\eta$ values for lower-mass, lower-SFR systems, and scale with stellar mass as $\eta \propto M_*^{-0.43}$, consistently with the momentum-driven hypothesis. Our model suggests that outflows are widespread phenomena in high-z galaxies. However, in low-mass systems the halo extended [CII] emission is likely too faint to be detected with the current levels of sensitivity., Comment: Accepted for publication in MNRAS. 11 pages, 6 figures

Recent attempts to detect [OIII] 88$\mu$m emission from super-early ($z>10$) galaxy candidates observed by JWST have been unsuccessful. By using zoom-in simulations, we show that these galaxies are faint, and mostly fall below the local metal-poor $\rm [OIII]-SFR$ relation as a result of their low ionization parameter, $U_{\rm ion}\lesssim 10^{-3}$. Such low $U_{\rm ion}$ values are found in galaxies that are in an early assembly stage, and whose stars are still embedded in high-density natal clouds. However, the most luminous galaxy in our sample ($\rm{log}[L_{\rm{[OIII]}}/L_\odot] = 8.4$, $U_{\rm ion} \approx 0.1$) could be detected by ALMA in only $2.8$ hrs., Comment: 6 pages, 3 figures, Submitted to MNRAS (comments are welcome)

The characterization of the dynamical state of galaxies up to z~7 is crucial for constraining the mechanisms driving the mass assembly in the early Universe. However, it is unclear whether the data quality of current and future observations is sufficient to perform a solid dynamical analysis. This paper defines the angular resolution and S/N required for a robust characterization of the dynamical state of galaxies up to the EoR. The final aim is to help design spatially-resolved surveys targeting emission lines of primeval galaxies. We investigate the [CII]-158um emission from z~6-7 LBGs from the SERRA cosmological simulation, covering a range of dynamical states: from disks to major mergers. We create ALMA mock observations with various data quality and apply the kinematic classification methods used in the literature. These tests allow us to quantify the performances of such methods as a function of angular resolution and S/N. We find that barely-resolved observations do not allow the correct dynamical characterization of a galaxy, resulting in the misclassification of all disks in our sample. However, even when using spatially-resolved observations with data quality typical of high-z galaxies, the standard kinematic classification methods, based on the analysis of the moment maps, fail to distinguish a merger from a disk. The high angular resolution and S/N needed to apply these standard methods successfully can be achieved with current data only for a handful of bright galaxies. We propose a new classification method, called PVsplit, that quantifies the asymmetries and morphological features in position-velocity diagrams using three empirical parameters. We test PVsplit on our mock data concluding that it can predict whether a galaxy is a disk or a merger provided that S/N $\gtrsim10$, and the major axis is covered by $\gtrsim3$ independent resolution elements., Comment: Submitted to Astronomy and Astrophysics (A&A) Journal. Comments are welcome

We introduce SERRA, a suite of zoom-in high-resolution ($\sim 10\,\rm pc$) cosmological simulations including non-equilibrium chemistry and on-the-fly radiative transfer. The outputs are post-processed to derive galaxy UV+FIR continuum and emission line properties. Results are compared with available multi-wavelength data to constrain the physical properties (e.g., star formation rates, stellar/gas/dust mass, metallicity) of high-redshift $6 \lesssim z \lesssim 15$ galaxies. This flagship paper focuses on the $z=7.7$ sub-sample, including 202 galaxies with stellar mass $10^7 M_\odot \lesssim M_\star \lesssim 5\times 10^{10}M_\odot$, and specific star formation ranging from ${\rm sSFR} \sim 100\,{\rm Gyr}^{-1}$ in young, low-mass galaxies to $\sim 10\,{\rm Gyr}^{-1}$ for older, massive ones. At this redshift, SERRA galaxies are typically bursty, i.e. they are located above the Schmidt-Kennicutt relation by a factor $\kappa_s = 3.03^{+4.9}_{-1.8}$, consistent with recent findings for [OIII] and [CII] emitters at high-$z$. They also show relatively large ${\rm IRX} = L_{\rm FIR}/L_{\rm UV}$ values as a result of their compact/clumpy morphology effectively blocking the stellar UV luminosity. Note that this conclusion might be affected by insufficient spatial resolution at the molecular cloud level. We confirm that early galaxies lie on the standard $\rm [CII]-SFR$ relation; their observed $L_{\rm [OIII]}/L_{\rm [CII]} \simeq 1-10$ ratios can be reproduced by a part of the SERRA galaxies without the need of a top-heavy IMF and/or anomalous C/O abundances. [OI] line intensities are similar to local ones, making ALMA high-$z$ detections challenging but feasible ($\sim 6\,\rm hr$ for a SFR of $50\,M_\odot\,{\rm yr}^{-1}$)., Comment: 23 pages, 14 figures, accepted by MNRAS; part of the data used for this study is available at the website http://cosmology.sns.it/data_access.html

We study the structure of spatially resolved, line-of-sight velocity dispersion for galaxies in the Epoch of Reionization (EoR) traced by [CII] $158\mu\rm{m}$ line emission. Our laboratory is a simulated prototypical Lyman-break galaxy, "Freesia", part of the SERRA suite. The analysis encompasses the redshift range 6 < z < 8, when Freesia is in a very active assembling phase. We build velocity dispersion maps for three dynamically distinct evolutionary stages (Spiral Disk at z=7.4, Merger at z=8.0, and Disturbed Disk at z=6.5) using [CII] hyperspectral data cubes. We find that, at a high spatial resolution of 0.005" ($\simeq 30 pc$), the luminosity-weighted average velocity dispersion is $\sigma_{\rm{CII}}$~23-38 km/s with the highest value belonging to the highly-structured Disturbed Disk stage. Low resolution observations tend to overestimate $\sigma_{\rm CII}$ values due to beam smearing effects that depend on the specific galaxy structure. For an angular resolution of 0.02" (0.1"), the average velocity dispersion is 16-34% (52-115%) larger than the actual one. The [CII] emitting gas in Freesia has a Toomre parameter $\mathcal{Q}$~0.2 and a rotational-to-dispersion ratio of $v_{\rm c}/\sigma$~ 7 similar to that observed in z=2-3 galaxies. The primary energy source for the velocity dispersion is due to gravitational processes, such as merging/accretion events; energy input from stellar feedback is generally subdominant (< 10%). Finally, we find that the resolved $\sigma_{\rm{CII}} - {\Sigma}_{\rm SFR}$ relation is relatively flat for $0.02<{\Sigma}_{\rm SFR}/{{\rm M}_{\odot}} \mathrm{yr}^{-1} {\mathrm kpc}^{-2} < 30$, with the majority of data lying on the derived analytical relation $\sigma \propto \Sigma_{\rm SFR}^{5/7}$. At high SFR, the increased contribution from stellar feedback steepens the relation, and $\sigma_{\rm{CII}}$ rises slightly., Comment: Accepted for publication in MNRAS; 16 pages (plus appendix), 8 figures, 3 tables

We investigate the spatially-resolved morphology of galaxies in the early Universe. We consider a typical redshift z = 6 Lyman Break galaxy, "Althaea" from the SERRA hydrodynamical simulations. We create mock rest-frame ultraviolet, optical, and far-infrared observations, and perform a two-dimensional morphological analysis to de-blend the galaxy disk from substructures (merging satellites or star-forming regions). We find that the [CII]158um emitting region has an effective radius 1.5 - 2.5 times larger than the optical one, consistent with recent observations. This [CII] halo in our simulated galaxy arises as the joint effect of stellar outflows and carbon photoionization by the galaxy UV field, rather than from the emission of unresolved nearby satellites. At the typical angular resolution of current observations (> 0.15") only merging satellites can be detected; detection of star-forming regions requires resolutions of < 0.05". The [CII]-detected satellite has a 2.5 kpc projected distance from the galaxy disk, whereas the star-forming regions are embedded in the disk itself (distance < 1 kpc). This suggests that multi-component systems reported in the literature, which have separations > 2 kpc, are merging satellites, rather than galactic substructures. Finally, the star-forming regions found in our mock maps follow the local L[CII] - SFR_UV relation of galaxy disks, although sampling the low-luminosity, low-SFR tail of the distribution. We show that future JWST observations, bridging UV and [CII] datasets, will be exceptionally suited to characterize galaxy substructures thanks to their exquisite spatial resolution and sensitivity to both low-metallicity and dust-obscured regions that are bright at infrared wavelengths., Comment: Accepted for publication in MNRAS; 17 pages (plus appendix), 7 figures, 4 tables

ALMA observations have revealed that [CII] 158$\mu$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($\Delta^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations.

A tight relation between the [CII]158$\mu$m line luminosity and star formation rate is measured in local galaxies. At high redshift ($z>5$), though, a much larger scatter is observed, with a considerable (15-20\%) fraction of the outliers being [CII]-deficient. Moreover, the [CII] surface brightness ($\Sigma_{\rm CII}$) of these sources is systematically lower than expected from the local relation. To clarify the origin of such [CII]-deficiency we have developed an analytical model that fits local [CII] data, and has been validated against radiative transfer simulations performed with CLOUDY. The model predicts an overall increase of $\Sigma_{\rm CII}$ with the surface star formation rate ($\Sigma_*$). However, for $\Sigma_* > 1 M_\odot~{\rm yr}^{-1}~{\rm kpc}^{-2}$, $\Sigma_{\rm CII}$ saturates. We conclude that underluminous [CII] systems can result from a combination of three factors: (a) large upward deviations from the Kennicutt-Schmidt relation ($\kappa_s \gg 1$), parameterized by the "burstiness" parameter $\kappa_s$; (b) low metallicity; (c) low gas density, at least for the most extreme sources (e.g. CR7). Observations of [CII] emission alone cannot break the degeneracy among the above three parameters; this requires additional information coming from other emission lines (e.g. [OIII]88$\mu$m, CIII]1909A, CO lines). Simple formulae are given to interpret available data for low and high-$z$ galaxies., Comment: Accepted for publication in MNRAS. 12 pages, 6 Figures. Comments welcome

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies., Comment: 20 pages (plus appendix), 8 figures, 2 tables. Accepted for publication in MNRAS

We study the kinematical properties of galaxies in the Epoch of Reionization via the [CII] 158$\mu$m line emission. The line profile provides information on the kinematics as well as structural properties such as the presence of a disk and satellites. To understand how these properties are encoded in the line profile, first we develop analytical models from which we identify disk inclination and gas turbulent motions as the key parameters affecting the line profile. To gain further insights, we use "Althaea", a highly-resolved ($30\, \rm pc$) simulated prototypical Lyman Break Galaxy, in the redshift range $z = 6-7$, when the galaxy is in a very active assembling phase. Based on morphology, we select three main dynamical stages: I) Merger , II) Spiral Disk, and III) Disturbed Disk. We identify spectral signatures of merger events, spiral arms, and extra-planar flows in I), II), and III), respectively. We derive a generalised dynamical mass vs. [CII]-line FWHM relation. If precise information on the galaxy inclination is (not) available, the returned mass estimate is accurate within a factor $2$ ($4$). A Tully-Fisher relation is found for the observed high-$z$ galaxies, i.e. $L_{\rm[CII]}\propto (FWHM)^{1.80\pm 0.35}$ for which we provide a simple, physically-based interpretation. Finally, we perform mock ALMA simulations to check the detectability of [CII]. When seen face-on, Althaea is always detected at $> 5\sigma$; in the edge-on case it remains undetected because the larger intrinsic FWHM pushes the line peak flux below detection limit. This suggests that some of the reported non-detections might be due to inclination effects., Comment: 14 pages, 12 figures, accepted for publication in MNRAS

We study the formation and evolution of a sample of Lyman Break Galaxies in the Epoch of Reionization by using high-resolution ($\sim 10 \,{\rm pc}$), cosmological zoom-in simulations part of the SERRA suite. In SERRA, we follow the interstellar medium (ISM) thermo-chemical non-equilibrium evolution, and perform on-the-fly radiative transfer of the interstellar radiation field (ISRF). The simulation outputs are post-processed to compute the emission of far infrared lines ([CII], [NII], and [OIII]). At $z=8$, the most massive galaxy, `Freesia', has an age $t_\star \simeq 409\,{\rm Myr}$, stellar mass $M_{\star} \simeq 4.2\times 10^9 {\rm M}_{\odot}$, and a star formation rate ${\rm SFR} \simeq 11.5\,{\rm M}_{\odot}{\rm yr}^{-1}$, due to a recent burst. Freesia has two stellar components (A and B) separated by $\simeq 2.5\, {\rm kpc}$; other 11 galaxies are found within $56.9 \pm 21.6 \, {\rm kpc}$. The mean ISRF in the Habing band is $G = 7.9\, G_0$ and is spatially uniform; in contrast, the ionisation parameter is $U = 2^{+20}_{-2} \times 10^{-3}$, and has a patchy distribution peaked at the location of star-forming sites. The resulting ionising escape fraction from Freesia is $f_{\rm esc}\simeq 2\%$. While [CII] emission is extended (radius 1.54 kpc), [OIII] is concentrated in Freesia-A (0.85 kpc), where the ratio $\Sigma_{\rm [OIII]}/\Sigma_{\rm [CII]} \simeq 10$. As many high-$z$ galaxies, Freesia lies below the local [CII]-SFR relation. We show that this is the general consequence of a starburst phase (pushing the galaxy above the Kennicutt-Schmidt relation) which disrupts/photodissociates the emitting molecular clouds around star-forming sites. Metallicity has a sub-dominant impact on the amplitude of [CII]-SFR deviations., Comment: 22 pages, 14 figures, accepted by MNRAS

In this paper we study anisotropic spherical polytropes within the framework of general relativity. Using the anisotropic Tolman-Oppenheimer-Volkov (TOV) equations, we explore the relativistic anisotropic Lane-Emden equations. We find how the anisotropic pressure affects the boundary conditions of these equations. Also we argue that the behaviour of physical quantities near the center of star changes in the presence of anisotropy. For constant density, a class of exact solution is derived with the aid of a new ansatz and its physical properties are discussed., Comment: 4 figures

In this paper we are concerned with the effects of anisotropic pressure on the boundary conditions of anisotropic Lane-Emden equation and homology theorem. Some new exact solutions of this equation are derived. Then some of the theorems governing the Newtonian perfect fluid star are extended taking the anisotropic pressure into account.

The tumour control probability (TCP) is a formalism derived to compare various treatment regimens of radiation therapy, defined as the probability that given a prescribed dose of radiation, a tumour has been eradicated or controlled. In the traditional view of cancer, all cells share the ability to divide without limit and thus have the potential to generate a malignant tumour. However, an emerging notion is that only a sub-population of cells, the so-called cancer stem cells (CSCs), are responsible for the initiation and maintenance of the tumour. A key implication of the CSC hypothesis is that these cells must be eradicated to achieve cures, thus we define TCP_S as the probability of eradicating CSCs for a given dose of radiation. A cell surface protein expression profile, such as CD44high/CD24low for breast cancer, is often used as a biomarker to monitor CSCs enrichment. However, it is increasingly recognized that not all cells bearing this expression profile are necessarily CSCs, and in particular early generations of progenitor cells may share the same phenotype. Thus, due to the lack of a perfect biomarker for CSCs, we also define a novel measurable TCP_CD+, that is the probability of eliminating or controlling biomarker positive cells. Based on these definitions, we use stochastic methods and numerical simulations to compare the theoretical TCP_S and the measurable TCP_CD+. We also use the measurable TCP to compare the effect of various radiation protocols., Comment: 18 pages, 5 figures

The combination of cytotoxic therapies and anti-angiogenic agents is emerging as a most promising strategy in the treatment of malignant tumors. However, the timing and sequencing of these treatments seem to play essential roles in achieving a synergic outcome. Using a mathematical modeling approach that is grounded on available experimental data, we investigate the spatial and temporal targeting of tumor cells and neovasculature with a nanoscale delivery system. Our model suggests that the experimental success of the nanoscale delivery system depends crucially on the trapping of chemotherapeutic agents within the tumor tissue. The numerical results also indicate that substantial further improvements in the efficiency of the nanoscale delivery system can be achieved through an adjustment of the temporal targeting mechanism., Comment: 17 pages, 5 figures