Refine
Year of publication
Document Type
- Article (4132) (remove)
Keywords
- stars: massive (53)
- gamma rays: general (45)
- stars: early-type (43)
- stars: winds, outflows (41)
- cosmic rays (38)
- diffusion (37)
- radiation mechanisms: non-thermal (37)
- Magellanic Clouds (36)
- X-rays: stars (36)
- ISM: supernova remnants (33)
Institute
- Institut für Physik und Astronomie (4132) (remove)
Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
In recent years, complex network analysis facilitated the identification of universal and unexpected patterns in complex climate systems. However, the analysis and representation of a multiscale complex relationship that exists in the global climate system are limited. A logical first step in addressing this issue is to construct multiple networks over different timescales. Therefore, we propose to apply the wavelet multiscale correlation (WMC) similarity measure, which is a combination of two state-of-the-art methods, viz. wavelet and Pearson’s correlation, for investigating multiscale processes through complex networks. Firstly we decompose the data over different timescales using the wavelet approach and subsequently construct a corresponding network by Pearson’s correlation. The proposed approach is illustrated and tested on two synthetics and one real-world example. The first synthetic case study shows the efficacy of the proposed approach to unravel scale-specific connections, which are often undiscovered at a single scale. The second synthetic case study illustrates that by dividing and constructing a separate network for each time window we can detect significant changes in the signal structure. The real-world example investigates the behavior of the global sea surface temperature (SST) network at different timescales. Intriguingly, we notice that spatial dependent structure in SST evolves temporally. Overall, the proposed measure has an immense potential to provide essential insights on understanding and extending complex multivariate process studies at multiple scales.
We study two coupled spatially extended dynamical systems which exhibit space-time chaos. The transition to the synchronized state is treated as a nonequilibrium phase transition, where the average synchronization error is the order parameter. The transition in one-dimensional systems is found to be generically in the universality class of the Kardar- Parisi-Zhang equation with a growth-limiting term ("bounded KPZ"). For systems with very strong nonlinearities in the local dynamics, however, the transition is found to be in the universality class of directed percolation.
The behavior of the Lyapunov exponents (LEs) of a disordered system consisting of mutually coupled chaotic maps with different parameters is studied. The LEs are demonstrated to exhibit avoided crossing and level repulsion, qualitatively similar to the behavior of energy levels in quantum chaos. Recent results for the coupling dependence of the LEs of two coupled chaotic systems are used to explain the phenomenon and to derive an approximate expression for the distribution functions of LE spacings. The depletion of the level spacing distribution is shown to be exponentially strong at small values. The results are interpreted in terms of the random matrix theory.
Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios. Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature. Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays). Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties.
Aims. We aim to characterize the multiwavelength emission from Markarian 501 (Mrk 501), quantify the energy-dependent variability, study the potential multiband correlations, and describe the temporal evolution of the broadband emission within leptonic theoretical scenarios. Methods. We organized a multiwavelength campaign to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Results. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of similar to 0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was similar to 3 CU, and the peak of the high-energy spectral component was found to be at similar to 2 TeV. Both the X-ray and VHE gamma-ray spectral slopes were measured to be extremely hard, with spectral indices <2 during most of the observing campaign, regardless of the X-ray and VHE flux. This study reports the hardest Mrk 501 VHE spectra measured to date. The fractional variability was found to increase with energy, with the highest variability occurring at VHE. Using the complete data set, we found correlation between the X-ray and VHE bands; however, if the June 9 flare is excluded, the correlation disappears (significance <3 sigma) despite the existence of substantial variability in the X-ray and VHE bands throughout the campaign. Conclusions. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency-peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The data set acquired shows that the broadband spectral energy distribution (SED) of Mrk 501, and its transient evolution, is very complex, requiring, within the framework of synchrotron self-Compton (SSC) models, various emission regions for a satisfactory description. Nevertheless the one-zone SSC scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behavior seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.
Context. The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes above the sensitivity limits of current Cherenkov telescopes have been predicted for both the central X-ray binary system and the interaction regions of SS 433 jets with the surrounding W50 nebula. Non-thermal emission at lower energies has been previously reported, indicating that efficient particle acceleration is taking place in the system. Aims. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses (P-orb similar to 13.1 days) and precession of the circumstellar disk (P-pre similar to 162 days) periodically covering the central binary system is expected to be at its minimum. The eastern and western SS 433/W50 interaction regions are also examined using the whole data set available. We aim to constrain some theoretical models previously developed for this system with our observations. Methods. We made use of dedicated observations from the Major Atmospheric Gamma Imaging Cherenkov telescopes (MAGIC) and High Energy Spectroscopic System (H.E.S.S.) of SS 433 taken from 2006 to 2011. These observation were combined for the first time and accounted for a total effective observation time of 16.5 h, which were scheduled considering the expected phases of minimum absorption of the putative VHE emission. Gamma-ray attenuation does not affect the jet/medium interaction regions. In this case, the analysis of a larger data set amounting to similar to 40-80 h, depending on the region, was employed. Results. No evidence of VHE gamma-ray emission either from the central binary system or from the eastern/western interaction regions was found. Upper limits were computed for the combined data set. Differential fluxes from the central system are found to be less than or similar to 10(-12)-10(-13) TeV-1 cm(-2) s(-1) in an energy interval ranging from similar to few x 100 GeV to similar to few TeV. Integral flux limits down to similar to 10(-12)-10(-13) ph cm(-2) s(-1) and similar to 10(-13)-10(-14) ph cm(-2) s(-1) are obtained at 300 and 800 GeV, respectively. Our results are used to place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Conclusions. Our findings suggest that the fraction of the jet kinetic power that is transferred to relativistic protons must be relatively small in SS 433, q(p) <= 2.5 x 10(-5), to explain the lack of TeV and neutrino emission from the central system. At the SS 433/W50 interface, the presence of magnetic fields greater than or similar to 10 mu G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with E-e up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.
We study localized traveling waves and chaotic states in strongly nonlinear one-dimensional Hamiltonian lattices. We show that the solitary waves are superexponentially localized and present an accurate numerical method allowing one to find them for an arbitrary nonlinearity index. Compactons evolve from rather general initially localized perturbations and collide nearly elastically. Nevertheless, on a long time scale for finite lattices an extensive chaotic state is generally observed. Because of the system's scaling, these dynamical properties are valid for any energy.
Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed.
New amphiphilic compounds 1-9 that feature a construction with dendronized hydrophilic and hydrophobic segment groups connected to a specific aromatic or aliphatic spacer unit have been synthesized, following a modular building block strategy. The hydrophilic dendrons are typically branched elements with peripheral carboxylic groups, unlike the hydrophobic dendrons that contain peripheral alkyl chains as part of respective amide functions. The hydrophilic dendrons are in different generations of branching, while the hydrophobic dendrons are all in the first generation of branching (three terminal branching), but differ in the length of the alkyl chains, thus giving rise to designed structure and amphiphilic properties in the new compounds. The resulting surfactants are capable of forming well-defined Langmuir films of remarkable stability when spread from a solution onto an aqueous subphase. Nevertheless, specific packing behaviour and orientation of the amphiphilic molecules were found, depending on the molecular structure, as determined using analysis of the surface pressure-area (pi-A) isotherms. Langmuir-Blodgett transfer of the first monolayer from a pure water subphase to a clean silicon wafer proved possible for the amphiphiles of peripheral alkyl chain length C-12, while the amphiphiles with the longer alkyl chains failed, possibly due to the more rigid monolayers they form, impeding the transfer.
A comprehensive numerical device simulation of the electrical and optical characteristics accompanied with experimental measurements of a new highly efficient system for polymer-based light-emitting diodes doped with phosphorescent dyes is presented. The system under investigation comprises an electron transporter attached to a polymer backbone blended with an electronically inert small molecule and an iridium-based green phosphorescent dye which serves as both emitter and hole transporter. The device simulation combines an electrical and an optical model. Based on the known highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) levels of all components as well as the measured electrical and optical characteristics of the devices, we model the emissive layer as an effective medium using the dye's HOMO as hole transport level and the polymer LUMO as electron transport level. By fine-tuning the injection barriers at the electron and hole-injecting contact, respectively, in simulated devices, unipolar device characteristics were fitted to the experimental data. Simulations using the so-obtained set of parameters yielded very good agreement to the measured currentvoltage, luminancevoltage characteristics, and the emission profile of entire bipolar light-emitting diodes, without additional fitting parameters. The simulation was used to gain insight into the physical processes and the mechanisms governing the efficiency of the organic light-emitting diode, including the position and extent of the recombination zone, carrier concentration profiles, and field distribution inside the device. The simulations show that the device is severely limited by hole injection, and that a reduction of the hole-injection barrier would improve the device efficiency by almost 50%.
Suppression of the TeV Pair-beam-Plasma Instability by a Tangled Weak Intergalactic Magnetic Field
(2022)
We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kiloparsec, we find that weak magnetic fields increase the transverse momentum of the pair-beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission.
In planetary rings, binary collisions and mutual gravity are the predominant particle interactions. Based on a viscoelastic contact model we implement the concept of static adhesion. We discuss the collision dynamics and obtain a threshold velocity for restitution or agglomeration to occur. The latter takes place within a range of a few cm s(-1) for icy grains at low temperatures. The stability of such two-body agglomerates bound by adhesion and gravity in a tidal environment is discussed and applied to the saturnian system. A maximal agglomerate size for a given orbit location is obtained. In this way we are able to resolve the borderline of the zone where agglomerates can exist as a function of the agglomerate size and thus gain an alternative to the classical Roche limit. An increasing ring grain size with distance to Saturn as observed by the VIMS-experiment on board the Cassini spacecraft can be found by our estimates and implications for the saturnian system will be addressed.
On the basis of a stochastic model of hopping transport in disordered solids we present results of simulations of the dark discharge of surface charged thin films considering the energetic distribution of the localised states within the sheet. A non- linear differential equation with suitable boundary and initial conditions describes the time evolution of the space charge. We suppose a Gaussian distribution of energies for the carrier transporting states with a standard deviation s. The arrival time of carriers at the rear electrode of the sandwich sample is studied in dependence on film thickness, initial surface charge, and energetic disorder. Our calculations confirm that the transit time increases indirect proportional with initial surface charge and proportional with the square thickness. For standard deviations of energetic distribution of 0.05 eV up to 0.25 eV the normalised transit time grows as s1.44. We discuss this in terms of the stochastic transport model.
In this work, the authors present a 7.5% efficient hybrid tandem solar cell with the bottom cell made of amorphous silicon and a Si-PCPDTBT:PC70BM bulk heterojunction top cell. Loss-free recombination contacts were realized by combing Al-doped ZnO with either the conducting polymer composite PEDOT:PSS or with a bilayer of ultrathin Al and MoO3. Optimization of these contacts results in tandem cells with high fill factors of 70% and an open circuit voltage close to the sum of those of the sub-cells. This is the best efficiency reported for this type of hybrid tandem cell so far. Optical and electrical device modeling suggests that the efficiency can be increased to similar to 12% on combining a donor polymer with suitable absorption onset with PCBM. We also describe proof-of-principle studies employing light trapping in hybrid tandem solar cells, suggesting that this device architecture has the potential to achieve efficiencies well above 12%. (C) 2014 Elsevier B.V. All rights reserved.
A novel fluorinated copolymer (F-PCPDTBT) is introduced and shown to exhibit significantly higher power conversion efficiency in bulk heterojunction solar cells with PC70BM compared to the well-known low-band-gap polymer PCPDTBT. Fluorination lowers the polymer HOMO level, resulting in high open-circuit voltages well exceeding 0.7 V. Optical spectroscopy and morphological studies with energy-resolved transmission electron microscopy reveal that the fluorinated polymer aggregates more strongly in pristine and blended layers, with a smaller amount of additives needed to achieve optimum device performance. Time-delayed collection field and charge extraction by linearly increasing voltage are used to gain insight into the effect of fluorination on the field dependence of free charge-carrier generation and recombination. F-PCPDTBT is shown to exhibit a significantly weaker field dependence of free charge-carrier generation combined with an overall larger amount of free charges, meaning that geminate recombination is greatly reduced. Additionally, a 3-fold reduction in non-geminate recombination is measured compared to optimized PCPDTBT blends. As a consequence of reduced non-geminate recombination, the performance of optimized blends of fluorinated PCPDTBT with PC70BM is largely determined by the field dependence of free-carrier generation, and this field dependence is considerably weaker compared to that of blends comprising the non-fluorinated polymer. For these optimized blends, a short-circuit current of 14 mA/cm(2), an open-circuit voltage of 0.74 V, and a fill factor of 58% are achieved, giving a highest energy conversion efficiency of 6.16%. The superior device performance and the low band-gap render this new polymer highly promising for the construction of efficient polymer-based tandem solar cells.
We compare standard and inverted bulk heterojunction solar cells composed of PCPDTBT:PC70BM blends. Inverted devices comprising 100 nm thick active layers exhibited short circuit currents of 15 mA/cm(2), 10% larger than in corresponding standard devices. Modeling of the optical field distribution in the different device stacks proved that this enhancement originates from an increased absorption of incident light in the active layer. Internal quantum efficiencies (IQEs) were obtained from the direct comparison of experimentally derived and modeled currents for different layer thicknesses, yielding IQEs of similar to 70% for a layer thickness of 100 nm. Simulations predict a significant increase of the light harvesting efficiency upon increasing the layer thickness to 270 nm. However, a continuous deterioration of the photovoltaic properties with layer thickness was measured for both device architectures, attributed to incomplete charge extraction. On the other hand, our optical modeling suggests that inverted devices based on PCPDTBT should be able to deliver high power conversion efficiencies (PCEs) of more than 7% provided that recombination losses can be reduced.
We have applied time-delayed collection field (TDCF) and charge extraction by linearly increasing voltage (CELIV) to investigate the photogeneration, transport, and recombination of charge carriers in blends composed of PCPDTBT/PC70BM processed with and without the solvent additive diiodooctane. The results suggest that the solvent additive has severe impacts on the elementary processes involved in the photon to collected electron conversion in these blends. First, a pronounced field dependence of the free carrier generation is found for both blends, where the field dependence is stronger without the additive. Second, the fate of charge carriers in both blends can be described with a rather high bimolecular recombination coefficients, which increase with decreasing internal field. Third, the mobility is three to four times higher with the additive. Both blends show a negative field dependence of mobility, which we suggest to cause bias-dependent recombination coefficients.
We introduce a new and simple method to quantify the effective extraction mobility in organic solar cells at low electric fields and charge carrier densities comparable to operation conditions under one sun illumination. By comparing steady-state carrier densities at constant illumination intensity and under open-circuit conditions, the gradient of the quasi-Fermi potential driving the current is estimated as a function of external bias and charge density. These properties are then related to the respective steady-state current to determine the effective extraction mobility. The new technique is applied to different derivatives of the well-known low-band-gap polymer PCPDTBT blended with PC70BM. We show that the slower average extraction due to lower mobility accounts for the moderate fill factor when solar cells are fabricated with mono- or difluorinated PCPDTBT. This lower extraction competes with improved generation and reduced nongeminate recombination, rendering the monofluorinated derivative the most efficient donor polymer.
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.
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.
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.
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 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.
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.
Professional and amateur astronomers around the world contributed to a 4-month long campaign in 2013, mainly in spectroscopy but also in photometry, interferometry and polarimetry, to observe the first 3 Wolf-Rayet stars discovered: WR 134 (WN6b), WR 135 (WC8) and WR 137 (WC7pd+O9). Each of these stars are interesting in their own way, showing a variety of stellar wind structures. The spectroscopic data from this campaign were reduced and analyzed for WR 134 in order to better understand its behavior and long-term periodicity in the context of CIRs in the wind. We will be presenting the results of these spectroscopic data, which include the confirmation of the CIR variability and a time-coherency of ∼ 40 days (half-life of ∼ 20 days).
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.
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.
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.
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.
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.
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.
Bacterial chemotaxis-a fundamental example of directional navigation in the living world-is key to many biological processes, including the spreading of bacterial infections. Many bacterial species were recently reported to exhibit several distinct swimming modes-the flagella may, for example, push the cell body or wrap around it. How do the different run modes shape the chemotaxis strategy of a multimode swimmer? Here, we investigate chemotactic motion of the soil bacterium Pseudomonas putida as a model organism. By simultaneously tracking the position of the cell body and the configuration of its flagella, we demonstrate that individual run modes show different chemotactic responses in nutrition gradients and, thus, constitute distinct behavioral states. On the basis of an active particle model, we demonstrate that switching between multiple run states that differ in their speed and responsiveness provides the basis for robust and efficient chemotaxis in complex natural habitats.
The 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.
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.
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 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.
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 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 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.
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.
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.
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.
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 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.
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 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.
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 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 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.
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.
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.
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.
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.
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
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.
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
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⊙.
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.
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.
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 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.
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.
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.
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.
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.
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.
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.
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.
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.
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
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.
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.
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.
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
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 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
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.
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.
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.