Calcium-rich gap transients: Tidal detonations of white dwarfs?

We hypothesize that at least some of the recently discovered class of calcium-rich gap transients are tidal detonation events of white dwarfs (WDs) by black holes (BHs) or possibly neutron stars. We show that the properties of the calcium-rich gap transients agree well with the predictions of the tidal detonation model. Under the predictions of this model, we use a follow-up X-ray observation of one of these transients, SN 2012hn, to place weak upper limits on the detonator mass of this system that include all intermediate-mass BHs (IMBHs). As these transients are preferentially in the stellar haloes of galaxies, we discuss the possibility that these transients are tidal detonations of WDs caused by random flyby encounters with IMBHs in dwarf galaxies or globular clusters. This possibility has been already suggested in the literature but without connection to the calcium-rich gap transients. In order for the random flyby cross-section to be high enough, these events would have to be occurring inside these dense stellar associations. However, there is a lack of evidence for IMBHs in these systems, and recent observations have ruled out all but the very faintest dwarf galaxies and globular clusters for a few of these transients. Another possibility is that these are tidal detonations caused by three-body interactions, where a WD is perturbed towards the detonator in isolated multiple star systems. We highlight a number of ways this could occur, even in lower mass systems with stellar-mass BHs or neutron stars. Finally, we outline several new observational tests of this scenario, which are feasible with current instrumentation.

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Context: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are, however, found in multiple systems. 

Aims: By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary subpopulations with one another and with that of presumed-single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the spin rates of massive stars. 

Methods: We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (ν_esini) for components of 114 spectroscopic binaries in 30 Doradus. The ν_esini values are derived from the full width at half maximum (FWHM) of a set of spectral lines, using a FWHM vs. ν_esini calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. 

Results: The overall ν_esini distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at ν_esini<200 kms^-1) and a shoulder at intermediate velocities (200 <ν_esini<300 kms^-1). The distributions of binaries and single stars, however, differ in two ways. First, the main peak at ν_esini ~ 100kms^-1 is broader and slightly shifted towards higher spin rates in the binary distribution than that of the presumed-single stars. This shift is mostly due to short-period binaries (P_orb~<10 d). Second, the ν_esini distribution of primaries lacks a significant population of stars spinning faster than 300 kms^-1, while such a population is clearly present in the single-star sample. The ν_esini distribution of binaries with amplitudes of radial velocity variation in the range of 20 to 200 kms^-1 (mostly binaries with P_orb ~ 10-1000 d and/or with q<0.5) is similar to that of single O stars below ν_esini_~<170kms^-1. 

Conclusions: Our results are compatible with the assumption that binary components formed with the same spin distribution as single stars, and that this distribution contains few or no fast-spinning stars. The higher average spin rate of stars in short-period binaries may either be explained by spin-up through tides in such tight binary systems, or by spin-down of a fraction of the presumed-single stars and long-period binaries through magnetic braking (or by a combination of both mechanisms). Most primaries and secondaries of SB2 systems with P_orb~<10 d appear to have similar rotational velocities. This is in agreement with tidal locking in close binaries where the components have similar radii. The lack of very rapidly spinning stars among binary systems supports the idea that most stars with ν_esini_~> 300kms^-1 in the single-star sample are actually spun-up post-binary interaction products. Finally, the overall similarities (low-velocity peak and intermediate-velocity shoulder) of the spin distribution of binary and single stars argue for a massive star formation process in which the initial spin is set independently of whether stars are formed as single stars or as components of a binary system.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

The VLT-FLAMES tarantula survey:XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations - in terms of spectral type, luminosity class, and spatial location - in the field of view. Methods. We measured projected rotational velocities, 3e sin i, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(3e), of the equatorial rotational velocity, 3e. Results. The distribution of 3e sin i shows a two-component structure: a peak around 80 km s1 and a high-velocity tail extending up to 600 km s^-1 This structure is also present in the inferred distribution P(3e) with around 80% of the sample having 0 <3e ≤ 300 km s^-1 and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the highvelocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts. © 2013 Author(s).

