Institut für Physik und Astronomie
Refine
Year of publication
- 2024 (16)
- 2023 (48)
- 2022 (298)
- 2021 (173)
- 2020 (246)
- 2019 (312)
- 2018 (312)
- 2017 (281)
- 2016 (258)
- 2015 (221)
- 2014 (237)
- 2013 (185)
- 2012 (218)
- 2011 (187)
- 2010 (90)
- 2009 (115)
- 2008 (29)
- 2007 (81)
- 2006 (171)
- 2005 (191)
- 2004 (205)
- 2003 (102)
- 2002 (138)
- 2001 (160)
- 2000 (162)
- 1999 (208)
- 1998 (168)
- 1997 (149)
- 1996 (116)
- 1995 (177)
- 1994 (84)
- 1993 (10)
- 1992 (8)
- 1991 (2)
- 1989 (1)
- 1972 (1)
Document Type
- Article (4152)
- Doctoral Thesis (776)
- Postprint (126)
- Monograph/Edited Volume (120)
- Other (85)
- Review (36)
- Habilitation Thesis (24)
- Preprint (24)
- Master's Thesis (10)
- Conference Proceeding (6)
Is part of the Bibliography
- yes (5362) (remove)
Keywords
- diffusion (58)
- stars: massive (57)
- gamma rays: general (47)
- anomalous diffusion (46)
- stars: early-type (45)
- cosmic rays (43)
- stars: winds, outflows (42)
- Magellanic Clouds (39)
- radiation mechanisms: non-thermal (38)
- ISM: supernova remnants (37)
Institute
- Institut für Physik und Astronomie (5362)
- Extern (61)
- Institut für Chemie (18)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (18)
- Institut für Geowissenschaften (10)
- Mathematisch-Naturwissenschaftliche Fakultät (7)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (7)
- Department Psychologie (6)
- Institut für Biochemie und Biologie (4)
- Zentrum für Lehrerbildung und Bildungsforschung (ZeLB) (4)
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.
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 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.
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.
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.
Under the assumption of spherical symmetry, the run of intensity with impact parameter for a spatially resolved and optically thin bubble can be inverted for an "effective emissivity" as a function of radius. The effective emissivity takes into account instrumental sensitivity and even interstellar absorption. This work was supported by a grant from NASA (G03-14008X).
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.
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.
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.
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 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.
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).
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.
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.
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.
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.
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.
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.
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.
Function by structure
(2015)
Forcing Earth’s sea level
(2015)
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.
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.
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.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
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.
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.
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.
Atmospheric interactions with land surface in the arctic based on regional climate model solutions
(2014)
Scientific inquiry requires that we formulate not only what we know, but also what we do not know and by how much. In climate data analysis, this involves an accurate specification of measured quantities and a consequent analysis that consciously propagates the measurement errors at each step. The dissertation presents a thorough analytical method to quantify errors of measurement inherent in paleoclimate data. An additional focus are the uncertainties in assessing the coupling between different factors that influence the global mean temperature (GMT).
Paleoclimate studies critically rely on `proxy variables' that record climatic signals in natural archives. However, such proxy records inherently involve uncertainties in determining the age of the signal. We present a generic Bayesian approach to analytically determine the proxy record along with its associated uncertainty, resulting in a time-ordered sequence of correlated probability distributions rather than a precise time series. We further develop a recurrence based method to detect dynamical events from the proxy probability distributions. The methods are validated with synthetic examples and
demonstrated with real-world proxy records. The proxy estimation step reveals the interrelations between proxy variability and uncertainty. The recurrence analysis of the East Asian Summer Monsoon during the last 9000 years confirms the well-known `dry' events at 8200 and 4400 BP, plus an additional significantly dry event at 6900 BP.
We also analyze the network of dependencies surrounding GMT. We find an intricate, directed network with multiple links between the different factors at multiple time delays. We further uncover a significant feedback from the GMT to the El Niño Southern Oscillation at quasi-biennial timescales. The analysis highlights the need of a more nuanced formulation of influences between different climatic factors, as well as the limitations in trying to estimate such dependencies.
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.
It is generally agreed upon that stars typically form in open clusters and stellar associations, but little is known about the structure of the open cluster system. Do open clusters and stellar associations form isolated or do they prefer to form in groups and complexes? Open cluster groups and complexes could verify star forming regions to be larger than expected, which would explain the chemical homogeneity over large areas in the Galactic disk. They would also define an additional level in the hierarchy of star formation and could be used as tracers for the scales of fragmentation in giant molecular clouds? Furthermore, open cluster groups and complexes could affect Galactic dynamics and should be considered in investigations and simulations on the dynamical processes, such as radial migration, disc heating, differential rotation, kinematic resonances, and spiral structure.
In the past decade there were a few studies on open cluster pairs (de La Fuente Marcos & de La Fuente Marcos 2009a,b,c) and on open cluster groups and complexes (Piskunov et al. 2006). The former only considered spatial proximity for the identification of the pairs, while the latter also required tangential velocities to be similar for the members. In this work I used the full set of 6D phase-space information to draw a more detailed picture on these structures. For this purpose I utilised the most homogeneous cluster catalogue available, namely the Catalogue of Open Cluster Data (COCD; Kharchenko et al. 2005a,b), which contains parameters for 650 open clusters and compact associations, as well as for their uniformly selected members. Additional radial velocity (RV) and metallicity ([M/H]) information on the members were obtained from the RAdial Velocity Experiment (RAVE; Steinmetz et al. 2006; Kordopatis et al. 2013) for 110 and 81 clusters, respectively. The RAVE sample was cleaned considering quality parameters and flags provided by RAVE (Matijevič et al. 2012; Kordopatis et al. 2013). To ensure that only real members were included for the mean values, also the cluster membership, as provided by Kharchenko et al. (2005a,b), was considered for the stars cross-matched in RAVE.
6D phase-space information could be derived for 432 out of the 650 COCD objects and I used an adaption of the Friends-of-Friends algorithm, as used in cosmology, to identify potential groupings. The vast majority of the 19 identified groupings were pairs, but I also found four groups of 4-5 members and one complex with 15 members. For the verification of the identified structures, I compared the results to a randomly selected subsample of the catalogue for the Milky Way global survey of Star Clusters (MWSC; Kharchenko et al. 2013), which became available recently, and was used as reference sample. Furthermore, I implemented Monte-Carlo simulations with randomised samples created from two distinguished input distributions for the spatial and velocity parameters. On the one hand, assuming a uniform distribution in the Galactic disc and, on the other hand, assuming the COCD data distributions to be representative for the whole open cluster population.
The results suggested that the majority of identified pairs are rather by chance alignments, but the groups and the complex seemed to be genuine. A comparison of my results to the pairs, groups and complexes proposed in the literature yielded a partial overlap, which was most likely because of selection effects and different parameters considered. This is another verification for the existence of such structures.
The characteristics of the found groupings favour that members of an open cluster grouping originate from a common giant molecular cloud and formed in a single, but possibly sequential, star formation event. Moreover, the fact that the young open cluster population showed smaller spatial separations between nearest neighbours than the old cluster population indicated that the lifetime of open cluster groupings is most likely comparable to that of the Galactic open cluster population itself. Still even among the old open clusters I could identify groupings, which suggested that the detected structure could be in some cases more long lived as one might think.
In this thesis I could only present a pilot study on structures in the Galactic open cluster population, since the data sample used was highly incomplete. For further investigations a far more complete sample would be required. One step in this direction would be to use data from large current surveys, like SDSS, RAVE, Gaia-ESO and VVV, as well as including results from studies on individual clusters. Later the sample can be completed by data from upcoming missions, like Gaia and 4MOST. Future studies using this more complete open cluster sample will reveal the effect of open cluster groupings on star formation theory and their significance for the kinematics, dynamics and evolution of the Milky Way, and thereby of spiral galaxies.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
The Epoch of Reionization marks after recombination the second major change in the ionization state of the universe, going from a neutral to an ionized state. It starts with the appearance of the first stars and galaxies; a fraction of high-energy photons emitted from galaxies permeate into the intergalactic medium (IGM) and gradually ionize the hydrogen, until the IGM is completely ionized at z~6 (Fan et al., 2006). While the progress of reionization is driven by galaxy evolution, it changes the ionization and thermal state of the IGM substantially and affects subsequent structure and galaxy formation by various feedback mechanisms.
Understanding this interaction between reionization and galaxy formation is further impeded by a lack of understanding of the high-redshift galactic properties such as the dust distribution and the escape fraction of ionizing photons. Lyman Alpha Emitters (LAEs) represent a sample of high-redshift galaxies that are sensitive to all these galactic properties and the effects of reionization.
In this thesis we aim to understand the progress of reionization by performing cosmological simulations, which allows us to investigate the limits of constraining reionization by high-redshift galaxies as LAEs, and examine how galactic properties and the ionization state of the IGM affect the visibility and observed quantities of LAEs and Lyman Break galaxies (LBGs).
In the first part of this thesis we focus on performing radiative transfer calculations to simulate reionization. We have developed a mapping-sphere-scheme, which, starting from spherically averaged temperature and density fields, uses our 1D radiative transfer code and computes the effect of each source on the IGM temperature and ionization (HII, HeII, HeIII) profiles, which are subsequently mapped onto a grid. Furthermore we have updated the 3D Monte-Carlo radiative transfer pCRASH, enabling detailed reionization simulations which take individual source characteristics into account.
In the second part of this thesis we perform a reionization simulation by post-processing a smoothed-particle hydrodynamical (SPH) simulation (GADGET-2) with 3D radiative transfer (pCRASH), where the ionizing sources are modelled according to the characteristics of the stellar populations in the hydrodynamical simulation. Following the ionization fractions of hydrogen (HI) and helium (HeII, HeIII), and temperature in our simulation, we find that reionization starts at z~11 and ends at z~6, and high density regions near sources are ionized earlier than low density regions far from sources.
In the third part of this thesis we couple the cosmological SPH simulation and the radiative transfer simulations with a physically motivated, self-consistent model for LAEs, in order to understand the importance of the ionization state of the IGM, the escape fraction of ionizing photons from galaxies and dust in the interstellar medium (ISM) on the visibility of LAEs. Comparison of our models results with the LAE Lyman Alpha (Lya) and UV luminosity functions at z~6.6 reveals a three-dimensional degeneracy between the ionization state of the IGM, the ionizing photons escape fraction and the ISM dust distribution, which implies that LAEs act not only as tracers of reionization but also of the ionizing photon escape fraction and of the ISM dust distribution. This degeneracy does not even break down when we compare simulated with observed clustering of LAEs at z~6.6. However, our results show that reionization has the largest impact on the amplitude of the LAE angular correlation functions, and its imprints are clearly distinguishable from those of properties on galactic scales. These results show that reionization cannot be constrained tightly by exclusively using LAE observations. Further observational constraints, e.g. tomographies of the redshifted hydrogen 21cm line, are required.
In addition we also use our LAE model to probe the question when a galaxy is visible as a LAE or a LBG. Within our model galaxies above a critical stellar mass can produce enough luminosity to be visible as a LBG and/or a LAE. By finding an increasing duty cycle of LBGs with Lya emission as the UV magnitude or stellar mass of the galaxy rises, our model reveals that the brightest (and most massive) LBGs most often show Lya emission.
Predicting the Lya equivalent width (Lya EW) distribution and the fraction of LBGs showing Lya emission at z~6.6, we reproduce the observational trend of the Lya EWs with UV magnitude. However, the Lya EWs of the UV brightest LBGs exceed observations and can only be reconciled by accounting for an increased Lya attenuation of massive galaxies, which implies that the observed Lya brightest LAEs do not necessarily coincide with the UV brightest galaxies. We have analysed the dependencies of LAE observables on the properties of the galactic and intergalactic medium and the LAE-LBG connection, and this enhances our understanding of the nature of LAEs.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.