Refine
Year of publication
- 2017 (276) (remove)
Document Type
- Article (215)
- Doctoral Thesis (22)
- Other (17)
- Postprint (14)
- Review (6)
- Habilitation Thesis (2)
Language
- English (276) (remove)
Is part of the Bibliography
- yes (276)
Keywords
- stars: massive (11)
- stars: early-type (8)
- stars: winds, outflows (7)
- Magellanic Clouds (6)
- stars: mass-loss (6)
- X-rays: stars (5)
- acceleration of particles (5)
- cosmic rays (5)
- stars: magnetic field (5)
- X-rays: binaries (4)
Institute
- Institut für Physik und Astronomie (276) (remove)
Studying the temporal variability of BL Lac objects at the highest energies provides unique insights into the extreme physical processes occurring in relativistic jets and in the vicinity of super-massive black holes. To this end, the long-term variability of the BL Lac object PKS 2155 304 is analyzed in the high (HE, 100MeV < E < 300 GeV) and very high energy (VHE, E > 200 GeV) gamma-ray domain. Over the course of similar to 9 yr of H. E. S. S. observations the VHE light curve in the quiescent state is consistent with a log-normal behavior. The VHE variability in this state is well described by flicker noise (power-spectral-density index beta(VHE) = 1 .10(+ 0 : 10) (0 : 13)) on timescales larger than one day. An analysis of similar to 5.5 yr of HE Fermi-LAT data gives consistent results (beta(HE) = 1 : 20(+ 0 : 21) (0 : 23), on timescales larger than 10 days) compatible with the VHE findings. The HE and VHE power spectral densities show a scale invariance across the probed time ranges. A direct linear correlation between the VHE and HE fluxes could neither be excluded nor firmly established. These long-term-variability properties are discussed and compared to the red noise behavior (beta similar to 2) seen on shorter timescales during VHE-flaring states. The difference in power spectral noise behavior at VHE energies during quiescent and flaring states provides evidence that these states are influenced by different physical processes, while the compatibility of the HE and VHE long-term results is suggestive of a common physical link as it might be introduced by an underlying jet-disk connection.
Measurement of the EBL spectral energy distribution using the VHE gamma-ray spectra of HESS blazars
(2017)
Very high-energy gamma rays (VHE, E greater than or similar to 100 GeV) propagating over cosmological distances can interact with the low-energy photons of the extragalactic background light (EBL) and produce electron-positron pairs. The transparency of the Universe to VHE gamma rays is then directly related to the spectral energy distribution (SED) of the EBL. The observation of features in the VHE energy spectra of extragalactic sources allows the EBL to be measured, which otherwise is very difficult. An EBL model-independent measurement of the EBL SED with the H.E.S.S. array of Cherenkov telescopes is presented. It was obtained by extracting the EBL absorption signal from the reanalysis of high-quality spectra of blazars. From H.E.S.S. data alone the EBL signature is detected at a significance of 9.5 sigma, and the intensity of the EBL obtained in different spectral bands is presented together with the associated gamma-ray horizon.
We search for high-energy gamma-ray emission from the binary neutron star merger GW170817 with the H.E.S.S. Imaging Air Cherenkov Telescopes. The observations presented here have been obtained starting only 5.3 hr after GW170817. The H.E.S.S. target selection identified regions of high probability to find a counterpart of the gravitational-wave event. The first of these regions contained the counterpart SSS17a that has been identified in the optical range several hours after our observations. We can therefore present the first data obtained by a ground-based pointing instrument on this object. A subsequent monitoring campaign with the H.E.S.S. telescopes extended over several days, covering timescales from 0.22 to 5.2 days and energy ranges between 270 GeV to 8.55 TeV. No significant gamma-ray emission has been found. The derived upper limits on the very-high-energy gamma-ray flux for the first time constrain non-thermal, high-energy emission following the merger of a confirmed binary neutron star system.
Discovery of Very-high-energy Emission from RGB J2243+203 and Derivation of Its Redshift Upper Limit
(2017)
Very-high-energy (VHE; > 100 GeV) gamma-ray emission from the blazar RGB J2243+203 was discovered with the VERITAS Cherenkov telescope array, during the period between 2014 December 21 and 24. The VERITAS energy spectrum from this source can be fitted by a power law with a photon index of 4.6 +/- 0.5, and a flux normalization at 0.15 TeV of (6.3 +/- 1.1) x 10(-10) cm(-2) s(-1) TeV-1. The integrated Fermi-LAT flux from 1 to 100 GeV during the VERITAS detection is (4.1 +/- 0.8) x 10(-8) cm(-2) s(-1), which is an order of magnitude larger than the four-year-averaged flux in the same energy range reported in the 3FGL catalog, (4.0 +/- 0.1 x 10(-9) cm(-2) s(-1)). The detection with VERITAS triggered observations in the X-ray band with the Swift-XRT. However, due to scheduling constraints Swift-XRT observations were performed 67 hr after the VERITAS detection, rather than simultaneously with the VERITAS observations. The observed X-ray energy spectrum between 2 and 10 keV can be fitted with a power law with a spectral index of 2.7 +/- 0.2, and the integrated photon flux in the same energy band is (3.6 +/- 0.6) x 10(-13) cm(-2) s(-1). EBL-model-dependent upper limits of the blazar redshift have been derived. Depending on the EBL model used, the upper limit varies in the range from z < 0.9 to z < 1.1.
B2 1215+30 is a BL-Lac-type blazar that was first detected at TeV energies by the MAGIC atmospheric Cherenkov telescopes and subsequently confirmed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) observatory with data collected between 2009 and 2012. In 2014 February 08, VERITAS detected a large-amplitude flare from B2. 1215+30 during routine monitoring observations of the blazar 1ES. 1218+304, located in the same field of view. The TeV flux reached 2.4 times the Crab Nebula flux with a variability timescale of <3.6 hr. Multiwavelength observations with Fermi-LAT, Swift, and the Tuorla Observatory revealed a correlated high GeV flux state and no significant optical counterpart to the flare, with a spectral energy distribution where the gamma-ray luminosity exceeds the synchrotron luminosity. When interpreted in the framework of a onezone leptonic model, the observed emission implies a high degree of beaming, with Doppler factor delta > 10, and an electron population with spectral index p < 2.3.
Prospects for Cherenkov Telescope Array Observations of the Young Supernova Remnant RX J1713.7-3946
(2017)
We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX J1713.7-3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very high energy (VHE) gamma rays. Special attention is paid to exploring possible spatial (anti) correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H I emission. We present a series of simulated images of RX J1713.7-3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the nonthermal X-ray emission observed by XMM-Newton, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H I observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic versus leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.
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.
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.
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⊙.
Galaxies are among the most complex systems that can currently be modelled with a computer. A realistic simulation must take into account cosmology and gravitation as well as effects of plasma, nuclear, and particle physics that occur on very different time, length, and energy scales. The Milky Way is the ideal test bench for such simulations, because we can observe millions of its individual stars whose kinematics and chemical composition are records of the evolution of our Galaxy. Thanks to the advent of multi-object spectroscopic surveys, we can systematically study stellar populations in a much larger volume of the Milky Way. While the wealth of new data will certainly revolutionise our picture of the formation and evolution of our Galaxy and galaxies in general, the big-data era of Galactic astronomy also confronts us with new observational, theoretical, and computational challenges.
This thesis aims at finding new observational constraints to test Milky-Way models, primarily based on infra-red spectroscopy from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and asteroseismic data from the CoRoT mission. We compare our findings with chemical-evolution models and more sophisticated chemodynamical simulations. In particular we use the new powerful technique of combining asteroseismic and spectroscopic observations that allows us to test the time dimension of such models for the first time. With CoRoT and APOGEE (CoRoGEE) we can infer much more precise ages for distant field red-giant stars, opening up a new window for Galactic archaeology.
Another important aspect of this work is the forward-simulation approach that we pursued when interpreting these complex datasets and comparing them to chemodynamical models.
The first part of the thesis contains the first chemodynamical study conducted with the APOGEE survey. Our sample comprises more than 20,000 red-giant stars located within 6 kpc from the Sun, and thus greatly enlarges the Galactic volume covered with high-resolution spectroscopic observations. Because APOGEE is much less affected by interstellar dust extinction, the sample covers the disc regions very close to the Galactic plane that are typically avoided by optical surveys. This allows us to investigate the chemo-kinematic properties of the Milky Way's thin disc outside the solar vicinity. We measure, for the first time with high-resolution data, the radial metallicity gradient of the disc as a function of distance from the Galactic plane, demonstrating that the gradient flattens and even changes its sign for mid-plane distances greater than 1 kpc.
Furthermore, we detect a gap between the high- and low-[$\alpha$/Fe] sequences in the chemical-abundance diagram (associated with the thin and thick disc) that unlike in previous surveys can hardly be explained by selection effects. Using 6D kinematic information, we also present chemical-abundance diagrams cleaned from stars on kinematically hot orbits. The data allow us to confirm without doubt that the scale length of the (chemically-defined) thick disc is significantly shorter than that of the thin disc.
In the second part, we present our results of the first combination of asteroseismic and spectroscopic data in the context of Galactic Archaeology. We analyse APOGEE follow-up observations of 606 solar-like oscillating red giants in two CoRoT fields close to the Galactic plane. These stars cover a large radial range of the Galactic disc (4.5 kpc $\lesssim R_{\rm Gal}\lesssim15$ kpc) and a large age baseline (0.5 Gyr $\lesssim \tau\lesssim$ 13 Gyr), allowing us to study the age- and radius-dependence of the [$\alpha$/Fe] vs. [Fe/H] distributions. We find that the age distribution of the high-[$\alpha$/Fe] sequence appears to be broader than expected from a monolithically-formed old thick disc that stopped to form stars 10 Gyr ago. In particular, we discover a significant population of apparently young, [$\alpha$/Fe]-rich stars in the CoRoGEE data whose existence cannot be explained by standard chemical-evolution models. These peculiar stars are much more abundant in the inner CoRoT field LRc01 than in the outer-disc field LRc01, suggesting that at least part of this population has a chemical-evolution rather than a stellar-evolution origin, possibly due to a peculiar chemical-enrichment history of the inner disc. We also find that strong radial migration is needed to explain the abundance of super-metal-rich stars in the outer disc.
Finally, we use the CoRoGEE sample to study the time evolution of the radial metallicity gradient in the thin disc, an observable that has been the subject of observational and theoretical debate for more than 20 years. By dividing the CoRoGEE dataset into six age bins, performing a careful statistical analysis of the radial [Fe/H], [O/H], and [Mg/Fe] distributions, and accounting for the biases introduced by the observation strategy, we obtain reliable gradient measurements. The slope of the radial [Fe/H] gradient of the young red-giant population ($-0.058\pm0.008$ [stat.] $\pm0.003$ [syst.] dex/kpc) is consistent with recent Cepheid data. For the age range of $1-4$ Gyr, the gradient steepens slightly ($-0.066\pm0.007\pm0.002$ dex/kpc), before flattening again to reach a value of $\sim-0.03$ dex/kpc for stars with ages between 6 and 10 Gyr. This age dependence of the [Fe/H] gradient can be explained by a nearly constant negative [Fe/H] gradient of $\sim-0.07$ dex/kpc in the interstellar medium over the past 10 Gyr, together with stellar heating and migration. Radial migration also offers a new explanation for the puzzling observation that intermediate-age open clusters in the solar vicinity (unlike field stars) tend to have higher metallicities than their younger counterparts. We suggest that non-migrating clusters are more likely to be kinematically disrupted, which creates a bias towards high-metallicity migrators from the inner disc and may even steepen the intermediate-age cluster abundance gradient.
HESS J1826-130
(2017)
HESS J1826-130 is an unidentified hard spectrum source discovered by H.E.S.S. along the Galactic plane, the spectral index being Gamma = 1.6 with an exponential cut-off at about 12 TeV. While the source does not have a clear counterpart at longer wavelengths, the very hard spectrum emission at TeV energies implies that electrons or protons accelerated up to several hundreds of TeV are responsible for the emission. In the hadronic case, the VHE emission can be produced by runaway cosmic-rays colliding with the dense molecular clouds spatially coincident with the H.E.S.S. source.
Chemotactic motion in a chemical gradient is an essential cellular function that controls many processes in the living world. For a better understanding and more detailed modelling of the underlying mechanisms of chemotaxis, quantitative investigations in controlled environments are needed. We developed a setup that allows us to separately address the dependencies of the chemotactic motion on the average background concentration and on the gradient steepness of the chemoattractant. In particular, both the background concentration and the gradient steepness can be kept constant at the position of the cell while it moves along in the gradient direction. This is achieved by generating a well-defined chemoattractant gradient using flow photolysis. In this approach, the chemoattractant is released by a light-induced reaction from a caged precursor in a microfluidic flow chamber upstream of the cell. The flow photolysis approach is combined with an automated real-time cell tracker that determines changes in the cell position and triggers movement of the microscope stage such that the cell motion is compensated and the cell remains at the same position in the gradient profile. The gradient profile can be either determined experimentally using a caged fluorescent dye or may be alternatively determined by numerical solutions of the corresponding physical model. To demonstrate the function of this adaptive microfluidic gradient generator, we compare the chemotactic motion of Dictyostelium discoideum cells in a static gradient and in a gradient that adapts to the position of the moving cell. Published by AIP Publishing.
Dark matter constraints from a joint analysis of dwarf Spheroidal galaxy observations with VERITAS
(2017)
We present constraints on the annihilation cross section of weakly interacting massive particles dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of similar to 230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is 1.35 x 10(-23) cm(3) s(-1) at 1 TeV for the bottom quark (b (b) over bar) final state, 2.85 x 10(-24) cm(3) s(-1) at 1 TeV for the tau lepton (tau+tau(-)) final state and 1.32 x 10-25 cm(3) s(-1) at 1 TeV for the gauge boson (gamma gamma) final state.
Imaging atmospheric Cherenkov telescopes (IACTs) are equipped with sensitive photomultiplier tube (PMT) cameras. Exposure to high levels of background illumination degrades the efficiency of and potentially destroys these photo-detectors over time, so IACTs cannot be operated in the same configuration in the presence of bright moonlight as under dark skies. Since September 2012, observations have been carried out with the VERITAS IACTs under bright moonlight (defined as about three times the night-sky-background (NSB) of a dark extragalactic field, typically occurring when Moon illumination > 35%) in two observing modes, firstly by reducing the voltage applied to the PMTs and, secondly, with the addition of ultra-violet (UV) bandpass filters to the cameras. This has allowed observations at up to about 30 times previous NSB levels (around 80% Moon illumination), resulting in 30% more observing time between the two modes over the course of a year. These additional observations have already allowed for the detection of a flare from the 1ES 1727 + 502 and for an observing program targeting a measurement of the cosmic-ray positron fraction. We provide details of these new observing modes and their performance relative to the standard VERITAS observations. (C) 2017 Elsevier B.V. All rights reserved.
We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGNs), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when multi-TeV gamma-rays from an AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background photons via inverse-Compton scattering to produce gamma-rays. Due to the deflection of the electron- positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218 vertical bar 304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around 10(-14) G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES. 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.
Optical control of magnetization using femtosecond laser without applying any external magnetic field offers the advantage of switching magnetic states at ultrashort time scales. Recently, all-optical helicity-dependent switching (AO-HDS) has drawn a significant attention for potential information and data storage device applications. In this work, we employ element and magnetization sensitive photoemission electron microscopy (PEEM) to investigate the role of heating in AO-HDS for thin films of the rare-earth transition-metal alloy TbFe. Spatially resolved measurements in a 3–5 μm sized stationary laser spot demonstrate that AO-HDS is a local phenomenon in the vicinity of thermal demagnetization in a ‘ring’ shaped region. The efficiency of AO-HDS further depends on a local temperature profile around the demagnetized region and thermally activated domain wall motion. We also demonstrate that the thickness of the film determines the preferential switching direction for a particular helicity.
Up until recently, signatures of the ultrarelativistic electron loss driven by electromagnetic ion cyclotron (EMIC) waves in the Earth's outer radiation belt have been limited to direct or indirect measurements of electron precipitation or the narrowing of normalized pitch angle distributions in the heart of the belt. In this study, we demonstrate additional observational evidence of ultrarelativistic electron loss that can be driven by resonant interaction with EMIC waves. We analyzed the profiles derived from Van Allen Probe particle data as a function of time and three adiabatic invariants between 9 October and 29 November 2012. New local minimums in the profiles are accompanied by the narrowing of normalized pitch angle distributions and ground‐based detection of EMIC waves. Such a correlation may be indicative of ultrarelativistic electron precipitation into the Earth's atmosphere caused by resonance with EMIC waves.
The work done during the PhD studies has been focused on measurements of distribution functions of rotating galaxies using integral field spectroscopy observations.
Throughout the main body of research presented here we have been using CALIFA (Calar Alto Legacy Integral Field Area) survey stellar velocity fields to obtain robust measurements of circular velocities for rotating galaxies of all morphological types. A crucial part of the work was enabled by well-defined CALIFA sample selection criteria: it enabled reconstructing sample-independent distributions of galaxy properties.
In Chapter 2, we measure the distribution in absolute magnitude - circular velocity space for a well-defined sample of 199 rotating CALIFA galaxies using their stellar kinematics. Our aim in this analysis is to avoid subjective selection criteria and to take volume and large-scale structure factors into account. Using stellar velocity fields instead of gas emission line kinematics allows including rapidly rotating early type galaxies. Our initial sample contains 277 galaxies with available stellar velocity fields and growth curve r-band photometry. After rejecting 51 velocity fields that could not be modelled due to the low number of bins, foreground contamination or significant interaction we perform Markov Chain Monte Carlo (MCMC) modelling of the velocity fields, obtaining the rotation curve and kinematic parameters and their realistic uncertainties. We perform an extinction correction and calculate the circular velocity v_circ accounting for pressure support a given galaxy has. The resulting galaxy distribution on the M_r - v_circ plane is then modelled as a mixture of two distinct populations, allowing robust and reproducible rejection of outliers, a significant fraction of which are slow rotators. The selection effects are understood well enough that the incompleteness of the sample can be corrected and the 199 galaxies can be weighted by volume and large-scale structure factors enabling us to fit a volume-corrected Tully-Fisher relation (TFR). More importantly, we also provide the volume-corrected distribution of galaxies in the M_r - v_circ plane, which can be compared with cosmological simulations. The joint distribution of the luminosity and circular velocity space densities, representative over the range of -20 > M_r > -22 mag, can place more stringent constraints on the galaxy formation and evolution scenarios than linear TFR fit parameters or the luminosity function alone.
In Chapter 3, we measure one of the marginal distributions of the M_r - v_circ distribution: the circular velocity function of rotating galaxies. The velocity function is a fundamental observable statistic of the galaxy population, being of a similar importance as the luminosity function, but much more difficult to measure. We present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km s^-1 for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population 60 < v_circ < 170 km s^-1, we use the HIPASS velocity function. For the massive galaxy population 170 < v_circ < 320 km s^-1, we use stellar circular velocities from CALIFA. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies. It has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic Illustris simulations fare better, but still do not fully reproduce observations.
In Chapter 4, we present some other work done during the PhD studies, namely, a method that improves the precision of specific angular measurements by combining simultaneous Markov Chain Monte Carlo modelling of ionised gas 2D velocity fields and HI linewidths. To test the method we use a sample of 25 galaxies from the Sydney-AAO Multi-object Integral field (SAMI) survey that had matching ALFALFA HI linewidths. Such a method allows constraining the rotation curve both in the inner regions of a galaxy and in its outskirts, leading to increased precision of specific angular momentum measurements. It could be used to further constrain the observed relation between galaxy mass, specific angular momentum and morphology (Obreschkow & Glazebrook 2014).
Mathematical and computational methods are presented in the appendices.
Organic solar cells demonstrate external quantum efficiencies and fill factors approaching those of conventional photovoltaic technologies. However, as compared with the optical gap of the absorber materials, their open-circuit voltage is much lower, largely due to the presence of significant non-radiative recombination. Here, we study a large data set of published and new material combinations and find that non-radiative voltage losses decrease with increasing charge-transfer-state energies. This observation is explained by considering non-radiative charge-transfer-state decay as electron transfer in the Marcus inverted regime, being facilitated by a common skeletal molecular vibrational mode. Our results suggest an intrinsic link between non-radiative voltage losses and electron-vibration coupling, indicating that these losses are unavoidable. Accordingly, the theoretical upper limit for the power conversion efficiency of single-junction organic solar cells would be reduced to about 25.5% and the optimal optical gap increases to (1.45-1.65) eV, that is, (0.2-0.3) eV higher than for technologies with minimized non-radiative voltage losses.
Dynamic processes in living cells are highly organized in space and time. Unraveling the underlying molecular mechanisms of spatiotemporal pattern formation remains one of the outstanding challenges at the interface between physics and biology. A fundamental recurrent pattern found in many different cell types is that of self-sustained oscillations. They are involved in a wide range of cellular functions, including second messenger signaling, gene expression, and cytoskeletal dynamics. Here, we review recent developments in the field of cellular oscillations and focus on cases where concepts from physics have been instrumental for understanding the underlying mechanisms. We consider biochemical and genetic oscillators as well as oscillations that arise from chemo-mechanical coupling. Finally, we highlight recent studies of intracellular waves that have increasingly moved into the focus of this research field.
We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit the extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta beta p =5.10(-4) and 0.5 beta p =), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane, while keeping it at 90 degrees to the shock normal. Ion reflection off of the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all of the cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field, the fraction of suprathermal electrons is much higher than for other configurations, because only in this case are the Buneman modes fully captured by the 2D simulation grid. Shocks in plasma with moderate bp provide more efficient pre-acceleration. The relevance of our results to the physics of fully 3D systems is discussed.
Chimera states consisting of synchronous and asynchronous domains in a medium of nonlinearly coupled phase oscillators have been considered. Stationary inhomogeneous solutions of the Ott-Antonsen equation for a complex order parameter that correspond to fundamental chimeras have been constructed. The direct numerical simulation has shown that these structures under certain conditions are transformed to oscillatory (breathing) chimera regimes because of the development of instability.
Objective: The aims of this systematic review and meta-analysis are (i) to determine the impact of school-based interventions on objectively measured physical activity among adolescents and (ii) to examine accelerometer methods and decision rule reporting in previous interventions. Methods: A systematic search was performed to identify randomized controlled trials targeting adolescents (age: >= 10 years), conducted in the school setting, and reporting objectively measured physical activity. Random effects meta-analyses were conducted to determine the pooled effects of previous interventions on total and moderate-to-vigorous physical activity. Potential moderators of intervention effects were also explored. Results: Thirteen articles met the inclusion criteria, and twelve were included in the meta-analysis. The pooled effects were small and non-significant for both total physical activity (standardized mean difference = 0.02 [95% confidence interval = -0.13 to 0.18]) and moderate-to-vigorous physical activity (standardized mean difference = 0.24 [95% confidence interval = -0.08 to 0.56]). Sample age and accelerometer compliance were significant moderators for total physical activity, with a younger sample and higher compliance associated with larger effects. Conclusion: Previous school-based physical activity interventions targeting adolescents have been largely unsuccessful, particularly for older adolescents. There is a need for more high-quality research using objective monitoring in this population. Future interventions should comply with best-practice recommendations regarding physical activity monitoring protocols.
In a recent paper [U. Lorenz and P. Saalfrank, Chem. Phys. 482, 69 (2017)], we proposed a robust scheme to set up a system-bath model Hamiltonian, describing the coupling of adsorbate vibrations (system) to surface phonons (bath), from first principles. The method is based on an embedded cluster approach, using orthogonal coordinates for system and bath modes, and an anharmonic phononic expansion of the system-bath interaction up to second order. In this contribution, we use this model Hamiltonian to calculate vibrational relaxation rates of H–Si and D–Si bending modes, coupled to a fully H(D)-covered Si(100)-(2×1) surface, at zero temperature. The D–Si bending mode has an anharmonic frequency lying inside the bath frequency spectrum, whereas the H–Si bending mode frequency is outside the bath Debye band. Therefore, in the present calculations, we only take into account one-phonon system-bath couplings for the D–Si system and both one- and two-phonon interaction terms in the case of H–Si. The computation of vibrational lifetimes is performed with two different approaches, namely, Fermi’s golden rule, and a generalized Bixon-Jortner model built in a restricted vibrational space of the adsorbate-surface zeroth-order Hamiltonian. For D–Si, the Bixon-Jortner Hamiltonian can be solved by exact diagonalization, serving as a benchmark, whereas for H–Si, an iterative scheme based on the recursive residue generation method is applied, with excellent convergence properties. We found that the lifetimes obtained with perturbation theory, albeit having almost the same order of magnitude—a few hundred fs for D–Si and a couple of ps for H–Si—, are strongly dependent on the discretized numerical representation of the bath spectral density. On the other hand, the Bixon-Jortner model is free of such numerical deficiencies, therefore providing better estimates of vibrational relaxation rates, at a very low computational cost. The results obtained with this model clearly show a net exponential decay of the time-dependent survival probability for the H–Si initial vibrational state, allowing an easy extraction of the bending mode “lifetime.” This is in contrast with the D–Si system, whose survival probability exhibits a non-monotonic decay, making it difficult to define such a lifetime. This different behavior of the vibrational decay is rationalized in terms of the power spectrum of the adsorbate-surface system. In the case of D–Si, it consists of several, non-uniformly distributed peaks around the bending mode frequency, whereas the H–Si spectrum exhibits a single Lorentzian lineshape, whose width corresponds to the calculated lifetime. The present work gives some insight into mechanisms of vibration-phonon coupling at surfaces. It also serves as a benchmark for multidimensional system-bath quantum dynamics, for comparison with approximate schemes such as reduced, open-system density matrix theory (where the bath is traced out and a Liouville-von Neumann equation is solved) or approximate wavefunction methods to solve the combined system-bath Schrödinger equation.
Context. Supergiant fast X-ray (SFXT) transients are a peculiar class of supergiant X-ray binaries characterized by a remarkable variability in the X-ray domain, widely ascribed to accretion from a clumpy stellar wind. Aims. In this paper we performed a systematic and homogeneous analysis of the sufficiently bright X-ray flares observed with XMM-Newton from the supergiant fast X-ray transients to probe spectral variations on timescales as short as a few hundred seconds. Our ultimate goal is to investigate whether SFXT flares and outbursts are triggered by the presence of clumps, and to reveal whether strongly or mildly dense clumps are required. Methods. For all sources, we employ a technique developed by our group already exploited in a number of our previous papers, making use of an adaptive rebinned hardness ratio to optimally select the time intervals for the spectral extraction. A total of twelve observations performed in the direction of five SFXTs are reported, providing the largest sample of events available so far. Results. Using the original results reported here and those obtained with our technique from the analysis of two previously published XMM-Newton observations of IGR J17544-2619 and IGR J18410-0535, we show that both strongly and mildly dense clumps can trigger these events. In the former case, the local absorption column density may increase by a factor of >> 3, while in the latter case, the increase is only a factor of similar to 2-3 (or lower). An increase in the absorption column density is generally recorded during the rise of the flares/outbursts, while a drop follows when the source achieves peak flux. In a few cases, a re-increase of the absorption column density after the flare is also detected, and we discovered one absorption event related to the passage of an unaccreted clump in front of the compact object. Overall, there seems to be no obvious correlation between the dynamic ranges in the X-ray fluxes and absorption column densities in supergiant fast X-ray transients, with an indication that lower densities are recorded at the highest fluxes. Conclusions. The spectral variations measured in all sources are in agreement with the idea that the flares/outbursts are triggered by the presence of dense structures in the wind interacting with the X-rays from the compact object (leading to photoionization). The lack of correlation between the dynamic ranges in the X-ray fluxes and absorption column densities can be explained by the presence of accretion inhibition mechanism(s). Based on the knowledge acquired so far on the SFXTs, we propose a classification of the flares/outbursts from these sources in order to drive future observational investigations. We suggest that the difference between the classes of flares/outbursts is related to the fact that the mechanism(s) inhibiting accretion can be overcome more easily in some sources compared to others. We also investigate the possibility that different stellar wind structures, other than clumps, could provide the means to temporarily overcome the inhibition of accretion in supergiant fast X-ray transients.
The long-term X-ray light curves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of three-dimensional (3D) hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ 16493-4348, which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light curve are similar to the observed ones, thus providing support to our scenario. We propose that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.
The evolution of the crystal structure and crystallographic texture of porous synthetic cordierite was studied by in situ high-temperature neutron diffraction up to 1373 K, providing the first in situ high-temperature texture measurement of this technologically important material. It was observed that the crystal texture slightly weakens with increasing temperature, concurrently with subtle changes in the crystal structure. These changes are in agreement with previous work, leading the authors to the conclusion that high-temperature neutron diffraction allows reliable crystallographic characterization of materials with moderate texture. It was also observed that structural changes occur at about the glass transition temperature of the cordierite glass (between 973 and 1073 K). Crystal structure refinements were conducted with and without quantitative texture analysis being part of the Rietveld refinement, and a critical comparison of the results is presented, contributing to the sparse body of literature on combined texture and crystal structure refinements.
Over the last decades, the percentage of the age group choosing to pursue university studies has increased significantly across the world. At the same time, there are university teachers who believe that the standards have fallen. There is little research on whether students nowadays demonstrate knowledge or abilities similar to that of the preceding cohorts. However, in times of educational expansion, empirical evidence on student test performance is extremely helpful in evaluating how well educational systems cope with the increasing numbers of students. In this study, we compared a sample of 2322 physics freshmen from 2013 with another sample of 2718 physics freshmen from 1978 at universities in Germany with regard to their physics knowledge based on their results in the same entrance test. Previous results on mathematics knowledge and abilities in the same sample of students indicated that there was no severe decline in their average achievement. This paper compares the physics knowledge of the same two samples of students. Contrary to their mathematics results, their physics results showed a substantial decrease in physics knowledge as measured by the test.
Stress-induced damage evolution in cast AlSi12CuMgNi alloy with one- and two-ceramic reinforcements
(2017)
Two composites, consisting of an as-cast AlSi12CuMgNi alloy reinforced with 15 vol% Al2O3 short fibres and with 7 vol% Al2O3 short fibres + 15 vol% SiC particles, were studied. Synchrotron computed tomography disclosed distribution, orientation, and volume fraction of the different phases. In-situ compression tests during neutron diffraction in direction parallel to the fibres plane revealed the load partition between phases. Internal damage (fragmentation) of the Si phase and Al2O3 fibres was directly observed in CT reconstructions. Significant debonding between Al matrix and SiC particles was also found. Finally, based on the Maxwell scheme, a micromechanical model was utilized for the new composite with two-ceramic reinforcements; it rationalizes the experimental data and predicts the evolution of all internal stress components in each phase.
How different are the properties of critical adsorption of polyampholytes and polyelectrolytes onto charged surfaces? How important are the details of polyampholyte charge distribution on the onset of critical adsorption transition? What are the scaling relations governing the dependence of critical surface charge density on salt concentration in the surrounding solution? Here, we employ Metropolis Monte Carlo simulations and uncover the scaling relations for critical adsorption for quenched periodic and random charge distributions along the polyampholyte chains. We also evaluate and discuss the dependence of the adsorbed layer width on solution salinity and details of the charge distribution. We contrast our findings to the known results for polyelectrolyte adsorption onto oppositely charged surfaces, in particular, their dependence on electrolyte concentration.
Electron flux in the Earth’s outer radiation belt is highly variable due to a delicate balance between competing acceleration and loss processes. It has been long recognized that Electromagnetic Ion Cyclotron (EMIC) waves may play a crucial role in the loss of radiation belt electrons. Previous theoretical studies proposed that EMIC waves may account for the loss of the relativistic electron population. However, recent observations showed that while EMIC waves are responsible for the significant loss of ultra-relativistic electrons, the relativistic electron population is almost unaffected. In this study, we provide a theoretical explanation for this discrepancy between previous theoretical studies and recent observations. We demonstrate that EMIC waves mainly contribute to the loss of ultra-relativistic electrons. This study significantly improves the current understanding of the electron dynamics in the Earth’s radiation belt and also can help us understand the radiation environments of the exoplanets and outer planets.
Metal halide perovskites have emerged as materials of high interest for solar energy-to-electricity conversion, and in particular, the use of mixed-ion structures has led to high power conversion efficiencies and improved stability. For this reason, it is important to develop means to obtain atomic level understanding of the photoinduced behavior of these materials including processes such as photoinduced phase separation and ion migration. In this paper, we implement a new methodology combining visible laser illumination of a mixed-ion perovskite ((FAP-bI(3))(0.85)(MAPbBr(3))(0.15)) with the element specificity and chemical sensitivity of core-level photoelectron spectroscopy. By carrying out measurements at a synchrotron beamline optimized for low X-ray fluxes, we are able to avoid sample changes due to X-ray illumination and are therefore able to monitor what sample changes are induced by visible illumination only. We find that laser illumination causes partially reversible chemistry in the surface region, including enrichment of bromide at the surface, which could be related to a phase separation into bromide- and iodide-rich phases. We also observe a partially reversible formation of metallic lead in the perovskite structure. These processes occur on the time scale of minutes during illumination. The presented methodology has a large potential for understanding light-induced chemistry in photoactive materials and could specifically be extended to systematically study the impact of morphology and composition on the photostability of metal halide perovskites.
Luminous and high-frequency peaked blazars: the origin of the gamma-ray emission from PKS 1424+240
(2017)
Context. The current generation of ground-based Cherenkov telescopes, together with the LAT instrument on-board the Fermi satellite, have greatly increased our knowledge of gamma-ray blazars. Among them, the high-frequency-peaked BL Lacertae object (HBL) PKS 1424+240 (z similar or equal to 0.6) is the farthest persistent emitter of very-high-energy (VHE; E >= 100 GeV) gamma-ray photons. Current emission models can satisfactorily reproduce typical blazar emission assuming that the dominant emission process is synchrotron-self-Compton (SSC) in HBLs; and external-inverse-Compton (EIC) in low-frequency-peaked BL Lacertae objects and flat-spectrum-radio-quasars. Alternatively, hadronic models are also able to correctly reproduce the gamma-ray emission from blazars, although they are in general disfavored for bright quasars and rapid flares. Aims. The blazar PKS 1424+240 is a rare example of a luminous HBL, and we aim to determine which is the emission process most likely responsible for its gamma-ray emission. This will impact more generally our comprehension of blazar emission models, and how they are related to the luminosity of the source and the peak frequency of the spectral energy distribution. Methods. We have investigated different blazar emission models applied to the spectral energy distribution of PKS 1424+240. Among leptonic models, we study a one-zone SSC model (including a systematic study of the parameter space), a two-zone SSC model, and an EIC model. We then investigated a blazar hadronic model, and finally a scenario in which the gamma-ray emission is associated with cascades in the line-of-sight produced by cosmic rays from the source. Results. After a systematic study of the parameter space of the one-zone SSC model, we conclude that this scenario is not compatible with gamma-ray observations of PKS 1424+240. A two-zone SSC scenario can alleviate this issue, as well as an EIC solution. For the latter, the external photon field is assumed to be the infra-red radiation from the dusty torus, otherwise the VHE gamma-ray emission would have been significantly absorbed. Alternatively, hadronic models can satisfactorily reproduce the gamma-ray emission from PKS 1424+240, both as in-source emission and as cascade emission.
We present an approach for reconstructing networks of pulse-coupled neuronlike oscillators from passive observation of pulse trains of all nodes. It is assumed that units are described by their phase response curves and that their phases are instantaneously reset by incoming pulses. Using an iterative procedure, we recover the properties of all nodes, namely their phase response curves and natural frequencies, as well as strengths of all directed connections.
In ageing systems physical observables explicitly depend on the time span elapsing between the original initiation of the system and the actual start of the recording of the particle motion. We here study the signatures of ageing in the framework of ultraslow continuous time random walk processes with super-heavy tailed waiting time densities. We derive the density for the forward or recurrent waiting time of the motion as function of the ageing time, generalise the Montroll-Weiss equation for this process, and analyse the ageing behaviour of the ensemble and time averaged mean squared displacements.
A growing number of biological, soft, and active matter systems are observed to exhibit normal diffusive dynamics with a linear growth of the mean-squared displacement, yet with a non-Gaussian distribution of increments. Based on the Chubinsky-Slater idea of a diffusing diffusivity, we here establish and analyze a minimal model framework of diffusion processes with fluctuating diffusivity. In particular, we demonstrate the equivalence of the diffusing diffusivity process with a superstatistical approach with a distribution of diffusivities, at times shorter than the diffusivity correlation time. At longer times, a crossover to a Gaussian distribution with an effective diffusivity emerges. Specifically, we establish a subordination picture of Brownian but non-Gaussian diffusion processes, which can be used for a wide class of diffusivity fluctuation statistics. Our results are shown to be in excellent agreement with simulations and numerical evaluations.
Cost-efficient, visible-light-driven hydrogen production from water is an attractive potential source of clean, sustainable fuel. Here, it is shown that thermal solid state reactions of traditional carbon nitride precursors (cyanamide, melamine) with NaCl, KCl, or CsCl are a cheap and straightforward way to prepare poly(heptazine imide) alkali metal salts, whose thermodynamic stability decreases upon the increase of the metal atom size. The chemical structure of the prepared salts is confirmed by the results of X-ray photoelectron and infrared spectroscopies, powder X-ray diffraction and electron microscopy studies, and, in the case of sodium poly(heptazine imide), additionally by atomic pair distribution function analysis and 2D powder X-ray diffraction pattern simulations. In contrast, reactions with LiCl yield thermodynamically stable poly(triazine imides). Owing to the metastability and high structural order, the obtained heptazine imide salts are found to be highly active photo-catalysts in Rhodamine B and 4-chlorophenol degradation, and Pt-assisted sacrificial water reduction reactions under visible light irradiation. The measured hydrogen evolution rates are up to four times higher than those provided by a benchmark photocatalyst, mesoporous graphitic carbon nitride. Moreover, the products are able to photocatalytically reduce water with considerable reaction rates, even when glycerol is used as a sacrificial hole scavenger.
What are the features of partitioning of crystalline materials on the surface of a two-component icosahedral vesicle? We model the response of the rigid hardly stretchable crystalline icosahedra upon addition of a softer component on its surface. We demonstrate how the soft phase "invades" the shell regions with the highest elastic energy density around 12 5-fold topological defects. We explore the phase diagram of these inhomogeneous shells as a function of the soft material fraction, shell radius, and elastic moduli of the two phases. The findings are compared with the recent computer simulation findings, and their biological relevance, for example, for the structure of icosahedral viruses, is also discussed.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black-Scholes-Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.
In this Letter, we study the role of the donor:acceptor interface nanostructure upon charge separation and recombination in organic photovoltaic devices and blend films, using mixtures of PBTTT and two different fullerene derivatives (PC70BM and ICTA) as models for intercalated and nonintercalated morphologies, respectively. Thermodynamic simulations show that while the completely intercalated system exhibits a large free-energy barrier for charge separation, this barrier is significantly lower in the nonintercalated system and almost vanishes when energetic disorder is included in the model. Despite these differences, both femtosecond-resolved transient absorption spectroscopy (TAS) and time-delayed collection field (TDCF) exhibit extensive first-order losses in both systems, suggesting that geminate pairs are the primary product of photoexcitation. In contrast, the system that comprises a combination of fully intercalated polymer:fullerene areas and fullerene-aggregated domains (1:4 PBTTT:PC70BM) is the only one that shows slow, second-order recombination of free charges, resulting in devices with an overall higher short-circuit current and fill factor. This study therefore provides a novel consideration of the role of the interfacial nanostructure and the nature of bound charges and their impact upon charge generation and recombination.
A constitutive model for the nonlinear or "pseudoplastic" mechanical behavior in a linear-elastic solid with thermally induced microcracks is developed and applied to experimental results. The model is termed strain dependent microcrack density approximation (SDMDA) and is an extension of the modified differential scheme that describes the slope of the stress-strain curves of microcracked solids. SDMDA allows a continuous variation in the microcrack density with tensile loading. Experimental uniaxial tensile response of beta-eucryptite glass and ceramics with controlled levels of microcracking is reported. It is demonstrated that SDMDA can well describe the extent of non-linearity in the experimental uniaxial tensile response of beta-eucryptite with varying levels of microcracking. The advantages of the SDMDA are discussed in regard to tensile loading.
The dynamics of fragmentation and vibration of molecular systems with a large number of coupled degrees of freedom are key aspects for understanding chemical reactivity and properties. Here we present a resonant inelastic X-ray scattering (RIXS) study to show how it is possible to break down such a complex multidimensional problem into elementary components. Local multimode nuclear wave packets created by X-ray excitation to different core-excited potential energy surfaces (PESs) will act as spatial gates to selectively probe the particular ground-state vibrational modes and, hence, the PES along these modes. We demonstrate this principle by combining ultra-high resolution RIXS measurements for gas-phase water with state-of-the-art simulations.
The dynamics of fragmentation and vibration of molecular systems with a large number of coupled degrees of freedom are key aspects for understanding chemical reactivity and properties. Here we present a resonant inelastic X-ray scattering (RIXS) study to show how it is possible to break down such a complex multidimensional problem into elementary components. Local multimode nuclear wave packets created by X-ray excitation to different core-excited potential energy surfaces (PESs) will act as spatial gates to selectively probe the particular ground-state vibrational modes and, hence, the PES along these modes. We demonstrate this principle by combining ultra-high resolution RIXS measurements for gas-phase water with state-of-the-art simulations.
Arctic Amplification of climate warming is caused by various feedback processes in the atmosphere-ocean-ice system and yields the strongest temperature increase during winter in the Arctic North Atlantic region. In our study, we attempt to quantify the advective contribution to the observed atmospheric warming in the Svalbard area. Based on radiosonde measurements from Ny-Ålesund, a strong dependence of the tropospheric temperature on the synoptic flow direction is revealed. Using FLEXTRA backward trajectories, an increase of advection from the lower latitude Atlantic region towards Ny-Ålesund is found that is attributed to a change in atmospheric circulation patterns. We find that about one-quarter (0.45 K per decade) of the observed atmospheric winter near surface warming trend in the North Atlantic region of the Arctic (2 K per decade) is due to increased advection of warm and moist air from the lower latitude Atlantic region, affecting the entire troposphere.
We study triggering processes in triaxial compression experiments under a constant displacement rate on sandstone and granite samples using spatially located acoustic emission events and their focal mechanisms. We present strong evidence that event-event triggering plays an important role in the presence of large-scale or macrocopic imperfections, while such triggering is basically absent if no significant imperfections are present. In the former case, we recover all established empirical relations of aftershock seismicity including the Gutenberg-Richter relation, a modified version of the Omori-Utsu relation and the productivity relation-despite the fact that the activity is dominated by compaction-type events and triggering cascades have a swarmlike topology. For the Gutenberg-Richter relations, we find that the b value is smaller for triggered events compared to background events. Moreover, we show that triggered acoustic emission events have a focal mechanism much more similar to their associated trigger than expected by chance.
e-ASTROGAM (‘enhanced ASTROGAM’) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV – the lower energy limit can be pushed to energies as low as 150 keV, albeit with rapidly degrading angular resolution, for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and the promise of eLISA.
Massive runaway stars produce bow shocks through the interaction of their winds with the interstellar medium, with the prospect for particle acceleration by the shocks. These objects are consequently candidates for non-thermal emission. Our aim is to investigate the X-ray emission from these sources. We observed with XMM-Newton a sample of five bow shock runaways, which constitutes a significant improvement of the sample of bow shock runaways studied in X-rays so far. A careful analysis of the data did not reveal any X-ray emission related to the bow shocks. However, X-ray emission from the stars is detected, in agreement with the expected thermal emission from stellar winds. On the basis of background measurements we derive conservative upper limits between 0.3 and 10 keV on the bow shocks emission. Using a simple radiation model, these limits together with radio upper limits allow us to constrain some of the main physical quantities involved in the non-thermal emission processes, such as the magnetic field strength and the amount of incident infrared photons. The reasons likely responsible for the non-detection of non-thermal radiation are discussed. Finally, using energy budget arguments, we investigate the detectability of inverse Compton X-rays in a more extended sample of catalogued runaway star bow shocks. From our analysis we conclude that a clear identification of non-thermal X-rays from massive runaway bow shocks requires one order of magnitude (or higher) sensitivity improvement with respect to present observatories.
We present an effective dynamical model for the onset of bacterial bioluminescence, one of the most studied quorum sensing-mediated traits. Our model is built upon simple equations that describe the growth of the bacterial colony, the production and accumulation of autoinducer signal molecules, their sensing within bacterial cells, and the ensuing quorum activation mechanism that triggers bioluminescent emission. The model is directly tested to quantitatively reproduce the experimental distributions of photon emission times, previously measured for bacterial colonies of Vibrio jasicida, a luminescent bacterium belonging to the Harveyi clade, growing in a highly drying environment. A distinctive and novel feature of the proposed model is bioluminescence ‘quenching’ after a given time elapsed from activation. Using an advanced fitting procedure based on the simulated annealing algorithm, we are able to infer from the experimental observations the biochemical parameters used in the model. Such parameters are in good agreement with the literature data. As a further result, we find that, at least in our experimental conditions, light emission in bioluminescent bacteria appears to originate from a subtle balance between colony growth and quorum activation due to autoinducers diffusion, with the two phenomena occurring on the same time scale. This finding is consistent with a negative feedback mechanism previously reported for Vibrio harveyi.
We study the effect of common noise on coupled active rotators. While such a noise always facilitates synchrony, coupling may be attractive (synchronizing) or repulsive (desynchronizing). We develop an analytical approach based on a transformation to approximate angle-action variables and averaging over fast rotations. For identical rotators, we describe a transition from full to partial synchrony at a critical value of repulsive coupling. For nonidentical rotators, the most nontrivial effect occurs at moderate repulsive coupling, where a juxtaposition of phase locking with frequency repulsion (anti-entrainment) is observed. We show that the frequency repulsion obeys a nontrivial power law.
Electromagnetic ion cyclotron (EMIC) waves play an important role in the dynamics of ultrarelativistic electron population in the radiation belts. However, as EMIC waves are very sporadic, developing a parameterization of such wave properties is a challenging task. Currently, there are no dynamic, activity-dependent models of EMIC waves that can be used in the long-term (several months) simulations, which makes the quantitative modeling of the radiation belt dynamics incomplete. In this study, we investigate Kp, Dst, and AE indices, solar wind speed, and dynamic pressure as possible parameters of EMIC wave presence. The EMIC waves are included in the long-term simulations (1year, including different geomagnetic activity) performed with the Versatile Electron Radiation Belt code, and we compare results of the simulation with the Van Allen Probes observations. The comparison shows that modeling with EMIC waves, parameterized by solar wind dynamic pressure, provides a better agreement with the observations among considered parameterizations. The simulation with EMIC waves improves the dynamics of ultrarelativistic fluxes and reproduces the formation of the local minimum in the phase space density profiles.
The diffusion in the comb structures is a popular model of geometrically induced anomalous diffusion. In the present work we concentrate on the diffusion along the backbone in a system where sidebranches are planes, and the diffusion thereon is anomalous and described by continuous time random walks (CTRW). We show that the mean squared displacement (MSD) in the backbone of the comb behaves differently depending on whether the waiting time periods in the sidebranches are reset after the step in the backbone is done (a rejuvenating junction model), or not (a non-rejuvenating junction model). In the rejuvenating case the subdiffusion in the sidebranches only changes the prefactor in the ultra-slow (logarithmic) diffusion along the backbone, while in the non-rejuvenating case the ultraslow, logarithmic subdiffusion is changed to a much faster power-law subdiffusion (with a logarithmic correction) as it was found earlier by Iomin and Mendez [25]. Moreover, in the first case the result does not change if the diffusion in the backbone is itself anomalous, while in the second case it does. Two of the special cases of the considered models (the non-rejuvenating junction under normal diffusion in the backbone, and rejuvenating junction for the same waiting time distribution in the sidebranches and in junction points) were also investigated within the approach based on the corresponding generalized Fokker-Planck equations. (c) 2017 Elsevier Ltd. All rights reserved.
The valence orbitals of aqueous histidine under basic, neutral and acidic conditions and their X-ray induced transformations have been monitored through N 1s resonant inelastic X-ray scattering. Using density functional ab initio molecular dynamics simulations in the core-hole state within the Z + 1 approximation, core-excitation-induced molecular transformations are quantified. Spectroscopic evidence for a highly directional X-ray-induced local N-H dissociation within the scattering duration is presented for acidic histidine. Our report demonstrates a protonation-state and chemical-environment dependent propensity for a molecular dissociation, which is induced by the absorption of high energy photons. This case study indicates that structural deformations in biomolecules under exposure to ionizing radiation, yielding possible alteration or loss of function, is highly dependent on the physiological state of the molecule upon irradiation.
The femtosecond excited-state dynamics following resonant photoexcitation enable the selective deformation of N-H and N-C chemical bonds in 2-thiopyridone in aqueous solution with optical or X-ray pulses. In combination with multiconfigurational quantum-chemical calculations, the orbital-specific electronic structure and its ultrafast dynamics accessed with resonant inelastic X-ray scattering at the N 1s level using synchrotron radiation and the soft X-ray free-electron laser LCLS provide direct evidence for this controlled photoinduced molecular deformation and its ultrashort time-scale.
The femtosecond excited-state dynamics following resonant photoexcitation enable the selective deformation of N-H and N-C chemical bonds in 2-thiopyridone in aqueous solution with optical or X-ray pulses. In combination with multiconfigurational quantum-chemical calculations, the orbital-specific electronic structure and its ultrafast dynamics accessed with resonant inelastic X-ray scattering at the N 1s level using synchrotron radiation and the soft X-ray free-electron laser LCLS provide direct evidence for this controlled photoinduced molecular deformation and its ultrashort time-scale.
Shape change is a fundamental process occurring in biological tissues during embryonic development and regeneration of tissues and organs. This process is regulated by cells that are constrained within a complex environment of biochemical and physical cues. The spatial constraint due to geometry has a determining role on tissue mechanics and the spatial distribution of force patterns that, in turn, influences the organization of the tissue structure. An understanding of the underlying principles of tissue organization may have wide consequences for the understanding of healing processes and the development of organs and, as such, is of fundamental interest for the tissue engineering community.
This thesis aims to further our understanding of how the collective behaviour of cells is influenced by the 3D geometry of the environment. Previous research studying the role of geometry on tissue growth has mainly focused either on flat surfaces or on substrates where at least one of the principal curvatures is zero. In the present work, tissue growth from MC3T3-E1 pre-osteoblasts was investigated on surfaces of controlled mean curvature.
One key aspect of this thesis was the development of substrates of controlled mean curvature and their visualization in 3D. It was demonstrated that substrates of controlled mean curvature suitable for cell culture can be fabricated using liquid polymers and surface tension effects.
Using these substrates, it was shown that the mean surface curvature has a strong impact on the rate of tissue growth and on the organization of the tissue structure. It was thereby not only demonstrated that the amount of tissue produced (i.e. growth rates) by the cells depends on the mean curvature of the substrate but also that the tissue surface behaves like a viscous fluid with an equilibrium shape governed by the Laplace-Young-law. It was observed that more tissue was formed on highly concave surfaces compared to flat or convex surfaces.
Motivated by these observations, an analytical model was developed, where the rate of tissue growth is a function of the mean curvature, which could successfully describe the growth kinetics. This model was also able to reproduce the growth kinetics of previous experiments where tissues have been cultured in straight-sided prismatic pores.
A second part of this thesis focuses on the tissue structure, which influences the mechanical properties of the mature bone tissue. Since the extracellular matrix is produced by the cells, the cell orientation has a strong impact on the direction of the tissue fibres. In addition, it was recently shown that some cell types exhibit collective alignment similar to liquid crystals.
Based on this observation, a computational model of self-propelled active particles was developed to explore in an abstract manner how the collective behaviour of cells is influenced by 3D curvature. It was demonstrated that the 3D curvature has a strong impact on the self-organization of active particles and gives, therefore, first insights into the principles of self-organization of cells on curved surfaces.
This study introduces a thermally responsive, polymer-based electrode system. The key component is a surface-attached, temperature-responsive poly(oligoethylene glycol) methacrylate (poly(OEGMA)) type polymer bearing photoreactive benzophenone and carboxy groups containing side chains. The responsive behavior of the polymer in aqueous media has been investigated by turbidimetry measurements. Polymer films are formed on gold substrates by means of the photoreactive 2(dicyclohexylphosphino)benzophenone (DPBP) through photocrosslinking. The electrochemical behavior of the resulting polymer-substrate interface has been investigated in buffered [Fe(CN)6](3-)/[Fe (CN)6](4-)solutions at room temperature and under temperature variation by cyclic voltammetry (CV). The CV experiments show that with increasing temperature structural changes of the polymer layer occur, which alter the output of the electrochemical measurement. Repeated heating/cooling cycles analyzed by CV measurements and pH changes analyzed by quartz crystal microbalance with dissipation monitoring (QCM-D) reveal the reversible nature of the restructuring process. The immobilized films are further modified by covalent coupling of two small biomolecules - a hydrophobic peptide and a more hydrophilic one. These attached components influence the hydrophobicity of the layer in a different way the resulting change of the temperature-caused behavior has been studied by CV indicating a different state of the polymer after coupling of the hydrophobic peptide.
During muscle contractions, radial-force distributions are generated on muscle surfaces due to muscle-volume changes, from which the corresponding body motions can be recorded by means of so-called force myography (FMG). Piezo- or ferroelectrets are flexible piezoelectric materials with attractive materials and sensing properties. In addition to several other applications, they are suitable for detecting force variations by means of wearable devices. In this paper, we prepared piezoelectrets from cellular polypropylene films by optimizing the fabrication procedures, and developed an FMG-recording system based on piezoelectret sensors. Different hand and wrist movements were successfully detected on able-bodied subjects with the FMG system. The FMG patterns were evaluated and identified by means of linear discriminant analysis and artificial neural network algorithms, and average motion-classification accuracies of 96.1% and 94.8%, respectively, were obtained. This paper demonstrates the feasibility of using piezoelectret-film sensors for FMG and may thus lead to alternative methods for detecting body motion and to related applications, e.g., in biomedical engineering or structural-health monitoring.
Here we report on a cyclic, physical ice-discharge instability in the Parallel Ice Sheet Model, simulating the flow of a three-dimensional, inherently buttressed ice-sheet-shelf system which periodically surges on a millennial timescale. The thermomechanically coupled model on 1 km horizontal resolution includes an enthalpy-based formulation of the thermodynamics, a nonlinear stress-balance-based sliding law and a very simple subglacial hydrology. The simulated unforced surging is characterized by rapid ice streaming through a bed trough, resulting in abrupt discharge of ice across the grounding line which is eventually calved into the ocean. We visualize the central feedbacks that dominate the subsequent phases of ice buildup, surge and stabilization which emerge from the interaction between ice dynamics, thermodynamics and the subglacial till layer. Results from the variation of surface mass balance and basal roughness suggest that ice sheets of medium thickness may be more susceptible to surging than relatively thin or thick ones for which the surge feedback loop is damped. We also investigate the influence of different basal sliding laws (ranging from purely plastic to nonlinear to linear) on possible surging. The presented mechanisms underlying our simulations of self-maintained, periodic ice growth and destabilization may play a role in large-scale ice-sheet surging, such as the surging of the Laurentide Ice Sheet, which is associated with Heinrich events, and ice-stream shutdown and reactivation, such as observed in the Siple Coast region of West Antarctica.
Bifurcations of dynamos in rotating and buoyancy-driven spherical Rayleigh-Benard convection in an electrically conducting fluid are investigated numerically. Both nonmagnetic and magnetic solution branches comprised of rotating waves are traced by path-following techniques, and their bifurcations and interconnections for different Ekman numbers are determined. In particular, the question of whether the dynamo branches bifurcate super- or sub-critically and whether a direct link to the primary pure convective states exists is answered.
Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using very deep observations (27 h) of the Hubble Deep Field South with the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M-star = 8 x 10(9) M-circle dot, star formation rate SFR = 77(-25)(+40) M-circle dot yr(-1), and star formation rate surface brightness Sigma(SFR) = 1.6 M-circle dot kpc(-2) within the [OII] lambda lambda 3727, 3729 half-light radius R-1/2, ([OII]) = 2.76 +/- 0.17 kpc. From a component of the strong resonant Mg II and Fe II absorptions at -350 km s(-1), we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant Fe II* emission, at lambda 2365, lambda 2396, lambda 2612, and lambda 2626, at 1.2-2.4-1.5-2.7 x 10-(18) erg s(-1) cm(-2) respectively. The flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant Fe II* emission lines, we spatially map the Fe II* emission from an individual galaxy for the first time. The Fe II* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius, R-1/2, (Fe II*) = 4.1 +/- 0.4 kpc, is 70% larger than the stellar continuum (R-1/2,(star) similar or equal to 2.34 +/- 0.17) or the [O II] nebular line. Moreover, the Fe II* emission shows a blue wing extending up to -400 km s(-1), which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p - v diagram along that axis. These features are consistent with a bi-conical outflow.
We present a setup combining a liquid flatjet sample delivery and a MHz laser system for time-resolved soft X-ray absorption measurements of liquid samples at the high brilliance undulator beamline UE52-SGM at Bessy II yielding unprecedented statistics in this spectral range. We demonstrate that the efficient detection of transient absorption changes in transmission mode enables the identification of photoexcited species in dilute samples. With iron(II)-trisbipyridine in aqueous solution as a benchmark system, we present absorption measurements at various edges in the soft X-ray regime. In combination with the wavelength tunability of the laser system, the set-up opens up opportunities to study the photochemistry of many systems at low concentrations, relevant to materials sciences, chemistry, and biology. (C) 2017 Author(s).
We develop a simple one-zone model of the steady-state Crab nebula spectrum encompassing both the radio/soft X-ray and the GeV/multi-TeV observations. By solving the transport equation for GeV-TeV electrons injected at the wind termination shock as a log-parabola momentum distribution and evolved via energy losses, we determine analytically the resulting differential energy spectrum of photons. We find an impressive agreement with the observed spectrum of synchrotron emission, and the synchrotron self-Compton component reproduces the previously unexplained broad 200-GeV peak that matches the Fermi/Large Area Telescope (LAT) data beyond 1 GeV with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) data. We determine the parameters of the single log-parabola electron injection distribution, in contrast with multiple broken power-law electron spectra proposed in the literature. The resulting photon differential spectrum provides a natural interpretation of the deviation from power law customarily fitted with empirical multiple broken power laws. Our model can be applied to the radio-to-multi-TeV spectrum of a variety of astrophysical outflows, including pulsar wind nebulae and supernova remnants, as well as to interplanetary shocks.
We develop a simple two-zone interpretation of the broadband baseline Crab nebula spectrum between 10(-5) eV and similar to 100 TeV by using two distinct log-parabola energetic electrons distributions. We determine analytically the very-high energy photon spectrum as originated by inverse-Compton scattering of the far-infrared soft ambient photons within the nebula off a first population of electrons energized at the nebula termination shock. The broad and flat 200 GeV peak jointly observed by Fermi/LAT and MAGIC is naturally reproduced. The synchrotron radiation from a second energetic electron population explains the spectrum from the radio range up to similar to 10 keV. We infer from observations the energy dependence of the microscopic probability of remaining in proximity of the shock of the accelerating electrons.
We study the frequently used assumption in multi-messenger astrophysics that the gamma-ray and neutrino fluxes are directly connected because they are assumed to be produced by the same photohadronic production chain. An interesting candidate source for this test is the flat-spectrum radio quasar PKS B1424-418, which recently called attention to a potential correlation between an IceCube PeV neutrino event and its burst phase. We simulate both the multi-waveband photon and the neutrino emission from this source using a self-consistent radiation model. We demonstrate that a simple hadronic model cannot adequately describe the spectral energy distribution for this source, but a lepto-hadronic model with a subdominant hadronic component can reproduce the multi-waveband photon spectrum observed during various activity phases of the blazar. As a conclusion, up to about 0.3 neutrino events may coincide with the burst, which implies that the leptonic contribution dominates in the relevant energy band. We also demonstrate that the time-wise correlation between the neutrino event and burst phase is weak.
The influence of different parameters on the silicification procedure using lysozyme is reported. When polyethoxysiloxane (PEOS), an internally crosslinked silica reservoir, is used, regular structures with a narrow size distribution could be obtained only via introducing the silica precursor in two steps including initial dropping and subsequent addition of residual oil phase in one portion. We found that mixing sequence of mineralizing agents in the presence of a positively charged surfactant plays a key role in terms of silica precipitation when tetraethoxyorthosilicate (TEOS) is the oil phase. In contrast, well mineralized crumpled features with high specific surface area could be synthesized in the presence of PEOS as a silica precursor polymer, regardless of mixing sequence. Moreover, introducing sodium dodecyl sulfate (SDS) as a negatively charged surfactant resulted in regular silica sphere formation only in combination with hexylene glycol (MPD) as a specific co-solvent. Finally, it is demonstrated that by inclusion of different nanoparticles even more sophisticated hybrid materials can be generated.
Recently a multitude of empirically derived damage models have been applied to project future tropical cyclone (TC) losses for the United States. In their study (Geiger et al 2016 Environ. Res. Lett. 11 084012) compared two approaches that differ in the scaling of losses with socio-economic drivers: the commonly-used approach resulting in a sub-linear scaling of historical TC losses with a nation's affected gross domestic product (GDP), and the disentangled approach that shows a sub-linear increase with affected population and a super-linear scaling of relative losses with per capita income. Statistics cannot determine which approach is preferable but since process understanding demands that there is a dependence of the loss on both GDP per capita and population, an approach that accounts for both separately is preferable to one which assumes a specific relation between the two dependencies. In the accompanying comment, Rybski et al argued that there is no rigorous evidence to reach the conclusion that high-income does not protect against hurricane losses. Here we affirm that our conclusion is drawn correctly and reply to further remarks raised in the comment, highlighting the adequateness of our approach but also the potential for future extension of our research.
We present rigorous results for the mean first passage time and first passage time statistics for two-channel Markov additive diffusion in a 3-dimensional spherical domain. Inspired by biophysical examples we assume that the particle can only recognise the target in one of the modes, which is shown to effect a non-trivial first passage behaviour. We also address the scenario of intermittent immobilisation. In both cases we prove that despite the perfectly non-recurrent motion of two-channel Markov additive diffusion in 3 dimensions the first passage statistics at long times do not display Poisson-like behaviour if none of the phases has a vanishing diffusion coefficient. This stands in stark contrast to the standard (one-channel) Markov diffusion counterpart. We also discuss the relevance of our results in the context of cellular signalling.
We describe analytically synchronization and desynchronization effects in an ensemble of phase oscillators driven by common noise and by global coupling. Adopting the Ott-Antonsen ansatz, we reduce the dynamics to closed stochastic equations for the order parameters, and study these equations for the cases of populations of identical and nonidentical oscillators. For nonidentical oscillators we demonstrate a counterintuitive effect of divergence of individual frequencies for moderate repulsive coupling, while the order parameter remains large.
We study a theoretical model for the toxin-antitoxin (hok/sok) mechanism for plasmid maintenance in bacteria. Toxin-antitoxin systems enforce the maintenance of a plasmid through post-segregational killing of cells that have lost the plasmid. Key to their function is the tight regulation of expression of a protein toxin by an sRNA antitoxin. Here, we focus on the nonlinear nature of the regulatory circuit dynamics of the toxin-antitoxin mechanism. The mechanism relies on a transient increase in protein concentration rather than on the steady state of the genetic circuit. Through a systematic analysis of the parameter dependence of this transient increase, we confirm some known design features of this system and identify new ones: for an efficient toxin-antitoxin mechanism, the synthesis rate of the toxin’s mRNA template should be lower that of the sRNA antitoxin, the mRNA template should be more stable than the sRNA antitoxin, and the mRNA-sRNA complex should be more stable than the sRNA antitoxin. Moreover, a short half-life of the protein toxin is also beneficial to the function of the toxin-antitoxin system. In addition, we study a therapeutic scenario in which a competitor mRNA is introduced to sequester the sRNA antitoxin, causing the toxic protein to be expressed.
Statistical properties of a free-electron laser revealed by Hanbury Brown-Twiss interferometry
(2017)
We present a comprehensive experimental analysis of statistical properties of the self-amplified spontaneous emission free-electron laser (FEL) FLASH by means of Hanbury Brown and Twiss interferometry. The experiments were performed at FEL wavelengths of 5.5, 13.4, and 20.8 nm. We determined the second-order intensity correlation function for all wavelengths and different operation conditions of FLASH. In all experiments a high degree of spatial coherence (above 50%) was obtained. Our analysis performed in spatial and spectral domains provided us with the independent measurements of an average pulse duration of the FEL that were below 60 fs. To explain the complicated behavior of the second-order intensity correlation function we developed an advanced theoretical model that includes the presence of multiple beams and external positional jitter of the FEL pulses. By this analysis we determined that in one of the experiments external positional jitter was about 25% of the beam size. We envision that methods developed in our study will be used widely for analysis and diagnostics of FEL radiation.
Networks of coupled dynamical systems provide a powerful way to model systems with enormously complex dynamics, such as the human brain. Control of synchronization in such networked systems has far-reaching applications in many domains, including engineering and medicine. In this paper, we formulate the synchronization control in dynamical systems as an optimization problem and present a multi-objective genetic programming-based approach to infer optimal control functions that drive the system from a synchronized to a non-synchronized state and vice versa. The genetic programming-based controller allows learning optimal control functions in an interpretable symbolic form. The effectiveness of the proposed approach is demonstrated in controlling synchronization in coupled oscillator systems linked in networks of increasing order complexity, ranging from a simple coupled oscillator system to a hierarchical network of coupled oscillators. The results show that the proposed method can learn highly effective and interpretable control functions for such systems.
The linear Boltzmann equation approach is generalized to describe fractional superdiffusive transport of the Levy walk type in external force fields. The time distribution between scattering events is assumed to have a finite mean value and infinite variance. It is completely characterized by the two scattering rates, one fractional and a normal one, which defines also the mean scattering rate. We formulate a general fractional linear Boltzmann equation approach and exemplify it with a particularly simple case of the Bohm and Gross scattering integral leading to a fractional generalization of the Bhatnagar, Gross and Krook kinetic equation. Here, at each scattering event the particle velocity is completely randomized and takes a value from equilibrium Maxwell distribution at a given fixed temperature. We show that the retardation effects are indispensable even in the limit of infinite mean scattering rate and argue that this novel fractional kinetic equation provides a viable alternative to the fractional Kramers-Fokker-Planck (KFP) equation by Barkai and Silbey and its generalization by Friedrich et al. based on the picture of divergent mean time between scattering events. The case of divergent mean time is also discussed at length and compared with the earlier results obtained within the fractional KFP. Also a phenomenological fractional BGK equation without retardation effects is proposed in the limit of infinite scattering rates. It cannot be, however, rigorously derived from a scattering model, being rather clever postulated. It this respect, this retardationless equation is similar to the fractional KFP by Barkai and Silbey. However, it corresponds to the opposite, much more physical limit and, therefore, also presents a viable alternative.
Can the statistical properties of single-electron transfer events be correctly predicted within a common equilibrium ensemble description? This fundamental in nanoworld question of ergodic behavior is scrutinized within a very basic semi-classical curve-crossing problem. It is shown that in the limit of non-adiabatic electron transfer (weak tunneling) well-described by the Marcus-Levich-Dogonadze (MLD) rate the answer is yes. However, in the limit of the so-called solvent-controlled adiabatic electron transfer, a profound breaking of ergodicity occurs. Namely, a common description based on the ensemble reduced density matrix with an initial equilibrium distribution of the reaction coordinate is not able to reproduce the statistics of single-trajectory events in this seemingly classical regime. For sufficiently large activation barriers, the ensemble survival probability in a state remains nearly exponential with the inverse rate given by the sum of the adiabatic curve crossing (Kramers) time and the inverse MLD rate. In contrast, near to the adiabatic regime, the single-electron survival probability is clearly non-exponential, even though it possesses an exponential tail which agrees well with the ensemble description. Initially, it is well described by a Mittag-Leffler distribution with a fractional rate. Paradoxically, the mean transfer time in this classical on the ensemble level regime is well described by the inverse of the nonadiabatic quantum tunneling rate on a single particle level. An analytical theory is developed which perfectly agrees with stochastic simulations and explains our findings.
Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains
(2017)
We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry-characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. A similar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA. We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters. We analyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.
Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains
(2017)
We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry—characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. Asimilar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA.We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters.Weanalyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.
Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains
(2017)
Westudy the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry—characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. Asimilar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA.We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters.Weanalyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.
Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase similar to 0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium, and neon whose strengths and presence depend on the overall level of absorption. These features imply a co-existence of cool and hot gas phases in the system, which we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.
Observational and computational extragalactic astrophysics are two fields of research that study a similar subject from different perspectives. Observational extragalactic astrophysics aims, by recovering the spectral energy distribution of galaxies at different wavelengths, to reliably measure their properties at different cosmic times and in a large variety of environments. Analyzing the light collected by the instruments, observers try to disentangle the different processes occurring in galaxies at the scales of galactic physics, as well as the effect of larger scale processes such as mergers and accretion, in order to obtain a consistent picture of galaxy formation and evolution. On the other hand, hydrodynamical simulations of galaxy formation in cosmological context are able to follow the evolution of a galaxy along cosmic time, taking into account both external processes such as mergers, interactions and accretion, and internal mechanisms such as feedback from Supernovae and Active Galactic Nuclei. Due to the great advances in both fields of research, we have nowadays available spectral and photometric information for a large number of galaxies in the Universe at different cosmic times, which has in turn provided important knowledge about the evolution of the Universe; at the same time, we are able to realistically simulate galaxy formation and evolution in large volumes of the Universe, taking into account the most relevant physical processes occurring in galaxies.
As these two approaches are intrinsically different in their methodology and in the information they provide, the connection between simulations and observations is still not fully established, although simulations are often used in galaxies' studies to interpret observations and assess the effect of the different processes acting on galaxies on the observable properties, and simulators usually test the physical recipes implemented in their hydrodynamical codes through the comparison with observations. In this dissertation we aim to better connect the observational and computational approaches in the study of galaxy formation and evolution, using the methods and results of one field to test and validate the methods and results of the other.
In a first work we study the biases and systematics in the derivation of the galaxy properties in observations. We post-process hydrodynamical cosmological simulations of galaxy formation to calculate the galaxies' Spectral Energy Distributions (SEDs) using different approaches, including radiative transfer techniques. Comparing the direct results of the simulations with the quantities obtained applying observational techniques to these synthetic SEDs, we are able to make an analysis of the biases intrinsic in the observational algorithms, and quantify their accuracy in recovering the galaxies' properties, as well as estimating the uncertainties affecting a comparison between simulations and observations when different approaches to obtain the observables are followed. Our results show that for some quantities such as the stellar ages, metallicities and gas oxygen abundances large differences can appear, depending on the technique applied in the derivation.
In a second work we compare a set of fifteen galaxies similar in mass to the Milky Way and with a quiet merger history in the recent past (hence expected to have properties close to spiral galaxies), simulated in a cosmological context, with data from the Sloan Digital Sky Survey (SDSS). We use techniques to obtain the observables as similar as possible to the ones applied in SDSS, with the aim of making an unbiased comparison between our set of hydrodynamical simulations and SDSS observations. We quantify the differences in the physical properties when these are obtained directly from the simulations without post-processing, or mimicking the SDSS observational techniques. We fit linear relations between the values derived directly from the simulations and following SDSS observational procedures, which in most of the cases have relatively high correlation, that can be easily used to more reliably compare simulations with SDSS data. When mimicking SDSS techniques, these simulated galaxies are photometrically similar to galaxies in the SDSS blue sequence/green valley, but have in general older ages, lower SFRs and metallicities compared to the majority of the spirals in the observational dataset.
In a third work, we post-process hydrodynamical simulations of galaxies with radiative transfer techniques, to generate synthetic data that mimic the properties of the CALIFA Integral Field Spectroscopy (IFS) survey. We reproduce the main characteristics of the CALIFA observations in terms of field of view and spaxel physical size, data format, point spread functions and detector noise. This 3-dimensional dataset is suited to be analyzed by the same algorithms applied to the CALIFA dataset, and can be used as a tool to test the ability of the observational algorithms in recovering the properties of the CALIFA galaxies. To this purpose, we also generate the resolved maps of the simulations' properties, calculated directly from the hydrodynamical snapshots, or from the simulated spectra prior to the addition of the noise.
Our work shows that a reliable connection between the models and the data is of crucial importance both to judge the output of galaxy formation codes and to accurately test the observational algorithms used in the analysis of galaxy surveys' data. A correct interpretation of observations will be particularly important in the future, in light of the several ongoing and planned large galaxy surveys that will provide the community with large datasets of properties of galaxies (often spatially-resolved) at different cosmic times, allowing to study galaxy formation physics at a higher level of detail than ever before. We have shown that neglecting the observational biases in the comparison between simulations and an observational dataset may move the simulations to different regions in the planes of the observables, strongly affecting the assessment of the correctness of the sub-resolution physical models implemented in galaxy formation codes, as well as the interpretation of given observational results using simulations.
The influence of relaxation processes on the thermal electret stability of high-impact polystyrene (HIPS) free-standing films filled with titanium dioxide (TiO2) of the rutile modification are investigated by means of a combination of dielectric methods (dielectric relaxation spectroscopy (DRS), thermally stimulated depolarization current (TSDC) and thermally stimulated surface-potential decay (TSSPD)), supplemented by differential scanning calorimetry (DSC) and dynamic mechanical analysis (DMA). Films with 2, 4, 6 and 8 vol.% TiO2 are compared to each other and to non-filled samples. Filling HIPS with up to 8 vol.% of TiO2 enhances the elastic modulus below the glass transition and increases the thermal electret stability above the glass transition without significantly increasing the DC conductivity. The improvement of the electret stability is caused by the build-up of an interface polarization which decays only gradually if the glass transition is exceeded. Two kinds of Arrhenius processes are considered in order to explain the decay of the composite-polymer electrets: (1) charge release from chemical traps located at the phenyl rings of the polymer chain with an activation energy of E-a = 1.1 eV after passing the glass transition at about 100 degrees C and (2) charge release from traps formed by the TiO2 particles with E-a = 2.4 eV at temperatures above 130 degrees C. Finally, the activation energies are discussed with respect to their significance.
Aims. Changes of the magnetic field and the line-of-sight velocities in the photosphere are being reported for an M-class flare that originated at a delta-spot belonging to active region NOAA 11865. Methods. High-resolution ground-based near-infrared spectropolarimetric observations were acquired simultaneously in two photospheric spectral lines, Fe I 10783 angstrom and Si I 10786 angstrom, with the Tenerife Infrared Polarimeter at the Vacuum Tower Telescope (VTT) in Tenerife on 2013 October 15. The observations covered several stages of the M-class flare. Inversions of the full-Stokes vector of both lines were carried out and the results were put into context using (extreme)-ultraviolet filtergrams from the Solar Dynamics Observatory (SDO). Results. The active region showed high flaring activity during the whole observing period. After the M-class flare, the longitudinal magnetic field did not show significant changes along the polarity inversion line (PIL). However, an enhancement of the transverse magnetic field of approximately 550G was found that bridges the PIL and connects umbrae of opposite polarities in the delta-spot. At the same time, a newly formed system of loops appeared co-spatially in the corona as seen in 171 angstrom filtergrams of the Atmospheric Imaging Assembly (AIA) on board SDO. However, we cannot exclude that the magnetic connection between the umbrae already existed in the upper atmosphere before the M-class flare and became visible only later when it was filled with hot plasma. The photospheric Doppler velocities show a persistent upflow pattern along the PIL without significant changes due to the flare. Conclusions. The increase of the transverse component of the magnetic field after the flare together with the newly formed loop system in the corona support recent predictions of flare models and flare observations.
Galaxies are surrounded by sizeable gas reservoirs which host a significant amount of metals: the circum-galactic medium (CGM). The CGM acts as a mediator between the galaxy and the extragalactic medium. However, our understanding of how galaxy mergers, a major evolutionary transformation, impact the CGM remains deficient. We present a theoretical study of the effect of galaxy mergers on the CGM. We use hydrodynamical cosmological zoom-in simulations of a major merger selected from the Illustris project such that the z = 0 descendant has a halo mass and stellar mass comparable to the Milky Way. To study the CGM we then re-simulated this system at a 40 times better mass resolution, and included detailed post-processing ionization modelling. Our work demonstrates the effect the merger has on the characteristic size of the CGM, its metallicity, and the predicted covering fraction of various commonly observed gas-phase species, such as H I, C IV, and O VI. We show that merger-induced outflows can increase the CGM metallicity by 0.2-0.3 dex within 0.5 Gyr post-merger. These effects last up to 6 Gyr post-merger. While the merger increases the total metal covering fractions by factors of 2-3, the covering fractions of commonly observed UV ions decrease due to the hard ionizing radiation from the active galactic nucleus, which we model explicitly. Our study of the single simulated major merger presented in this work demonstrates the significant impact that a galaxy interaction can have on the size, metallicity, and observed column densities of the CGM.
Background: Inferring regulatory interactions between genes from transcriptomics time-resolved data, yielding reverse engineered gene regulatory networks, is of paramount importance to systems biology and bioinformatics studies. Accurate methods to address this problem can ultimately provide a deeper insight into the complexity, behavior, and functions of the underlying biological systems. However, the large number of interacting genes coupled with short and often noisy time-resolved read-outs of the system renders the reverse engineering a challenging task. Therefore, the development and assessment of methods which are computationally efficient, robust against noise, applicable to short time series data, and preferably capable of reconstructing the directionality of the regulatory interactions remains a pressing research problem with valuable applications.
Results: Here we perform the largest systematic analysis of a set of similarity measures and scoring schemes within the scope of the relevance network approach which are commonly used for gene regulatory network reconstruction from time series data. In addition, we define and analyze several novel measures and schemes which are particularly suitable for short transcriptomics time series. We also compare the considered 21 measures and 6 scoring schemes according to their ability to correctly reconstruct such networks from short time series data by calculating summary statistics based on the corresponding specificity and sensitivity. Our results demonstrate that rank and symbol based measures have the highest performance in inferring regulatory interactions. In addition, the proposed scoring scheme by asymmetric weighting has shown to be valuable in reducing the number of false positive interactions. On the other hand, Granger causality as well as information-theoretic measures, frequently used in inference of regulatory networks, show low performance on the short time series analyzed in this study.
Conclusions: Our study is intended to serve as a guide for choosing a particular combination of similarity measures and scoring schemes suitable for reconstruction of gene regulatory networks from short time series data. We show that further improvement of algorithms for reverse engineering can be obtained if one considers measures that are rooted in the study of symbolic dynamics or ranks, in contrast to the application of common similarity measures which do not consider the temporal character of the employed data. Moreover, we establish that the asymmetric weighting scoring scheme together with symbol based measures (for low noise level) and rank based measures (for high noise level) are the most suitable choices.
Nanoscale Thermal Transfer
(2017)
We analyze the cross-over of a homogeneous, weakly interacting Bose gas in one dimension from the ideal gas into the dense quasi-condensate phase. We review a number of mean-field theories, perturbative or self-consistent, and provide accurate evaluations of equation of state, density fluctuations, and correlation functions. A smooth crossover is reproduced by classical-field simulations based on the stochastic Gross-Pitaevskii equation and the Yang-Yang solution to the one-dimensional Bose gas.
Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.
Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.
Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.