Search

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6., Comment: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

We present results from the James Webb Space Telescope (JWST) Director's Discretionary Time Early Release Science (ERS) program 1328 targeting the nearby, Luminous Infrared Galaxy (LIRG), VV 114. We use the MIRI and NIRSpec instruments to obtain integral-field spectroscopy of the heavily obscured Eastern nucleus (V114E) and surrounding regions. The spatially resolved, high-resolution, spectra reveal the physical conditions in the gas and dust over a projected area of 2-3 kpc that includes the two brightest IR sources, the NE and SW cores. Our observations show for the first time spectroscopic evidence that the SW core hosts an AGN as evidenced by its very low 6.2 {\mu}m and 3.3 {\mu}m PAH equivalent widths (0.12 and 0.017 {\mu}m respectively) and mid and near-IR colors. Our observations of the NE core show signs of deeply embedded star formation including absorption features due to aliphatic hydrocarbons, large quantities of amorphous silicates, as well as HCN due to cool gas along the line of sight. We detect elevated [Fe II]/Pf{\alpha} consistent with extended shocks coincident with enhanced emission from warm H$_{2}$, far from the IR-bright cores and clumps. We also identify broadening and multiple kinematic components in both H$_{2}$ and fine structure lines caused by outflows and previously identified tidal features., Comment: 12 pages, 4 figures, accepted for publication in ApJL

We present mid-infrared spectroscopic observations of the nucleus of the nearby Seyfert galaxy NGC 7469 taken with the MIRI instrument on the James Webb Space Telescope (JWST) as part of Directors Discretionary Time Early Release Science (ERS) program 1328. The high resolution nuclear spectrum contains 19 emission lines covering a wide range of ionization. The high ionization lines show broad, blueshifted emission reaching velocities up to 1700 km s$^{-1}$ and FWHM ranging from $\sim500 - 1100$ km s$^{-1}$. The width of the broad emission and the broad to narrow line flux ratios correlate with ionization potential. The results suggest a decelerating, stratified, AGN driven outflow emerging from the nucleus. The estimated mass outflow rate is one to two orders of magnitude larger than the current black hole accretion rate needed to power the AGN. Eight pure rotational H$_{2}$ emission lines are detected with intrinsic widths ranging from FWHM $\sim 125-330$ km s$^{-1}$. We estimate a total mass of warm H$_{2}$ gas of $\sim1.2\times10^{7}$M$_{\odot}$ in the central 100 pc. The PAH features are extremely weak in the nuclear spectrum, but a $6.2\mu$m PAH feature with an equivalent width $\sim0.07\mu$m and a flux of $2.7\times10^{-17}$ W m$^{-2}$ is detected. The spectrum is steeply rising in the mid-infrared, with a silicate strength $\sim0.02$, significantly smaller than seen in most PG QSOs, but comparable to other Seyfert 1's. These early MIRI mid-infrared IFU data highlight the power of JWST to probe the multi-phase interstellar media surrounding actively accreting supermassive black holes.

We present the analysis of $\sim 100$pc-scale compact radio continuum sources detected in 63 local (Ultra) Luminous Infrared Galaxies (U/LIRGs; $L_{\rm IR} \ge 10^{11} L_\odot$), using FWHM $\lesssim 0''.1 - 0''.2$ resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8 - 170pc, which are classified into four main categories -- "AGN" (AGN), "AGN/SBnuc" (AGN-starburst composite nucleus), "SBnuc" (starburst nucleus) and "SF" (star-forming clumps) -- based on ancillary datasets and the literature. We find that "AGN" and "AGN/SBnuc" more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with "SBnuc" and "SF", which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities ($\Sigma_{\rm SFR}$) are measured for "SF" and "SBnuc" using both the total 33 GHz continuum emission (SFR $\sim 0.14 - 13$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR} \sim 13 - 1600$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$) and the thermal free-free emission from HII regions (median SFR$_{\rm th} \sim 0.4$ M$_\odot$ yr$^{-1}$, $\Sigma_{\rm SFR_{th}} \sim 44$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$). These values are 1 - 2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have much flatter median 15 - 33 GHz spectral index ($\sim -0.08$) compared with "SBnuc" and "SF" ($\sim -0.46$), which may reflect higher non-thermal contribution from supernovae and/or ISM densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100\,pc scales., Comment: Accepted for publication in ApJ

