Refine
Document Type
- Article (14)
- Doctoral Thesis (1)
Language
- English (15)
Is part of the Bibliography
- yes (15)
Keywords
- methods: data analysis (15) (remove)
Institute
The Galactic center is an interesting region for high-energy (0.1-100 GeV) and very-high-energy (E > 100 GeV) gamma-ray observations. Potential sources of GeV/TeV gamma-ray emission have been suggested, e.g., the accretion of matter onto the supermassive black hole, cosmic rays from a nearby supernova remnant (e.g., Sgr A East), particle acceleration in a plerion, or the annihilation of dark matter particles. The Galactic center has been detected by EGRET and by Fermi/LAT in the MeV/GeV energy band. At TeV energies, the Galactic center was detected with moderate significance by the CANGAROO and Whipple 10 m telescopes and with high significance by H.E.S.S., MAGIC, and VERITAS. We present the results from three years of VERITAS observations conducted at large zenith angles resulting in a detection of the Galactic center on the level of 18 standard deviations at energies above similar to 2.5 TeV. The energy spectrum is derived and is found to be compatible with hadronic, leptonic, and hybrid emission models discussed in the literature. Future, more detailed measurements of the high-energy cutoff and better constraints on the high-energy flux variability will help to refine and/or disentangle the individual models.
Thermophysical modelling and parameter estimation of small solar system bodies via data assimilation
(2020)
Deriving thermophysical properties such as thermal inertia from thermal infrared observations provides useful insights into the structure of the surface material on planetary bodies. The estimation of these properties is usually done by fitting temperature variations calculated by thermophysical models to infrared observations. For multiple free model parameters, traditional methods such as least-squares fitting or Markov chain Monte Carlo methods become computationally too expensive. Consequently, the simultaneous estimation of several thermophysical parameters, together with their corresponding uncertainties and correlations, is often not computationally feasible and the analysis is usually reduced to fitting one or two parameters. Data assimilation (DA) methods have been shown to be robust while sufficiently accurate and computationally affordable even for a large number of parameters. This paper will introduce a standard sequential DA method, the ensemble square root filter, for thermophysical modelling of asteroid surfaces. This method is used to re-analyse infrared observations of the MARA instrument, which measured the diurnal temperature variation of a single boulder on the surface of near-Earth asteroid (162173) Ryugu. The thermal inertia is estimated to be 295 +/- 18 Jm(-2) K-1 s(-1/2), while all five free parameters of the initial analysis are varied and estimated simultaneously. Based on this thermal inertia estimate the thermal conductivity of the boulder is estimated to be between 0.07 and 0.12,Wm(-1) K-1 and the porosity to be between 0.30 and 0.52. For the first time in thermophysical parameter derivation, correlations and uncertainties of all free model parameters are incorporated in the estimation procedure that is more than 5000 times more efficient than a comparable parameter sweep.
We present K-s-band light curves for 299 Cepheids in the Small Magellanic Cloud (SMC) of which 288 are new discoveries that we have identified using multi-epoch near-infrared photometry obtained by the VISTA survey of the Magellanic Clouds system (VMC). The new Cepheids have periods in the range from 0.34 to 9.1 d and cover the magnitude interval 12.9 <= currency sign < K-s > <= currency sign 17.6 mag. Our method was developed using variable stars previously identified by the optical microlensing survey OGLE. We focus on searching new Cepheids in external regions of the SMC for which complete VMC K-s-band observations are available and no comprehensive identification of different types of variable stars from other surveys exists yet.
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.
Aims. We study the evolution of an arch filament system (AFS) and of its individual arch filaments to learn about the processes occurring in them. Methods. We observed the AFS at the GREGOR solar telescope on Tenerife at high cadence with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) in the He I 10 830 angstrom spectral range. The He I triplet profiles were fitted with analytic functions to infer line-of-sight (LOS) velocities to follow plasma motions within the AFS. Results. We tracked the temporal evolution of an individual arch filament over its entire lifetime, as seen in the He I 10 830 angstrom triplet. The arch filament expanded in height and extended in length from 13 ' to 21 '. The lifetime of this arch filament is about 30 min. About 11 min after the arch filament is seen in He I, the loop top starts to rise with an average Doppler velocity of 6 km s(-1). Only two minutes later, plasma drains down with supersonic velocities towards the footpoints reaching a peak velocity of up to 40 km s(-1) in the chromosphere. The temporal evolution of He I 10 830 angstrom profiles near the leading pore showed almost ubiquitous dual red components of the He I triplet, indicating strong downflows, along with material nearly at rest within the same resolution element during the whole observing time.
We present results of the analysis of 70 RR Lyrae stars located in the bar of the Large Magellanic Cloud (LMC). Combining the spectroscopically determined metallicity of these stars from the literature with precise periods from the OGLE III catalog and multi-epoch K-s photometry from the VISTA survey of the Magellanic Clouds system, we derive a new near-infrared period-luminosity-metallicity (PLKsZ) relation for RR Lyrae variables. In order to fit the relation we use a fitting method developed specifically for this study. The zero-point of the relation is estimated two different ways: by assuming the value of the distance to the LMC and by using Hubble Space Telescope parallaxes of five RR Lyrae stars in the Milky Way (MW). The difference in distance moduli derived by applying these two approaches is similar to 0.2 mag. To investigate this point further we derive the PL(Ks)Z relation based on 23 MW RR Lyrae stars that had been analyzed in Baade-Wesselink studies. We compared the derived PL(Ks)Z relations for RR Lyrae stars in the MW and LMC. Slopes and zero-points are different, but still consistent within the errors. The shallow slope of the metallicity term is confirmed by both LMC and MW variables. The astrometric space mission Gaia is expected to provide a huge contribution to the determination of the RR Lyrae PL(Ks)Z relation; however, calculating an absolute magnitude from the trigonometric parallax of each star and fitting a PL(Ks)Z relation directly to period and absolute magnitude leads to biased results. We present a tool to achieve an unbiased solution by modeling the data and inferring the slope and zero-point of the relation via statistical methods.
Aims. We demonstrate the high multiplex advantage of crowded field 3D spectroscopy with the new integral field spectrograph MUSE by means of a spectroscopic analysis of more than 12 000 individual stars in the globular cluster NGC 6397. Methods. The stars are deblended with a point spread function fitting technique, using a photometric reference catalogue from HST as prior, including relative positions and brightnesses. This catalogue is also used for a first analysis of the extracted spectra, followed by an automatic in-depth analysis via a full-spectrum fitting method based on a large grid of PHOENIX spectra. Results. We analysed the largest sample so far available for a single globular cluster of 18 932 spectra from 12 307 stars in NGC 6397. We derived a mean radial velocity of v(rad) = 17.84 +/- 0.07 km s(-1) and a mean metallicity of [Fe/H] = -2.120 +/- 0.002, with the latter seemingly varying with temperature for stars on the red giant branch (RGB). We determine Teff and [Fe/H] from the spectra, and log g from HST photometry. This is the first very comprehensive Hertzsprung-Russell diagram (HRD) for a globular cluster based on the analysis of several thousands of stellar spectra, ranging from the main sequence to the tip of the RGB. Furthermore, two interesting objects were identified; one is a post-AGB star and the other is a possible millisecond-pulsar companion.
Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature. Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays). Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.
Context. The solar magnetic field is responsible for all aspects of solar activity. Thus, emergence of magnetic flux at the surface is the first manifestation of the ensuing solar activity. Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. The small active region NOAA 12118 emerged on 2014 July 17 and was observed one day later with the 1.5-m GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-Perot Interferometer (GFPI) were complemented by synoptic line-of-sight magnetograms and continuum images obtained with the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator (DAVE), respectively. Morphological image processing was employed to measure the photometric and magnetic area, magnetic flux, and the separation profile of the emerging flux region during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days provides a comprehensive view of growth and decay. It traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km s(-1) is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of dark pores. We detected many elongated rapidly expanding granules between the two major polarities, with dimensions twice as large as the normal granules.