Institut für Physik und Astronomie
Refine
Year of publication
Document Type
- Postprint (136) (remove)
Language
- English (136) (remove)
Keywords
- diffusion (16)
- anomalous diffusion (9)
- models (3)
- nonergodicity (3)
- transport (3)
- azobenzene (2)
- brushes (2)
- dynamics (2)
- films (2)
- first passage time (2)
Institute
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
We use ultrafast x-ray diffraction to investigate the effect of expansive phononic and contractive magnetic stress driving the picosecond strain response of a metallic perovskite SrRuO3 thin film upon femtosecond laser excitation. We exemplify how the anisotropic bulk equilibrium thermal expansion can be used to predict the response of the thin film to ultrafast deposition of energy. It is key to consider that the laterally homogeneous laser excitation changes the strain response compared to the near-equilibrium thermal expansion because the balanced in-plane stresses suppress the Poisson stress on the picosecond timescale. We find a very large negative Grüneisen constant describing the large contractive stress imposed by a small amount of energy in the spin system. The temperature and fluence dependence of the strain response for a double-pulse excitation scheme demonstrates the saturation of the magnetic stress in the high-fluence regime.
Using time-resolved x-ray diffraction, we demonstrate the manipulation of the picosecond strain response of a metallic heterostructure consisting of a dysprosium (Dy) transducer and a niobium (Nb) detection layer by an external magnetic field. We utilize the first-order ferromagnetic–antiferromagnetic phase transition of the Dy layer, which provides an additional large contractive stress upon laser excitation compared to its zerofield response. This enhances the laser-induced contraction of the transducer and changes the shape of the picosecond strain pulses driven in Dy and detected within the buried Nb layer. Based on our experiment with rare-earth metals we discuss required properties for functional transducers, which may allow for novel field-control of the emitted picosecond strain pulses.
Pre-service physics teachers often do not recognise the relevance for their future career in their university content knowledge courses. A lower perceived relevance can, however, have a negative effect on their motivation and on their academic success. Several intervention studies have been undertaken with the goal to increase this perceived relevance. A previous study shows that conceptual physics problems used in university physics courses are perceived by pre-service physics teachers as more relevant for their future career than regular, quantitative problems. It is however not clear, what the students' meaning of the construct 'relevance' is: what makes a problem more relevant to them than another problem? To answer this question, N = 7 pre-service teachers were interviewed using the repertory grid technique, based on the personal construct theory. Nine physics problems were discussed with regards to their perceived relevance and with regards to problem properties that distinguish these problems from each other. We are able to identify six problem properties that have a positive influence on the perceived relevance. Physics problems that are based on these properties should therefore potentially have a higher perceived relevance, which can have a positive effect on the motivation of the pre-service teachers who solve these problems.
We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T−h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.
The Ornstein–Uhlenbeck process is a stationary and ergodic Gaussian process, that is fully determined by its covariance function and mean. We show here that the generic definitions of the ensemble- and time-averaged mean squared displacements fail to capture these properties consistently, leading to a spurious ergodicity breaking. We propose to remedy this failure by redefining the mean squared displacements such that they reflect unambiguously the statistical properties of any stochastic process. In particular we study the effect of the initial condition in the Ornstein–Uhlenbeck process and its fractional extension. For the fractional Ornstein–Uhlenbeck process representing typical experimental situations in crowded environments such as living biological cells, we show that the stationarity of the process delicately depends on the initial condition.
Persistent episodes of extreme weather in the Northern Hemisphere summer have been associated with high-amplitude quasi-stationary atmospheric Rossby waves, with zonal wave numbers 6 to 8 resulting from the phenomenon of quasi-resonant amplification (QRA). A fingerprint for the occurrence of QRA can be defined in terms of the zonally averaged surface temperature field. Examining state-of-the-art [Coupled Model Intercomparison Project Phase 5 (CMIP5)] climate model projections, we find that QRA events are likely to increase by similar to 50% this century under business-as-usual carbon emissions, but there is considerable variation among climate models. Some predict a near tripling of QRA events by the end of the century, while others predict a potential decrease. Models with amplified Arctic warming yield the most pronounced increase in QRA events. The projections are strongly dependent on assumptions regarding the nature of changes in radiative forcing associated with anthropogenic aerosols over the next century. One implication of our findings is that a reduction in midlatitude aerosol loading could actually lead to Arctic de-amplification this century, ameliorating potential increases in persistent extreme weather events.
We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes-namely, the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving while varying the Alfven Mach number (M-A), plasma beta, and the sonic Mach number from subsonic to transsonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven-mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high M-A. The slow-mode properties are similar to the Alfven mode. On the other hand, the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast-mode behavior does not show any transition in going from low to high M-A. We find indications that there is some interaction between the different modes, and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic-ray transport and magnetic reconnection are discussed.
Despite the ongoing progress in nanotechnology and its applications, the development of strategies for connecting nano-scale systems to micro- or macroscale elements is hampered by the lack of structural components that have both, nano- and macroscale dimensions. The production of nano-scale wires with macroscale length is one of the most interesting challenges here. There are a lot of strategies to fabricate long nanoscopic stripes made of metals, polymers or ceramics but none is suitable for mass production of ordered and dense arrangements of wires at large numbers. In this paper, we report on a technique for producing arrays of ordered, flexible and free-standing polymer nano-wires filled with different types of nano-particles. The process utilizes the strong response of photosensitive polymer brushes to irradiation with UV-interference patterns, resulting in a substantial mass redistribution of the polymer material along with local rupturing of polymer chains. The chains can wind up in wires of nano-scale thickness and a length of up to several centimeters. When dispersing nano-particles within the film, the final arrangement is similar to a core-shell geometry with mainly nano-particles found in the core region and the polymer forming a dielectric jacket.
We demonstrate in situ recorded motion of nano-objects adsorbed on a photosensitive polymer film. The motion is induced by a mass transport of the underlying photoresponsive polymer material occurring during irradiation with interference pattern. The polymer film contains azobenzene molecules that undergo reversible photoisomerization reaction from trans- to cis-conformation. Through a multi-scale chain of physico-chemical processes, this finally results in the macro-deformations of the film due to the changing elastic properties of polymer. The topographical deformation of the polymer surface is sensitive to a local distribution of the electrical field vector that allows for the generation of dynamic changes in the surface topography during irradiation with different light interference patterns. Polymer film deformation together with the motion of the adsorbed nano-particles are recorded using a homemade set-up combining an optical part for the generation of interference patterns and an atomic force microscope for acquiring the surface deformation. The particles undergo either translational or rotational motion. The direction of particle motion is towards the topography minima and opposite to the mass transport within the polymer film. The ability to relocate particles by photo-induced dynamic topography fluctuation offers a way for a non-contact simultaneous manipulation of a large number of adsorbed particles just in air at ambient conditions.
Herein, we report the chain-growth tin-free room temperature polymerization method to synthesize n-type perylene diimide-dithiophene-based conjugated polymers (PPDIT2s) suitable for solar cell and transistor applications. The palladium/electron-rich tri-tert-butylphosphine catalyst is effective to enable the chain-growth polymerization of anion-radical monomer Br-TPDIT-Br/Zn to PPDIT2 with a molecular weight up to Mw ≈ 50 kg mol−1 and moderate polydispersity. This is the second example of the polymerization of unusual anion-radical aromatic complexes formed in a reaction of active Zn and electron-deficient diimide-based aryl halides. As such, the discovered polymerization method is not a specific reactivity feature of the naphthalene-diimide derivatives but is rather a general polymerization tool. This is an important finding, given the significantly higher maximum external quantum efficiency that can be reached with PDI-based copolymers (32–45%) in all-polymer solar cells compared to NDI-based materials (15–30%). Our studies revealed that PPDIT2 synthesized by the new method and the previously published polymer prepared by step-growth Stille polycondensation show similar electron mobility and all-polymer solar cell performance. At the same time, the polymerization reported herein has several technological advantages as it proceeds relatively fast at room temperature and does not involve toxic tin-based compounds. Because several chain-growth polymerization reactions are well-suited for the preparation of well-defined multi-functional polymer architectures, the next target is to explore the utility of the discovered polymerization in the synthesis of end-functionalized polymers and block copolymers. Such materials would be helpful to improve the nanoscale morphology of polymer blends in all-polymer solar cells.
This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the Itô calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.
In this paper, we report X-ray absorption and core-level electron spectra of the nucleobase derivative 2-thiouracil at the sulfur L1- and L2,3-edges. We used soft X-rays from the free-electron laser FLASH2 for the excitation of isolated molecules and dispersed the outgoing electrons with a magnetic bottle spectrometer. We identified photoelectrons from the 2p core orbital, accompanied by an electron correlation satellite, as well as resonant and non-resonant Coster–Kronig and Auger–Meitner emission at the L1- and L2,3-edges, respectively. We used the electron yield to construct X-ray absorption spectra at the two edges. The experimental data obtained are put in the context of the literature currently available on sulfur core-level and 2-thiouracil spectroscopy.
Recombination of free charge is a key process limiting the performance of solar cells. For low mobility materials, such as organic semiconductors, the kinetics of non-geminate recombination (NGR) is strongly linked to the motion of charges. As these materials possess significant disorder, thermalization of photogenerated carriers in the inhomogeneously broadened density of state distribution is an unavoidable process. Despite its general importance, knowledge about the kinetics of NGR in complete organic solar cells is rather limited. We employ time delayed collection field (TDCF) experiments to study the recombination of photogenerated charge in the high-performance polymer:fullerene blend PCDTBT:PCBM. NGR in the bulk of this amorphous blend is shown to be highly dispersive, with a continuous reduction of the recombination coefficient throughout the entire time scale, until all charge carriers have either been extracted or recombined. Rapid, contact-mediated recombination is identified as an additional loss channel, which, if not properly taken into account, would erroneously suggest a pronounced field dependence of charge generation. These findings are in stark contrast to the results of TDCF experiments on photovoltaic devices made from ordered blends, such as P3HT:PCBM, where non-dispersive recombination was proven to dominate the charge carrier dynamics under application relevant conditions.
Open-circuit voltages of lead-halide perovskite solar cells are improving rapidly and are approaching the thermodynamic limit. Since many different perovskite compositions with different bandgap energies are actively being investigated, it is not straightforward to compare the open-circuit voltages between these devices as long as a consistent method of referencing is missing. For the purpose of comparing open-circuit voltages and identifying outstanding values, it is imperative to use a unique, generally accepted way of calculating the thermodynamic limit, which is currently not the case. Here a meta-analysis of methods to determine the bandgap and a radiative limit for open-circuit voltage is presented. The differences between the methods are analyzed and an easily applicable approach based on the solar cell quantum efficiency as a general reference is proposed.
We study the experimentally measured ciprofloxacin antibiotic diffusion through a gel-like artificial sputum medium (ASM) mimicking physiological conditions typical for a cystic fibrosis layer, in which regions occupied by Pseudomonas aeruginosa bacteria are present. To quantify the antibiotic diffusion dynamics we employ a phenomenological model using a subdiffusion-absorption equation with a fractional time derivative. This effective equation describes molecular diffusion in a medium structured akin Thompson’s plumpudding model; here the ‘pudding’ background represents the ASM and the ‘plums’ represent the bacterial biofilm. The pudding is a subdiffusion barrier for antibiotic molecules that can affect bacteria found in plums. For the experimental study we use an interferometric method to determine the time evolution of the amount of antibiotic that has diffused through the biofilm. The theoretical model shows that this function is qualitatively different depending on whether or not absorption of the antibiotic in the biofilm occurs. We show that the process can be divided into three successive stages: (1) only antibiotic subdiffusion with constant biofilm parameters, (2) subdiffusion and absorption of antibiotic molecules with variable biofilm transport parameters, (3) subdiffusion and absorption in the medium but the biofilm parameters are constant again. Stage 2 is interpreted as the appearance of an intensive defence build–up of bacteria against the action of the antibiotic, and in the stage 3 it is likely that the bacteria have been inactivated. Times at which stages change are determined from the experimentally obtained temporal evolution of the amount of antibiotic that has diffused through the ASM with bacteria. Our analysis shows good agreement between experimental and theoretical results and is consistent with the biologically expected biofilm response. We show that an experimental method to study the temporal evolution of the amount of a substance that has diffused through a biofilm is useful in studying the processes occurring in a biofilm. We also show that the complicated biological process of antibiotic diffusion in a biofilm can be described by a fractional subdiffusion-absorption equation with subdiffusion and absorption parameters that change over time.
Many cellular processes require decision making mechanisms, which must act reliably even in the unavoidable presence of substantial amounts of noise. However, the multistable genetic switches that underlie most decision-making processes are dominated by fluctuations that can induce random jumps between alternative cellular states. Here we show, via theoretical modeling of a population of noise-driven bistable genetic switches, that reliable timing of decision-making processes can be accomplished for large enough population sizes, as long as cells are globally coupled by chemical means. In the light of these results, we conjecture that cell proliferation, in the presence of cell-cell communication, could provide a mechanism for reliable decision making in the presence of noise, by triggering cellular transitions only when the whole cell population reaches a certain size. In other words , the summation performed by the cell population would average out the noise and reduce its detrimental impact.
Photoluminescence spectroscopy is a widely applied characterization technique for semiconductor materials in general and halide perovskite solar cell materials in particular. It can give direct information on the recombination kinetics and processes as well as the internal electrochemical potential of free charge carriers in single semiconductor layers, layer stacks with transport layers, and complete solar cells. The correct evaluation and interpretation of photoluminescence requires the consideration of proper excitation conditions, calibration and application of the appropriate approximations to the rather complex theory, which includes radiative recombination, non-radiative recombination, interface recombination, charge transfer, and photon recycling. In this article, an overview is given of the theory and application to specific halide perovskite compositions, illustrating the variables that should be considered when applying photoluminescence analysis in these materials.
Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.
Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used.
In this paper we report on photoswitchable polymer surfaces with dynamically and reversibly fluctuating topographies. It is well known that when azobenzene containing polymer films are irradiated with optical interference patterns the film topography changes to form a surface relief grating. In the simplest case, the film shape mimics the intensity distribution and deforms into a wave like, sinusoidal manner with amplitude that may be as large as the film thickness. This process takes place in the glassy state without photo-induced softening. Here we report on an intriguing discovery regarding the formation of reliefs under special illumination conditions. We have developed a novel setup combining the optical part for creating interference patterns, an AFM for in situ acquisition of topography changes and diffraction efficiency signal measurements. In this way we demonstrate that these gratings can be “set in motion” like water waves or dunes in the desert. We achieve this by applying repetitive polarization changes to the incoming interference pattern. Such light responsive surfaces represent the prerequisite for providing practical applications ranging from conveyer or transport systems for adsorbed liquid objects and colloidal particles to generation of adaptive and dynamic optical devices.
Estimating global mean sea-level rise and its uncertainties by 2100 and 2300 from an expert survey
(2020)
Sea-level rise projections and knowledge of their uncertainties are vital to make informed mitigation and adaptation decisions. To elicit projections from members of the scientific community regarding future global mean sea-level (GMSL) rise, we repeated a survey originally conducted five years ago. Under Representative Concentration Pathway (RCP) 2.6, 106 experts projected a likely (central 66% probability) GMSL rise of 0.30-0.65 m by 2100, and 0.54-2.15 m by 2300, relative to 1986-2005. Under RCP 8.5, the same experts projected a likely GMSL rise of 0.63-1.32 m by 2100, and 1.67-5.61 m by 2300. Expert projections for 2100 are similar to those from the original survey, although the projection for 2300 has extended tails and is higher than the original survey. Experts give a likelihood of 42% (original survey) and 45% (current survey) that under the high-emissions scenario GMSL rise will exceed the upper bound (0.98 m) of the likely range estimated by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, which is considered to have an exceedance likelihood of 17%. Responses to open-ended questions suggest that the increases in upper-end estimates and uncertainties arose from recent influential studies about the impact of marine ice cliff instability on the meltwater contribution to GMSL rise from the Antarctic Ice Sheet.
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.
A novel atomic beam splitter, using reflection of atoms off an evanescent light wave, is investigated theoretically. The intensity or frequency of the light is modulated in order to create sidebands on the reflected de Broglie wave. The weights and phases of the various sidevands are calculated using three different approaches: the Born approximation, a semiclassical path integral approach, and a numerical solution of the time-dependent Schrdinger equation. We show how this modulated mirror could be used to build practical atomic interferometers.
A detailed theoretical investigation of the reflection of an atomic de Broglie wave at an evanescent wave mirror is presented. The classical and the semiclassical descriptions of the reflection process are reviewed, and a full wave-mechanical approach based on the analytical soution of the corresponding Schrödinger equation is presented. The phase shift at reflection is calculated exactly and interpreted in terms of instantaneous reflection of the atom at an effective mirror. Besides the semiclassical regime of reflection describable by the WKB method, a pure quantum regime of reflection is identified in the limit where the incident de Broglie wavelength is large compared to the evanescent wave decay length.
We present a semiclassical perturbation method for the description of atomic diffraction by a weakly modulated potential. It proceeds in a way similar to the treatment of light diffraction by a thin phase grating, and consists in calculating the atomic wavefunction by means of action integrals along the classical trajectories of the atoms in the absence of the modulated part of the potential. The capabilities and the validity condition of the method are illustrated on the well-known case of atomic diffraction by a Gaussian standing wave. We prove that in this situation the perturbation method is equivalent to the Raman-Nath approximation, and we point out that the usually-considered Raman-Nath validity condition can lead to inaccuracies in the evaluation of the phases of the diffraction amplitudes. The method is also applied to the case of an evanescent wave reflection grating, and an analytical expression for the diffraction pattern at any incidence angle is obtained for the first time. Finally, the application of the method to other situations is briefly discussed.
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.
The molecular ability to selectively and efficiently convert sunlight into other forms of energy like heat, bond change, or charge separation is truly remarkable. The decisive steps in these transformations often happen on a femtosecond timescale and require transitions among different electronic states that violate the Born-Oppenheimer approximation (BOA). Non-BOA transitions pose challenges to both theory and experiment. From a theoretical point of view, excited state dynamics and nonadiabatic transitions both are difficult problems (see Figure 1(a)). However, the theory on non-BOA dynamics has advanced significantly over the last two decades. Full dynamical simulations for molecules of the size of nucleobases have been possible for a couple of years and allow predictions of experimental observables like photoelectron energy or ion yield. The availability of these calculations for isolated molecules has spurred new experimental efforts to develop methods that are sufficiently different from all optical techniques. For determination of transient molecular structure, femtosecond X-ray diffraction and electron diffraction have been implemented on optically excited molecules.
Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.
We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.
We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.
We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.
In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.
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.
We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary.
We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed.
Stochastic Wilson
(2015)
We consider a simple Markovian class of the stochastic Wilson–Cowan type models of neuronal network dynamics, which incorporates stochastic delay caused by the existence of a refractory period of neurons. From the point of view of the dynamics of the individual elements, we are dealing with a network of non-Markovian stochastic two-state oscillators with memory, which are coupled globally in a mean-field fashion. This interrelation of a higher-dimensional Markovian and lower-dimensional non-Markovian dynamics is discussed in its relevance to the general problem of the network dynamics of complex elements possessing memory. The simplest model of this class is provided by a three-state Markovian neuron with one refractory state, which causes firing delay with an exponentially decaying memory within the two-state reduced model. This basic model is used to study critical avalanche dynamics (the noise sustained criticality) in a balanced feedforward network consisting of the excitatory and inhibitory neurons. Such avalanches emerge due to the network size dependent noise (mesoscopic noise). Numerical simulations reveal an intermediate power law in the distribution of avalanche sizes with the critical exponent around −1.16. We show that this power law is robust upon a variation of the refractory time over several orders of magnitude. However, the avalanche time distribution is biexponential. It does not reflect any genuine power law dependence.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.
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.
What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive–repulsive rod–rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed.
Fluorination of organic spacer impacts on the structural and optical response of 2D perovskites
(2020)
Low-dimensional hybrid perovskites have triggered significant research interest due to their intrinsically tunable optoelectronic properties and technologically relevant material stability. In particular, the role of the organic spacer on the inherent structural and optical features in two-dimensional (2D) perovskites is paramount for material optimization. To obtain a deeper understanding of the relationship between spacers and the corresponding 2D perovskite film properties, we explore the influence of the partial substitution of hydrogen atoms by fluorine in an alkylammonium organic cation, resulting in (Lc)(2)PbI4 and (Lf)(2)PbI4 2D perovskites, respectively. Consequently, optical analysis reveals a clear 0.2 eV blue-shift in the excitonic position at room temperature. This result can be mainly attributed to a band gap opening, with negligible effects on the exciton binding energy. According to Density Functional Theory (DFT) calculations, the band gap increases due to a larger distortion of the structure that decreases the atomic overlap of the wavefunctions and correspondingly bandwidth of the valence and conduction bands. In addition, fluorination impacts the structural rigidity of the 2D perovskite, resulting in a stable structure at room temperature and the absence of phase transitions at a low temperature, in contrast to the widely reported polymorphism in some non-fluorinated materials that exhibit such a phase transition. This indicates that a small perturbation in the material structure can strongly influence the overall structural stability and related phase transition of 2D perovskites, making them more robust to any phase change. This work provides key information on how the fluorine content in organic spacer influence the structural distortion of 2D perovskites and their optical properties which possess remarkable importance for future optoelectronic applications, for instance in the field of light-emitting devices or sensors.
The strong adhesion of sub-micron sized particles to surfaces is a nuisance, both for removing contaminating colloids from surfaces and for conscious manipulation of particles to create and test novel micro/nano-scale assemblies. The obvious idea of using detergents to ease these processes suffers from a lack of control: the action of any conventional surface-modifying agent is immediate and global. With photosensitive azobenzene containing surfactants we overcome these limitations. Such photo-soaps contain optical switches (azobenzene molecules), which upon illumination with light of appropriate wavelength undergo reversible trans-cis photo-isomerization resulting in a subsequent change of the physico-chemical molecular properties. In this work we show that when a spatial gradient in the composition of trans- and cis- isomers is created near a solid-liquid interface, a substantial hydrodynamic flow can be initiated, the spatial extent of which can be set, e.g., by the shape of a laser spot. We propose the concept of light induced diffusioosmosis driving the flow, which can remove, gather or pattern a particle assembly at a solid-liquid interface. In other words, in addition to providing a soap we implement selectivity: particles are mobilized and moved at the time of illumination, and only across the illuminated area.
Ground-based gamma-ray astronomy has had a major breakthrough with the impressive results obtained using systems of imaging atmospheric Cherenkov telescopes. Ground-based gamma-ray astronomy has a huge potential in astrophysics, particle physics and cosmology. CTA is an international initiative to build the next generation instrument, with a factor of 5-10 improvement in sensitivity in the 100 GeV-10 TeV range and the extension to energies well below 100 GeV and above 100 TeV. CTA will consist of two arrays (one in the north, one in the south) for full sky coverage and will be operated as open observatory. The design of CTA is based on currently available technology. This document reports on the status and presents the major design concepts of CTA.
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.
We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems.
Based on extensive Monte Carlo simulations and analytical considerations we study the electrostatically driven adsorption of flexible polyelectrolyte chains onto charged Janus nanospheres. These net-neutral colloids are composed of two equally but oppositely charged hemispheres. The critical binding conditions for polyelectrolyte chains are analysed as function of the radius of the Janus particle and its surface charge density, as well as the salt concentration in the ambient solution. Specifically for the adsorption of finite-length polyelectrolyte chains onto Janus nanoparticles, we demonstrate that the critical adsorption conditions drastically differ when the size of the Janus particle or the screening length of the electrolyte are varied. We compare the scaling laws obtained for the adsorption–desorption threshold to the known results for uniformly charged spherical particles, observing significant disparities. We also contrast the changes to the polyelectrolyte chain conformations close to the surface of the Janus nanoparticles as compared to those for simple spherical particles. Finally, we discuss experimentally relevant physico-chemical systems for which our simulations results may become important. In particular, we observe similar trends with polyelectrolyte complexation with oppositely but heterogeneously charged proteins.
We have investigated the electrochemical, spectroscopic and electroluminescent properties of a family of aza-aromatic complexes of ruthenium of type [RuII(bpy/phen)2(L)]2+ (4d6) with three isomeric L ligands, where, bpy = 2,2′-bipyridine, phen = 1,10-phenanthroline and the L ligands are 3-(2-pyridyl)[1,2,4]triazolo[1,5-a]pyridine (L1), 3-(2-pyridyl[1,2,3])triazolo[1,5-a]pyridine (L2) and 2-(2-pyridyl)[1,2,4]triazolo[1,5-a]pyridine (L3). The complexes display two bands in the visible region near 410–420 and 440–450 nm. The complexes are diamagnetic and show well defined 1H NMR lines. They are electroactive in acetonitrile solution and exhibit a well defined RuII/RuIII couple near 1.20 to 1.30 V and −1.40 to −1.50 V due to ligand reduction versus Saturated Calomel Electrode (SCE). The solutions are also luminescent, with peaks are near 600 nm. All the complexes are electroluminescent in nature with peaks lying near 580 nm. L1 and L3 ligated complexes with two bpy co-ligands show weak photoluminescence (PL) but stronger electroluminescence (EL) compared to corresponding L2 ligated analogues.
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.
Droughts in tropical South America have an imminent and severe impact on the Amazon rainforest and affect the livelihoods of millions of people. Extremely dry conditions in Amazonia have been previously linked to sea surface temperature (SST) anomalies in the adjacent tropical oceans. Although the sources and impacts of such droughts have been widely studied, establishing reliable multi-year lead statistical forecasts of their occurrence is still an ongoing challenge. Here, we further investigate the relationship between SST and rainfall anomalies using a complex network approach. We identify four ocean regions which exhibit the strongest overall SST correlations with central Amazon rainfall, including two particularly prominent regions in the northern and southern tropical Atlantic. Based on the time-dependent correlation between SST anomalies in these two regions alone, we establish a new early-warning method for droughts in the central Amazon basin and demonstrate its robustness in hindcasting past major drought events with lead-times up to 18 months.
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.