Spectral analysis of four 'hypervariable' AGN: a micro-needle in the haystack?

We analyze four extreme AGN transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 magnitudes in the g-band on a timescale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein timescale of 7.5/10.8 years respectively. In one case the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects there is spectroscopic evidence for an intervening absorber but the extra structure seen in the lightcurves requires a more complex lensing scenario to adequately explain.

Photoionization of Co+ and electron-impact excitation of Co2 + using the Dirac R-matrix method

Modelling of massive stars and supernovae (SNe) plays a crucial role in understanding galaxies. From this modelling we can derive fundamental constraints on stellar evolution, mass-loss processes, mixing, and the products of nucleosynthesis. Proper account must be taken of all important processes that populate and depopulate the levels (collisional excitation, de-excitation, ionization, recombination, photoionization, bound–bound processes). For the analysis of Type Ia SNe and core collapse SNe (Types Ib, Ic and II) Fe group elements are particularly important. Unfortunately little data is currently available and most noticeably absent are the photoionization cross-sections for the Fe-peaks which have high abundances in SNe. Important interactions for both photoionization and electron-impact excitation are calculated using the relativistic Dirac atomic R-matrix codes (DARC) for low-ionization stages of Cobalt. All results are calculated up to photon energies of 45 eV and electron energies up to 20 eV. The wavefunction representation of Co III has been generated using GRASP0 by including the dominant 3d⁷, 3d⁶[4s, 4p], 3p⁴3d⁹ and 3p⁶3d⁹ configurations, resulting in 292 fine structure levels. Electron-impact collision strengths and Maxwellian averaged effective collision strengths across a wide range of astrophysically relevant temperatures are computed for Co III. In addition, statistically weighted level-resolved ground and metastable photoionization cross-sections are presented for Co II and compared directly with existing work.

The microphysics of collisionless shock waves

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebulæ, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in situ observations, analytical and numerical developments. A particular emphasis is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics.

A sub-kiloparsec scale view of star formation in M31

Thesis (Ph.D.)--University of Washington, 2016-08

Resolved Stars' Insights into Galaxy Physics

Thesis (Ph.D.)--University of Washington, 2016-08

Histoire évolutive de la galaxie spirale NGC 5430

La galaxie spirale barrée NGC 5430 est particulière en ce sens qu’elle présente un noeud Wolf-Rayet très lumineux et des bras asymétriques. Des spectres longue-fente le long de la barre et dans le bras déformé ainsi que des données SpIOMM couvrant l’ensemble de la galaxie ont été analysées. L’absorption stellaire sous-jacente a été soustraite des spectres longue-fente à l’aide d’un ajustement de modèles théoriques de populations stellaires fait avec le programme GANDALF. L’absorption a un impact très important sur le calcul de l’extinction ainsi que sur les différents diagnostics propres aux régions HII et aux populations stellaires jeunes. Enfin, cette étude montre que NGC 5430 comporte une composante gazeuse ionisée diffuse sur toute son étendue et qu’il est important d’en tenir compte afin d’appliquer correctement les diagnostics. Un des scénarios évolutifs proposés au terme de cette étude est que le noeud Wolf-Rayet constitue le restant d’une petite galaxie ou d’un nuage intergalactique qui serait entré en collision avec NGC 5430. Une structure englobant le noeud Wolf-Rayet se déplace à une vitesse considérablement inférieure (50 - 70 km s-1) à celle attendue à une telle distance du centre de la galaxie (200 - 220 km s-1). De plus, le noeud Wolf-Rayet semble très massif puisque l’intensité maximale du continu stellaire de cette région est semblable à celle du noyau et est de loin supérieure à celle de l’autre côté de la barre. Le nombre d’étoiles Wolf-Rayet (2150) est aussi considérable. Il n’est toutefois pas exclu que la différence de vitesses observée témoigne d’un écoulement de gaz le long de la barre, qui alimenterait la formation stellaire du noeud Wolf-Rayet ou du noyau.

