Supernova spectral time series contain a wealth of information about the progenitor and explosion process of these energetic events. The modeling of these data requires the exploration of very high dimensional posterior probabilities with expensive radiative transfer codes. Even modest parametrizations of supernovae contain more than ten parameters and a detailed exploration demands at least several million function evaluations. Physically realistic models require at least tens of CPU minutes per evaluation putting a detailed reconstruction of the explosion out of reach of traditional methodology. The advent of widely available libraries for the training of neural networks combined with their ability to approximate almost arbitrary functions with high precision allows for a new approach to this problem. Instead of evaluating the radiative transfer model itself, one can build a neural network proxy trained on the simulations but evaluating orders of magnitude faster. Such a framework is called an emulator or surrogate model. In this work, we present an emulator for the TARDIS supernova radiative transfer code applied to Type Ia supernova spectra. We show that we can train an emulator for this problem given a modest training set of a hundred thousand spectra (easily calculable on modern supercomputers). The results show an accuracy on the percent level (that are dominated by the Monte Carlo nature of TARDIS and not the emulator) with a speedup of several orders of magnitude. This method has a much broader set of applications and is not limited to the presented problem.

While ancient scientists often had patrons to fund their work, peer review of proposals for the allocation of resources is a foundation of modern science. A very common method is that proposals are evaluated by a small panel of experts (due to logistics and funding limitations) nominated by the grant-giving institutions. The expert panel process introduces several issues - most notably: 1) biases introduced in the selection of the panel. 2) experts have to read a very large number of proposals. Distributed Peer Review promises to alleviate several of the described problems by distributing the task of reviewing among the proposers. Each proposer is given a limited number of proposals to review and rank. We present the result of an experiment running a machine-learning enhanced distributed peer review process for allocation of telescope time at the European Southern Observatory. In this work, we show that the distributed peer review is statistically the same as a `traditional' panel, that our machine learning algorithm can predict expertise of reviewers with a high success rate, and we find that seniority and reviewer expertise have an influence on review quality. The general experience has been overwhelmingly praised from the participating community (using an anonymous feedback mechanism).

We present stellar metallicity measurements of more than 600 late-type stars in the central 10 pc of the Galactic centre. Together with our previously published KMOS data, this data set allows us to investigate, for the first time, spatial variations of the nuclear star cluster's metallicity distribution. Using the integral-field spectrograph KMOS (VLT) we observed almost half of the area enclosed by the nuclear star cluster's effective radius. We extract spectra at medium spectral resolution, and apply full spectral fitting utilising the PHOENIX library of synthetic stellar spectra. The stellar metallicities range from [M/H]=-1.25 dex to [M/H]> +0.3 dex, with most of the stars having super-solar metallicity. We are able to measure an anisotropy of the stellar metallicity distribution. In the Galactic North, the portion of sub-solar metallicity stars with [M/H]<0.0 dex is more than twice as high as in the Galactic South. One possible explanation for different fractions of sub-solar metallicity stars in different parts of the cluster is a recent merger event. We propose to test this hypothesis with high-resolution spectroscopy, and by combining the metallicity information with kinematic data.

There are now hundreds of publicly available supernova spectral time series. Radiative transfer modeling of this data gives insights into the physical properties of these explosions such as the composition, the density structure, or the intrinsic luminosity---this is invaluable for understanding the supernova progenitors, the explosion mechanism, or for constraining the supernova distance. However, a detailed parameter study of the available data has been out of reach due to the high dimensionality of the problem coupled with the still significant computational expense. We tackle this issue through the use of machine-learning emulators, which are algorithms for high-dimensional interpolation. These use a pre-calculated training dataset to mimic the output of a complex code but with run times orders of magnitude shorter. We present the application of such an emulator to synthetic type II supernova spectra generated with the TARDIS radiative transfer code. The results show that with a relatively small training set of 780 spectra we can generate emulated spectra with interpolation uncertainties of less than one percent. We demonstrate the utility of this method by automatic spectral fitting of two well-known type IIP supernovae; as an exemplary application, we determine the supernova distances from the spectral fits using the tailored-expanding-photosphere method. We compare our results to previous studies and find good agreement. This suggests that emulation of TARDIS spectra can likely be used to perform automatic and detailed analysis of many transient classes putting the analysis of large data repositories within reach.

We have used two methods to search for surviving companions of Type Ia supernova progenitors in three Balmer-dominated supernova remnants (SNRs) in the Large Magellanic Cloud: 0519-69.0, 0505-67.9 (DEM L71), and 0548-70.4. In the first method, we use the Hubble Space Telescope photometric measurements of stars to construct color-magnitude diagrams (CMDs), and compare positions of stars in the CMDs with those expected from theoretical post-impact evolution of surviving main sequence or helium star companions. No obvious candidates of surviving companion are identified in this photometric search. Future models for surviving red giant companions or with different explosion mechanisms are needed for thorough comparisons with these observations in order to make more definitive conclusions. In the second method, we use Multi-Unit Spectroscopic Explorer (MUSE) observations of 0519-69.0 and DEM L71 to carry out spectroscopic analyses of stars in order to use large peculiar radial velocities as diagnostics of surviving companions. We find a star in 0519-69.0 and a star in DEM L71 moving at radial velocities of 182 $\pm$ 0 km s$^{-1}$ and 213 $\pm$ 0 km s$^{-1}$, more than 2.5$\sigma$ from the mean radial velocity of the underlying stellar population, 264 km s$^{-1}$ and 270 km s$^{-1}$, respectively. These stars need higher-quality spectra to investigate their abundances and rotation velocities to determine whether they are indeed surviving companions of the SN progenitors.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

General Relativity predicts that a star passing close to a supermassive black hole should exhibit a relativistic redshift. We test this using observations of the Galactic center star S0-2. We combine existing spectroscopic and astrometric measurements from 1995-2017, which cover S0-2's 16-year orbit, with measurements in 2018 March to September which cover three events during its closest approach to the black hole. We detect the combination of special relativistic- and gravitational-redshift, quantified using a redshift parameter, $\Upsilon$. Our result, $\Upsilon=0.88 \pm 0.17$, is consistent with General Relativity ($\Upsilon=1$) and excludes a Newtonian model ($\Upsilon=0$ ) with a statistical significance of 5 $\sigma$.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

How an otherwise inert carbon-oxygen white dwarf can be made to explode as a Type Ia supernova remains unknown. A promising test of theoretical models is to constrain the distribution of material that is left unburned, in particular of carbon. So far, most investigations used line identification codes to detect carbon in the ejecta, a method that cannot be readily compared against model predictions because it requires assumed opacities and temperatures. Here, we instead use tomographic techniques to investigate the amount of carbon in the inner layers of SN~2011fe, starting from the previously published tomographic analysis of Mazzali et al. From the presence of the carbon feature in the optical at early epochs and its disappearance later on, we derive an average carbon mass fraction between 0.001 and 0.05 for velocities in the range $13500 \lesssim v \lesssim 16000\ \rm{km\ s^{-1}}$, and an upper limit of 0.005 inside that region. Based on our models and the assumed density profile, only small amounts of carbon should be in the neutral state, too little to be responsible for features seen in near-infrared spectra that were previously identified as due to neutral carbon; We discuss possible reasons for this discrepancy. We compare our results against a suite of explosion models, although uncertainties in both the models and our simulations make it difficult to draw definitive conclusions.

Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as "stripped-envelope supernovae", showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-rich progenitors are too large to fit. The physical processes responsible for the stripping process and the fate of the companion are still very poorly understood. Aiming to find new clues, we investigate Cas~A, which is a very young ($\sim$340 \,yr) and near ($\sim$3.4\,kpc) remnant of a core collapse supernova. Cas~A has been subject to several searches for possible companions, all unsuccessfully. We present new measurements of the proper motions and photometry of stars in the vicinity based on deep HST ACS/WFC and WFC3-IR data. We identify stellar sources that are close enough in projection, but using their proper motions we show that none are compatible with being at the location of center at the time of explosion, in agreement with earlier findings. Our photometric measurements allow us to place much deeper (order of magnitude) upper limits on the brightness of possible undetected companions. We systematically compare them with model predictions for a wide variety of scenarios. We can confidently rule out the presence of any stellar companion of any reasonable mass and age (main sequence, pre main sequence or stripped). Although this finding is not in direct conflict with model predictions, it does rule out what many considered to be the most likely evolutionary scenario for type IIb supernova. (abstract abbreviated)

We present substantial extensions to the Monte Carlo radiative transfer code TARDIS to perform spectral synthesis for type II supernovae. By incorporating a non-LTE ionization and excitation treatment for hydrogen, a full account of free-free and bound-free processes, a self-consistent determination of the thermal state and by improving the handling of relativistic effects, the improved code version includes the necessary physics to perform spectral synthesis for type II supernovae to high precision as required for the reliable inference of supernova properties. We demonstrate the capabilities of the extended version of TARDIS by calculating synthetic spectra for the prototypical type II supernova SN1999em and by deriving a new and independent set of dilution factors for the expanding photosphere method. We have investigated in detail the dependence of the dilution factors on photospheric properties and, for the first time, on changes in metallicity. We also compare our results with two previously published sets of dilution factors by Eastman et al. (1996) and by Dessart & Hillier (2005), and discuss the potential sources of the discrepancies between studies.

We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a sampling algorithm allowing us to probe a broad parameter space. We derive velocity shifts, line widths, and abundance ratios for iron and cobalt. The measured line widths and velocity shifts of the singly ionised ions suggest a shared emitting region. Our data are fully compatible with radioactive $^{56}$Ni decay as the origin for cobalt and iron. We compare the measured abundance ratios of iron and cobalt to theoretical predictions of various SN$\,$Ia explosion models. These models include, in addition to $^{56}$Ni, different amounts of $^{57}$Ni and stable $^{54,56}$Fe. We can exclude models that produced only $^{54,56}$Fe or only $^{57}$Ni in addition to $^{56}$Ni. If we consider a model that has $^{56}$Ni, $^{57}$Ni, and $^{54,56}$Fe then our data imply that these ratios are $^{54,56}$Fe / $^{56}$Ni $=0.272\pm0.086$ and $^{57}$Ni / $^{56}$Ni $=0.032\pm0.011$.

Recent studies have demonstrated the diversity in type Ia supernovae (SNe Ia) at early times and highlighted a need for a better understanding of the explosion physics as manifested by observations soon after explosion. To this end, we present a Monte Carlo code designed to model the light curves of radioactively driven, hydrogen-free transients from explosion to approximately maximum light. In this initial study, we have used a parametrised description of the ejecta in SNe Ia, and performed a parameter study of the effects of the $^{56}$Ni distribution on the observed colours and light curves for a fixed $^{56}$Ni mass of 0.6 $M_\odot$. For a given density profile, we find that models with $^{56}$Ni extending throughout the entirety of the ejecta are typically brighter and bluer shortly after explosion. Additionally, the shape of the density profile itself also plays an important role in determining the shape, rise time, and colours of observed light curves. We find that the multi-band light curves of at least one SNe Ia (SN 2009ig) are inconsistent with less extended $^{56}$Ni distributions, but show good agreement with models that incorporate $^{56}$Ni throughout the entire ejecta. We further demonstrate that comparisons with full $UVOIR$ colour light curves are powerful tools in discriminating various $^{56}$Ni distributions, and hence explosion models.

The increase in the number of researchers coupled with the ease of publishing and distribution of scientific papers (due to technological advancements) has resulted in a dramatic increase in astronomy literature. This has likely led to the predicament that the body of the literature is too large for traditional human consumption and that related and crucial knowledge is not discovered by researchers. In addition to the increased production of astronomical literature, recent decades have also brought several advancements in computational linguistics. Especially, the machine-aided processing of literature dissemination might make it possible to convert this stream of papers into a coherent knowledge set. In this paper, we present the application of computational linguistics techniques to astronomy literature. In particular, we developed a tool that will find similar articles purely based on text content from an input paper. We find that our technique performs robustly in comparison with other tools recommending articles given a reference paper (known as recommender system). Our novel tool shows the great power in combining computational linguistics with astronomy literature and suggests that additional research in this endeavor will likely produce even better tools that will help researchers cope with the vast amounts of knowledge being produced.

We present direct spectroscopic modeling of five Type Iax supernovae (SNe) with the one dimensional Monte Carlo radiative transfer code TARDIS. The abundance tomography technique is used to map the chemical structure and physical properties of the SN atmosphere. Through via fitting of multiple spectral epochs with self-consistent ejecta models, we can then constrain the location of some elements within the ejecta. The synthetic spectra of the best-fit models are able to reproduce the flux continuum and the main absorption features in the whole sample. We find that the mass fractions of IGEs and IMEs show a decreasing trend toward the outer regions of the atmospheres using density profiles similar to those of deflagration models in the literature. Oxygen is the only element, which could be dominant at higher velocities. The stratified abundance structure contradicts the well-mixed chemical profiles predicted by pure deflagration models. Based on the derived densities and abundances, a template model atmosphere is created for the SN Iax class and compared to the observed spectra. Free parameters are the scaling of the density profile, the velocity shift of the abundance template, and the peak luminosity. The results of this test support the idea that all SNe Iax can be described by a similar internal structure, which argues for a common origin of this class of explosions.

If some or all Type Ia SNe arise from accretion onto a massive WD from a companion, then the companion will remain in some form after the SN explosion. Tycho-B is an unusual, relatively hot star along the line of sight to Tycho's SNR -- conclusively shown to be a Type Ia -- and has been suggested as such a companion. If the interior of Tycho's SNR contains unshocked Fe, and if Tycho-B is either within the SNR shell or in the background, then one might hope to see evidence of this in the UV spectrum. Such is the case for SN\,1006, where spectra of the background Schweizer-Middleditch star, as well as two AGNs, show broad absorption lines of Fe II. To test this idea, we have used STIS on HST to obtain a UV spectrum of Tycho-B. The observed spectrum, however, shows no evidence of Fe II absorption. Furthermore, a luminosity distance estimate using UV and optical spectra of Tycho-B suggests that the star is consistent with a foreground interloper. We conclude either that Tycho B is nearer than Tycho's SNR, or that all of the Fe in the interior of Tycho's SNR is more highly ionized.

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.

We present near infrared (NIR) spectroscopy of the nearby supernova 2014J obtained $\sim$450 d after explosion. We detect the [Ni II] 1.939 $\mu$m line in the spectra indicating the presence of stable $^{58}$Ni in the ejecta. The stable nickel is not centrally concentrated but rather distributed as the iron. The spectra are dominated by forbidden [Fe II] and [Co II] lines. We use lines, in the NIR spectra, arising from the same upper energy levels to place constraints on the extinction from host galaxy dust. We find that that our data are in agreement with the high $A_V$ and low $R_V$ found in earlier studies from data near maximum light. Using a $^{56}$Ni mass prior from near maximum light $\gamma$-ray observations, we find $\sim$0.05 M$_\odot$ of stable nickel to be present in the ejecta. We find that the iron group features are redshifted from the host galaxy rest frame by $\sim$600 km s$^{-1}$.

One of the big challenges in astrophysics is the comparison of complex simulations to observations. As many codes do not directly generate observables (e.g. hydrodynamic simulations), the last step in the modelling process is often a radiative-transfer treatment. For this step, the community relies increasingly on Monte Carlo radiative transfer due to the ease of implementation and scalability with computing power. We consider simulations in which the number of photon packets is Poisson distributed, while the weight assigned to a single photon packet follows any distribution of choice. We show how to estimate the statistical uncertainty of the sum of weights in each bin from the output of a single radiative-transfer simulation. Our Bayesian approach produces a posterior distribution that is valid for any number of packets in a bin, even zero packets, and is easy to implement in practice. Our analytic results for large number of packets show that we generalise existing methods that are valid only in limiting cases. The statistical problem considered here appears in identical form in a wide range of Monte Carlo simulations including particle physics and importance sampling. It is particularly powerful in extracting information when the available data are sparse or quantities are small.

We present adaptive-optics assisted near-infrared high-spectral resolution observations of late-type giants in the nuclear star cluster of the Milky Way. The metallicity and elemental abundance measurements of these stars offer us an opportunity to understand the formation and evolution of the nuclear star cluster. In addition, their proximity to the supermassive black hole ($\sim 0.5$ pc) offers a unique probe of the star formation and chemical enrichment in this extreme environment. We observed two stars identified by medium spectral-resolution observations as potentially having very high metallicities. We use spectral-template fitting with the PHOENIX grid and Bayesian inference to simultaneously constrain the overall metallicity, [M/H], alpha-element abundance [$\alpha$/Fe], effective temperature, and surface gravity of these stars. We find that one of the stars has very high metallicity ([M/H] $> 0.6$) and the other is slightly above solar metallicity. Both Galactic center stars have lines from scandium (Sc), vanadium (V), and yttrium (Y) that are much stronger than allowed by the PHOENIX grid. We find, using the spectral synthesis code Spectroscopy Made Easy, that [Sc/Fe] may be an order of magnitude above solar. For comparison, we also observed an empirical calibrator in NGC6791, the highest metallicity cluster known ([M/H] $\sim 0.4$). Most lines are well matched between the calibrator and the Galactic center stars, except for Sc, V, and Y, which confirms that their abundances must be anomalously high in these stars. These unusual abundances, which may be a unique signature of nuclear star clusters, offer an opportunity to test models of chemical enrichment in this region.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mechanism, and flame speed). Theory predicts that at very late times (beyond 1000 d) it might be possible to distinguish between explosion models. Few very nearby supernovae can be observed that long after the explosion. The Type Ia supernova SN 2011fe located in M101 and along a line of sight with negligible extinction, provides us with the once-in-a-lifetime chance to obtain measurements that may distinguish between theoretical models. In this work, we present the analysis of photometric data of SN 2011fe taken between 900 and 1600 days after explosion with Gemini and HST. At these extremely late epochs theory suggests that the light curve shape might be used to measure isotopic abundances which is a useful model discriminant. However, we show in this work that there are several currently not well constrained physical processes introducing large systematic uncertainties to the isotopic abundance measurement. We conclude that without further detailed knowledge of the physical processes at this late stage one cannot reliably exclude any models on the basis of this dataset.

Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surveys, however, have only looked for stars as faint as $\approx 0.1 L_\odot$ and thus would have missed a surviving white dwarf companion. In this work, we present very deep DECAM imaging $(u, g, r, z)$ of the Type Ia supernova remnant SN 1006 specifically to search for a potential surviving white dwarf companion. We find no object within the inner third of the SN 1006 remnant that is consistent with a relatively young cooling white dwarf. We find that if there is a companion white dwarf, it must have formed long ago and cooled undisturbed for $> 10^8$ yr to be consistent with the redder objects in our sample. We conclude that our findings are consistent with the complete destruction of the secondary (such as in a merger) or an unexpectedly cool and thus very dim surviving companion white dwarf.

Type Ia supernovae are generally agreed to arise from thermonuclear explosions of carbon-oxygen white dwarfs. The actual path to explosion, however, remains elusive, with numerous plausible parent systems and explosion mechanisms suggested. Observationally, type Ia supernovae have multiple subclasses, distinguished by their lightcurves and spectra. This raises the question whether these reflect that multiple mechanisms occur in nature, or instead that explosions have a large but continuous range of physical properties. We revisit the idea that normal and 91bg-like supernovae can be understood as part of a spectral sequence, in which changes in temperature dominate. Specifically, we find that a single ejecta structure is sufficient to provide reasonable fits of both the normal type Ia supernova SN~2011fe and the 91bg-like SN~2005bl, provided that the luminosity and thus temperature of the ejecta are adjusted appropriately. This suggests that the outer layers of the ejecta are similar, thus providing some support of a common explosion mechanism. Our spectral sequence also helps to shed light on the conditions under which carbon can be detected in pre-maximum SN~Ia spectra -- we find that emission from iron can "fill in" the carbon trough in cool SN~Ia. This may indicate that the outer layers of the ejecta of events in which carbon is detected are relatively metal poor compared to events where carbon is not detected.

We present a detailed spectral analysis of Type Iax SN 2011ay. Our spectra cover epochs between -3 and +19 days with respect to the maximum light in B-band. This time range allows us to employ a so-called abundance tomography technique. The synthetic spectral fitting was made with the 1D Monte Carlo radiative transfer code TARDIS. In this paper, we describe our method to fit multiple epochs with a self-consistent, stratified atmospheric model. We compare our results to previously published SYN++ models and the predictions of different explosion scenarios. Using a fixed density profile (exponential fit of W7), we find that a uniform abundance model cannot reproduce the spectral features before maximum light because of the emergence of excessively strong Fe lines. In our best-fit TARDIS model, we find an abundance profile that separated into two different regimes: a well-mixed region under 10,000 km s$^{-1}$ and a stratified region with decreasing IGE abundances above 10,000 km s$^{-1}$. Based on a detailed comparative analysis, our conclusion is that the available pure deflagration models cannot fully explain either the observed properties of SN 2011ay or the results of our TARDIS modeling. Further examinations are necessary to find an adequate explanation for the origin of this object.

We present late-time optical $R$-band imaging data from the Palomar Transient Factory (PTF) for the nearby type Ia supernova SN 2011fe. The stacked PTF light curve provides densely sampled coverage down to $R\simeq22$ mag over 200 to 620 days past explosion. Combining with literature data, we estimate the pseudo-bolometric light curve for this event from 200 to 1600 days after explosion, and constrain the likely near-infrared contribution. This light curve shows a smooth decline consistent with radioactive decay, except over ~450 to ~600 days where the light curve appears to decrease faster than expected based on the radioactive isotopes presumed to be present, before flattening at around 600 days. We model the 200-1600d pseudo-bolometric light curve with the luminosity generated by the radioactive decay chains of $^{56}$Ni, $^{57}$Ni and $^{55}$Co, and find it is not consistent with models that have full positron trapping and no infrared catastrophe (IRC); some additional energy escape other than optical/near-IR photons is required. However, the light curve is consistent with models that allow for positron escape (reaching 75% by day 500) and/or an IRC (with 85% of the flux emerging in non-optical wavelengths by day 600). The presence of the $^{57}$Ni decay chain is robustly detected, but the $^{55}$Co decay chain is not formally required, with an upper mass limit estimated at 0.014 M$_{\odot}$. The measurement of the $^{57}$Ni/$^{56}$Ni mass ratio is subject to significant systematic uncertainties, but all of our fits require a high ratio >0.031 (>1.3 in solar abundances).

The double-detonation explosion model has been considered a candidate for explaining astrophysical transients with a wide range of luminosities. In this model, a carbon-oxygen white dwarf star explodes following detonation of a surface layer of helium. One potential signature of this explosion mechanism is the presence of unburned helium in the outer ejecta, left over from the surface helium layer. In this paper we present simple approximations to estimate the optical depths of important He I lines in the ejecta of double-detonation models. We use these approximations to compute synthetic spectra, including the He I lines, for double-detonation models obtained from hydrodynamical explosion simulations. Specifically, we focus on photospheric-phase predictions for the near-infrared 10830 \AA~and 2 $\mu$m lines of He I. We first consider a double detonation model with a luminosity corresponding roughly to normal SNe Ia. This model has a post-explosion unburned He mass of 0.03 $M_{\odot}$ and our calculations suggest that the 2 $\mu$m feature is expected to be very weak but that the 10830 \AA~feature may have modest opacity in the outer ejecta. Consequently, we suggest that a moderate-to-weak He I 10830 \AA~feature may be expected to form in double-detonation explosions at epochs around maximum light. However, the high velocities of unburned helium predicted by the model ($\sim 19,000$~km~s$^{-1}$) mean that the He I 10830 \AA~feature may be confused or blended with the C I 10690~\AA~line forming at lower velocities. We also present calculations for the He I 10830 \AA~and 2 $\mu$m lines for a lower mass (low luminosity) double detonation model, which has a post-explosion He mass of 0.077 $M_{\odot}$. In this case, both the He I features we consider are strong and can provide a clear observational signature of the double-detonation mechanism.

Due to its proximity, the Milky Way nuclear star cluster provides us with a wealth of data not available in other galactic nuclei. In particular, with adaptive optics, we can observe the detailed properties of individual stars, which can offer insight into the origin and evolution of the cluster. We summarize work on the central parsec of the Galactic center based on imaging and spectroscopic observations at the Keck and Gemini telescopes. These observations include stellar positions in two dimension and the velocity in three dimensions. Spectroscopic observations also enable measurements of the physical properties of individual stars, such as the spectral type and in some cases the effective temperature, metallicity, and surface gravity. We present a review of our latest measurements of the density profiles and luminosity functions of the young and old stars in this region. These observations show a complex stellar population with a young (4-6 Myr) compact star cluster in the central 0.5 pc embedded in an older and much more massive nuclear star cluster. Surprisingly, the old late-type giants do not show a cusp profile as long been expected from theoretical work. The solution to the missing cusp problem may offer us insight into the dynamical evolution of the cluster. Finally, we also discuss recent work on the metallicity of stars in this region and how they might be used to trace their origin. The nuclear star cluster shows a large variation in metallicity ([M/Fe]). The majority of the stars have higher than solar metallicity, with about 6\% having [M/Fe] $< -0.5$. These observations indicate that the NSC was not built from the globular clusters that we see today. The formation of the nuclear star cluster is more likely from the inward migration of gas originating in the disk of the Milky Way.

Knowing the metallicity distribution of stars in the Galactic Centre has important implications for the formation history of the Milky Way nuclear star cluster. However, this distribution is not well known, and is currently based on a small sample of fewer than 100 stars. We obtained near-infrared K-band spectra of more than 700 late-type stars in the central 4 pc^2 of the Milky Way nuclear star cluster with the integral-field spectrograph KMOS (VLT). We analyse the medium-resolution spectra using a full-spectral fitting method employing the G\"ottingen Spectral library of synthetic PHOENIX spectra. The derived stellar metallicities range from metal-rich [M/H]>+0.3 dex to metal-poor [M/H]<-1.0 dex, with a fraction of 5.2(^{+6.0}+{-3.1}) per cent metal-poor ([M/H]<-0.5 dex) stars. The metal-poor stars are distributed over the entire observed field. The origin of metal-poor stars remains unclear. They could originate from infalling globular clusters. For the metal-rich stellar population ([M/H]>0 dex) a globular cluster origin can be ruled out. As there is only a very low fraction of metal-poor stars in the central 4 pc^2 of the Galactic Centre, we believe that our data can discard a scenario in which the Milky Way nuclear star cluster is purely formed from infalling globular clusters.

Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we find that the background subtraction method adopted by the authors does not adequately remove contaminating stars in the small number Poisson limit. Hence, we conclude that their results do not provide evidence of young generations of stars within these clusters.

We present a metallicity analysis of 83 late-type giants within the central 1 pc of the Milky Way. K-band spectroscopy of these stars were obtained with the medium-spectral resolution integral-field spectrograph NIFS on Gemini North using laser-guide star adaptive optics. Using spectral template fitting with the MARCS synthetic spectral grid, we find that there is large variation in metallicity, with stars ranging from [M/H] $<$ -1.0 to above solar metallicity. About 6\% of the stars have [M/H] $<$ -0.5. This result is in contrast to previous observations, with smaller samples, that show stars at the Galactic center have approximately solar metallicity with only small variations. Our current measurement uncertainties are dominated by systematics in the model, especially at [M/H] $>$ 0, where there are stellar lines not represented in the model. However, the conclusion that there are low metallicity stars, as well as large variations in metallicity is robust. The metallicity may be an indicator of the origin of these stars. The low-metallicity population is consistent with that of globular clusters in the Milky Way, but their small fraction likely means that globular cluster infall is not the dominant mechanism for forming the Milky Way nuclear star cluster. The majority of stars are at or above solar metallicity, which suggests they were formed closer to the Galactic center or from the disk. In addition, our results indicate that it will be important for star formation history analyses using red giants at the Galactic center to consider the effect of varying metallicity.

In this letter we present an optical spectrum of SN 2011fe taken 1034 d after the explosion, several hundred days later than any other spectrum of a Type Ia supernova (disregarding light-echo spectra and local-group remnants). The spectrum is still dominated by broad emission features, with no trace of a light echo or interaction of the supernova ejecta with surrounding interstellar material. Comparing this extremely late spectrum to an earlier one taken 331 d after the explosion, we find that the most prominent feature at 331 d - [Fe III] emission around 4700 A - has entirely faded away, suggesting a significant change in the ionisation state. Instead, [Fe II] lines are probably responsible for most of the emission at 1034 d. An emission feature at 6300-6400 A has newly developed at 1034 d, which we tentatively identify with Fe I {\lambda}6359, [Fe I] {\lambda}{\lambda}6231, 6394 or [O I] {\lambda}{\lambda}6300, 6364. Interestingly, the features in the 1034-d spectrum seem to be collectively redshifted, a phenomenon that we currently have no convincing explanation for. We discuss the implications of our findings for explosion models, but conclude that sophisticated spectral modelling is required for any firm statement.

The Type Ia supernova SN 2011fe is one of the closest supernovae of the past decades. Due to its proximity and low dust extinction, this object provides a very rare opportunity to study the extremely late time evolution ($>900$ days) of thermonuclear supernovae. In this Letter, we present our photometric data of SN 2011fe taken at an unprecedented late epoch of $\approx 930$ days with GMOS-N mounted on the Gemini North telescope ($g=23.43 \pm 0.28$, $r=24.14 \pm 0.14$, $i=23.91 \pm 0.18$, and $z=23.90 \pm 0.17$) to study the energy production and retention in the ejecta of SN 2011fe. Together with previous measurements by other groups, our result suggests that the optical supernova light curve can still be explained by the full thermalization of the decay positrons of Co56. This is in spite of theoretical predicted effects (e.g. infrared catastrophe, positron escape, and dust) that advocate a substantial energy redistribution and/or loss via various processes that result in a more rapid dimming at these very late epochs.

The helium-rich hot subdwarf LS IV -14 116 shows remarkably high surface abundances of zirconium, yttrium, strontium, and germanium, indicative of strong chemical stratification in the photosphere. It also shows photometric behaviour indicative of non-radial g-mode pulsations, despite having surface properties inconsistent with any known pulsational instability zone. We have conducted a search for radial velocity variability. This has demonstrated that at least one photometric period is observable in several absorption lines as a radial velocity variation with a semi-amplitude in excess of 5 km s$^{-1}$. A correlation between line strength and pulsation amplitude provides evidence that the photosphere pulsates differentially. The ratio of light to velocity amplitude is too small to permit the largest amplitude oscillation to be radial.

We present TARDIS - an open-source code for rapid spectral modelling of supernovae (SNe). Our goal is to develop a tool that is sufficiently fast to allow exploration of the complex parameter spaces of models for SN ejecta. This can be used to analyse the growing number of high-quality SN spectra being obtained by transient surveys. The code uses Monte Carlo methods to obtain a self-consistent description of the plasma state and to compute a synthetic spectrum. It has a modular design to facilitate the implementation of a range of physical approximations that can be compared to asses both accuracy and computational expediency. This will allow users to choose a level of sophistication appropriate for their application. Here, we describe the operation of the code and make comparisons with alternative radiative transfer codes of differing levels of complexity (SYN++, PYTHON, and ARTIS). We then explore the consequence of adopting simple prescriptions for the calculation of atomic excitation, focussing on four species of relevance to Type Ia supernova spectra - Si II, S II, Mg II, and Ca II. We also investigate the influence of three methods for treating line interactions on our synthetic spectra and the need for accurate radiative rate estimates in our scheme.

The identity of Type Ia supernova progenitors remains a mystery, with various lines of evidence pointing towards either accretion from a non-degenerate companion, or the rapid merger of two degenerate stars leading to the thermonuclear destruction of a white dwarf. In this paper we spectroscopically scrutinize 24 of the brightest stars residing in the central 38" x 38" of the SN 1604 (Kepler) supernova remnant to search for a possible surviving companion star. We can rule out, with high certainty, a red giant companion star - a progenitor indicated by some models of the supernova remnant. Furthermore, we find no star that exhibits properties uniquely consistent with those expected of a donor star down to L>10Lsun. While the distribution of star properties towards the remnant are consistent with unrelated stars, we identify the most promising candidates for further astrometric and spectroscopic follow-up. Such a program would either discover the donor star, or place strong limits on progenitor systems to luminosities with L<<Lsun.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

In the SD scenario of SNe Ia the companion survives the SN explosion and thus should be visible near the center of the SN remnant and may show some unusual features. A promising approach to test progenitor models of SNe Ia is to search for the companion in SNRs. Here we present the results of 3D hydrodynamics simulations of the interaction between the SN Ia blast wave and a MS companion taking into consideration its orbital motion and spin. The primary goal of this work is to investigate the rotation of surviving companions after SN Ia explosions in the WD+MS scenario. We use Eggleton's code including the optically thick accretion wind model to obtain realistic models of companions. The impact of the SN blast wave on these companions is followed in 3D hydrodynamic simulations employing the SPH code GADGET3. We find that the rotation of the companion does not significantly affect the amount of stripped mass and the kick velocity caused by the SN impact. However, in our simulations, the rotational velocity of the companion is significantly reduced to about 14% to 32% of its pre-explosion value due to the expansion of the companion and the fact that 55%-89% of the initial angular momentum is carried away by the stripped matter. Compared with the observed rotational velocity of the presumed companion star of Tycho's SN, Tycho G, of 6 km/s the final rotational velocity we obtain is still higher by at least a factor of two. Whether this difference is significant, and may cast doubts on the suggestion that Tycho G is the companion of SN 1572, has to be investigated in future studies. Based on binary population synthesis results we present, for the first time, the expected distribution of rotational velocities of companions after the explosion which may provide useful information for the identification of the surviving companion in observational searches in other historical SNRs.

Located within 10 degrees of the North Pole, northern Ellesmere Island offers continuous darkness in the winter months. This capability can greatly enhance the detection efficiency of planetary transit surveys and other time domain astronomy programs. We deployed two wide-field cameras at 80 degrees North, near Eureka, Nunavut, for a 152-hour observing campaign in February 2012. The 16-megapixel-camera systems were based on commercial f/1.2 lenses with 70mm and 42mm apertures, and they continuously imaged 504 and 1,295 square degrees respectively. In total, the cameras took over 44,000 images and produced better-than-1% precision light curves for approximately 10,000 stars. We describe a new high-speed astrometric and photometric data reduction pipeline designed for the systems, test several methods for the precision flat-fielding of images from very-wide-angle cameras, and evaluate the cameras' image qualities. We achieved a scintillation-limited photometric precision of 1-2% in each 10s exposure. Binning the short exposures into 10-min chunks provided a photometric stability of 2-3 millimagnitudes, sufficient for the detection of transiting exoplanets around the bright stars targeted by our survey. We estimate that the cameras, when operated over the full arctic winter, will be capable of discovering several transiting exoplanets around bright (V<9) stars.

In this paper, we report on our analysis using Hubble Space Telescope astrometry and Keck-I HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With these data, we measured the proper motions, radial velocities, rotational velocities, and chemical abundances of these objects. Regarding the chemical abundances, we do not confirm the unusu- ally high [Ni/Fe] ratio previously reported for Tycho-G. Rather, we find that for all metrics in all stars, none exhibit the characteristics expected from traditional SN Ia single-degenerate-scenario calculations. The only possible exception is Tycho-B, a rare, metal-poor A-type star; however, we are unable to find a suitable scenario for it. Thus, we suggest that SN 1572 cannot be explained by the standard single-degenerate model.

Type Ia supernovae play a significant role in the evolution of the Universe and have a wide range of applications. It is widely believed that these events are the thermonuclear explosions of carbon-oxygen white dwarfs close to the Chandrasekhar mass (1.38 M\odot). However, CO white dwarfs are born with masses much below the Chandrasekhar limit and thus require mass accretion to become Type Ia supernovae. There are two main scenarios for accretion. First, the merger of two white dwarfs and, second, a stable mass accretion from a companion star. According to predictions, this companion star (also referred to as donor star) survives the explosion and thus should be visible in the center of Type Ia remnants. In this paper we scrutinize the central stars (79 in total) of the SN 1006 remnant to search for the surviving donor star as predicted by this scenario. We find no star consistent with the traditional accretion scenario in SN1006.

The element abundance distributions in the ejecta of Type Ia supernova (SN) is studied by modelling a time series of optical spectra of SN 2003du until ~1 year after the explosion. Since SN 2003du is a very normal Type Ia SN both photometrically and spectroscopically, the abundance distribution derived for it can be considered as representative of normal Type Ia SNe. We find that the innermost layers are dominated by stable Fe-group elements, with a total mass of ~ 0.2 Msun, which are synthesized through electron capture. Above the core of stable elements there are thick 56Ni-rich layers. The total mass of 56Ni is 0.65 Msun. The Si- and S-rich layers are located above the 56Ni-rich layers. The dominant element in the outermost layers (M_r > 1.1 Msun, v > 13000 km/s) is O, with a small amount of Si. Little unburned C remains, with an upper limit of 0.016 Msun. The element distributions in the ejecta are moderately mixed, but not fully mixed as seen in three-dimensional deflagration models.

It is widely believed that Type Ia supernovae (SN Ia) originate in binary systems where a white dwarf accretes material from a companion star until its mass approaches the Chandrasekhar mass and carbon is ignited in the white dwarf's core. This scenario predicts that the donor star should survive the supernova explosion, providing an opportunity to understand the progenitors of Type Ia supernovae.In this paper we argue that rotation is a generic signature expected of most non-giant donor stars that is easily measurable. \citep{2004Natur.431.1069R} examined stars in the center of the remnant of SN 1572 (Tycho's SN) and showed evidence that a subgiant star (Star G by their naming convention) near the remnant's centre was the system's donor star. We present high-resolution (R \simeq 40000) spectra taken with the High Dispersion Spectrograph on Subaru of this candidate donor star and measure the star's radial velocity as $79\pm 2$ \kms with respect to the LSR and put an upper limit on the star's rotation of 7.5 \kms. In addition, by comparing images that were taken in 1970 and 2004, we measure the proper motion of Star G to be $\mu_l = -1.6 \pm 2.1$ \masyr and $\mu_b = -2.7 \pm 1.6$ \masyr. We demonstrate that all of the measured properties of Star G presented in this paper are consistent with those of a star in the direction of Tycho's SN that is not associated with the supernova event. However, we discuss an unlikely, but still viable scenario for Star G to be the donor star, and suggest further observations that might be able to confirm or refute it.