We investigate mass ejection from accretion disks formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH-NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disk configurations expected from BH-NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of black hole accretion disks that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower-limit to the fraction of the initial disk mass ejected. We find a nearly linear inverse dependence of this fraction on disk compactness (BH mass over initial disk radius). The dependence is related to the fraction of the disk mass accreted at the time when weak interactions freeze out. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disk mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disk masses. We provide estimates of the radioactive luminosity from the disk outflow alone available to power kilonovae over the range of BH and NS combinations studied, finding a spread of two orders of magnitude. For most of the BH-NS parameter space, the disk outflow contribution is well below the kilonova mass upper limits for S190814bv.

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

Neutron star mergers are unique laboratories of accretion, ejection, and r-process nucleosynthesis. We used 3D general relativistic magnetohydrodynamic simulations to study the role of the post-merger magnetic geometry in the evolution of merger remnant discs around stationary Kerr black holes. Our simulations fully capture mass accretion, ejection, and jet production, owing to their exceptionally long duration exceeding $4$ s. Poloidal post-merger magnetic field configurations produce jets with energies $E_\mathrm{jet} \sim (4{-}30)\times10^{50}$ erg, isotropic equivalent energies $E_\mathrm{iso}\sim(4{-}20)\times10^{52}$ erg, opening angles $\theta_\mathrm{jet}\sim6{-}13^\circ$, and durations $t_j\lesssim1$ s. Accompanying the production of jets is the ejection of $f_\mathrm{ej}\sim30{-}40\%$ of the post-merger disc mass, continuing out to times $> 1$ s. We discover that a more natural, purely toroidal post-merger magnetic field geometry generates large-scale poloidal magnetic flux of alternating polarity and striped jets. The first stripe, of $E_\mathrm{jet}\simeq2\times10^{48}\,\mathrm{erg}$, $E_\mathrm{iso}\sim10^{51}$ erg, $\theta_\mathrm{jet}\sim3.5{-}5^\circ$, and $t_j\sim0.1$ s, is followed by $\gtrsim4$ s of striped jet activity with $f_\mathrm{ej}\simeq27\%$. The dissipation of such stripes could power the short gamma-ray burst (sGRB) prompt emission. Our simulated jet energies and durations span the range of sGRBs. We find that although the blue kilonova component is initially hidden from view by the red component, it expands faster, outruns the red component, and becomes visible to off-axis observers. In comparison to GW 170817/GRB 170817A, our simulations under-predict the mass of the blue relative to red component by a factor of few. Including the dynamical ejecta and neutrino absorption may reduce this tension.

We study mass ejection from accretion disks formed in the merger of a white dwarf with a neutron star or black hole. These disks are mostly radiatively-inefficient and support nuclear fusion reactions, with ensuing outflows and electromagnetic transients. Here we perform time-dependent, axisymmetric hydrodynamic simulations of these disks including a physical equation of state, viscous angular momentum transport, a coupled $19$-isotope nuclear network, and self-gravity. We find no detonations in any of the configurations studied. Our global models extend from the central object to radii much larger than the disk. We evolve these global models for several orbits, as well as alternate versions with an excised inner boundary to much longer times. We obtain robust outflows, with a broad velocity distribution in the range $10^2-10^4$ km s$^{-1}$. The outflow composition is mostly that of the initial white dwarf, with burning products mixed in at the $\lesssim 10-30\%$ level by mass, including up to $\sim 10^{-2}M_\odot$ of ${}^{56}$Ni. These heavier elements (plus ${}^{4}$He) are ejected within $\lesssim 40^\circ$ of the rotation axis, and should have higher average velocities than the lighter elements that make up the white dwarf. These results are in broad agreement with previous one- and two-dimensional studies, and point to these systems as progenitors of rapidly-rising ($\sim $ few day) transients. If accretion onto the central BH/NS powers a relativistic jet, these events could be accompanied by high energy transients with peak luminosities $\sim 10^{47}-10^{50}$ erg s$^{-1}$ and peak durations of up to several minutes, possibly accounting for events like CDF-S XT2.

We study mass ejection from accretion disks around newly-formed hypermassive neutron stars (HMNS). Standard kilonova model fits to GW170817 require at least a lanthanide-poor ('blue') and lanthanide-rich ('red') component. The existence of a blue component has been used as evidence for a HMNS remnant of finite lifetime, but average disk outflow velocities from existing long-term HMNS simulations fall short of the inferred value ($\sim 0.25c$) by a factor of $\sim 2$. Here we use time-dependent, axisymmetric hydrodynamic simulations of HMNS disks to explore the limits of the model and its ability to account for observations. For physically plausible parameter choices compatible with GW170817, we find that hydrodynamic models that use shear viscosity to transport angular momentum cannot eject matter with mass-averaged velocities larger than $\sim 0.15c$. While outflow velocities in our simulations can exceed the asymptotic value for a steady-state neutrino-driven wind, the increase in the average velocity due to viscosity is not sufficient. Therefore, viscous HMNS disk winds cannot reproduce by themselves the ejecta properties inferred from multi-component fits to kilonova light curves from GW170817. Three possible resolutions remain feasible within standard merger ejecta channels: more sophisticated radiative transfer models that allow for photon reprocessing between ejecta components, inclusion of magnetic stresses, or enhancement of the dynamical ejecta. We provide fits to our disk outflow models once they reach homologous expansion.

We investigate the long-term evolution of black hole accretion disks formed in neutron star mergers. These disks expel matter that contributes to an $r$-process kilonova, and can produce relativistic jets powering short gamma-ray bursts. Here we report the results of a three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulation of such a disk which is evolved for long enough ($\sim 9$s, or $\sim 6\times 10^5 r_{\rm g}/c$) to achieve completion of mass ejection far from the disk. Our model starts with a poloidal field, and fully resolves the most unstable mode of the magnetorotational instability. We parameterize the dominant microphysics and neutrino cooling effects, and compare with axisymmetric hydrodynamic models with shear viscosity. The GRMHD model ejects mass in two ways: a prompt MHD-mediated outflow and a late-time, thermally-driven wind once the disk becomes advective. The total amount of unbound mass ejected ($0.013M_\odot$, or $\simeq 40\%$ of the initial torus mass) is twice as much as in hydrodynamic models, with higher average velocity ($0.1c$) and a broad electron fraction distribution with a lower average value ($0.16$). Scaling the ejected fractions to a disk mass of $\sim 0.1M_\odot$ can account for the red kilonova from GW170817 but underpredicts the blue component. About $\sim 10^{-3}M_\odot$ of material should undergo neutron freezout and could produce a bright kilonova precursor in the first few hours after the merger. With our idealized initial magnetic field configuration, we obtain a robust jet and sufficient ejecta with Lorentz factor $\sim 1-10$ to (over)produce the non-thermal emission from GW1708107.

We examine the excitation of unstable magnetosonic waves in the radiative envelopes of intermediate- and high-mass stars with a magnetic field of ~kG strength. Wind clumping close to the star and microturbulence can often be accounted for when including small-scale, sub-photospheric density or velocity perturbations. Compressional waves - with wavelengths comparable to or shorter than the gas pressure scale height - can be destabilized by the radiative flux in optically-thick media when a magnetic field is present, in a process called the Radiation-Driven Magneto-Acoustic Instability (RMI). The instability does not require radiation or magnetic pressure to dominate over gas pressure, and acts independently of sub-surface convection zones. Here we evaluate the conditions for the RMI to operate on a grid of stellar models covering a mass range $3-40M_\odot$ at solar metallicity. For a uniform 1kG magnetic field, fast magnetosonic modes are unstable down to an optical depth of a few tens, while unstable slow modes extend beyond the depth of the iron convection zone. The qualitative behavior is robust to magnetic field strength variations by a factor of a few. When combining our findings with previous results for the saturation amplitude of the RMI, we predict velocity fluctuations in the range ~0.1-10 km/s. These amplitudes are a monotonically increasing function of the ratio of radiation to gas pressure, or alternatively, of the zero-age main sequence mass.

