AT 2019abn was discovered in the nearby M51 galaxy by the Zwicky Transient Facility at more than two magnitudes and around three weeks prior to its optical peak. We aim to conduct a detailed photometric and spectroscopic follow-up campaign for AT 2019abn, with early discovery allowing for significant pre-maximum observations of an intermediate luminosity red transient (ILRT) for the first time. This work is based on the analysis of u'BVi'z'H photometry and low-resolution spectroscopy using the Liverpool Telescope, medium-resolution spectroscopy with the Gran Telescopio Canarias (GTC), and near-infrared imaging with the GTC and the Nordic Optical Telescope. We present the most detailed optical light curve of an ILRT to date, with multi-band photometry starting around three weeks before peak brightness. The transient peaked at an observed absolute magnitude of M_r' = -13.1, although it is subject to significant reddening from dust in M51, implying an intrinsic M_r' ~ -15.2. The initial light curve showed a linear, achromatic rise in magnitude before becoming bluer at peak. After peak brightness, the transient gradually cooled. This is reflected in our spectra, which at later times show absorption from such species as Fe I, Ni I and Li I. A spectrum taken around peak brightness shows narrow, low-velocity absorption lines, which we interpret as likely to originate from pre-existing circumstellar material. We conclude that while there are some peculiarities, such as the radius evolution, AT 2019abn fits in well overall with the ILRT class of objects and is the most luminous member of the class seen to date.

