We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood.

We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production.

We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($\eta < 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $\eta > 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa.

We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallicity $Z=0.01$, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including $^{51}$V and $^{52}$Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 $M_\odot$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 $M_\odot$ model at $Z=0.01$ in the yields interpolation scheme of our GCE model for stars in between 15 and 25 $M_\odot$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] $\approx$ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the pre-supernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si-shell could launch a convective-reactive global oscillation, if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,\gamma)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,\gamma)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $\gamma$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects.

(Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates.

We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2\sigma$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3\sigma$ ($\rm5\sigma$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN.

We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and significantly blue-shifted. The derived M_BH from the blueshift-corrected broad CIV line width agrees with the MgII result. From direct measurement using a well-sampled SED, the bolometric luminosity is 3.6 * 10^14 L_sun. The corresponding Eddington ratio for W2246-0526 is lambda_Edd = L_AGN / L_Edd = 2.8. This high Eddington ratio may reach the level where the luminosity is saturating due to photon trapping in the accretion flow, and be insensitive to the mass accretion rate. In this case, the M_BH growth rate in W2246-0526 would exceed the apparent accretion rate derived from the observed luminosity.

We provide here a significant extension of the NuGrid Set 1 models in mass coverage and toward lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M/Msun = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with alpha-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass loss, convective boundary mixing in low- and intermediate mass models and H and He core burning massive star models is included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K and Sc. In AGB models with hot dredge-up the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity H-ingestion events are observed and lead to i-process nucleosynthesis and substantial N-15 production. Complete yield data tables, derived data products and online analytic data access are provided.

The first peak s-process elements Rb, Sr, Y and Zr in the post-AGB star Sakurai's object (V4334 Sagittarii) have been proposed to be the result of i-process nucleosynthesis in a post-AGB very-late thermal pulse event. We estimate the nuclear physics uncertainties in the i-process model predictions to determine whether the remaining discrepancies with observations are significant and point to potential issues with the underlying astrophysical model. We find that the dominant source in the nuclear physics uncertainties are predictions of neutron capture rates on unstable neutron rich nuclei, which can have uncertainties of more than a factor 20 in the band of the i-process. We use a Monte Carlo variation of 52 neutron capture rates and a 1D multi-zone post-processing model for the i-process in Sakurai's object to determine the cumulative effect of these uncertainties on the final elemental abundance predictions. We find that the nuclear physics uncertainties are large and comparable to observational errors. Within these uncertainties the model predictions are consistent with observations. A correlation analysis of the results of our MC simulations reveals that the strongest impact on the predicted abundances of Rb, Sr, Y and Zr is made by the uncertainties in the (n,$\gamma$) reaction rates of $^{85}$Br, $^{86}$Br, $^{87}$Kr, $^{88}$Kr, $^{89}$Kr, $^{89}$Rb, $^{89}$Sr, and $^{92}$Sr. (abridged)

Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented.

Convective boundary mixing (CBM) in the advanced evolutionary stages of massive stars is not well understood. Structural changes caused by convection have an impact on the evolution as well as the subsequent supernova, or lack thereof. The effects of convectively driven mixing across convective boundaries during the post He core burning evolution of $25\mathrm{M}_{\odot}$, solar-metallicity, non-rotating stellar models is studied using the MESA stellar evolution code. CBM is modelled using the exponentially decaying diffusion coefficient equation, the free parameter of which, $f_{\mathrm{cbm}}$, is varied systematically throughout the course of the stellar model's evolution with values of $(0.002, 0.012, 0.022, 0.032)$. The effect of varying this parameter produces mass ranges at collapse in the ONe, Si, Fe cores of $(1.82\mathrm{M}_{\odot}, 4.36\mathrm{M}_{\odot})$, $(1.67\mathrm{M}_{\odot}, 1.99\mathrm{M}_{\odot})$ and $(1.46\mathrm{M}_{\odot}, 1.70\mathrm{M}_{\odot})$ respectively, with percent differences from the model with minimal CBM as large as 86.3%. At the presupernova stage, the compactness of the stellar cores from O'Connor & Ott (2011), $\xi_M$, exhibit a range of $(0.120, 0.354)$, suggesting that the extent of CBM in the advanced burning stages of massive stars is an important consideration for the explodability and type of compact remnant. The nucleosynthetic yields from the models, most notably C, O, Ne, Mg and Si are also significantly affected by the CBM assumptions, showing non-linear trends with increased mixing. The simulations show that interactions between convective C, Ne and O shells produce significant non-linear changes in the evolution, whereas from the end of Si burning, the structural changes attributed to the CBM are dominated by the growth of the convective C shell. Progenitor structures for all the models are available from HERE (link and DOI to appear).

We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.

Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $\mu m$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3\sigma$ level. Based on the IR SED decomposition approach, we derive the IR luminosities of AGN torus and cold dust components. Hot DOGs in our sample are extremely luminous with most of them having $L_{\rm IR}^{\rm tot}>10^{14} L_\odot$. The torus emissions dominate the total IR energy output. However, the cold dust contribution is still non-negligible, with the fraction of the cold dust contribution to the total IR luminosity $(\sim 8-24\%)$ being dependent on the choice of torus model. The derived cold dust temperatures in Hot DOGs are comparable to those in UV bright quasars with similar IR luminosity, but much higher than those in SMGs. Higher dust temperatures in Hot DOGs may be due to the more intense radiation field caused by intense starburst and obscured AGN activities. Fourteen and five submillimeter serendipitous sources in the 10 SCUBA-2 fields around Hot DOGs have been detected at $>3\sigma$ and $>3.5\sigma$ levels, respectively. By estimating their cumulative number counts, we confirm the previous argument that Hot DOGs lie in dense environments. Our results support the scenario in which Hot DOGs are luminous, dust-obscured quasars lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars.

We have reported previously on a new method we are developing for using image-based information to improve global coronal magnetic field models. In that work we presented early tests of the method which proved its capability to improve global models based on flawed synoptic magnetograms, given excellent constraints on the field in the model volume. In this follow-up paper we present the results of similar tests given field constraints of a nature that could realistically be obtained from quality white-light coronagraph images of the lower corona. We pay particular attention to difficulties associated with the line-of-sight projection of features outside of the assumed coronagraph image plane, and the effect on the outcome of the optimization of errors in localization of constraints. We find that substantial improvement in the model field can be achieved with this type of constraints, even when magnetic features in the images are located outside of the image plane.

We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$.

We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space densities in both samples to be overdense compared to normal star-forming galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy Survey (S2CLS). Both samples of companion sources have consistent mid-IR colours and mid-IR to submm ratios as SMGs. The brighter population around WISE/radio AGNs could be responsible for the higher overdensity reported. We also find the star formation rate density (SFRDs) are higher than the field, but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be good signposts for protoclusters at high redshift on arcmin scales. The results reported here provide an upper limit to the strength of angular clustering using the two-point correlation function. Monte Carlo simulations show no angular correlation, which could indicate protoclusters on scales larger than the SCUBA-2 1.5arcmin scale maps.

The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation.

Based on stellar evolution simulations, we demonstrate that rapidly accreting white dwarfs in close binary systems are an astrophysical site for the intermediate neutron-capture process. During recurrent and very strong He-shell flashes in the stable H-burning accretion regime H-rich material enters the He-shell flash convection zone. $^{12}$C(p,$\gamma)^{13}$N reactions release enough energy to potentially impact convection, and i process is activated through the $^{13}$C($\alpha$,n)$^{16}$O reaction. The H-ingestion flash may not cause a split of the convection zone as it was seen in simulations of He-shell flashes in post-AGB and low-Z AGB stars. We estimate that for the production of first-peak heavy elements this site can be of similar importance for galactic chemical evolution as the s-process production by low-mass AGB stars. The He-shell flashes result in the expansion and, ultimately, ejection of the accreted and then i-process enriched material, via super-Eddington luminosity winds or Roche-lobe overflow. The white dwarf models do not retain any significant amount of the accreted mass, with a He retention efficiency of $\leq 10\%$ depending on mass and convective boundary mixing assumptions. This makes the evolutionary path of such systems to supernova Ia explosion highly unlikely.

This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4\pi$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2\sigma$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles.

The origin of the low luminosity radio emission in radio-quiet AGN, is unknown. The detection of a positive correlation between the radio and X-ray emission would imply a jet-like origin, similar to that seen in `hard state' X-ray binary systems. In our previous work, we found no believable radio variability in the well known X-ray bright Seyfert 1 galaxy NGC 4051, despite large amplitude X-ray variability. In this study we have carefully re-analysed radio and X-ray observations using the same methods as our previous work, we again find no evidence for core radio variability. In direct contrast to our findings, another study claim significant radio variability and a distinctive anti-correlation between radio and X-ray data for the same source. The other study report only integral flux values and do not consider the effect of the changing array on the synthesised beam. In both our studies of NGC 4051 we have taken great care to account for the effect that the changing beam size has on the measured radio flux and as a result we are confident that our method gives more accurate values for the intrinsic core radio flux. However, the lack of radio variability we find is hard to reconcile because radio images of NGC 4051 do show jet-like structure. We suggest that the radio structures observed are likely the result of a previous period of higher radio activity and that the current level of radio emission from a compact nuclear jet is low.

Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the JINA-NuGrid chemical evolution pipeline. It is designed to probe the impact of nuclear astrophysics uncertainties on galactic chemical evolution, to improve our knowledges regarding the origin of the elements in a cosmological context, and to create the required interdisciplinary connections.

In the ONeMg cores of $8.8-9.5~{\rm M}_\odot$ stars, neon and oxygen burning is ignited off-center. Whether the neon-oxygen flame propagates to the center is critical to determine whether these stars undergo Fe core collapse or electron capture induced ONeMg core collapse. We present more details of stars that ignite neon and oxygen burning off-center. The neon flame is established in a similar manner to the carbon flame of super-AGB stars, albeit with a narrower flame width. The criteria for establishing a flame are able to be met if the strict Schwarzschild criterion for convective instability is adopted. Mixing across the interface of the convective shell disrupts the conditions for the propagation of the burning front and instead the shell burns as a series of inward-moving flashes. While this may not directly affect whether the burning will reach the center (as in super-AGB stars), the core is allowed to contract between each shell flash. Reduction of the electron fraction in the shell reduces the Chandrasekhar mass and the center reaches the threshold density for the URCA process to activate and steer the remaining evolution of the core. This highlights the importance of a more accurate treatment of mixing in the stellar interior for yet another important question in stellar astrophysics - determining the properties of stellar evolution and supernova progenitors at the boundary between electron capture supernova and iron core-collapse supernova.

In the classical picture, electron-capture supernovae and the accretion-induced collapse of oxygen-neon white dwarfs (ONeWDs) undergo an oxygen deflagration phase before gravitational collapse produces a neutron star (NS). These types of core collapse events are postulated to explain several astronomical phenomena. In this work, the deflagration phase is simulated for the first time using multidimensional hydrodynamics, with the aim of gaining new insight into the explosive deaths of $8-10~M_\odot$ stars and ONeWDs that accrete material from a binary companion star. The main aim is to determine whether these events are thermonuclear or core-collapse supernova explosions, and hence whether NSs are formed by such phenomena. The deflagration is simulated in ONe cores with three different central ignition densities. The intermediate density case is perhaps the most realistic, being based on recent nuclear physics calculations and 1D stellar models. The 3D hydrodynamic simulations presented in this work begin from a centrally confined flame structure using a level-set-based flame approach and are performed in $256^3$ and $512^3$ numerical resolutions. In the simulations with intermediate and low ignition density, the cores do not appear to collapse into NSs. Instead, almost a solar mass of material becomes unbound from the cores, leaving bound remnants. These simulations represent the case in which semiconvective mixing during the electron-capture phase preceding the deflagration is inefficient. The masses of the bound remnants double when Coulomb corrections are included in the EoS, however they still do not exceed the effective Chandrasekhar mass and, hence, would not collapse into NSs. The simulations with the highest ignition density ($\log_{10}\rho_{\rm c}=10.3$), representing the case where semiconvective mixing is very efficient, show clear signs that the core will collapse into a NS.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

The coronal magnetic field directly or indirectly affects a majority of the phenomena studied in space physics. It provides energy for coronal heating, controls the release of coronal mass ejections (CMEs), and drives heliospheric and magnetospheric activity, yet the coronal magnetic field itself has proven difficult to measure. This difficulty has prompted a decades-long effort to develop accurate, timely, models of the field - an effort that continues today. We have developed a method for improving global coronal magnetic field models by incorporating the type of morphological constraints which could be derived from coronal images. Here we report promising initial tests of this approach on two theoretical problems, and discuss opportunities for application.

We investigate the evolution of super-AGB thermal pulse (TP) stars for a range of metallicities (Z) and explore the effect of convective boundary mixing (CBM). With decreasing metallicity and evolution along the TP phase, the He-shell flash and the third dredge-up (TDU) occur closer together in time. After some time (depending upon the CBM parameterisation), efficient TDU begins while the pulse-driven convection zone (PDCZ) is still present, causing a convective exchange of material between the PDCZ and the convective envelope. This results in the ingestion of protons into the convective He-burning pulse. Even small amounts of CBM encourage the interaction of the convection zones leading to transport of protons from the convective envelope into the He layer. H-burning luminosities exceed $10^9$ (in some cases $10^{10}$) $\mathrm{L}_\odot$. We also calculate models of dredge-out in the most massive super-AGB stars and show that the dredge-out phenomenon is another likely site of convective-reactive H-$^{12}$C combustion. We discuss the substantial uncertainties of stellar evolution models under these conditions. Nevertheless, the simulations suggest that in the convective-reactive H-combustion regime of H ingestion the star may encounter conditions for the intermediate neutron capture process (i process). We speculate that some CEMP-s/r stars could originate in i-process conditions in the H-ingestion phases of low-Z SAGB stars. This scenario would however suggest a very low electron-capture supernova rate from super-AGB stars. We also simulate potential outbursts triggered by such H-ingestion events, present their light curves and briefly discuss their transient properties.

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-s+r (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars, super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found that single-exposure one-zone simulations of the i-process reproduce the abundance patterns in some of the CEMP-r/s stars much better than the model that assumes a superposition of yields from s- and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-process revealed that they could have a significant impact on the i-process yields obtained in our idealized one-zone calculations, leading, for example, to ~0.7dex uncertainty in our predicted [Ba/La] ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGB Sakurai's object have discovered a new mode of non-radial instabilities: the Global Oscillation of Shell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulations cannot be used to accurately model physical conditions for the i-process.

Submillimetre (submm) observations of WISE-selected, dusty, luminous, high-redshift galaxies have revealed intriguing overdensities around them on arcmin scales. They could be the best signposts of overdense environments on the sky.

We present JCMT SCUBA-2 850microns submillimetre (submm) observations of 30 mid-infrared (mid-IR) luminous AGN, detected jointly by the WISE all-sky IR survey and the NVSS/FIRST radio survey. These rare sources are selected by their extremely red mid-infrared spectral energy distributions (SEDs) and compact radio counterparts. Further investigations show that they are highly obscured, have abundant warm AGN-heated dust and are thought to be experiencing intense AGN feedback. These galaxies appear to be consistent with an AGN-dominated galaxy, and could be a transient phase of merging galaxies. When comparing the number of submm galaxies (SMGs) detected serendipitously in the surrounding 1.5-arcmin to those in blank-field submm surveys, there is a very significant overdensity, of order 5, but no sign of radial clustering centred at our primary objects. The WISE/radio-selected AGN thus reside in 10-Mpc-scale overdense environments, that could be forming in pre-viralised clusters of galaxies. WISE/radio-selected AGNs appear to be the strongest signposts of high-density regions of active, luminous and dusty galaxies. SCUBA-2 850microns observations indicate that their submm fluxes are low compared to many popular AGN SED templates, hence the WISE/radio-selected AGNs have either less cold and/or more warm dust emission than normally assumed for typical AGN. Most of the targets are not detected, only four targets are detected at SCUBA-2 850microns, and have total IR luminosities >= 10^13 L_solar, if their redshifts are consistent with the subset of the 10 SCUBA-2 undetected targets with known redshifts, z ~ 0.44 - 2.86.

Massive stars are key sources of radiative, kinetic, and chemical feedback in the universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes---GENEC, KEPLER and MESA---to compute non-rotating stellar models of $15~\mathrm{M}_\odot$, $20~\mathrm{M}_\odot$, and $25~\mathrm{M}_\odot$ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The GENEC and KEPLER models hold physics assumptions used in large grids of published models. The MESA code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by GENEC. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool MPPNP. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than $30\,\%$---smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g., helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models.

When carbon is ignited off-centre in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the center. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-center C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the center. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

Context. The high intensities and narrow bandwidths exhibited by some astronomical masers make them ideal tools for studying star-forming giant molecular clouds. The water maser transition $J_{K_{-1}K_{1}}=5_{32}-4_{41}$ at 620.701 GHz can only be observed from above Earth's strongly absorbing atmosphere; its emission has recently been detected from space. Aims. We sought to further characterize the star-forming environment of Orion KL by investigating the linear polarisation of a source emitting a narrow 620.701 GHz maser feature with the heterodyne spectrometer HIFI on board the Herschel Space Observatory. Methods. High-resolution spectral datasets were collected over a thirteen month period beginning in 2011 March, to establish not only the linear polarisation but also the temporal variability of the source. Results. Within a $3\sigma$ uncertainty, no polarisation was detected to an upper limit of approximately 2%. These results are compared with coeval linear polarisation measurements of the 22.235 GHz $J_{K_{-1}K_{1}}=6_{16}-5_{23}$ maser line from the Effelsberg 100-m radio telescope, typically a much stronger maser transition. Although strongly polarised emission is observed for one component of the 22.235 GHz maser at 7.2 km s$^{-1}$, a weaker component at the same velocity as the 620.701 GHz maser at 11.7 km s$^{-1}$ is much less polarised.

We present SCUBA-2 850um submillimetre (submm) observations of the fields of 10 dusty, luminous galaxies at z ~ 1.7 - 4.6, detected at 12um and/or 22um by the WISE all-sky survey, but faint or undetected at 3.4um and 4.6um; dubbed hot, dust-obscured galaxies (Hot DOGs). The six detected targets all have total infrared luminosities greater than 10^13 L_sun, with one greater than 10^14 L_sun. Their spectral energy distributions (SEDs) are very blue from mid-infrared to submm wavelengths and not well fitted by standard AGN SED templates, without adding extra dust extinction to fit the WISE 3.4um and 4.6um data. The SCUBA-2 850um observations confirm that the Hot DOGs have less cold and/or more warm dust emission than standard AGN templates, and limit an underlying extended spiral or ULIRG-type galaxy to contribute less than about 2% or 55% of the typical total Hot DOG IR luminosity, respectively. The two most distant and luminous targets have similar observed submm to mid-infrared ratios to the rest, and thus appear to have even hotter SEDs. The number of serendipitous submm galaxies (SMGs) detected in the 1.5-arcmin-radius SCUBA-2 850um maps indicates there is a significant over-density of serendipitous sources around Hot DOGs. These submm observations confirm that the WISE-selected ultra-luminous galaxies have very blue mid-infrared to submm SEDs, suggesting that they contain very powerful AGN, and are apparently located in unusual arcmin-scale overdensities of very luminous dusty galaxies.

The stellar mass range 8<M/Mo<12 corresponds to the most massive AGB stars and the most numerous massive stars. It is host to a variety of supernova progenitors and is therefore very important for galactic chemical evolution and stellar population studies. In this paper, we study the transition from super-AGB star to massive star and find that a propagating neon-oxygen burning shell is common to both the most massive electron capture supernova (EC-SN) progenitors and the lowest mass iron-core collapse supernova (FeCCSN) progenitors. Of the models that ignite neon burning off-center, the 9.5Mo model would evolve to an FeCCSN after the neon-burning shell propagates to the center, as in previous studies. The neon-burning shell in the 8.8Mo model, however, fails to reach the center as the URCA process and an extended (0.6 Mo) region of low Ye (0.48) in the outer part of the core begin to dominate the late evolution; the model evolves to an EC-SN. This is the first study to follow the most massive EC-SN progenitors to collapse, representing an evolutionary path to EC-SN in addition to that from SAGB stars undergoing thermal pulses. We also present models of an 8.75Mo super-AGB star through its entire thermal pulse phase until electron captures on 20Ne begin at its center and of a 12Mo star up to the iron core collapse. We discuss key uncertainties and how the different pathways to collapse affect the pre-supernova structure. Finally, we compare our results to the observed neutron star mass distribution.

We present measurements of circular polarization from rotational spectral lines of molecular species in Orion KL, most notably 12CO (J=2 - 1), obtained at the Caltech Submillimeter Observatory with the Four-Stokes-Parameter Spectra Line Polarimeter. We find levels of polarization of up to 1 to 2% in general, for 12CO (J=2 - 1) this level is comparable to that of linear polarization also measured for that line. We present a physical model based on resonant scattering in an attempt to explain our observations. We discuss how slight differences in scattering amplitudes for radiation polarized parallel and perpendicular to the ambient magnetic field, responsible for the alignment of the scattering molecules, can lead to the observed circular polarization. We also show that the effect is proportional to the square of the magnitude of the plane of the sky component of the magnetic field, and therefore opens up the possibility of measuring this parameter from circular polarization measurements of Zeeman insensitive molecules.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

The Census of High- and Medium-mass Protostars (CHaMP) is the first large-scale, unbiased, uniform mapping survey at sub-parsec scale resolution of 90 GHz line emission from massive molecular clumps in the Milky Way. We present the first Mopra (ATNF) maps of the CHaMP survey region (300{\deg}>l>280{\deg}) in the HCO+ J=1-0 line, which is usually thought to trace gas at densities up to 10^11 m-3. In this paper we introduce the survey and its strategy, describe the observational and data reduction procedures, and give a complete catalogue of moment maps of the HCO+ J=1-0 emission from the ensemble of 301 massive molecular clumps. From these maps we also derive the physical parameters of the clumps, using standard molecular spectral-line analysis techniques. This analysis yields the following range of properties: integrated line intensity 1-30 K km s-1, peak line brightness 1-7 K, linewidth 1-10 km s-1, integrated line luminosity 0.5-200 K km s-1 pc^2, FWHM size 0.2-2.5 pc, mean projected axial ratio 2, optical depth 0.08-2, total surface density 30-3000 M{\sun} pc-2, number density 0.2-30 x 10^9 m-3, mass 15-8000 M{\sun}, virial parameter 1-55, and total gas pressure 0.3-700 pPa. We find that the CHaMP clumps do not obey a Larson-type size-linewidth relation. Among the clumps, there exists a large population of subthermally excited, weakly-emitting (but easily detectable) dense molecular clumps, confirming the prediction of Narayanan et al. (2008). These weakly-emitting clumps comprise 95% of all massive clumps by number, and 87% of the molecular mass, in this portion of the Galaxy; their properties are distinct from the brighter massive star-forming regions that are more typically studied. If the clumps evolve by slow contraction, the 95% of fainter clumps may represent a long-lived stage of pressure-confined, gravitationally stable massive clump evolution, while the CHaMP ... (abridged)

We present intensive quasi-simultaneous X-ray and radio monitoring of the narrow line Seyfert 1 galaxy NGC 4051, over a 16 month period in 2000-2001. Observations were made with the Rossi Timing X-ray Explorer (RXTE) and the Very Large Array (VLA) at 8.4 and 4.8 GHz. In the X-ray band NGC 4051 behaves much like a Galactic black hole binary (GBH) system in a `soft-state'. In such systems, there has so far been no firm evidence for an active, radio-emitting jet like those found in `hard state' GBHs. VLBI observations of NGC 4051 show three co-linear compact components. This structure resembles the core and outer hot spots seen in powerful, jet-dominated, extragalactic radio sources and suggests the existence of a weak jet. Radio monitoring of the core of NGC 4051 is complicated by the presence of surrounding extended emission and by the changing array configurations of the VLA. Only in the A configuration is the core reasonably resolved. We have carefully removed the contaminations of the core by extended emission in the various arrays. The resulting lightcurve shows no sign of large amplitude variability (i.e. factor 50 %) over the 16 month period. Within the most sensitive configuration (A array) we see marginal evidence for radio core variability of ~25% (~0.12 mJy at 8.4GHz) on a 2-week timescale, correlated with X-ray variations. Even if the radio variations in NGC 4051 are real, the percentage variability is much less than in the X-ray band. Within the B configuration observations, where sensitivity is reduced, there is no sign of correlated X-ray/radio variability. The lack of radio variability in NGC 4051, which we commonly see in `hard state' GBHs, may be explained by orientation effects. Another possibility is that the radio emission arises from the X-ray corona, although the linear structure of the compact radio components here is hard to explain.

SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {\mu}m allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

An extended data set of extreme ultraviolet images of the solar corona provided by the SOHO spacecraft are analyzed using statistical methods common to studies of self-organized criticality (SOC) and intermittent turbulence (IT). The data exhibits simultaneous hallmarks of both regimes, namely power law avalanche statistics as well as multiscaling of structure functions for spatial activity. This implies that both SOC and IT may be manifestations of a single complex dynamical process entangling avalanches of magnetic energy dissipation with turbulent particle flows.

Contemporary interstellar dust has never been analyzed in the laboratory, despite its obvious astronomical importance and its potential as a probe of stellar nucleosynthesis and galactic chemical evolution. Here we report the discovery of a novel fluorescent aerogel which is capable of capturing hypervelocity dust grains and passively recording their kinetic energies. An array of these "calorimetric" aerogel collectors in low earth orbit would lead to the capture and identification of large numbers of interstellar dust grains.

Using the Very Large Array (VLA), we have detected radio emission from the site of SN 1961V in the Sc galaxy NGC 1058. With a peak flux density of 0.063 +/- 0.008 mJy/beam at 6 cm and 0.147 +/- 0.026 mJy/beam at 18 cm, the source is non-thermal, with a spectral index of -0.79 +/- 0.23. Within errors, this spectral index is the same value reported for previous VLA observations taken in 1984 and 1986. The radio emission at both wavelengths has decayed since the mid 1980's observations with power-law indices of beta(20cm) = -0.69 +/- 0.23 and beta(6cm) = -1.75 +/- 0.16. We discuss the radio properties of this source and compare them with those of Type II radio supernovae and luminous blue variables.

We present results from X-ray observations of the quadruply lensed quasar RX J0911.4+0551 using data obtained with the Advanced CCD Imaging Spectrometer (ACIS) on board the Chandra X-ray Observatory. The 29 ks observation detects a total of ~404 X-ray photons (0.3 to 7.0 keV) from the four images of the lensed quasar. Deconvolution of the aspect corrected data resolves all four lensed images, with relative positions in good agreement with optical measurements. When compared to contemporaneous optical data, one of the lensed images (component A3) is dimmer by a factor of ~6 in X-rays with respect to the 2 brighter images (components A1 and A2). Spectral fitting for the combined images shows significant intrinsic absorption in the soft (0.2 to 2.4 keV) energy band, consistent with the mini-BAL nature of this quasar, while a comparison with ROSAT PSPC observations from 1990 shows a drop of ~6.5 in the total soft bandpass flux. The observations also detect ~157 X-ray photons arising from extended emission of the nearby cluster (peaked ~42" SW of RXJ0911.4+0551) responsible for the large external shear present in the system. The Chandra observation reveals the cluster emission to be complex and non-spherical, and yields a cluster temperature of kT = 2.3^{+1.8}_{-0.8} keV and a 2.0 to 10 keV cluster luminosity within a 1 Mpc radius of L_X = 7.6_{-0.2}^{+0.6} x 10^{43} ergs/s (error bars denote 90% confidence limits). Our mass estimate of the cluster within its virial radius is 2.3^{+1.8}_{-0.7} x 10^{14} solar, and is a factor of 2 smaller than, although consistent with, previous mass estimates based on the observed cluster velocity dispersion.