We investigate the cosmic evolution of the absolute and specific star formation rate (SFR, sSFR) of galaxies as derived from a spatially-resolved study of the stellar populations in a set of 366 nearby galaxies from the CALIFA survey. The analysis combines GALEX and SDSS images with the 4000 break, H_beta, and [MgFe] indices measured from the datacubes, to constrain parametric models for the SFH, which are then used to study the cosmic evolution of the star formation rate density (SFRD), the sSFR, the main sequence of star formation (MSSF), and the stellar mass density (SMD). A delayed-tau model, provides the best results, in good agreement with those obtained from cosmological surveys. Our main results from this model are: a) The time since the onset of the star formation is larger in the inner regions than in the outer ones, while tau is similar or smaller in the inner than in the outer regions. b) The sSFR declines rapidly as the Universe evolves, and faster for early than for late type galaxies, and for the inner than for the outer regions of galaxies. c) SFRD and SMD agree well with results from cosmological surveys. At z< 0.5, most star formation takes place in the outer regions of late spiral galaxies, while at z>2 the inner regions of the progenitors of the current E and S0 are the major contributors to SFRD. d) The inner regions of galaxies are the major contributor to SMD at z> 0.5, growing their mass faster than the outer regions, with a lookback time at 50% SMD of 9 and 6 Gyr for the inner and outer regions. e) The MSSF follows a power-law at high redshift, with the slope evolving with time, but always being sub-linear. f) In agreement with galaxy surveys at different redshifts, the average SFH of CALIFA galaxies indicates that galaxies grow their mass mainly in a mode that is well represented by a delayed-tau model, with the peak at z~2 and an e-folding time of 3.9 Gyr.

We study the ejection of mass during stellar core-collapse when the stalled shock does not revive and a black hole forms. Neutrino emission during the protoneutron star phase causes a decrease in the gravitational mass of the core, resulting in an outward going sound pulse that steepens into a shock as it travels out through the star. We explore the properties of this mass ejection mechanism over a range of stellar progenitors using spherically-symmetric, time-dependent hydrodynamic simulations that treat neutrino mass loss parametrically and follow the shock propagation over the entire star. We find that all types of stellar progenitor can eject mass through this mechanism. The ejected mass is a decreasing function of the surface gravity of the star, ranging from several $M_\odot$ for red supergiants to $\sim 0.1M_\odot$ for blue supergiants and $\sim 10^{-3} M_\odot $ for Wolf-Rayet stars. We find that the final shock energy at the surface is a decreasing function of the core-compactness, and is $\lesssim 10^{47}-10^{48}$ erg in all cases. In progenitors with a sufficiently large envelope, high core-compactness, or a combination of both, the sound pulse fails to unbind mass. Successful mass ejection is accompanied by significant fallback accretion that can last from hours to years. We predict the properties of shock breakout and thermal plateau emission produced by the ejection of the outer envelope of blue supergiant and Wolf-Rayet progenitors in otherwise failed supernovae.

During the core collapse of massive stars, the formation of the protoneutron star is accompanied by the emission of a significant amount of mass-energy ($\sim 0.3 \, M_{\odot}$) in the form of neutrinos. This mass-energy loss generates an outward-propagating pressure wave that steepens into a shock near the stellar surface, potentially powering a weak transient associated with an otherwise-failed supernova. We analytically investigate this mass-loss-induced wave generation and propagation. Heuristic arguments provide an accurate estimate of the amount of energy contained in the outgoing sound pulse. We then develop a general formalism for analyzing the response of the star to centrally concentrated mass loss in linear perturbation theory. To build intuition, we apply this formalism to polytropic stellar models, finding qualitative and quantitative agreement with simulations and heuristic arguments. We also apply our results to realistic pre-collapse massive star progenitors (both giants and compact stars). Our analytic results for the sound pulse energy, excitation radius, and steepening in the stellar envelope are in good agreement with full time-dependent hydrodynamic simulations. We show that {prior} to the sound pulses arrival at the stellar photosphere, the photosphere has already reached velocities $\sim 20-100 \%$ of the local sound speed, thus likely modestly decreasing the stellar effective temperature prior to the star disappearing. Our results provide important constraints on the physical properties and observational appearance of failed supernovae.

We investigate the nucleosynthesis of heavy elements in the winds ejected by accretion disks formed in neutron star mergers. We compute the element formation in disk outflows from hypermassive neutron star (HMNS) remnants of variable lifetime, including the effect of angular momentum transport in the disk evolution. We employ long-term axisymmetric hydrodynamic disk simulations to model the ejecta, and compute r-process nucleosynthesis with tracer particles using a nuclear reaction network containing $\sim 8000$ species. We find that the previously known strong correlation between HMNS lifetime, ejected mass, and average electron fraction in the outflow is directly related to the amount of neutrino irradiation on the disk, which dominates mass ejection at early times in the form of a neutrino-driven wind. Production of lanthanides and actinides saturates at short HMNS lifetimes ($\lesssim 10$ ms), with additional ejecta contributing to a blue optical kilonova component for longer-lived HMNSs. We find good agreement between the abundances from the disk outflow alone and the solar r-process distribution only for short HMNS lifetimes ($\lesssim 10$ ms). For longer lifetimes, the rare-earth and third r-process peaks are significantly under-produced compared to the solar pattern, requiring additional contributions from the dynamical ejecta. The nucleosynthesis signature from a spinning black hole (BH) can only overlap with that from a HMNS of moderate lifetime ($\lesssim 60$ ms). Finally, we show that angular momentum transport not only contributes with a late-time outflow component, but that it also enhances the neutrino-driven component by moving material to shallower regions of the gravitational potential, in addition to providing additional heating.

This paper presents the mass assembly time scales of nearby galaxies observed by CALIFA at the 3.5m telescope in Calar Alto. We apply the fossil record method of the stellar populations to the complete sample of the 3rd CALIFA data release, with a total of 661 galaxies, covering stellar masses from 10$^{8.4}$ to 10$^{12}$ M$_{\odot}$ and a wide range of Hubble types. We apply spectral synthesis techniques to the datacubes and process the results to produce the mass growth time scales and mass weighted ages, from which we obtain temporal and spatially resolved information in seven bins of galaxy morphology and six bins of stellar mass (M$_{\star}$) and stellar mass surface density ($\Sigma_{\star}$). We use three different tracers of the spatially resolved star formation history (mass assembly curves, ratio of half mass to half light radii, and mass-weighted age gradients) to test if galaxies grow inside-out, and its dependence with galaxy stellar mass, $\Sigma_{\star}$, and morphology. Our main results are as follows: (a) The innermost regions of galaxies assemble their mass at an earlier time than regions located in the outer parts; this happens at any given M$_{\star}$, $\Sigma_{\star}$, or Hubble type, including the lowest mass systems. (b) Galaxies present a significant diversity in their characteristic formation epochs for lower-mass systems. This diversity shows a strong dependence of the mass assembly time scales on $\Sigma_{\star}$ and Hubble type in the lower-mass range (10$^{8.4}$ to 10$^{10.4}$), but a very mild dependence in higher-mass bins. (c) All galaxies show negative $\langle$log age$\rangle_{M}$ gradients in the inner 1 HLR. The profile flattens with increasing values of $\Sigma_{\star}$. There is no significant dependence on M$_{\star}$ within a particular $\Sigma_{\star}$ bin, except for the lowest bin, where the gradients becomes steeper.

