AE Aur Research Articles

Spectral investigation of the interstellar cloud towards a runaway star AE Aur

Abstract The goal of this study is tracking of possible changes in the spectrum of the interstellar cloud seen along the sightline towards a runaway star AE Aur. The study considers variations in the radial velocities of the interstellar molecules and atoms, dwelling in the cloud, for the period ≈17 years. The most recent results of 6 year monitoring of equivalent widths variability for the CH and CH+ absorption lines are also presented. The study was done using the CH B-X(0,0) 3886 Å, A-X (0,0) 4300 Å and the CH+ A-X (1-0) 3958 Å, A-X (0,0) 4232 Å lines as well as the interstellar potassium K i 7699 Å line. Moreover a short discussion concerning the distance to the star is included. The present study is a continuation of our previous works and works of other researchers dedicated to the problem of small-scale structure in diffuse interstellar clouds. High resolution optical spectra from different spectrographs have been used in this work. The results of the study do not demonstrate significant changes in the radial velocities of the spectral features mentioned. The intensities of the CH and CH+ lines show slight tendency to vary at 1–3% level.

The origin of early-type runaway stars from open clusters

Runaway stars are ejected from their place of birth in the Galactic disk, with some young B-type runaways found several tens of kiloparsecs from the plane traveling at speeds beyond the escape velocity, which calls for violent ejection processes. Young open clusters are a likely place of origin, and ejection may be either through N-body interactions or in binary supernova explosions. The most energetic events may require dynamical interaction with massive black holes. The excellent quality of Gaia astrometry opens up the path to study the kinematics of young runaway stars to such a high precision that the place of origin in open stellar clusters can be identified uniquely even when the star is a few kiloparsecs away. We developed an efficient minimization method to calculate whether two or more objects may come from the same place, which we tested against samples of Orion runaways. Our fitting procedure was then used to calculate trajectories for known runaway stars, which have previously been studied from HIPPARCOS astrometry as well as known open clusters. For runaways in our sample we used Gaia data and updated radial velocities, and found that only half of the sample could be classified as runaways. The other half of the sample moves so slowly (< 30 km s−1) that they have to be considered as walkaway stars. Most of the latter stars turned out to be binaries. We identified parent clusters for runaways based on their trajectories. We then used cluster age and flight time of the stars to investigate whether the ejection was likely due to a binary supernova or due to a dynamical ejection. In particular we show that the classical runaways AE Aurigae and μ Columbae might not have originated together, with μ Columbae having an earlier ejection from Collinder 69, a cluster near the ONC. The second sample investigated comprises a set of distant runaway B stars in the halo which have been studied carefully by quantitative spectral analyses. We are able to identify candidate parent clusters for at least four stars including the hyper-runaway candidate HIP 60350. The ejection events had to be very violent, ejecting stars at velocities as large as 150–400 km s−1.

Radial velocity measurements of the runaway stars AE Aur and ι Ori

We present the results of our spectroscopic monitoring observations of the late O‐type runaway stars AE Aur and ι Ori, carried out between September 2016 and May 2018 using the Échelle spectrograph FLECHAS at the University Observatory Jena. For both stars, we obtained current radial velocity (RV) measurements with better precision than before. These results can be used for further advanced calculations of their space motion, needed to determine, for example, their place of origin in our Galaxy. In this paper, we describe the observations and the data reduction, as well as the analysis of all spectral data. We find a stable RV for AE Aur and confirm ι Ori as a double‐lined spectroscopic binary and derive its current Keplerian orbital elements.

Rise of the intensity of interstellar CH+ in the spectrum of AE Aur

AbstractUsing many spectra of the star AE Aur (HD 34078), covering the period 1997–2019, we demonstrate temporal variations of the intensities of spectral features carried by simple interstellar radicals. The drop of the CH+ and CH line intensities was demonstrated by Krełowski et al. (2016)—Paper I. More recent data suggest that the CH+ interstellar lines reach the minimum of intensity in 2015–2016 being then followed by the evident rise, which continues. The intensities of CH lines show a similar drop to CH+. However, the further (after 2016) rise seems marginal if any. Thus, the CH+ features grow since 2016 while those of CH remain almost constant, likely due to the fact that both species are formed in different environments. No variations in the radial velocities of the observed interstellar features were detected during the whole covered time period.

Runaway O-star Bow Shocks as Particle Accelerators? The Case of AE Aur Revisited

