Institut für Physik und Astronomie
Refine
Year of publication
Document Type
- Article (4156)
- Doctoral Thesis (776)
- Postprint (126)
- Monograph/Edited Volume (120)
- Other (86)
- Review (36)
- Habilitation Thesis (24)
- Preprint (24)
- Master's Thesis (10)
- Conference Proceeding (6)
Is part of the Bibliography
- yes (5367) (remove)
Keywords
- diffusion (58)
- stars: massive (57)
- gamma rays: general (47)
- anomalous diffusion (46)
- stars: early-type (45)
- cosmic rays (43)
- stars: winds, outflows (42)
- Magellanic Clouds (39)
- radiation mechanisms: non-thermal (39)
- ISM: supernova remnants (37)
Institute
- Institut für Physik und Astronomie (5367)
- Extern (61)
- Institut für Chemie (18)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (18)
- Institut für Geowissenschaften (10)
- Mathematisch-Naturwissenschaftliche Fakultät (7)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (7)
- Department Psychologie (6)
- Institut für Biochemie und Biologie (4)
- Zentrum für Lehrerbildung und Bildungsforschung (ZeLB) (4)
- Institut für Informatik und Computational Science (2)
- Institut für Mathematik (2)
- Department Erziehungswissenschaft (1)
- Department Linguistik (1)
- Department Sport- und Gesundheitswissenschaften (1)
- Institut für Ernährungswissenschaft (1)
- Institut für Romanistik (1)
- Institut für Umweltwissenschaften und Geographie (1)
- Strukturbereich Bildungswissenschaften (1)
- Wirtschaftswissenschaften (1)
We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGNs), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma-rays from an AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background photons via inverse-Compton scattering to produce gamma-rays. Due to the deflection of the electron- positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218 vertical bar 304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around 10(-14) G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES. 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.
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.
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.
Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E > 100 GeV) gamma-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of similar to 570 hr. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog that are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have spectroscopic distance estimates. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample, we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data set, which shows a 4s excess.
The TeV binary system LS I +61 degrees 303 is known for its regular, non-thermal emission pattern that traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I + 61 degrees. 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I + 61 degrees 303 during the flares, provides constraints on the properties of the accelerator in the source.
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.
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.
The IR-based method of infrared transition moment orientational analysis (IR-TMOA) is employed to unravel molecular order in thin layers of the semiconducting polymer poly[N,N'-bis(2-octyldodecyl),-1,4,5,8-naphthalene-diimide-2,6-diyl]-alt-5-5'-(2,2'-bithiophene) (P(NDI2OD-T2)). Structure-specific vibrational bands are analyzed in dependence On polarization and inclination of the sample-With respect to the optical axis. By that the molecular Order parameter tensor for the respective molecular moieties with regard to the sample: coordinate system is deduced. Making use of the specificity of the IR spectral range, we are able to determine separately the orientation of atomistic planes defined through the naphthalenediimide (NDI) and bithiophene (T2) units relative to the substrate, and hence, relative to each other. A pronounced solvent effect is observed While chlorobenzene causes the T2 planes to align preferentially parallel to the substrate at an angle of 29 degrees, using a 1:1 chloronaphthalene:xylene mixture results in a reorientation of the T2 units from a face on into an edge on arrangement. In contrast the NDI unit remains unaffected. Additionally, for both solvents evidence is observed for the aggregation of chains in accord With recently published results obtained by UV-vis absorption spectroscopy.
We present results of physical experiments where we measure the autocorrelation function (ACF) and the spectral linewidth of the basic frequency of a spiral chaotic attractor in a generator with inertial nonlinearity both without and in the presence of external noise. It is shown that the ACF of spiral attractors decays according to an exponential law with a decrement which is defined by the phase diffusion coefficient. It is also established that the evolution of the instantaneous phase can be approximated by a Wiener random process
We study numerically the behavior of the autocorrelation function (ACF) and the power spectrum of spiral attractors without and in the presence of noise. It is shown that the ACF decays exponentially and has two different time scales. The rate of the ACF decrease is defined by the amplitude fluctuations on small time intervals, i.e., when tau < tau(cor), and by the effective diffusion coefficient of the instantantaneous phase on large time intervals. it is also demonstrated that the ACF in the Poincare map also decreases according to the exponential law exp(-lambda(+)k), where lambda(+) is the positive Lyapunov exponent. The obtained results are compared with the theory of fluctuations for the Van der Pol oscillator
We propose a new autonomous dynamical system of dimension N=4 that demonstrates the regime of stable two- frequency motions and period-doubling bifurcations of a two-dimensional torus. It is shown that the period-doubling bifurcation of the two-dimensional torus is not followed by the resonance phenomenon, and the two-dimensional ergodic torus undergoes a period-doubling bifurcation. The interaction of two generators is also analyzed. The phenomenon of external and mutual synchronization of two-frequency oscillations is observed, for which winding number locking on a two- dimensional torus takes place
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.
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.
Changes of the molecular arrangenemt that can be induced by means of the LB technique in the multilayers of a disc-shaped multialkynyl amphiphile are monitored by means of small angle X-ray diffraction. Studies of the monolayers at the air-water interface reveal "edge-on" orientation of the discs. Specific effects of the counter- ions (Na+, Cd²+, Pb²+, and Ba²+) and sub-solution pH on the monolayer collapse pressure, transfer efficiency and molecular order in the multilayers are found. A correlation between the monolayer properties and the ability for formation of periodic discotic structures in the presence of divalent counterions is established. The discotic molecules retain their "edge-on" arrangement in the highly compressed transferred films with slight irregular interdigtation of the flexible wings and inclination to the substrate normal. The tilt and the inter- digitation are reduced when the discotic monolayers are deposited in alternating LB films with barium arachidate spacer layers.
In this paper a perturbation-theory study of vibrational lifetimes for the bending and stretching modes of hydrogen adsorbed on a Si(100) surface is presented. The hydrogen-silicon interaction is treated with a semiempirical bond-order potential. Calculations are performed for H-Si clusters of different sizes. The finite lifetime is due to vibration-phonon coupling, which is assumed to be linear or bilinear in the phonon and nonlinear in the H-Si stretching and bending modes. Lifetimes and vibrational transition rates are evaluated with one- and two-phonon processes taken into account. Temperature effects are also discussed. In agreement with the experiment and previous theoretical treatment it is found that the H-Si (upsilon(s)=1) stretching vibration decays on a nanosecond timescale, whereas for the H-Si (upsilon(b)=1) bending mode a picosecond decay is predicted. For higher-excited vibrations, simple scaling laws are found if the excitation energies are not too large. The relaxation mechanisms for the excited H-Si stretching and the H-Si bending modes are analyzed in detail.
Recent experiments using time- and angle-resolved two-photon photoemission (2PPE) spectroscopy at metal/polar adsorbate interfaces succeeded in time-dependent analysis of the process of electron solvation. A fully quantum mechanical, two-dimensional simulation of this process, which explicitly includes laser excitation, is presented here, confirming the origin of characteristic features, such as the experimental observation of an apparently negative dispersion. The inference of the spatial extent of the localized electron states from the angular dependence of the 2PPE spectra has been found to be non-trivial and system-dependent. (C) 2005 American Institute of Physics
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.
Vortex ripples in sand are studied experimentally in a one-dimensional setup with periodic boundary conditions. The nonlinear evolution, far from the onset of instability, is analyzed in the framework of a simple model developed for homogeneous patterns. The interaction function describing the mass transport between neighboring ripples is extracted from experimental runs using a recently proposed method for data analysis, and the predictions of the model are compared to the experiment. An analytic explanation of the wavelength selection mechanism in the model is provided, and the width of the stable band of ripples is measured.
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.
We present a catalogue of 362 million stellar parameters, distances, and extinctions derived from Gaia's Early Data Release (EDR3) cross-matched with the photometric catalogues of Pan-STARRS1, SkyMapper, 2MASS, and All WISE. The higher precision of the Gaia EDR3 data, combined with the broad wavelength coverage of the additional photometric surveys and the new stellar-density priors of the StarHorse code, allows us to substantially improve the accuracy and precision over previous photo-astrometric stellar-parameter estimates. At magnitude G = 14 (17), our typical precisions amount to 3% (15%) in distance, 0.13 mag (0.15 mag) in V-band extinction, and 140 K (180 K) in effective temperature. Our results are validated by comparisons with open clusters, as well as with asteroseismic and spectroscopic measurements, indicating systematic errors smaller than the nominal uncertainties for the vast majority of objects. We also provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps, and extensive stellar density maps that reveal detailed substructures in the Milky Way and beyond. The new density maps now probe a much greater volume, extending to regions beyond the Galactic bar and to Local Group galaxies, with a larger total number density. We publish our results through an ADQL query interface (gaia . aip . de) as well as via tables containing approximations of the full posterior distributions. Our multi-wavelength approach and the deep magnitude limit render our results useful also beyond the next Gaia release, DR3.
Galaxies are among the most complex systems that can currently be modelled with a computer. A realistic simulation must take into account cosmology and gravitation as well as effects of plasma, nuclear, and particle physics that occur on very different time, length, and energy scales. The Milky Way is the ideal test bench for such simulations, because we can observe millions of its individual stars whose kinematics and chemical composition are records of the evolution of our Galaxy. Thanks to the advent of multi-object spectroscopic surveys, we can systematically study stellar populations in a much larger volume of the Milky Way. While the wealth of new data will certainly revolutionise our picture of the formation and evolution of our Galaxy and galaxies in general, the big-data era of Galactic astronomy also confronts us with new observational, theoretical, and computational challenges.
This thesis aims at finding new observational constraints to test Milky-Way models, primarily based on infra-red spectroscopy from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and asteroseismic data from the CoRoT mission. We compare our findings with chemical-evolution models and more sophisticated chemodynamical simulations. In particular we use the new powerful technique of combining asteroseismic and spectroscopic observations that allows us to test the time dimension of such models for the first time. With CoRoT and APOGEE (CoRoGEE) we can infer much more precise ages for distant field red-giant stars, opening up a new window for Galactic archaeology.
Another important aspect of this work is the forward-simulation approach that we pursued when interpreting these complex datasets and comparing them to chemodynamical models.
The first part of the thesis contains the first chemodynamical study conducted with the APOGEE survey. Our sample comprises more than 20,000 red-giant stars located within 6 kpc from the Sun, and thus greatly enlarges the Galactic volume covered with high-resolution spectroscopic observations. Because APOGEE is much less affected by interstellar dust extinction, the sample covers the disc regions very close to the Galactic plane that are typically avoided by optical surveys. This allows us to investigate the chemo-kinematic properties of the Milky Way's thin disc outside the solar vicinity. We measure, for the first time with high-resolution data, the radial metallicity gradient of the disc as a function of distance from the Galactic plane, demonstrating that the gradient flattens and even changes its sign for mid-plane distances greater than 1 kpc.
Furthermore, we detect a gap between the high- and low-[$\alpha$/Fe] sequences in the chemical-abundance diagram (associated with the thin and thick disc) that unlike in previous surveys can hardly be explained by selection effects. Using 6D kinematic information, we also present chemical-abundance diagrams cleaned from stars on kinematically hot orbits. The data allow us to confirm without doubt that the scale length of the (chemically-defined) thick disc is significantly shorter than that of the thin disc.
In the second part, we present our results of the first combination of asteroseismic and spectroscopic data in the context of Galactic Archaeology. We analyse APOGEE follow-up observations of 606 solar-like oscillating red giants in two CoRoT fields close to the Galactic plane. These stars cover a large radial range of the Galactic disc (4.5 kpc $\lesssim R_{\rm Gal}\lesssim15$ kpc) and a large age baseline (0.5 Gyr $\lesssim \tau\lesssim$ 13 Gyr), allowing us to study the age- and radius-dependence of the [$\alpha$/Fe] vs. [Fe/H] distributions. We find that the age distribution of the high-[$\alpha$/Fe] sequence appears to be broader than expected from a monolithically-formed old thick disc that stopped to form stars 10 Gyr ago. In particular, we discover a significant population of apparently young, [$\alpha$/Fe]-rich stars in the CoRoGEE data whose existence cannot be explained by standard chemical-evolution models. These peculiar stars are much more abundant in the inner CoRoT field LRc01 than in the outer-disc field LRc01, suggesting that at least part of this population has a chemical-evolution rather than a stellar-evolution origin, possibly due to a peculiar chemical-enrichment history of the inner disc. We also find that strong radial migration is needed to explain the abundance of super-metal-rich stars in the outer disc.
Finally, we use the CoRoGEE sample to study the time evolution of the radial metallicity gradient in the thin disc, an observable that has been the subject of observational and theoretical debate for more than 20 years. By dividing the CoRoGEE dataset into six age bins, performing a careful statistical analysis of the radial [Fe/H], [O/H], and [Mg/Fe] distributions, and accounting for the biases introduced by the observation strategy, we obtain reliable gradient measurements. The slope of the radial [Fe/H] gradient of the young red-giant population ($-0.058\pm0.008$ [stat.] $\pm0.003$ [syst.] dex/kpc) is consistent with recent Cepheid data. For the age range of $1-4$ Gyr, the gradient steepens slightly ($-0.066\pm0.007\pm0.002$ dex/kpc), before flattening again to reach a value of $\sim-0.03$ dex/kpc for stars with ages between 6 and 10 Gyr. This age dependence of the [Fe/H] gradient can be explained by a nearly constant negative [Fe/H] gradient of $\sim-0.07$ dex/kpc in the interstellar medium over the past 10 Gyr, together with stellar heating and migration. Radial migration also offers a new explanation for the puzzling observation that intermediate-age open clusters in the solar vicinity (unlike field stars) tend to have higher metallicities than their younger counterparts. We suggest that non-migrating clusters are more likely to be kinematically disrupted, which creates a bias towards high-metallicity migrators from the inner disc and may even steepen the intermediate-age cluster abundance gradient.
Chemotaxis, the directed motion of a cell toward a chemical source, plays a key role in many essential biological processes. Here, we derive a statistical model that quantitatively describes the chemotactic motion of eukaryotic cells in a chemical gradient. Our model is based on observations of the chemotactic motion of the social ameba Dictyostelium discoideum, a model organism for eukaryotic chemotaxis. A large number of cell trajectories in stationary, linear chemoattractant gradients is measured, using microfluidic tools in combination with automated cell tracking. We describe the directional motion as the interplay between deterministic and stochastic contributions based on a Langevin equation. The functional form of this equation is directly extracted from experimental data by angle-resolved conditional averages. It contains quadratic deterministic damping and multiplicative noise. In the presence of an external gradient, the deterministic part shows a clear angular dependence that takes the form of a force pointing in gradient direction. With increasing gradient steepness, this force passes through a maximum that coincides with maxima in both speed and directionality of the cells. The stochastic part, on the other hand, does not depend on the orientation of the directional cue and remains independent of the gradient magnitude. Numerical simulations of our probabilistic model yield quantitative agreement with the experimental distribution functions. Thus our model captures well the dynamics of chemotactic cells and can serve to quantify differences and similarities of different chemotactic eukaryotes. Finally, on the basis of our model, we can characterize the heterogeneity within a population of chemotactic cells.
The chemotaxis of eukaryotic cells depends both on the average concentration of the chemoattractant and on the steepness of its gradient. For the social amoeba Dictyostelium discoideum, we test quantitatively the prediction by Ueda and Shibata [Biophys. J. 93, 11 (2007)] that the efficacy of chemotaxis depends on a single control parameter only, namely, the signal-to-noise ratio (SNR), determined by the stochastic fluctuations of (i) the binding of the chemoattractant molecule to the transmembrane receptor and (ii) the intracellular activation of the effector of the signaling cascade. For SNR less than or similar to 1, the theory captures the experimental findings well, while for larger SNR noise sources further downstream in the signaling pathway need to be taken into account.
The PAC2MAN mission
(2015)
An accurate forecast of flare and coronal mass ejection (CME) initiation requires precise measurements of the magnetic energy buildup and release in the active regions of the solar atmosphere. We designed a new space weather mission that performs such measurements using new optical instruments based on the Hanle and Zeeman effects. The mission consists of two satellites, one orbiting the L1 Lagrangian point (Spacecraft Earth, SCE) and the second in heliocentric orbit at 1AU trailing the Earth by 80 degrees (Spacecraft 80, SC80). Optical instruments measure the vector magnetic field in multiple layers of the solar atmosphere. The orbits of the spacecraft allow for a continuous imaging of nearly 73% of the total solar surface. In-situ plasma instruments detect solar wind conditions at 1AU and ahead of our planet. Earth-directed CMEs can be tracked using the stereoscopic view of the spacecraft and the strategic placement of the SC80 satellite. Forecasting of geoeffective space weather events is possible thanks to an accurate surveillance of the magnetic energy buildup in the Sun, an optical tracking through the interplanetary space, and in-situ measurements of the near-Earth environment.
Gravity dictates the structure of the whole Universe and, although it is triumphantly described by the theory of General Relativity, it is the force that we least understand in nature. One of the cardinal predictions of this theory are black holes. Massive, dark objects are found in the majority of galaxies. Our own galactic center very contains such an object with a mass of about four million solar masses. Are these objects supermassive black holes (SMBHs), or do we need alternatives? The answer lies in the event horizon, the characteristic that defines a black hole. The key to probe the horizon is to model the movement of stars around a SMBH, and the interactions between them, and look for deviations from real observations. Nuclear star clusters harboring a massive, dark object with a mass of up to ~ ten million solar masses are good testbeds to probe the event horizon of the potential SMBH with stars. The channel for interactions between stars and the central MBH are the fact that (a) compact stars and stellar-mass black holes can gradually inspiral into the SMBH due to the emission of gravitational radiation, which is known as an “Extreme Mass Ratio Inspiral” (EMRI), and (b) stars can produce gases which will be accreted by the SMBH through normal stellar evolution, or by collisions and disruptions brought about by the strong central tidal field. Such processes can contribute significantly to the mass of the SMBH. These two processes involve different disciplines, which combined will provide us with detailed information about the fabric of space and time. In this habilitation I present nine articles of my recent work directly related with these topics.
Simulating the space weather in the AU Mic system: stellar winds and extreme coronal mass ejections
(2022)
Two close-in planets have been recently found around the M-dwarf flare star AU Microscopii (AU Mic). These Neptune-sized planets (AU Mic b and c) seem to be located very close to the so-called "evaporation valley" in the exoplanet population, making this system an important target for studying atmospheric loss on exoplanets. This process, while mainly driven by high-energy stellar radiation, will be strongly mediated by the space environment surrounding the planets. Here we present an investigation of this last area, performing 3D numerical modeling of the quiescent stellar wind from AU Mic, as well as time-dependent simulations describing the evolution of a highly energetic coronal mass ejection (CME) event in this system. Observational constraints on the stellar magnetic field and properties of the eruption are incorporated in our models. We carry out qualitative and quantitative characterizations of the stellar wind, the emerging CMEs, as well as the expected steady and transient conditions along the orbit of both exoplanets. Our results predict extreme space weather for AU Mic and its planets. This includes sub-Alfvenic regions for the large majority of the exoplanet orbits, very high dynamic and magnetic pressure values in quiescence (varying within 10(2)-10(5) times the dynamic pressure experienced by Earth), and an even harsher environment during the passage of any escaping CME associated with the frequent flaring observed in AU Mic. These space weather conditions alone pose an immense challenge for the survival of exoplanetary atmospheres (if any) in this system.
A template-based lamination technique for the manufacture of ferroelectrets from uniform electret films was recently reported. In the present work, this technique is used to prepare similar ferroelectret structures from low-density polyethylene (LDPE) films and from fluoro-ethylene-propylene (FEP) copolymer films. A comparative analysis of the pressure-, temperature-, and frequency-dependent piezoelectric properties has been performed on the two ferroelectret systems. It is observed that the FEP ferroelectrets exhibit better piezoelectric responses and are thermally more stable. The difference between the piezoelectric d(33) coefficients of the two ferroelectret systems is partially explained here by their different elastic moduli. The anti-resonance peaks of both structures have been investigated by means of dielectric resonance spectroscopy and electroacoustic sound-pressure measurements. A difference of more than 10 kHz is observed between the anti-resonance frequencies of the two ferroelectret systems.
We describe the concept, the fabrication, and the most relevant properties of a piezoelectric-polymer system: Two fluoroethylenepropylene (FEP) films with good electret properties are laminated around a specifically designed and prepared polytetrafluoroethylene (PTFE) template at 300 degrees C. After removing the PTFE template, a two-layer FEP film with open tubular channels is obtained. For electric charging, the two-layer FEP system is subjected to a high electric field. The resulting dielectric barrier discharges inside the tubular channels yield a ferroelectret with high piezoelectricity. d(33) coefficients of up to 160 pC/N have already been achieved on the ferroelectret films. After charging at suitable elevated temperatures, the piezoelectricity is stable at temperatures of at least 130 degrees C. Advantages of the transducer films include ease of fabrication at laboratory or industrial scales, a wide range of possible geometrical and processing parameters, straightforward control of the uniformity of the polymer system, flexibility, and versatility of the soft ferroelectrets, and a large potential for device applications e.g., in the areas of biomedicine, communications, production engineering, sensor systems, environmental monitoring, etc.
The acetone extracts of the root bark and stem bark of Erythrina sacleuxii showed antiplasmodial activities against the chloroquine-sensitive (D6) and chloroquine-resistant (W2) strains of Plasmodium falciparum. Chromatographic separation of the acetone extract of the root bark afforded a new isoflavone, 7-hydroxy-4 -methoxy-3'- prenylisoflavone (trivial name 5-deoxy-3' - prenylbiochanin A) along with known isoflavonoids as the antiplasmodial principles. Flavonoids and isoflavonoids isolated from the stem bark of E. sucleuxii were also tested and showed antiplasmodial activities. The structures were determined on the basis of spectroscopic evidence.
Fluoropolymer piezoelectrets with tubular channels resonance behavior controlled by channel geometry
(2012)
Ferro- or piezoelectrets are dielectric materials with two elastically very different macroscopic phases and electrically charged interfaces between them. One of the newer piezoelectret variants is a system of two fluoroethylenepropylene (FEP) films that are first laminated around a polytetrafluoroethylene (PTFE) template. Then, by removing the PTFE template, a two-layer FEP structure with open tubular channels is obtained. After electrical charging, the channels form easily deformable macroscopic electric dipoles whose changes under mechanical or electrical stress lead to significant direct or inverse piezoelectricity, respectively. Here, different PTFE templates are employed to generate channel geometries that vary in height or width. It is shown that the control of the channel geometry allows a direct adjustment of the resonance frequencies in the tubular-channel piezoelectrets. By combining several different channel widths in a single ferroelectret, it is possible to obtain multiple resonance peaks that may lead to a rather flat frequency-response region of the transducer material. A phenomenological relation between the resonance frequency and the geometrical parameters of a tubular channel is also presented. This relation may help to design piezoelectrets with a specific frequency response.
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.
Modeling random crawling, membrane deformation and intracellular polarity of motile amoeboid cells
(2018)
Amoeboid movement is one of the most widespread forms of cell motility that plays a key role in numerous biological contexts. While many aspects of this process are well investigated, the large cell-to-cell variability in the motile characteristics of an otherwise uniform population remains an open question that was largely ignored by previous models. In this article, we present a mathematical model of amoeboid motility that combines noisy bistable kinetics with a dynamic phase field for the cell shape. To capture cell-to-cell variability, we introduce a single parameter for tuning the balance between polarity formation and intracellular noise. We compare numerical simulations of our model to experiments with the social amoeba Dictyostelium discoideum. Despite the simple structure of our model, we found close agreement with the experimental results for the center-of-mass motion as well as for the evolution of the cell shape and the overall intracellular patterns. We thus conjecture that the building blocks of our model capture essential features of amoeboid motility and may serve as a starting point for more detailed descriptions of cell motion in chemical gradients and confined environments.
Modeling random crawling, membrane deformation and intracellular polarity of motile amoeboid cells
(2018)
Amoeboid movement is one of the most widespread forms of cell motility that plays a key role in numerous biological contexts. While many aspects of this process are well investigated, the large cell-to-cell variability in the motile characteristics of an otherwise uniform population remains an open question that was largely ignored by previous models. In this article, we present a mathematical model of amoeboid motility that combines noisy bistable kinetics with a dynamic phase field for the cell shape. To capture cell-to-cell variability, we introduce a single parameter for tuning the balance between polarity formation and intracellular noise. We compare numerical simulations of our model to experiments with the social amoeba Dictyostelium discoideum. Despite the simple structure of our model, we found close agreement with the experimental results for the center-of-mass motion as well as for the evolution of the cell shape and the overall intracellular patterns. We thus conjecture that the building blocks of our model capture essential features of amoeboid motility and may serve as a starting point for more detailed descriptions of cell motion in chemical gradients and confined environments.
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⊙.
We present constraints on the mean dark energy density, Omega(X) and dark energy equation of state parameter, w(X), based on Chandra measurements of the X-ray gas mass fraction in 26 X-ray luminous, dynamically relaxed galaxy clusters spanning the redshift range 0.07 < z < 0.9. Under the assumption that the X-ray gas mass fraction measured within r(2500) is constant with redshift and using only weak priors on the Hubble constant and mean baryon density of the Universe, we obtain a clear detection of the effects of dark energy on the distances to the clusters, confirming (at comparable significance) previous results from Type la supernovae studies. For a standard Lambda cold dark matter (CDM) cosmology with the curvature Omega(K) included as a free parameter, we find Omega(Lambda) = 0.94(-0.23)(+0.21) (68 per cent confidence limits). We also examine extended XCDM dark energy models. Combining the Chandra data with independent constraints from cosmic microwave background experiments, we find Omega(X) = 0.75 +/- 0.04, Omega(m) = 0.26(- 0.04)(+0.06) and w(X) =-1.26 +/- 0.24. Imposing the prior constraint w(X) > -1, the same data require w(X) < -0.7 at 95 per cent confidence. Similar results on the mean matter density and dark energy equation of state parameter, &UOmega;(m) = 0.24 ± 0.04 and w(X) 1.20(-0.28)(+0.24), are obtained by replacing the cosmic microwave background data with standard priors on the Hubble constant and mean baryon density and assuming a flat geometry
We present very-high-energy gamma-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above the background of 20.8 sigma in 47.2 h (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations, the temporal properties of 1ES 2344+514 are studied on short and long times-scales. We fit a constant-flux model to nightly and seasonally binned light curves and apply a fractional variability test to determine the stability of the source on different time-scales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly binned light curves and for the long-term seasonally binned light curve at the > 3 sigma level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (x(2)/NDF = 7.89/6) by a power-law function with an index Gamma = 2.46 +/- 0.06(stat) +/- 0.20(sys) and extends to at least 8 TeV. The extragalactic-backgroundlight-deabsorbed spectrum is adequately fit (x(2)/NDF = 6.73/6) by a power-law function with an index Gamma = 2.15 +/- 0.06(stat) +/- 0.20(sys) while an F-test indicates that the power law with an exponential cut-off function provides a marginally better fit (x(2)/NDF = 2.56/5) at the 2.1 sigma level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source.
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
We present different tests for phase synchronization which improve the procedures currently used in the literature. This is accomplished by using a two-sample 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
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.
We report on the VERITAS discovery of very high energy (VHE) gamma-ray emission above 200 GeV from the high-frequency-peaked BL Lac (HBL) object RX J0648.7+1516 (GB J0648+1516), associated with 1FGL J0648.8+1516. The photon spectrum above 200 GeV is fitted by a power law dN/dE = F-0(E/E-0)(-Gamma) with a photon index Gamma of 4.4 +/- 0.8(stat) +/- 0.3(syst) and a flux normalization F-0 of (2.3 +/- 0.5(stat) +/- 1.2(sys)) x 10(-11) TeV-1 cm(-2) s(-1) with E-0 = 300 GeV. No VHE variability is detected during VERITAS observations of RX J0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identification and spectroscopic redshift were obtained using the Shane 3 m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to subclassify the blazar as an HBL object, including data from the MDM observatory, Swift-UVOT, and X-Ray Telescope, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton (SSC) model overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The spectral energy distribution can be parameterized satisfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.
We present very high energy (VHE) imaging of MGRO J2019+37 obtained with the VERITAS observatory. The bright extended (similar to 2 degrees) unidentified Milagro source is located toward the rich star formation region Cygnus-X. MGRO J2019+37 is resolved into two VERITAS sources. The faint, point-like source VER J2016+371 overlaps CTB 87, a filled-center remnant (SNR) with no evidence of a supernova remnant shell at the present time. Its spectrum is well fit in the 0.65-10 TeV energy range by a power-law model with photon index 2.3 +/- 0.4. VER J2019+378 is a bright extended (similar to 1 degrees) source that likely accounts for the bulk of the Milagro emission and is notably coincident with PSR J2021+3651 and the star formation region Sh 2-104. Its spectrum in the range 1-30 TeV is well fit with a power-law model of photon index 1.75 +/- 0.3, among the hardest values measured in the VHE band, comparable to that observed near Vela-X. We explore the unusual spectrum and morphology in the radio and X-ray bands to constrain possible emission mechanisms for this source.
Prompt emission from the very fluent and nearby (z = 0.34) gamma-ray burst GRB130427A was detected by several orbiting telescopes and by ground-based, wide-field-of-view optical transient monitors. Apart from the intensity and proximity of this GRB, it is exceptional due to the extremely long-lived high-energy (100 MeV to 100 GeV) gamma-ray emission, which was detected by the Large Area Telescope on the Fermi Gamma-Ray Space Telescope for similar to 70 ks after the initial burst. The persistent, hard-spectrum, high-energy emission suggests that the highest-energy gamma rays may have been produced via synchrotron self-Compton processes though there is also evidence that the high-energy emission may instead be an extension of the synchrotron spectrum. VERITAS, a ground-based imaging atmospheric Cherenkov telescope array, began follow-up observations of GRB130427A similar to 71 ks (similar to 20 hr) after the onset of the burst. The GRB was not detected with VERITAS; however, the high elevation of the observations, coupled with the low redshift of the GRB, make VERITAS a very sensitive probe of the emission from GRB130427A for E > 100 GeV. The non-detection and consequent upper limit derived place constraints on the synchrotron self-Compton model of high-energy gamma-ray emission from this burst.
VERITAS has been monitoring the very-high-energy (VHE; > 100 GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +/- 0.16) x 10(-11) cm(-2) s(-1) at energies above 350 GeV. In 2010 April, VERITAS detected a flare from M 87 with peak flux of (2.71 +/- 0.68) x 10(-11) cm(-2) s(-1) for E > 350 GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good-quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90(-0.15)(+0.22) days. The shortest detected exponential rise time is three times as long, at 2.87(+1.65)(-0.99) days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +/- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M 87 and the constraints they provide on the physical size of the emitting region.
We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 giga-electron volts (GeV) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 mega-electron volts and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.
We report on very high energy (>100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the Swift Observatory and later detected by multiple radio, optical, and X-ray observatories. A total exposure of 28 hr was obtained on Swift J164449.3+573451 with the Very Energetic Radiation Imaging Telescope Array System ( VERITAS) during 2011 March 28-April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of 1.4 x 10(-12) erg cm(-2) s(-1) (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.
We present results from VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560 GeV, is well described by a power law with a spectral index of 4.33 +/- 0.09. The time-averaged integral flux above 200 GeV measured for this period was (1.69 +/- 0.06) x 10(-11) photons cm(-2) s(-1), corresponding to 6.9% of the Crab Nebula flux. We also present the combined gamma-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100 MeV to 560 GeV. The data are well fit by a power law with an exponential cutoff at 101.9 +/- 3.2 GeV. The origin of the cutoff could be intrinsic to PG 1553+113 or be due to the gamma-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of z > 0.395 based on optical/UV observations of PG 1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of z <= 0.62. A strongly elevated mean flux of (2.50 +/- 0.14) x10(-11) photons cm(-2) s(-1) (10.3% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as (4.44 +/- 0.71) x10(-11) photons cm(-2) s(-1) (18.3% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a chi(2) probability for a steady flux of 0.03%.
We present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 degrees 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift/XRT (2-10 keV). The data were taken from 2011 December through 2012 January and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hr of quality selected livetime data in which LS I +61 degrees 303 was detected at a statistical significance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 degrees 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 degrees 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.
A search for enhanced very high energy GAMMA-RAY emission from the 2013 march crab nebula flare
(2014)
In 2013 March, a flaring episode from the Crab Nebula lasting similar to 2 weeks was detected by Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). The Very Energetic Radiation Imaging Telescope Array System (VERITAS) provides simultaneous observations throughout this period. During the flare, Fermi-LAT detected a 20 fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very high energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 x 10(-3) with 95% confidence.
We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10(stat) +/- 0.20(sys). The TeV emission is extended, covering the region near PSR J1907+0602 and also extending toward SNR G40.5-0.5. When fitted with a two-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma(src) = 0 degrees.44 +/- 0 degrees.02. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5-0.5, 0 degrees.33 away.
We present the results of a joint observational campaign between the Green Bank radio telescope and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of very-high-energy (VHE) gamma rays (E-gamma > 150 GeV) and giant radio pulses (GRPs) from the Crab pulsar at 8.9 GHz. A total of 15,366 GRPs were recorded during 11.6 hr of simultaneous observations, which were made across four nights in 2008 December and in 2009 November and December. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed around the arrival time of the GRP events. In total, eight different time windows with durations ranging from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Furthermore, we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on the flux simultaneous with interpulse GRPs on single-rotation-period timescales. On similar to 8 s timescales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP emission correlations are below the derived limits.
We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 x 10(-12) erg cm(-2) s(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.
The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant gamma-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are <sigma upsilon >(95%) (CL) less than or similar to 10(-23) cm(3) s(-1), improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of 2 for dark matter particle masses m(chi) greater than or similar to 300 GeV. The lower limits on the decay lifetime are at the level of tau(95%) (CL) greater than or similar to 10(24) s. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.
Multiwavelength observations of the AGN 1ES 0414+009 with veritas, FERMI-LAT, SWIFT-XRT, AND MDM
(2012)
We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between 2008 January and 2011 February, resulting in 56.2 hr of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as (5.2 +/- 1.1(stat) +/- 2.6(sys)) x 10(-12) photons cm(-2) s(-1) above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with a photon index of Gamma = 3.4 +/- 0.5(stat) +/- 0.3(sys) and a flux normalization of (1.6 +/- 0.3(stat) +/- 0.8(sys)) x 10(-11) photons cm(-2) s(-1) at 300 GeV. We also present multiwavelength results taken in the optical (MDM), x-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.
We report on VERITAS very high energy (VHE; E >= 100 GeV) observations of six blazars selected from the Fermi Large Area Telescope First Source Catalog (1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the VHE band, taking gamma-ray absorption by the extragalactic background light into account. This allowed the selection of six bright, hard-spectrum blazars that were good candidate TeV emitters. Spectroscopic redshift measurements were attempted with the Keck Telescope for the targets without Sloan Digital Sky Survey spectroscopic data. No VHE emission is detected during the observations of the six sources described here. Corresponding TeV upper limits are presented, along with contemporaneous Fermi observations and non-concurrent Swift UVOT and X-Ray Telescope data. The blazar broadband spectral energy distributions (SEDs) are assembled and modeled with a single-zone synchrotron self-Compton model. The SED built for each of the six blazars shows a synchrotron peak bordering between the intermediate-and high-spectrum-peak classifications, with four of the six resulting in particle-dominated emission regions.
We report on the discovery of high-energy (HE; E > 0.1 GeV) and very high energy (VHE; E > 100 GeV) gamma-ray emission from the high-frequency-peaked BL Lac object RBS 0413. VERITAS, a ground-based gamma-ray observatory, detected VHE. rays from RBS 0413 with a statistical significance of 5.5 standard deviations (sigma) and a gamma-ray flux of (1.5 +/- 0.6(stat) +/- 0.7(syst)) x 10(-8) photons m(-2) s(-1) (similar to 1% of the Crab Nebula flux) above 250 GeV. The observed spectrum can be described by a power law with a photon index of 3.18 +/- 0.68(stat) +/- 0.30(syst). Contemporaneous observations with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope detected HE gamma rays from RBS 0413 with a (stat)istical significance of more than 9 sigma, a power-law photon index of 1.57 +/- 0.12(stat-0.12sys')(+0.11) and a gamma-ray flux between 300 MeV and 300 GeV of (1.64 +/- 0.43(stat-0.22sys)(+ 0.31)) x 10(-5) photons m(-2) s(-1). We present the results from Fermi-LAT and VERITAS, including a spectral energy distribution modeling of the gamma-ray, quasi-simultaneous X-ray (Swift-XRT), ultraviolet (Swift-UVOT), and R-band optical (MDM) data. We find that, if conditions close to equipartition are required, both the combined synchrotron self-Compton/external-Compton and the lepto-hadronic models are preferred over a pure synchrotron self-Compton model.
We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hr of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of the order of less than or similar to 2 in the high energy (E > 1MeV) and very high energy (E > 100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10 m telescope light curves.
We summarize broadband observations of the TeV-emitting blazar 1ES 1959+650, including optical R-band observations by the robotic telescopes Super-LOTIS and iTelescope, UV observations by Swift Ultraviolet and Optical Telescope, X-ray observations by the Swift X-ray Telescope, high-energy gamma-ray observations with the Fermi Large Area Telescope, and very-high-energy (VHE) gamma-ray observations by VERITAS above 315 GeV, all taken between 2012 April 17 and 2012 June 1 (MJD 56034 and 56079). The contemporaneous variability of the broadband spectral energy distribution is explored in the context of a simple synchrotron self Compton (SSC) model. In the SSC emission scenario, we find that the parameters required to represent the high state are significantly different than those in the low state. Motivated by possible evidence of gas in the vicinity of the blazar, we also investigate a reflected emission model to describe the observed variability pattern. This model assumes that the non-thermal emission from the jet is reflected by a nearby cloud of gas, allowing the reflected emission to re-enter the blob and produce an elevated gamma-ray state with no simultaneous elevated synchrotron flux. The model applied here, although not required to explain the observed variability pattern, represents one possible scenario which can describe the observations. As applied to an elevated VHE state of 66% of the Crab Nebula flux, observed on a single night during the observation period, the reflected emission scenario does not support a purely leptonic non-thermal emission mechanism. The reflected emission model does, however, predict a reflected photon field with sufficient energy to enable elevated gamma-ray emission via pion production with protons of energies between 10 and 100 TeV.
A SEARCH FOR VERY HIGH ENERGY GAMMA RAYS FROM THE MISSING LINK BINARY PULSAR J1023+0038 WITH VERITAS
(2016)
The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than similar to 2 G before the disappearance of the radio pulsar and greater than similar to 10 G afterward.
We present the results of 71.6 hr of observations of the Geminga pulsar (PSR J0633+1746) with the VERITAS very-high-energy gamma-ray telescope array. Data taken with VERITAS between 2007 November and 2013 February were phase-folded using a Geminga pulsar timing solution derived from data recorded by the XMM-Newton and Fermi-LAT space telescopes. No significant pulsed emission above 100 GeV is observed, and we report upper limits at the 95% confidence level on the integral flux above 135 GeV (spectral analysis threshold) of 4.0x10(-13) s(-1) cm(-2) and 1.7 x 10(-13) s(-1) cm(-2) for the two principal peaks in the emission profile. These upper limits, placed in context with phase-resolved spectral energy distributions determined from 5 yr of data from the Fermi-Large Area Telescope (LAT), constrain possible hardening of the Geminga pulsar emission spectra above similar to 50 GeV.
The very high energy (VHE; E > 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54 938-54 956), 2009, as part of an extensive multiwavelength campaign from radio to VHE. Strong VHE yray activity was detected on May 1st with Whipple and VERITAS, when the flux (E > 400 GeV) increased to 10 times the preflare baseline flux (3.9 x 10(-11) ph cm(-2) s(-1)), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15. This VHE flare showed a fast flux variation with an increase of a factor similar to 4 in 25 min, and a falling time of similar to 50 min. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.
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.
Context
Transmission spectroscopy is a promising tool for the atmospheric characterization of transiting exoplanets. Because the planetary signal is faint, discrepancies have been reported regarding individual targets.
Aims
We investigate the dependence of the estimated transmission spectrum on deviations of the orbital parameters of the star-planet system that are due to the limb-darkening effects of the host star. We describe how the uncertainty on the orbital parameters translates into an uncertainty on the planetary spectral slope.
Methods
We created synthetic transit light curves in seven different wavelength bands, from the near-ultraviolet to the near-infrared, and fit them with transit models parameterized by fixed deviating values of the impact parameter b. First, we performed a qualitative study to illustrate the effect by presenting the changes in the transmission spectrum slope with different deviations of b. Then, we quantified these variations by creating an error envelope (for centrally transiting, off-center, and grazing systems) based on a derived typical uncertainty on b from the literature. Finally, we compared the variations in the transmission spectra for different spectral types of host stars.
Results
Our simulations show a wavelength-dependent offset that is more pronounced at the blue wavelengths where the limb-darkening effect is stronger. This offset introduces a slope in the planetary transmission spectrum that becomes steeper with increasing b values. Variations of b by positive or negative values within its uncertainty interval introduce positive or negative slopes, thus the formation of an error envelope. The amplitude from blue optical to near-infrared wavelength for a typical uncertainty on b corresponds to one atmospheric pressure scale height and more. This impact parameter degeneracy is confirmed for different host types; K stars present prominently steeper slopes, while M stars indicate features at the blue wavelengths.
Conclusions
We demonstrate that transmission spectra can be hard to interpret, basically because of the limitations in defining a precise impact parameter value for a transiting exoplanet. This consequently limits a characterization of its atmosphere.
We revisited 10 known exoplanetary systems using publicly available data provided by the transiting exoplanet survey satellite (TESS). The sample presented in this work consists of short period transiting exoplanets, with inflated radii and large reported uncertainty on their planetary radii. The precise determination of these values is crucial in order to develop accurate evolutionary models and understand the inflation mechanisms of these systems. Aiming to evaluate the planetary radius measurement, we made use of the planet-to-star radii ratio, a quantity that can be measured during a transit event. We fit the obtained transit light curves of each target with a detrending model and a transit model. Furthermore, we used emcee, which is based on a Markov chain Monte Carlo approach, to assess the best fit posterior distributions of each system parameter of interest. We refined the planetary radius of WASP-140 b by approximately 12%, and we derived a better precision on its reported asymmetric radius uncertainty by approximately 86 and 67%. We also refined the orbital parameters of WASP-120 b by 2 sigma. Moreover, using the high-cadence TESS datasets, we were able to solve a discrepancy in the literature, regarding the planetary radius of the exoplanet WASP-93 b. For all the other exoplanets in our sample, even though there is a tentative trend that planetary radii of (near-) grazing systems have been slightly overestimated in the literature, the planetary radius estimation and the orbital parameters were confirmed with independent observations from space, showing that TESS and ground-based observations are overall in good agreement.
Planets outside our solar system, so-called "exoplanets", can be detected with different methods, and currently more than 5000 exoplanets have been confirmed, according to NASA Exoplanet Archive. One major highlight of the studies on exoplanets in the past twenty years is the characterization of their atmospheres usingtransmission spectroscopy as the exoplanet transits. However, this characterization is a challenging process and sometimes there are reported discrepancies in the literature regarding the atmosphere of the same exoplanet. One potential reason for the observed atmospheric inconsistencies is called impact parameter degeneracy, and it is highly driven by the limb darkening effect of the host star. A brief introductionto those topics in presented in chapter 1, while the motivation and objectives of thiswork are described in chapter 2.The first goal is to clarify the origin of the transmission spectrum, which is anindicator of an exoplanet’s atmosphere; whether it is real or influenced by the impactparameter degeneracy. A second goal is to determine whether photometry from space using the Transiting Exoplanet Survey Satellite (TESS), could improve on the major parameters, which are responsible for the aforementioned degeneracy, of known exoplanetary systems. Three individual projects were conducted in order toaddress those goals. The three manuscripts are presented, in short, in the manuscriptoverview in chapter 3.More specifically, in chapter 4, the first manuscript is presented, which is an ex-tended investigation on the impact parameter degeneracy and its application onsynthetic transmission spectra. Evidently, the limb darkening of the host star isan important driver for this effect. It keeps the degeneracy persisting through different groups of exoplanets, based on the uncertainty of their impact parameter and on the type of their host star. The second goal, was addressed in the second and third manuscripts (chapter 5 and chapter 6 respectively). Using observationsfrom the TESS mission, two samples of exoplanets were studied; 10 transiting inflated hot-Jupiters and 43 transiting grazing systems. Potentially, the refinement or confirmation of their major system parameters’ measurements can assist in solving current or future discrepancies regarding their atmospheric characterization.In chapter 7 the conclusions of this work are discussed, while in chapter 8 itis proposed how TESS’s measurements can be able to discern between erroneousinterpretations of transmission spectra, especially on systems where the impact parameter degeneracy is likely not applicable.
Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010
(2015)
Context. Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED).
Aims. The goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail.
Methods. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, Fermi-LAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population.
Results. The activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2-10 keV band and the VHE gamma-ray flux above 200 GeV, consistent with the gamma-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE gamma-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 x 10(4) < gamma < 6 x 10(5), where. is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands.
Conclusions. Both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low-and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (delta, B and R), in addition to the parameters related to the particle population.
Aims. We perform an extensive characterization of the broadband emission of Mrk 421, as well as its temporal evolution, during the non-flaring (low) state. The high brightness and nearby location (z = 0.031) of Mrk 421 make it an excellent laboratory to study blazar emission. The goal is to learn about the physical processes responsible for the typical emission of Mrk 421, which might also be extended to other blazars that are located farther away and hence are more difficult to study.
Methods. We performed a 4.5-month multi-instrument campaign on Mrk 421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. This extensive radio to very-high-energy (VHE; E > 100 GeV) gamma-ray dataset provides excellent temporal and energy coverage, which allows detailed studies of the evolution of the broadband spectral energy distribution.
Results. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical /UV and X-rays extending over the duration of the campaign.
Conclusions. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multiwavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk 421 during non-flaring activity. Such a temporally extended X-ray /VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.
Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. However, most of the gamma-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available gamma-ray instrumentation was easier to perform.
Aims. Our goal is to characterize the source gamma-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity.
Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign.
Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density.
Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
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.
We present VR observations of QSO 2237+0305 conducted by the Gravitational Lensing International Time Project collaboration from 1999 October 1 to 2000 February 3. The observations were made with the 2.56 m Nordic Optical Telescope at Roque de los Muchachos Observatory, La Palma (Spain). The point-spread function (PSF) fitting method and an adapted version of the ISIS subtraction method have been used to derive the VR light curves of the four components (A-D) of the quasar. The mean errors range in the intervals 0.01-0.04 mag (PSF fitting) and 0.01-0.02 mag (ISIS subtraction), with the faintest component (D) having the largest uncertainties. We address the relatively good agreement between the A and D light curves derived using different filters, photometric techniques, and telescopes. The new VR light curves of component A extend the time coverage of a high-magnification microlensing peak, which was discovered by the OGLE team.
We calculate the shift in the critical temperature of Bose-Einstein condensation for a dilute Bose-Fermi mixture confined by a harmonic potential, to lowest order in both the Bose-Bose and Bose-Fermi coupling constans. The relativ importance of the effect on the critical temperature of the boson-fermion interactions is investigated as a function of the parameters of mixture. The possible relevance of the shift of the transition temperature in current experiments on trapped Bose-Fermi mixtures is discussed.
We consider a dilute homogeneous mixture of bosons and spin-polarized fermions at zero temperature. We first construct the formal scheme for carrying out systematic perturbation theory in terms of single particle Green's functions. We introduce a new relevant object, the renormalized boson-fermion T-matrix which we determine to second order in the boson-fermion s-wave scattering length. We also discuss how to incorporate the usual boson-boson T-matrix in mean-field approximation to obtain the total ground state properties of the system. The next order term beyond mean- field stems from the boson-fermion interaction and is proportional to $a_{scriptsize BF}k_{scriptsize F}$. The total ground-state energy-density reads $E/V =epsilon_{scriptsize F} + epsilon_{scriptsize B} + (2pihbar^{2}a_{
m BF}n_{scriptsize B}n_{scriptsize F}/m) [1 + a_{scriptsize BF}k_{scriptsize F}f(delta)/pi]$. The first term is the kinetic energy of the free fermions, the second term is the boson-boson mean-field interaction, the pre-factor to the additional term is the usual mean-field contribution to the boson-fermion interaction energy, and the second term in the square brackets is the second-order correction, where $f(delta)$ is a known function of $delta= (m_{scriptsize B} - m_{scriptsize F})/(m_{scriptsize B} + m_{scriptsize F})$. We discuss the relevance of this new term, how it can be incorporated into existing theories of boson-fermion mixtures, and its importance in various parameter regimes, in particular considering mixtures of $^{6}$Li and $^{7}$Li and of $^{3}$He and $^{4}$He.
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.
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.
Recent observations and modeling studies emphasize the crucial role of fracture mechanics for the stability of ice shelves and thereby the evolution of ice sheets. Here we introduce a macroscopic fracture-density field into a prognostic continuum ice-flow model and compute its evolution incorporating the initiation and growth of fractures as well as their advection with two-dimensional ice flow. To a first approximation, fracture growth is assumed to depend on the spreading rate only, while fracture initiation is defined in terms of principal stresses. The inferred fracture-density fields compare well with observed elongate surface structures. Since crevasses and other deep-reaching fracture structures have been shown to influence the overall ice-shelf dynamics, we propose the fracture-density field introduced here be used as a measure for ice softening and decoupling of the ice flow in fracture-weakened zones. This may yield more accurate and realistic velocity patterns in prognostic simulations. Additionally, the memory of past fracture events links the calving front to the upstream dynamics. Thus the fracture-density field proposed here may be employed in fracture-based calving parameterizations. The aim of this study is to introduce the field and investigate which of the observed surface structures can be reproduced by the simplest physically motivated fracture source terms.
Antarctic ice-discharge constitutes the largest uncertainty in future sea-level projections. Floating ice shelves, fringing most of Antarctica, exert retentive forces onto the ice flow. While abrupt ice-shelf retreat has been observed, it is generally considered a localized phenomenon. Here we show that the disintegration of an ice shelf may induce the spontaneous retreat of its neighbor. As an example, we reproduce the spontaneous but gradual retreat of the Larsen B ice front as observed after the disintegration of the adjacent Larsen A ice shelf. We show that the Larsen A collapse yields a change in spreading rate in Larsen B via their connecting ice channels and thereby causes a retreat of the ice front to its observed position of the year 2000, prior to its collapse. This mechanism might be particularly relevant for the role of East Antarctica and the Antarctic Peninsula in future sea level.
Floating ice shelves can exert a retentive and hence stabilizing force onto the inland ice sheet of Antarctica. However, this effect has been observed to diminish by the dynamic effects of fracture processes within the protective ice shelves, leading to accelerated ice flow and hence to a sea-level contribution. In order to account for the macroscopic effect of fracture processes on large-scale viscous ice dynamics (i.e., ice-shelf scale) we apply a continuum representation of fractures and related fracture growth into the prognostic Parallel Ice Sheet Model (PISM) and compare the results to observations. To this end we introduce a higher order accuracy advection scheme for the transport of the two-dimensional fracture density across the regular computational grid. Dynamic coupling of fractures and ice flow is attained by a reduction of effective ice viscosity proportional to the inferred fracture density. This formulation implies the possibility of non-linear threshold behavior due to self-amplified fracturing in shear regions triggered by small variations in the fracture-initiation threshold. As a result of prognostic flow simulations, sharp across-flow velocity gradients appear in fracture-weakened regions. These modeled gradients compare well in magnitude and location with those in observed flow patterns. This model framework is in principle expandable to grounded ice streams and provides simple means of investigating climate-induced effects on fracturing (e. g., hydro fracturing) and hence on the ice flow. It further constitutes a physically sound basis for an enhanced fracture-based calving parameterization.