The Calar Alto Legacy Integral Field Area (CALIFA) survey, a pioneer in integral field spectroscopy legacy projects, has fostered many studies exploring the information encoded on the spatially resolved data on gaseous and stellar features in the optical range of galaxies. We describe a value-added catalogue of stellar population properties for CALIFA galaxies analysed with the spectral synthesis code STARLIGHT and processed with the PyCASSO platform. Our public data base (http://pycasso.ufsc.br/, mirror at http://pycasso.iaa.es/) comprises 445 galaxies from the CALIFA Data Release 3 with COMBO data. The catalogue provides maps for the stellar mass surface density, mean stellar ages and metallicities, stellar dust attenuation, star formation rates, and kinematics. Example applications both for individual galaxies and for statistical studies are presented to illustrate the power of this data set. We revisit and update a few of our own results on mass density radial profiles and on the local mass-metallicity relation. We also show how to employ the catalogue for new investigations, and show a pseudo Schmidt-Kennicutt relation entirely made with information extracted from the stellar continuum. Combinations to other databases are also illustrated. Among other results, we find a very good agreement between star formation rate surface densities derived from the stellar continuum and the $\mathrm{H}\alpha$ emission. This public catalogue joins the scientific community's effort towards transparency and reproducibility, and will be useful for researchers focusing on (or complementing their studies with) stellar properties of CALIFA galaxies.

This paper presents the spatially resolved star formation history (2D-SFH) of a small sample of four local mergers: the early-stage mergers IC1623, NGC6090, and the Mice, and the more advanced merger NGC2623, by analyzing IFS data from the CALIFA survey and PMAS in LArr mode. Full spectral fitting techniques are applied to the datacubes to obtain the spatially resolved mass growth histories, the time evolution of the star formation rate intensity ($\Sigma_{SFR}$), and the local specific star formation rate, over three different timescales (30 Myr, 300 Myr, and 1 Gyr). The results are compared with non-interacting Sbc--Sc galaxies, to quantify if there is an enhancement of the star formation and to trace its time scale and spatial extent. Our results for the three LIRGs (IC1623W, NGC6090, and NGC2623) show that a major phase of star formation is occurring in time scales of 10$^{7}$ yr to few 10$^{8}$ yr, with global SFR enhancements of $\sim$2--6 with respect to main-sequence star forming (MSSF) galaxies. In the two early-stage mergers IC1623W and NGC6090, which are between first pericenter passage and coalescence, the most remarkable increase of the SFR with respect to non-interacting spirals occurred in the last 30 Myr, and it is spatially extended, with enhancements of factors 2--7 both in the centres ($r <$ 0.5 half light radius, HLR), and in the disks ($r >$ 1 HLR). In the more advanced merger NGC 2623 an extended phase of star formation occurred on a longer time-scale of $\sim$1 Gyr, with a SFR enhancement of a factor $\sim$2--3 larger than the one in Sbc--Sc MSSF galaxies over the same period, probably relic of the first pericenter passage epoch. A SFR enhancement in the last 30 Myr is also present, but only in NGC2623 centre, by a factor 3. In general, the spatially resolved SFHs of the LIRG-mergers are consistent with the predictions from high spatial resolution simulations.

We investigate the ejecta from black hole - neutron star mergers by modeling the formation and interaction of mass ejected in a tidal tail and a disk wind. The outflows are neutron-rich, giving rise to optical/infrared emission powered by the radioactive decay of $r$-process elements (a kilonova). Here we perform an end-to-end study of this phenomenon, where we start from the output of a fully-relativistic merger simulation, calculate the post-merger hydrodynamical evolution of the ejecta and disk winds including neutrino physics, determine the final nucleosynthetic yields using post-processing nuclear reaction network calculations, and compute the kilonova emission with a radiative transfer code. We study the effects of the tail-to-disk mass ratio by scaling the tail density. A larger initial tail mass results in fallback matter becoming mixed into the disk and ejected in the subsequent disk wind. Relative to the case of a disk without dynamical ejecta, the combined outflow has lower mean electron fraction, faster speed, larger total mass, and larger absolute mass free of high-opacity Lanthanides or Actinides. In most cases, the nucleosynthetic yield is dominated by the heavy $r$-process contribution from the unbound part of the tidal tail. A Solar-like abundance distribution can however be obtained when the total mass of the dynamical ejecta is comparable to the mass of the disk outflows. The kilonova has a characteristic duration of 1 week and a luminosity of ~$10^{41}$ erg/s, with orientation effects leading to variations of a factor ~2 in brightness. At early times (< 1 day) the emission includes an optical component from the (hot) Lanthanide-rich material, but the spectrum evolves quickly to the infrared thereafter.

This paper presents the spatially resolved star formation history (SFH) of nearby galaxies with the aim of furthering our understanding of the different processes involved in the formation and evolution of galaxies. To this end, we apply the fossil record method of stellar population synthesis to a rich and diverse data set of 436 galaxies observed with integral field spectroscopy in the CALIFA survey. The sample covers a wide range of Hubble types, with stellar masses ranging from $M_\star \sim 10^9$ to $7 \times 10^{11} M_\odot$. Spectral synthesis techniques are applied to the datacubes to retrieve the spatially resolved time evolution of the star formation rate (SFR), its intensity ($\Sigma_{\rm SFR}$), and other descriptors of the 2D-SFH in seven bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc, and Sd), and five bins of stellar mass. Our main results are: a) Galaxies form very fast independently of their current stellar mass, with the peak of star formation at high redshift ($z > 2$). Subsequent star formation is driven by $M_\star$ and morphology, with less massive and later type spirals showing more prolonged periods of star formation. b) At any epoch in the past the SFR is proportional to $M_\star$, with most massive galaxies having the highest absolute (but lowest specific) SFRs. c) While nowadays $\Sigma_{\rm SFR}$ is similar for all spirals, and significantly lower in early type galaxies (ETG), in the past $\Sigma_{\rm SFR}$ scales well with morphology. The central regions of today's ETGs are where $\Sigma_{\rm SFR}$ reached the highest values ($> 10^3 \,M_\odot\,$Gyr$^{-1}\,$pc$^{-2}$), similar to those measured in high redshift star forming galaxies. d) The evolution of $\Sigma_{\rm SFR}$ in Sbc systems matches that of models for Milky-Way-like galaxies, suggesting that the formation of a thick disk may be a common phase in spirals at early epochs.