We ask how the inclusion of various physical heating processes due to the metal content of gas affect the evolution of massive galaxies and compute a suite of cosmological hydrodynamical simulations that follow these systems and their supermassive black holes. We use a smoothed particle hydrodynamics code with a pressure-entropy formulation and a more accurate treatment of the metal production, turbulent diffusion and cooling rate based on individual element abundances. The feedback models include (1) AGN feedback via high velocity BAL winds and Compton/photoionization heating, (2) explicit stellar feedback from multiple processes including powerful winds from supernova events, stellar winds from young massive stars and AGB stars as well as radiative heating within Stromgren spheres around massive stars, and (3) additional heating effects due to the presence of metals including grain photoelectric heating, metallicity dependent X-ray heating by nearby accreting black holes and from the cosmic X-ray background, which are the major improvement in our feedback model. With a suite of zoom-in simulations of 30 halos with $M_{vir} \sim 10^{12-13.4}$, we show that energy and momentum budget from all feedback effects generate realistic galaxy properties. We explore the detailed role of each feedback model with three additional sets of simulations with varying input physics. We show that the metal induced heating mechanisms reduce the fraction of accreted stellar material by mainly suppressing the growth of diffuse small stellar systems at high redshift but overall have a relatively minor effect on the final stellar and gas properties of massive galaxies. The inclusion of AGN feedback significantly improves the ability of our cosmological simulations to yield realistic gas and stellar properties of massive galaxies with reasonable fraction of the final stellar mass accreted from other galaxies.

Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{\mathrm{hot}}$, as $(1 - f_{\mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Str\"omgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$\times 10^4$ solar masses.

We present the results of a spectroscopic survey to characterize chromospheric activity, as measured by H$\alpha$ emission, in low-mass members of the 500-Myr-old open cluster M37. Combining our new measurements of H$\alpha$ luminosities ($L_{H\alpha}$) with previously cataloged stellar properties, we identify saturated and unsaturated regimes in the dependence of the $L_{H\alpha}$-to-bolometric-luminosity ratio, $L_{H\alpha}/L_{bol}$, on the Rossby number $R_o$. All rotators with Ro smaller than 0.03$\pm$0.01 converge to an activity level of $L_{H\alpha}/L_{bol}$ = (1.27$\pm$0.02) x 10$^{-4}$. This saturation threshold ($R_{o,sat}$ = 0.03$\pm$0.01) is statistically smaller than that found in most studies of the rotation-activity relation. In the unsaturated regime, slower rotators have lower levels of chromospheric activity, with $L_{H\alpha}/L_{bol}$($R_o$) following a power-law of index $\beta$ = -0.51$\pm$0.02, slightly shallower than the one found for a combined $\approx$650-Myr-old sample of Hyades and Praesepe stars. By comparing this unsaturated behavior to that previously found for coronal activity in M37 (as measured via the X-ray luminosity, $L_X$), we confirm that chromospheric activity decays at a much slower rate than coronal activity with increasing $R_o$. While a comparison of $L_{H\alpha}$ and $L_X$ for M37 members with measurements of both reveals a nearly 1:1 relation, removing the mass-dependencies by comparing instead $L_{H\alpha}/L_{bol}$ and $L_X/L_{bol}$ does not provide clear evidence for such a relation. Finally, we find that $R_{o,sat}$ is smaller for our chromospheric than for our coronal indicator of activity ($R_{o,sat}$ = 0.03$\pm$0.01 versus 0.09$\pm$0.01). We interpret this as possible evidence for coronal stripping.

Empirical calibrations of the stellar age-rotation-activity relation (ARAR) rely on observations of the co-eval populations of stars in open clusters. We used the Chandra X-ray Observatory to study M37, a 500-Myr-old open cluster that has been extensively surveyed for rotation periods ($P_{\rm rot}$). M37 was observed almost continuously for five days, for a total of 440.5 ksec, to measure stellar X-ray luminosities ($L_{\mathrm{X}}$), a proxy for coronal activity, across a wide range of masses. The cluster's membership catalog was revisited to calculate updated membership probabilities from photometric data and each star's distance to the cluster center. The result is a comprehensive sample of 1699 M37 members: 426 with $P_{\rm rot}$, 278 with X-ray detections, and 76 with both. We calculate Rossby numbers, $R_o = P_{\rm rot}/\tau$, where $\tau$ is the convective turnover time, and ratios of the X-ray-to-bolometric luminosity, $L_{\rm X}/L_{\rm bol}$, to minimize mass dependencies in our characterization of the rotation-coronal activity relation at 500 Myr. We find that fast rotators, for which $R_o<0.09\pm0.01$, show saturated levels of activity, with log($L_{\rm X}/L_{\rm bol}$)$=-3.06\pm0.04$. For $R_o\geq0.09\pm0.01$, activity is unsaturated and follows a power law of the form $R_o^\beta$, where $\beta$=$-2.03_{-0.14}^{+0.17}$. This is the largest sample available for analyzing the dependence of coronal emission on rotation for a single-aged population, covering stellar masses in the range 0.4$-$1.3 $M_{\odot}$, $P_{\rm rot}$ in the range 0.4$-$12.8 d, and $L_{\rm X}$ in the range 10$^{28.4-30.5}$ erg s$^{-1}$. Our results make M37 a new benchmark open cluster for calibrating the ARAR at ages of $\approx$500 Myr.

We use a 440.5 ks Chandra observation of the $\approx$500-Myr-old open cluster M37 to derive the X-ray luminosity functions of its $\leq1.2$ $M_{\odot}$ stars. Combining detections of 162 M37 members with upper limits for 160 non-detections, we find that its G, K, and M stars have a similar median (0.5$-$7 keV) X-ray luminosity L$_X =10^{29.0}$ erg/s, whereas the L$_X$-to-bolometric-luminosity ratio (L$_X$/L$_{bol}$) indicates that M stars are more active than G and K stars by $\approx$1 order of magnitude at 500 Myr. To characterize the evolution of magnetic activity in low-mass stars over their first $\approx$600 Myr, we consolidate X-ray and optical data from the literature for stars in six other open clusters: from youngest to oldest, the Orion Nebula Cluster (ONC), NGC 2547, NGC 2516, the Pleiades, NGC 6475, and the Hyades. For these, we homogenize the conversion of instrumental count rates to L$_X$ by applying the same one-temperature emission model as for M37, and obtain masses using the same empirical mass-absolute magnitude relation (except for the ONC). We find that for G and K stars X-ray activity decreases $\approx$2 orders of magnitude over their first 600 Myr, and for M stars, $\approx$1.5. The decay rate of the median L$_X$ follows the relation L$_X \propto~t^b$, where $b=-0.61\pm0.12$ for G, $-0.82\pm0.16$ for K, and $-0.40\pm0.17$ for M stars. In L$_X$/L$_{bol}$ space, the slopes are $-0.68\pm0.12$, $-0.81\pm0.19$, and $-0.61\pm0.12$, respectively. These results suggest that for low-mass stars the age-activity relation steepens after $\approx$625 Myr, consistent with the faster decay in activity observed in solar analogs at $t>1$ Gyr.

We present parallax and proper motion measurements, near-infrared spectra, and WISE photometry for the low surface gravity L5gamma dwarf 2MASSJ035523.37+113343.7 (2M0355). We use these data to evaluate photometric, spectral, and kinematic signatures of youth as 2M0355 is the reddest isolated L dwarf yet classified. We confirm its low-gravity spectral morphology and find a strong resemblance to the sharp triangular shaped $H$-band spectrum of the 10 Myr planetary-mass object 2M1207b. We find that 2M0355 is underluminous compared to a normal field L5 dwarf in the optical and MKO J,H, and K bands and transitions to being overluminous from 3-12 microns, indicating that enhanced photospheric dust shifts flux to longer wavelengths for young, low-gravity objects, creating a red spectral energy distribution. Investigating the near-infrared color magnitude diagram for brown dwarfs confirms that 2M0355 is redder and underluminous compared to the known brown dwarf population, similar to the peculiarities of directly imaged exoplanets 2M1207b and HR8799bcd. We calculate UVW space velocities and find that the motion of 2M0355 is consistent with young disk objects (< 2-3 Gyr) and it shows a high likelihood of membership in the AB Doradus association.