Traitement de données observationnelles et développement d’un code de simulation pour SpIOMM et SITELLE

L’arrivée du spectromètre imageur à transformée de Fourier SITELLE au télescope Canada-France-Hawaï souligne la nécessité d’un calculateur de temps d’exposition permettant aux utilisateurs de l’instrument de planifier leurs observations et leurs demandes de temps de télescope. Une grande partie de mon projet est ainsi le développement d’un code de simulation capable de reproduire les résultats de SITELLE et de son prédecesseur SpIOMM, installé à l’Observatoire du Mont-Mégantic. La précision des simulations est confirmée par une comparaison avec des données SpIOMM et les premières observations de SITELLE. La seconde partie de mon projet consiste en une analyse spectrale de données observationelles. Prenant avantage du grand champ de vue de SpIOMM, les caractéristiques du gaz ionisé (vitesse radiale et intensité) sont étudiées pour l’ensemble de la paire de galaxies en interaction Arp 72. La courbe de rotation dans le visible ainsi que le gradient de métallicité de NGC 5996, la galaxie principale d’Arp 72, sont obtenues ici pour la première fois. La galaxie spirale NGC 7320 est également étudiée à partir d’observations faites à la fois avec SpIOMM et SITELLE.

Aperture effects on the oxygen abundance determinations from Califa data