We report on a detailed study of the stellar populations and ionized gas properties in the merger LIRG NGC 2623, analysing optical Integral Field Spectroscopy from the CALIFA survey and PMAS LArr, multiwavelength HST imaging, and OSIRIS narrow band H$\alpha$ and [NII]$\lambda$6584 imaging. The spectra were processed with the STARLIGHT full spectral fitting code, and the results compared with those for two early-stage merger LIRGs (IC 1623 W and NGC 6090), together with CALIFA Sbc/Sc galaxies. We find that NGC 2623 went through two periods of increased star formation (SF), a first and widespread episode, traced by intermediate-age stellar populations ISP (140 Myr-1.4 Gyr), and a second one, traced by young stellar populations YSP ($<$140 Myr), which is concentrated in the central regions ($<$1.4 kpc). Our results are in agreement with the epochs of the first peri-center passage ($\sim$200 Myr ago) and coalescence ($<$100 Myr ago) predicted by dynamical models, and with high resolution merger simulations in the literature, consistent with NGC 2623 representing an evolved version of the early-stage mergers. Most ionized gas is concentrated within $<$2.8 kpc, where LINER-like ionization and high velocity dispersion ($\sim$220 km/s) are found, consistent with the previously reported outflow. As revealed by the highest resolution OSIRIS and HST data, a collection of HII regions is also present in the plane of the galaxy, which explains the mixture of ionization mechanisms in this system. It is unlikely that the outflow in NGC 2623 will escape from the galaxy, given the low SFR intensity ($\sim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), the fact that the outflow rate is 3 times lower than the current SFR, and the escape velocity in the central areas higher than the outflow velocity.

The role of major mergers in galaxy evolution is investigated through a detailed characterization of the stellar populations, ionized gas properties, and star formation rates (SFR) in the early-stage merger LIRGs IC 1623 W and NGC 6090, by analysing optical Integral Field Spectroscopy (IFS) and high resolution HST imaging. The spectra were processed with the Starlight full spectral fitting code, and the emission lines measured in the residual spectra. The results are compared with control non-interacting spiral galaxies from the CALIFA survey. Merger-induced star formation is extended and recent, as revealed by the young ages (50-80 Myr) and high contributions to light of young stellar populations (50-90$\%$), in agreement with merger simulations in the literature. These early-stage mergers have positive central gradients of the stellar metallicity, with an average $\sim$0.6 Z$_{\odot}$. Compared to non-interacting spirals, they have lower central nebular metallicity, and flatter profiles, in agreement with the gas inflow scenario. We find that they are dominated by star formation, although shock excitation cannot be discarded in some regions, where high velocity dispersion is found (170-200 km s$^{-1}$). The average SFR in these early-stage mergers ($\sim$23-32 M$_{\odot}$ yr$^{-1}$) is enhanced with respect to main-sequence Sbc galaxies by factors of 6-9, slightly above the predictions from classical merger simulations, but still possible in about 15$\%$ of major galaxy mergers, where U/LIRGs belong.

We consider $r$-process nucleosynthesis in outflows from black hole accretion discs formed in double neutron star and neutron star -- black hole mergers. These outflows, powered by angular momentum transport processes and nuclear recombination, represent an important -- and in some cases dominant -- contribution to the total mass ejected by the merger. Here we calculate the nucleosynthesis yields from disc outflows using thermodynamic trajectories from hydrodynamic simulations, coupled to a nuclear reaction network. We find that outflows produce a robust abundance pattern around the second $r$-process peak (mass number $A \sim 130$), independent of model parameters, with significant production of $A < 130$ nuclei. This implies that dynamical ejecta with high electron fraction may not be required to explain the observed abundances of $r$-process elements in metal poor stars. Disc outflows reach the third peak ($ A \sim 195$) in most of our simulations, although the amounts produced depend sensitively on the disc viscosity, initial mass or entropy of the torus, and nuclear physics inputs. Some of our models produce an abundance spike at $A = 132$ that is absent in the Solar system $r$-process distribution. The spike arises from convection in the disc and depends on the treatment of nuclear heating in the simulations. We conclude that disc outflows provide an important -- and perhaps dominant -- contribution to the $r$-process yields of compact binary mergers, and hence must be included when assessing the contribution of these systems to the inventory of $r$-process elements in the Galaxy.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

The aim of this paper is to characterize the radial structure of the star formation rate (SFR) in galaxies in the nearby Universe as represented by the CALIFA survey. The sample under study contains 416 galaxies observed with IFS, covering a wide range of Hubble types and stellar masses. Spectral synthesis techniques are applied to obtain radial profiles of the intensity of the star formation rate in the recent past, and the local sSFR. To emphasize the behavior of these properties for galaxies that are on and off the main sequence of star formation (MSSF) we stack the individual radial profiles in bins of galaxy morphology and stellar masses. Our main results are: a) The intensity of SFR shows declining profiles that exhibit very little differences between spirals. The dispersion between the profiles is significantly smaller in late type spirals. This confirms that the MSSF is a sequence of galaxies with nearly constant intensity of SFR b) sSFR values scale with Hubble type and increase radially outwards, with a steeper slope in the inner 1 HLR. This behavior suggests that galaxies are quenched inside-out, and that this process is faster in the central, bulge-dominated part than in the disks. c) As a whole, and at all radii, E and S0 are off the MSSF. d) Applying the volume-corrections for the CALIFA sample, we obtain a density of star formation in the local Universe of 0.0105 Msun/yr/Mpc^{-3}. Most of the star formation is occurring in the disks of spirals. e) The volume averaged birthrate parameter, b'=0.39, suggests that the present day Universe is forming stars at 1/3 of its past average rate. E, S0, and the bulge of early type spirals contribute little to the recent SFR of the Universe, which is dominated by the disks of later spirals. f) There is a tight relation between the intensity of the SFR and stellar mass, defining a local MSSF relation with a logarithmic slope of 0.8.

We present an extended version of the spectral synthesis code STARLIGHT designed to incorporate both $\lambda$-by-$\lambda$ spectra and photometric fluxes in the estimation of stellar population properties of galaxies. The code is tested with simulations and data for 260 galaxies culled from the CALIFA survey, spatially matching the 3700--7000 \AA\ optical datacubes to GALEX near and far UV images. The sample spans E--Sd galaxies with masses from $10^9$ to $10^{12} M_\odot$ and stellar populations all the way from star-forming to old, passive systems. Comparing results derived from purely optical fits with those which also consider the NUV and FUV data we find that: (1) The new code is capable of matching the input UV data within the errors while keeping the quality of the optical fit essentially unchanged. (2) Despite being unreliable predictors of the UV fluxes, purely optical fits yield stellar population properties which agree well with those obtained in optical+UV fits for nearly 90% of our sample. (3) The addition of UV constraints has little impact on properties such as stellar mass and dust optical depth. Mean stellar ages and metallicities also remain nearly the same for most galaxies, the exception being low-mass, late-type galaxies, which become older and less enriched due to rearrangements of their youngest populations. (4) The revised ages are better correlated with observables such as the 4000 \AA\ break index, and the $NUV - r$ and $u - r$ colours, an empirical indication that the addition of UV constraints helps mitigating the effects of age-metallicity-extinction degeneracies.

The mergers of binaries containing neutron stars and stellar-mass black holes are the most promising sources for direct detection in gravitational waves by the interferometers Advanced LIGO and Virgo over the next few years. The concurrent detection of electromagnetic emission from these events would greatly enhance the scientific return of these discoveries. Here we review the state of the art in modeling the electromagnetic signal of neutron star binary mergers across different phases of the merger and multiple wavelengths. We focus on those observables which provide the most sensitive diagnostics of the merger physics and the contribution to the synthesis of rapid neutron capture ($r$-process) elements in the Galaxy. We also outline expected future developments on the observational and theoretical sides of this rapidly evolving field.

