530 Physik
Refine
Year of publication
Document Type
- Article (904)
- Doctoral Thesis (373)
- Postprint (123)
- Preprint (50)
- Other (48)
- Habilitation Thesis (23)
- Master's Thesis (10)
- Review (10)
- Monograph/Edited Volume (4)
- Conference Proceeding (3)
Keywords
- diffusion (43)
- anomalous diffusion (33)
- gamma rays: general (20)
- synchronization (19)
- Synchronisation (16)
- organic solar cells (15)
- stochastic processes (15)
- cosmic rays (14)
- ISM: supernova remnants (13)
- data analysis (12)
Institute
- Institut für Physik und Astronomie (1417)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (52)
- Extern (42)
- Institut für Chemie (40)
- Mathematisch-Naturwissenschaftliche Fakultät (25)
- Institut für Mathematik (22)
- Institut für Geowissenschaften (15)
- Institut für Biochemie und Biologie (7)
- Institut für Umweltwissenschaften und Geographie (7)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (6)
The unidentified very-high-energy (VHE; E > 0.1 TeV) gamma -ray source, HESS J1826-130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady gamma -ray flux from HESS J1826-130, which appears extended with a half-width of 0.21 degrees +/- 0.02 <br /> (stat)degrees <br /> stat degrees +/- 0.05 <br /> (sys)degrees sys degrees . The source spectrum is best fit with either a power-law function with a spectral index Gamma = 1.78 +/- 0.10(stat) +/- 0.20(sys) and an exponential cut-off at 15.2 <br /> (+5.5)(-3.2) -3.2+5.5 TeV, or a broken power-law with Gamma (1) = 1.96 +/- 0.06(stat) +/- 0.20(sys), Gamma (2) = 3.59 +/- 0.69(stat) +/- 0.20(sys) for energies below and above E-br = 11.2 +/- 2.7 TeV, respectively. The VHE flux from HESS J1826-130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula, HESS J1825-137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826-130 VHE emission related to PSR J1826-1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826-130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to greater than or similar to 200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants, molecular clouds, and counterparts detected in radio, X-rays, and TeV energies.
A search for new supernova remnants (SNRs) has been conducted using TeV gamma-ray data from the H.E.S.S. Galactic plane survey. As an identification criterion, shell morphologies that are characteristic for known resolved TeV SNRs have been used. Three new SNR candidates were identified in the H.E.S.S. data set with this method. Extensive multiwavelength searches for counterparts were conducted. A radio SNR candidate has been identified to be a counterpart to HESS J1534-571. The TeV source is therefore classified as a SNR. For the other two sources, HESS J1614-518 and HESS J1912 + 101, no identifying counterparts have been found, thus they remain SNR candidates for the time being. TeV-emitting SNRs are key objects in the context of identifying the accelerators of Galactic cosmic rays. The TeV emission of the relativistic particles in the new sources is examined in view of possible leptonic and hadronic emission scenarios, taking the current multiwavelength knowledge into account.
The diffuse very high-energy (VHE; > 100 GeV) gamma-ray emission observed in the central 200 pc of the Milky Way by H.E.S.S. was found to follow dense matter distribution in the central molecular zone (CMZ) up to a longitudinal distance of about 130 pc to the Galactic centre (GC), where the flux rapidly decreases. This was initially interpreted as the result of a burst-like injection of energetic particles 104 yr ago, but a recent more sensitive H.E.S.S. analysis revealed that the cosmic-ray (CR) density profile drops with the distance to the centre, making data compatible with a steady cosmic PeVatron at the GC. In this paper, we extend this analysis to obtain, for the first time, a detailed characterisation of the correlation with matter and to search for additional features and individual gamma-ray sources in the inner 200 pc. Taking advantage of 250 h of H.E.S.S. data and improved analysis techniques, we perform a detailed morphology study of the diffuse VHE emission observed from the GC ridge and reconstruct its total spectrum. To test the various contributions to the total gamma-ray emission, we used an iterative 2D maximum-likelihood approach that allows us to build a phenomenological model of the emission by summing a number of different spatial components. We show that the emission correlated with dense matter covers the full CMZ and that its flux is about half the total diffuse emission flux. We also detect some emission at higher latitude that is likely produced by hadronic collisions of CRs in less dense regions of the GC interstellar medium. We detect an additional emission component centred on the GC and extending over about 15 pc that is consistent with the existence of a strong CR density gradient and confirms the presence of a CR accelerator at the very centre of our Galaxy. We show that the spectrum of full ridge diffuse emission is compatible with that previously derived from the central regions, suggesting that a single population of particles fills the entire CMZ. Finally, we report the discovery of a VHE gamma-ray source near the GC radio arc and argue that it is produced by the pulsar wind nebula candidate G0.13-0.11.
Aims. We study gamma-ray emission from the shell-type supernova remnant (SNR) RXJ0852.0-4622 to better characterize its spectral properties and its distribution over the SNR. Methods. The analysis of an extended High Energy Spectroscopic System (H.E.S.S.) data set at very high energies (E > 100 GeV) permits detailed studies, as well as spatially resolved spectroscopy, of the morphology and spectrum of the whole RXJ0852.0-4622 region. The H.E.S.S. data are combined with archival data from other wavebands and interpreted in the framework of leptonic and hadronic models. The joint Fermi-LAT-H.E.S.S. spectrum allows the direct determination of the spectral characteristics of the parent particle population in leptonic and hadronic scenarios using only GeV-TeV data. Results. An updated analysis of the H.E.S.S. data shows that the spectrum of the entire SNR connects smoothly to the high-energy spectrum measured by Fermi-LAT. The increased data set makes it possible to demonstrate that the H.E.S.S. spectrum deviates significantly from a power law and is well described by both a curved power law and a power law with an exponential cutoff at an energy of E-cut = (6.7 +/- 1.2(stat) +/- 1.2(syst)) TeV. The joint Fermi-LAT-H.E.S.S. spectrum allows the unambiguous identification of the spectral shape as a power law with an exponential cutoff. No significant evidence is found for a variation of the spectral parameters across the SNR, suggesting similar conditions of particle acceleration across the remnant. A simple modeling using one particle population to model the SNR emission demonstrates that both leptonic and hadronic emission scenarios remain plausible. It is also shown that at least a part of the shell emission is likely due to the presence of a pulsar wind nebula around PSR J0855-4644.
Context. Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E >= 0.1 TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the Fermi-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Aims. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Methods. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. Results. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Conclusions. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1% of the kinetic wind energy.
The nine-year H.E.S.S. Galactic Plane Survey (HGPS) has yielded the most uniform observation scan of the inner Milky Way in the TeV gamma-ray band to date. The sky maps and source catalogue of the HGPS allow for a systematic study of the population of TeV pulsar wind nebulae found throughout the last decade. To investigate the nature and evolution of pulsar wind nebulae, for the first time we also present several upper limits for regions around pulsars without a detected TeV wind nebula. Our data exhibit a correlation of TeV surface brightness with pulsar spindown power (E) over dot. This seems to be caused both by an increase of extension with decreasing (E) over dot, and hence with time, compatible with a power law R-PWN((E) over dot) similar to(E) over dot(0.65 +/- 0.20), and by a mild decrease of TeV gamma-ray luminosity with decreasing (E) over dot, compatible with L-1 (10 TeV) similar to (E) over dot(0.59 +/- 0.21). We also find that the off sets of pulsars with respect to the wind nebula centre with ages around 10 kyr are frequently larger than can be plausibly explained by pulsar proper motion and could be due to an asymmetric environment. In the present data, it seems that a large pulsar off set is correlated with a high apparent TeV efficiency L1- 10 TeV / (E) over dot. In addition to 14 HGPS sources considered firmly identified pulsar wind nebulae and 5 additional pulsar wind nebulae taken from literature, we find 10 HGPS sources that are likely TeV pulsar wind nebula candidates. Using a model that subsumes the present common understanding of the very high-energy radiative evolution of pulsar wind nebulae, we find that the trends and variations of the TeV observables and limits can be reproduced to a good level, drawing a consistent picture of present-day TeV data and theory.
The supernova remnant (SNR) W49B originated from a core-collapse supernova that occurred between one and four thousand years ago, and subsequently evolved into a mixed-morphology remnant, which is interacting with molecular clouds (MC). Gamma-ray observations of SNR-MC associations are a powerful tool to constrain the origin of Galactic cosmic rays, as they can probe the acceleration of hadrons through their interaction with the surrounding medium and subsequent emission of non-thermal photons. We report the detection of a gamma-ray source coincident with W49B at very high energies (VHE; E > 100 GeV) with the H.E.S.S. Cherenkov telescopes together with a study of the source with five years of Fermi-LAT high-energy gamma-ray (0.06-300 GeV) data. The smoothly connected, combined source spectrum, measured from 60 MeV to multi-TeV energies, shows two significant spectral breaks at 304 +/- 20 MeV and 8.4(-2.5)(+2.5) GeV; the latter is constrained by the joint fit from the two instruments. The detected spectral features are similar to those observed in several other SNR-MC associations and are found to be indicative of gamma-ray emission produced through neutral-pion decay.
Context. Microquasars are potential gamma-ray emitters. Indications of transient episodes of gamma-ray emission were recently reported in at least two systems: Cyg X-1 and Cyg X-3. The identification of additional gamma-ray-emitting microquasars is required to better understand how gamma-ray emission can be produced in these systems. Aims. Theoretical models have predicted very high-energy (VHE) gamma-ray emission from microquasars during periods of transient outburst. Observations reported herein were undertaken with the objective of observing a broadband flaring event in the gamma-ray and X-ray bands. Methods. Contemporaneous observations of three microquasars, GRS 1915+105, Circinus X-1, and V4641 Sgr, were obtained using the High Energy Spectroscopic System (H.E.S.S.) telescope array and the Rossi X-ray Timing Explorer (RXTE) satellite. X-ray analyses for each microquasar were performed and VHE gamma-ray upper limits from contemporaneous H.E.S.S. observations were derived. Results. No significant gamma-ray signal has been detected in any of the three systems. The integral gamma-ray photon flux at the observational epochs is constrained to be I(>560 GeV) < 7.3 x 10(-13) cm(-2) S-1, I(>560 GeV) < 1.2 x 10-(12) cm s(-1), and I(>240 GeV) < 4.5 x 10(-12) cm(-2) s(-1) for GRS 1915+105, Circinus X-1, and V4641 Sgr, respectively. Conclusions. The gamma-ray upper limits obtained using H.E.S.S. are examined in the context of previous Cherenkov telescope observations of microquasars. The effect of intrinsic absorption is modelled for each target and found to have negligible impact on the flux of escaping gamma-rays. When combined with the X-ray behaviour observed using RXTE, the derived results indicate that if detectable VHE gamma-ray emission from microquasars is commonplace, then it is likely to be highly transient.
Extended VHE gamma-ray emission towards SGR1806-20, LBV 1806-20, and stellar cluster Cl*1806-20
(2018)
Using the High Energy Spectroscopic System (H.E.S.S.) telescopes we have discovered a steady and extended very high-energy (VHE) gamma-ray source towards the luminous blue variable candidate LBV 1806-20, massive stellar cluster Cl* 1806-20, and magnetar SGR 1806-20. The new VHE source, HESS J1808-204, was detected at a statistical significance of >6 sigma (post-trial) with a photon flux normalisation (2.9 +/- 0.4(stat) +/- 0.5(sys)) x 10(-13) ph cm(-2) s(-1) TeV-1 at 1 TeV and a power-law photon index of 2.3 +/- 0.2(stat) +/- 0.3(sys). The luminosity of this source (0.2 to 10 TeV; scaled to distance d = 8 : 7 kpc) is L-VHE similar to 1.6 x 10(34)(d = 8.7 kpc)(2) erg s(-1). The VHE gamma-ray emission is extended and is well fit by a single Gaussian with statistical standard deviation of 0.095 degrees +/- 0.015 degrees. This extension is similar to that of the synchrotron radio nebula G10.0-0.3, which is thought to be powered by LBV 1806-20. The VHE gamma-ray luminosity could be provided by the stellar wind luminosity of LBV 1806-20 by itself and/or the massive star members of Cl* 1806-20. Alternatively, magnetic dissipation (e.g. via reconnection) from SGR 1806-20 can potentially account for the VHE luminosity. The origin and hadronic and/or leptonic nature of the accelerated particles responsible for HESS J1808-204 is not yet clear. If associated with SGR 1806 20, the potentially young age of the magnetar (650 yr) can be used to infer the transport limits of these particles to match the VHE source size. This discovery provides new interest in the potential for high-energy particle acceleration from magnetars, massive stars, and/or stellar clusters.
Shell-type supernova remnants (SNRs) are considered prime candidates for the acceleration of Galactic cosmic rays (CRs) up to the knee of the CR spectrum at E approximate to 3 x 10(15) eV. Our MilkyWay galaxy hosts more than 350 SNRs discovered at radio wavelengths and at high energies, of which 220 fall into the H.E.S.S. Galactic Plane Survey (HGPS) region. Of those, only 50 SNRs are coincident with a H.E.S.S source and in 8 cases the very high-energy (VHE) emission is firmly identified as an SNR. The H.E.S.S. GPS provides us with a legacy for SNR population study in VHE gamma-rays and we use this rich data set to extract VHE flux upper limits from all undetected SNRs. Overall, the derived flux upper limits are not in contradiction with the canonical CR paradigm. Assuming this paradigm holds true, we can constrain typical ambient density values around shell-type SNRs to n <= 7 cm(-3) and electron-to-proton energy fractions above 10 TeV to epsilon(ep) <= 5 x 10(-3). Furthermore, comparisons of VHE with radio luminosities in non-interacting SNRs reveal a behaviour that is in agreement with the theory of magnetic field amplification at shell-type SNRs.
We present the results of the most comprehensive survey of the Galactic plane in very high-energy (VHE) gamma-rays, including a public release of Galactic sky maps, a catalog of VHE sources, and the discovery of 16 new sources of VHE gamma-rays. The High Energy Spectroscopic System (H.E.S.S.) Galactic plane survey (HGPS) was a decade-long observation program carried out by the H.E.S.S. I array of Cherenkov telescopes in Namibia from 2004 to 2013. The observations amount to nearly 2700 h of quality-selected data, covering the Galactic plane at longitudes from l = 250 degrees to 65 degrees and latitudes vertical bar b vertical bar <= 3 degrees. In addition to the unprecedented spatial coverage, the HGPS also features a relatively high angular resolution (0.08 degrees approximate to 5 arcmin mean point spread function 68% containment radius), sensitivity (less than or similar to 1.5% Crab flux for point-like sources), and energy range (0.2-100 TeV). We constructed a catalog of VHE gamma-ray sources from the HGPS data set with a systematic procedure for both source detection and characterization of morphology and spectrum. We present this likelihood-based method in detail, including the introduction of a model component to account for unresolved, large-scale emission along the Galactic plane. In total, the resulting HGPS catalog contains 78 VHE sources, of which 14 are not reanalyzed here, for example, due to their complex morphology, namely shell-like sources and the Galactic center region. Where possible, we provide a firm identification of the VHE source or plausible associations with sources in other astronomical catalogs. We also studied the characteristics of the VHE sources with source parameter distributions. 16 new sources were previously unknown or unpublished, and we individually discuss their identifications or possible associations. We firmly identified 31 sources as pulsar wind nebulae (PWNe), supernova remnants (SNRs), composite SNRs, or gamma-ray binaries. Among the 47 sources not yet identified, most of them (36) have possible associations with cataloged objects, notably PWNe and energetic pulsars that could power VHE PWNe.
The H.E.S.S. Collaboration has discovered a new very high energy (VHE, E > 0.1 TeV) gamma-ray source, HESS J1741-302, located in the Galactic plane. Despite several attempts to constrain its nature, no plausible counterpart has been found so far at X-ray and MeV/GeV gamma-ray energies, and the source remains unidentified. An analysis of 145-h of observations of HESS J1741-302 at VHEs has revealed a steady and relatively weak TeV source (similar to 1% of the Crab Nebula flux), with a spectral index of Gamma = 2.3 +/- 0.2(stat) +/- 0.2(sys), extending to energies up to 10 TeV without any clear signature of a cut-off. In a hadronic scenario, such a spectrum implies an object with particle acceleration up to energies of several hundred TeV. Contrary to most H.E.S.S. unidentified sources, the angular size of HESS J1741-302 is compatible with the H.E.S.S. point spread function at VHEs, with an extension constrained to be below 0.068 degrees at a 99% confidence level. The gamma-ray emission detected by H.E.S.S. can be explained both within a hadronic scenario, due to collisions of protons with energies of hundreds of TeV with dense molecular clouds, and in a leptonic scenario, as a relic pulsar wind nebula, possibly powered by the middle-aged (20 kyr) pulsar PSR B1737-30. A binary scenario, related to the compact radio source 1LC 358.266+0.038 found to be spatially coincident with the best fit position of HESS J1741-302, is also envisaged.
Spectral lines are among the most powerful signatures for dark matter (DM) annihilation searches in very-high-energy gamma rays. The central region of the Milky Way halo is one of the most promising targets given its large amount of DM and proximity to Earth. We report on a search for a monoenergetic spectral line from self-annihilations of DM particles in the energy range from 300 GeV to 70 TeV using a two-dimensional maximum likelihood method taking advantage of both the spectral and spatial features of the signal versus background. The analysis makes use of Galactic center observations accumulated over ten years (2004-2014) with the H.E.S.S. array of ground-based Cherenkov telescopes. No significant gamma-ray excess above the background is found. We derive upper limits on the annihilation cross section (sigma v) for monoenergetic DM lines at the level of 4 x 10(-28) cm(3) s(-1) at 1 TeV, assuming an Einasto DM profile for the Milky Way halo. For a DM mass of 1 TeV, they improve over the previous ones by a factor of 6. The present constraints are the strongest obtained so far for DM particles in the mass range 300 GeV-70 TeV. Ground-based gamma-ray observations have reached sufficient sensitivity to explore relevant velocity-averaged cross sections for DM annihilation into two gamma-ray photons at the level expected from the thermal relic density for TeV DM particles.
Gamma-ray bursts (GRBs) are brief flashes of gamma-rays and are considered to be the most energetic explosive phenomena in the Universe(1). The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow-produced by the interaction between the ejected matter and the circumburst medium-slows down, and a gradual decrease in brightness is observed(2). GRBs typically emit most of their energy via.-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments(3). However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elusive(4). Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow-ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.
The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy (100MeV < E < 100 GeV) gamma-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 target-of-opportunity observations were undertaken with the High Energy Stereoscopic System (H.E.S.S.) in the very-high-energy (VHE, E > 100 GeV) gamma-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with 8 : 7 sigma significance above an energy threshold of 66 GeV. No VHE variability was detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of 4.2 +/- 0.3. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole was estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at r greater than or similar to 1 : 7 X 1017 cm from the black hole, that is beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models were used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data were used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279.
We report on the detection of very high energy (VHE; E > 100 GeV) gamma-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multiwavelength observations with Fermi/LAT, XRT and UVOT onboard the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE gamma-ray regime, we deduce a 95 per cent confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98 and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results, the redshift of KUV 00311-1938 is constrained to 0.51 <= z < 0.98 and of PKS 1440-389 to 0.14 (sic) z < 0.53.
We report on the detection of very high energy (VHE; E > 100 GeV) gamma-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multiwavelength observations with Fermi/LAT, XRT and UVOT onboard the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE gamma-ray regime, we deduce a 95 per cent confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98 and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results, the redshift of KUV 00311-1938 is constrained to 0.51 <= z < 0.98 and of PKS 1440-389 to 0.14 (sic) z < 0.53.
Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95% confidence level of about <sigma upsilon > less than or similar to 3 x 10(-25) cm(3) s(-1) are obtained in the mass range of 400 GeV to 1TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.
Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten super- novae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 to 53 h. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 h. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the >1 TeV gamma-ray flux of the order of similar to 10(-13) cm(-)(2)s(-1) are established, corresponding to upper limits on the luminosities in the range similar to 2 x 10(39) to similar to 1 x 10(42) erg s(-1). These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between similar to 2 x 10(-5) and similar to 2 x 10(-3) M-circle dot yr(-1) under reasonable assumptions on the particle acceleration parameters.
Context. NGC 253 is one of only two starburst galaxies found to emit gamma-rays from hundreds of MeV to multi-TeV energies. Accurate measurements of the very-high-energy (VHE; E> 100 GeV) and high-energy (HE; E > 60 MeV) spectra are crucial to study the underlying particle accelerators, probe the dominant emission mechanism(s) and to study cosmic-ray interaction and transport. Aims. The measurement of the VHE gamma-ray emission of NGC 253 published in 2012 by H.E.S.S. was limited by large systematic uncertainties. Here, the most up to date measurement of the gamma-ray spectrum of NGC 253 is investigated in both HE and VHE gamma-rays. Assuming a hadronic origin of the gamma-ray emission, the measurement uncertainties are propagated into the interpretation of the accelerated particle population. Methods. The data of H.E.S.S. observations are reanalysed using an updated calibration and analysis chain. The improved Fermi-LAT analysis employs more than 8 yr of data processed using pass 8. The cosmic-ray particle population is evaluated from the combined HE-VHE gamma-ray spectrum using NAIMA in the optically thin case. Results. The VHE gamma-ray energy spectrum is best fit by a power-law distribution with a flux normalisation of (1.34 +/- 0.14(stat) +/- 0.27(sys)) x 10(-13) cm(-2) s(-1) TeV-1 at 1 TeV - about 40% above, but compatible with the value obtained in Abramowski et al. (2012). The spectral index Gamma = 2.39 +/- 0.14(stat) +/- 0.25(sys) is slightly softer than but consistent with the previous measurement within systematic errors. In the Fermi energy range an integral flux of F(E > 60 MeV) = (1.56 +/- 0.28(stat) +/- 0.15(sys)) x 10(-8) cm(-2) s(-1) is obtained. At energies above similar to 3 GeV the HE spectrum is consistent with a power-law ranging into the VHE part of the spectrum measured by H.E.S.S. with an overall spectral index Gamma = 2.22 +/- 0.06(stat). Conclusions. Two scenarios for the starburst nucleus are tested, in which the gas in the starburst nucleus acts as either a thin or a thick target for hadronic cosmic rays accelerated by the individual sources in the nucleus. In these two models, the level to which NGC 253 acts as a calorimeter is estimated to a range of f(cal) = 0.1 to 1 while accounting for the measurement uncertainties. The presented spectrum is likely to remain the most accurate measurements until the Cherenkov Telescope Array (CTA) has collected a substantial set of data towards NGC 253.
The blazar Mrk 501 (z = 0.034) was observed at very-high-energy (VHE, E greater than or similar to 100 GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an energy coverage extending significantly up to 20 TeV. This data set is used to constrain Lorentz invariance violation (LIV) using two independent channels: a temporal approach considers the possibility of an energy dependence in the arrival time of gamma-rays, whereas a spectral approach considers the possibility of modifications to the interaction of VHE gamma-rays with extragalactic background light (EBL) photons. The non-detection of energy-dependent time delays and the non-observation of deviations between the measured spectrum and that of a supposed power-law intrinsic spectrum with standard EBL attenuation are used independently to derive strong constraints on the energy scale of LIV (E-QG) in the subluminal scenario for linear and quadratic perturbations in the dispersion relation of photons. For the case of linear perturbations, the 95% confidence level limits obtained are E-QG,E-1 > 3.6 x 10(17) GeV using the temporal approach and E-QG,E-1 > 2.6 x 10(19) GeV using the spectral approach. For the case of quadratic perturbations, the limits obtained are E-QG,E-2 > 8.5 x 10(10) GeV using the temporal approach and E-QG,E-2 > 7.8 x 10(11) GeV using the spectral approach.
PKS 1830-211 is a known macrolensed quasar located at a redshift of z = 2.5. Its highenergy gamma-ray emission has been detected with the Fermi-Large Area Telescope (LAT) instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the High Energy Stereoscopic System (H.E.S.S.) array of Imaging Atmospheric Cherenkov Telescopes in 2014 August, following a flare alert by the Fermi-LAT Collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 d from the alert flare, as expected from observations at other wavelengths. More than 12 h of good-quality data were taken with an analysis threshold of similar to 67 GeV. The significance of a potential signal is computed as a function of the date and the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum.
Context. We present a detailed view of the pulsar wind nebula (PWN) HESS J1825-137. We aim to constrain the mechanisms dominating the particle transport within the nebula, accounting for its anomalously large size and spectral characteristics. Aims. The nebula was studied using a deep exposure from over 12 years of H.E.S.S. I operation, together with data from H.E.S.S. II that improve the low-energy sensitivity. Enhanced energy-dependent morphological and spatially resolved spectral analyses probe the very high energy (VHE, E > 0.1 TeV) gamma-ray properties of the nebula. Methods. The nebula emission is revealed to extend out to 1.5 degrees from the pulsar, similar to 1.5 times farther than previously seen, making HESS J1825-137, with an intrinsic diameter of similar to 100 pc, potentially the largest gamma-ray PWN currently known. Characterising the strongly energy-dependent morphology of the nebula enables us to constrain the particle transport mechanisms. A dependence of the nebula extent with energy of R proportional to E alpha with alpha = -0.29 +/- 0.04(stat) +/- 0.05(sys) disfavours a pure diffusion scenario for particle transport within the nebula. The total gamma-ray flux of the nebula above 1 TeV is found to be (1.12 +/- 0.03(stat) +/- 0.25(sys)) +/- 10(-11) cm(-2) s(-1), corresponding to similar to 64% of the flux of the Crab nebula. Results. HESS J1825-137 is a PWN with clearly energy-dependent morphology at VHE gamma-ray energies. This source is used as a laboratory to investigate particle transport within intermediate-age PWNe. Based on deep observations of this highly spatially extended PWN, we produce a spectral map of the region that provides insights into the spectral variation within the nebula.
Context. Pulsar wind nebulae (PWNe) represent the most prominent population of Galactic very-high-energy gamma-ray sources and are thought to be an efficient source of leptonic cosmic rays. Vela X is a nearby middle-aged PWN, which shows bright X-ray and TeV gamma-ray emission towards an elongated structure called the cocoon. Aims. Since TeV emission is likely inverse-Compton emission of electrons, predominantly from interactions with the cosmic microwave background, while X-ray emission is synchrotron radiation of the same electrons, we aim to derive the properties of the relativistic particles and of magnetic fields with minimal modelling. Methods. We used data from the Suzaku XIS to derive the spectra from three compact regions in Vela X covering distances from 0.3 to 4 pc from the pulsar along the cocoon. We obtained gamma-ray spectra of the same regions from H.E.S.S. observations and fitted a radiative model to the multi-wavelength spectra. Results. The TeV electron spectra and magnetic field strengths are consistent within the uncertainties for the three regions, with energy densities of the order 10(-12) erg cm(-3). The data indicate the presence of a cutoff in the electron spectrum at energies of similar to 100 TeV and a magnetic field strength of similar to 6 mu G. Constraints on the presence of turbulent magnetic fields are weak. Conclusions. The pressure of TeV electrons and magnetic fields in the cocoon is dynamically negligible, requiring the presence of another dominant pressure component to balance the pulsar wind at the termination shock. Sub-TeV electrons cannot completely account for the missing pressure, which may be provided either by relativistic ions or from mixing of the ejecta with the pulsar wind. The electron spectra are consistent with expectations from transport scenarios dominated either by advection via the reverse shock or by diffusion, but for the latter the role of radiative losses near the termination shock needs to be further investigated in the light of the measured cutoff energies. Constraints on turbulent magnetic fields and the shape of the electron cutoff can be improved by spectral measurements in the energy range greater than or similar to 10 keV.
Context. Recently, the high-energy (HE, 0.1-100 GeV) gamma-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic gamma-ray binary. Aims. This work aims at the detection of very-high-energy (VHE, >100 GeV) gamma-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Results. VHE gamma-ray emission is detected with a statistical significance of 6.4 sigma. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the 1-10 TeV energy range is (1.4 +/- 0.2) x 10(35) erg s(-1). A luminosity of (5 +/- 1) x 10(35) erg s(-1) is reached during 20% of the orbit. HE and VHE gamma-ray emissions are anti-correlated. LMC P3 is the most luminous gamma-ray binary known so far.
One of the major challenges in chemical synthesis is to trigger and control a specific reaction route leading to a specific final product, while side products are avoided. Methodologies based on resonant processes at the molecular level, for example, photochemistry, offer the possibility of inducing selective reactions. Electrons at energies below the molecular ionization potential (<10 eV) are known to dissociate molecules via resonant processes with higher cross sections and specificity than photons. Here we show that even subexcitation electrons with energies as low as 1 eV produce ethylene and acetylene from dimethyl sulfide in competing reactions. However, the production of ethylene can specifically be targeted by controlling the energy of electrons (similar to 3 to 4 eV). Finally, pure ethylene is selectively desorbed by heating the substrate from 90 to 105 K. Beyond the synthesis of these versatile hydrocarbons for various industrial applications from a biogenic sulfur compound, our findings demonstrate the feasibility of electron induced selective chemistry applicable on the nanoscale.
B2 1215+30 is a BL-Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes and subsequently confirmed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2. 1215+30 during routine monitoring observations of the blazar 1ES. 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of <3.6 hr. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla Observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a onezone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor delta > 10, and an electron population with spectral index p < 2.3.
We present results from deep observations toward the Cygnus region using 300 hr of very high energy (VHE)gamma-ray data taken with the VERITAS Cerenkov telescope array and over 7 yr of high-energy.-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse gamma-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130) and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL J2021.0+4031e (the Gamma Cygni supernova remnant) was examined, and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE gamma-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula flux. We interpret these observations in a multiwavelength context and present the most detailed gamma-ray view of the region to date.
On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event IC 170922A, of potential astrophysical origin. It was soon determined that the neutrino direction was consistent with the location of the gamma-ray blazar TXS 0506+056. (3FGL J0509.4+ 0541), which was in an elevated gamma-ray emission state as measured by the Fermi satellite. Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E > 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+ 056 was detected by VERITAS with a significance of 5.8 standard deviations (sigma) in the full 35 hr data set. The average photon flux of the source during this period was (8.9 +/- 1.6). x. 10(-12) cm(-2) s(-1), or 1.6% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of 4.8. +/-. 1.3.
VERITAS and Fermi-LAT Observations of TeV Gamma-Ray Sources Discovered by HAWC in the 2HWC Catalog
(2018)
The High Altitude Water Cherenkov (HAWC) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100 GeV) gamma-ray sources based on 507 days of observation. Among these, 19 sources are not associated with previously known teraelectronvolt (TeV) gamma-ray sources. We have studied 14 of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1 TeV-30 TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected 14 new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected gigaelectronvolt (GeV) gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC, and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.
On 2015 March 23, the Very Energetic Radiation Imaging Telescope Array System (VERITAS) responded to a Swift-Burst Alert Telescope (BAT) detection of a gamma-ray burst, with observations beginning 270 s after the onset of BAT emission, and only 135 s after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40-minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant because the very-high-energy (VHE) observation started only similar to 2 minutes after the prompt emission peaked, and Fermi-Large Area Telescope observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB 150323A (z = 0.593) limits the attenuation by the extragalactic background light to similar to 50% at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below similar to 100 GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be A greater than or similar to 3 x 10(11) g . cm(-1), consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the interstellar medium, which therefore cannot be ruled out as the environment of GRB 150323A.
Combined with measurements made by very-long-baseline interferometry, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL. Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL. Lacertae observed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS). with a rise time of similar to 2.3 hr and a decay time of similar to 36 min. The peak flux above 200 GeV is (4.2 +/- 0.6) x 10(-6) photon m(-2) s(-1) measured with a 4-minute-binned light curve, corresponding to similar to 180% of the flux that is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in Very Long Baseline Array observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models that invoke relativistic plasma passing stationary shocks.
Aims. We aim for an understanding of the morphological and spectral properties of the supernova remnant RCW 86 and for insights into the production mechanism leading to the RCW 86 very high-energy gamma-ray emission. Methods. We analyzed High Energy Spectroscopic System (H.E.S.S.) data that had increased sensitivity compared to the observations presented in the RCW 86 H.E.S.S. discovery publication. Studies of the morphological correlation between the 0.5-1 keV X-ray band, the 2-5 keV X-ray band, radio, and gamma-ray emissions have been performed as well as broadband modeling of the spectral energy distribution with two different emission models. Results. We present the first conclusive evidence that the TeV gamma-ray emission region is shell-like based on our morphological studies. The comparison with 2-5 keV X-ray data reveals a correlation with the 0.4-50 TeV gamma-ray emission. The spectrum of RCW 86 is best described by a power law with an exponential cutoff at E-cut = (3.5 +/- 1.2(stat)) TeV and a spectral index of Gamma approximate to 1.6 +/- 0.2. A static leptonic one-zone model adequately describes the measured spectral energy distribution of RCW 86, with the resultant total kinetic energy of the electrons above 1 GeV being equivalent to similar to 0.1% of the initial kinetic energy of a Type Ia supernova explosion (10(51) erg). When using a hadronic model, a magnetic field of B approximate to 100 mu G is needed to represent the measured data. Although this is comparable to formerly published estimates, a standard E-2 spectrum for the proton distribution cannot describe the gamma-ray data. Instead, a spectral index of Gamma(p) approximate to 1.7 would be required, which implies that similar to 7 x 10(49)/n(cm-3) erg has been transferred into high-energy protons with the effective density n(cm-3) = n/1 cm(-3). This is about 10% of the kinetic energy of a typical Type Ia supernova under the assumption of a density of 1 cm(-3).
Prospects for Cherenkov Telescope Array Observations of the Young Supernova Remnant RX J1713.7-3946
(2017)
We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX J1713.7-3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very high energy (VHE) gamma rays. Special attention is paid to exploring possible spatial (anti) correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H I emission. We present a series of simulated images of RX J1713.7-3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the nonthermal X-ray emission observed by XMM-Newton, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H I observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic versus leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.
Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
In recent years, complex network analysis facilitated the identification of universal and unexpected patterns in complex climate systems. However, the analysis and representation of a multiscale complex relationship that exists in the global climate system are limited. A logical first step in addressing this issue is to construct multiple networks over different timescales. Therefore, we propose to apply the wavelet multiscale correlation (WMC) similarity measure, which is a combination of two state-of-the-art methods, viz. wavelet and Pearson’s correlation, for investigating multiscale processes through complex networks. Firstly we decompose the data over different timescales using the wavelet approach and subsequently construct a corresponding network by Pearson’s correlation. The proposed approach is illustrated and tested on two synthetics and one real-world example. The first synthetic case study shows the efficacy of the proposed approach to unravel scale-specific connections, which are often undiscovered at a single scale. The second synthetic case study illustrates that by dividing and constructing a separate network for each time window we can detect significant changes in the signal structure. The real-world example investigates the behavior of the global sea surface temperature (SST) network at different timescales. Intriguingly, we notice that spatial dependent structure in SST evolves temporally. Overall, the proposed measure has an immense potential to provide essential insights on understanding and extending complex multivariate process studies at multiple scales.
Subject of this work is the investigation of universal scaling laws which are observed in coupled chaotic systems. Progress is made by replacing the chaotic fluctuations in the perturbation dynamics by stochastic processes. First, a continuous-time stochastic model for weakly coupled chaotic systems is introduced to study the scaling of the Lyapunov exponents with the coupling strength (coupling sensitivity of chaos). By means of the the Fokker-Planck equation scaling relations are derived, which are confirmed by results of numerical simulations. Next, the new effect of avoided crossing of Lyapunov exponents of weakly coupled disordered chaotic systems is described, which is qualitatively similar to the energy level repulsion in quantum systems. Using the scaling relations obtained for the coupling sensitivity of chaos, an asymptotic expression for the distribution function of small spacings between Lyapunov exponents is derived and compared with results of numerical simulations. Finally, the synchronization transition in strongly coupled spatially extended chaotic systems is shown to resemble a continuous phase transition, with the coupling strength and the synchronization error as control and order parameter, respectively. Using results of numerical simulations and theoretical considerations in terms of a multiplicative noise partial differential equation, the universality classes of the observed two types of transition are determined (Kardar-Parisi-Zhang equation with saturating term, directed percolation).
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.
In the present work, we study wave phenomena in strongly nonlinear lattices. Such lattices are characterized by the absence of classical linear waves. We demonstrate that compactons – strongly localized solitary waves with tails decaying faster than exponential – exist and that they play a major role in the dynamics of the system under consideration. We investigate compactons in different physical setups. One part deals with lattices of dispersively coupled limit cycle oscillators which find various applications in natural sciences such as Josephson junction arrays or coupled Ginzburg-Landau equations. Another part deals with Hamiltonian lattices. Here, a prominent example in which compactons can be found is the granular chain. In the third part, we study systems which are related to the discrete nonlinear Schrödinger equation describing, for example, coupled optical wave-guides or the dynamics of Bose-Einstein condensates in optical lattices. Our investigations are based on a numerical method to solve the traveling wave equation. This results in a quasi-exact solution (up to numerical errors) which is the compacton. Another ansatz which is employed throughout this work is the quasi-continuous approximation where the lattice is described by a continuous medium. Here, compactons are found analytically, but they are defined on a truly compact support. Remarkably, both ways give similar qualitative and quantitative results. Additionally, we study the dynamical properties of compactons by means of numerical simulation of the lattice equations. Especially, we concentrate on their emergence from physically realizable initial conditions as well as on their stability due to collisions. We show that the collisions are not exactly elastic but that a small part of the energy remains at the location of the collision. In finite lattices, this remaining part will then trigger a multiple scattering process resulting in a chaotic state.
Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed.
New amphiphilic compounds 1-9 that feature a construction with dendronized hydrophilic and hydrophobic segment groups connected to a specific aromatic or aliphatic spacer unit have been synthesized, following a modular building block strategy. The hydrophilic dendrons are typically branched elements with peripheral carboxylic groups, unlike the hydrophobic dendrons that contain peripheral alkyl chains as part of respective amide functions. The hydrophilic dendrons are in different generations of branching, while the hydrophobic dendrons are all in the first generation of branching (three terminal branching), but differ in the length of the alkyl chains, thus giving rise to designed structure and amphiphilic properties in the new compounds. The resulting surfactants are capable of forming well-defined Langmuir films of remarkable stability when spread from a solution onto an aqueous subphase. Nevertheless, specific packing behaviour and orientation of the amphiphilic molecules were found, depending on the molecular structure, as determined using analysis of the surface pressure-area (pi-A) isotherms. Langmuir-Blodgett transfer of the first monolayer from a pure water subphase to a clean silicon wafer proved possible for the amphiphiles of peripheral alkyl chain length C-12, while the amphiphiles with the longer alkyl chains failed, possibly due to the more rigid monolayers they form, impeding the transfer.
Suppression of the TeV Pair-beam-Plasma Instability by a Tangled Weak Intergalactic Magnetic Field
(2022)
We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kiloparsec, we find that weak magnetic fields increase the transverse momentum of the pair-beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission.
Since their discovery in 1610 by Galileo Galilei, Saturn's rings continue to fascinate both experts and amateurs. Countless numbers of icy grains in almost Keplerian orbits reveal a wealth of structures such as ringlets, voids and gaps, wakes and waves, and many more. Grains are found to increase in size with increasing radial distance to Saturn. Recently discovered "propeller" structures in the Cassini spacecraft data, provide evidence for the existence of embedded moonlets. In the wake of these findings, the discussion resumes about origin and evolution of planetary rings, and growth processes in tidal environments. In this thesis, a contact model for binary adhesive, viscoelastic collisions is developed that accounts for agglomeration as well as restitution. Collisional outcomes are crucially determined by the impact speed and masses of the collision partners and yield a maximal impact velocity at which agglomeration still occurs. Based on the latter, a self-consistent kinetic concept is proposed. The model considers all possible collisional outcomes as there are coagulation, restitution, and fragmentation. Emphasizing the evolution of the mass spectrum and furthermore concentrating on coagulation alone, a coagulation equation, including a restricted sticking probability is derived. The otherwise phenomenological Smoluchowski equation is reproduced from basic principles and denotes a limit case to the derived coagulation equation. Qualitative and quantitative analysis of the relevance of adhesion to force-free granular gases and to those under the influence of Keplerian shear is investigated. Capture probability, agglomerate stability, and the mass spectrum evolution are investigated in the context of adhesive interactions. A size dependent radial limit distance from the central planet is obtained refining the Roche criterion. Furthermore, capture probability in the presence of adhesion is generally different compared to the case of pure gravitational capture. In contrast to a Smoluchowski-type evolution of the mass spectrum, numerical simulations of the obtained coagulation equation revealed, that a transition from smaller grains to larger bodies cannot occur via a collisional cascade alone. For parameters used in this study, effective growth ceases at an average size of centimeters.
A dramatic efficiency improvement of bulk heterojunction solar cells based on electron-donating conjugated polymers in combination with soluble fullerene derivatives has been achieved over the past years. Certified and reported power conversion efficiencies now reach over 9% for single junctions and exceed the 10% benchmark for tandem solar cells. This trend brightens the vision of organic photovoltaics becoming competitive with inorganic solar cells including the realization of low-cost and large-area organic photovoltaics. For the best performing organic materials systems, the yield of charge generation can be very efficient. However, a detailed understanding of the free charge carrier generation mechanisms at the donor acceptor interface and the energy loss associated with it needs to be established. Moreover, organic solar cells are limited by the competition between charge extraction and free charge recombination, accounting for further efficiency losses. A conclusive picture and the development of precise methodologies for investigating the fundamental processes in organic solar cells are crucial for future material design, efficiency optimization, and the implementation of organic solar cells into commercial products.
In order to advance the development of organic photovoltaics, my thesis focuses on the comprehensive understanding of charge generation, recombination and extraction in organic bulk heterojunction solar cells summarized in 6 chapters on the cumulative basis of 7 individual publications.
The general motivation guiding this work was the realization of an efficient hybrid inorganic/organic tandem solar cell with sub-cells made from amorphous hydrogenated silicon and organic bulk heterojunctions. To realize this project aim, the focus was directed to the low band-gap copolymer PCPDTBT and its derivatives, resulting in the examination of the charge carrier dynamics in PCPDTBT:PC70BM blends in relation to by the blend morphology. The phase separation in this blend can be controlled by the processing additive diiodooctane, enhancing domain purity and size. The quantitative investigation of the free charge formation was realized by utilizing and improving the time delayed collection field technique. Interestingly, a pronounced field dependence of the free carrier generation for all blends is found, with the field dependence being stronger without the additive. Also, the bimolecular recombination coefficient for both blends is rather high and increases with decreasing internal field which we suggest to be caused by a negative field dependence of mobility. The additive speeds up charge extraction which is rationalized by the threefold increase in mobility.
By fluorine attachment within the electron deficient subunit of PCPDTBT, a new polymer F-PCPDTBT is designed. This new material is characterized by a stronger tendency to aggregate as compared to non-fluorinated PCPDTBT. Our measurements show that for F-PCPDTBT:PCBM blends the charge carrier generation becomes more efficient and the field-dependence of free charge carrier generation is weakened. The stronger tendency to aggregate induced by the fluorination also leads to increased polymer rich domains, accompanied in a threefold reduction in the non-geminate recombination coefficient at conditions of open circuit. The size of the polymer domains is nicely correlated to the field-dependence of charge generation and the Langevin reduction factor, which highlights the importance of the domain size and domain purity for efficient charge carrier generation. In total, fluorination of PCPDTBT causes the PCE to increase from 3.6 to 6.1% due to enhanced fill factor, short circuit current and open circuit voltage. Further optimization of the blend ratio, active layer thickness, and polymer molecular weight resulted in 6.6% efficiency for F-PCPDTBT:PC70BM solar cells.
Interestingly, the double fluorinated version 2F-PCPDTBT exhibited poorer FF despite a further reduction of geminate and non-geminate recombination losses. To further analyze this finding, a new technique is developed that measures the effective extraction mobility under charge carrier densities and electrical fields comparable to solar cell operation conditions. This method involves the bias enhanced charge extraction technique. With the knowledge of the carrier density under different electrical field and illumination conditions, a conclusive picture of the changes in charge carrier dynamics leading to differences in the fill factor upon fluorination of PCPDTBT is attained. The more efficient charge generation and reduced recombination with fluorination is counterbalanced by a decreased extraction mobility. Thus, the highest fill factor of 60% and efficiency of 6.6% is reached for F-PCPDTBT blends, while 2F-PCPDTBT blends have only moderate fill factors of 54% caused by the lower effective extraction mobility, limiting the efficiency to 6.5%.
To understand the details of the charge generation mechanism and the related losses, we evaluated the yield and field-dependence of free charge generation using time delayed collection field in combination with sensitive measurements of the external quantum efficiency and absorption coefficients for a variety of blends. Importantly, both the yield and field-dependence of free charge generation is found to be unaffected by excitation energy, including direct charge transfer excitation below the optical band gap. To access the non-detectable absorption at energies of the relaxed charge transfer emission, the absorption was reconstructed from the CT emission, induced via the recombination of thermalized charges in electroluminescence. For a variety of blends, the quantum yield at energies of charge transfer emission was identical to excitations with energies well above the optical band-gap. Thus, the generation proceeds via the split-up of the thermalized charge transfer states in working solar cells. Further measurements were conducted on blends with fine-tuned energy levels and similar blend morphologies by using different fullerene derivatives. A direct correlation between the efficiency of free carrier generation and the energy difference of the relaxed charge transfer state relative to the energy of the charge separated state is found. These findings open up new guidelines for future material design as new high efficiency materials require a minimum energetic offset between charge transfer and the charge separated state while keeping the HOMO level (and LUMO level) difference between donor and acceptor as small as possible.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
The theory of atomic Boson-Fermion mixtures in the dilute limit beyond mean-field is considered in this thesis. Extending the formalism of quantum field theory we derived expressions for the quasi-particle excitation spectra, the ground state energy, and related quantities for a homogenous system to first order in the dilute gas parameter. In the framework of density functional theory we could carry over the previous results to inhomogeneous systems. We then determined to density distributions for various parameter values and identified three different phase regions: (i) a stable mixed regime, (ii) a phase separated regime, and (iii) a collapsed regime. We found a significant contribution of exchange-correlation effects in the latter case. Next, we determined the shift of the Bose-Einstein condensation temperature caused by Boson-Fermion interactions in a harmonic trap due to redistribution of the density profiles. We then considered Boson-Fermion mixtures in optical lattices. We calculated the criterion for stability against phase separation, identified the Mott-insulating and superfluid regimes both, analytically within a mean-field calculation, and numerically by virtue of a Gutzwiller Ansatz. We also found new frustrated ground states in the limit of very strong lattices. ----Anmerkung: Der Autor ist Träger des durch die Physikalische Gesellschaft zu Berlin vergebenen Carl-Ramsauer-Preises 2004 für die jeweils beste Dissertation der vier Universitäten Freie Universität Berlin, Humboldt-Universität zu Berlin, Technische Universität Berlin und Universität Potsdam.
Whereas stable homogenous states of aqueous hydrocarbon solutions are typically observed at high temperatures and pressures far beyond the critical values corresponding to individual components, the stability of such system may be preserved upon transition into the region of metastable states at low temperatures and low pressures. This work is dedicated to the study of the thermal stability of a water-methane mixture formed by cryogenic vapor phase deposition. The obtained thin films were studied using vibrational spectroscopy in the temperature range of 16-180 K. During thermal annealing of the samples, characteristic vibrational C-H modes of methane were monitored alongside the chamber pressure to register both structural changes and desorption of the film material. The obtained results reveal that upon the co-deposition of methane and water, methane molecules appear both in non-bound and trapped states. The observed broadening of the characteristic C-H stretching mode at 3010 cm(-1) upon an increase in temperature of the sample from 16 to 90 K, followed by narrowing of the peak as the temperature is reduced back to 16 K, indicates localization of methane molecules within the water matrix at lower temperatures.
Bacterial chemotaxis-a fundamental example of directional navigation in the living world-is key to many biological processes, including the spreading of bacterial infections. Many bacterial species were recently reported to exhibit several distinct swimming modes-the flagella may, for example, push the cell body or wrap around it. How do the different run modes shape the chemotaxis strategy of a multimode swimmer? Here, we investigate chemotactic motion of the soil bacterium Pseudomonas putida as a model organism. By simultaneously tracking the position of the cell body and the configuration of its flagella, we demonstrate that individual run modes show different chemotactic responses in nutrition gradients and, thus, constitute distinct behavioral states. On the basis of an active particle model, we demonstrate that switching between multiple run states that differ in their speed and responsiveness provides the basis for robust and efficient chemotaxis in complex natural habitats.
Bacterial chemotaxis-a fundamental example of directional navigation in the living world-is key to many biological processes, including the spreading of bacterial infections. Many bacterial species were recently reported to exhibit several distinct swimming modes-the flagella may, for example, push the cell body or wrap around it. How do the different run modes shape the chemotaxis strategy of a multimode swimmer? Here, we investigate chemotactic motion of the soil bacterium Pseudomonas putida as a model organism. By simultaneously tracking the position of the cell body and the configuration of its flagella, we demonstrate that individual run modes show different chemotactic responses in nutrition gradients and, thus, constitute distinct behavioral states. On the basis of an active particle model, we demonstrate that switching between multiple run states that differ in their speed and responsiveness provides the basis for robust and efficient chemotaxis in complex natural habitats.
The topic of synchronization forms a link between nonlinear dynamics and neuroscience. On the one hand, neurobiological research has shown that the synchronization of neuronal activity is an essential aspect of the working principle of the brain. On the other hand, recent advances in the physical theory have led to the discovery of the phenomenon of phase synchronization. A method of data analysis that is motivated by this finding - phase synchronization analysis - has already been successfully applied to empirical data. The present doctoral thesis ties up to these converging lines of research. Its subject are methodical contributions to the further development of phase synchronization analysis, as well as its application to event-related potentials, a form of EEG data that is especially important in the cognitive sciences. The methodical contributions of this work consist firstly in a number of specialized statistical tests for a difference in the synchronization strength in two different states of a system of two oscillators. Secondly, in regard of the many-channel character of EEG data an approach to multivariate phase synchronization analysis is presented. For the empirical investigation of neuronal synchronization a classic experiment on language processing was replicated, comparing the effect of a semantic violation in a sentence context with that of the manipulation of physical stimulus properties (font color). Here phase synchronization analysis detects a decrease of global synchronization for the semantic violation as well as an increase for the physical manipulation. In the latter case, by means of the multivariate analysis the global synchronization effect can be traced back to an interaction of symmetrically located brain areas.<BR> The findings presented show that the method of phase synchronization analysis motivated by physics is able to provide a relevant contribution to the investigation of event-related potentials in the cognitive sciences.
In Allefeld & Kurths [2004], we introduced an approach to multivariate phase synchronization analysis in the form of a Synchronization Cluster Analysis (SCA). A statistical model of a synchronization cluster was described, and an abbreviated instruction on how to apply this model to empirical data was given, while an implementation of the corresponding algorithm was (and is) available from the authors. In this letter, the complete details on how the data analysis algorithm is to be derived from the model are filled in.
Phase synchronization analysis, including our recently introduced multivariate approach, is applied to event-related EEG data from an experiment on language processing, following a classic psycholinguistic paradigm. For the two types of experimental manipulation distinct effects in overall synchronization are found; for one of them they can also be localized. The synchronization effects occur earlier than those found by the conventional analysis method, indicating that the new approach provides additional information on the underlying neuronal process.
In order to investigate the temporal characteristics of cognitive processing, we apply multivariate phase synchronization analysis to event-related potentials. The experimental design combines a semantic incongruity in a sentence context with a physical mismatch (color change). In the ERP average, these result in an N400 component and a P300-like positivity, respectively. The synchronization analysis shows an effect of global desynchronization in the theta band around 288ms after stimulus presentation for the semantic incongruity, while the physical mismatch elicits an increase of global synchronization in the alpha band around 204ms. Both of these effects clearly precede those in the ERP average. Moreover, the delay between synchronization effect and ERP component correlates with the complexity of the cognitive processes.
We present different tests for phase synchronization which improve the procedures currently used in the literature. This is accomplished by using a two-samples test setup and by utilizing insights and methods from directional statistics and bootstrap theory. The tests differ in the generality of the situation in which they can be applied as well as in their complexity, including computational cost. A modification of the resampling technique of the bootstrap is introduced, making it possible to fully utilize data from time series.
A method for the multivariate analysis of statistical phase synchronization phenomena in empirical data is presented. A first statistical approach is complemented by a stochastic dynamic model, to result in a data analysis algorithm which can in a specific sense be shown to be a generic multivariate statistical phase synchronization analysis. The method is applied to EEG data from a psychological experiment, obtaining results which indicate the relevance of this method in the context of cognitive science as well as in other fields.
Context:
Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way.
Aims:
The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region.
Methods:
In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries.
Results:
Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus’s core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2–O7) tend to have shorter orbital periods than later spectral types (O9.2–O9.7).
Conclusions:
Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe’s peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙.
The interplay between nanomorphology and efficiency of polymer-fullerene bulk-heterojunction (BHJ) solar cells has been the subject of intense research, but the generality of these concepts for small-molecule (SM) BHJs remains unclear. Here, the relation between performance; charge generation, recombination, and extraction dynamics; and nanomorphology achievable with two SM donors benzo[1,2-b:4,5-b]dithiophene-pyrido[3,4-b]-pyrazine BDT(PPTh2)(2), namely SM1 and SM2, differing by their side-chains, are examined as a function of solution additive composition. The results show that the additive 1,8-diiodooctane acts as a plasticizer in the blends, increases domain size, and promotes ordering/crystallinity. Surprisingly, the system with high domain purity (SM1) exhibits both poor exciton harvesting and severe charge trapping, alleviated only slightly with increased crystallinity. In contrast, the system consisting of mixed domains and lower crystallinity (SM2) shows both excellent exciton harvesting and low charge recombination losses. Importantly, the onset of large, pure crystallites in the latter (SM2) system reduces efficiency, pointing to possible differences in the ideal morphologies for SM-based BHJ solar cells compared with polymer-fullerene devices. In polymer-based systems, tie chains between pure polymer crystals establish a continuous charge transport network, whereas SM-based active layers may in some cases require mixed domains that enable both aggregation and charge percolation to the electrodes.
Spin precession in magnetic materials is commonly modelled with the classical phenomenological Landau-Lifshitz-Gilbert (LLG) equation. Based on a quantized three-dimensional spin + environment Hamiltonian, we here derive a spin operator equation of motion that describes precession and includes a general form of damping that consistently accounts for memory, coloured noise and quantum statistics. The LLG equation is recovered as its classical, Ohmic approximation. We further introduce resonant Lorentzian system-reservoir couplings that allow a systematic comparison of dynamics between Ohmic and non-Ohmic regimes. Finally, we simulate the full non-Markovian dynamics of a spin in the semi-classical limit. At low temperatures, our numerical results demonstrate a characteristic reduction and flattening of the steady state spin alignment with an external field, caused by the quantum statistics of the environment. The results provide a powerful framework to explore general three-dimensional dissipation in quantum thermodynamics.
Development of functional and stable solid polymer electrolytes (SPEs) for battery applications is an important step towards both safer batteries and for the realization of lithium-based or anode-less batteries. The interface between the lithium and the solid polymer electrolyte is one of the bottlenecks, where severe degradation is expected. Here, the stability of three different SPEs - poly(ethylene oxide) (PEO), poly(epsilon-caprolactone) (PCL) and poly(trimethylene carbonate) (PTMC) - together with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt, is investigated after they have been exposed to lithium metal under UHV conditions. Degradation compounds, e.g. Li-O-R, LiF and LixSyOz, are identified for all SPEs using soft X-ray photoelectron spectroscopy. A competing degradation between polymer and salt is identified in the outermost surface region (<7 nm), and is dependent on the polymer host. PTMC:LiTFSI shows the most severe decomposition of both polymer and salt followed by PCL:LiTFSI and PEO:LiTFSI. In addition, the movement of lithium species through the decomposed interface shows large variation depending on the polymer electrolyte system.
HESS J1826-130
(2017)
HESS J1826-130 is an unidentified hard spectrum source discovered by H.E.S.S. along the Galactic plane, the spectral index being Gamma = 1.6 with an exponential cut-off at about 12 TeV. While the source does not have a clear counterpart at longer wavelengths, the very hard spectrum emission at TeV energies implies that electrons or protons accelerated up to several hundreds of TeV are responsible for the emission. In the hadronic case, the VHE emission can be produced by runaway cosmic-rays colliding with the dense molecular clouds spatially coincident with the H.E.S.S. source.
Chemotactic motion in a chemical gradient is an essential cellular function that controls many processes in the living world. For a better understanding and more detailed modelling of the underlying mechanisms of chemotaxis, quantitative investigations in controlled environments are needed. We developed a setup that allows us to separately address the dependencies of the chemotactic motion on the average background concentration and on the gradient steepness of the chemoattractant. In particular, both the background concentration and the gradient steepness can be kept constant at the position of the cell while it moves along in the gradient direction. This is achieved by generating a well-defined chemoattractant gradient using flow photolysis. In this approach, the chemoattractant is released by a light-induced reaction from a caged precursor in a microfluidic flow chamber upstream of the cell. The flow photolysis approach is combined with an automated real-time cell tracker that determines changes in the cell position and triggers movement of the microscope stage such that the cell motion is compensated and the cell remains at the same position in the gradient profile. The gradient profile can be either determined experimentally using a caged fluorescent dye or may be alternatively determined by numerical solutions of the corresponding physical model. To demonstrate the function of this adaptive microfluidic gradient generator, we compare the chemotactic motion of Dictyostelium discoideum cells in a static gradient and in a gradient that adapts to the position of the moving cell. Published by AIP Publishing.
Gravitational-wave (GW) astrophysics is a field in full blossom. Since the landmark detection of GWs from a binary black hole on September 14th 2015, fifty-two compact-object binaries have been reported by the LIGO-Virgo collaboration. Such events carry astrophysical and cosmological information ranging from an understanding of how black holes and neutron stars are formed, what neutron stars are composed of, how the Universe expands, and allow testing general relativity in the highly-dynamical strong-field regime. It is the goal of GW astrophysics to extract such information as accurately as possible. Yet, this is only possible if the tools and technology used to detect and analyze GWs are advanced enough. A key aspect of GW searches are waveform models, which encapsulate our best predictions for the gravitational radiation under a certain set of parameters, and that need to be cross-correlated with data to extract GW signals. Waveforms must be very accurate to avoid missing important physics in the data, which might be the key to answer the fundamental questions of GW astrophysics. The continuous improvements of the current LIGO-Virgo detectors, the development of next-generation ground-based detectors such as the Einstein Telescope or the Cosmic Explorer, as well as the development of the Laser Interferometer Space Antenna (LISA), demand accurate waveform models. While available models are enough to capture the low spins, comparable-mass binaries routinely detected in LIGO-Virgo searches, those for sources from both current and next-generation ground-based and spaceborne detectors must be accurate enough to detect binaries with large spins and asymmetry in the masses. Moreover, the thousands of sources that we expect to detect with future detectors demand accurate waveforms to mitigate biases in the estimation of signals’ parameters due to the presence of a foreground of many sources that overlap in the frequency band. This is recognized as one of the biggest challenges for the analysis of future-detectors’ data, since biases might hinder the extraction of important astrophysical and cosmological information from future detectors’ data. In the first part of this thesis, we discuss how to improve waveform models for binaries with high spins and asymmetry in the masses. In the second, we present the first generic metrics that have been proposed to predict biases in the presence of a foreground of many overlapping signals in GW data.
For the first task, we will focus on several classes of analytical techniques. Current models for LIGO and Virgo studies are based on the post-Newtonian (PN, weak-field, small velocities) approximation that is most natural for the bound orbits that are routinely detected in GW searches. However, two other approximations have risen in prominence, the post-Minkowskian (PM, weak- field only) approximation natural for unbound (scattering) orbits and the small-mass-ratio (SMR) approximation typical of binaries in which the mass of one body is much bigger than the other. These are most appropriate to binaries with high asymmetry in the masses that challenge current waveform models. Moreover, they allow one to “cover” regions of the parameter space of coalescing binaries, thereby improving the interpolation (and faithfulness) of waveform models. The analytical approximations to the relativistic two-body problem can synergically be included within the effective-one-body (EOB) formalism, in which the two-body information from each approximation can be recast into an effective problem of a mass orbiting a deformed Schwarzschild (or Kerr) black hole. The hope is that the resultant models can cover both the low-spin comparable-mass binaries that are routinely detected, and the ones that challenge current models. The first part of this thesis is dedicated to a study about how to best incorporate information from the PN, PM, SMR and EOB approaches in a synergistic way. We also discuss how accurate the resulting waveforms are, as compared against numerical-relativity (NR) simulations. We begin by comparing PM models, whether alone or recast in the EOB framework, against PN models and NR simulations. We will show that PM information has the potential to improve currently-employed models for LIGO and Virgo, especially if recast within the EOB formalism. This is very important, as the PM approximation comes with a host of new computational techniques from particle physics to exploit. Then, we show how a combination of PM and SMR approximations can be employed to access previously-unknown PN orders, deriving the third subleading PN dynamics for spin-orbit and (aligned) spin1-spin2 couplings. Such new results can then be included in the EOB models currently used in GW searches and parameter estimation studies, thereby improving them when the binaries have high spins. Finally, we build an EOB model for quasi-circular nonspinning binaries based on the SMR approximation (rather than the PN one as usually done). We show how this is done in detail without incurring in the divergences that had affected previous attempts, and compare the resultant model against NR simulations. We find that the SMR approximation is an excellent approximation for all (quasi-circular nonspinning) binaries, including both the equal-mass binaries that are routinely detected in GW searches and the ones with highly asymmetric masses. In particular, the SMR-based models compare much better than the PN models, suggesting that SMR-informed EOB models might be the key to model binaries in the future. In the second task of this thesis, we work within the linear-signal ap- proximation and describe generic metrics to predict inference biases on the parameters of a GW source of interest in the presence of confusion noise from unfitted foregrounds and from residuals of other signals that have been incorrectly fitted out. We illustrate the formalism with simple (yet realistic) LISA sources, and demonstrate its validity against Monte-Carlo simulations. The metrics we describe pave the way for more realistic studies to quantify the biases with future ground-based and spaceborne detectors.
Active matter broadly covers the dynamics of self-propelled particles.
While the onset of collective behavior in homogenous active systems is relatively well understood, the effect of inhomogeneities such as obstacles and traps lacks overall clarity.
Here, we study how interacting, self-propelled particles become trapped and released from a trap.
We have found that captured particles aggregate into an orbiting condensate with a crystalline structure. As more particles are added, the trapped condensates escape as a whole.
Our results shed light on the effects of confinement and quenched disorder in active matter.
Dark matter constraints from a joint analysis of dwarf Spheroidal galaxy observations with VERITAS
(2017)
We present constraints on the annihilation cross section of weakly interacting massive particles dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of similar to 230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is 1.35 x 10(-23) cm(3) s(-1) at 1 TeV for the bottom quark (b (b) over bar) final state, 2.85 x 10(-24) cm(3) s(-1) at 1 TeV for the tau lepton (tau+tau(-)) final state and 1.32 x 10-25 cm(3) s(-1) at 1 TeV for the gauge boson (gamma gamma) final state.
Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727 + 502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations. (C) 2017 Elsevier B.V. All rights reserved.
Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy, and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance, there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, is primarily utilized for gamma-ray astronomy but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 +/- 40(stat) +/- 140(syst) GeV.
We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 TeV to 500 TeV; The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law dF/dE = f(0) center dot (E/E-0)(-gamma) over the full energy range, with gamma = 2.82 +/- 0.30(stat)(-0.27)(+0.24)(syst) and f(0) = (4.82 +/- 0.98(stat)(-2.70)(+2.12)(syst)) x 10(-7) m(-2) s(-1) TeV-1 sr(-1) at E-0 = 50 TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies.
A Search for Pulsed Very High-energy Gamma-Rays from 13 Young Pulsars in Archival VERITAS Data
(2019)
We conduct a search for periodic emission in the very high-energy (VHE) gamma-ray band (E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hr. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in nondetections of pulsed VHE gamma-rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar.
HESS J1943+213
(2018)
HESS J1943+213 is a very high energy (VHE; > 100 GeV) gamma-ray source in the direction of the Galactic plane. Studies exploring the classification of the source are converging toward its identification as an extreme synchrotron BL Lac object. Here we present 38 hr of VERITAS observations of HESS J1943+213 taken over 2 yr. The source is detected with a significance of similar to 20 standard deviations, showing a remarkably stable flux and spectrum in VHE gamma-rays. Multifrequency Very Long Baseline Array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with the European VLBI Network and detect this component in the 4.6 and 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, Fermi-LAT, Swift-XRT, the FLWO 48 ' telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron self-Compton model. The well-measured gamma-ray peak of the SED with VERITAS and Fermi-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary gamma-rays from ultra-high-energy cosmic-ray-initiated electromagnetic cascades to the gamma-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and gamma-ray bands. No statistically significant flux or spectral variability is detected.
Kosmologie beschreibt die Entwicklung des Universums als Ganzes. Kosmologische Entdeckungen in Theorie und Praxis haben daher unser modernes wissenschaftliches Weltbild entscheidend geprägt. Die Vermittlung eines modernen Weltbildes durch Unterricht ist ein häufiger Wunsch in der naturwissenschaftlichen Bildungsdiskussion. Dennoch existieren weiterhin Forschungs- und Entwicklungsbedarfe. Kosmologische Themen finden sich häufig in den Medien und sind gleichzeitig weiter vom Alltag entfernt, so dass sich hier besonders leicht wissenschaftlich inkorrekte Vorstellungen entwickeln können, die zu Problemen im Unterricht führen können.
Das Ziel dieser wissenschaftlichen Arbeit ist es, zu diesem Forschungsgebiet beizutragen und die Voraussetzungen hinsichtlich vorhandener Vorkenntnisse und Präkonzepte in Kosmologie, mit denen Schülerinnen und Schüler in den Unterricht kommen, zu untersuchen und anschließend mit denen anderer Länder zu vergleichen. Dies erfolgt anhand einer qualitativen Inhaltsanalyse eines offenen Fragebogens. Auf dieser Grundlage wird schließlich ein Multiple-Choice Fragebogen entwickelt, angewendet und evaluiert.
Die Ergebnisse zeigen große Wissenslücken im Bereich der Kosmologie auf und geben erste Hinweise auf vorhandene Unterschiede zwischen den Ländern. Es existieren ebenfalls einige teils weit verbreitete wissenschaftlich inkorrekte Vorstellungen wie beispielsweise die Assoziation des Urknalls mit einer Explosion, der Urknall verursacht durch eine Kollision von Teilchen oder größeren Objekten, oder die Vorstellung der Ausdehnung des Universums als neue Entdeckungen und/oder Wissen. Des Weiteren gab nur etwa jeder Fünfte das korrekte Alter des Universums oder die Ausdehnung des Universums als einen der drei Belege der Urknalltheorie an, während fast 40% keinen einzigen Beleg nennen konnten. Für den geschlossenen Fragebogen konnten gute Hinweise für verschiedene Validitätsaspekte herausgearbeitet werden und es existieren erste Hinweise darauf, dass der Fragebogen Wissenszuwachs messen kann und damit wahrscheinlich zur Untersuchung der Wirksamkeit von Lerneinheiten eingesetzt werden kann. Auch ein entsprechendes Modell zur Verständnisentwicklung der Ausdehnung des Universums zeigte sich vielversprechend.
Diese Arbeit liefert insgesamt einen Forschungsbeitrag zum Schülervorwissen und Vorstellungen in der Kosmologie und deren Large Scale Assessment. Dies eröffnet die Möglichkeit zukünftiger Forschungen im Bereich von Gruppenvergleichen insbesondere hinsichtlich objektiver Ländervergleiche sowie der Untersuchungen der Wirksamkeit von einzelnen Lerneinheiten als auch Vergleiche verschiedener Lerneinheiten untereinander.
A Shockley-Type polymer
(2018)
Charge extraction rate in solar cells made of blends of electron donating/accepting organic semiconductors is typically slow due to their low charge carrier mobility. This sets a limit on the active layer thickness and has hindered the industrialization of organic solar cells (OSCs). Herein, charge transport and recombination properties of an efficient polymer (NT812):fullerene blend are investigated. This system delivers power conversion efficiency of >9% even when the junction thickness is as large as 800 nm. Experimental results indicate that this material system exhibits exceptionally low bimolecular recombination constant, 800 times smaller than the diffusion-controlled electron and hole encounter rate. Comparing theoretical results based on a recently introduced modified Shockley model for fill factor, and experiments, clarifies that charge collection is nearly ideal in these solar cells even when the thickness is several hundreds of nanometer. This is the first realization of high-efficiency Shockley-type organic solar cells with junction thicknesses suitable for scaling up.
Future magnetic recording industry needs a high-density data storage technology. However, switching the magnetization of small bits requires high magnetic fields that cause excessive heat dissipation. Therefore, controlling magnetism without applying external magnetic field is an important research topic for potential applications in data storage devices with low power consumption. Among the different approaches being investigated, two of them stand out, namely i) all-optical helicity dependent switching (AO-HDS) and ii) ferroelectric control of magnetism. This thesis aims to contribute towards a better understanding of the physical processes behinds these effects as well as reporting new and exciting possibility for the optical and/or electric control of magnetic properties. Hence, the thesis contains two differentiated chapters of results; the first devoted to AO-HDS on TbFe alloys and the second to the electric field control of magnetism in an archetypal Fe/BaTiO3 system.
In the first part, the scalability of the AO-HDS to small laser spot-sizes of few microns in the ferrimagnetic TbFe alloy is investigated by spatially resolving the magnetic contrast with photo-emission electron microscopy (PEEM) and X-ray magnetic circular dichroism (XMCD). The results show that the AO-HDS is a local effect within the laser spot size that occurs in the ring-shaped region in the vicinity of thermal demagnetization. Within the ring region, the helicity dependent switching occurs via thermally activated domain wall motion. Further, the thesis reports on a novel effect of thickness dependent inversion of the switching orientation. It addresses some of the important questions like the role of laser heating and the microscopic mechanism driving AO-HDS.
The second part of the thesis focuses on the electric field control of magnetism in an artificial multiferroic heterostructure. The sample consists of an Fe wedge with thickness varying between 0:5 nm and 3 nm, deposited on top of a ferroelectric and ferroelastic BaTiO3 [001]-oriented single crystal substrate. Here, the magnetic contrast is imaged via PEEM and XMCD as a function of out-of-plane voltage. The results show the evidence of the electric field control of superparamagnetism mediated by a ferroelastic modification of the magnetic anisotropy. The changes in the magnetoelastic anisotropy drive the transition from the superparamagnetic to superferromagnetic state at localized sample positions.
In this paper, we introduce the phenomenon of light driven diffusioosmotic long-range attraction and repulsion of porous particles under irradiation with UV light. The change in the inter-particle interaction potential is governed by flow patterns generated around single colloids and results in reversible aggregation or separation of the mesoporous silica particles that are trapped at a solid surface. The range of the interaction potential extends to several times the diameter of the particle and can be adjusted by varying the light intensity. The "fuel" of the process is a photosensitive surfactant undergoing photo-isomerization from a more hydrophobic trans-state to a rather hydrophilic cis-state. The surfactant has different adsorption affinities to the particles depending on the isomerization state. The trans-isomer, for example, tends to accumulate in the negatively charged pores of the particles, while the cis-isomer prefers to remain in the solution. This implies that when under UV irradiation cis-isomers are being formed within the pores, they tend to diffuse out readily and generate an excess concentration near the colloid's outer surface, ultimately resulting in the initiation of diffusioosmotic flow. The direction of the flow depends strongly on the dynamic redistribution of the fraction of trans- and cis-isomers near the colloids due to different kinetics of photo-isomerization within the pores as compared to the bulk. The unique feature of the mechanism discussed in the paper is that the long-range mutual repulsion but also the attraction can be tuned by convenient external optical stimuli such as intensity so that a broad variety of experimental situations for manipulation of a particle ensemble can be realized.
We report on photoisomerization kinetics of azobenzene containing surfactants in aqueous solution. The surfactant molecule consists of a positively charged trimethylammonium bromide head group, a hydrophobic spacer connecting via 6 to 10 CH2 groups to the azobenzene unit, and the hydrophobic tail of 1 and 3CH(2) groups. Under exposure to light, the azobenzene photoisomerizes from more stable trans- to metastable cis-state, which can be switched back either thermally in dark or by illumination with light of a longer wavelength. The surfactant isomerization is described by a kinetic model of a pseudo first order reaction approaching equilibrium, where the intensity controls the rate of isomerization until the equilibrated state. The rate constants of the trans-cis and cis-trans photoisomerization are calculated as a function of several parameters such as wavelength and intensity of light, the surfactant concentration, and the length of the hydrophobic tail. The thermal relaxation rate from cis- to trans-state is studied as well. The surfactant isomerization shows a different kinetic below and above the critical micellar concentration of the trans isomer due to steric hindrance within the densely packed micelle but does not depend on the spacer length.
The Earth's inner magnetosphere is a very dynamic system, mostly driven by the external solar wind forcing exerted upon the magnetic field of our planet. Disturbances in the solar wind, such as coronal mass ejections and co-rotating interaction regions, cause geomagnetic storms, which lead to prominent changes in charged particle populations of the inner magnetosphere - the plasmasphere, ring current, and radiation belts. Satellites operating in the regions of elevated energetic and relativistic electron fluxes can be damaged by deep dielectric or surface charging during severe space weather events. Predicting the dynamics of the charged particles and mitigating their effects on the infrastructure is of particular importance, due to our increasing reliance on space technologies.
The dynamics of particles in the plasmasphere, ring current, and radiation belts are strongly coupled by means of collisions and collisionless interactions with electromagnetic fields induced by the motion of charged particles. Multidimensional numerical models simplify the treatment of transport, acceleration, and loss processes of these particles, and allow us to predict how the near-Earth space environment responds to solar storms. The models inevitably rely on a number of simplifications and assumptions that affect model accuracy and complicate the interpretation of the results. In this dissertation, we quantify the processes that control electron dynamics in the inner magnetosphere, paying particular attention to the uncertainties of the employed numerical codes and tools.
We use a set of convenient analytical solutions for advection and diffusion equations to test the accuracy and stability of the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. We show that numerical schemes implemented in the code converge to the analytical solutions and that the VERB-4D code demonstrates stable behavior independent of the assumed time step. The order of the numerical scheme for the convection equation is demonstrated to affect results of ring current and radiation belt simulations, and it is crucially important to use high-order numerical schemes to decrease numerical errors in the model.
Using the thoroughly tested VERB-4D code, we model the dynamics of the ring current electrons during the 17 March 2013 storm. The discrepancies between the model and observations above 4.5 Earth's radii can be explained by uncertainties in the outer boundary conditions. Simulation results indicate that the electrons were transported from the geostationary orbit towards the Earth by the global-scale electric and magnetic fields.
We investigate how simulation results depend on the input models and parameters. The model is shown to be particularly sensitive to the global electric field and electron lifetimes below 4.5 Earth's radii. The effects of radial diffusion and subauroral polarization streams are also quantified.
We developed a data-assimilative code that blends together a convection model of energetic electron transport and loss and Van Allen Probes satellite data by means of the Kalman filter. We show that the Kalman filter can correct model uncertainties in the convection electric field, electron lifetimes, and boundary conditions. It is also demonstrated how the innovation vector - the difference between observations and model prediction - can be used to identify physical processes missing in the model of energetic electron dynamics.
We computed radial profiles of phase space density of ultrarelativistic electrons, using Van Allen Probes measurements. We analyze the shape of the profiles during geomagnetically quiet and disturbed times and show that the formation of new local minimums in the radial profiles coincides with the ground observations of electromagnetic ion-cyclotron (EMIC) waves. This correlation indicates that EMIC waves are responsible for the loss of ultrarelativistic electrons from the heart of the outer radiation belt into the Earth's atmosphere.
Here, a promising approach for producing piezo-polymer transducers in a one-step process is presented. Using 3D-printing technology and polypropylene (PP) filaments, we are able to print a two-layered film structure with regular cavities of precisely controlled size and shape. It is found that the 3D-printed samples exhibit piezoelectric coefficients up to 200 pC/N, similar to those of other PP ferroelectrets, and their temporal and thermal behavior is in good agreement with those known of PP ferroelectrets. The piezoelectric response strongly decreases for applied pressures above 20 kPa, as the pressure in the air-filled cavities strongly influences the overall elastic modulus of ferroelectrets.
Recurrence plot analyses suggest a novel reference system involved in newborn spontaneous movements
(2006)
The movements of newborns have been thoroughly studied in terms of reflexes, muscle synergies, leg coordination, and target-directed arm/hand movements. Since these approaches have concentrated mainly on separate accomplishments, there has remained a clear need for more integrated investigations. Here, we report an inquiry in which we explicitly concentrated on taking such a perspective and, additionally, were guided by the methodological concept of home base behavior, which Ilan Golard developed for studies of exploratory behavior in animals. Methods from nonlinear dynamics, such as symbolic dynamics and recurrence plot analyses of kinematic data received from audiovisual newborn recordings, yielded new insights into the spatial and temporal organization of limb movements. In the framework of home base behavior, our approach uncovered a novel reference system of spontaneous newborn movements.
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.
The Arctic is the hot spot of the ongoing, global climate change. Over the last decades, near-surface temperatures in the Arctic have been rising almost four times faster than on global average. This amplified warming of the Arctic and the associated rapid changes of its environment are largely influenced by interactions between individual components of the Arctic climate system. On daily to weekly time scales, storms can have major impacts on the Arctic sea-ice cover and are thus an important part of these interactions within the Arctic climate. The sea-ice impacts of storms are related to high wind speeds, which enhance the drift and deformation of sea ice, as well as to changes in the surface energy budget in association with air mass advection, which impact the seasonal sea-ice growth and melt.
The occurrence of storms in the Arctic is typically associated with the passage of transient cyclones. Even though the above described mechanisms how storms/cyclones impact the Arctic sea ice are in principal known, there is a lack of statistical quantification of these effects. In accordance with that, the overarching objective of this thesis is to statistically quantify cyclone impacts on sea-ice concentration (SIC) in the Atlantic Arctic Ocean over the last four decades. In order to further advance the understanding of the related mechanisms, an additional objective is to separate dynamic and thermodynamic cyclone impacts on sea ice and assess their relative importance. Finally, this thesis aims to quantify recent changes in cyclone impacts on SIC. These research objectives are tackled utilizing various data sets, including atmospheric and oceanic reanalysis data as well as a coupled model simulation and a cyclone tracking algorithm.
Results from this thesis demonstrate that cyclones are significantly impacting SIC in the Atlantic Arctic Ocean from autumn to spring, while there are mostly no significant impacts in summer. The strength and the sign (SIC decreasing or SIC increasing) of the cyclone impacts strongly depends on the considered daily time scale and the region of the Atlantic Arctic Ocean. Specifically, an initial decrease in SIC (day -3 to day 0 relative to the cyclone) is found in the Greenland, Barents and Kara Seas, while SIC increases following cyclones (day 0 to day 5 relative to the cyclone) are mostly limited to the Barents and Kara Seas.
For the cold season, this results in a pronounced regional difference between overall (day -3 to day 5 relative to the cyclone) SIC-decreasing cyclone impacts in the Greenland Sea and overall SIC-increasing cyclone impacts in the Barents and Kara Seas. A cyclone case study based on a coupled model simulation indicates that both dynamic and thermodynamic mechanisms contribute to cyclone impacts on sea ice in winter. A typical pattern consisting of an initial dominance of dynamic sea-ice changes followed by enhanced thermodynamic ice growth after the cyclone passage was found. This enhanced ice growth after the cyclone passage most likely also explains the (statistical) overall SIC-increasing effects of cyclones in the Barents and Kara Seas in the cold season.
Significant changes in cyclone impacts on SIC over the last four decades have emerged throughout the year. These recent changes are strongly varying from region to region and month to month. The strongest trends in cyclone impacts on SIC are found in autumn in the Barents and Kara Seas. Here, the magnitude of destructive cyclone impacts on SIC has approximately doubled over the last four decades. The SIC-increasing effects following the cyclone passage have particularly weakened in the Barents Sea in autumn. As a consequence, previously existing overall SIC-increasing cyclone impacts in this region in autumn have recently disappeared. Generally, results from this thesis show that changes in the state of the sea-ice cover (decrease in mean sea-ice concentration and thickness) and near-surface air temperature are most important for changed cyclone impacts on SIC, while changes in cyclone properties (i.e. intensity) do not play a significant role.
Corvino, Corvino and Schoen, Chruściel and Delay have shown the existence of a large class of asymptotically flat vacuum initial data for Einstein's field equations which are static or stationary in a neighborhood of space-like infinity, yet quite general in the interior. The proof relies on some abstract, non-constructive arguments which makes it difficult to calculate such data numerically by using similar arguments. A quasilinear elliptic system of equations is presented of which we expect that it can be used to construct vacuum initial data which are asymptotically flat, time-reflection symmetric, and asymptotic to static data up to a prescribed order at space-like infinity. A perturbation argument is used to show the existence of solutions. It is valid when the order at which the solutions approach staticity is restricted to a certain range. Difficulties appear when trying to improve this result to show the existence of solutions that are asymptotically static at higher order. The problems arise from the lack of surjectivity of a certain operator. Some tensor decompositions in asymptotically flat manifolds exhibit some of the difficulties encountered above. The Helmholtz decomposition, which plays a role in the preparation of initial data for the Maxwell equations, is discussed as a model problem. A method to circumvent the difficulties that arise when fast decay rates are required is discussed. This is done in a way that opens the possibility to perform numerical computations. The insights from the analysis of the Helmholtz decomposition are applied to the York decomposition, which is related to that part of the quasilinear system which gives rise to the difficulties. For this decomposition analogous results are obtained. It turns out, however, that in this case the presence of symmetries of the underlying metric leads to certain complications. The question, whether the results obtained so far can be used again to show by a perturbation argument the existence of vacuum initial data which approach static solutions at infinity at any given order, thus remains open. The answer requires further analysis and perhaps new methods.
The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function.
Anomalous diffusion with a power-law time dependence vertical bar R vertical bar(2)(t) similar or equal to t(alpha i) of the mean squared displacement occurs quite ubiquitously in numerous complex systems. Often, this anomalous diffusion is characterised by crossovers between regimes with different anomalous diffusion exponents alpha(i). Here we consider the case when such a crossover occurs from a first regime with alpha(1) to a second regime with alpha(2) such that alpha(2) > alpha(1), i.e., accelerating anomalous diffusion. A widely used framework to describe such crossovers in a one-dimensional setting is the bi-fractional diffusion equation of the so-called modified type, involving two time-fractional derivatives defined in the Riemann-Liouville sense. We here generalise this bi-fractional diffusion equation to higher dimensions and derive its multidimensional propagator (Green's function) for the general case when also a space fractional derivative is present, taking into consideration long-ranged jumps (Levy flights). We derive the asymptotic behaviours for this propagator in both the short- and long-time as well the short- and long-distance regimes. Finally, we also calculate the mean squared displacement, skewness and kurtosis in all dimensions, demonstrating that in the general case the non-Gaussian shape of the probability density function changes.
We study by Monte Carlo simulations a kinetic exchange trading model for both fixed and distributed saving propensities of the agents and rationalize the person and wealth distributions. We show that the newly introduced wealth distribution – that may be more amenable in certain situations – features a different power-law exponent, particularly for distributed saving propensities of the agents. For open agent-based systems, we analyze the person and wealth distributions and find that the presence of trap agents alters their amplitude, leaving however the scaling exponents nearly unaffected. For an open system, we show that the total wealth – for different trap agent densities and saving propensities of the agents – decreases in time according to the classical Kohlrausch–Williams–Watts stretched exponential law. Interestingly, this decay does not depend on the trap agent density, but rather on saving propensities. The system relaxation for fixed and distributed saving schemes are found to be different.
We study buckling instabilities of filaments in biological systems. Filaments in a cell are the building blocks of the cytoskeleton. They are responsible for the mechanical stability of cells and play an important role in intracellular transport by molecular motors, which transport cargo such as organelles along cytoskeletal filaments. Filaments of the cytoskeleton are semiflexible polymers, i.e., their bending energy is comparable to the thermal energy such that they can be viewed as elastic rods on the nanometer scale, which exhibit pronounced thermal fluctuations. Like macroscopic elastic rods, filaments can undergo a mechanical buckling instability under a compressive load. In the first part of the thesis, we study how this buckling instability is affected by the pronounced thermal fluctuations of the filaments. In cells, compressive loads on filaments can be generated by molecular motors. This happens, for example, during cell division in the mitotic spindle. In the second part of the thesis, we investigate how the stochastic nature of such motor-generated forces influences the buckling behavior of filaments. In chapter 2 we review briefly the buckling instability problem of rods on the macroscopic scale and introduce an analytical model for buckling of filaments or elastic rods in two spatial dimensions in the presence of thermal fluctuations. We present an analytical treatment of the buckling instability in the presence of thermal fluctuations based on a renormalization-like procedure in terms of the non-linear sigma model where we integrate out short-wavelength fluctuations in order to obtain an effective theory for the mode of the longest wavelength governing the buckling instability. We calculate the resulting shift of the critical force by fluctuation effects and find that, in two spatial dimensions, thermal fluctuations increase this force. Furthermore, in the buckled state, thermal fluctuations lead to an increase in the mean projected length of the filament in the force direction. As a function of the contour length, the mean projected length exhibits a cusp at the buckling instability, which becomes rounded by thermal fluctuations. Our main result is the observation that a buckled filament is stretched by thermal fluctuations, i.e., its mean projected length in the direction of the applied force increases by thermal fluctuations. Our analytical results are confirmed by Monte Carlo simulations for buckling of semiflexible filaments in two spatial dimensions. We also perform Monte Carlo simulations in higher spatial dimensions and show that the increase in projected length by thermal fluctuations is less pronounced than in two dimensions and strongly depends on the choice of the boundary conditions. In the second part of this work, we present a model for buckling of semiflexible filaments under the action of molecular motors. We investigate a system in which a group of motors moves along a clamped filament carrying a second filament as a cargo. The cargo-filament is pushed against the wall and eventually buckles. The force-generating motors can stochastically unbind and rebind to the filament during the buckling process. We formulate a stochastic model of this system and calculate the mean first passage time for the unbinding of all linking motors which corresponds to the transition back to the unbuckled state of the cargo filament in a mean-field model. Our results show that for sufficiently short microtubules the movement of kinesin-I-motors is affected by the load force generated by the cargo filament. Our predictions could be tested in future experiments.
The H.E.S.S. array is a third generation Imaging Atmospheric Cherenkov Telescope (IACT) array. It is located in the Khomas Highland in Namibia, and measures very high energy (VHE) gamma-rays. In Phase I, the array started data taking in 2004 with its four identical 13 m telescopes. Since then, H.E.S.S. has emerged as the most successful IACT experiment to date. Among the almost 150 sources of VHE gamma-ray radiation found so far, even the oldest detection, the Crab Nebula, keeps surprising the scientific community with unexplained phenomena such as the recently discovered very energetic flares of high energy gamma-ray radiation. During its most recent flare, which was detected by the Fermi satellite in March 2013, the Crab Nebula was simultaneously observed with the H.E.S.S. array for six nights. The results of the observations will be discussed in detail during the course of this work. During the nights of the flare, the new 24 m × 32 m H.E.S.S. II telescope was still being commissioned, but participated in the data taking for one night. To be able to reconstruct and analyze the data of the H.E.S.S. Phase II array, the algorithms and software used by the H.E.S.S. Phase I array had to be adapted. The most prominent advanced shower reconstruction technique developed by de Naurois and Rolland, the template-based model analysis, compares real shower images taken by the Cherenkov telescope cameras with shower templates obtained using a semi-analytical model. To find the best fitting image, and, therefore, the relevant parameters that describe the air shower best, a pixel-wise log-likelihood fit is done. The adaptation of this advanced shower reconstruction technique to the heterogeneous H.E.S.S. Phase II array for stereo events (i.e. air showers seen by at least two telescopes of any kind), its performance using MonteCarlo simulations as well as its application to real data will be described.
Over the past decades, there has been a growing interest in ‘extreme events’ owing to the increasing threats that climate-related extremes such as floods, heatwaves, droughts, etc., pose to society. While extreme events have diverse definitions across various disciplines, ranging from earth science to neuroscience, they are characterized mainly as dynamic occurrences within a limited time frame that impedes the normal functioning of a system. Although extreme events are rare in occurrence, it has been found in various hydro-meteorological and physiological time series (e.g., river flows, temperatures, heartbeat intervals) that they may exhibit recurrent behavior, i.e., do not end the lifetime of the system. The aim of this thesis to develop some
sophisticated methods to study various properties of extreme events.
One of the main challenges in analyzing such extreme event-like time series is that they have large temporal gaps due to the paucity of the number of observations of extreme events. As a result, existing time series analysis tools are usually not helpful to decode the underlying
information. I use the edit distance (ED) method to analyze extreme event-like time series in their unaltered form. ED is a specific distance metric, mainly designed to measure the similarity/dissimilarity between point process-like data. I combine ED with recurrence plot techniques to identify the recurrence property of flood events in the Mississippi River in the United States. I also use recurrence quantification analysis to show the deterministic properties
and serial dependency in flood events.
After that, I use this non-linear similarity measure (ED) to compute the pairwise dependency in extreme precipitation event series. I incorporate the similarity measure within the framework of complex network theory to study the collective behavior of climate extremes. Under this architecture, the nodes are defined by the spatial grid points of the given spatio-temporal climate dataset. Each node is associated with a time series corresponding to the temporal evolution
of the climate observation at that grid point. Finally, the network links are functions of the pairwise statistical interdependence between the nodes. Various network measures, such as degree, betweenness centrality, clustering coefficient, etc., can be used to quantify the network’s topology. We apply the methodology mentioned above to study the spatio-temporal coherence pattern of extreme rainfall events in the United States and the Ganga River basin, which reveals its relation to various climate processes and the orography of the region.
The identification of precursors associated with the occurrence of extreme events in the near future is extremely important to prepare the masses for an upcoming disaster and mitigate the potential risks associated with such events. Under this motivation, I propose an in-data prediction recipe for predicting the data structures that typically occur prior to extreme events using the Echo state network, a type of Recurrent Neural Network which is a part of the reservoir
computing framework. However, unlike previous works that identify precursory structures in the same variable in which extreme events are manifested (active variable), I try to predict these structures by using data from another dynamic variable (passive variable) which does not show large excursions from the nominal condition but carries imprints of these extreme events. Furthermore, my results demonstrate that the quality of prediction depends on the magnitude
of events, i.e., the higher the magnitude of the extreme, the better is its predictability skill. I show quantitatively that this is because the input signals collectively form a more coherent pattern for an extreme event of higher magnitude, which enhances the efficiency of the machine to predict the forthcoming extreme events.
Classical semiconductor physics has been continuously improving electronic components such as diodes, light-emitting diodes, solar cells and transistors based on highly purified inorganic crystals over the past decades. Organic semiconductors, notably polymeric, are a comparatively young field of research, the first light-emitting diode based on conjugated polymers having been demonstrated in 1990. Polymeric semiconductors are of tremendous interest for high-volume, low-cost manufacturing ("printed electronics"). Due to their rather simple device structure mostly comprising only one or two functional layers, polymeric diodes are much more difficult to optimize compared to small-molecular organic devices. Usually, functions such as charge injection and transport are handled by the same material which thus needs to be highly optimized. The present work contributes to expanding the knowledge on the physical mechanisms determining device performance by analyzing the role of charge injection and transport on device efficiency for blue and white-emitting devices, based on commercially relevant spiro-linked polyfluorene derivatives. It is shown that such polymers can act as very efficient electron conductors and that interface effects such as charge trapping play the key role in determining the overall device efficiency. This work contributes to the knowledge of how charges drift through the polymer layer to finally find neutral emissive trap states and thus allows a quantitative prediction of the emission color of multichromophoric systems, compatible with the observed color shifts upon driving voltage and temperature variation as well as with electrical conditioning effects. In a more methodically oriented part, it is demonstrated that the transient device emission observed upon terminating the driving voltage can be used to monitor the decay of geminately-bound species as well as to determine trapped charge densities. This enables direct comparisons with numerical simulations based on the known properties of charge injection, transport and recombination. The method of charge extraction under linear increasing voltages (CELIV) is investigated in some detail, correcting for errors in the published approach and highlighting the role of non-idealized conditions typically present in experiments. An improved method is suggested to determine the field dependence of charge mobility in a more accurate way. Finally, it is shown that the neglect of charge recombination has led to a misunderstanding of experimental results in terms of a time-dependent mobility relaxation.
We analysed recent K2 data of the short-period eclipsing binary system HW Vir, which consists of a hot subdwarf-B type primary with an M-dwarf companion. We determined the mid-times of eclipses, calculated O-C diagrams, and an average shift of the secondary minimum. Our results show that the orbital period is stable within the errors over the course of the 70 days of observations. Interestingly, the offset from mid-orbital phase between the primary and the secondary eclipses is found to be 1.62 s. If the shift is explained solely by light-travel time, the mass of the sdB primary must be 0.26 M-circle dot, which is too low for the star to be core-helium burning. However, we argue that this result is unlikely to be correct and that a number of effects caused by the relative sizes of the stars conspire to reduce the effective light-travel time measurement. After removing the flux variation caused by the orbit, we calculated the amplitude spectrum to search for pulsations. The spectrum clearly shows periodic signal from close to the orbital frequency up to 4600 mu Hz, with the majority of peaks found below 2600 mu Hz. The amplitudes are below 0.1 part-per-thousand, too low to be detected with ground-based photometry. Thus, the high-precision data from the Kepler spacecraft has revealed that the primary of the HW Vir system is a pulsating sdBV star. We argue that the pulsation spectrum of the primary in HW Vir differs from that in other sdB stars due to its relatively fast rotation that is (nearly) phase-locked with the orbit.
We present an analysis of two pulsating subdwarf B stars PHL 457 and EQ Psc observed during the K2 mission. The K2 light curves of both stars show variation consistent with irradiation of a cooler companion by the hot subdwarf. They also show higher frequency oscillations consistent with pulsation. Using new spectroscopic data, we measured the radial velocity, effective temperature, surface gravity, and helium abundance of both hot subdwarfs as a function of orbital phase. We confirm the previously published spectroscopic orbit of PHL 457, and present the first spectroscopic orbit of EQ Psc. The orbital periods are 0.313 and 0.801 d, respectively. For EQPsc, we find a strong correlation between T-eff and orbital phase, due to contribution of light from the irradiated companion. We calculated amplitude spectra, identified significant pulsation frequencies, and searched for multiplets and asymptotic period spacings. By means of multiplets and period spacing, we identified the degrees of several pulsation modes in each star. The g-mode multiplets indicate subsynchronous core rotation with periods of 4.6 d (PHL 457) and 9.4 d (EQ Psc). We made spectral energy disctribution (SED) fits of PHL 457 and EQ Psc using available broad-band photometry and Gaia data. While the SED of PHL 457 shows no evidence of a cool companion, the SED for EQPsc clearly shows an infrared (IR) access consistent with a secondary with a temperature of about 6800K and a radius of 0.23 R-circle dot. This is the first detection of an IR access in any sdB + dM binary.
Carbonfasern haben sich in der Luft- und Raumfahrt etabliert und gewinnen in Alltagsanwendungen wie dem Automobilbereich, Windkraft- und Sportbereich durch ihre hohen Zugfestigkeiten, insbesondere ihrer hohen E-Moduli, und ihrer geringen Dichte immer mehr an Bedeutung. Auf Grund ihrer hohen Kosten, welche sich zur Hälfte aus der Precursorherstellung, inklusive seiner Synthese und seinem Verspinnprozess, dem Lösungsspinnverfahren, ergeben, erhalten zunehmend alternative und schmelzspinnbare Precursoren Interesse. Für die Carbonfaserherstellung wird fast ausschließlich Polyacrylnitril (PAN) verwendet, das vor dem Schmelzen irreversible exotherme Zyklisierungsreaktionen aufweist, welchen sich seine Zersetzung anschließt. Eine Möglichkeit der Reduzierung der Schmelztemperatur von Polymeren ist die Einbringung von Comonomeren zur Erhöhung des freien Volumens und die Reduzierung der intermolekularen Wechselwirkungen als interne Weichmacher. Wie am Fraunhofer IAP gezeigt wurde, kann mittels 2-Methoxyethylacrylat (MEA) die Schmelztemperatur zu neuartigen PAN-basierten Precursoren verringert werden. Um den PAN-co-MEA-Precursor für die nachfolgenden Prozessschritte der Carbonfaserherstellung zu verwenden, müssen die thermoplastischen Fasern in thermisch stabile Fasern ohne thermoplastisches Verhalten überführt werden. Es wurde ein neuer Prozessschritt (Prästabilisierung) eingeführt, welcher unter alkalischen Bedingungen zur Abspaltung der Comonomerseitenkette führt. Neben der Esterhydrolyse finden Reaktionen statt, welche an diesem Material noch nicht hinreichend untersucht wurden. Weiterhin stellt sich die Frage nach der Kinetik der Prästabilisierung und der Ermittlung einer geeigneten Prozessführung.
Hierzu wurde die Prästabilisierung in den Labormaßstab überführt und die möglichen Zusammensetzungen des aus DMSO und einer KOH-Lösung bestehenden Reaktionsmediums evaluiert. Weiterhin wurde die Behandlung bei verschiedenen Prästabilisierungszeiten von maximal 30 min und Temperaturen von 40, 50 und 60 °C durchgeführt, um primär mittels NMR-Spektroskopie die chemischen Strukturänderungen aufzuklären. Die Esterhydrolyse des Comonomers, welche zur Abspaltung des 2-Methoxyethanols führt, wurde mittels 1H-NMR-spektroskopischer Untersuchungen detektiert.
Es wurde ein Modell aufgestellt, das die chemisch-physikalischen Strukturänderungen während der Prästabilisierung aufzeigt. Die zuerst ablaufende Reaktion ist die Esterhydrolyse am Comonomer, welche vom Faserrand nach innen verläuft und durch die Präsenz des DMSO in Kombination mit der KOH-Lösung (Superbase) initiiert wird. Der zeitliche Reaktionsverlauf der Esterhydrolyse kann in drei Bereiche eingeteilt werden. Der erste Bereich ab dem Prästabilisierungsbeginn wird durch die Diffusion der basischen Anionen in die Faser, der zweite Bereich durch die Reaktion an der Estergruppe des Comonomers und der dritte Bereich durch letzte Reaktionen im Faserinneren und diffusiven Prozessen der Produkte und Edukte charakterisiert. Der zweite Bereich kann mit einer Reaktion pseudo 1. Ordnung abgebildet werden, da in diesem Bereich bereits eine ausreichende Diffusion der Edukte in die Faser stattgefunden hat. Bei 50 °C spielt die Diffusion im ersten Bereich im Vergleich zur Reaktion eine untergeordnete Rolle. Mit Erhöhung der Temperatur auf 60 °C kann eine im Verhältnis geringere Diffusions- als Reaktionsgeschwindigkeit beobachtet werden. Die Nebenreaktionen wurden mittels 13C-CP/MAS-NMR-spektroskopischen, elementaranlaytischen Untersuchungen sowie Doppelbrechungsmessungen charakterisiert. Während der alkalischen Esterhydrolyse beginnt die Reduzierung der Nitrilgruppen unter der Bildung von primären Carbonsäureamiden und Carbonsäuren. Zur Beschreibung dieser Umsetzung wurde eine Methode entwickelt, welche die Addition von 13C-CP/MAS-NMR-Spektren der Modellsubstanzen PAN, PAM und PAA beinhaltet. Weitere stattfindende Reaktionen sind die Bildung von konjugierten Doppelbindungen, welche insbesondere auf eine Zyklisierung der Nitrile hinweisen. Die nasschemisch initiierte Zyklisierung der Nitrilgruppen kann zu kürzeren Stabilisierungszeiten und einem besser kontrollierbaren Stabilisierungsprozess durch geringere Wärmefreisetzung und schlussendlich zu einer Kostenersparnis des gesamten Verfahrens führen. Die Umsetzung der Nitrilgruppen konnte mit einer Reaktion pseudo 1. Ordnung gut abgebildet werden. DMSO initiiert die Esterhydrolyse, wobei die KOH-Konzentration einen höheren Einfluss auf die Reaktionsgeschwindigkeit der Ester- und Nitrilhydrolyse als die DMSO-Konzentration besitzt. Beide Reaktionen zeigen eine vergleichbare Abhängigkeit von der Temperatur. Die Erhöhung der Prästabilisierungszeit und der KOH- bzw. DMSO-Konzentration führt zur Migration niedermolekularer Bestandteile des Fasermaterials an die Oberfläche und der Bildung punktueller Ablagerungen bis hin zu miteinander verbundenen Einzelfasern. Eine weitere Erhöhung der Prästabilisierungszeit bzw. der Konzentration führt zu einem steigenden Carbonsäureanteil und zur Quellung des Fasermaterials, wodurch die Ablagerungen in das Reaktionsmedium diffundieren. Die Ablagerungen enthalten Chlor, welches durch den Waschvorgang mit HCl in das Materialsystem gelangt ist und durch Parameteranpassungen reduziert wurde. Die schmelzbaren Fasern konnten durch die Prästabilisierung erfolgreich über eine Kern-Mantel-Struktur in nicht-thermoplastische Fasern überführt werden.
Zur Ermittlung eines geeigneten Prozessfensters für nachfolgende thermische Beanspruchungen der prästabilisierten Fasern wurden drei Kriterien identifiziert, anhand welcher die Evaluation erfolgte. Das erste Kriterium beinhaltet die Notwendigkeit der vollständigen Aufhebung der thermoplastischen Eigenschaft der Fasern. Als zweites Kriterium diente die Fasermorphologie. Anhand von REM-Aufnahmen wurden Faserbündel mit separierten Einzelfasern ohne Ablagerungen für die nachfolgende Stabilisierung ausgewählt. Das dritte Kriterium bezieht sich auf eine möglichst geringe Umsetzung der Nitrilgruppen, um Prästabilisierungsbedingungen mit Nebenreaktionen zu vermeiden.
Aus den Untersuchungen konnte eine Prästabilisierungstemperatur von 60 °C als geeignet identifiziert werden. Weiterhin führen hoch alkalische Zusammensetzungen des Reaktionsmediums mit KOH-Konzentrationen von 1, 1,5 und 2 M, vorzugsweise 1,5 M und 50 vol% DMSO mit Reaktionszeiten von unter 10 min zu geeigneten Fasern. Ein MEA-Anteil unterhalb von 2 mol% bewirkt eine Überführung in die Unschmelzbarkeit. Thermisch stabile und für die nachfolgende Stabilisierung geeignete Fasern besitzen weiterhin 68 – 80 mol% Nitrilgruppen, 20 – 25 mol% Carbonsäuren, bis zu 15 mol% primäre Carbonsäureamide und zyklisierte Strukturen.
We introduce non-Markovian SIR epidemic spreading model inspired by the characteristics of the COVID-19, by considering discrete-and continuous-time versions. The distributions of infection intensity and recovery period may take an arbitrary form. By taking corresponding choice of these functions, it is shown that the model reduces to the classical Markovian case. The epidemic threshold is analytically determined for arbitrary functions of infectivity and recovery and verified numerically. The relevance of the model is shown by modeling the first wave of the epidemic in Italy, Spain and the UK, in the spring, 2020.
Seasonal forecasts are of great interest in many areas. Knowing the amount of precipitation for the upcoming season in regions of water scarcity would facilitate a better water management. If farmers knew the weather conditions of the upcoming summer at sowing time, they could select those cereal species that are best adapted to these conditions. This would allow farmers to improve the harvest and potentially even reduce the amount of pesticides used. However, the undoubted advantages of seasonal forecasts are often opposed by their high degree of uncertainty. The great challenge of generating seasonal forecasts with lead times of several months mainly originates from the chaotic nature of the earth system. In a chaotic system, even tiny differences in the initial conditions can lead to strong deviations in the system’s state in the long run.
In this dissertation we propose an emergent machine learning approach for seasonal forecasting, called the AnlgModel. The AnlgModel combines the analogue method with myopic feature selection and bootstrapping. To benchmark the abilities of the AnlgModel we apply it to seasonal cyclone activity forecasts in the North Atlantic and Northwest Pacific. The AnlgModel demonstrates competitive hindcast skills with two operational forecasts and even outperforms these for long lead times.
In the second chapter we comprehend the forecasting strategy of the Anlg-Model. We thereby analyse the analogue selection process for the 2017 North Atlantic and the 2018 Northwest Pacific seasonal cyclone activity. The analysis shows that those climate indices which are known to influence the seasonal cyclone activity, such as the Niño 3.4 SST, are correctly represented among the selected analogues. Furthermore the selected analogues reflect large-scale climate patterns that were identified by expert reports as being determinative for these particular seasons.
In the third chapter we analyse the features that are used by the AnlgModel for its predictions. We therefore inspect the feature relevance (FR). The FR patterns learned by the AnlgModel show a high congruence with the predictor regions used by the operational forecasts. However, the AnlgModel also discovered new features, such as the SST anomaly in the Gulf of Guinea during November. This SST pattern exhibits a remarkably high predictive potential for the upcoming Atlantic hurricane activity.
In the final chapter we investigate potential mechanisms, that link two of these regions with high feature relevance to the Atlantic hurricane activity. We mainly focus on ocean surface transport. The ocean surface flow paths are calculated using Lagrangian particle analysis. We demonstrate that the FR patterns in the region of the Canary islands do not correspond with ocean surface transport. It is instead likely that these FR patterns fingerprint a wind transport of latent heat. The second region to be studied is situated in the Gulf of Guinea. Our analysis shows that the FR patterns seen there do fingerprint ocean surface transport. However, our simulations also show that at least one other mechanism is involved in linking the Gulf of Guinea SST anomaly in November to the hurricane activity of the upcoming season.
In this work the AnlgModel does not only demonstrate its outstanding forecast skills but also shows its capabilities as research tool for detecting oceanic and atmospheric mechanisms.
Bioinspired confinement of upconversion nanoparticles for improved performance in aqueous solution
(2020)
The resonance energy transfer (RET) from NaYF4:Yb,Er upconverting nanoparticles (UNCPs) to a dye (5-carboxytetramethylrhodamine (TAMRA)) was investigated by photoluminescence experiments and microscale thermophoresis (MST). The dye was excited via RET from the UCNPs which was excited in the near-infrared (NIR). The change of the dye diffusion speed (free vs coupled) was investigated by MST. RET shows significant changes in the decay times of the dye as well as of the UCNPs. MST reveals significant changes in the diffusion speed. A unique amphiphilic coating polymer (customized mussel protein (CMP) polymer) for UCNP surface coating was used, which mimics blood protein adsorption and mussel food protein adhesion to transfer the UCNP into the aqueous phase and to allow surface functionalization. The CMP provides very good water dispersibility to the UCNPs and minimizes ligand exchange and subsequent UCNP aging reactions because of the interlinkage of the CMP on the UCNP surface. Moreover, CMP provides N-3-functional groups for dick chemistry-based functionalization demonstrated with the dye 5-carboxytetramethylrhodamine (TAMRA). This establishes the principle coupling scheme for suitable biomarkers such as antibodies. The CMP provides very stable aqueous UCNP dispersions that are storable up to 3 years in a fridge at 5 degrees C without dissolution or coagulation. The outstanding properties of CMP in shielding the UCNP from unwanted solvent effects is reflected in the distinct increase of the photoluminescence decay times after UCNP functionalization. The UCNP-to-TAMRA energy transfer is also spectroscopically investigated at low temperatures (4-200 K), revealing that one of the two green Er(III) emission bands contributes the major part to the energy transfer. The TAMRA fluorescence decay time increases by a factor of 9500 from 2.28 ns up to 22 mu s due to radiationless energy transfer from the UCNP after NIR excitation of the latter. This underlines the unique properties of CMP as a versatile capping ligand for distinctly improving the UCNPs' performance in aqueous solutions, for coupling of biomolecules, and for applications for in vitro and in vivo experiments using UCNPs as optical probes in life science applications.
Moving spiral wave chimeras
(2021)
We consider a two-dimensional array of heterogeneous nonlocally coupled phase oscillators on a flat torus and study the bound states of two counter-rotating spiral chimeras, shortly two-core spiral chimeras, observed in this system. In contrast to other known spiral chimeras with motionless incoherent cores, the two-core spiral chimeras typically show a drift motion. Due to this drift, their incoherent cores become spatially modulated and develop specific fingerprint patterns of varying synchrony levels. In the continuum limit of infinitely many oscillators, the two-core spiral chimeras can be studied using the Ott-Antonsen equation. Numerical analysis of this equation allows us to reveal the stability region of different spiral chimeras, which we group into three main classes-symmetric, asymmetric, and meandering spiral chimeras.