Abstract We present results of our Chandra/ACIS observations of the field centered on the fast, runaway O star AE Aur and its bow shock. Previous XMM-Newton observations revealed an X-ray “blob” near the IR arc tracing the bow shock, possibly a nonthermal source consistent with models of Inverse Compton scattering of dust IR photons by electrons accelerated at the shock. The new, subarcsecond-resolution Chandra data, while confirming the presence of the XMM-Newton source, clearly indicate that the latter is neither extended nor coincident with the IR arc and strongly suggest it is a background active galactic nucleus. Motivated by results published for the bow shock of BD+43°3654, we extended our study to the radio domain by analyzing archival EVLA data. We find no radio emission from the AE Aur bow shock either. The corresponding upper limits for the absorbed (unabsorbed) X-ray flux of 5.9(7.8) × 10−15 erg cm−2 s−1 (3σ) and, in the radio range of 2 mJy (1.4 GHz) and 0.4 mJy (5.0 GHz), are used to put constraints on model predictions for particle acceleration within the bow shock. In the “classical” framework of diffusive shock acceleration, we find that the predicted X-ray and radio emission by the bow shock is at least two orders of magnitude below the current upper limits, consistent with the systematic nondetections of up to 60 stellar bow shocks. The only exception so far remains that of BD+43°3654, which is probably the result of its very large mass-loss rate among runaway O stars.

Evidence for feedback and stellar-dynamically regulated bursty star cluster formation: the case of the Orion Nebula Cluster

A scenario for the formation of multiple co-eval populations separated in age by about 1 Myr in very young clusters (VYCs, ages less than 10 Myr) and with masses in the range 600–20 000 M⊙ is outlined. It rests upon a converging inflow of molecular gas building up a first population of pre-main sequence stars. The associated just-formed O stars ionise the inflow and suppress star formation in the embedded cluster. However, they typically eject each other out of the embedded cluster within 106 yr, that is before the molecular cloud filament can be ionised entirely. The inflow of molecular gas can then resume forming a second population. This sequence of events can be repeated maximally over the life-time of the molecular cloud (about 10 Myr), but is not likely to be possible in VYCs with mass <300 M⊙, because such populations are not likely to contain an O star. Stellar populations heavier than about 2000 M⊙ are likely to have too many O stars for all of these to eject each other from the embedded cluster before they disperse their natal cloud. VYCs with masses in the range 600–2000 M⊙ are likely to have such multi-age populations, while VYCs with masses in the range 2000–20 000 M⊙ can also be composed solely of co-eval, mono-age populations. More massive VYCs are not likely to host sub-populations with age differences of about 1 Myr. This model is applied to the Orion Nebula Cluster (ONC), in which three well-separated pre-main sequences in the colour–magnitude diagram of the cluster have recently been discovered. The mass-inflow history is constrained using this model and the number of OB stars ejected from each population are estimated for verification using Gaia data. As a further consequence of the proposed model, the three runaway O star systems, AE Aur, μ Col and ι Ori, are considered as significant observational evidence for stellar-dynamical ejections of massive stars from the oldest population in the ONC. Evidence for stellar-dynamical ejections of massive stars in the currently forming population is also discussed.

26. Note on the variable or new nebulae possibly connected with FU orionis, AE aurigae and nova aurigae 1891

26. Note on the variable or new nebulae possibly connected with FU orionis, AE aurigae and nova aurigae 1891

Formation of runaway stars in a star-cluster potential

We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $\mu$ Col and binary $\iota$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $\mu$ Col.

An analytic method for identifying dynamically formed runaway stars

In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $\mu$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.

On the Absence of Non-thermal X-Ray Emission around Runaway O Stars

Abstract Theoretical models predict that the compressed interstellar medium around runaway O stars can produce high-energy non-thermal diffuse emission, in particular, non-thermal X-ray and γ-ray emission. So far, detection of non-thermal X-ray emission was claimed for only one runaway star, AE Aur. We present a search for non-thermal diffuse X-ray emission from bow shocks using archived XMM-Newton observations for a clean sample of six well-determined runaway O stars. We find that none of these objects present diffuse X-ray emission associated with their bow shocks, similarly to previous X-ray studies toward ζ Oph and BD+43°3654. We carefully investigated multi-wavelength observations of AE Aur and could not confirm previous findings of non-thermal X-rays. We conclude that so far there is no clear evidence of non-thermal extended emission in bow shocks around runaway O stars.

Modeling nonthermal emission from stellar bow shocks

Runaway O- and early B-type stars passing throughout the interstellar medium at supersonic velocities and characterized by strong stellar winds may produce bow shocks that can serve as particle acceleration sites. Previous theoretical models predict the production of high energy photons by non-thermal radiative processes, but their efficiency is still debated. We aim to test and explain the possibility of emission from the bow shocks formed by runaway stars traveling through the interstellar medium by using previous theoretical models. We apply our model to AE Aurigae, the first reported star with an X-ray detected bow shock, to BD+43 3654, in which the observations failed in detecting high energy emission, and to the transition phase of a supergiant star in the late stages of its life.From our analysis, we confirm that the X-ray emission from the bow shock produced by AE Aurigae can be explained by inverse Compton processes involving the infrared photons of the heated dust. We also predict low high energy flux emission from the bow shock produced by BD+43 3654, and the possibility of high energy emission from the bow shock formed by a supergiant star during the transition phase from blue to red supergiant.Bow shock formed by different type of runaway stars are revealed as a new possible source of high energy photons in our neighbourhood.