We investigate the effect of dimensionality on the transition to explosion in neutrino-driven core-collapse supernovae. Using parameterized hydrodynamic simulations of the stalled supernova shock in one-, two- (2D), and three spatial dimensions (3D), we systematically probe the extent to which hydrodynamic instabilities alone can tip the balance in favor of explosion. In particular, we focus on systems that are well into the regimes where the Standing Accretion Shock Instability (SASI) or neutrino-driven convection dominate the dynamics, and characterize the difference between them. We find that SASI-dominated models can explode with up to ~20% lower neutrino luminosity in 3D than in 2D, with the magnitude of this difference decreasing with increasing resolution. This improvement in explosion conditions is related to the ability of spiral modes to generate more non-radial kinetic energy than a single sloshing mode, increasing the size of the average shock radius, and hence generating better conditions for the formation of large-scale, high-entropy bubbles. In contrast, convection-dominated explosions show a smaller difference in their critical heating rate between 2D and 3D (<8%), in agreement with previous studies. The ability of our numerical implementation to maintain arbitrary symmetries is quantified with a set of SASI-based tests. We discuss implications for the diversity of explosion paths in a realistic supernova environment.

This paper characterizes the radial structure of stellar population properties of galaxies in the nearby universe, based on 300 galaxies from the CALIFA survey. The sample covers a wide range of Hubble types, and galaxy stellar mass. We apply the spectral synthesis techniques to recover the stellar mass surface density, stellar extinction, light and mass-weighted ages, and mass-weighted metallicity, for each spatial resolution element in our target galaxies. To study mean trends with overall galaxy properties, the individual radial profiles are stacked in seven bins of galaxy morphology. We confirm that more massive galaxies are more compact, older, more metal rich, and less reddened by dust. Additionally, we find that these trends are preserved spatially with the radial distance to the nucleus. Deviations from these relations appear correlated with Hubble type: earlier types are more compact, older, and more metal rich for a given mass, which evidences that quenching is related to morphology, but not driven by mass. Negative gradients of ages are consistent with an inside-out growth of galaxies, with the largest ages gradients in Sb-Sbc galaxies. Further, the mean stellar ages of disks and bulges are correlated, with disks covering a wider range of ages, and late type spirals hosting younger disks. The gradients in stellar mass surface density depend mostly on stellar mass, in the sense that more massive galaxies are more centrally concentrated. There is a secondary correlation in the sense that at the same mass early type galaxies have steeper gradients. We find mildly negative metallicity gradients, shallower than predicted from models of galaxy evolution in isolation. The largest gradients occur in Sb galaxies. Overall we conclude that quenching processes act in manners that are independent of mass, while metallicity and galaxy structure are influenced by mass-dependent processes.

We characterize in detail the radial structure of the stellar population properties of 300 galaxies in the nearby universe, observed with integral field spectroscopy in the CALIFA survey. The sample covers a wide range of Hubble types, from spheroidal to spiral galaxies, ranging in stellar masses from $M_\star \sim 10^9$ to $7 \times 10^{11}$ $M_\odot$. We derive the stellar mass surface density ($\mu_\star$), light-weighted and mass-weighted ages ($\langle {\rm log}\,age\rangle _L$, $\langle {\rm log}\,age\rangle _M$), and mass-weighted metallicity ($\langle {\rm log}\,Z_\star\rangle _M$), applying the spectral synthesis technique. We study the mean trends with galaxy stellar mass, $M_\star$, and morphology (E, S0, Sa, Sb, Sbc, Sc and Sd). We confirm that more massive galaxies are more compact, older, more metal rich, and less reddened by dust. Additionally, we find that these trends are preserved spatially with the radial distance to the nucleus. Deviations from these relations appear correlated with Hubble type: earlier types are more compact, older, and more metal rich for a given M$_\star$, which evidences that quenching is related to morphology, but not driven by mass.

We explore the evolution of the different ejecta components generated during the merger of a neutron star (NS) and a black hole (BH). Our focus is the interplay between material ejected dynamically during the merger, and the wind launched on a viscous timescale by the remnant accretion disk. These components are expected to contribute to an electromagnetic transient and to produce r-process elements, each with a different signature when considered separately. Here we introduce a two-step approach to investigate their combined evolution, using two- and three-dimensional hydrodynamic simulations. Starting from the output of a merger simulation, we identify each component in the initial condition based on its phase space distribution, and evolve the accretion disk in axisymmetry. The wind blown from this disk is injected into a three-dimensional computational domain where the dynamical ejecta is evolved. We find that the wind can suppress fallback accretion on timescales longer than ~100 ms. Due to self-similar viscous evolution, the disk accretion at late times nevertheless approaches a power-law time dependence $\propto t^{-2.2}$. This can power some late-time GRB engine activity, although the available energy is significantly less than in traditional fallback models. Inclusion of radioactive heating due to the r-process does not significantly affect the fallback accretion rate or the disk wind. We do not find any significant modification to the wind properties at large radius due to interaction with the dynamical ejecta. This is a consequence of the different expansion velocities of the two components.

We analyze the spatially resolved star formation history of 300 nearby galaxies from the CALIFA integral field spectroscopic survey to investigate the radial structure and gradients of the present day stellar populations properties as a function of Hubble type and galaxy stellar mass. A fossil record method based on spectral synthesis techniques is used to recover spatially and temporally resolved maps of stellar population properties of spheroidal and spiral galaxies with masses $10^9$ to $7 \times 10^{11}$ M$_\odot$. The results show that galaxy-wide spatially averaged stellar population properties (stellar mass, mass surface density, age, metallicity, and extinction) match those obtained from the integrated spectrum, and that these spatially averaged properties match those at $R = 1$ HLR (half light radius), proving that the effective radii are really effective. Further, the individual radial profiles of the stellar mass surface density ($\mu_\star$), luminosity weighted ages ($< {\rm log}\,age>_L$), and mass weighted metallicity ($< \log Z_\star >_M$) are stacked in bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc and Sd). All these properties show negative gradients as a sign of the inside-out growth of massive galaxies. However, the gradients depend on the Hubble type in different ways. For the same galaxy mass, E and S0 galaxies show the largest inner gradients in $\mu_\star$; while MW-like galaxies (Sb with $M_\star \sim 10^{11} M_\odot$) show the largest inner age and metallicity gradients. The age and metallicity gradients suggest that major mergers have a relevant role in growing the center (within 3 HLR) of massive early type galaxies; and radial mixing may play a role flattening the radial metallicity gradient in MW-like disks.

