Refine
Year of publication
- 2016 (258) (remove)
Document Type
- Article (203)
- Doctoral Thesis (21)
- Postprint (18)
- Other (10)
- Review (4)
- Habilitation Thesis (2)
Is part of the Bibliography
- yes (258)
Keywords
- Magellanic Clouds (8)
- Sun: magnetic fields (6)
- intergalactic medium (6)
- methods: data analysis (6)
- stars: Wolf-Rayet (6)
- techniques: spectroscopic (6)
- galaxies: formation (5)
- stars: atmospheres (5)
- turbulence (5)
- X-rays: binaries (4)
Institute
- Institut für Physik und Astronomie (258) (remove)
We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.
We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.
Particle diffusion and localized acceleration in inhomogeneous AGN jets - II. Stochastic variation
(2016)
We study the stochastic variation of blazar emission under a 2D spatially resolved leptonic jet model we previously developed. Random events of particle acceleration and injection in small zones within the emission region are assumed to be responsible for flux variations. In addition to producing spectral energy distributions that describe the observed flux of Mrk 421, we further analyse the timing properties of the simulated light curves, such as the power spectral density (PSD) at different bands, flux-flux correlations, aswell as the cross-correlation function between X-rays and TeV gamma-rays. We find spectral breaks in the PSD at a time-scale comparable to the dominant characteristic time-scale in the system, which is usually the predefined decay time-scale of an acceleration event. Cooling imposes a delay, and so PSDs taken at lower energy bands in each emission component (synchrotron or inverse Compton) generally break at longer time-scales. The flux-flux correlation between X-rays and TeV gamma-rays can be either quadratic or linear, depending on whether or not there are large variation of the injection into the particle acceleration process. When the relationship is quadratic, the TeV flares lag the X-ray flares, and the optical and GeV flares are large enough to be comparable to the ones in X-ray. When the relationship is linear, the lags are insignificant, and the optical and GeV flares are small.
Adsorption of amino acids on the magnetite-(111)-surface: a force field study (vol 19, 851, 2013)
(2016)
Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback.
This work reports about new high-resolution imaging and spectroscopic observations of solar type III radio bursts at low radio frequencies in the range from 30 to 80 MHz. Solar type III radio bursts are understood as result of the beam-plasma interaction of electron beams in the corona. The Sun provides a unique opportunity to study these plasma processes of an active star. Its activity appears in eruptive events like flares, coronal mass ejections and radio bursts which are all accompanied by enhanced radio emission. Therefore solar radio emission carries important information about plasma processes associated with the Sun’s activity. Moreover, the Sun’s atmosphere is a unique plasma laboratory with plasma processes under conditions not found in terrestrial laboratories. Because of the Sun’s proximity to Earth, it can be studied in greater detail than any other star but new knowledge about the Sun can be transfer to them. This “solar stellar connection” is important for the understanding of processes on other stars.
The novel radio interferometer LOFAR provides imaging and spectroscopic capabilities to study these processes at low frequencies. Here it was used for solar observations.
LOFAR, the characteristics of its solar data and the processing and analysis of the latter with the Solar Imaging Pipeline and Solar Data Center are described. The Solar Imaging Pipeline is the central software that allows using LOFAR for solar observations. So its development was necessary for the analysis of solar LOFAR data and realized here. Moreover a new density model with heat conduction and Alfvén waves was developed that provides the distance of radio bursts to the Sun from dynamic radio spectra.
Its application to the dynamic spectrum of a type III burst observed on March 16, 2016 by LOFAR shows a nonuniform radial propagation velocity of the radio emission. The analysis of an imaging observation of type III bursts on June 23, 2012 resolves a burst as bright, compact region localized in the corona propagating in radial direction along magnetic field lines with an average velocity of 0.23c. A nonuniform propagation velocity is revealed. A new beam model is presented that explains the nonuniform motion of the radio source as a propagation effect of an electron ensemble with a spread velocity distribution and rules out a monoenergetic electron distribution. The coronal electron number density is derived in the region from 1.5 to 2.5 R☉ and fitted with the newly developed density model. It determines the plasma density for the interplanetary space between Sun and Earth. The values correspond to a 1.25- and 5-fold Newkirk model for harmonic and fundamental emission, respectively. In comparison to data from other radio instruments the LOFAR data shows a high sensitivity and resolution in space, time and frequency.
The new results from LOFAR’s high resolution imaging spectroscopy are consistent with current theories of solar type III radio bursts and demonstrate its capability to track fast moving radio sources in the corona. LOFAR solar data is found to be a valuable source for solar radio physics and opens a new window for studying plasma processes associated with highly energetic electrons in the solar corona.
The Milky Way is only one out of billions of galaxies in the universe. However, it is a special galaxy because it allows to explore the main mechanisms involved in its evolution and formation history by unpicking the system star-by-star. Especially, the chemical fingerprints of its stars provide clues and evidence of past events in the Galaxy’s lifetime. These information help not only to decipher the current structure and building blocks of the Milky Way, but to learn more about the general formation process of galaxies.
In the past decade a multitude of stellar spectroscopic Galactic surveys have scanned millions of stars far beyond the rim of the solar neighbourhood. The obtained spectroscopic information provide unprecedented insights to the chemo-dynamics of the Milky Way. In addition analytic models and numerical simulations of the Milky Way provide necessary descriptions and predictions suited for comparison with observations in order to decode the physical properties that underlie the complex system of the Galaxy.
In the thesis various approaches are taken to connect modern theoretical modelling of galaxy formation and evolution with observations from Galactic stellar surveys. With its focus on the chemo-kinematics of the Galactic disk this work aims to determine new observational constraints on the formation of the Milky Way providing also proper comparisons with two different models. These are the population synthesis model TRILEGAL based on analytical distribution functions, which aims to simulate the number and distribution of stars in the Milky Way and its different components, and a hybrid model (MCM) that combines an N-body simulation of a Milky Way like galaxy in the cosmological framework with a semi-analytic chemical evolution model for the Milky Way. The major observational data sets in use come from two surveys, namely the “Radial Velocity Experiment” (RAVE) and the “Sloan Extension for Galactic Understanding and Exploration” (SEGUE).
In the first approach the chemo-kinematic properties of the thin and thick disk of the Galaxy as traced by a selection of about 20000 SEGUE G-dwarf stars are directly compared to the predictions by the MCM model. As a necessary condition for this, SEGUE's selection function and its survey volume are evaluated in detail to correct the spectroscopic observations for their survey specific selection biases. Also, based on a Bayesian method spectro-photometric distances with uncertainties below 15% are computed for the selection of SEGUE G-dwarfs that are studied up to a distance of 3 kpc from the Sun.
For the second approach two synthetic versions of the SEGUE survey are generated based on the above models. The obtained synthetic stellar catalogues are then used to create mock samples best resembling the compiled sample of observed SEGUE G-dwarfs. Generally, mock samples are not only ideal to compare predictions from various models. They also allow validation of the models' quality and improvement as with this work could be especially achieved for TRILEGAL. While TRILEGAL reproduces the statistical properties of the thin and thick disk as seen in the observations, the MCM model has shown to be more suitable in reproducing many chemo-kinematic correlations as revealed by the SEGUE stars. However, evidence has been found that the MCM model may be missing a stellar component with the properties of the thick disk that the observations clearly show. While the SEGUE stars do indicate a thin-thick dichotomy of the stellar Galactic disk in agreement with other spectroscopic stellar studies, no sign for a distinct metal-poor disk is seen in the MCM model.
Usually stellar spectroscopic surveys are limited to a certain volume around the Sun covering different regions of the Galaxy’s disk. This often prevents to obtain a global view on the chemo-dynamics of the Galactic disk. Hence, a suitable combination of stellar samples from independent surveys is not only useful for the verification of results but it also helps to complete the picture of the Milky Way. Therefore, the thesis closes with a comparison of the SEGUE G-dwarfs and a sample of RAVE giants. The comparison reveals that the chemo-kinematic relations agree in disk regions where the samples of both surveys show a similar number of stars. For those parts of the survey volumes where one of the surveys lacks statistics they beautifully complement each other. This demonstrates that the comparison of theoretical models on the one side, and the combined observational data gathered by multiple surveys on the other side, are key ingredients to understand and disentangle the structure and formation history of the Milky Way.
The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the nonstationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total timescale of several hours), the transitions of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the nonstationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favored in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.
In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting approximately 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high, we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detections of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousand seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.
We present a statistical analysis of phase space density data from the first 26 months of the Van Allen Probes mission. In particular, we investigate the relationship between the tens and hundreds of keV seed electrons and >1 MeV core radiation belt electron population. Using a cross-correlation analysis, we find that the seed and core populations are well correlated with a coefficient of approximate to 0.73 with a time lag of 10-15 h. We present evidence of a seed population threshold that is necessary for subsequent acceleration. The depth of penetration of the seed population determines the inner boundary of the acceleration process. However, we show that an enhanced seed population alone is not enough to produce acceleration in the higher energies, implying that the seed population of hundreds of keV electrons is only one of several conditions required for MeV electron radiation belt acceleration.
We use a background quasar to detect the presence of circumgalactic gas around a z = 0.91 low-mass star-forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope show that the galaxy has a dust-corrected star formation rate (SFR) of 4.7 +/- 2.0. M-circle dot yr(-1), with no companion down to 0.22 M-circle dot yr(-1) (5 sigma) within 240 h(-1) kpc ("30"). Using a high-resolution spectrum of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle alpha of only 15 degrees), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold-flow disk" extending at least 12 kpc (3 x R-1/2). We estimate the mass accretion rate M-in to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the H (I) column density of log N-H (I)/cm(-2) similar or equal to 20.4 obtained from a Hubble Space Telescope/COS near-UV spectrum. From a detailed analysis of the low-ionization lines (e.g., Zn II, Cr II, Ti II, MnII, Si II), the accreting material appears to be enriched to about 0.4 Z(circle dot) (albeit with large uncertainties: log Z/Z(circle dot) = -0.4 +/- 0.4), which is comparable to the galaxy metallicity (12 + log O/H = 8.7 +/- 0.2), implying a large recycling fraction from past outflows. Blueshifted Mg II and Fe II absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The Mg II and Fe II absorption line ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent Fe II* emission.
We consider two coupled populations of leaky integrate-and-fire neurons. Depending on the coupling strength, mean fields generated by these populations can have incommensurate frequencies or become frequency locked. In the observed 2:1 locking state of the mean fields, individual neurons in one population are asynchronous with the mean fields, while in another population they have the same frequency as the mean field. These synchronous neurons form a chimera state, where part of them build a fully synchronized cluster, while other remain scattered. We explain this chimera as a marginal one, caused by a self-organized neutral dynamics of the effective circle map.
This publications-based thesis summarizes my contribution to the scientific field of ultrafast structural dynamics. It consists of 16 publications, about the generation, detection and coupling of coherent gigahertz longitudinal acoustic phonons, also called hypersonic waves. To generate such high frequency phonons, femtosecond near infrared laser pulses were used to heat nanostructures composed of perovskite oxides on an ultrashort timescale. As a consequence the heated regions of such a nanostructure expand and a high frequency acoustic phonon pulse is generated. To detect such coherent acoustic sound pulses I use ultrafast variants of optical Brillouin and x-ray scattering. Here an incident optical or x-ray photon is scattered by the excited sound wave in the sample. The scattered light intensity measures the occupation of the phonon modes.
The central part of this work is the investigation of coherent high amplitude phonon wave packets which can behave nonlinearly, quite similar to shallow water waves which show a steepening of wave fronts or solitons well known as tsunamis. Due to the high amplitude of the acoustic wave packets in the solid, the acoustic properties can change significantly in the vicinity of the sound pulse. This may lead to a shape change of the pulse. I have observed by time-resolved Brillouin scattering, that a single cycle hypersound pulse shows a wavefront steepening. I excited hypersound pulses with strain amplitudes until 1% which I have calibrated by ultrafast x-ray diffraction (UXRD).
On the basis of this first experiment we developed the idea of the nonlinear mixing of narrowband phonon wave packets which we call "nonlinear phononics" in analogy with the nonlinear optics, which summarizes a kaleidoscope of surprising optical phenomena showing up at very high electric fields. Such phenomena are for instance Second Harmonic Generation, four-wave-mixing or solitons. But in case of excited coherent phonons the wave packets have usually very broad spectra which make it nearly impossible to look at elementary scattering processes between phonons with certain momentum and energy.
For that purpose I tested different techniques to excite narrowband phonon wave packets which mainly consist of phonons with a certain momentum and frequency. To this end epitaxially grown metal films on a dielectric substrate were excited with a train of laser pulses. These excitation pulses drive the metal film to oscillate with the frequency given by their inverse temporal displacement and send a hypersonic wave of this frequency into the substrate. The monochromaticity of these wave packets was proven by ultrafast optical Brillouin and x-ray scattering.
Using the excitation of such narrowband phonon wave packets I was able to observe the Second Harmonic Generation (SHG) of coherent phonons as a first example of nonlinear wave mixing of nanometric phonon wave packets.
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion
(2016)
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
The concept of bonding and antibonding orbitals is fundamental in chemistry. The population of those orbitals and the energetic difference between the two reflect the strength of the bonding interaction. Weakening the bond is expected to reduce this energetic splitting, but the transient character of bond-activation has so far prohibited direct experimental access. Here we apply time-resolved soft X-ray spectroscopy at a free electron laser to directly observe the decreased bonding antibonding splitting following bond-activation using an ultrashort optical laser pulse.
To turn or not to turn?
(2016)
Bacteria typically swim in straight runs, interruped by sudden turning events. In particular, some species are limited to a reversal in the swimming direction as the only turning maneuver at their disposal. In a recent article, Grossmann et al (2016 New J. Phys. 18 043009) introduce a theoretical framework to analyze the diffusive properties of active particles following this type of run-and-reverse pattern. Based on a stochastic clock model to mimic the regulatory pathway that triggers reversal events, they show that a run-and-reverse swimmer can optimize its diffusive spreading by tuning the reversal rate according to the level of rotational noise. With their approach, they open up promising new perspectives of how to incorporate the dynamics of intracellular signaling into coarse-grained active particle descriptions.
In this Comment, we review the results of pattern formation in a reaction-diffusion-advection system following the kinetics of the Gray-Scott model. A recent paper by Das [Phys. Rev. E 92, 052914 (2015)] shows that spatiotemporal chaos of the intermittency type can disappear as the advective flow is increased. This study, however, refers to a single point in the space of kinetic parameters of the original Gray-Scott model. Here we show that the wealth of patterns increases substantially as some of these parameters are changed. In addition to spatiotemporal intermittency, defect-mediated turbulence can also be found. In all cases, however, the chaotic behavior is seen to disappear as the advective flow is increased, following a scenario similar to what was reported in our earlier work [I. Berenstein and C. Beta, Phys. Rev. E 86, 056205 (2012)] as well as by Das. We also point out that a similar phenomenon can be found in other reaction-diffusion-advection models, such as the Oregonator model for the Belousov-Zhabotinsky reaction under flow conditions.
We do magnetohydrodynamic (MHD) simulations of local box models of turbulent Interstellar Medium (ISM) and analyse the process of amplification and saturation of mean magnetic fields with methods of mean field dynamo theory. It is shown that the process of saturation of mean fields can be partially described by the prolonged diffusion time scales in presence of the dynamically significant magnetic fields. However, the outward wind also plays an essential role in the saturation in higher SN rate case. Algebraic expressions for the back reaction of the magnetic field onto the turbulent transport coefficients are derived, which allow a complete description of the nonlinear dynamo. We also present the effects of dynamically significant mean fields on the ISM configuration and pressure distribution. We further add the cosmic ray component in the simulations and investigate the kinematic growth of mean fields with a dynamo perspective.
A very general cosmological consideration suggests that, along with galactic dark matter halos, much smaller dark matter structures may exist. These structures are usually called `clumps', and their mass extends to 10−6 M ⊙ or even lower. The clumps should give the main contribution into the signal of dark matter annihilation, provided that they have survived until the present time. Recent observations favor a cored profile for low-mass astrophysical halos. We consider cored clumps and show that they are significantly less firm than the standard NFW ones. In contrast to the standard scenario, the cored clumps should have been completely destroyed inside ~ 20 kpc from the Milky Way center. The dwarf spheroidals should not contain any dark matter clumps. On the other hand, even under the most pessimistic assumption about the clump structure, the clumps should have survived in the Milky Way at a distance exceeding 50 kpc from the center, as well as in low-density cosmic structures. There they significantly boost the dark matter annihilation. We show that at least 70% of the clumps endured the primordial structure formation should still exist untouched in the present-day Universe.
Traditional economic theory could not explain, much less predict, the near collapse of the financial system and its long-lasting effects on the global economy. Since the 2008 crisis, there has been increasing interest in using ideas from complexity theory to make sense of economic and financial markets. Concepts, such as tipping points, networks, contagion, feedback, and resilience have entered the financial and regulatory lexicon, but actual use of complexity models and results remains at an early stage. Recent insights and techniques offer potential for better monitoring and management of highly interconnected economic and financial systems and, thus, may help anticipate and manage future crises.
Arch filament systems occur in active sunspot groups, where a fibril structure connects areas of opposite magnetic polarity, in contrast to active region filaments that follow the polarity inversion line. We used the GREGOR Infrared Spectrograph (GRIS) to obtain the full Stokes vector in the spectral lines SiI lambda 1082.7 nm, He I lambda 1083.0 nm, and Ca I lambda 1083.9 nm. We focus on the near-infrared calcium line to investigate the photospheric magnetic field and velocities, and use the line core intensities and velocities of the helium line to study the chromospheric plasma. The individual fibrils of the arch filament system connect the sunspot with patches of magnetic polarity opposite to that of the spot. These patches do not necessarily coincide with pores, where the magnetic field is strongest. Instead, areas are preferred not far from the polarity inversion line. These areas exhibit photospheric downflows of moderate velocity, but significantly higher downflows of up to 30 km s(-1) in the chromospheric helium line. Our findings can be explained with new emerging flux where the matter flows downward along the field lines of rising flux tubes, in agreement with earlier results. (C) 2016 WILEY-VCH Verlag GmbH& Co. KGaA, Weinheim
In the current paradigm of cosmology, the formation of large-scale structures is mainly driven by non-radiating dark matter, making up the dominant part of the matter budget of the Universe. Cosmological observations however, rely on the detection of luminous galaxies, which are biased tracers of the underlying dark matter. In this thesis I present cosmological reconstructions of both, the dark matter density field that forms the cosmic web, and cosmic velocities, for which both aspects of my work are delved into, the theoretical formalism and the results of its applications to cosmological simulations and also to a galaxy redshift survey.The foundation of our method is relying on a statistical approach, in which a given galaxy catalogue is interpreted as a biased realization of the underlying dark matter density field. The inference is computationally performed on a mesh grid by sampling from a probability density function, which describes the joint posterior distribution of matter density and the three dimensional velocity field. The statistical background of our method is described in Chapter ”Implementation of argo”, where the introduction in sampling methods is given, paying special attention to Markov Chain Monte-Carlo techniques. In Chapter ”Phase-Space Reconstructions with N-body Simulations”, I introduce and implement a novel biasing scheme to relate the galaxy number density to the underlying dark matter, which I decompose into a deterministic part, described by a non-linear and scale-dependent analytic expression, and a stochastic part, by presenting a negative binomial (NB) likelihood function that models deviations from Poissonity. Both bias components had already been studied theoretically, but were so far never tested in a reconstruction algorithm. I test these new contributions againstN-body simulations to quantify improvements and show that, compared to state-of-the-art methods, the stochastic bias is inevitable at wave numbers of k≥0.15h Mpc^−1 in the power spectrum in order to obtain unbiased results from the reconstructions. In the second part of Chapter ”Phase-Space Reconstructions with N-body Simulations” I describe and validate our approach to infer the three dimensional cosmic velocity field jointly with the dark matter density. I use linear perturbation theory for the large-scale bulk flows and a dispersion term to model virialized galaxy motions, showing that our method is accurately recovering the real-space positions of the redshift-space distorted galaxies. I analyze the results with the isotropic and also the two-dimensional power spectrum.Finally, in Chapter ”Phase-space Reconstructions with Galaxy Redshift Surveys”, I show how I combine all findings and results and apply the method to the CMASS (for Constant (stellar) Mass) galaxy catalogue of the Baryon Oscillation Spectroscopic Survey (BOSS). I describe how our method is accounting for the observational selection effects inside our reconstruction algorithm. Also, I demonstrate that the renormalization of the prior distribution function is mandatory to account for higher order contributions in the structure formation model, and finally a redshift-dependent bias factor is theoretically motivated and implemented into our method. The various refinements yield unbiased results of the dark matter until scales of k≤0.2 h Mpc^−1in the power spectrum and isotropize the galaxy catalogue down to distances of r∼20h^−1 Mpc in the correlation function. We further test the results of our cosmic velocity field reconstruction by comparing them to a synthetic mock galaxy catalogue, finding a strong correlation between the mock and the reconstructed velocities. The applications of both, the density field without redshift-space distortions, and the velocity reconstructions, are very broad and can be used for improved analyses of the baryonic acoustic oscillations, environmental studies of the cosmic web, the kinematic Sunyaev-Zel’dovic or integrated Sachs-Wolfe effect.
In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.
Methods. In most cases, the given orientation of the ecliptic is used to set up the heliographic coordinate system for the drawings. Positions and sizes are measured manually on screen. Very early drawings have no indication of their orientation. A rotational matching using common spots of adjacent days is used in some cases, while in other cases, the assumption that images were aligned with a zenith-horizon coordinate system appeared to be the most probable.
The Galactic Center ridge has been observed extensively in the past by both GeV and TeV gamma-ray instruments revealing a wealth of structure, including a diffuse component and the point sources G0.9+0.1 (a composite supernova remnant) and Sgr A* (believed to be associated with the supermassive black hole located at the center of our Galaxy). Previous very high energy (VHE) gamma-ray observations with the H.E.S.S.. experiment have also detected an extended TeV gamma-ray component along the Galactic plane in the >300 GeV gamma-ray regime. Here we report on observations of the Galactic Center ridge from 2010 to 2014 by the VERITAS telescope array in the >2 TeV energy range. From these observations we (1) provide improved measurements of the differential energy spectrum for Sgr A* in the >2 TeV gamma-ray regime, (2) provide a detection in the >2 TeV gamma-ray emission from the composite SNR G0.9+0.1 and an improved determination of its multi-TeV gamma-ray energy spectrum, and. (3) report on the detection of VER J1746-289, a localized enhancement of >2 TeV gamma-ray emission along the Galactic plane.
Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E > 100 GeV) gamma-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of similar to 570 hr. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog that are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have spectroscopic distance estimates. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample, we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data set, which shows a 4s excess.
The TeV binary system LS I +61 degrees 303 is known for its regular, non-thermal emission pattern that traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I + 61 degrees. 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I + 61 degrees 303 during the flares, provides constraints on the properties of the accelerator in the source.
Gravity dictates the structure of the whole Universe and, although it is triumphantly described by the theory of General Relativity, it is the force that we least understand in nature. One of the cardinal predictions of this theory are black holes. Massive, dark objects are found in the majority of galaxies. Our own galactic center very contains such an object with a mass of about four million solar masses. Are these objects supermassive black holes (SMBHs), or do we need alternatives? The answer lies in the event horizon, the characteristic that defines a black hole. The key to probe the horizon is to model the movement of stars around a SMBH, and the interactions between them, and look for deviations from real observations. Nuclear star clusters harboring a massive, dark object with a mass of up to ~ ten million solar masses are good testbeds to probe the event horizon of the potential SMBH with stars. The channel for interactions between stars and the central MBH are the fact that (a) compact stars and stellar-mass black holes can gradually inspiral into the SMBH due to the emission of gravitational radiation, which is known as an “Extreme Mass Ratio Inspiral” (EMRI), and (b) stars can produce gases which will be accreted by the SMBH through normal stellar evolution, or by collisions and disruptions brought about by the strong central tidal field. Such processes can contribute significantly to the mass of the SMBH. These two processes involve different disciplines, which combined will provide us with detailed information about the fabric of space and time. In this habilitation I present nine articles of my recent work directly related with these topics.
A SEARCH FOR VERY HIGH ENERGY GAMMA RAYS FROM THE MISSING LINK BINARY PULSAR J1023+0038 WITH VERITAS
(2016)
The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than similar to 2 G before the disappearance of the radio pulsar and greater than similar to 10 G afterward.
The very high energy (VHE; E > 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54 938-54 956), 2009, as part of an extensive multiwavelength campaign from radio to VHE. Strong VHE yray activity was detected on May 1st with Whipple and VERITAS, when the flux (E > 400 GeV) increased to 10 times the preflare baseline flux (3.9 x 10(-11) ph cm(-2) s(-1)), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15. This VHE flare showed a fast flux variation with an increase of a factor similar to 4 in 25 min, and a falling time of similar to 50 min. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
Galactic cosmic rays reach energies of at least a few petaelectronvolts (of the order of 1015 electronvolts). This implies that our Galaxy contains petaelectronvolt accelerators (‘PeVatrons’), but all proposed models of Galactic cosmic-ray accelerators encounter difficulties at exactly these energies. Dozens of Galactic accelerators capable of accelerating particles to energies of tens of teraelectronvolts (of the order of 1013 electronvolts) were inferred from recent γ-ray observations3. However, none of the currently known accelerators—not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays—has shown the characteristic tracers of petaelectronvolt particles, namely, power-law spectra of γ-rays extending without a cut-off or a spectral break to tens of teraelectronvolts4. Here we report deep γ-ray observations with arcminute angular resolution of the region surrounding the Galactic Centre, which show the expected tracer of the presence of petaelectronvolt protons within the central 10 parsecs of the Galaxy. We propose that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts5and an outflow from the Galactic Centre6. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last 106–107 years, and therefore should be considered as a viable alternative to supernova remnants as a source of petaelectronvolt Galactic cosmic rays.
The F-type star KIC. 8462852 has recently been identified as an exceptional target for search for extraterrestrial intelligence (SETI) observations. We describe an analysis methodology for optical SETI, which we have used to analyze nine hours of serendipitous archival observations of KIC. 8462852 made with the VERITAS gamma-ray observatory between 2009 and 2015. No evidence of pulsed optical beacons, above a pulse intensity at the Earth of approximately 1 photon m(-2), is found. We also discuss the potential use of imaging atmospheric Cherenkov telescope arrays in searching for extremely short duration optical transients in general.
We present results from multiwavelength observations of the BL Lacertae object 1ES 1741 + 196, including results in the very high energy gamma-ray regime using the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is well modelled by a power law with a spectral index of 2.7 +/- 0.7(stat) +/- 0.2(syst). The integral flux above 180 GeV is (3.9 +/- 0.8(stat) +/- 1.0(syst)) x 10(-8) m(-2) s(-1), corresponding to 1.6 per cent of the Crab nebula flux on average. The multiwavelength spectral energy distribution of the source suggests that 1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The observations analysed in this paper extend over a period of six years, during which time no strong flares were observed in any band. This analysis is therefore one of the few characterizations of a blazar in a non-flaring state.
The inner region of the Milky Way halo harbors a large amount of dark matter (DM). Given its proximity, it is one of the most promising targets to look for DM. We report on a search for the annihilations of DM particles using gamma-ray observations towards the inner 300 pc of the Milky Way, with the H.E.S.S. array of ground-based Cherenkov telescopes. The analysis is based on a 2D maximum likelihood method using Galactic Center (GC) data accumulated by H.E.S.S. over the last 10 years (2004-2014), and does not show any significant gamma-ray signal above background. Assuming Einasto and Navarro-Frenk-White DM density profiles at the GC, we derive upper limits on the annihilation cross section <sigma nu >. These constraints are the strongest obtained so far in the TeV DM mass range and improve upon previous limits by a factor 5. For the Einasto profile, the constraints reach <sigma nu > values of 6 x 10(-26) cm(3) s(-1) in the W+W- channel for a DM particle mass of 1.5 TeV, and 2 x 10(-26) cm(3) s(-1) in the tau(+)tau(-) channel for a 1 TeV mass. For the first time, ground-based gamma-ray observations have reached sufficient sensitivity to probe <sigma nu > values expected from the thermal relic density for TeV DM particles.
A search for dark matter linelike signals iss performed in the vicinity of the Galactic Center by the H.E.S.S. experiment on observational data taken in 2014. An unbinned likelihood analysis iss developed to improve the sensitivity to linelike signals. The upgraded analysis along with newer data extend the energy coverage of the previous measurement down to 100 GeV. The 18 h of data collected with the H.E.S.S. array allow one to rule out at 95% C.L. the presence of a 130 GeV line (at l = -1.5 degrees, b = 0 degrees and for a dark matter profile centered at this location) previously reported in Fermi-LAT data. This new analysis overlaps significantly in energy with previous Fermi-LAT and H.E.S.S. results. No significant excess associated with dark matter annihilations was found in the energy range of 100 GeV to 2 TeV and upper limits on the gamma-ray flux and the velocity weighted annihilation cross section are derived adopting an Einasto dark matter halo profile. Expected limits for present and future large statistics H.E.S.S. observations are also given.