Nuclear rings are excellent laboratories for studying intense star formation. We present results from a study of nuclear star-forming rings in five nearby normal galaxies from the Star Formation in Radio Survey (SFRS) and four local LIRGs from the Great Observatories All-sky LIRG Survey (GOALS) at sub-kpc resolutions using VLA high-frequency radio continuum observations. We find that nuclear ring star formation (NRSF) contributes 49 - 60\% of the total star formation of the LIRGs, compared to 7 - 40\% for the normal galaxies. We characterize a total of 58 individual star-forming regions in these rings, and find that with measured sizes of 10 - 200 pc, NRSF regions in the LIRGs have SFR and $\Sigma_\mathrm{SFR}$ up to 1.7 M$_\odot$yr$^{-1}$ and 402 M$_\odot$yr$^{-1}$kpc$^{-2}$, respectively, which are about 10 times higher than NRSF regions in the normal galaxies with similar sizes, and comparable to lensed high-$z$ star-forming regions. At $\sim 100 - 300$ pc scales, we estimate low contributions ($< 50\%$) of thermal free-free emission to total radio continuum emission at 33 GHz in the NRSF regions in the LIRGs, but large variations possibly exist at smaller physical scales. Finally, using archival sub-kpc resolution CO (J=1-0) data of nuclear rings in the normal galaxies and NGC 7469 (LIRG), we find a large scatter in gas depletion times at similar molecular gas surface densities, which tentatively points to a multi-modal star formation relation on sub-kpc scales., Comment: Accepted for publication in ApJ

We present a numerical method which evolves a two-temperature, magnetized, radiative, accretion flow around a black hole, within the framework of general relativistic radiation magnetohydrodynamics. As implemented in the code KORAL, the gas consists of two sub-components -- ions and electrons -- which share the same dynamics but experience independent, relativistically consistent, thermodynamical evolution. The electrons and ions are heated independently according to a standard prescription from the literature for magnetohydrodynamical turbulent dissipation. Energy exchange between the particle species via Coulomb collisions is included. In addition, electrons gain and lose energy and momentum by absorbing and emitting synchrotron and bremsstrahlung radiation, and through Compton scattering. All evolution equations are handled within a fully covariant framework in the relativistic fixed-metric spacetime of the black hole. Numerical results are presented for five models of low luminosity black hole accretion. In the case of a model with a mass accretion rate $\dot{M}\sim10^{-8} \dot M_{\rm Edd}$, we find that radiation has a negligible effect on either the dynamics or the thermodynamics of the accreting gas. In contrast, a model with a larger $\dot{M}\sim 4\times 10^{-4} \dot M_{\rm Edd}$ behaves very differently. The accreting gas is much cooler and the flow is geometrically less thick, though it is not quite a thin accretion disk., Comment: 21 pages, improved numerical method, MNRAS, in press

A new general relativistic radiation magnetohydrodynamical code KORAL, is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a=0.0) and a spinning black hole (a=0.9) are presented. The accretion rates in the two models are \dot M = 100-200 \dot M_Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a=0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a=0) and 33% (a=0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a=0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only \sim 1-10 L_Edd, which corresponds to a radiative efficiency \sim 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington., Comment: 19 pages, 13 figures, submitted to MNRAS

The expulsion of axisymmetric magnetic field from the event horizons of rapidly rotating black holes has been seen as an astrophysically important effect that may significantly reduce or even nullify the efficiency of the Blandford-Znajek mechanism of powering the relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts. However, this Meissner-like effect is seen in vacuum solutions of black hole electrodynamics whereas the Blandford-Znajek mechanism is concerned with plasma-filled magnetospheres. In this paper we argue that conductivity dramatically changes the properties of axisymmetric electromagnetic solutions -- even for a maximally rotating Kerr black hole the magnetic field is pulled inside the event horizon. Moreover, the conditions resulting in outgoing Poynting flux in the Blandford-Znajek mechanism exist not on the event horizon but everywhere within the black hole ergosphere. Thus, the Meissner effect is unlikely to be of interest in astrophysics of black holes, at least not in the way this has been suggested so far. These conclusions are supported by the results of time-dependent numerical simulations with three different computer codes. The test problems involve black holes with the rotation parameter ranging from $a=0.999$ to $a=1$. The pure electrodynamic simulations deal with the structure of conductive magnetospheres of black holes placed in a uniform-at-infinity magnetic field (Wald's problem) and the magnetohydrodynamic simulations are used to study the magnetospheres arising in the problem of disk accretion., Comment: Accepted by MNRAS Letters

Results indicated going steady among high school sophomores was associated with lower self-esteem. Lowest scores were held by girls who had gone steady in the past and were still going steady. For boys the difference in self-esteem between those who were going steady and those who were not was unrelated to their past dating practices. (BF)