The accretion disk that forms after a neutron star merger is a source of neutron-rich ejecta. The ejected material contributes to a radioactively-powered electromagnetic transient, with properties that depend sensitively on the composition of the outflow. Here we investigate how the spin of the black hole remnant influences mass ejection on the thermal and viscous timescales. We carry out two-dimensional, time-dependent hydrodynamic simulations of merger remnant accretion disks including viscous angular momentum transport and approximate neutrino self-irradiation. The gravity of the spinning black hole is included via a pseudo-Newtonian potential. We find that a disk around a spinning black hole ejects more mass, up to a factor of several, relative to the non-spinning case. The enhanced mass loss is due to energy release by accretion occurring deeper in the gravitational potential, raising the disk temperature and hence the rate of viscous heating in regions where neutrino cooling is ineffective. The mean electron fraction of the outflow increases moderately with BH spin due to a highly-irradiated (though not neutrino-driven) wind component. While the bulk of the ejecta is still very neutron-rich, thus generating heavy r-process elements, the leading edge of the wind contains a small amount of Lanthanide-free material. This component can give rise to a ~1 day blue optical `bump' in a kilonova light curve, even in the case of prompt BH formation, which may facilitate its detection.

Methods to recover the fossil record of galaxy evolution encoded in their optical spectra have been instrumental in processing the avalanche of data from mega-surveys along the last decade, effectively transforming observed spectra onto a long and rich list of physical properties: from stellar masses and mean ages to full star formation histories. This promoted progress in our understanding of galaxies as a whole. Yet, the lack of spatial resolution introduces undesirable aperture effects, and hampers advances on the internal physics of galaxies. This is now changing with 3D surveys. The mapping of stellar populations in data-cubes allows us to figure what comes from where, unscrambling information previously available only in integrated form. This contribution uses our starlight-based analysis of 300 CALIFA galaxies to illustrate the power of spectral synthesis applied to data-cubes. The selected results highlighted here include: (a) The evolution of the mass-metallicity and mass-density-metallicity relations, as traced by the mean stellar metallicity. (b) A comparison of star formation rates obtained from H{\alpha} to those derived from full spectral fits. (c) The relation between star formation rate and dust optical depth within galaxies, which turns out to mimic the Schmidt-Kennicutt law. (d) PCA tomography experiments.

We use spatially and temporally resolved maps of stellar population properties of 300 galaxies from the CALIFA integral field survey to investigate how the stellar metallicity (Z*) relates to the total stellar mass (M*) and the local mass surface density ($\mu$*) in both spheroidal and disk dominated galaxies. The galaxies are shown to follow a clear stellar mass-metallicity relation (MZR) over the whole 10$^9$ to 10$^{12}$ M$_{\odot}$ range. This relation is steeper than the one derived from nebular abundances, which is similar to the flatter stellar MZR derived when we consider only young stars. We also find a strong relation between the local values of $\mu$* and Z* (the $\mu$ZR), betraying the influence of local factors in determining Z*. This shows that both local ($\mu$*-driven) and global (M*-driven) processes are important in determining the metallicity in galaxies. We find that the overall balance between local and global effects varies with the location within a galaxy. In disks, $\mu$* regulates Z*, producing a strong $\mu$ZR whose amplitude is modulated by M*. In spheroids it is M* who dominates the physics of star formation and chemical enrichment, with $\mu$* playing a minor, secondary role. These findings agree with our previous analysis of the star formation histories of CALIFA galaxies, which showed that mean stellar ages are mainly governed by surface density in galaxy disks and by total mass in spheroids.

Mergers of binary neutron stars (NSs) usually result in the formation of a hypermassive neutron star (HMNS). Whether- and when this remnant collapses to a black hole (BH) depends primarily on the equation of state and on angular momentum transport processes, both of which are uncertain. Here we show that the lifetime of the merger remnant may be directly imprinted in the radioactively powered kilonova emission following the merger. We employ axisymmetric, time-dependent hydrodynamic simulations of remnant accretion disks orbiting a HMNS of variable lifetime, and characterize the effect of this delay to BH formation on the disk wind ejecta. When BH formation is relatively prompt (~ 100 ms), outflows from the disk are sufficiently neutron rich to form heavy r-process elements, resulting in ~ week-long emission with a spectral peak in the near-infrared (NIR), similar to that produced by the dynamical ejecta. In contrast, delayed BH formation allows neutrinos from the HMNS to raise the electron fraction in the polar direction to values such that potentially Lanthanide-free outflows are generated. The lower opacity would produce a brighter, bluer, and shorter-lived ~ day-long emission (a `blue bump') prior to the late NIR peak from the dynamical ejecta and equatorial wind. This new diagnostic of BH formation should be useful for events with a signal to noise lower than that required for direct detection of gravitational waveform signatures.

In the collapsing core of massive stars, the standing accretion shock instability (SASI) can drive spiral modes that efficiently redistribute angular momentum. This process can impart a spin to the forming neutron star even when the progenitor star is non-rotating. Here we develop the first analytical description of the angular momentum redistribution driven by a spiral mode of the SASI. Our analysis, valid in the limit of small mode amplitude, shows that the angular momentum separation is driven by the Reynolds stress generated by the spiral mode. The resulting solutions compare favorably with previous three-dimensional hydrodynamic simulations of the SASI in the linear and weakly non-linear phases. Reasonable agreement is also found when extrapolating the solutions into the fully non-linear phase. A Reynolds-decomposition of the flow is performed in the saturated state of these simulations, showing that outward angular momentum transport by the Reynolds stress and the fluctuating component of the mass flux balance inward transport by advection. We derive an approximate analytic expression for the maximum angular momentum deposited in the neutron star as a function of the mass accretion rate, shock radius, shock compression ratio, and amplitude of the spiral mode at the time of explosion. Implications for the birth spin periods of neutron stars are discussed.

The success of the neutrino mechanism of core-collapse supernovae relies on the supporting action of two hydrodynamic instabilities: neutrino-driven convection and the Standing Accretion Shock Instability (SASI). Depending on the structure of the stellar progenitor, each of these instabilities can dominate the evolution of the gain region prior to the onset of explosion, with implications for the ensuing asymmetries. Here we examine the flow dynamics in the neighborhood of explosion by means of parametric two-dimensional, time-dependent hydrodynamic simulations for which the linear stability properties are well understood. We find that systems for which the convection parameter is sub-critical (SASI-dominated) develop explosions once large-scale, high-entropy bubbles are able to survive for several SASI oscillation cycles. These long-lived structures are seeded by the SASI during shock expansions. Finite-amplitude initial perturbations do not alter this outcome qualitatively, though they can lead to significant differences in explosion times. Supercritical systems (convection-dominated) also explode by developing large-scale bubbles, though the formation of these structures is due to buoyant activity. Non-exploding systems achieve a quasi-steady state in which the time-averaged flow adjusts itself to be convectively sub-critical. We characterize the turbulent flow using a spherical Fourier-Bessel decomposition, identifying the relevant scalings and connecting temporal and spatial components. Finally, we verify the applicability of these principles on the general relativistic, radiation-hydrodynamic simulations of Mueller, Janka, & Heger (2012), and discuss implications for the three-dimensional case.

Type Ia supernovae (SNe Ia), thermonuclear explosions of carbon-oxygen white dwarfs (CO-WDs), are currently the best cosmological "standard candles", but the triggering mechanism of the explosion is unknown. It was recently shown that the rate of head-on collisions of typical field CO-WDs in triple systems may be comparable to the SNe Ia rate. Here we provide evidence supporting a scenario in which the majority of SNe Ia are the result of such head-on collisions of CO-WDs. In this case, the nuclear detonation is due to a well understood shock ignition, devoid of commonly introduced free parameters such as the deflagration velocity or transition to detonation criteria. By using two-dimensional hydrodynamical simulations with a fully resolved ignition process, we show that zero-impact-parameter collisions of typical CO-WDs with masses $0.5-1\,M_{\odot}$ result in explosions that synthesize $^{56}$Ni masses in the range of $\sim0.1-1\,M_{\odot}$, spanning the wide distribution of yields observed for the majority of SNe Ia. All collision models yield the same late-time ($>60$ days since explosion) bolometric light curve when normalized by $^{56}$Ni masses (to better than $30\%$), in agreement with observations. The calculated widths of the $^{56}$Ni-mass-weighted-line-of-sight velocity distributions are correlated with the calculated $^{56}$Ni yield, agreeing with the observed correlation. The strong correlation, shown here for the first time, between $^{56}$Ni yield and total mass of the colliding CO-WDs (insensitive to their mass ratio), is suggestive as the source for the continuous distribution of observed SN Ia features, possibly including the Philips relation.

