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)
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.