Microscopic processes in global relativistic jets containing helical magnetic fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron–proton (e−–p+) and electron–positron (e±) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of “global” jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the e−–p+ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the e± jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

(Table T1) Trace-element abundances for igneous basement at ODP Site 206-1256

In this report, I present trace element data for basement samples at Ocean Drilling Program (ODP) Site 1256. The samples analyzed represent a subset of the group ("pool") samples from ODP Leg 206, and these trace element data are part of a more comprehensive data suite for the same samples, with analyses of stable and radiogenic isotopes (e.g., Sr, Li, and O) in progress or recently completed that will be presented elsewhere. The trace element analyses were performed in the GeoAnalytical Lab at Washington State University. The following elements were analyzed: La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Ba, Th, Nb, Y, Hf, Ta, U, Pb, Rb, Cs, Sr, Sc, and Zr. Trace element data indicate that the igneous basement at Site 1256 is geochemically normal mid-ocean-ridge basalt. A massive ponded flow sampled in both Holes 1256C and 1256D is distinguished by higher abundances of rare earth elements (REE) and most of the other trace elements analyzed. One interval of highly altered basalt has significantly higher concentrations of Cs, Rb, and Ba and lower concentrations of Sr, Pb, Zr, Hf, Sc, and most REE than the samples of background alteration or halos. No correlation is obvious between trace element abundance and macroscopic type of alteration within the background alteration or halos.

(Table 1) Sulfur and oxygen isotopes from two anhydrite samples recovered from DSDP Hole 70-504B

The first anhydrite reported from oceanic basalts occurs in altered basalts drilled during DSDP Leg 70 from Hole 504B. Anhydrite has been identified in several samples, two of which were studied in detail. Anhydrite in Sample 504B-40-3 (130-135 cm), which was acquired at 310 meters sub-basement, occurs in a dolerite at the center of a vug rimmed by saponite and calcite. Red iron-hydroxide-rich alteration halos occur from 0 to 310 meters sub-basement; primary sulfides in these halos are oxidized, and the rocks have lost large amounts of sulfur. The anhydrite in this sample has a d34S value of 18.5 per mil, and it is interpreted to have formed from a fluid containing a mixture of seawater sulfate (20.9 per mil) and basaltic sulfur (0 per mil) released during the oxidation of primary sulfides. Anhydrite in Sample 504B-48-3 (14-18 cm), which was found at 376 meters sub-basement, occurs intergrown with gyrolite at the center of a 1-cm-wide vein that is rimmed by saponite and quartz. At sub-basement depths below 310 meters to the bottom of the Leg 70 section (562 m sub-basement), the rocks exhibit the effects of anoxic alteration with common secondary pyrite. Anhydrite in Sample 504B-48-3 (14-18 cm) has a d34S value of 36.7 per mil, and it is interpreted to have formed from seawater-derived fluids enriched in 34S through sulfate reduction. Temperatures of alteration calculated from oxygen isotope data range from 60 to 100°C. Sulfate reduction may have occurred in situ, or elsewhere at higher temperature, possibly deeper in the crust. The secondary mineral paragenetic sequence indicates a progressive decrease in Mg and increase in Ca in the circulating fluids. This eventually led to anhydrite formation late in the alteration process.

Geochemistry of sheeted dike complex minerals of ODP Hole 140-504B

Diabases were recovered during Legs 137 and 140 at Hole 504B from depths between 1621.5 and 2000.4 meters below seafloor in the lower sheeted dike complex. The samples contain multiple generations of millimetric to centimetric veins. The orientation of the measured veins suggests that two main vein sets exist: one characterized by shallow dipping and the other by random trend. Thermal contraction during rock cooling is considered the main mechanism responsible for fracture formation. Vein infill is related to the circulation of hydrothermal fluids near the spreading axis. Some veins are surrounded by millimeter-sized alteration halos due to fluid percolation from the fractures through the host rock. Vein-filling minerals are essentially amphibole, chlorite, and zeolites. Amphibole composition is controlled by the microstructural site of the rock. Actinolite is the main amphibole occurring in the veins and also in the groundmass away from the halos. In the alteration halos, amphibole shows composition of actinolitic hornblende and Mg-hornblende. Late-stage tension gashes and interstitial spaces in some amphibole-bearing veins are filled with zeolites, suggesting that the veins likely suffered multiple opening stages that record the cooling history of the circulating fluids. Evidence of deformation recorded by the recovered samples seems to be restricted to veins that clearly represent elements of weakness of the rock. On the basis of vein geometry and microstructure we infer structural interpretations for the formation mechanism and for deformation of veins.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

Supernova 2013fc in a circumnuclear ring of a luminous infrared galaxy: The big brother of SN 1998S

We present photometric and spectroscopic observations of SN 2013fc, a bright type II supernova (SN) in a circumnuclear star-forming ring in the luminous infrared galaxy ESO 154-G010, observed as part of the Public ESO Spectroscopic Survey of Transient Objects. SN 2013fc is both photometrically and spectroscopically similar to the well-studied type IIn SN 1998S and to the bright type II-L SN 1979C. It exhibits an initial linear decline, followed by a short plateau phase and a tail phase with a decline too fast for ⁵⁶Co decay with full γ -ray trapping. Initially, the spectrum was blue and featureless. Later on, a strong broad (~8000 km s^-1) H α emission profile became prominent. We apply a STARLIGHT stellar population model fit to the SN location (observed when the SN had faded) to estimate a high extinction of A_V = 2.9 ± 0.2 mag and an age of 10⁺³ _-2 Myr for the underlying cluster.We compare the SN to SNe 1998S and 1979C and discuss its possible progenitor star considering the similarities to these events. With a peak brightness of B = -20.46 ± 0.21 mag, SN 2013fc is 0.9 mag brighter than SN 1998S and of comparable brightness to SN 1979C.We suggest that SN 2013fc was consistent with a massive red supergiant (RSG) progenitor. Recent mass loss probably due to a strong RSG wind created the circumstellar matter illuminated through its interaction with the SN ejecta. We also observe a near-infrared excess, possibly due to newly condensed dust.

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).

The effects of transients on photospheric and chromospheric power distributions

We have observed a quiet Sun region with the Swedish 1-meter Solar Telescope (SST) equipped with CRISP Imaging SpectroPolarimeter. High-resolution, high-cadence, H-alpha line scanning images were taken to observe different layers of the solar atmosphere from the photosphere to upper chromosphere. We study the distribution of power in different period-bands at different heights. Power maps of the upper photosphere and the lower chromosphere show suppressed power surrounding the magnetic-network elements, known as "magnetic shadows". These also show enhanced power close to the photosphere, traditionally referred to as"power halos". The interaction between acoustic waves and inclined magnetic fields is generally believed to be responsible for these two effects. In this study we explore if small-scale transients can influence the distribution of power at different heights. We show that the presence of transients, like mottles, Rapid Blueshifted Excursions (RBEs) and Rapid Redshifted Excursions (RREs), can strongly influence the power-maps. The short and finite lifetime of these events strongly affects all powermaps, potentially influencing the observed power distribution. We show that Doppler-shifted transients like RBEs and RREs that occur ubiquitously, can have a dominant effect on the formation ofthe power halos in the quiet Sun. For magnetic shadows, transients like mottles do not seem to have a significant effect in the power suppression around 3 minutes and wave interaction may play a key role here. Our high cadence observations reveal that flows, waves and shocks manifest in presence of magnetic fields to form a non-linear magnetohydrodynamic system.

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