This paper aims to provide aperture corrections for emission lines in a sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey (CALIFA) database. In particular, we explore the behavior of the log([O III] λ5007/Hβ)/([N II] λ6583/Hα) (O3N2) and log[N II] lambda 6583/Hα (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star-forming galaxies. We compute the median growth curves of Hα, Hα/Hβ, O3N2, and N-2 up to 2.5R(50) and 1.5 disk R-eff. These distances cover most of the optical spatial extent of the CALIFA galaxies. The growth curves simulate the effect of observing galaxies through apertures of varying radii. We split these growth curves by morphological types and stellar masses to check if there is any dependence on these properties. The median growth curve of the Hα flux shows a monotonous increase with radius with no strong dependence on galaxy inclination, morphological type, and stellar mass. The median growth curve of the Hα/HβH ratio monotonically decreases from the center toward larger radii, showing for small apertures a maximum value of ≈10% larger than the integrated one. It does not show any dependence on inclination, morphological type, and stellar mass. The median growth curve of N-2 shows a similar behavior, decreasing from the center toward larger radii. No strong dependence is seen on the inclination, morphological type, and stellar mass. Finally, the median growth curve of O3N2 increases monotonically with radius, and it does not show dependence on the inclination. However, at small radii it shows systematically higher values for galaxies of earlier morphological types and for high stellar mass galaxies. Applying our aperture corrections to a sample of galaxies from the SDSS survey at 0.02 ≤ z ≤ 0.3 shows that the average difference between fiber-based and aperture-corrected oxygen abundances, for different galaxy stellar mass and redshift ranges, reaches typically to ≈11%, depending on the abundance calibration used. This average difference is found to be systematically biased, though still within the typical uncertainties of oxygen abundances derived from empirical calibrations. Caution must be exercised when using observations of galaxies for small radii (e.g., below 0.5 R_eff) given the high dispersion shown around the median growth curves. Thus, the application of these median aperture corrections to derive abundances for individual galaxies is not recommended when their fluxes come from radii much smaller than either R_50 or R_eff.

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, 6 x 10(19) eV. The anisotropy was measured by the fraction of arrival directions that are less than 3.1 degrees from the position of an active galactic nucleus within 75 Mpc (using the Veron-Cetty and Veron 12th catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is (38(-6)(+7))%, compared with 21% expected for isotropic cosmic rays. This is down from the early estimate of (69-(+11)(13))%. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation. (C) 2010 Elsevier B.V. All rights reserved.

Near-infrared Instrumentation For Rapid-response Astronomy

Ɣ-ray bursts (GRBs) are the Universe's most luminous transient events. Since the discovery of GRBs was announced in 1973, efforts have been ongoing to obtain data over a broader range of the electromagnetic spectrum at the earliest possible times following the initial detection. The discovery of the theorized ``afterglow'' emission in radio through X-ray bands in the late 1990s confirmed the cosmological nature of these events. At present, GRB afterglows are among the best probes of the early Universe (z ≳ 9). In addition to informing theories about GRBs themselves, observations of afterglows probe the circum-burst medium (CBM), properties of the host galaxies and the progress of cosmic reionization. To explore the early-time variability of afterglows, I have developed a generalized analysis framework which models near-infrared (NIR), optical, ultra-violet (UV) and X-ray light curves without assuming an underlying model. These fits are then used to construct the spectral energy distribution (SED) of afterglows at arbitrary times within the observed window. Physical models are then used to explore the evolution of the SED parameter space with time. I demonstrate that this framework produces evidence of the photodestruction of dust in the CBM of GRB 120119A, similar to the findings from a previous study of this afterglow. The framework is additionally applied to the afterglows of GRB 140419A and GRB 080607. In these cases the evolution of the SEDs appears consistent with the standard fireball model. Having introduced the scientific motivations for early-time observations, I introduce the Rapid Infrared Imager-Spectrometer (RIMAS). Once commissioned on the 4.3 meter Discovery Channel Telescope (DCT), RIMAS will be used to study the afterglows of GRBs through photometric and spectroscopic observations beginning within minutes of the initial burst. The instrument will operate in the NIR, from 0.97 μm to 2.37 μm, permitting the detection of very high redshift (z ≳ 7) afterglows which are attenuated at shorter wavelengths by Lyman-ɑ absorption in the intergalactic medium (IGM). A majority of my graduate work has been spent designing and aligning RIMAS's cryogenic (~80 K) optical systems. Design efforts have included an original camera used to image the field surrounding spectroscopic slits, tolerancing and optimizing all of the instrument's optics, thermal modeling of optomechanical systems, and modeling the diffraction efficiencies for some of the dispersive elements. To align the cryogenic optics, I developed a procedure that was successfully used for a majority of the instrument's sub-assemblies. My work on this cryogenic instrument has necessitated experimental and computational projects to design and validate designs of several subsystems. Two of these projects describe simple and effective measurements of optomechanical components in vacuum and at cryogenic temperatures using an 8-bit CCD camera. Models of heat transfer via electrical harnesses used to provide current to motors located within the cryostat are also presented.

The telescopic eye: aqua optica

The Telescopic Eye: Aqua Optica was an installation work which explores concepts and philosophies related to instrumental, lens based observations of heavens as experienced when seeing through a telescope. In particular it considers how the eye of the observer and the lens of the telescope must act in unison to extend both the perception and conception of that which is remote and beyond direct naked eye sensing - an act of transcendence. A series of projected images (the video) of celestial objects are observed by the viewer through a large spherical aquatic lens which is mounted on a large wooden tripod. These images slowly dissolve and morph into each other. The sequence in which these images are presented is based on approximate chronological order in which these celestial objects were historically observed through the telescope – such as the local bodies of our the solar system starting with the Moon, Jupiter and Saturn (as first observed by Galileo in 1610) and eventually extending out to distant nebula, star clusters and galaxies. For the Transduction exhibition the Telescopic Eye: Aqua Optica was installed outside in a small quadrangle next to the Tech Bar and as the night drew darker it was able to also embody the lighting and other visual elements from the surrounding buildings and streets. The Aqua Optica system often encompasses such local lighting and images as a added layer to the video projection work. In this case (at night and outside at Federation Square) the extent of the layering was very successful and increased dramatically as the night became darker. This multi layering of images invokes (for the viewer) a connection between the immediate environment of the Federation square and cityscape with the vastly distant images of the cosmos. For me such a multilayered imaging represents how scientific instruments of observation can generate a complex and augmented visualisation of reality, one which has historically come to re-define our conception of the cosmos and our place in it both spatially and temporally.