Refine
Year of publication
- 2015 (295) (remove)
Document Type
- Article (250)
- Doctoral Thesis (28)
- Postprint (7)
- Other (3)
- Conference Proceeding (2)
- Review (2)
- Habilitation Thesis (1)
- Master's Thesis (1)
- Preprint (1)
Keywords
- anomalous diffusion (8)
- stars: early-type (8)
- Magellanic Clouds (7)
- stars: massive (7)
- stars: atmospheres (6)
- stars: mass-loss (6)
- stars: winds, outflows (6)
- ISM: supernova remnants (5)
- X-rays: binaries (5)
- X-rays: stars (5)
Institute
- Institut für Physik und Astronomie (295) (remove)
We analyze star-type networks of phase oscillators by virtue of two methods. For identical oscillators we adopt the Watanabe-Strogatz approach, which gives full analytical description of states, rotating with constant frequency. For nonidentical oscillators, such states can be obtained by virtue of the self-consistent approach in a parametric form. In this case stability analysis cannot be performed, however with the help of direct numerical simulations we show which solutions are stable and which not. We consider this system as a model for a drum orchestra, where we assume that the drummers follow the signal of the leader without listening to each other and the coupling parameters are determined by a geometrical organization of the orchestra. (C) 2015 AIP Publishing LLC.
We analyze quasiperiodic partially synchronous states in an ensemble of Stuart-Landau oscillators with global nonlinear coupling. We reveal two types of such dynamics: in the first case the time-averaged frequencies of oscillators and of the mean field differ, while in the second case they are equal, but the motion of oscillators is additionally modulated. We describe transitions from the synchronous state to both types of quasiperiodic dynamics, and a transition between two different quasiperiodic states. We present an example of a bifurcation diagram, where we show the borderlines for all these transitions, as well as domain of bistability.
We study numerically secondary modes on top of a chaotic state in disordered nonlinear lattices. Two basic models are considered, with or without a local on-site potential. By performing periodic spatial modulation of displacement and kinetic energy, and following the temporal evolution of the corresponding spatial profiles, we reveal different modes which can be interpreted as first and second sound.
We study how coherence of noisy oscillations can be optimally enhanced by external locking. Based on the condition of minimizing the phase diffusion constant, we find the optimal forcing explicitly in the limits of small and large noise, in dependence of the phase sensitivity of the oscillator. We show analytically that the form of the optimal force bifurcates with the noise intensity; this is confirmed by the analysis of an optimal locking forcing for an experimentally obtained phase sensitivity of a neural cell. In the limit of small noise, the results are compared with purely deterministic conditions of optimal locking.
Partial synchronization in networks of non-linearly coupled oscillators: The Deserter Hubs Model
(2015)
We study the Deserter Hubs Model: a Kuramoto-like model of coupled identical phase oscillators on a network, where attractive and repulsive couplings are balanced dynamically due to nonlinearity of interactions. Under weak force, an oscillator tends to follow the phase of its neighbors, but if an oscillator is compelled to follow its peers by a sufficient large number of cohesive neighbors, then it actually starts to act in the opposite manner, i.e., in anti-phase with the majority. Analytic results yield that if the repulsion parameter is small enough in comparison with the degree of the maximum hub, then the full synchronization state is locally stable. Numerical experiments are performed to explore the model beyond this threshold, where the overall cohesion is lost. We report in detail partially synchronous dynamical regimes, like stationary phase-locking, multistability, periodic and chaotic states. Via statistical analysis of different network organizations like tree, scale-free, and random ones, we found a measure allowing one to predict relative abundance of partially synchronous stationary states in comparison to time-dependent ones. (C) 2015 AIP Publishing LLC.
We describe synchronization transitions in an ensemble of globally coupled phase oscillators with a bi-harmonic coupling function, and two sources of disorder-diversity of the intrinsic oscillators' frequencies, and external independent noise forces. Based on the self-consistent formulation, we derive analytic solutions for different synchronous states. We report on various non-trivial transitions from incoherence to synchrony, with the following possible scenarios: simple supercritical transition (similar to classical Kuramoto model); subcritical transition with large area of bistability of incoherent and synchronous solutions; appearance of a symmetric two-cluster solution which can coexist with the regular synchronous state. We show that the interplay between relatively small white noise and finite-size fluctuations can lead to metastability of the asynchronous solution.
We generalize the Kuramoto model of globally coupled oscillators to multifrequency communities. A situation when mean frequencies of two subpopulations are close to the resonance 2 : 1 is considered in detail. We construct uniformly rotating solutions describing synchronization inside communities and between them. Remarkably, cross coupling across the frequencies can promote synchrony even when ensembles are separately asynchronous. We also show that the transition to synchrony due to the cross coupling is accompanied by a huge multiplicity of distinct synchronous solutions, which is directly related to a multibranch entrainment. On the other hand, for synchronous populations, the cross-frequency coupling can destroy phase locking and lead to chaos of mean fields.
Paper-based microfluidics provide an inexpensive, easy to use technology for point-of-care diagnostics in developing countries. Here, we combine paper-based microfluidic devices with responsive hydrogels to add an entire new class of functions to these versatile low-cost fluidic systems. The hydrogels serve as fluid reservoirs. In response to an external stimulus, e.g. an increase in temperature, the hydrogels collapse and release fluid into the structured paper substrate. In this way, chemicals that are either stored on the paper substrate or inside the hydrogel pads can be dissolved, premixed, and brought to reaction to fulfill specific analytic tasks. We demonstrate that multi-step sequences of chemical reactions can be implemented in a paper-based system and operated without the need for external precision pumps. We exemplify this technology by integrating an antibody-based E. coli test on a small and easy to use paper device.
X-ray spectroscopy is a powerful tool to study the local charge distribution of chemical systems. Together with the liquid jet it becomes possible to probe chemical systems in their natural environment, the liquid phase. In this work, we present X-ray absorption (XA), X-ray emission (XE) and resonant inelastic X-ray scattering (RIXS) data of pure water and various salt solutions and show the possibilities these methods offer to elucidate the nature of ion-water interaction.
Liquid water molecules interact strongly with each other, forming a fluctuating hydrogen bond network and thereby giving rise to the anomalous phase diagram of liquid water. Consequently, symmetric and asymmetric water molecules have been found in the picosecond time average with IR and optical Raman spectroscopy. With subnatural linewidth resonant inelastic x-ray scattering (RIXS) at vibrational resolution, we take sub-femtosecond snapshots of the electronic and structural properties of water molecules in the hydrogen bond network. We derive a strong dominance of nonsymmetric molecules in liquid water in contrast to the gas phase on the sub-femtosecond timescale of RIXS and determine the fraction of highly asymmetrically distorted molecules.
Rydberg-Resolved Resonant Inelastic Soft X-Ray Scattering: Dynamics at Core Ionization Thresholds
(2015)
Resonant inelastic x-ray scattering spectra excited in the immediate vicinity of the core-level ionization thresholds of N-2 have been recorded. Final states of well-resolved symmetry-selected Rydberg series converging to valence-level ionization thresholds with vibrational excitations are observed. The results are well described by a quasi-two-step model which assumes that the excited electron is unaffected by the radiative decay. This threshold dynamics simplifies the interpretation of resonant inelastic x-ray scattering spectra considerably and facilitates characterization of low-energy excited final states in molecular systems.
Wolf-Rayet Stars
(2015)
Nearly 150 years ago, the French astronomers Charles Wolf and Georges Rayet described stars with very conspicuous spectra that are dominated by bright and broad emission lines. Meanwhile termed Wolf-Rayet Stars after their discoverers, those objects turned out to represent important stages in the life of massive stars.
As the first conference in a long time that was specifically dedicated to Wolf-Rayet stars, an international workshop was held in Potsdam, Germany, from 1.-5. June 2015. About 100 participants, comprising most of the leading experts in the field as well as as many young scientists, gathered for one week of extensive scientific exchange and discussions. Considerable progress has been reported throughout, e.g. on finding such stars, modeling and analyzing their spectra, understanding their evolutionary context, and studying their circumstellar nebulae. While some major questions regarding Wolf-Rayet stars still remain open 150 years after their discovery, it is clear today that these objects are not just interesting stars as such, but also keystones in the evolution of galaxies.
These proceedings summarize the talks and posters presented at the Potsdam Wolf-Rayet workshop. Moreover, they also include the questions, comments, and discussions emerging after each talk, thereby giving a rare overview not only about the research, but also about the current debates and unknowns in the field. The Scientific Organizing Committee (SOC) included Alceste Bonanos (Athens), Paul Crowther (Sheffield), John Eldridge (Auckland), Wolf-Rainer Hamann (Potsdam, Chair), John Hillier (Pittsburgh), Claus Leitherer (Baltimore), Philip Massey (Flagstaff), George Meynet (Geneva), Tony Moffat (Montreal), Nicole St-Louis (Montreal), and Dany Vanbeveren (Brussels).
A significant number of the central stars of planetary nebulae (CSPNe) are hydrogen-deficient, showing a chemical composition of helium, carbon, and oxygen. Most of them exhibit Wolf-Rayet-like emission line spectra, similar to those of the massive WC Pop I stars, and are therefore classified as of spectral type [WC]. In the last years, CSPNe of other Wolf-Rayet spectral subtypes have been identified, namely PB 8, which is of spectral type [WN/C], and IC 4663 and Abell 48, which are of spectral type [WN]. We review spectral analyses of Wolf-Rayet type central stars of different evolutionary stages and discuss the results in the context of stellar evolution. Especially we consider the question of a common evolutionary channel for [WC] stars. The constraints on the formation of [WN] or [WC/N] subtype stars will also be addressed.
The distribution of angular momentum in massive stars is a critical component of their evolution, yet not much is known on the rotation velocities of Wolf-Rayet stars. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars by implementing a 3D integration scheme in the PoWR code. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may imply rapid rotation, infer the corresponding rotation parameters, and discuss the implications of our results. We find that rotation helps to reproduce the unique spectra analyzed here. However, if rotation is indeed involved, the inferred rotational velocities at the stellar surface are large (∼ 200 km/s), and the implied co-rotation radii (∼ 10R∗) suggest the existence of very strong photospheric magnetic fields (∼ 20 kG).
In the last decades, stellar atmosphere codes have become a key tool in understanding massive stars, including precise calculations of stellar and wind parameters, such as temperature, massloss rate, and terminal wind velocity. Nevertheless, for these models the hydrodynamic equation is not solved in the wind. Motivated by the results of the CAK theory, the models typically use a beta velocity law, which however turns out not to be adequate for stars with very strong winds, and treat the mass-loss rate as a free parameter. In a new branch of the Potsdam Wolf-Rayet model atmosphere (PoWR) code, we solve the hydrodynamic equation consistently throughout the stellar atmosphere. The PoWR code performs the calculation of the radiative force without approximations (e.g. Sobolev). We show the impact of hydrodynamically consistent modelling on OB and WR stars in comparison to conventional models and discuss the obtained velocity fields and their impact on the observed spectral lines.
Macroclumping in WR 136
(2015)
Macroclumping proved to resolve the discordance between different mass-loss rate diagnostics for O-type stars, in particular between Hα and the P v resonance lines. In this paper, we report first results from a corresponding investigation for WR stars. We apply our detailed 3-D Monte Carlo (MC) line formation code to the P v resonance doublet and show, for the Galactic WNL star WR136, that macroclumping is require to bring this line in accordance with the mass-loss rate derived from the emission-line spectrum.
Context. Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars.
Aims. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars.
Methods. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations.
Results. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10(5.5) to 10(6.1) L-circle dot. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC.
Conclusions. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.
Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres.
Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra.
Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses.
Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.
We present new grids of Potsdam Wolf-Rayet (PoWR) model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN stars). The models have been calculated with the latest version of the PoWR stellar atmosphere code for spherical stellar winds. The WN model atmospheres include the non-LTE solutions of the statistical equations for complex model atoms, as well as the radiative transfer equation in the co-moving frame. Iron-line blanketing is treated with the help of the superlevel approach, while wind inhomogeneities are taken into account via optically thin clumps. Three of our model grids are appropriate for Galactic metallicity. The hydrogen mass fraction of these grids is 50%, 20%, and 0%, thus also covering the hydrogen-rich late-type WR stars that have been discovered in recent years. Three grids are adequate for LMC WN stars and have hydrogen fractions of 40%, 20%, and 0%. Recently, additional grids with SMC metallicity and with 60%, 40%, 20%, and 0% hydrogen have been added. We provide contour plots of the equivalent widths of spectral lines that are usually used for classification and diagnostics.
As WR 6 is a putatively single WN4 star, and is relatively bright (V = 6.9), it is an ideal case for studying the wind mechanisms in these extremely luminous stars. To obtain higher resolution spectra at higher energy (above 1 keV) than previously obtained with the XMM/Newton RGS, we have observed WR 6 with the Chandra High Energy Transmission Grating Spectrometer for 450 ks. We have resolved emission lines of S, Si, Mg, Ne, and Fe, which all show a “fin"-shaped prole, characteristic of a self-absorbed uniformly expanding shell. Steep blue edges gives robust maximal expansion velocities of about 2000 km/s, somewhat larger than the 1700km/s derived from UV lines. The He-like lines all indicate that X-ray emitting plasmas are far from the photosphere – even at the higher energies where opacity is lowest { as was also the case for the longer wavelength lines observed with XMM-Newton/RGS. Abundances determined from X-ray spectral modeling indicate enhancements consistent with nucleosynthesis. The star was also variable in X-rays and in simultaneous optical photometry obtained with Chandra aspect camera, but not coherently with the optically known period of 3.765 days.
We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system delta Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P > 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.
We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeK alpha emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeK alpha in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeK alpha spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in gamma Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeK alpha is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeK alpha (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeK alpha is narrow (sigma(line) < 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V similar to 1000 km s(-1)). (6) The equivalent hydrogen column (N-H) directly correlates to the EW of FeK alpha, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object - wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor's wind in the FeK alpha emission and the absorption when the donor is a supergiant massive star.
While the majority of very massive stars is clearly found in clusters, there are also very massive objects not associated with any cluster, suggesting they may have been born in isolation. In order to gain more insights, we studied the regions around two WR stars in the Galactic Center region. To understand the nature of the potential cluster around massive stars, photometry alone is not sufficient. We therefore used the ESO VLT/SINFONI integral field spectrograph to obtain photometry and spectra for the whole region around our two candidate stars. In total, more than 60 stars have been found and assigned a spectral type.
We present an analysis of a 78 ks Chandra high-energy transmission gratings observation of the B0I star QV Nor, the massive donor of the wind-accreting pulsar 4U1538-52. The neutron star (NS) orbits its companion in a very close orbit (r < 1.4R(*), in units of the stellar radii), thereby allowing probing of the innermost wind regions. The flux of the Fe K alpha line during eclipse reduces to only similar to 30% of the flux measured out of eclipse. This indicates that the majority of Fe fluorescence must be produced in regions close to the NS, at distances smaller than 1R(*) from its surface. The fact that the flux of the continuum decreases to only similar to 3% during eclipse allows for a high contrast of the Fe Ka line fluorescence during eclipse. The line is not resolved and centered at lambda = 1.9368(-0.0018)(+0.0032) angstrom. From the inferred plasma speed limit of v < c Delta lambda/lambda < 800 km s(-1) and range of ionization parameters of log xi =[-1, 2], together with the stellar density profile, we constrain the location of the cold, dense material in the stellar wind of QV Nor using simple geometrical considerations. We then use the Fe K alpha line fluorescence as a tracer of wind clumps and determine that these clumps in the stellar wind of QV Nor (B0I) must already be present at radii r < 1.25R(*), close to the photosphere of the star.
The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS 2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33 000 +/- 1000K and 4.0 +/- 0.1 dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H alpha shows a variable emission. The H alpha emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.
Only a small fraction of massive stars seem to host a measurable structured magnetic field, whose origin is still unknown and whose implications for stellar evolution still need to be assessed. Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD 44743; B1 II/III) and epsilon CMa (HD 52089; B1.5II) in December 2013 and April 2014. For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field, approximately constant with time. We determined the physical parameters of both stars and characterise their X-ray spectrum. For the beta Cep star beta CMa, our mode identification analysis led to determining a rotation period of 13.6 +/- 1.2 days and of an inclination angle of the rotation axis of 57.6 +/- 1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of about 22 degrees and a dipolar magnetic field strength (B-d) of about 100 G (60 < B-d < 230 G within the 1 sigma level), below what is typically found for other magnetic massive stars. This conclusion is strengthened further by considerations of the star's X-ray spectrum. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Our results imply that both stars are expected to have a dynamical magnetosphere, so the magnetic field is not able to support a circumstellar disk. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.
Aims. Massive B-type stars with strong magnetic fields and fast rotation are very rare and pose a mystery for theories of star formation and magnetic field evolution. Only two such stars, called sigma Ori E analogues, were known until recently. A team involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating magnetosphere stars, HD 23478 and HD 345439. The magnetic fields in these newly discovered sOri E analogues have not been investigated so far.
Methods. In the framework of our ESO Large Programme and one normal ESO programme, we carried out low-resolution FORS 2 spectropolarimetric observations of HD 23478 and HD 345439.
Results. In the measurements of hydrogen lines, we discover a rather strong longitudinal magnetic field of up to 1.5 kG in HD 23478 and up to 1.3 kG using the entire spectrum. The analysis of HD 345439 using four subsequent spectropolarimetric subexposures does not reveal a magnetic field at a significance level of 3 sigma. On the other hand, individual subexposures indicate that HD 345439 may host a strong magnetic field that rapidly varies over 88 min. The fast rotation of HD 345439 is also indicated by the behaviour of several metallic and He I lines in the low-resolution FORS 2 spectra that show profile variations already on this short time-scale.
Context. beta Cep-type variables are early B-type stars that are characterized by oscillations observable in their optical light curves. At least one beta Cep-variable also shows periodic variability in X-rays.
Aims. Here we study the X-ray light curves in a sample of beta Cep-variables to investigate how common X-ray pulsations are for this type of stars.
Methods. We searched the Chandra and XMM-Newton X-ray archives and selected stars that were observed by these telescopes for at least three optical pulsational periods. We retrieved and analyzed the X-ray data for kappa Sco, beta Cru, and alpha Vir. The X-ray light curves of these objects were studied to test for their variability and periodicity.
Results. While there is a weak indication for X-ray variability in beta Cru, we find no statistically significant evidence of X-ray pulsations in any of our sample stars. This might be due either to the insufficient data quality or to the physical lack of modulations. New, more sensitive observations should settle this question.
Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system delta Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the Hipparcos parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if delta Ori lies at about twice the Hipparcos distance, in the vicinity of the sigma-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be Delta V approximate to 2.(m)8. The inferred parameters suggest that the secondary is an early B-type dwarf (approximate to B1 V), while the tertiary is an early B-type subgiant (approximate to B0 IV). We find evidence for rapid turbulent velocities (similar to 200 km s(-1)) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind (logL(X)/L-Bol approximate to -6.85), initiating close to the stellar surface at R-0 similar to 1.1 R-*. Accounting for clumping, the mass-loss rate of the primary is found to be log (M) over dot approximate to -6.4 (M-circle dot yr(-1))., which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical, and radio spectral domains.
With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear-processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars.
Projections of changes in Antarctic Ice Sheet (AIS) surface mass balance indicate a negative contribution to sea level because of the expected increase in precipitation due to the higher moisture holding capacity of warmer air(1). Observations over the past decades, however, are unable to constrain the relation between temperature and accumulation changes because both are dominated by strong natural variability(2-5). Here we derive a consistent continental-scale increase in accumulation of approximately 5 +/- 1% K-1, through the assessment of ice-core data (spanning the large temperature change during the last deglaciation, 21,000 to 10,000 years ago), in combination with palaeo-simulations, future projections by 35 general circulation models (GCMs), and one high-resolution future simulation. The ice-core data and modelling results for the last deglaciation agree, showing uniform local sensitivities of similar to 6% K-1. The palaeo-simulation allows for a continental-scale aggregation of accumulation changes reaching 4.3% K-1. Despite the different timescales, these sensitivities agree with the multi-model mean of 6.1 +/- 2.6% K-1 (GCMprojections) and the continental-scale sensitivity of 4.9% K-1 (high-resolution future simulation). Because some of the mass gain of the AIS is offset by dynamical losses induced by accumulation(6,7), we provide a response function allowing projections of sea-level fall in terms of continental-scale accumulation changes that compete with surface melting and dynamical losses induced by other mechanisms(6,8,9).
We report on rendering polyelectrolyte brushes photosensitive by loading them with azobenzene-containing cationic surfactants. Planar poly(methacrylic acid) (PMAA) brushes are synthesized using the “grafting from” free-radical polymerization scheme followed by exposure to a solution of photosensitive surfactants consisting of positively-charged head groups and hydrophobic tails into which azobenzene moieties are inserted. In this study the length of the hydrophobic methylene spacer connecting the azobenzene and the charged head group ranges from 4 to 10 CH2 groups. Under irradiation with UV light, the photo-isomerization of azobenzene integrated into a surfactant results in a change in size, geometry, dipole moment and free volume of the whole molecule. When the brush loaded with photosensitive surfactants is exposed to irradiation with UV interference patterns, the topography of the brush deforms following the distribution of the light intensity, exhibiting surface relief gratings (SRG). Since SRG formation is accompanied by a local rupturing of polymer chains in areas from which the polymer material is receding, most of the polymer material is removed from the surface during treatment with good solvent, leaving behind characteristic patterns of lines or dots. The azobenzene molecules still integrated within the polymer film can be removed by washing the brush with water. The remaining nano-structured brush can then be re-used for further functionalization. Although the opto-mechanically induced rupturing occurs for all surfactants, larger species do not penetrate deep enough into the brush such that after rupturing a leftover layer of polymer material remains on the substrate. This indicates that rupturing occurs predominantly in regions of high surfactant density.
Here, we report on two photosensitive amorphous polymers showing opposite behavior upon exposure to illumination. The first polymer (PAZO) consists of linear backbone to which azobenzene-containing side chains are covalently attached, while in the second polymer (azo-PEI), the azobenzene side chains are attached ionically to a polyelectrolyte backbone. When irradiated through a mask, the PAZO goes away from the intensity maxima, leaving behind topography trenches, while the direction of the mass transport of the azo-PEI polymer points towards the intensity maxima. This kind of behavior has been reported only for certain liquid crystalline polymers that exhibit in-phase reaction on illumination, that is, topography maxima coincides with the intensity maxima. Furthermore, flat nanocrystals placed on top of azo-PEI film was found to be moved together with the mass transport of the underlying polymer film as visualized using in situ atomic force microscopy (AFM) measurements. It was also demonstrated that the two polymer films respond differently on irradiation with the polarization and intensity interference patterns (IPs). To record the kinetic of the surface relief grating formation within two polymers during irradiation with different IPs, we utilized a homemade setup combining the optical part for the generation of IP and AFM. A possible mechanism explaining different responses on the irradiation of amorphous polymers is discussed in the frame of a theoretical model proposed by Saphiannikova et al. (J. Phys. Chem. B 113, 5032-5045 (2009)).
We report on the experimental and theoretical investigation of a considerable increase in the rate for thermal cis -> trans isomerization of azobenzene-containing molecules in the presence of gold nanopartides. Experimentally, by means of UV vis spectroscopy, we studied a series of azobenzene-containing surfactants and 4-nitroazobenzene. We found that in the presence of gold,nanoparticles the thermal lifetime of the cis isomer of the azobenzenecontaining molecules was decreased by up to 3 orders of magnitude in comparison to the lifetime in solution without nanoparticles. The electron transfer between azobenzene-containing molecules and a surface of gold nanopartides is a possible reason to promote the thermal cis trans switching. To investigate the effect of electron attachment to, and withdrawal from, the azobenzene-containing molecules on the isomerization rate, we performed density functional theory calculations of activation energy barriers of the reaction together with Eyring's transition state theory calculations of the rates for azobenzene derivatives with donor and acceptor groups in para position of one of the phenyl rings, as well as for one of the azobenzene-containing surfactants. We found that activation barriers are greatly lowered for azobenzene-containing molecules, both upon electron attachment and withdrawal, which leads, in turn, to a dramatic increase in the thermal isomerization rate.
DNA binding to trans- and cis-isomers of azobenzene containing cationic surfactant in 5 mM NaCl solution was investigated by the methods of dynamic light scattering (DLS), low-gradient viscometry (LGV), atomic force microscopy (AFM), circular dichroism (CD), gel electrophoresis (GE), flow birefringence (FB), UV-Vis spectrophotometry. Light-responsive conformational transitions of DNA in complex with photosensitive surfactant, changes in DNA optical anisotropy and persistent length, phase transition of DNA into nanoparticles induced by high surfactant concentration, as well as transformation of surfactant conformation under its binding to macromolecule were studied. Computer simulations of micelles formation for cis- and trans-isomers of azobenzene containing surfactant, as well as DNA-surfactant interaction, were carried out. Phase diagram for DNA-surfactant solutions was designed. The possibility to reverse the DNA packaging induced by surfactant binding with the dilution and light irradiation was shown. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 109-122, 2015.
Many chemical reactions in biological cells occur at very low concentrations of constituent molecules. Thus, transcriptional gene-regulation is often controlled by poorly expressed transcription-factors, such as E.coli lac repressor with few tens of copies. Here we study the effects of inherent concentration fluctuations of substrate-molecules on the seminal Michaelis-Menten scheme of biochemical reactions. We present a universal correction to the Michaelis-Menten equation for the reaction-rates. The relevance and validity of this correction for enzymatic reactions and intracellular gene-regulation is demonstrated. Our analytical theory and simulation results confirm that the proposed variance-corrected Michaelis-Menten equation predicts the rate of reactions with remarkable accuracy even in the presence of large non-equilibrium concentration fluctuations. The major advantage of our approach is that it involves only the mean and variance of the substrate-molecule concentration. Our theory is therefore accessible to experiments and not specific to the exact source of the concentration fluctuations.
We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral data set of the delta Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of approximate to 479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range of 5-25 is is confirmed, with a maximum amplitude of about +/- 15% within a single approximate to 125 ks observation. Periods of 4.76 and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in the flux level throughout the nine-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate the variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi = 0.0 when the secondary delta Ori Aa2 is at the inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.
A rather strong mean longitudinal magnetic field of the order of a few hundred gauss was detected a few years ago in the Of?p star CPD -28 degrees 2561 using FORS2 (FOcal Reducer low dispersion Spectrograph 2) low-resolution spectropolarimetric observations. In this work, we present additional low-resolution spectropolarimetric observations obtained during several weeks in 2013 December using FORS 2 mounted at the 8-m Antu telescope of the Very Large Telescope (VLT). These observations cover a little less than half of the stellar rotation period of 73.41 d mentioned in the literature. The behaviour of the mean longitudinal magnetic field is consistent with the assumption of a single-wave variation during the stellar rotation cycle, indicating a dominant dipolar contribution to the magnetic field topology. The estimated polar strength of the surface dipole B-d is larger than 1.15 kG. Further, we compared the behaviour of the line profiles of various elements at different rotation phases associated with different magnetic field strengths. The strongest contribution of the emission component is observed at the phases when the magnetic field shows a negative or positive extremum. The comparison of the spectral behaviour of CPD -28 degrees 2561 with that of another Of?p star, HD 148937 of similar spectral type, reveals remarkable differences in the degree of variability between both stars. Finally, we present new X-ray observations obtained with the Suzaku X-ray Observatory. We report that the star is X-ray bright with log L-X/L-bol approximate to -5.7. The low-resolution X-ray spectra reveal the presence of a plasma heated up to 24 MK. We associate the 24 MK plasma in CPD -28 degrees 2561 with the presence of a kG strong magnetic field capable to confine stellar wind.
In this review I briefly summarize our knowledge of the X-ray emission from single WN, WC, and WO stars. These stars have relatively modest X-ray luminosities, typically not exceeding 1L⊙. The analysis of X-ray spectra usually reveals thermal plasma with temperatures reaching a few x10 MK. X-ray variability is detected in some WN stars. At present we don't fully understand how X-ray radiation in produced in WR stars, albeit there are some promising research avenues, such as the presence of CIRs in the winds of some stars. To fully understand WR stars we need to unravel mechanisms of X-ray production in their winds.
We present the results of a broadband X-ray study of the enigmatic Be star Gamma Cassiopeiae (herein gamma Cas) based on observations made with both the Suzaku and INTEGRAL observatories.. Cas has long been recognized as the prototypical example of a small subclass of Be stars with moderately strong X-ray emission dominated by a hot thermal component in the 0.5-12 keV energy range (L-x approximate to 10(32)-10(33) erg s(-1)). This places them at the high end of the known luminosity distribution for stellar emission, but several orders of magnitude below typical accretion-powered Be X-ray binaries. The INTEGRAL observations spanned an eight-year baseline and represent the deepest measurement to date at energies above similar to 50 keV. We find that the INTEGRAL data are consistent within statistics to a constant intensity source above 20 keV, with emission extending up to similar to 100 keV, and that searches for all of the previously reported periodicities of the system at lower energies led to null results. We further find that our combined Suzaku and INTEGRAL spectrum, which we suggest is the most accurate broadband X-ray measurement of gamma Cas to date, is fitted extremely well with a thermal plasma emission model with a single absorption component. We found no compelling need for an additional non-thermal high-energy component. We discuss these results in the context of a currently favored models for gamma Cas and its analogs.
The structure of bulk liquid water was recently probed by x-ray scattering below the temperature limit of homogeneous nucleation (T-H) of similar to 232 K [J. A. Sellberg et al., Nature 510, 381-384 (2014)]. Here, we utilize a similar approach to study the structure of bulk liquid water below T-H using oxygen K-edge x-ray emission spectroscopy (XES). Based on previous XES experiments [T. Tokushima et al., Chem. Phys. Lett. 460, 387-400 (2008)] at higher temperatures, we expected the ratio of the 1b(1)' and 1b(1)" peaks associated with the lone-pair orbital in water to change strongly upon deep supercooling as the coordination of the hydrogen (H-) bonds becomes tetrahedral. In contrast, we observed only minor changes in the lone-pair spectral region, challenging an interpretation in terms of two interconverting species. A number of alternative hypotheses to explain the results are put forward and discussed. Although the spectra can be explained by various contributions from these hypotheses, we here emphasize the interpretation that the line shape of each component changes dramatically when approaching lower temperatures, where, in particular, the peak assigned to the proposed disordered component would become more symmetrical as vibrational interference becomes more important. (C) 2015 AIP Publishing LLC.
The discovery of ultrafast X-ray induced optical reflectivity changes enabled the development of X-ray/optical cross correlation techniques at X-ray free electron lasers worldwide. We have now linked through experiment and theory the fundamental excitation and relaxation steps with the transient optical properties in finite solid samples. Therefore, we gain a thorough interpretation and an optimized detection scheme of X-ray induced changes to the refractive index and the X-ray/optical cross correlation response. (C) 2015 AIP Publishing LLC.
The initiation of a marine ice-sheet instability (MISI) is generally discussed from the ocean side of the ice sheet. It has been shown that the reduction in ice-shelf buttressing and softening of the coastal ice can destabilize a marine ice sheet if the bedrock is sloping upward towards the ocean. Using a conceptional flow-line geometry, we investigate the possibility of whether a MISI can be triggered from the direction of the ice divide as opposed to coastal forcing and explore the interaction between connected basins. We find that the initiation of a MISI in one basin can induce a destabilization in the other. The underlying mechanism of basin interaction is based on dynamic thinning and a consecutive motion of the ice divide which induces a thinning in the adjacent basin and a successive initiation of the instability. Our simplified and symmetric topographic setup allows scaling both the geometry and the transition time between both instabilities. We find that the ice profile follows a universal shape that is scaled with the horizontal extent of the ice sheet and that the same exponent of 1/2 applies for the scaling relation between central surface elevation and horizontal extent as in the pure shallow ice approximation (Vialov profile). Altering the central bed elevation, we find that the extent of grounding-line retreat in one basin determines the degree of interaction with the other. Different scenarios of basin interaction are discussed based on our modeling results as well as on a conceptual flux-balance analysis. We conclude that for the three-dimensional case, the possibility of drainage basin interaction on timescales on the order of 1 kyr or larger cannot be excluded and hence needs further investigation.
Anthropogenic carbon emissions lock in long-term sea-level rise that greatly exceeds projections for this century, posing profound challenges for coastal development and cultural legacies. Analysis based on previously published relationships linking emissions to warming and warming to rise indicates that unabated carbon emissions up to the year 2100 would commit an eventual global sea-level rise of 4.3-9.9 m. Based on detailed topographic and population data, local high tide lines, and regional long-term sea-level commitment for different carbon emissions and ice sheet stability scenarios, we compute the current population living on endangered land at municipal, state, and national levels within the United States. For unabated climate change, we find that land that is home to more than 20 million people is implicated and is widely distributed among different states and coasts. The total area includes 1,185-1,825 municipalities where land that is home to more than half of the current population would be affected, among them at least 21 cities exceeding 100,000 residents. Under aggressive carbon cuts, more than half of these municipalities would avoid this commitment if the West Antarctic Ice Sheet remains stable. Similarly, more than half of the US population-weighted area under threat could be spared. We provide lists of implicated cities and state populations for different emissions scenarios and with and without a certain collapse of the West Antarctic Ice Sheet. Although past anthropogenic emissions already have caused sea-level commitment that will force coastal cities to adapt, future emissions will determine which areas we can continue to occupy or may have to abandon.
The future evolution of the Antarctic Ice Sheet represents the largest uncertainty in sea-level projections of this and upcoming centuries. Recently, satellite observations and high-resolution simulations have suggested the initiation of an ice-sheet instability in the Amundsen Sea sector of West Antarctica, caused by the last decades' enhanced basal ice-shelf melting. Whether this localized destabilization will yield a full discharge of marine ice from West Antarctica, associated with a global sea-level rise of more than 3 m, or whether the ice loss is limited by ice dynamics and topographic features, is unclear. Here we show that in the Parallel Ice Sheet Model, a local destabilization causes a complete disintegration of the marine ice in West Antarctica. In our simulations, at 5-km horizontal resolution, the region disequilibrates after 60 y of currently observed melt rates. Thereafter, the marine ice-sheet instability fully unfolds and is not halted by topographic features. In fact, the ice loss in Amundsen Sea sector shifts the catchment's ice divide toward the Filchner-Ronne and Ross ice shelves, which initiates grounding-line retreat there. Our simulations suggest that if a destabilization of Amundsen Sea sector has indeed been initiated, Antarctica will irrevocably contribute at least 3 m to global sea-level rise during the coming centuries to millennia.
Regional and sectoral disaggregation of multi-regional input-output tables - a flexible algorithm
(2015)
A common shortcoming of available multi-regional input-output (MRIO) data sets is their lack of regional and sectoral detail required for many research questions (e.g. in the field of disaster impact analysis). We present a simple algorithm to refine MRIO tables regionally and/or sectorally. By the use of proxy data, each MRIO flow in question is disaggregated into the corresponding sub-flows. This downscaling procedure is complemented by an adjustment rule ensuring that the sub-flows match the superordinate flow in sum. The approximation improves along several iteration steps. The algorithm unfolds its strength through the flexible combination of multiple, possibly incomplete proxy data sources. It is also flexible in a sense that any target sector and region resolution can be chosen. As an exemplary case we apply the algorithm to a regional and sectoral refinement of the Eora MRIO database.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
In this paper, we discuss recent progress in research of ensembles of mean field coupled oscillators. Without an ambition to present a comprehensive review, we outline most interesting from our viewpoint results and surprises, as well as interrelations between different approaches. (c) 2015 AIP Publishing LLC.
We examine the non-ergodic properties of scaled Brownian motion (SBM), a non-stationary stochastic process with a time dependent diffusivity of the form D(t) similar or equal to t(alpha-1). We compute the ergodicity breaking parameter EB in the entire range of scaling exponents a, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths T and short lag times Delta the EB parameter as function of the scaling exponent a has no divergence at alpha - 1/2 and present the asymptotes for EB in different limits. We generalize the analytical and simulations results for the time averaged and ergodic properties of SBM in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion.
Aging scaled Brownian motion
(2015)
Scaled Brownian motion (SBM) is widely used to model anomalous diffusion of passive tracers in complex and biological systems. It is a highly nonstationary process governed by the Langevin equation for Brownian motion, however, with a power-law time dependence of the noise strength. Here we study the aging properties of SBM for both unconfined and confined motion. Specifically, we derive the ensemble and time averaged mean squared displacements and analyze their behavior in the regimes of weak, intermediate, and strong aging. A very rich behavior is revealed for confined aging SBM depending on different aging times and whether the process is sub- or superdiffusive. We demonstrate that the information on the aging factorizes with respect to the lag time and exhibits a functional form that is identical to the aging behavior of scale-free continuous time random walk processes. While SBM exhibits a disparity between ensemble and time averaged observables and is thus weakly nonergodic, strong aging is shown to effect a convergence of the ensemble and time averaged mean squared displacement. Finally, we derive the density of first passage times in the semi-infinite domain that features a crossover defined by the aging time.
We present 3D numerical simulations of the NGC6888 nebula considering the proper motion and the evolution of the star, from the red supergiant (RSG) to the Wolf-Rayet (WR) phase. Our simulations reproduce the limb-brightened morphology observed in [OIII] and X-ray emission maps. The synthetic maps computed by the numerical simulations show filamentary and clumpy structures produced by instabilities triggered in the interaction between the WR wind and the RSG shell.
The origin of cosmic rays was the subject of several studies for over a century. The investigations done within this dissertation are one small step to shed some more light on this mystery.
Locating the sources of cosmic rays is not trivial due to the interstellar magnetic field. However, the Hillas criterion allows us to arrive at the conclusion that supernova remnants are our main suspect for the origin of galactic cosmic rays. The mechanism by which they are accelerating particles is found within the field of shock physics as diffusive shock acceleration. To allow particles to enter this process also known as Fermi acceleration pre-acceleration processes like shock surfing acceleration and shock drift acceleration are necessary. Investigating the processes happening in the plasma shocks of supernova remnants is possible by utilising a simplified model which can be simulated on a computer using Particle-in-Cell simulations.
We developed a new and clean setup to simulate the formation of a double shock, i.e., consisting of a forward and a reverse shock and a contact discontinuity, by the collision of two counter-streaming plasmas, in which a magnetic field can be woven into. In a previous work, we investigated the processes at unmagnetised and at magnetised parallel shocks, whereas in the current work, we move our investigation on to magnetised perpendicular shocks.
Due to a much stronger confinement of the particles to the collision region the perpendicular shock develops much faster than the parallel shock. On the other hand, this leads to much weaker turbulence. We are able to find indications for shock surfing acceleration and shock drift acceleration happening at the two shocks leading to populations of pre-accelerated particles that are suitable as a seed population to be injected into further diffusive shock acceleration to be accelerated to even higher energies. We observe the development of filamentary structures in the shock ramp of the forward shock, but not at the reverse shock. This leads to the conclusion that the development of such structures in the shock ramp of quasi-perpendicular collisionless shocks might not necessarily be determined by the existence of a critical sonic Mach number but by a critical shock speed.
The results of the investigations done within this dissertation might be useful for further studies of oblique shocks and for studies using hybrid or magnetohydrodynamic simulations. Together with more sophisticated observational methods, these studies will help to bring us closer to an answer as to how particles can be accelerated in supernova remnants and eventually become cosmic rays that can be detected on Earth.
G349.7+0.2 is a young Galactic supernova remnant (SNR) located at the distance of 11.5 kpc and observed across the entire electromagnetic spectrum from radio to high energy (HE; 0.1 GeV < E < 100 GeV) gamma-rays. Radio and infrared observations indicate that the remnant is interacting with a molecular cloud. In this paper, the detection of very high energy (VHE, E > 100 GeV) gamma-ray emission coincident with this SNR with the High Energy Stereoscopic System (HESS.) is reported. This makes it one of the farthest Galactic SNR ever detected in this domain. An integral flux F(E > 400 GeV) = (6.5 +/- 1.1(stat) +/- 1.3(syst)) x 10-11 ph cm(-2) s(-1) corresponding to similar to 0.7% of that of the Crab Nebula and to a luminosity of similar to 10(34) erg s(-1) above the same energy threshold, and a steep photon index Gamma(VHE) = 2.8 +/- 0.27(stat) +/- 0.20(syst) are measured. The analysis of more than 5 yr of Fermi-LAT data towards this source shows a power-law like spectrum with a best-fit photon index Gamma(HE) = 2.2 +/- 0.04.2(stat-0.31sys)(+0.13), The combined gamma-ray spectrum of 0349.7+0.2 can be described by either a broken power law (I3PL) or a power law with exponential (or sub exponential) cutoff (PLC). In the former case, the photon break energy is found at E-br,E-gamma = 551(-30)(+70) GeV, slightly higher than what is usually observed in the HE/VHE gamma-ray emitting middle-aged SNRs known to be interacting with molecular clouds. In the latter case. the exponential (respectively sub-exponential) cutoff energy is measured at E-cat,E-gamma = 1.4(-0.55)(+1.6) (respectively 0.35(-0.21)(+0.75)) TeV. A pion decay process resulting from the interaction of the accelerated protons and nuclei with the dense surrounding medium is clearly the preferred scenario to explain the gamma-ray emission. The BPL with a spectral steepening of 0.5-1 and the PLC provide equally good fits to the data. The product or the average gas density and the total energy content of accelerated protons and nuclei amounts to nu W-p similar to 5 x 10(51) erg cm(-3)
The gamma-ray spectrum of the low-frequency-peaked BL Lac (LBL) object AP Librae is studied, following the discovery of very-high-energy (VHE; E > 100 GeV) gamma-ray emission up to the TeV range by the H.E.S.S. experiment. Thismakes AP Librae one of the few VHE emitters of the LBL type. The measured spectrum yields a flux of (8.8 +/- 1.5(stat) +/- 1.8(sys)) x 10(-12) cm(-2) s(-1) above 130 GeV and a spectral index of Gamma = 2.65 +/- 0.19(stat) +/- 0.20(sys). This study also makes use of Fermi-LAT observations in the high energy (HE, E > 100 MeV) range, providing the longest continuous light curve (5 years) ever published on this source. The source underwent a flaring event between MJD 56 306-56 376 in the HE range, with a flux increase of a factor of 3.5 in the 14 day bin light curve and no significant variation in spectral shape with respect to the low-flux state. While the H.E.S.S. and (low state) Fermi-LAT fluxes are in good agreement where they overlap, a spectral curvature between the steep VHE spectrum and the Fermi-LAT spectrum is observed. The maximum of the gamma-ray emission in the spectral energy distribution is located below the GeV energy range.
The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known, N 157B; the radio-loud supernova remnant N 132D; and the largest nonthermal x-ray shell, the superbubble 30 Dor C. The unique object SN 1987A is, unexpectedly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a g-ray source population in an external galaxy and provide via 30 Dor C the unambiguous detection of g-ray emission from a superbubble.
Context. RX J1713.7-3946 is the brightest shell-type supernova remnant (SNR) of the TeV gamma-ray sky. Earlier Fermi-LAT results on low energy gamma-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent with a hadronic origin.
Aims. We update the GeV-band spectra using improved estimates for the diffuse Galactic gamma-ray emission and more than double the volume of data. We further investigate the viability of hadronic emission models for RX J1713.7-3946.
Methods. We produced a high-resolution map of the diffuse Galactic gamma-ray background corrected for the HI self-absorption and used it in the analysis of more than five years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport equation to calculate the acceleration and propagation of cosmic rays in SNR. We then determined spectra of hadronic gamma-ray emission from RX J1713.7-3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux variations that would allow us to test the model with observations.
Results. We find that RX J1713.7-3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphology is best described by that seen in the TeV band. The measured spectrum of RX J1713.7-3946 is hard with index gamma = 1.53 +/- 0.07, and the integral flux above 500 MeV is F = (5 : 5 +/- 1 : 1) x 10(-9) photons cm(-2) s(-1). We demonstrate that scenarios based on hadronic emission from the cosmic-ray precursor region are acceptable for RX J1713.7-3946, and we predict a secular flux increase at a few hundred GeV at the level of around 15% over ten years, which may be detectable with the upcoming Cherenkov Telescope Array (CTA) observatory.
Particle diffusion and localized acceleration in inhomogeneous AGN jets - I. Steady-state spectra
(2015)
We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton (SSC) emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axisymmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard. gamma-ray spectra observed in some blazars.
We present results from VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560 GeV, is well described by a power law with a spectral index of 4.33 +/- 0.09. The time-averaged integral flux above 200 GeV measured for this period was (1.69 +/- 0.06) x 10(-11) photons cm(-2) s(-1), corresponding to 6.9% of the Crab Nebula flux. We also present the combined gamma-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100 MeV to 560 GeV. The data are well fit by a power law with an exponential cutoff at 101.9 +/- 3.2 GeV. The origin of the cutoff could be intrinsic to PG 1553+113 or be due to the gamma-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of z > 0.395 based on optical/UV observations of PG 1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of z <= 0.62. A strongly elevated mean flux of (2.50 +/- 0.14) x10(-11) photons cm(-2) s(-1) (10.3% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as (4.44 +/- 0.71) x10(-11) photons cm(-2) s(-1) (18.3% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a chi(2) probability for a steady flux of 0.03%.
An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated on-off observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of similar to 9 h of on-off observations. Upper limits on the velocity averaged cross section, <sigma upsilon >, for the annihilation of dark matter particles with masses in the range of similar to 300 GeV to similar to 10 TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of <sigma upsilon > that are larger than 3 x 10(-24) cm(3)/s are excluded for dark matter particles with masses between similar to 1 and similar to 4 TeV at 95% C.L. if the radius of the central dark matter density core does not exceed 500 pc. This is the strongest constraint that is derived on <sigma upsilon > for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.
There is a recurring question in solar physics regarding whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Other sources of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net versus neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both direct and return currents, (2) induce very weak compression currents-not observed in 2.5D-in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line-a special condition that is rarely observed. We conclude that. photospheric flows,. as magnetic flux emergence, can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields that possess a net coronal current.
The majority of the matter in the universe is still unidentified and under investigation by both direct and indirect means. Many experiments searching for the recoil of dark-matter particles off target nuclei in underground laboratories have established increasingly strong constraints on the mass and scattering cross sections of weakly interacting particles, and some have even seen hints at a possible signal. Other experiments search for a possible mixing of photons with light scalar or pseudo-scalar particles that could also constitute dark matter. Furthermore, annihilation or decay of dark matter can contribute to charged cosmic rays, photons at all energies, and neutrinos. Many existing and future ground-based and satellite experiments are sensitive to such signals. Finally, data from the Large Hadron Collider at CERN are scrutinized for missing energy as a signature of new weakly interacting particles that may be related to dark matter. In this review article we summarize the status of the field with an emphasis on the complementarity between direct detection in dedicated laboratory experiments, indirect detection in the cosmic radiation, and searches at particle accelerators. (C) 2015 Elsevier B.V. All rights reserved.
Re-observations with the HESS telescope array of the very high-energy (VHE) source HESS J1018-589A that is coincident with the Fermi-LAT gamma-ray binary 1FGL J1018.6-5856 have resulted in a source detection significance of more than 9 sigma and the detection of variability (chi(2)/nu of 238.3/155) in the emitted gamma-ray flux. This variability confirms the association of HESS J1018-589A with the high-energy gamma-ray binary detected Fermi-LAT and also confirms the point-like source as a new VHE binary system. The spectrum of HESS J1018-589A is best fit with a power-law function with photon index Gamma = 2.20 +/- 0.14(stat) +/- 0.2(sys). Emission is detected up to similar to 20 TeV. The mean differential flux level is (2.9 +/- 0.4) x 10(-13) TeV-1 cm(-2) s(-1) at 1 TeV, equivalent to similar to 1% of the flux from the Crab Nebula at the same energy. Variability is clearly detected the night-by-night light curve. When folded on the orbital period of 16.58 days, the rebinned light curve peaks in phase with the observed X-ray high-energy phaseograms. The fit of the HESS phaseogram to a constant flux provides evidence of periodicity at the level of N-sigma > 3 sigma. The of the VHE phaseogram and measured spectrum suggest a low-inclination, low-eccentricity system with a modest impact from VHE gamma-ray due to pair production (tau less than or similar to 1 at 300 GeV).
We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 2013 April 1 and August 10, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope, Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsahovi, and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a Light Detection and Ranging (LIDAR) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution (SED) between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) show evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton (SSC) model to five simultaneous broadband SEDs. We find that the SSC model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.
Context. Puppis A is an interesting similar to 4 kyr-old supernova remnant (SNR) that shows strong evidence of interaction between the forward shock and a molecular cloud. It has been studied in detail from radio frequencies to high-energy (HE, 0.1-100 GeV) gamma-rays. An analysis of the Fermi-LAT data has shown extended HE gamma-ray emission with a 0.2-100 GeV spectrum exhibiting no significant deviation from a power law, unlike most of the GeV-emitting SNRs known to be interacting with molecular clouds. This makes it a promising target for imaging atmospheric Cherenkov telescopes (IACTs) to probe the gamma-ray emission above 100 GeV.
Aims. Very-high-energy (VHE, E >= 0.1 TeV) gamma-ray emission from Puppis A has been, for the first time, searched for with the High Energy Stereoscopic System (HESS.).
Methods. Stereoscopic imaging of Cherenkov radiation from extensive air showers is used to reconstruct the direction and energy of the incident gamma-rays in order to produce sky images and source spectra. The profile likelihood method is applied to find constraints on the existence of a potential break or cutoff in the photon spectrum.
Results. The analysis of the HESS. data does not reveal any significant emission towards Puppis A. The derived upper limits on the differential photon flux imply that its broadband gamma-ray spectrum must exhibit a spectral break or cutoff. By combining Fermi-LAT and HESS. measurements, the 99% confidence-level upper limits on such a cutoff are found to be 450 and 280 GeV, assuming a power law with a simple exponential and a sub-exponential cutoff, respectively. It is concluded that none of the standard limitations (age, size, radiative losses) on the particle acceleration mechanism, assumed to be continuing at present, can explain the lack of VHE signal. The scenario in which particle acceleration has ceased some time ago is considered as an alternative explanation. The HE/VHE spectrum of Puppis A could then exhibit a break of non-radiative origin (as observed in several other interacting SNRs, albeit at somewhat higher energies), owing to the interaction with dense and neutral material, in particular towards the NE region.
Supernova remnants (SNRs) are among the most important targets for gamma-ray observatories. Being prominent non-thermal sources, they are very likely responsible for the acceleration of the bulk of Galactic cosmic rays (CRS). To firmly establish the SNR paradigm for the origin of cosmic rays, it should be confirmed that protons are indeed accelerated in, and released from, SNRs with the appropriate flux and spectrum. This can be done by detailed theoretical models which account for microphysics of acceleration and various radiation processes of hadrons and leptons. The current generation of Cherenkov telescopes has insufficient sensitivity to constrain theoretical models. A new facility, the Cherenkov Telescope Array (CTA), will have superior capabilities and may finally resolve this long standing issue of high-energy astrophysics. We want to assess the capabilities of CTA to reveal the physics of various types of SNRs in the initial 2000 years of their evolution. During this time, the efficiency to accelerate cosmic rays is highest. We perform time-dependent simulations of the hydrodynamics, the magnetic fields, the cosmic-ray acceleration, and the non-thermal emission for type Ia, Ic and IIP SNRs. We calculate the CTA response to the y-ray emission from these SNRs for various ages and distances, and we perform a realistic analysis of the simulated data. We derive distance limits for the detectability and resolvability of these SNR types at several ages. We test the ability of CTA to reconstruct their morphological and spectral parameters as a function of their distance. Finally, we estimate how well CTA data will constrain the theoretical models. (C) 2014 Elsevier B.V. All rights reserved.
Very high energy (VHE, E > 100 GeV)gamma-ray flaring activity of the high-frequency peaked BL Lac object PG 1553 + 113 has been detected by the H.E.S.S. telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with a hint of intra-night variability. No counterpart of this event has been detected in the Fermi-Large Area Telescope data. This pattern is consistent with VHE gamma(-)ray flaring being caused by the injection of ultrarelativistic particles, emitting.-rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at z = 0.49 +/- 0.04 using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz invariance violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are E-QG,E- 1 > 4.10 x 10(17) GeV and E-QG,E- 2 > 2.10 x 10(10) GeV for linear and quadratic LIV effects, respectively.
Aims. Previous observations with the High Energy Stereoscopic System (H.E.S.S.) have revealed an extended very-high-energy (VHE; E > 100 GeV) gamma-ray source, HESS J1834-087, coincident with the supernova remnant (SNR) W41. The origin of the gamma-ray emission was investigated in more detail with the H.E.S.S. array and the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope.
Methods. The gamma-ray data provided by 61 h of observations with H.E.S.S., and four years with the Fermi LAT were analyzed, covering over five decades in energy from 1.8 GeV up to 30 TeV. The morphology and spectrum of the TeV and GeV sources were studied and multiwavelength data were used to investigate the origin of the gamma-ray emission toward W41.
Results. The TeV source can be modeled with a sum of two components: one point-like and one significantly extended (sigma(TeV) = 0.17 degrees +/- 0.01 degrees), both centered on SNR W41 and exhibiting spectra described by a power law with index Gamma(TeV) similar or equal to 2.6. The GeV source detected with Fermi LAT is extended (sigma(GeV) = 0.15 degrees +/- 0.03 degrees) and morphologically matches the VHE emission. Its spectrum can be described by a power-law model with an index Gamma(GeV) = 2.15 +/- 0.12 and smoothly joins the spectrum of the whole TeV source. A break appears in the gamma-ray spectra around 100 GeV. No pulsations were found in the GeV range.
Conclusions. Two main scenarios are proposed to explain the observed emission: a pulsar wind nebula (PWN) or the interaction of SNR W41 with an associated molecular cloud. X-ray observations suggest the presence of a point-like source (a pulsar candidate) near the center of the remnant and nonthermal X-ray diffuse emission that could arise from the possibly associated PWN. The PWN scenario is supported by the compatible positions of the TeV and GeV sources with the putative pulsar. However, the spectral energy distribution from radio to gamma-rays is reproduced by a one-zone leptonic model only if an excess of low-energy electrons is injected following a Maxwellian distribution by a pulsar with a high spin-down power (> 10(37) erg s(-1)). This additional low-energy component is not needed if we consider that the point-like TeV source is unrelated to the extended GeV and TeV sources. The interacting SNR scenario is supported by the spatial coincidence between the gamma-ray sources, the detection of OH (1720 MHz) maser lines, and the hadronic modeling.
We discuss our most recent findings on the diffuse X-ray emission within Wolf-Rayet (WR) nebulae. The best-quality X-ray observations of these objects are those performed by XMM- Newton and Chandra towards S 308, NGC 2359, and NGC 6888. Even though these three WR nebulae might have different formation scenarios, they all share similar characteristics: i) the main plasma temperatures of the X-ray-emitting gas is found to be T =[1–2]×^K, ii) the diffuse X-ray emission is confined inside the [O iii] shell, and iii) their X-ray luminosities and electron densities in the 0.3–2.0 keV energy range are LX ≈10^33–10^34 erg s-1 and ne ≈0.1–1 cm^-3 . These properties and the nebular-like abundances of the hot gas suggest mixing and/or thermal conduction is taking an important rôle reducing the temperature of the hot bubble.
Two types of X-ray sources are mostly found in planetary nebulae (PNe): point sources at their central stars and diffuse emission inside hot bubbles. Here we describe these two types of sources based on the most recent observations obtained in the framework of the Chandra Planetary Nebula Survey, ChanPlaNS, an X-ray survey targeting a volume-limited sample of PNe. Diffuse X-ray emission is found preferentially in young PNe with sharp, closed inner
shells. Point sources of X-ray emission at the central stars reveal magnetically active binary companions and shock-in stellar winds.
Detection and Characterization of Wolf-Rayet stars in M81 with GTC/OSIRIS spectra and HST images
(2015)
Here we investigate a sample of young star clusters (YSCs) and other regions of recent star formation with Wolf-Rayet (W-R) features detected in the relatively nearby spiral galaxy M81 by analysing long-slit (LS) and Multi-Object Spectroscopy (MOS) spectra obtained with the OSIRIS instrument at the 10.4-m Gran Telescopio Canarias (GTC). We take advantage of the synergy between GTC spectra and Hubble Space Telescope (HST) images to also reveal their spatial localization and the environments hosting these stars. We finally discuss and comment on the next steps of our study.
Under the assumption of spherical symmetry, the run of intensity with impact parameter for a spatially resolved and optically thin bubble can be inverted for an "effective emissivity" as a function of radius. The effective emissivity takes into account instrumental sensitivity and even interstellar absorption. This work was supported by a grant from NASA (G03-14008X).
We present the XMM-Newton discovery of X-ray emission from the planetary nebula (PN) A78, the second born-again PN detected in X-rays apart from A30. These two PNe share similar spectral and morphological characteristics: they harbor diffuse soft X-ray emission associated with the interaction between the H-poor ejecta and the current fast stellar wind and a point-like source at the position of the central star (CSPN). We present the spectral analysis of the CSPN, using for the first time an NLTE code for expanding atmospheres that takes line blanketing into account for the UV and optical spectra. The wind abundances are used for the X-ray spectral analysis of the CSPN and the diffuse emission. The X-ray emission from the CSPN in A78 can be modeled by a single C VI emission line, while the X-ray emission from its diffuse component is better described by an optically thin plasma emission model with a temperature of kT = 0.088 keV (T approximate to 1.0 x 10(6) K). We estimate X-ray luminosities in the 0.2-2.0 keV energy band of L-X,L-CSPN =(1.2 +/- 0.3) x 10(31) erg s(-1) and L-X,L-DIFF =(9.2 +/- 2.3) x 10(30) erg s(-1) for the CSPN and diffuse components, respectively.
We study distributed-order time fractional diffusion equations characterized by multifractal memory kernels, in contrast to the simple power-law kernel of common time fractional diffusion equations. Based on the physical approach to anomalous diffusion provided by the seminal Scher-Montroll-Weiss continuous time random walk, we analyze both natural and modified-form distributed-order time fractional diffusion equations and compare the two approaches. The mean squared displacement is obtained and its limiting behavior analyzed. We derive the connection between the Wiener process, described by the conventional Langevin equation and the dynamics encoded by the distributed-order time fractional diffusion equation in terms of a generalized subordination of time. A detailed analysis of the multifractal properties of distributed-order diffusion equations is provided.
We consider anomalous stochastic processes based on the renewal continuous time random walk model with different forms for the probability density of waiting times between individual jumps. In the corresponding continuum limit we derive the generalized diffusion and Fokker-Planck-Smoluchowski equations with the corresponding memory kernels. We calculate the qth order moments in the unbiased and biased cases, and demonstrate that the generalized Einstein relation for the considered dynamics remains valid. The relaxation of modes in the case of an external harmonic potential and the convergence of the mean squared displacement to the thermal plateau are analyzed.
The electrical conductivity of organic semiconductors can be enhanced by orders of magnitude via doping with strong molecular electron acceptors or donors. Ground-state integer charge transfer and charge-transfer complex formation between organic semiconductors and molecular dopants have been suggested as the microscopic mechanisms causing these profound changes in electrical materials properties. Here, we study charge-transfer interactions between the common molecular p-dopant 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane and a systematic series of thiophene-based copolymers by a combination of spectroscopic techniques and electrical measurements. Subtle variations in chemical structure are seen to significantly impact the nature of the charge-transfer species and the efficiency of the doping process, underlining the need for a more detailed understanding of the microscopic doping mechanism in organic semiconductors to reliably guide targeted chemical design.
Organic bulk heterojunction (BHJ) solar cells based on polymer:fullerene blends are a promising alternative for a low-cost solar energy conversion. Despite significant improvements of the power conversion efficiency in recent years, the fundamental working principles of these devices are yet not fully understood. In general, the current output of organic solar cells is determined by the generation of free charge carriers upon light absorption and their transport to the electrodes in competition to the loss of charge carriers due to recombination.
The object of this thesis is to provide a comprehensive understanding of the dynamic processes and physical parameters determining the performance. A new approach for analyzing the characteristic current-voltage output was developed comprising the experimental determination of the efficiencies of charge carrier generation, recombination and transport, combined with numerical device simulations.
Central issues at the beginning of this work were the influence of an electric field on the free carrier generation process and the contribution of generation, recombination and transport to the current-voltage characteristics. An elegant way to directly measure the field dependence of the free carrier generation is the Time Delayed Collection Field (TDCF) method. In TDCF charge carriers are generated by a short laser pulse and subsequently extracted by a defined rectangular voltage pulse. A new setup was established with an improved time resolution compared to former reports in literature. It was found that charge generation is in general independent of the electric field, in contrast to the current view in literature and opposed to the expectations of the Braun-Onsager model that was commonly used to describe the charge generation process. Even in cases where the charge generation was found to be field-dependend, numerical modelling showed that this field-dependence is in general not capable to account for the voltage dependence of the photocurrent. This highlights the importance of efficient charge extraction in competition to non-geminate recombination, which is the second objective of the thesis.
Therefore, two different techniques were combined to characterize the dynamics and efficiency of non-geminate recombination under device-relevant conditions. One new approach is to perform TDCF measurements with increasing delay between generation and extraction of charges. Thus, TDCF was used for the first time to measure charge carrier generation, recombination and transport with the same experimental setup. This excludes experimental errors due to different measurement and preparation conditions and demonstrates the strength of this technique. An analytic model for the description of TDCF transients was developed and revealed the experimental conditions for which reliable results can be obtained. In particular, it turned out that the $RC$ time of the setup which is mainly given by the sample geometry has a significant influence on the shape of the transients which has to be considered for correct data analysis.
Secondly, a complementary method was applied to characterize charge carrier recombination under steady state bias and illumination, i.e. under realistic operating conditions. This approach relies on the precise determination of the steady state carrier densities established in the active layer. It turned out that current techniques were not sufficient to measure carrier densities with the necessary accuracy. Therefore, a new technique {Bias Assisted Charge Extraction} (BACE) was developed. Here, the charge carriers photogenerated under steady state illumination are extracted by applying a high reverse bias. The accelerated extraction compared to conventional charge extraction minimizes losses through non-geminate recombination and trapping during extraction. By performing numerical device simulations under steady state, conditions were established under which quantitative information on the dynamics can be retrieved from BACE measurements.
The applied experimental techniques allowed to sensitively analyse and quantify geminate and non-geminate recombination losses along with charge transport in organic solar cells. A full analysis was exemplarily demonstrated for two prominent polymer-fullerene blends.
The model system P3HT:PCBM spincast from chloroform (as prepared) exhibits poor power conversion efficiencies (PCE) on the order of 0.5%, mainly caused by low fill factors (FF) and currents. It could be shown that the performance of these devices is limited by the hole transport and large bimolecular recombination (BMR) losses, while geminate recombination losses are insignificant. The low polymer crystallinity and poor interconnection between the polymer and fullerene domains leads to a hole mobility of the order of 10^-7 cm^2/Vs which is several orders of magnitude lower than the electron mobility in these devices. The concomitant build up of space charge hinders extraction of both electrons and holes and promotes bimolecular recombination losses.
Thermal annealing of P3HT:PCBM blends directly after spin coating improves crystallinity and interconnection of the polymer and the fullerene phase and results in comparatively high electron and hole mobilities in the order of 10^-3 cm^2/Vs and 10^-4 cm^2/Vs, respectively. In addition, a coarsening of the domain sizes leads to a reduction of the BMR by one order of magnitude. High charge carrier mobilities and low recombination losses result in comparatively high FF (>65%) and short circuit current (J_SC ≈ 10 mA/cm^2). The overall device performance (PCE ≈ 4%) is only limited by a rather low spectral overlap of absorption and solar emission and a small V_OC, given by the energetics of the P3HT.
From this point of view the combination of the low bandgap polymer PTB7 with PCBM is a promising approach. In BHJ solar cells, this polymer leads to a higher V_OC due to optimized energetics with PCBM. However, the J_SC in these (unoptimized) devices is similar to the J_SC in the optimized blend with P3HT and the FF is rather low (≈ 50%). It turned out that the unoptimized PTB7:PCBM blends suffer from high BMR, a low electron mobility of the order of 10^-5 cm^2/Vs and geminate recombination losses due to field dependent charge carrier generation.
The use of the solvent additive DIO optimizes the blend morphology, mainly by suppressing the formation of very large fullerene domains and by forming a more uniform structure of well interconnected donor and acceptor domains of the order of a few nanometers. Our analysis shows that this results in an increase of the electron mobility by about one order of magnitude (3 x 10^-4 cm^2/Vs), while BMR and geminate recombination losses are significantly reduced. In total these effects improve the J_SC (≈ 17 mA/cm^2) and the FF (> 70%). In 2012 this polymer/fullerene combination resulted in a record PCE for a single junction OSC of 9.2%.
Remarkably, the numerical device simulations revealed that the specific shape of the J-V characteristics depends very sensitively to the variation of not only one, but all dynamic parameters. On the one hand this proves that the experimentally determined parameters, if leading to a good match between simulated and measured J-V curves, are realistic and reliable. On the other hand it also emphasizes the importance to consider all involved dynamic quantities, namely charge carrier generation, geminate and non-geminate recombination as well as electron and hole mobilities. The measurement or investigation of only a subset of these parameters as frequently found in literature will lead to an incomplete picture and possibly to misleading conclusions.
Importantly, the comparison of the numerical device simulation employing the measured parameters and the experimental $J-V$ characteristics allows to identify loss channels and limitations of OSC. For example, it turned out that inefficient extraction of charge carriers is a criticical limitation factor that is often disobeyed. However, efficient and fast transport of charges becomes more and more important with the development of new low bandgap materials with very high internal quantum efficiencies. Likewise, due to moderate charge carrier mobilities, the active layer thicknesses of current high-performance devices are usually limited to around 100 nm. However, larger layer thicknesses would be more favourable with respect to higher current output and robustness of production. Newly designed donor materials should therefore at best show a high tendency to form crystalline structures, as observed in P3HT, combined with the optimized energetics and quantum efficiency of, for example, PTB7.
Among the parameters that characterize a solar cell and define its power-conversion efficiency, the fill factor is the least well understood, making targeted improvements difficult. Here we quantify the competition between charge extraction and recombination by using a single parameter theta, and we demonstrate that this parameter is directly related to the fill factor of many different bulk-heterojunction solar cells. Our finding is supported by experimental measurements on 15 different donor: acceptor combinations, as well as by drift-diffusion simulations of organic solar cells in which charge-carrier mobilities, recombination rate, light intensity, energy levels and active-layer thickness are all varied over wide ranges to reproduce typical experimental conditions. The results unify the fill factors of several very different donor: acceptor combinations and give insight into why fill factors change so much with thickness, light intensity and materials properties. To achieve fill factors larger than 0.8 requires further improvements in charge transport while reducing recombination.
We investigate hybrid charge transfer states (HCTS) at the planar interface between a-NPD and ZnO by spectrally resolved electroluminescence (EL) and external quantum efficiency (EQE) measurements. Radiative decay of HCTSs is proven by distinct emission peaks in the EL spectra of such bilayer devices in the NIR at energies well below the bulk a-NPD or ZnO emission. The EQE spectra display low energy contributions clearly red-shifted with respect to the a-NPD photocurrent and partially overlapping with the EL emission. Tuning of the energy gap between the ZnO conduction band and a-NPD HOMO level (E-int) was achieved by modifying the ZnO surface with self-assembled monolayers based on phosphonic acids. We find a linear dependence of the peak position of the NIR EL on E-int, which unambiguously attributes the origin of this emission to radiative recombination between an electron on the ZnO and a hole on a-NPD. In accordance with this interpretation, we find a strictly linear relation between the open-circuit voltage and the energy of the charge state for such hybrid organicinorganic interfaces.
The IR-based method of infrared transition moment orientational analysis (IR-TMOA) is employed to unravel molecular order in thin layers of the semiconducting polymer poly[N,N'-bis(2-octyldodecyl),-1,4,5,8-naphthalene-diimide-2,6-diyl]-alt-5-5'-(2,2'-bithiophene) (P(NDI2OD-T2)). Structure-specific vibrational bands are analyzed in dependence On polarization and inclination of the sample-With respect to the optical axis. By that the molecular Order parameter tensor for the respective molecular moieties with regard to the sample: coordinate system is deduced. Making use of the specificity of the IR spectral range, we are able to determine separately the orientation of atomistic planes defined through the naphthalenediimide (NDI) and bithiophene (T2) units relative to the substrate, and hence, relative to each other. A pronounced solvent effect is observed While chlorobenzene causes the T2 planes to align preferentially parallel to the substrate at an angle of 29 degrees, using a 1:1 chloronaphthalene:xylene mixture results in a reorientation of the T2 units from a face on into an edge on arrangement. In contrast the NDI unit remains unaffected. Additionally, for both solvents evidence is observed for the aggregation of chains in accord With recently published results obtained by UV-vis absorption spectroscopy.
An approach is presented to modify the work function of solution-processed sol-gel derived zinc oxide (ZnO) over an exceptionally wide range of more than 2.3 eV. This approach relies on the formation of dense and homogeneous self-assembled monolayers based on phosphonic acids with different dipole moments. This allows us to apply ZnO as charge selective bottom electrodes in either regular or inverted solar cell structures, using poly(3-hexylthiophene): phenyl-C71-butyric acid methyl ester as the active layer. These devices compete with or even surpass the performance of the reference on indium tin oxide/poly(3,4-ethylenedioxythiophene) polystyrene sulfonate. Our findings highlight the potential of properly modified ZnO as electron or hole extracting electrodes in hybrid optoelectronic devices. (C) 2015 AIP Publishing LLC.
Blending the conjugated polymer poly(3-hexylthiophene) (P3HT) with the insulating electret polystyrene (PS), we show that the threshold voltage V-t of organic field-effect transistors (OFETs) can be easily and reversely tuned by applying a gate bias stress at 130 degrees C. It is proposed that this phenomenon is caused by thermally activated charge injection from P3HT into PS matrix, and that this charge is immobilized within the PS matrix after cooling down to room temperature. Therefore, room-temperature hysteresis-free FETs with desired V-t can be easily achieved. The approach is applied to reversely tune the OFET mode of operation from accumulation to depletion, and to build inverters. (C) 2015 AIP Publishing LLC.
When arranged in a proper nanoaggregate architecture, gold nanoparticles can offer controllable plasmon-related absorption/scattering, yielding distinct color effects that depend critically on the relative orientation and distance between nanoparticle constituents. Herein, we report on the implementation of novel plasmonic nanoarchitectures based on complexes between gold nanoparticles and an azobenzene-modified cationic surfactant that can exhibit a light-tunable plasmonic response. The formation of such complexes becomes possible through the use of strongly negatively charged bare gold nanoparticles (similar to 10-nm diameter) prepared by the method of laser ablation in deionized water. Driven by electrostatic interactions, the cationic surfactant molecules attach and form a shell around the negatively charged nanoparticles, resulting in neutralization of the particle charge or even overcompensation beyond which the nanoparticles become positively charged. At low and high surfactant concentrations, Au nanoparticles are negatively and positively charged, respectively, and are represented by single species due to electric repulsion effects having absorption peaks around 523-527 nm, whereas at intermediate concentrations, the Au nanoparticles become neutral, forming nanoscale 100-nm clusterlike aggregates and exhibiting an additional absorption peak at gimel > 600 nm and a visible change in the color of the solution from red to blue. Because of the presence of the photosensitive azobenzene unit in the surfactant tail that undergoes trans-to-cis isomerization under irradiation with UV light, we then demonstrate a light-controlled nanoclustering of nanoparticles, yielding a switch in the plasmonic absorption band and a related change in the solution color. The formed hybrid architectures with a light-controlled plasmonic response could be important for a variety of tasks, including biomedical, surface-enhanced Raman spectroscopy (SERS), data transmission, and storage applications.
Time-delayed collection field (TDCF), bias-assisted charge extraction (BACE), and space charge-limited current (SCLC) measurements are combined with complete numerical device simulations to unveil the effect of the solvent additive 1,8-diiodooctane (DIO) on the performance of PTB7:PCBM bulk heterojunction solar cells. DIO is shown to increase the charge generation rate, reduce geminate and bimolecular recombination, and increase the electron mobility. In total, the reduction of loss currents by processing with the additive raises the power conversion efficiency of the PTB7:PCBM blend by a factor of almost three. The lower generation rates and higher geminate recombination losses in devices without DIO are consistent with a blend morphology comprising large fullerene clusters embedded within a PTB7-rich matrix, while the low electron mobility suggests that these fullerene clusters are themselves composed of smaller pure fullerene aggregates separated by disordered areas. Our device simulations show unambiguously that the effect of the additive on the shape of the currentvoltage curve (J-V) cannot be ascribed to the variation of only the mobility, the recombination, or the field dependence of generation. It is only when the changes of all three parameters are taken into account that the simulation matches the experimental J-V characteristics under all illumination conditions and for a wide range of voltages.
PCPDTBT, a marginally crystallizable polymer, is crystallized into a new crystal structure using solvent-vapor annealing. Highly ordered areas with three different polymer-chain orientations are identified using TEM/ED, GIWAXS, and polarized Raman spectroscopy. The optical and structural properties differ significantly from films prepared by standard device preparation protocols. Bilayer solar cells, however, show similar performance.
This work elucidates the impact of charge transport on the photovoltaic properties of organic solar cells. Here we show that the analysis of current-voltage curves of organic solar cells under illumination with the Shockley equation results in values for ideality factor, photo-current and parallel resistance, which lack physical meaning. Drift-diffusion simulations for a wide range of charge-carrier mobilities and illumination intensities reveal significant carrier accumulation caused by poor transport properties, which is not included in the Shockley equation. As a consequence, the separation of the quasi Fermi levels in the organic photoactive layer (internal voltage) differs substantially from the external voltage for almost all conditions. We present a new analytical model, which considers carrier transport explicitly. The model shows excellent agreement with full drift-diffusion simulations over a wide range of mobilities and illumination intensities, making it suitable for realistic efficiency predictions for organic solar cells.
Hybrid multijunction solar cells comprising hydrogenated amorphous silicon and an organic bulk heterojunction are presented, reaching 11.7% power conversion efficiency. The benefits of merging inorganic and organic subcells are pointed out, the optimization of the cells, including optical modeling predictions and tuning of the recombination contact are described, and an outlook of this technique is given.