Refine
Has Fulltext
- yes (108) (remove)
Year of publication
- 2015 (108) (remove)
Document Type
- Article (83)
- Doctoral Thesis (15)
- Postprint (7)
- Conference Proceeding (1)
- Habilitation Thesis (1)
- Master's Thesis (1)
Language
- English (108)
Keywords
- regional climate model (2)
- Aerosole (1)
- Aktuation (1)
- Arktis (1)
- Asteroiden (1)
- Astronomie (1)
- Astrophysik (1)
- Biomechanik (1)
- Biophysik (1)
- Bogoliubov theory (1)
Institute
- Institut für Physik und Astronomie (108) (remove)
Optical frequency combs (OFC) constitute an array of phase-correlated equidistant spectral lines with nearly equal intensities over a broad spectral range. The adaptations of combs generated in mode-locked lasers proved to be highly efficient for the calibration of high-resolution (resolving power > 50000) astronomical spectrographs. The observation of different galaxy structures or the studies of the Milky Way are done using instruments in the low- and medium resolution range. To such instruments belong, for instance, the Multi Unit Spectroscopic Explorer (MUSE) being developed for the Very Large Telescope (VLT) of the European Southern Observatory (ESO) and the 4-metre Multi-Object Spectroscopic Telescope (4MOST) being in development for the ESO VISTA 4.1 m Telescope. The existing adaptations of OFC from mode-locked lasers are not resolvable by these instruments.
Within this work, a fibre-based approach for generation of OFC specifically in the low- and medium resolution range is studied numerically. This approach consists of three optical fibres that are fed by two equally intense continuous-wave (CW) lasers. The first fibre is a conventional single-mode fibre, the second one is a suitably pumped amplifying Erbium-doped fibre with anomalous dispersion, and the third one is a low-dispersion highly nonlinear optical fibre. The evolution of a frequency comb in this system is governed by the following processes: as the two initial CW-laser waves with different frequencies propagate through the first fibre, they generate an initial comb via a cascade of four-wave mixing processes. The frequency components of the comb are phase-correlated with the original laser lines and have a frequency spacing that is equal to the initial laser frequency separation (LFS), i.e. the difference in the laser frequencies. In the time domain, a train of pre-compressed pulses with widths of a few pico-seconds arises out of the initial bichromatic deeply-modulated cosine-wave. These pulses undergo strong compression in the subsequent amplifying Erbium-doped fibre: sub-100 fs pulses with broad OFC spectra are formed. In the following low-dispersion highly nonlinear fibre, the OFC experience a further broadening and the intensity of the comb lines are fairly equalised. This approach was mathematically modelled by means of a Generalised Nonlinear Schrödinger Equation (GNLS) that contains terms describing the nonlinear optical Kerr effect, the delayed Raman response, the pulse self-steepening, and the linear optical losses as well as the wavelength-dependent Erbium gain profile for the second fibre. The initial condition equation being a deeply-modulated cosine-wave mimics the radiation of the two initial CW lasers. The numerical studies are performed with the help of Matlab scripts that were specifically developed for the integration of the GNLS and the initial condition according to the proposed approach for the OFC generation. The scripts are based on the Fourth-Order Runge-Kutta in the Interaction Picture Method (RK4IP) in combination with the local error method.
This work includes the studies and results on the length optimisation of the first and the second fibre depending on different values of the group-velocity dispersion of the first fibre. Such length optimisation studies are necessary because the OFC have the biggest possible broadband and exhibit a low level of noise exactly at the optimum lengths. Further, the optical pulse build-up in the first and the second fibre was studied by means of the numerical technique called Soliton Radiation Beat Analysis (SRBA). It was shown that a common soliton crystal state is formed in the first fibre for low laser input powers. The soliton crystal continuously dissolves into separated optical solitons as the input power increases. The pulse formation in the second fibre is critically dependent on the features of the pulses formed in the first fibre. I showed that, for low input powers, an adiabatic soliton compression delivering low-noise OFC occurs in the second fibre. At high input powers, the pulses in the first fibre have more complicated structures which leads to the pulse break-up in the second fibre with a subsequent degradation of the OFC noise performance. The pulse intensity noise studies that were performed within the framework of this thesis allow making statements about the noise performance of an OFC. They showed that the intensity noise of the whole system decreases with the increasing value of LFS.
The high-latitudinal thermospheric processes driven by the solar wind and Interplanetary Magnetic Field (IMF) interaction with the Earth magnetosphere are highly variable parts of the complex dynamic plasma environment, which represent the coupled Magnetosphere – Ionosphere – Thermosphere (MIT) system. The solar wind and IMF interactions transfer energy to the MIT system via reconnection processes at the magnetopause. The Field Aligned Currents (FACs) constitute the energetic links between the magnetosphere and the Earth ionosphere. The MIT system depends on the highly variable solar wind conditions, in particular on changes of the strength and orientation of the IMF.
In my thesis, I perform an investigation on the physical background of the complex MIT system using the global physical - numerical, three-dimensional, time-dependent and self-consistent Upper Atmosphere Model (UAM). This model describes the thermosphere, ionosphere, plasmasphere and inner magnetosphere as well as the electrodynamics of the coupled MIT system for the altitudinal range from 80 (60) km up to the 15 Earth radii.
In the present study, I developed and investigated several variants of the high-latitudinal electrodynamic coupling by including the IMF dependence of FACs into the UAM model. For testing, the various variants were applied to simulations of the coupled MIT system for different seasons, geomagnetic activities, various solar wind and IMF conditions. Additionally, these variants of the theoretical model with the IMF dependence were compared with global empirical models. The modelling results for the most important thermospheric parameters like neutral wind and mass density were compared with satellite measurements. The variants of the UAM model with IMF dependence show a good agreement with the satellite observations. In comparison with the empirical models, the improved variants of the UAM model reproduce a more realistic meso-scale structures and dynamics of the coupled MIT system than the empirical models, in particular at high latitudes. The new configurations of the UAM model with IMF dependence contribute to the improvement of space weather prediction.
Most of the baryonic matter in the Universe resides in a diffuse gaseous phase in-between galaxies consisting mostly of hydrogen and helium. This intergalactic medium (IGM) is distributed in large-scale filaments as part of the overall cosmic web. The luminous extragalactic objects that we can observe today, such as galaxies and quasars, are surrounded by the IGM in the most dense regions within the cosmic web. The radiation of these objects contributes to the so-called ultraviolet background (UVB) which keeps the IGM highly ionized ever since the epoch of reionization.
Measuring the amount of absorption due to intergalactic neutral hydrogen (HI) against extragalactic background sources is a very useful tool to constrain the energy input of ionizing sources into the IGM. Observations suggest that the HI Lyman-alpha effective optical depth, τ_eff, decreases with decreasing redshift, which is primarily due to the expansion of the Universe. However, some studies find a smaller value of the effective optical depth than expected at the specific redshift z~3.2, possibly related to the complete reionization of helium in the IGM and a hardening of the UVB. The detection and possible cause of a decrease in τ_eff at z~3.2 is controversially debated in the literature and the observed features need further explanation.
To better understand the properties of the mean absorption at high redshift and to provide an answer for whether the detection of a τ_eff feature is real we study 13 high-resolution, high signal-to-noise ratio quasar spectra observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) at the Very Large Telescope (VLT). The redshift evolution of the effective optical depth, τ_eff(z), is measured in the redshift range 2.7≤z≤3.6. The influence of metal absorption features is removed by performing a comprehensive absorption-line-fitting procedure.
In the first part of the thesis, a line-parameter analysis of the column density, N, and Doppler parameter, b, of ≈7500 individually fitted absorption lines is performed. The results are in good agreement with findings from previous surveys.
The second (main) part of this thesis deals with the analysis of the redshift evolution of the effective optical depth. The τ_eff measurements vary around the empirical power law τ_eff(z)~(1+z)^(γ+1) with γ=2.09±0.52. The same analysis as for the observed spectra is performed on synthetic absorption spectra. From a comparison between observed and synthetic spectral data it can be inferred that the uncertainties of the τ_eff values are likely underestimated and that the scatter is probably caused by high-column-density absorbers with column densities in the range 15≤logN≤17. In the real Universe, such absorbers are rarely observed, however. Hence, the difference in τ_eff from different observational data sets and absorption studies is most likely caused by cosmic variance. If, alternatively, the disagreement between such data is a result of an too optimistic estimate of the (systematic) errors, it is also possible that all τ_eff measurements agree with a smooth evolution within the investigated redshift range. To explore in detail the different analysis techniques of previous studies an extensive literature comparison to the results of this work is presented in this thesis.
Although a final explanation for the occurrence of the τ_eff deviation in different studies at z~3.2 cannot be given here, our study, which represents the most detailed line-fitting analysis of its kind performed at the investigated redshifts so far, represents another important benchmark for the characterization of the HI Ly-alpha effective optical depth at high redshift and its indicated unusual behavior at z~3.2.
Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.
Semi-empirical sea-level models (SEMs) exploit physically motivated empirical relationships between global sea level and certain drivers, in the following global mean temperature. This model class evolved as a supplement to process-based models (Rahmstorf (2007)) which were unable to fully represent all relevant processes. They thus failed to capture past sea-level change (Rahmstorf et al. (2012)) and were thought likely to underestimate future sea-level rise. Semi-empirical models were found to be a fast and useful tool for exploring the uncertainties in future sea-level rise, consistently giving significantly higher projections than process-based models.
In the following different aspects of semi-empirical sea-level modelling have been studied. Models were first validated using various data sets of global sea level and temperature. SEMs were then used on the glacier contribution to sea level, and to infer past global temperature from sea-level data via inverse modelling. Periods studied encompass the instrumental period, covered by tide gauges (starting 1700 CE (Common Era) in Amsterdam) and satellites (first launched in 1992 CE), the era from 1000 BCE (before CE) to present, and the full length of the Holocene (using proxy data). Accordingly different data, model formulations and implementations have been used. It could be shown in Bittermann et al. (2013) that SEMs correctly predict 20th century sea-level when calibrated with data until 1900 CE. SEMs also turned out to give better predictions than the Intergovernmental Panel on Climate Change (IPCC) 4th assessment report (AR4, IPCC (2007)) models, for the period from 1961–2003 CE.
With the first multi-proxy reconstruction of global sea-level as input, estimate of the human-induced component of modern sea-level change and projections of future sea-level rise were calculated (Kopp et al. (2016)). It turned out with 90% confidence that more than 40 % of the observed 20th century sea-level rise is indeed anthropogenic. With the new semi-empirical and IPCC (2013) 5th assessment report (AR5) projections the gap between SEM and process-based model projections closes, giving higher credibility to both. Combining all scenarios, from strong mitigation to business as usual, a global sea-level rise of 28–131 cm relative to 2000 CE, is projected with 90% confidence. The decision for a low carbon pathway could halve the expected global sea-level rise by 2100 CE.
Present day temperature and thus sea level are driven by the globally acting greenhouse-gas forcing. Unlike that, the Milankovich forcing, acting on Holocene timescales, results mainly in a northern-hemisphere temperature change. Therefore a semi-empirical model can be driven with northernhemisphere temperatures, which makes it possible to model the main subcomponent of sea-level change over this period. It showed that an additional positive constant rate of the order of the estimated Antarctic sea-level contribution is then required to explain the sea-level evolution over the Holocene. Thus the global sea level, following the climatic optimum, can be interpreted as the sum of a temperature induced sea-level drop and a positive long-term contribution, likely an ongoing response to deglaciation coming from Antarctica.
We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T−h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.
The lives of more than 1/6 th of the world population is directly affected by the caprices of the South Asian summer monsoon rainfall. India receives around 78 % of the annual precipitation during the June-September months, the summer monsoon season of South Asia. But, the monsoon circulation is not consistent throughout the entire summer season. Episodes of heavy rainfall (active periods) and low rainfall (break periods) are inherent to the intraseasonal variability of the South Asian summer monsoon. Extended breaks or long-lasting dryness can result in droughts and hence trigger crop failures and in turn famines. Furthermore, India's electricity generation from renewable sources (wind and hydro-power), which is increasingly important in order to satisfy the rapidly rising demand for energy, is highly reliant on the prevailing meteorology. The major drought years 2002 and 2009 for the Indian summer monsoon during the last decades, which are results of the occurrence of multiple extended breaks, emphasise exemplary that the understanding of the monsoon system and its intraseasonal variation is of greatest importance. Although, numerous studies based on observations, reanalysis data and global model simulations have been carried out with the focus on monsoon active and break phases over India, the understanding of the monsoon intraseasonal variability is only in the infancy stage. Regional climate models could benefit the comprehension of monsoon breaks by its resolution advantage.
This study investigates moist dynamical processes that initiate and maintain breaks during the South Asian summer monsoon using the atmospheric regional climate model HIRHAM5 at a horizontal resolution of 25 km forced by the ECMWF ERA Interim reanalysis for the period 1979-2012. By calculating moisture and moist static energy budgets the various competing mechanisms leading to extended breaks are quantitatively estimated. Advection of dry air from the deserts of western Asia towards central India is the dominant moist dynamical process in initiating extended break conditions over South Asia. Once initiated, the extended breaks are maintained due to many competing mechanisms: (i) the anomalous easterlies at the southern flank of this anticyclonic anomaly weaken the low-level cross-equatorial jet and thus the moisture transport into the monsoon region, (ii) differential radiative heating over the continental and the oceanic tropical convergence zone induces a local Hadley circulation with anomalous rising over the equatorial Indian Ocean and descent over central India, and (iii) a cyclonic response to positive rainfall anomalies over the near-equatorial Indian Ocean amplifies the anomalous easterlies over India and hence contributes to the low-level divergence over central India.
A sensitivity experiment that mimics a scenario of higher atmospheric aerosol concentrations over South Asia addresses a current issue of large uncertainty: the role aerosols play in suppressing monsoon rainfall and hence in triggering breaks. To study the indirect aerosol effects the cloud droplet number concentration was increased to imitate the aerosol's function as cloud condensation nuclei. The sensitivity experiment with altered microphysical cloud properties shows a reduction in the summer monsoon precipitation together with a weakening of the South Asian summer monsoon. Several physical mechanisms are proposed to be responsible for the suppressed monsoon rainfall: (i) according to the first indirect radiative forcing the increase in the number of cloud droplets causes an increase in the cloud reflectivity of solar radiation, leading to a climate cooling over India which in turn reduces the hydrological cycle, (ii) a stabilisation of the troposphere induced by a differential cooling between the surface and the upper troposphere over central India inhibits the growth of deep convective rain clouds, (iii) an increase of the amount of low and mid-level clouds together with a decrease in high-level cloud amount amplify the surface cooling and hence the atmospheric stability, and (iv) dynamical changes of the monsoon manifested as a anomalous anticyclonic circulation over India reduce the moisture transport into the monsoon region. The study suggests that the changes in the total precipitation, which are dominated by changes in the convective precipitation, mainly result from the indirect radiative forcing. Suppression of rainfall due to the direct microphysical effect is found to be negligible over India. Break statistics of the polluted cloud scenario indicate an increase in the occurrence of short breaks (3 days), while the frequency of extended breaks (> 7 days) is clearly not affected. This disproves the hypothesis that more and smaller cloud droplets, caused by a high load of atmospheric aerosols trigger long drought conditions over central India.
This work investigates the influence of the Coriolis force on mass motion related to the Rheasilvia impact basin on asteroid (4) Vesta's southern hemisphere. The giant basin is 500km in diameter, with a centre which nearly coincides with the rotation axis of Vesta. The Rheasilvia basin partially overlaps an earlier, similarly large impact basin, Veneneia.
Mass motion within and in the vicinity of the Rheasilvia basin includes slumping and landslides, which, primarily due to their small linear extents, have not been noticeably affected by the Coriolis force. However, a series of ridges related to the basin exhibit significant curvature, which may record the effect of the Coriolis force on the mass motion which generated them.
In this thesis 32 of these curved ridges, in three geologically distinct regions, were examined. The mass motion velocities from which the ridge curvatures may have resulted during the crater modification stage were investigated. Velocity profiles were derived by fitting inertial circles along the curved ridges and considering both the current and past rotation states of Vesta. An iterative, statistical approach was used, whereby the radii of inertial circles were obtained through repeated fitting to triplets of points across the ridges. The most frequently found radius for each central point was then used for velocity derivation at that point.
The results of the velocity analysis are strongly supportive of a Coriolis force origin for the curved ridges. Derived velocities (29.6 ± 24.6 m/s) generally agree well with previously published predictions from numerical simulations of mass motion during the impact process. Topographical features such as local slope gradient and mass deposition regions on the curved ridges also independently agree with regions in which the calculated mass motion accelerates or decelerates.
Sections of constant acceleration, deceleration and constant velocity are found, showing that mass motion is being governed by varying conditions of topography, regolith structure and friction. Estimates of material properties such as the effective viscosities (1.9-9.0·10⁶ Pa·s) and coefficients of friction (0.02-0.81) are derived from the velocity profile information in these sections. From measured accelerations of mass motions on the crater wall, it is also shown that the crater walls must have been locally steeper at the time of the mass motion.
Together with these novel insights into the state and behaviour of material moving during the modification stage of Rheasilvia's formation, this work represents the first time that the Coriolis Effect on mass motions during crater formation has been shown to result in diagnostic features preserved until today.
The evolution of massive stars in very low metallicity galaxies is less well observationally
constrained than in environments more similar to the Milky Way, M33, or the LMC. We discuss
in this contribution the current state of our program to search for and characterize Wolf-Rayet stars (and other massive emission line stars) in low metallicity galaxies in the Local Volume.
PopIII-star siblings in IZw18 and metal-poor WR galaxies unveiled from integral field spectroscopy
(2015)
Here, we highlight our recent results from the IFS study of Mrk178, the closest metal-poor WR galaxy, and of IZw18, the most metal-poor star-forming galaxy known in the local Universe. The IFS data of Mrk178 show the importance of aperture effects on the search for WR features, and the extent to which physical variations in the ISM properties can be detected. Our IFS data of IZw18 reveal its entire nebular HeIIλ4686-emitting region, and indicate for the very first time that peculiar, hot (nearly) metal-free ionizing stars (called here PopIII-star siblings) might hold the key to the HeII-ionization in IZw18.
The main objective of this work is to investigate the evolution of massive stars, and the interplay between them and the ionized gas for a sample of local metal-poor Wolf-Rayet galaxies.
Optical integral field spectrocopy was used in combination with multi-wavelength radio data.
Combining optical and radio data, we locate Wolf-Rayet stars and supernova remnants across the Wolf-Rayet galaxies to study the spatial correlation between them. This study will shed light on the massive star formation and its feedback, and will help us to better understand
distant star-forming galaxies.
We highlight the basic physics that allows fundamental parameters, such as the effective
temperature, luminosity, abundances, and mass-loss rate, of Wolf-Rayet (W-R) stars to be
determined. Since the temperature deduced from the spectrum of a W-R star is an ionization
temperature, a detailed discussion of the ionization structure of W-R winds, and how it is set, is given. We also provide an overview of line and continuum formation in W-R stars. Mechanisms that contribute to the strength of different emission lines, such as collisional excitation, radiative recombination, dielectronic recombination, and continuum uorescence, are discussed.
For some years now, spectroscopic measurements of massive stars in the amateur domain have been fulfilling professional requirements. Various groups in the northern and southern hemispheres have been established, running successful professional-amateur (ProAm) collaborative campaigns, e.g., on WR, O and B type stars. Today high quality data (echelle and long-slit) are regularly delivered and corresponding results published. Night-to-night long-term observations over months to years open a new opportunity for massive-star research. We introduce recent and ongoing sample campaigns (e.g. ∊ Aur, WR 134, ζ Pup), show respective results and highlight the vast amount of data collected in various data bases. Ultimately it is in the time-dependent domain where amateurs can shine most.
Because most massive stars have been or will be affected by a companion during the course of their evolution, we cannot afford to neglect binaries when discussing the progenitors of supernovae and GRBs. Analyzing linear polarization in the emission lines of close binary systems allows us to probe the structures of these systems' winds and mass flows, making it possible to map the complex morphologies of the mass loss and mass transfer structures that shape their subsequent evolution. In Wolf-Rayet (WR) binaries, line polarization variations with orbital phase distinguish polarimetric signatures arising from lines that scatter near the stars from those that scatter far from the orbital plane. These far-scattering lines may form the basis for a "binary line-effect method" of identifying rapidly rotating WR stars (and hence GRB progenitor candidates) in binary systems.
Before GAIA improves the HIPPARCOS survey, direct determination of the distance via parallax is only possible for γ Vel, but the analysis of the cluster or association to which WR stars are associated can give distances with a 50% to a 10% accuracy. The list of Galactic clusters, associations and clusters/association candidates has grown significantly in the last decade with the numerous deep, high resolution surveys of the Milky Way. In this work, we revisit the fundamental parameters of known clusters with WR stars, and we present the search for new ones. All our work is based on the catalogs from the VVV (from the VISTA telescope) and the UKIDS (from the UKIRT telescope) near infrared surveys. Finally, the relations between the fundamental parameters of clusters with WR stars are explored.
The Westerlund 1 (Wd1) cluster hosts a rich and varied collection of massive stars. Its dynamical youth and the absence of ongoing star formation indicate a coeval population. As such, the simultaneous presence of both late-type supergiants and Wolf-Rayet stars has defied explanation in the context of single-star evolution. Observational evidence points to a high binary fraction, hence this stellar population offers a robust test for stellar models accounting for both single-star and binary evolution. We present an optical to near-IR (VLT & NTT) spectroscopic analysis of 22 WR stars in Wd 1, delivering physical properties for the WR stars.
We discuss how these differ from the Galactic field population, and how they may be reconciled with the predictions of single and binary evolutionary models.
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.
I review our current understanding of the interaction between a Wolf-Rayet star's fast wind and the surrounding medium, and discuss to what extent the predictions of numerical simulations coincide with multiwavelength observations of Wolf-Rayet nebulae. Through a series of examples, I illustrate how changing the input physics affects the results of the numerical simulations. Finally, I discuss how numerical simulations together with multiwavelength observations of these objects allow us to unpick the previous mass-loss history of massive stars.
We analyse whether a stellar atmosphere model computed with the code CMFGEN provides an optimal description of the stellar observations of WR 136 and simultaneously reproduces the nebular observations of NGC 6888, such as the ionization degree, which is modelled with the pyCloudy code. All the observational material available (far and near UV and optical spectra) were used to constrain such models. We found that the stellar temperature T∗, at τ = 20, can be in a range between 70 000 and 110 000 K, but when using the nebula as an additional restriction, we found that the stellar models with T∗ ∼ 70 000 K represent the best solution for both, the star and the nebula.
The Galactic Center (GC) hosts three of the most massive WR rich, resolved young clusters in the Local Group as well as a large number of apparently isolated massive stars. Therefore, it constitutes a test bed to study the star formation history of the region, to probe a possible top-heavy scenario and to address massive star formation (clusters vs isolation) in such a dense and harsh environment. We present results from our ongoing infrared spectroscopic studies of WRs and other massive stars at the Center of the Milky Way.
A detailed and comprehensive study of the Wolf-Rayet stars of the nitrogen sequence (WN
stars) in the Small Magellanic Cloud (SMC) and the Large Magellanic Cloud (LMC) is presented.
We derived the fundamental stellar and wind parameters for more than 100 massive stars, encompassing almost the whole WN population in the Magellanic Clouds (MCs). The observations are fitted with synthetic spectra, using the PotsdamWolf-Rayet model atmosphere
code (PoWR). For this purpose, large grids of line-blanket models for different metallicities have been calculated, covering a wide range of stellar temperatures, mass-loss rates, and hydrogen abundances. Our comprehensive sample facilitates statistical studies of the WN properties in the MCs without selection bias. To investigate the impact of the low LMC metallicity and the even lower SMC metallicity, we compare our new results to previous analyses of the Galactic WN population and the late type WN stars from M31. Based on these studies we derived an empirical relation between the WN mass-loss rates and the metallicity. Current stellar evolution tracks, even when accounting for rotationally induced mixing, partly fail to reproduce the observed ranges of luminosities and initial masses.
The super massive binary system, η Car, experienced periastron passage in the summer of 2014. We observed the star twice around the maximum (forb =0.97, 2014 June 6) and just before the minimum (ϕorb =0.99, 2014 July 28) of its wind-wind colliding (WWC) X-ray emis-sion using the XMM-Newton and NuSTAR observatories, the latter of which is equipped with extremely hard X-ray (>10 keV) focusing mirrors. In both observations, NuSTAR detected the thermal X-ray tail up to 40-50 keV. The hard slope is consistent with an electron tem- perature of ∼6 keV, which is significantly higher than the ionization temperature (kT ∼4 keV) measured from the Fe K emission lines, assuming collisional equilibrium plasma. The spectrum did not show a hard power-law component above this energy range, unlike earlier detections with INTEGRAL and Suzaku. In the second NuSTAR observation, the X-ray flux above 5 keV declined gradually in ∼1 day. This result suggests that the WWC apex was gradually hidden behind the optically thick primary wind around conjunction.
Helium stars
(2015)
There are outstanding problems in trying to reproduce the observed nature of Wolf–Rayet stars from theoretical stellar models. We have investigated the effects of uncertainties, such as composition and mass-loss rate, on the evolution and structure of Wolf–Rayet stars and their lower mass brethren. We find that the normal Conti scenario needs to be altered, with different WR types being due to different initial masses as well as different stages of evolution.
Wolf-Rayet stars are very hot stars close to the Eddington limit. In the conditions encountered in their radiation pressure dominated outer layers several instabilities are expected to arise. These instabilities could influence both the dynamic of their optically thick winds and the observed spectral lines introducing small and large scale variability. We investigate the conditions in the convective envelopes of our helium star models and relate them to the appearance of a high number of stochastic density inhomogeneities, i.e. clumping in the optically thick wind. We also investigate the pulsational stability of these envelope, considering the effect of the high stellar wind mass loss rates.
We suggest several ideas which when combined could lead to a new mechanism for long-term pulsations of very hot and luminous stars. These involve the interplay between convection, radiation, atmospheric clumping and winds, which collectively feed back to stellar expansion and contraction. We discuss these ideas and point out the future work required in order to fill in the blanks.
We first give a short historical overview with some key facts of massive star population synthesis with binaries. We then discuss binary population codes and focus on two ingredients which are important for massive star population synthesis and which may be different in different codes. Population simulations with binaries is the third part where we consider the initial massive binary frequency, the RSG/WR and WC/WN and SNII/SNIbc number ratio's, the probable initial rotational velocity distribution of massive stars.
We investigate the rarity of the Wolf-Rayet X-ray binaries (WRXRBs) in contrast to their predecessors, the high mass X-ray binaries (HMXRBs). Recent studies suggest that common envelope (CE) mergers during the evolution of a HMXRBs may be responsible (Linden et al. 2012). We conduct a binary population synthesis to generate a population of HMXRBs mimicking the Galactic sample and vary the efficiency parameter during the CE phase to match the current WRXRB to HMXRB ratio. We find that ∼50% of systems must merge to match observational constraints.
Wolf-Rayet (WR) stars, as they are advanced stages of the life of massive stars, provide a good test for various physical processes involved in the modelling of massive stars, such as rotation and mass loss. In this paper, we show the outputs of the latest grids of single massive stars computed with the Geneva stellar evolution code, and compare them with some observations. We present a short discussion on the shortcomings of single stars models and we also briefly discuss the impact of binarity on the WR populations.
We compute spectral libraries for populations of coeval stars using state-of-the-art massive-star evolutionary tracks that account for different astrophysics including rotation and close-binarity. Our synthetic spectra account for stellar and nebular contributions. We use our models to obtain E(B – V ), age, and mass for six clusters in spiral galaxy NGC 1566, which have ages of < 50 Myr and masses of > 5 x 104M⊙ according to standard models. NGC 1566 was observed from the NUV to the I-band as part of the imaging Treasury HST program LEGUS: Legacy Extragalactic UV Survey. We aim to establish i) if the models provide reasonable fits to the data, ii) how well the models and photometry are able to constrain the cluster properties, and iii) how different the properties obtained with different models are.
The interaction between massive star formation and gas is a key ingredient in galaxy evolution. Given the level of observational detail currently achievable in nearby starbursts, they constitute ideal laboratories to study interaction process that contribute to global evolution in all types of galaxies. Wolf-Rayet (WR) stars, as an observational marker of high mass star formation, play a pivotal role and their winds can strongly influence the surrounding gas. Imaging spectroscopy of two nearby (<4 Mpc) starbursts, both of which show multiple regions with WR stars, are discussed. The relation between the WR content and the physical and chemical properties of the surrounding ionized gas is explored.
Using ESPaDOnS optical spectra of WR6, we search variations on the stellar wind parameters during the different phases of the spectral variations. We use the radiative transfer code CMFGEN (Hillier & Miller 1998) to determine the wind parameters. Our work gives mean parameters for WR6, Teff = 55 kK, M = 2.7 × 10^-5 M⊙/yr and v∞ =1700 km/s. Furthermore the line profiles variations at different phases are the consequence of a variation of mass loss rate and temperature un the winds. Effective temperature reaches 59 kK at the highest intensity, whereas the mass-loss rate decreases to 2.5 × 10^-5 M⊙/yr in that case. On the other hand, effective temperature decreases to 52.5 kK and the mass-loss rate increases to 3 × 10^-5 M/⊙yr when the line profile reach its minimum intensity. Results confirm the variable nature of the stellar wind, presented in this case on two of its fundamental parameters: temperature and mass-loss; which could be used to constrain the nature of the instability at the basis of the wind.
The X-ray observations of the colliding wind binary WR 21a is reported. The first monitoring performed by Swift/XRT in order to reveal the phase-locked variation. Our observations cover 201 different epochs from 2013 October 1 to 2015 January 30 for a total exposure of about 306 ks. It is found for the first time that the luminosity varies roughly in inverse proportion to the separation of the two stars before the X-ray maximum but later drops rapidly toward periastron.
Spots on stellar surfaces are thought to be stellar analogues of sunspots. Thus, starspots are direct manifestations of strong magnetic fields. Their decay rate is directly related to the magnetic diffusivity, which itself is a key quantity for the deduction of an activity cycle length. So far, no single starspot decay has been observed, and thus no stellar activity cycle was inferred from its corresponding turbulent diffusivity.
We investigate the evolution of starspots on the rapidly-rotating K0 giant XX Triangulum. Continuous high-resolution and phase-resolved spectroscopy was obtained with the robotic 1.2-m STELLA telescope on Tenerife over a timespan of six years. With our line-profile inversion code iMap we reconstruct a total of 36 consecutive Doppler maps. To quantify starspot area decay and growth, we match the observed images with simplified spot models based on a Monte-Carlo approach.
It is shown that the surface of XX Tri is covered with large high-latitude and even polar spots and with occasional small equatorial spots. Just over the course of six years, we see a systematically changing spot distribution with various time scales and morphology such as spot fragmentation and spot merging as well as spot decay and formation.
For the first time, a starspot decay rate on another star than the Sun is determined. From our spot-decay analysis we determine an average linear decay rate of D = -0.067±0.006 Gm^2/day. From this decay rate, we infer a turbulent diffusivity of η_τ = (6.3±0.5) x 10^14 cm^2/s and consequently predict an activity cycle of 26±6 years. The obtained cycle length matches very well with photometric observations.
Our time-series of Doppler maps further enables to investigate the differential rotation of XX Tri. We therefore applied a cross-correlation analysis. We detect a weak solar-like differential rotation with a surface shear of α = 0.016±0.003. This value agrees with similar studies of other RS CVn stars.
Furthermore, we found evidence for active longitudes and flip-flops. Whereas the more active longitude is located in phase towards the (unseen) companion star, the weaker active longitude is located at the opposite stellar hemisphere. From their periodic appearance, we infer a flip-flop cycle of ~2 years. Both activity phenomena are common on late-type binary stars.
Last but not least we redetermine several astrophysical properties of XX Tri and its binary system, as large datasets of photometric and spectroscopic observations are available since its last determination in 1999. Additionally, we compare the rotational spot-modulation from photometric and spectroscopic studies.
We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.
Two of the main physical parameters that govern the massive star evolution, the mass and the mass-loss rate, are still poorly determined from the observational point of view. Only binary systems could provide well constrained masses and colliding-wind binaries could bring some constraints on the mass-loss rate. Therefore, colliding-wind binaries turn out to be very promising objects. In this framework, we present detailed studies of basic observational data obtained with the XMM-Newton facility and combined with ground-based observations and other data. We expose the results for two particularly interesting WR+O colliding-wind binaries: WR22 and WR21a.
75 WR stars and 164 RSGs are identified in a single WFC3 pointing of our M101 survey. We find that within it's large star-forming complex NGC 5462 WR stars are preferentially located in the core whilst RSGs are found in the halo, suggesting two bursts of star-formation. A review of our WR candidates reveals that only ∼30% are detected in the archival broad-band ACS imaging whilst only ∼50% are associated with HII regions.
Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then the shock breakout signals diffuse in them and can extend their rise times significantly. We show that our inflated, hydrogen-free, WR stellar models with a radius of ∼R⊙ can have shock breakout signals longer than ∼ 60 s. The puzzlingly long shock breakout signal observed in the Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events might argue in favour of existence of inflated envelopes in general.
Spectroscopy is the preferred way to study the physical and wind properties of Wolf-Rayet (WR) stars, but with decreasing brightness and increasing distance of the object spectroscopy become very expensive. However, photometry still delivers a high signal to noise ratio. Current and past astronomical surveys and space missions provide large data sets, that can be harvested to discover new WR stars and study them over a wide metallicity range with the help of state of the art stellar atmosphere and evolutionary models.
An overview of the known Wolf-Rayet (WR) population of the Milky Way is presented, including a brief overview of historical catalogues and recent advances based on infrared photometric and spectroscopic observations resulting in the current census of 642 (vl.13 online catalogue). The observed distribution of WR stars is considered with respect to known star clusters, given that ≤20% of WR stars in the disk are located in clusters. WN stars outnumber WC stars at all galactocentric radii, while early-type WC stars are strongly biased against the inner Milky Way. Finally, recent estimates of the global WR population in the Milky Way are reassessed, with 1,200±100 estimated, such that the current census may be 50% complete. A characteristic WR lifetime of 0.25 Myr is inferred for an initial mass threshold of 25 M⊙.
Luminous Blue Variables (LBVs) are stars is a transitional phase massive stars may enter while evolving from main-sequence to Wolf-Rayet stars. The to LBVs intrinsic photometric variability is based on the modulation of the stellar spectrum. Within a few years the spectrum shifts from OB to AF type and back. During their cool phase LBVs are close to the Humphreys-Davidson (equivalent to Eddington/Omega-Gamma) limit. LBVs have a rather high mass loss rate, with stellar winds that are fast in the hot and slower in the cool phase of an LBV. These alternating wind velocities lead to the formation of LBV nebulae by wind-wind interactions. A nebula can also be formed in a spontaneous giant eruption in which larger amounts of mass are ejected. LBV nebulae are generally small (< 5 pc) mainly gaseous circumstellar nebulae, with a rather large fraction of LBV nebulae being bipolar. After the LBV phase the star will turn into a Wolf-Rayet star, but note that not all WR stars need to have passed the LBV phase. Some follow from the RSG and the most massive directly from the MS phase. In general WRs have a large mass loss and really fast stellar winds. The WR wind may interact with winds of earlier phases (MS, RSG) to form WR nebulae. As for WR with LBV progenitors the scenario might be different, here no older wind is present but an LBV nebula! The nature of WR nebulae are therefore manifold and in particular the connection (or family ties) of WR to LBV nebulae is important to understand the transition between these two phases, the evolution of massive stars, their winds, wind-wind and wind-nebula interactions. Looking at the similarities and differences of LBV and WR nebula, figuring what is a genuine LBV and WR nebula are the basic question addressed in the analysis presented here.
In this thesis we utilize resolved stellar populations to improve our understanding of galaxy formation and evolution. In the first part we improve a method for metallicity determination of faint old stellar systems, in the second and third part we analyze the individual history of six nearby disk galaxies outside the Local Group.
A New Calibration of the Color Metallicity Relation of Red Giants for HST data:
It is well known, that the color distribution of stars on the the Red Giant Branch (RGB) can be used to determine metallicities of old stellar populations that have only shallow photometry. Based on the largest sample of globular clusters ever used for such studies, we quantify the relation between metallicity and color in the widely used HST ACS filters F606W and F814W.
We use a sample of globular clusters from the ACS Globular Cluster Survey and measure their RGB color at given absolute magnitudes to derive the color-metallicity relation. We find a clear relation between metallicity and RGB color; we investigate the scatter and the uncertainties in this relation and show its limitations. A comparison with isochrones shows reasonably good agreement with BaSTI models, a small offset to Dartmouth models, and a larger offset to Padua models.
Even for the best globular cluster data available, the metallicity of a simple stellar population can be determined from the RGB alone only with an accuracy of 0.3 dex for [M/H]<-1, and 0.15 dex for [M/H]>-1. For mixed populations, as they are observed in external galaxies, the uncertainties will be even larger due to uncertainties in extinction, age, etc. Therefore caution is necessary when interpreting photometric metallicities.
The Structural History of Nearby Low Mass Disk Galaxies:
We study the individual evolution histories of three nearby, low-mass, edge-on galaxies (IC5052, NGC4244, NGC5023).
Using the color magnitude diagrams of resolved stellar populations, we construct star count density maps for populations of different ages and analyze the change of structural parameters with stellar age within each galaxy.
The three galaxies show low vertical heating rates, which are much lower than the heating rate of the Milky Way. This indicates that heating agents, as giant molecular clouds and spiral structure are weak in low mass galaxies.
We do not detect a separate thick disk in any of the three galaxies, even though our observations cover a larger range in equivalent surface brightness than any integrated light study. While scaleheights increase with age, each population can be well described by a single disk. Only two of the galaxies contain a very weak additional component, which we identify as the faint halo. The mass of these faint halos is less than 1% of the mass of the disk.
All populations in the three galaxies exhibit no or only little flaring. While this finding is consistent with previous integrated light studies, it poses strong constraints on galaxy formation models, because most theoretical simulations often find strong flaring due to interactions or radial migration.
Furthermore, we find breaks in the radial profiles of all three galaxies. The radii of these breaks are independent of age, and the break strength is decreasing with age in two of the galaxies (NGC4244 and NGC5023). This is consistent with break formation models, that combine a star formation cutoff with radial migration. The differing behavior of IC5052 can be explained by a recent interaction or minor merger.
The Structural History of Massive Disk Galaxies:
We extend the structural analysis of stellar populations with distinct ages to three massive galaxies, NGC891, NGC4565 and NGC7814. While confusion effects due to the high stellar number densities in their central region, and the prominent dust lanes inhibit an detailed analysis of the radial profiles, we can study their vertical structure.
These massive galaxies also have a slower heating than the Milky Way, comparable to the low mass galaxies. This can be traced back to their already thick young populations and thick layers of their interstellar medium.
We do not find a clear separate thick disk in any of these three galaxies; all populations can be described by a single disk plus a S\'ersic bulge/halo component. In contrast to the low mass galaxies, we cannot rule out the presence of thick disks in the massive galaxies, because of the strong influence of the halo, that might hide the possible contribution of the thick disk to the vertical star count profiles. However, the faintness of the possible thick disks still points to problems in the earlier ubiquitous findings of thick disks in external galaxies.
We study the interaction of line-driven winds from massive stars with the magnetic field rooted in these stars by carrying out numerical simulations using the Nirvana MHD code in 2D in spherical polar coordinates. The code's adaptive mesh refinement feature allows high spatial resolution across the whole simulation box. We study both O and Wolf-Rayet stars for a range of magnetic field strengths from weak to strong as measured by the confinement parameter. For weak fields our simulations show that the initially dipolar field opens up far away from the star and a thin disk-like structure forms in the equatorial plane of the magnetic field. For stronger fields the disk is disrupted close to the stellar surface and closed field lines persist at low latitudes. For very strong fields a pronounced magnetosphere forms where the gas is forced to move along the field lines and eventually falls back to the stellar surface.
Magnetic fields, non-thermal radiation and particle acceleration in colliding winds of WR-O stars
(2015)
Non-thermal emission has been detected in WR-stars for many years at long wavelengths spectral range, in general attributed to synchrotron emission. Two key ingredients are needed to explain such emissions, namely magnetic fields and relativistic particles. Particles can be accelerated to relativistic speeds by Fermi processes at strong shocks. Therefore, strong synchrotron emission is usually attributed to WR binarity. The magnetic field may also be amplified at shocks, however the actual picture of the magnetic field geometry, intensity, and its role on the acceleration of particles at WR binary systems is still unclear. In this work we discuss the recent developments in MHD modelling of wind-wind collision regions by means of numerical simulations, and the coupled particle acceleration processes related.
Using a code that employs a self-consistent method for computing the effects of photoionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive Wolf-Rayet (W-R) stars. Our algorithm incorporates a simplified model of the photo-ionization source, computes the fractional ionization of hydrogen due to the photoionizing flux and recombination, and determines self-consistently the energy balance due to ionization, photo-heating and radiative cooling. We take into account changes in stellar properties and mass-loss over the star's evolution. Our multi-dimensional simulations clearly reveal the presence of strong ionization front instabilities. Using various X-ray emission models, and abundances consistent with those derived for W-R nebulae, we compute the X-ray flux and spectra from our wind bubble models. We show the evolution of the X-ray spectral features with time over the evolution of the star, taking the absorption of the X-rays by the ionized bubble into account. Our simulated X-ray spectra compare reasonably well with observed spectra of Wolf-Rayet bubbles. They suggest that X-ray nebulae around massive stars may not be easily detectable, consistent with observations.∗
The main goal of this cumulative thesis is the derivation of surface emissivity data in the infrared from radiance measurements of Venus. Since these data are diagnostic of the chemical composition and grain size of the surface material, they can help to improve knowledge of the planet’s geology. Spectrally resolved images of nightside emissions in the range 1.0-5.1 μm were recently acquired by the InfraRed Mapping channel of the Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS-M-IR) aboard ESA’s Venus EXpress (VEX). Surface and deep atmospheric thermal emissions in this spectral range are strongly obscured by the extremely opaque atmosphere, but three narrow spectral windows at 1.02, 1.10, and 1.18 μm allow the sounding of the surface. Additional windows between 1.3 and 2.6 μm provide information on atmospheric parameters that is required to interpret the surface signals. Quantitative data on surface and atmosphere can be retrieved from the measured spectra by comparing them to simulated spectra. A numerical radiative transfer model is used in this work to simulate the observable radiation as a function of atmospheric, surface, and instrumental parameters. It is a line-by-line model taking into account thermal emissions by surface and atmosphere as well as absorption and multiple scattering by gases and clouds. The VIRTIS-M-IR measurements are first preprocessed to obtain an optimal data basis for the subsequent steps. In this process, a detailed detector responsivity analysis enables the optimization of the data consistency. The measurement data have a relatively low spectral information content, and different parameter vectors can describe the same measured spectrum equally well. A usual method to regularize the retrieval of the wanted parameters from a measured spectrum is to take into account a priori mean values and standard deviations of the parameters to be retrieved. This decreases the probability to obtain unreasonable parameter values. The multi-spectrum retrieval algorithm MSR is developed to additionally consider physically realistic spatial and temporal a priori correlations between retrieval parameters describing different measurements. Neglecting geologic activity, MSR also allows the retrieval of an emissivity map as a parameter vector that is common to several spectrally resolved images that cover the same surface target. Even applying MSR, it is difficult to obtain reliable emissivity maps in absolute values. A detailed retrieval error analysis based on synthetic spectra reveals that this is mainly due to interferences from parameters that cannot be derived from the spectra themselves, but that have to be set to assumed values to enable the radiative transfer simulations. The MSR retrieval of emissivity maps relative to a fixed emissivity is shown to effectively avoid most emissivity retrieval errors. Relative emissivity maps at 1.02, 1.10, and 1.18 μm are finally derived from many VIRTIS-M-IR measurements that cover a surface target at Themis Regio. They are interpreted as spatial variations relative to an assumed emissivity mean of the target. It is verified that the maps are largely independent of the choice of many interfering parameters as well as the utilized measurement data set. These are the first Venus IR emissivity data maps based on a consistent application of a full radiative transfer simulation and a retrieval algorithm that respects a priori information. The maps are sufficiently reliable for future geologic interpretations.
The subject of the present thesis is the one-dimensional Bose gas. Since long-rang order is destroyed by infra-red fluctuations in one dimension, only the formation of a quasi-condensate is possible, which exhibits suppressed density fluctuations, but whose phase fluctuates strongly. It is shown that modified mean-field theories based on a symmetry-breaking approach can even characterise phase coherence properties of such a quasi-condensate properly. A correct description of the transition from the degenerate ideal Bose gas to the quasi-condensate, which is a smooth cross-over rather than a phase transition, is not possible though. Basic conditions for the applicability of the theories are not fulfilled in this regime, such that the existence of a critical point is predicted.
The theories are compared on the basis of their excitation sprectum, equation of state, density fluctuations and related correlation functions. High-temperature expansions of the corresponding integrals are derived analytically for the numerical evaluation of the self-consistent integral equations. Apart from that, the Stochastic Gross-Pitaevskii equation (SGPE), a non-linear Langevin equation, is analysed numerically by means of Monte-Carlo simulations and the results are compared to those of the mean-field theories. In this context, a lot of attention is payed to the appropriate choice of the parameters. The simulations prove that the SGPE is capable of describing the cross-over properly, but highlight the limitations of the widely used local density approximation as well.
Although we all use the name Wolf-Rayet to refer to specific groups of stars, “Wolf-Rayet” per se is really an astrophysical phenomenon of fast-moving, hot plasma, normally expanding around a hot star. However, expediency demands that we follow established traditions by referring to three specific kinds of WR stars: (1) cWR, “classical” He-burning descendants of massive, O-type stars, presumably all of which pass through a WR stage; (2) WNh, the most massive and luminous hydrogen-rich main-sequence stars with strong winds; and (3) [WR], the central stars of some 15 % of Planetary Nebulae. Wolf-Rayet stars are the epitome of relatively stable stars with the highest mass-loss rates for their kind. It behooves us to understand the what, how and why of this circumstance, along with its manyfold and fascinating consequences.
Biological materials, in addition to having remarkable physical properties, can also change shape and volume. These shape and volume changes allow organisms to form new tissue during growth and morphogenesis, as well as to repair and remodel old tissues. In addition shape or volume changes in an existing tissue can lead to useful motion or force generation (actuation) that may even still function in the dead organism, such as in the well known example of the hygroscopic opening or closing behaviour of the pine cone. Both growth and actuation of tissues are mediated, in addition to biochemical factors, by the physical constraints of the surrounding environment and the architecture of the underlying tissue. This habilitation thesis describes biophysical studies carried out over the past years on growth and swelling mediated shape changes in biological systems. These studies use a combination of theoretical and experimental tools to attempt to elucidate the physical mechanisms governing geometry controlled tissue growth and geometry constrained tissue swelling. It is hoped that in addition to helping understand fundamental processes of growth and morphogenesis, ideas stemming from such studies can also be used to design new materials for medicine and robotics.
The non-linear behaviour of the atmospheric dynamics is not well understood and makes the evaluation and usage of regional climate models (RCMs) difficult. Due to these non-linearities, chaos and internal variability (IV) within the RCMs are induced, leading to a sensitivity of RCMs to their initial conditions (IC). The IV is the ability of RCMs to realise different solutions of simulations that differ in their IC, but have the same lower and lateral boundary conditions (LBC), hence can be defined as the across-member spread between the ensemble members.
For the investigation of the IV and the dynamical and diabatic contributions generating the IV four ensembles of RCM simulations are performed with the atmospheric regional model HIRHAM5. The integration area is the Arctic and each ensemble consists of 20 members. The ensembles cover the time period from July to September for the years 2006, 2007, 2009 and 2012. The ensemble members have the same LBC and differ in their IC only. The different IC are arranged by an initialisation time that shifts successively by six hours. Within each ensemble the first simulation starts on 1st July at 00 UTC and the last simulation starts on 5th July at 18 UTC and each simulation runs until 30th September. The analysed time period ranges from 6th July to 30th September, the time period that is covered by all ensemble members. The model runs without any nudging to allow a free development of each simulation to get the full internal variability within the HIRHAM5.
As a measure of the model generated IV, the across-member standard deviation and the across-member variance is used and the dynamical and diabatic processes influencing the IV are estimated by applying a diagnostic budget study for the IV tendency of the potential temperature developed by Nikiema and Laprise [2010] and Nikiema and Laprise [2011]. The diagnostic budget study is based on the first law of thermodynamics for potential temperature and the mass-continuity equation. The resulting budget equation reveals seven contributions to the potential temperature IV tendency.
As a first study, this work analyses the IV within the HIRHAM5. Therefore, atmospheric circulation parameters and the potential temperature for all four ensemble years are investigated. Similar to previous studies, the IV fluctuates strongly in time. Further, due to the fact that all ensemble members are forced with the same LBC, the IV depends on the vertical level within the troposphere, with high values in the lower troposphere and at 500 hPa and low values in the upper troposphere and at the surface. By the same reason, the spatial distribution shows low values of IV at the boundaries of the model domain.
The diagnostic budget study for the IV tendency of potential temperature reveals that the seven contributions fluctuate in time like the IV. However, the individual terms reach different absolute magnitudes. The budget study identifies the horizontal and vertical ‘baroclinic’ terms as the main contributors to the IV tendency, with the horizontal ‘baroclinic’ term producing and the vertical ‘baroclinic’ term reducing the IV. The other terms fluctuate around zero, because they are small in general or are balanced due to the domain average.
The comparison of the results obtained for the four different ensembles (summers 2006, 2007, 2009 and 2012) reveals that on average the findings for each ensemble are quite similar concerning the magnitude and the general pattern of IV and its contributions. However, near the surface a weaker IV is produced with decreasing sea ice extent. This is caused by a smaller impact of the horizontal 'baroclinic' term over some regions and by the changing diabatic processes, particularly a more intense reducing tendency of the IV due to condensative heating. However, it has to be emphasised that the behaviour of the IV and its dynamical and diabatic contributions are influenced mainly by complex atmospheric feedbacks and large-scale processes and not by the sea ice distribution.
Additionally, a comparison with a second RCM covering the Arctic and using the same LBCs and IC is performed. For both models very similar results concerning the IV and its dynamical and diabatic contributions are found. Hence, this investigation leads to the conclusion that the IV is a natural phenomenon and is independent from the applied RCM.
We summarize past and current surveys for WRs among the Local Group galaxies, empha- sizing both the why and how. Such studies are invaluable for helping us learn about massive star evolution, and for providing sensitive tests of the stellar evolution models. But for such surveys to be useful, the completeness limits must be well understood. We illustrate that point by following the “evolution” of the observed WC/WN ratio in nearby galaxies. We end by examining our new survey for WR stars in the Magellanic Clouds, which has revealed a new type of WN star, never before seen.
The emission-line dominated spectra of Wolf-Rayet stars are formed in expanding layers of their atmosphere, i.e. in their strong stellar wind. Adequate modeling of such spectra has to face a couple of difficulties. Because of the supersonic motion, the radiative transfer is preferably formulated in the co-moving frame. The strong deviations from local thermodynamical equilibrium (LTE) require to solve the equations of statistical equilibrium for the population numbers, accounting for many hundred atomic energy levels and thousands of line transitions. Moreover, millions of lines from iron-group elements must be taken into account for their blanketing effect. Model atmospheres of the described kind can reproduce the observed WR spectra satisfyingly, and have been widely applied for corresponding spectral analyses.
The total population of Wolf-Rayet (WR) stars in the Galaxy is predicted by models to be as many as ~6000 stars, and yet the number of catalogued WR stars as a result of optical surveys was far lower than this (~200) at the turn of this century. When beginning our WR searches using infrared techniques it was not clear whether WR number predictions were too optimistic or whether there was more hidden behind interstellar and circumstellar extinction. During the last decade we pioneered a technique of exploiting the near- and mid-infrared continuum colours for individual point sources provided by large-format surveys of the Galaxy, including 2MASS and Spitzer/GLIMPSE, to pierce through the dust and reveal newly discovered WR stars throughout the Galactic Plane. The key item to the colour discrimination is via the characteristic infrared spectral index produced by the strong winds of the WR stars, combined with dust extinction, which place WR stars in a relatively depopulated area of infrared colour-colour diagrams. The use of the Spitzer/GLIMPSE 8µm and, more recently, WISE 22µm fluxes together with cross-referencing with X-ray measurements in selected Galactic regions have enabled improved candidate lists that increased our confirmation success rate, achieved via follow-up infrared and optical spectroscopy. To date a total of 102 new WR stars have been found with many more
candidates still available for follow-up. This constitutes an addition of ~16% to the current
inventory of 642 Galactic WR stars. In this talk we review our methods and provide some
new results and a preliminary review of their stellar and interstellar medium environments. We provide a roadmap for the future of this search, including statistical modeling, and what we can add to star formation and high mass star evolution studies.
In this contribution we present some preliminary results obtained from a SOAR-Goodman optical spectroscopic survey aimed to confirm the OIf* - OIf*/WN nature of a sample of Galactic candidates that were previously confirmed as massive stars based on near-infrared spectra taken with OSIRIS at SOAR. With only a few of such stars known in the Galaxy to date, our study significantly contributes to improve the number of known Galactic O2If* stars, as well as almost doubling the number of known members of the galactic sample of the rare type OIf*/WN.
The enigmatic oxygen-sequence Wolf-Rayet stars represent a rare stage in the evolution of massive stars. Their properties can provide unique constraints on the pre-supernova evolution of massive stars. This work presents the results of a quantitative spectroscopic analysis of the known single WO stars, with the aim to obtain the key stellar parameters and deduce their evolutionary state.X-Shooter spectra of the WO stars are modeled using the line-blanketed non-local thermal equilibrium atmosphere code cmfgen. The obtained stellar parameters show that the WO stars are very hot, with temperatures ranging from 150 kK to 210 kK. Their chemical composition is dominated by carbon (>50%), while the helium mass fraction is very low (down to 14%). Oxygen mass fractions reach as high as 25%. These properties can be reproduced with dedicated evolutionary models for helium stars, which show that the stars are post core-helium burning and very close to their eventual supernova explosion. The helium-star masses indicate initial masses or approximately 40 - 60M⊙.Thus, WO stars represent the final evolutionary stage of stars with estimated initial masses of 40 - 60M⊙. They are post core-helium burning and may explode as type Ic supernovae within a few thousand years.
Professional and amateur astronomers around the world contributed to a 4-month long campaign in 2013, mainly in spectroscopy but also in photometry, interferometry and polarimetry, to observe the first 3 Wolf-Rayet stars discovered: WR 134 (WN6b), WR 135 (WC8) and WR 137 (WC7pd+O9). Each of these stars are interesting in their own way, showing a variety of stellar wind structures. The spectroscopic data from this campaign were reduced and analyzed for WR 134 in order to better understand its behavior and long-term periodicity in the context of CIRs in the wind. We will be presenting the results of these spectroscopic data, which include the confirmation of the CIR variability and a time-coherency of ∼ 40 days (half-life of ∼ 20 days).
As part of our ongoing Wolf-Rayet (WR) Magellanic Cloud survey, we have discovered 13 new WRs. However, the most exciting outcome of our survey is not the number of new WRs, but their unique characteristics. Eight of our discoveries appear to belong to an entirely new class of WRs. While one might naively classify these stars as WN3+O3V binaries, such a pairing is unlikely. Preliminary CMFGN modeling suggests physical parameters similar to early-type WNs in the Large Magellanic Cloud except with mass-loss rates three to five times lower and slightly higher temperatures. The evolution status of these stars remains an open question.
Obtaining a complete census of massive, evolved stars in a galaxy would be a key ingredient for testing stellar evolution models. However, as the evolution of stars is also strongly dependent on their metallicity, it is inevitable to have this kind of data for a variety of galaxies with different metallicities. Between 2009 and 2011, we conducted the Magellanic Clouds Massive Stars and Feedback Survey (MSCF); a spatially complete, multi-epoch, broad- and narrow-band optical imaging survey of the Large and Small Magellanic Clouds. With the inclusion of shallow images, we are able to give a complete photometric catalog of stars between B ≈ 18 and B ≈ 19 mag.
These observations were augmented with additional photometric data of similar spatial res-
olution from UV to IR (e.g. from GALEX, 2MASS and Spitzer) in order to sample a large portion of the spectral energy distribution of the brightest stars (B < 16 mag) in the Magel- lanic Clouds. Using these data, were are able to train a machine learning algorithm that gives us a good estimate of the spectral type of tens of thousands of stars.
This method can be applied to the search for Wolf-Rayet-Stars to obtain a sample of candi- dates for follow-up observations. As this approach can, in principle, be adopted for any resolved galaxy as long as sufficient photometric data is available, it can form an effective alternative method to the classical strategies (e.g. He II filter imaging).
Colliding Wolf-Rayet (WR) winds produce thermal X-ray emission widely observed by X-ray telescopes. In wide WR+O binaries, such as WR 140, the X-ray flux is tied to the orbital phase, and is a direct probe of the winds’ properties. In the Galactic center, ~30 WRs orbit the super massive black hole (SMBH) within ~10”, leading to a smorgasbord of wind-wind collisions. To model the X-ray emission of WR 140 and the Galactic center, we perform 3D hydrodynamic simulations to trace the complex gaseous flows, and then carry out 3D radiative transfer calculations to compute the variable X-ray spectra. The model WR 140 RXTE light curve matches the data well for all phases except the X-ray minimum associated with periastron, while the model spectra agree with the RXTE hardness ratio and the shape of the Suzaku observations throughout the orbit. The Galactic center model of the Chandra flux and spectral shape match well in the region r ≤ 3”, but the model flux falls off too rapidly beyond this radius.
The feedback from massive stars is important to super star cluster (SSC) evolution and the timescales on which it occurs. SSCs form embedded in thick material, and eventually, the cluster is cleared out and revealed at optical wavelengths – however, this transition is not well understood. We are investigating this critical SSC evolutionary transition with a multi-wavelength observational campaign. Although previously thought to appear after the cluster has fully removed embedding natal material, we have found that SSCs may host large populations of Wolf-Rayet stars. These evolved stars provide ionization and mechanical feedback that we hypothesize is the tipping point in the combined feedback processes that drive a SSC to emerge. Utilizing optical spectra obtained with the 4m Mayall Telescope at Kitt Peak National Observatory and the 6.5m MMT, we have compiled a sample of embedded SSCs that are likely undergoing this short-lived evolutionary phase and in which we confirm the presence of Wolf-Rayet stars. Early results suggest that WRs may accelerate the cluster emergence.
Ring Nebulae
(2015)
Preliminary results are presented from spectroscopic data in the optical range of the Galactic ring nebulae NGC 6888, G2:4+1:4, RCW 58 and Sh2-308. Deep observations with long exposure times were carried out at the 6.5m Clay Telescope and at the 10.4m Gran Telescopio Canarias. In NGC 6888, recombination lines of C ii, O ii and N ii are detected with signal-to-noise ratios higher than 8. The chemical content of NGC 6888 is discussed within the chemical enrichment predicted by evolution models of massive stars. For all nebulae, a forthcoming work will content in-depth details about observations, analysis and final results (Esteban et al. 2015, in prep.).
Concluding Remarks
(2015)
The gas cloud G2 is currently being tidally disrupted by the Galactic Centre super-massive black hole, Sgr A*. The region around the black hole is populated by ∼ 30 Wolf-Rayet stars, which produce strong outflows. Here we explore the possibility that gas clumps like G2 originate from the collision of stellar winds via the non-linear thin shell instability.
We present the results of the new photometric observations of the famous hypergiant PCygni. New observations were obtained in 2014 using the 48 cm Cassegrain telescope of the Abastumani Astrophysical Observatory, Georgia. We reveal some interesting behaviors of the B,V,R,I light curves, and also report new results on the periodicity of PCygni's variation. The latter result is based on the analysis of the photometric data (U,B,V filters) collected at the Abastumani Observatory between 1937 and 1983.
We present results of investigation of spectral variability of one of the most interesting massive stars, Romano's star (M33/V532 or GR290), located in the M33 galaxy. Brightness of the star changes together with its spectral class, which varies from WN11 to WN8. Using CMFGEN code we estimated parameters of stellar atmosphere and found that during last ten years bolometric luminosity of the star changed synchronously with stellar magnitude. Our calculations argue in favor of the hypothesis of a post-LBV status of GR290.
In this review, I discuss the suitability of massive star progenitors, evolved in isolation or in interacting binaries, for the production of observed supernovae (SNe) IIb, Ib, Ic. These SN types can be explained through variations in composition. The critical need of non-thermal effects to produce He I lines favours low-mass He-rich ejecta (in which ^56 Ni can be more easily mixed with He) for the production of SNe IIb/Ib, which thus may arise preferentially from moderate-mass donors in interacting binaries. SNe Ic may instead arise from higher mass progenitors, He-poor or not, because their larger CO cores prevent efficient non-thermal excitation of He i lines. However, current single star evolution models tend to produce Wolf-Rayet (WR) stars at death that have a final mass of > 10 M⊙. Single WR star explosion models produce ejecta that are too massive to match the observed light curve widths and rise times of SNe IIb/Ib/Ic, unless their kinetic energy is systematically and far greater than the canonical value of 10^56 erg. Future work is needed to evaluate the energy/mass degeneracy in light curve properties. Alternatively, a greater mass loss during the WR phase, perhaps in the form of eruptions, as evidenced in SNe Ibn, may reduce the final WR mass. If viable, such explosions would nonetheless favour a SN Ic, not a Ib.
Nearly 50 post-common-envelope (post-CE) close binary central stars of planetary nebulae (CSPNe) are now known. Most contain either main sequence or white dwarf (WD) companions that orbit the WD primary in around 0.1–1.0 days. Only PN G222.8–04.2 and NGC 5189 have post-CE CSPNe with a Wolf-Rayet star primary (denoted [WR]), the low-mass analogues of massive Wolf-Rayet stars. It is not well understood how H-deficient [WR] CSPNe form, even though they are relatively common, appearing in over 100 PNe. The discovery and characterisation of post-CE [WR] CSPNe is essential to determine whether proposed binary formation scenarios are feasible to explain this enigmatic class of stars. The existence of post-CE [WR] binaries alone suggests binary mergers are not necessarily a pathway to form [WR] stars. Here we give an overview of the initial results of a radial velocity monitoring programme of [WR] CSPNe to search for new binaries. We discuss the motivation for the survey and the associated strong selection effects. The mass functions determined for PN G222.8–04.2 and NGC 5189, together with literature photometric variability data of other [WR] CSPNe, suggest that of the post-CE [WR] CSPNe yet to be found, most will have WD or subdwarf O/B-type companions in wider orbits than typical post-CE CSPNe (several days or months c.f. less than a day).
Carbon-rich Wolf-Rayet stars are efficient carbon dust makers. Despite the strong evidence for dust formation in these objects provided by infrared thermal emission from dust, the routes to nucleation and condensation and the physical conditions required for dust production are still poorly understood. We discuss here the potential routes to carbon dust and the possible locations conducive to dust formation in the colliding winds of WC binaries.
Observations of the WC9+OB system WR65 in the infrared show variations of its dust emission consistent with a period near 4.8 yr, suggesting formation in a colliding-wind binary (CWB) having an elliptical orbit. If we adopt the IR maximum as zero phase, the times of X-ray maximum count and minimum extinction to the hard component measured by Oskinova & Hamann fall at phases 0.4–0.5, when the separation of the WC9 and OB stars is greatest. We consider WR65 in the context of other WC8–9+OB stars showing dust emission.
WR Time Series Photometry
(2015)
We take a comprehensive look at Wolf Rayet photometric variability using the MOST satellite. This sample, consisting of 6 WR stars and 6 WC stars defies all typical photometric analysis. We do, however, confirm the presence of unusual periodic signals resembling sawtooth waves which are present in 11 out of 12 stars in this sample.
The stellar Eddington limit
(2015)
It is often assumed that when stars reach their Eddington limit, strong outflows are initiated, and that this happens only for extreme stellar masses. We discuss here that in models of up to 500 M⊙, the Eddington limit is never reached at the stellar surface. Instead, we argue that the Eddington limit is reached inside the stellar envelope in hydrogen-rich stars above ∼ 30 M⊙ and in Wolf-Rayet stars above ∼ 7 M⊙, with drastic effects for their struture and stability.
We look at how the dynamics of colliding wind binaries (CWB) can be investigated in 2D, and how several parameters influence the dynamics of the small scale structures inside the colliding wind and the shocked regions, as well as in how the dynamics influence the shape of the collision region at large distances. The parameters we adopt are based on the binary system WR98a, one of the few Wolf-Rayet (WR) dusty pinwheels known.
The morphological appearance of massive stars across their post-Main Sequence evolution and before the SN event is very uncertain, both from a theoretical and observational perspective. We recently developed coupled stellar evolution and atmospheric modeling of stars done with the Geneva and CMFGEN codes, for initial masses between 9 and 120 M⊙. We are able to predict the observables such as the high-resolution spectrum and broadband photometry. Here I discuss how the spectrum of a massive star changes across its evolution and before death, with focus on the WR stage. Our models indicate that single stars with initial masses larger than 30 M⊙ end their lives as WR stars. Depending on rotation, the spectrum of the star can either be that of a WN or WO subtype at the pre-SN stage. Our models allow, for the first time, direct comparison between predictions from stellar evolution models and observations of SN progenitors.
Eta Carinae
(2015)
Since Augusto Damineli's demonstration in 1996 that Eta Carinae is a binary with a 5.52 year period, many innovative observations and increasingly advanced three-dimensional models have led to considerable insight on this massive system that ejected at least ten, possibly forty, solar masses in the nineteenth century. Here we present a review of our current understanding of this complex system and point out continuing puzzles.
Key physical ingredients governing the evolution of massive stars are mass losses, convection and mixing in radiative zones. These effects are important both in the frame of single and close binary evolution. The present paper addresses two points: 1) the differences between two families of rotating models, i.e. the family of models computed with and without an efficient transport of angular momentum in radiative zones; 2) The impact of the mass losses in single and in close binary models.
I address uncertainties on the spatial distribution and mass of the dust formed in η Carinae's Homunculus nebula with data being combined from several space- and ground-based facilities spanning near-infrared to sub-mm wavelengths, in terms of observational constraints and modeling. Until these aspects are better understood, the mass loss history and mechanisms responsible for η Car's enormous eruption(s) remain poorly constrained.
The Wolf-Rayet (WR) phenomenon is widespread in astronomy. It involves classical WRs, very massive stars (VMS), WR central stars of planetary nebula CSPN [WRs], and supernovae (SNe). But what is the root cause for a certain type of object to turn into an emission-line star? In this contribution, I discuss the basic aspects of radiation-driven winds that might reveal the ultimate difference between WR stars and canonical O-type stars. I discuss the aspects of (i) self-enrichment via CNO elements, (ii) high effective temperatures (Tₑff), (iii) an increase in the helium abundance (Y ), and finally (iv) the Eddington factor Γₑ. Over the last couple of years, we have made a breakthrough in our understanding of Γₑ -dependent mass loss, which will have far-reaching consequences for the evolution and fate of the most massive stars in the Universe. Finally, I discuss the prospects for studies of the WR phenomenon in the highest redshift Lyα and He ii emitting galaxies.
We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system η Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on η Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty `pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.
We present results from our near-infrared spectroscopy with VLT/ISAAC of four, massive eclipsing binary systems in the young, heavily reddened, massive Danks clusters. We derive accurate fundamental parameters and the distance to these massive systems, which comprise of OIf+, WR and O-type stars. Our goal is to increase the sample of well-studied WR stars and constrain their physics by comparison with evolutionary models.
Wolf-Rayet stars are important sources for the enrichment of the ISM with nuclear processed elements, UV photons and momentum. They are descendants of high-mass stars for which short lifetimes and transition times can hamper the spectral classification of the stars in their different evolutionary phases. The expanded stellar atmospheres of Wolf-Rayet stars can show spectra which seem inconsistent with the anticipated underlying evolution phase, for example in late hydrogen-burning WN stars and Of/WN transition stars. We present a sequence of synthetic spectra of the Potsdam Wolf-Rayet models based on the latest Geneva stellar evolution models. This will visualize the changes in stellar spectra over a full stellar lifetime. Direct comparison with observed stellar spectra, as well as the evolution of diagnostic line ratios will improve the connection of spectral classification and evolution phase.
HD5980
(2015)
HD5980 is a multiple system containing at least 3 very massive and luminous stars. Located in the Small Magellanic Cloud, it is an ideal system for studying the massive star structure and evolutionary processes in low-metallicity environments. Intensely observed over the past few decades, HD5980 is a treasure trove of information on stellar wind structure, on wind-wind collisions and on the formation of wind-blown circumstellar structures. In addition, its characteristics suggest that the eclipsing WR+LBV stars of the system are the product of quasihomogeneous chemical evolution, thus making them candidate pair production supernovae or GRB progenitors. This paper summarizes some of the outstanding results derived from half a century of observations and recent theoretical studies.
We present the first physical characterization of the young open cluster VVVCL041. We spectroscopically observed the cluster main-sequence stellar population and a very-massive star candidate: WR62-2. CMFGEN modelling to our near-infrared spectra indicates that WR62-2 is a very luminous (10^6.4±0.2 L⊙)and massive (∼ 80M⊙) star.
We an optically-thick, transonic, steady wind model for a H-free Wolf-Rayet star. A bifurcation is found across a critical mass loss rate Mb. Slower winds M < Mb extend by several hydrostatic stellar radii, reproduce features of envelope in ation from Petrovic et al. (2006) and Gräfener et al. (2012), and are energetically unbound. This work is of particular interest for extended envelopes and winds, radiative hydrodynamic instabilities (eg. wind stagnation, clumping, etc.), and NLTE atmospheric models.
In this thesis, I study ultrafast dynamics in perovskite oxides using time resolved broadband spectroscopy. I focus on the observation of coherent phonon propagation by time resolved Brillouin scattering: following the excition of metal transducer films with a femtosecond infrared pump pulse, coherent phonon dynamics in the GHz frequency range are triggered. Their propagation is monitored using a delayed white light probe pulse. The technique is illustrated on various thin films and multilayered samples. I apply the technique to investigate the linear and nonlinear acoustic response in bulk SrTiO_3, which displays a ferroelastic phase transition from a cubic to a tetragonal structural phase at T_a=105 K. In the linear regime, I observe a coupling of the observed acoustic phonon mode to the softening optic modes describing the phase transition. In the nonlinear regime, I find a giant slowing down of the sound velocity in the low temperature phase that is only observable for a strain amplitude exceeding the tetragonality of the material. It is attributed to a coupling of the high frequency phonons to ferroelastic domain walls in the material. I propose a new mechanism for the coupling of strain waves to the domain walls that is only effective for high amplitude strain. A detailed study of the phonon attenuation across a wide temperature range shows that the phonon attenuation at low temperatures is influenced by the domain configuration, which is determined by interface strain. Preliminary measurements on magnetic-ferroelectric multilayers reveal that the excitation fluence needs to be carefully controlled when dynamics at phase transitions are studied.
We found original observations of PCygni by E. Kharadze and N. Magalashvili in the archives of the Abastumani Observatory. These observations were carried out in the period 1951–1983. Initially they used 29 Cygni as a comparison star, and all observations of PCygni were processed using this star. On the basis of their calculations, the authors decided that PCygni may be a WUMa type binary with an orbital period of 0.500565 d, but this hypothesis was not confirmed. The only observations that have been published in the Bulletin of the Abastumani Astrophysical Observatory were those of of 1951–1955. There are whole sets of observational data not only for PCygni and 29 Cygni, but in the majority of cases also for 36 Cygni in the archives. We recalculated all data (where it was possible) using 36 Cygni as a comparison star. We are presenting UBV light curves of the variable, and also observations made by V. Nikonov in Abastumani in the period 1935–1937
This thesis investigates the application of polyelectrolyte multilayers in plasmonics and picosecond acoustics. The observed samples were fabricated by the spin-assisted layer-by-layer deposition technique that allowed a precise tuning of layer thickness in the range of few nanometers.
The first field of interest deals with the interaction of light-induced localized surface plasmons (LSP) of rod-shaped gold nanoparticles with the particles' environment. The environment consists of an air phase and a phase of polyelectrolytes, whose ratio affects the spectral position of the LSP resonance.
Measured UV-VIS spectra showed the shift of the LSP absorption peak as a function of the cover layer thickness of the particles. The data are modeled using an average dielectric function instead of the dielectric functions of air and polyelectrolytes. In addition using a measured dielectric function of the gold nanoparticles, the position of the LSP absorption peak could be simulated with good agreement to the data.
The analytic model helps to understand the optical properties of metal nanoparticles in an inhomogeneous environment.
The second part of this work discusses the applicability of PAzo/PAH and dye-doped PSS/PAH polyelectrolyte multilayers as transducers to generate hypersound pulses. The generated strain pulses were detected by time-domain Brillouin scattering (TDBS) using a pump-probe laser setup. Transducer layers made of polyelectrolytes were compared qualitatively to common aluminum transducers in terms of measured TDBS signal amplitude, degradation due to laser excitation, and sample preparation.
The measurements proved that fast and easy prepared polyelectrolyte transducers provided stronger TDBS signals than the aluminum transducer. AFM topography measurements showed a degradation of the polyelectrolyte structures, especially for the PAzo/PAH sample.
To quantify the induced strain, optical barriers were introduced to separate the transducer material from the medium of the hypersound propagation. Difficulties in the sample preparation prohibited a reliable quantification. But the experiments showed that a coating with transparent polyelectrolytes increases the efficiency of aluminum transducers and modifies the excited phonon distribution.
The adoption of polyelectrolytes to the scientific field of picosecond acoustics enables a cheap and fast fabrication of transducer layers on most surfaces. In contrast to aluminum layers the polyelectrolytes are transparent over a wide spectral range. Thus, the strain modulation can be probed from surface and back.
Wolf-Rayet (WR) stars lose copious amounts of mass and momentum through dense stellar winds. The interaction of these outflows with their surroundings results in highly structured and complex circumstellar environments, often featuring knots, arcs, shells and spirals. Recent improvements in computational power and techniques have led to the development of detailed, multi-dimensional simulations that have given new insight into the origin of these structures, and better understanding of the physical mechanisms driving their formation. We review three of the main mechanisms that shape the outflows of WR stars:
• interaction with the interstellar medium (ISM), i.e., wind-ISM interactions;
• interaction with a stellar wind, either from a previous phase of evolution or the wind from a companion star, i.e., wind-wind interactions;
• and interaction with a companion star that has a weak or insignificant outflow (e.g., a compact companion such as a neutron star or black hole), i.e.,wind-companion interactions.
We also highlight the broader implications and impact of these circumstellar structures for related phenomena, e.g., for X-ray binaries and Gamma-ray bursts.
The plasmon resonance of metal nanoparticles determines their optical response in the visible spectral range. Many details such as the electronic properties of gold near the particle surface and the local environment of the particles influence the spectra. We show how the cheap but highly precise fabrication of composite nanolayers by spin-assisted layer-by-layer deposition of polyelectrolytes can be used to investigate the spectral response of gold nanospheres (GNS) and gold nanorods (GNR) in a self-consistent way, using the established Maxwell–Garnett effective medium (MGEM) theory beyond the limit of homogeneous media. We show that the dielectric function of gold nanoparticles differs from the bulk value and experimentally characterize the shape and the surrounding of the particles thoroughly by SEM, AFM and ellipsometry. Averaging the dielectric functions of the layered surrounding by an appropriate weighting with the electric field intensity yields excellent agreement for the spectra of several nanoparticles and nanorods with various cover-layer thicknesses.
We have investigated the electrochemical, spectroscopic and electroluminescent properties of a family of aza-aromatic complexes of ruthenium of type [RuII(bpy/phen)2(L)]2+ (4d6) with three isomeric L ligands, where, bpy = 2,2′-bipyridine, phen = 1,10-phenanthroline and the L ligands are 3-(2-pyridyl)[1,2,4]triazolo[1,5-a]pyridine (L1), 3-(2-pyridyl[1,2,3])triazolo[1,5-a]pyridine (L2) and 2-(2-pyridyl)[1,2,4]triazolo[1,5-a]pyridine (L3). The complexes display two bands in the visible region near 410–420 and 440–450 nm. The complexes are diamagnetic and show well defined 1H NMR lines. They are electroactive in acetonitrile solution and exhibit a well defined RuII/RuIII couple near 1.20 to 1.30 V and −1.40 to −1.50 V due to ligand reduction versus Saturated Calomel Electrode (SCE). The solutions are also luminescent, with peaks are near 600 nm. All the complexes are electroluminescent in nature with peaks lying near 580 nm. L1 and L3 ligated complexes with two bpy co-ligands show weak photoluminescence (PL) but stronger electroluminescence (EL) compared to corresponding L2 ligated analogues.
The evolution of massive stars is strongly influenced by their initial chemical composition. We have computed rapidly-rotating massive star models with low metallicity (∼1/50 Z⊙) that evolve chemically homogeneously and have optically-thin winds during the main sequence evolution. These luminous and hot stars are predicted to emit intense mid- and far-UV radiation, but without the broad emission lines that characterize WR stars with optically-thick winds. We show that such Transparent Wind UV-Intense (TWUIN) stars may be responsible for the high number of He ii ionizing photons observed in metal-poor dwarf galaxies, such as IZw 18. We find that these TWUIN stars are possible long-duration gamma-ray burst progenitors.
Synchronization of large ensembles of oscillators is an omnipresent phenomenon observed in different fields of science like physics, engineering, life sciences, etc. The most simple setup is that of globally coupled phase oscillators, where all the oscillators contribute to a global field which acts on all oscillators. This formulation of the problem was pioneered by Winfree and Kuramoto. Such a setup gives a possibility for the analysis of these systems in terms of global variables. In this work we describe nontrivial collective dynamics in oscillator populations coupled via mean fields in terms of global variables. We consider problems which cannot be directly reduced to standard Kuramoto and Winfree models.
In the first part of the thesis we adopt a method introduced by Watanabe and Strogatz. The main idea is that the system of identical oscillators of particular type can be described by a low-dimensional system of global equations. This approach enables us to perform a complete analytical analysis for a special but vast set of initial conditions. Furthermore, we show how the approach can be expanded for some nonidentical systems. We apply the Watanabe-Strogatz approach to arrays of Josephson junctions and systems of identical phase oscillators with leader-type coupling.
In the next parts of the thesis we consider the self-consistent mean-field theory method that can be applied to general nonidentical globally coupled systems of oscillators both with or without noise. For considered systems a regime, where the global field rotates uniformly, is the most important one. With the help of this approach such solutions of the self-consistency equation for an arbitrary distribution of frequencies and coupling parameters can be found analytically in the parametric form, both for noise-free and noisy cases.
We apply this method to deterministic Kuramoto-type model with generic coupling and an ensemble of spatially distributed oscillators with leader-type coupling. Furthermore, with the proposed self-consistent approach we fully characterize rotating wave solutions of noisy Kuramoto-type model with generic coupling and an ensemble of noisy oscillators with bi-harmonic coupling.
Whenever possible, a complete analysis of global dynamics is performed and compared with direct numerical simulations of large populations.
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.
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.
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.
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.