Position Displacement of Diffuse Interstellar Bands 1 1This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory (Chile).

We reconsider the already published phenomenon of the blue shift of diffuse interstellar bands, observed in spectra of HD34078 (AE Aur) and members of the Sco OB1 association, in particular HD152233. We have analyzed 29 diffuse bands. Some of them, already proven as blue-shifted in our earlier study, are now confirmed using another instrument: the 6.5 m Clay telescope equipped with the MIKE spectrograph. The high signal-to-noise ratio (over 600) of our spectra allowed us to reveal even small small-scale displacements of positions (both blue and redshifts) of diffuse bands along the considered lines of sight. In some cases, the magnitude of deviation exceeds 10 km s-1. Also, we prove that profiles of many diffuse bands in spectra of HD34078 suffer significant broadening. The origin of the observed phenomena is discussed.

Systematic search for high-energy gamma-ray emission from bow shocks of runaway stars

Context. It has been suggested that the bow shocks of runaway stars are sources of high-energy gamma rays (E > 100 MeV). Theoretical models predicting high-energy gamma-ray emission from these sources were followed by the first detection of non-thermal radio emission from the bow shock of BD+43$^\deg$ 3654 and non-thermal X-ray emission from the bow shock of AE Aurigae. Aims. We perform the first systematic search for MeV and GeV emission from 27 bow shocks of runaway stars using data collected by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (Fermi). Methods. We analysed 57 months of Fermi-LAT data at the positions of 27 bow shocks of runaway stars extracted from the Extensive stellar BOw Shock Survey catalogue (E-BOSS). A likelihood analysis was performed to search for gamma-ray emission that is not compatible with diffuse background or emission from neighbouring sources and that could be associated with the bow shocks. Results. None of the bow shock candidates is detected significantly in the Fermi-LAT energy range. We therefore present upper limits on the high-energy emission in the energy range from 100 MeV to 300 GeV for 27 bow shocks of runaway stars in four energy bands. For the three cases where models of the high-energy emission are published we compare our upper limits to the modelled spectra. Our limits exclude the model predictions for Zeta Ophiuchi by a factor $\approx$ 5.

AE AURIGAE: FIRST DETECTION OF NON-THERMAL X-RAY EMISSION FROM A BOW SHOCK PRODUCED BY A RUNAWAY STAR

Runaway stars produce shocks when passing through interstellar medium at supersonic velocities. Bow shocks have been detected in the mid-infrared for several high-mass runaway stars and in radio waves for one star. Theoretical models predict the production of high-energy photons by non-thermal radiative processes in a number sufficiently large to be detected in X-rays. To date, no stellar bow shock has been detected at such energies. We present the first detection of X-ray emission from a bow shock produced by a runaway star. The star is AE Aur, which was likely expelled from its birthplace by the encounter of two massive binary systems and now is passing through the dense nebula IC 405. The X-ray emission from the bow shock is detected at 30" to the northeast of the star, coinciding with an enhancement in the density of the nebula. From the analysis of the observed X-ray spectrum of the source and our theoretical emission model, we confirm that the X-ray emission is produced mainly by inverse Compton upscattering of infrared photons from dust in the shock front.

RUBIDIUM IN THE INTERSTELLAR MEDIUM

We present observations of interstellar rubidium toward o Per, zeta Per, AE Aur, HD 147889, chi Oph, zeta Oph, and 20 Aql. Theory suggests that stable 85Rb and long-lived 87Rb are produced predominantly by high-mass stars, through a combination of the weak s- and r-processes. The 85Rb/87Rb ratio was determined from measurements of the Rb I line at 7800 angstroms and was compared to the solar system meteoritic ratio of 2.59. Within 1-sigma uncertainties all directions except HD 147889 have Rb isotope ratios consistent with the solar system value. The ratio toward HD 147889 is much lower than the meteoritic value and similar to that toward rho Oph A (Federman et al. 2004); both lines of sight probe the Rho Ophiuchus Molecular Cloud. The earlier result was attributed to a deficit of r-processed 85Rb. Our larger sample suggests instead that 87Rb is enhanced in these two lines of sight. When the total elemental abundance of Rb is compared to the K elemental abundance, the interstellar Rb/K ratio is significantly lower than the meteoritic ratio for all the sight lines in this study. Available interstellar samples for other s- and r- process elements are used to help interpret these results.