We study the radial structure of the stellar mass surface density ($\mu$) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We use the fossil record to recover the star formation history (SFH) in spheroidal and disk dominated galaxies with masses from 10$^9$ to 10$^{12}$ M$_\odot$. We derive the half mass radius, and we find that galaxies are on average 15% more compact in mass than in light. HMR/HLR decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. We find that the galaxy-averaged stellar population age, stellar extinction, and $\mu$ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages support an inside-out formation. The larger inner age gradients occur in the most massive disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative $\mu$ gradients that steepen with stellar mass. In spheroidal galaxies $\mu$ saturates at a critical value that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local $\mu$ at the same radius (in HLR units). The SFH of the regions beyond 1 HLR are well correlated with their local $\mu$, and follow the same relation as the galaxy-averaged age and $\mu$; suggesting that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in $\mu$. Thus, galaxy mass is a more fundamental property in spheroidal systems while the local stellar mass surface density is more important in disks.

Expulsion of neutron-rich matter following the merger of neutron star (NS) binaries is crucial to the radioactively-powered electromagnetic counterparts of these events and to their relevance as sources of r-process nucleosynthesis. Here we explore the long-term (viscous) evolution of remnant black hole accretion disks formed in such mergers by means of two-dimensional, time-dependent hydrodynamical simulations. The evolution of the electron fraction due to charged-current weak interactions is included, and neutrino self-irradiation is modeled as a lightbulb that accounts for the disk geometry and moderate optical depth effects. Over several viscous times (~1s), a fraction ~10% of the initial disk mass is ejected as a moderately neutron-rich wind (Y_e ~ 0.2) powered by viscous heating and nuclear recombination, with neutrino self-irradiation playing a sub-dominant role. Although the properties of the outflow vary in time and direction, their mean values in the heavy-element production region are relatively robust to variations in the initial conditions of the disk and the magnitude of its viscosity. The outflow is sufficiently neutron-rich that most of the ejecta forms heavy r-process elements with mass number A >130, thus representing a new astrophysical source of r-process nucleosynthesis, distinct from that produced in the dynamical ejecta. Due to its moderately high entropy, disk outflows contain a small residual fraction ~1% of helium, which could produce a unique spectroscopic signature.

We examine the nonlinear development of unstable magnetosonic waves driven by a background radiative flux -- the Radiation-Driven Magneto-Acoustic Instability (RMI, a.k.a. the "photon bubble" instability). The RMI may serve as a persistent source of density, radiative flux, and magnetic field fluctuations in stably-stratified, optically-thick media. The conditions for instability are present in a variety of astrophysical environments, and do not require the radiation pressure to dominate or the magnetic field to be strong. Here we numerically study the saturation properties of the RMI, covering three orders of magnitude in the relative strength of radiation, magnetic field, and gas energies. Two-dimensional, time-dependent radiation-MHD simulations of local, stably-stratified domains are conducted with Zeus-MP in the optically-thick, highly-conducting limit. Our results confirm the theoretical expectations of Blaes and Socrates (2003) in that the RMI operates even in gas pressure-dominated environments that are weakly magnetized. The saturation amplitude is a monotonically increasing function of the ratio of radiation to gas pressure. Keeping this ratio constant, we find that the saturation amplitude peaks when the magnetic pressure is comparable to the radiation pressure. We discuss the implications of our results for the dynamics of magnetized stellar envelopes, where the RMI should act as a source of sub-photospheric perturbations.

The Calar Alto Legacy Integral Field Area (CALIFA) is an ongoing 3D spectroscopic survey of 600 nearby galaxies of all kinds. This pioneer survey is providing valuable clues on how galaxies form and evolve. Processed through spectral synthesis techniques, CALIFA datacubes allow us to, for the first time, spatially resolve the star formation history of galaxies spread across the color-magnitude diagram. The richness of this approach is already evident from the results obtained for the first 107 galaxies. Here we show how the different galactic spatial sub-components ("bulge" and "disk") grow their stellar mass over time. We explore the results stacking galaxies in mass bins, finding that, except at the lowest masses, galaxies grow inside-out, and that the growth rate depends on a galaxy's mass. The growth rate of inner and outer regions differ maximally at intermediate masses. We also find a good correlation between the age radial gradient and the stellar mass density, suggesting that the local density is a main driver of galaxy evolution.