Blueshifted diffuse interstellar bands in the spectrum of HD 34078

In this paper, we report the very first observation of diffuse interstellar bands (DIBs) that, in the spectrum of HD 34078 (AE Aur), are blueshifted with respect to the normal position that they have in other objects, where the rest-wavelength velocity frame is determined using very sharp interstellar atomic lines or molecular features. Only reasonably narrow DIBs seemingly show this effect, which is absent in broader ones. The result is confirmed independently using three different spectrographs attached to two different telescopes.

The Massive Close Binary in the δ Ori A Triple System

AbstractWe present an analysis of short-wave, high-dispersion ultraviolet spectra of the triple star δ Ori A from the International Ultraviolet Explorer Satellite’s (IUE) Final Archive. These spectra were cross-correlated against AE Aur to find the components’ radial velocities, which were then used to produce the system’s orbital elements. The long-period tertiary star in the δ Ori A system was not seen in the resulting cross-correlation functions (CCFs). The close binary’s eclipses allow the orbit’s inclination to be estimated by modeling of its Hipparcos light curve. The primary star appears to have a mass of 11.2 M⊙ and the secondary seems to have a mass of 5.6 M⊙, both of which are about 1/3 of the expected values for stars of their MK types. Although we expected the massive close binary in the δ Ori A system to be a pre-Roche lobe overflow (RLOF) system, these masses appear to require that it be a post-RLOF system. The full description of this work, including the tomographic separation of the spectra for the close binary’s components, appears in Harvin et al. (2002).

ι Orionis—Evidence for a Capture Origin Binary

The star ι Ori (HD 37043, HR 1899) is a 29.1338 day period, double-lined O-type spectroscopic binary with a highly eccentric orbit. The classification of the primary (O9 III) is well established. The secondary, roughly 2 mag fainter, is generally acknowledged as an early-type B star, but claims about its luminosity class have been made only from its estimated light contribution and not from its spectrum. We have collected visual data with the 1 m Multiple-Telescope Telescope (MTT; at the Center for High Angular Resolution at Georgia State University) in four wave bands well distributed in orbital phase, with the periastron passage particularly well covered. We present reconstructed component spectra from our tomographic reconstruction technique. According to three Morgan-Keenan (MK) luminosity criteria, we can unambiguously characterize the secondary as a B1 III-IV giant (MK type or B0.8 III-IV interpolated type). We also have constructed three line-strength indices from eighteen observed individual lines for the secondary and a grid of comparison stars. These results confirm the typing from the MK luminosity criteria. Furthermore, we argue, from the mass ratio (q = 0.5), evolutionary considerations, the orbital scale, and the unremarkable CNO element equivalent widths of the secondary, that this system is not coevolved nor has it undergone significant mass transfer. We postulate an origin for the system in a binary-binary collision and suspect that the runaway stars μ Col and AE Aur may have been the original partners for the current ι Ori secondary and primary, respectively.

The Origin of Runaway Stars

Milliarcsecond astrometry provided by Hipparcos and by radio observations makes it possible to retrace the orbits of some of the nearest runaway stars and pulsars to determine their site of origin. The orbits of the runaways AE Aurigae and μ Columbae and of the eccentric binary ι Orionis intersected each other ~2.5 Myr ago in the nascent Trapezium cluster, confirming that these runaways were formed in a binary-binary encounter. The path of the runaway star ζ Ophiuchi intersected that of the nearby pulsar PSR J1932+1059, ~1 Myr ago, in the young stellar group Upper Scorpius. We propose that this neutron star is the remnant of a supernova that occurred in a binary system that also contained ζ Oph and deduce that the pulsar received a kick velocity of ~350 km s-1 in the explosion. These two cases provide the first specific kinematic evidence that both mechanisms proposed for the production of runaway stars, the dynamical ejection scenario and the binary-supernova scenario, operate in nature.

Examination of the Interstellar Spectrum of AE Aur for Long‐Term Changes

ABSTRACTThe runaway O9.5 V star AE Aur is moving across the line of sight behind the curtain provided by the emission/reflection nebula IC 405. The interstellar line spectrum of AE Aur was examined on a HIRES spectrogram obtained in 1998 to determine if there had been any change since the original high‐resolution observations by Adams in the 1940s. No equivalent widths were published by Adams, so only a comparison of radial velocities was possible. No significant changes over the intervening half‐century were detected; during that time the cross motion of AE Aur in the plane of the sky amounted to about 1450 AU. Similarly, no convincing changes were found in the interstellar Ca ii lines of the runaway O star μ Col.