We explore the evolution of radiatively inefficient accretion disks in which nuclear reactions are dynamically important (`Nuclear Dominated Accretion Flows', or NuDAFs). Examples of such disks are those generated by the merger of a white dwarf with a neutron star or black hole, or by the collapse of a rotating star. Here we present two-dimensional hydrodynamic simulations that systematically explore the effect of adding a single nuclear reaction to a viscous torus. The equation of state, anomalous shear stress, and nuclear reactions are given parametric forms. Our results point to the existence of two qualitatively different regimes of NuDAF evolution: (1) steady accretion with quiescent burning; or (2) detonation of the disk. These outcomes are controlled primarily by the ratio of the nuclear energy released to the enthalpy at the burning radius. Disks detonate if this ratio exceeds a critical value (~1), and if burning occurs in regions where neutrino cooling is unimportant. Thermonuclear runaways are seeded by the turbulent mixing of hot ash with cold fuel at the burning front. Subcritical disks do not explode, but instead power a persistent collimated outflow of unbound material composed primarily of ash, with a mass-loss rate that increases with the nuclear binding energy released. We discuss the implications of our results for supernova-like counterparts from astrophysical events in the NuDAF regime. In particular, detonations following a white dwarf - neutron star merger could account for some subluminous Type Ia supernovae, such as the class defined by SN 2002cx.

We study the transition to runaway expansion of an initially stalled core-collapse supernova shock. The neutrino luminosity, mass accretion rate, and neutrinospheric radius are all treated as free parameters. In spherical symmetry, this transition is mediated by a global non-adiabatic instability that develops on the advection time and reaches non-linear amplitude. Here we perform high-resolution, time-dependent hydrodynamic simulations of stalled supernova shocks with realistic microphysics to analyze this transition. We find that radial instability is a sufficient condition for runaway expansion if the neutrinospheric parameters do not vary with time and if heating by the accretion luminosity is neglected. For a given unstable mode, transition to runaway occurs when fluid in the gain region reaches positive specific energy. We find approximate instability criteria that accurately describe the behavior of the system over a wide region of parameter space. The threshold neutrino luminosities are in general different than the limiting value for a steady-state solution. We hypothesize that multidimensional explosions arise from the excitation of unstable large-scale modes of the turbulent background flow, at threshold luminosities that are lower than in the laminar case.

In the magnetar model, the quiescent non-thermal soft X-ray emission from Anomalous X-ray Pulsars and Soft-Gamma Repeaters is thought to arise from resonant comptonization of thermal photons by charges moving in a twisted magnetosphere. Robust inference of physical quantities from observations is difficult, because the process depends strongly on geometry and current understanding of the magnetosphere is not very deep. The polarization of soft X-ray photons is an independent source of information, and its magnetospheric imprint remains only partially explored. In this paper we calculate how resonant cyclotron scattering would modify the observed polarization signal relative to the surface emission, using a multidimensional Monte Carlo radiative transfer code that accounts for the gradual coupling of polarization eigenmodes as photons leave the magnetosphere. We employ a globally-twisted, self-similar, force-free magnetosphere with a power-law momentum distribution, assume a blackbody spectrum for the seed photons, account for general relativistic light deflection close to the star, and assume that vacuum polarization dominates the dielectric properties of the magnetosphere. The latter is a good approximation if the pair multiplicity is not much larger than unity. Phase-averaged polarimetry is able to provide a clear signature of the magnetospheric reprocessing of thermal photons and to constrain mechanisms generating the thermal emission. Phase-resolved polarimetry, in addition, can characterize the spatial extent and magnitude of the magnetospheric twist angle at ~100 stellar radii, and discern between uni- or bidirectional particle energy distributions, almost independently of every other parameter in the system. We discuss prospects for detectability with GEMS.

A stalled spherical accretion shock, such as that arising in core-collapse supernovae, is unstable to non-spherical perturbations. In three dimensions, this Standing Accretion Shock Instability (SASI) can develop spiral modes that spin-up the protoneutron star. Here we study these non-axisymmetric modes by combining linear stability analysis and three-dimensional, time-dependent hydrodynamic simulations with Zeus-MP, focusing on characterizing their spatial structure and angular momentum content. We do not impose any rotation on the background accretion flow, and use simplified microphysics with no neutrino heating or nuclear dissociation. Spiral modes are examined in isolation by choosing flow parameters such that only the fundamental mode is unstable for a given polar index, leading to good agreement with linear theory. We find that any superposition of sloshing modes with non-zero relative phases survives in the nonlinear regime and leads to angular momentum redistribution. It follows that the range of perturbations required to obtain spin-up is broader than that needed to obtain the limiting case of a phase shift of pi/2. The bulk of the angular momentum redistribution occurs during a phase of exponential growth, and arises from internal torques that are second order in the perturbation amplitude. This redistribution gives rise to at least two counter rotating regions, with the maximum angular momentum of a given sign approaching a significant fraction of the mass accretion rate times the shock radius squared (Mdot * r_shock^2 ~ 1E+47 g/cm^2/s, spin period ~60 ms). Nonlinear mode coupling at saturation causes the angular momentum to fluctuate in all directions with much smaller amplitudes.

A hypothetical time-variation of the gravitational constant $G$ would make neutron stars expand or contract, so the matter in their interiors would depart from beta equilibrium. This induces non-equilibrium weak reactions, which release energy that is invested partly in neutrino emission and partly in internal heating. Eventually, the star arrives at a stationary state in which the temperature remains nearly constant, as the forcing through the change of $G$ is balanced by the ongoing reactions. Using the surface temperature of the nearest millisecond pulsar (PSR J0437$-$4715) inferred from ultraviolet observations and results from theoretical modelling of the thermal evolution, we estimate two upper limits for this variation: (1) $|\dot G/G| < 2 \times 10^{-10}\mathrm{yr}^{-1},$ if the fast, "direct Urca" reactions are allowed, and (2) $|\dot G/G|<4\times 10^{-12}\mathrm{yr}^{-1},$ considering only the slower, "modified Urca" reactions. The latter is among the most restrictive upper limits obtained by other methods.

We investigate the effects of neutrino heating and alpha-particle recombination on the hydrodynamics of core-collapse supernovae. Our focus is on the non-linear dynamics of the shock wave that forms in the collapse, and the assembly of positive energy material below it. To this end, we perform time-dependent hydrodynamic simulations with FLASH2.5 in spherical and axial symmetry. These generalize our previous calculations by allowing for bulk neutrino heating and for nuclear statistical equilibrium between n, p and alpha. The heating rate is freely tunable, as is the starting radius of the shock relative to the recombination radius of alpha-particles. An explosion in spherical symmetry involves the excitation of an overstable mode, which may be viewed as the L=0 version of the `Standing Accretion Shock Instability'. In 2D simulations, non-spherical deformations of the shock are driven by plumes of material with positive Bernoulli parameter, which are concentrated well outside the zone of strong neutrino heating. The non-spherical modes of the shock reach a large amplitude only when the heating rate is also high enough to excite convection below the shock. The critical heating rate that causes an explosion depends sensitively on the initial position of the shock relative to the recombination radius. Weaker heating is required to drive an explosion in 2D than in 1D, but the difference also depends on the size of the shock. Forcing the infalling heavy nuclei to break up into n and p below the shock only causes a slight increase in the critical heating rate, except when the shock starts out at a large radius. This shows that heating by neutrinos (or some other mechanism) must play a significant role in pushing the shock far enough out that recombination heating takes over.

We examine the stability of a standing shock wave within a spherical accretion flow onto a gravitating star, in the context of core-collapse supernova explosions. Our focus is on the effect of nuclear dissociation below the shock on the linear growth, and non-linear saturation, of non-radial oscillations of the shocked fluid. We combine two-dimensional, time-dependent hydrodynamic simulations using FLASH2.5 with a solution to the linear eigenvalue problem, and demonstrate the consistency of the two approaches. Previous studies of this `Standing Accretion Shock Instability' (SASI) have focused either on zero-energy accretion flows without nuclear dissociation, or made use of a detailed finite-temperature nuclear equation of state and included strong neutrino heating. Our main goal in this and subsequent papers is to introduce equations of state of increasing complexity, in order to isolate the various competing effects. In this work we employ an ideal gas equation of state with a constant rate of nuclear dissociation below the shock, and do not include neutrino heating. We find that a negative Bernoulli parameter below the shock significantly lowers the real frequency, growth rate, and saturation amplitude of the SASI. A decrease in the adiabatic index has similar effects. The non-linear development of the instability is characterized by an expansion of the shock driven by turbulent kinetic energy at nearly constant internal energy. Our results also provide further insight into the instability mechanism: the rate of growth of a particular mode is fastest when the radial advection time from the shock to the accretor overlaps with the period of a standing lateral sound wave. The fastest-growing mode can therefore be modified by nuclear dissociation.

(Abridged) The main sequence star beta Pictoris hosts the best studied circumstellar disk to date. Nonetheless, a long-standing puzzle has been around since the detection of metallic gas in the disk: radiation pressure from the star should blow the gas away, yet the observed motion is consistent with Keplerian rotation. In this work we search for braking mechanisms that can resolve this discrepancy. We find that all species affected by radiation force are heavily ionized and dynamically coupled into a single fluid by Coulomb collisions, reducing the radiation force on species feeling the strongest acceleration. For a gas of solar composition, the resulting total radiation force still exceeds gravity, while a gas of enhanced carbon abundance could be self-braking. We also explore two other braking agents: collisions with dust grains and neutral gas. Grains surrounding beta Pic are photoelectrically charged to a positive electrostatic potential. If a significant fraction of the grains are carbonaceous (10% in the midplane and larger at higher altitudes), ions can be slowed down to satisfy the observed velocity constraints. For neutral gas to brake the coupled ion fluid, we find the minimum required mass to be $\approx$ 0.03 $M_\earth$, consistent with observed upper limits of the hydrogen column density, and substantially reduced relative to previous estimates. Our results favor a scenario in which metallic gas is generated by grain evaporation in the disk, perhaps during grain-grain collisions. We exclude a primordial origin for the gas, but cannot rule out the possibility of its production by falling evaporating bodies near the star. We discuss the implications of this work for observations of gas in other debris disks.

Rotochemical heating originates in a departure from beta equilibrium due to spin-down compression in a rotating neutron star. The main consequence is that the star eventually arrives at a quasi-equilibrium state, in which the thermal photon luminosity depends only on the current value of the spin-down power, which is directly measurable. Only in millisecond pulsars the spin-down power remains high long enough for this state to be reached with a substantial luminosity. We report an extensive study of the effect of this heating mechanism on the thermal evolution of millisecond pulsars, developing a general formalism in the slow-rotation approximation of general relativity that takes the spatial structure of the star fully into account, and using a sample of realistic equations of state to solve the non-superfluid case numerically. We show that nearly all observed millisecond pulsars are very likely to be in the quasi-equilibrium state. Our predicted quasi-equilibrium temperatures for PSR J0437-4715 are only 20% lower than inferred from observations. Accounting for superfluidity should increase the predicted value.