530 Physik
Refine
Year of publication
Document Type
- Article (883)
- Doctoral Thesis (371)
- Postprint (123)
- Preprint (50)
- Other (48)
- Habilitation Thesis (23)
- Master's Thesis (10)
- Review (10)
- Monograph/Edited Volume (4)
- Course Material (3)
Keywords
- diffusion (43)
- anomalous diffusion (33)
- gamma rays: general (20)
- synchronization (19)
- Synchronisation (16)
- organic solar cells (15)
- stochastic processes (15)
- cosmic rays (14)
- ISM: supernova remnants (13)
- data analysis (12)
Institute
- Institut für Physik und Astronomie (1397)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (52)
- Extern (41)
- Institut für Chemie (38)
- Mathematisch-Naturwissenschaftliche Fakultät (25)
- Institut für Mathematik (21)
- Institut für Geowissenschaften (14)
- Institut für Biochemie und Biologie (7)
- Institut für Umweltwissenschaften und Geographie (7)
- Potsdam Institute for Climate Impact Research (PIK) e. V. (6)
Increasing Earth’s surface air temperature yields an intensification of its hydrological cycle. As a consequence, the risk of river floods will increase regionally within the next two decades due to the atmospheric warming caused by past anthropogenic greenhouse gas emissions. The direct economic losses caused by these floods can yield regionally heterogeneous losses and gains by propagation within the global trade and supply network. Here we show that, in the absence of large-scale structural adaptation, the total economic losses due to fluvial floods will increase in the next 20 years globally by 17% despite partial compensation through market adjustment within the global trade network. China will suffer the strongest direct losses, with an increase of 82%. The United States is mostly affected indirectly through its trade relations. By contrast to the United States, recent intensification of the trade relations with China leaves the European Union better prepared for the import of production losses in the future.
Earth’s surface temperature will continue to rise for another 20 to 30 years even with the strongest carbon emission reduction currently considered. The associated changes in rainfall patterns can result in an increased flood risk worldwide. We compute the required increase in flood protection to keep high-end fluvial flood risk at present levels. The analysis is carried out worldwide for subnational administrative units. Most of the United States, Central Europe, and Northeast and West Africa, as well as large parts of India and Indonesia, require the strongest adaptation effort. More than half of the United States needs to at least double their protection within the next two decades. Thus, the need for adaptation to increased river flood is a global problem affecting industrialized regions as much as developing countries.
Climate change affects societies across the globe in various ways. In addition to gradual changes in temperature and other climatic variables, global warming is likely to increase intensity and frequency of extreme weather events.
Beyond biophysical impacts, these also directly affect societal and economic activity. Additionally, indirect effects can occur; spatially, economic losses can spread along global supply-chains; temporally, climate impacts can change the economic development trajectory of countries.
This thesis first examines how climate change alters river flood risk and its local socio-economic implications. Then, it studies the global economic response to river floods in particular, and to climate change in general.
Changes in high-end river flood risk are calculated for the next three decades on a global scale with high spatial resolution. In order to account for uncertainties, this assessment makes use of an ensemble of climate and hydrological models as well as a river routing model, that is found to perform well regarding peak river discharge. The results show an increase in high-end flood risk in many parts of the world, which require profound adaptation efforts. This pressure to adapt is measured as the enhancement in protection level necessary to stay at historical high-end risk. In developing countries as well as in industrialized regions, a high pressure to adapt is observed - the former to increase low protection levels, the latter to maintain the low risk levels perceived in the past.
Further in this thesis, the global agent-based dynamic supply-chain model acclimate is developed. It models the cascading of indirect losses in the global supply network. As an anomaly model its agents - firms and consumers - maximize their profit locally to respond optimally to local perturbations. Incorporating quantities as well as prices on a daily basis, it is suitable to dynamically resolve the impacts of unanticipated climate extremes.
The model is further complemented by a static measure, which captures the inter-dependencies between sectors across regions that are only connected indirectly. These higher-order dependencies are shown to be important for a comprehensive assessment of loss-propagation and overall costs of local disasters.
In order to study the economic response to river floods, the acclimate model is driven by flood simulations. Within the next two decades, the increase in direct losses can only partially be compensated by market adjustments, and total losses are projected to increase by 17% without further adaptation efforts. The US and the EU are both shown to receive indirect losses from China, which is strongly affected directly. However, recent trends in the trade relations leave the EU in a better position to compensate for these losses.
Finally, this thesis takes a broader perspective when determining the investment response to the climate change damages employing the integrated assessment model DICE. On an optimal economic development path, the increase in damages is anticipated as emissions and consequently temperatures increase. This leads to a significant devaluation of investment returns and the income losses from climate damages almost double.
Overall, the results highlight the need to adapt to extreme weather events - local physical adaptation measures have to be combined with regional and global policy measures to prepare the global supply-chain network to climate change.
We investigate the ultrafast magnetization dynamics of FePt in the L1(0) phase after an optical heating pulse, as used in heat-assisted magnetic recording. We compare continuous and nano-granular thin films and emphasize the impact of the finite size on the remagnetization dynamics. The remagnetization speeds up significantly with increasing external magnetic field only for the continuous film, where domain-wall motion governs the dynamics. The ultrafast remagnetization dynamics in the continuous film are only dominated by heat transport in the regime of high magnetic fields, whereas the timescale required for cooling is prevalent in the granular film for all magnetic field strengths. These findings highlight the necessary conditions for studying the intrinsic heat transport properties in magnetic materials.
Ultrafast magnetisation dynamics have been investigated intensely for two decades. The recovery process after demagnetisation, however, was rarely studied experimentally and discussed in detail. The focus of this work lies on the investigation of the magnetisation on long timescales after laser excitation. It combines two ultrafast time resolved methods to study the relaxation of the magnetic and lattice system after excitation with a high fluence ultrashort laser pulse. The magnetic system is investigated by time resolved measurements of the magneto-optical Kerr effect. The experimental setup has been implemented in the scope of this work. The lattice dynamics were obtained with ultrafast X-ray diffraction. The combination of both techniques leads to a better understanding of the mechanisms involved in magnetisation recovery from a non-equilibrium condition. Three different groups of samples are investigated in this work: Thin Nickel layers capped with nonmagnetic materials, a continuous sample of the ordered L10 phase of Iron Platinum and a sample consisting of Iron Platinum nanoparticles embedded in a carbon matrix. The study of the remagnetisation reveals a general trend for all of the samples: The remagnetisation process can be described by two time dependences. A first exponential recovery that slows down with an increasing amount of energy absorbed in the system until an approximately linear time dependence is observed. This is followed by a second exponential recovery. In case of low fluence excitation, the first recovery is faster than the second. With increasing fluence the first recovery is slowed down and can be described as a linear function. If the pump-induced temperature increase in the sample is sufficiently high, a phase transition to a paramagnetic state is observed. In the remagnetisation process, the transition into the ferromagnetic state is characterised by a distinct transition between the linear and exponential recovery. From the combination of the transient lattice temperature Tp(t) obtained from ultrafast X-ray measurements and magnetisation M(t) gained from magneto-optical measurements we construct the transient magnetisation versus temperature relations M(Tp). If the lattice temperature remains below the Curie temperature the remagnetisation curve M(Tp) is linear and stays below the M(T) curve in equilibrium in the continuous transition metal layers. When the sample is heated above phase transition, the remagnetisation converges towards the static temperature dependence. For the granular Iron Platinum sample the M(Tp) curves for different fluences coincide, i.e. the remagnetisation follows a similar path irrespective of the initial laser-induced temperature jump.
We have investigated the structural dynamics in photoexcited 1,2-diiodotetrafluoroethane molecules (C2F4I2) in the gas phase experimentally using ultrafast electron diffraction and theoretically using FOMO-CASCI excited-state dynamics simulations. The molecules are excited by an ultraviolet femtosecond laser pulse to a state characterized by a transition from the iodine 5p perpendicular to orbital to a mixed 5p parallel to sigma hole and CF2 center dot antibonding orbital, which results in the cleavage of one of the carbon-iodine bonds. We have observed, with sub-Angstrom resolution, the motion of the nuclear wave packet of the dissociating iodine atom followed by coherent vibrations in the electronic ground state of the C2F4I radical. The radical reaches a stable classical (nonbridged) structure in less than 200 fs.
Concurrent observation technologies have made high-precision real-time data available in large quantities. Data assimilation (DA) is concerned with how to combine this data with physical models to produce accurate predictions. For spatial-temporal models, the ensemble Kalman filter with proper localisation techniques is considered to be a state-of-the-art DA methodology. This article proposes and investigates a localised ensemble Kalman Bucy filter for nonlinear models with short-range interactions. We derive dimension-independent and component-wise error bounds and show the long time path-wise error only has logarithmic dependence on the time range. The theoretical results are verified through some simple numerical tests.
Concurrent observation technologies have made high-precision real-time data available in large quantities. Data assimilation (DA) is concerned with how to combine this data with physical models to produce accurate predictions. For spatial-temporal models, the ensemble Kalman filter with proper localisation techniques is considered to be a state-of-the-art DA methodology. This article proposes and investigates a localised ensemble Kalman Bucy filter for nonlinear models with short-range interactions. We derive dimension-independent and component-wise error bounds and show the long time path-wise error only has logarithmic dependence on the time range. The theoretical results are verified through some simple numerical tests.
'Complex systems are information processors' is a statement that is frequently made. Here we argue for the distinction between information processing-in the sense of encoding and transmitting a symbolic representation-and the formation of correlations (pattern formation/self-organisation). The study of both uses tools from information theory, but the purpose is very different in each case: explaining the mechanisms and understanding the purpose or function in the first case, versus data analysis and correlation extraction in the latter. We give examples of both and discuss some open questions. The distinction helps focus research efforts on the relevant questions in each case.
The origin of cosmic rays was the subject of several studies for over a century. The investigations done within this dissertation are one small step to shed some more light on this mystery.
Locating the sources of cosmic rays is not trivial due to the interstellar magnetic field. However, the Hillas criterion allows us to arrive at the conclusion that supernova remnants are our main suspect for the origin of galactic cosmic rays. The mechanism by which they are accelerating particles is found within the field of shock physics as diffusive shock acceleration. To allow particles to enter this process also known as Fermi acceleration pre-acceleration processes like shock surfing acceleration and shock drift acceleration are necessary. Investigating the processes happening in the plasma shocks of supernova remnants is possible by utilising a simplified model which can be simulated on a computer using Particle-in-Cell simulations.
We developed a new and clean setup to simulate the formation of a double shock, i.e., consisting of a forward and a reverse shock and a contact discontinuity, by the collision of two counter-streaming plasmas, in which a magnetic field can be woven into. In a previous work, we investigated the processes at unmagnetised and at magnetised parallel shocks, whereas in the current work, we move our investigation on to magnetised perpendicular shocks.
Due to a much stronger confinement of the particles to the collision region the perpendicular shock develops much faster than the parallel shock. On the other hand, this leads to much weaker turbulence. We are able to find indications for shock surfing acceleration and shock drift acceleration happening at the two shocks leading to populations of pre-accelerated particles that are suitable as a seed population to be injected into further diffusive shock acceleration to be accelerated to even higher energies. We observe the development of filamentary structures in the shock ramp of the forward shock, but not at the reverse shock. This leads to the conclusion that the development of such structures in the shock ramp of quasi-perpendicular collisionless shocks might not necessarily be determined by the existence of a critical sonic Mach number but by a critical shock speed.
The results of the investigations done within this dissertation might be useful for further studies of oblique shocks and for studies using hybrid or magnetohydrodynamic simulations. Together with more sophisticated observational methods, these studies will help to bring us closer to an answer as to how particles can be accelerated in supernova remnants and eventually become cosmic rays that can be detected on Earth.
Ground ice and sedimentary records of a pingo exposure reveal insights into Holocene permafrost, landscape and climate dynamics. Early to mid-Holocene thermokarst lake deposits contain rich floral and faunal paleoassemblages, which indicate lake shrinkage and decreasing summer temperatures (chironomid-based T-July) from 10.5 to 3.5 cal kyr BP with the warmest period between 10.5 and 8 cal kyr BP. Talik refreezing and pingo growth started about 3.5 cal kyr BP after disappearance of the lake. The isotopic composition of the pingo ice (delta O-18 - 17.1 +/- 0.6 parts per thousand, delta D -144.5 +/- 3.4 parts per thousand, slope 5.85, deuterium excess -7.7 +/- 1.5 parts per thousand) point to the initial stage of closed-system freezing captured in the record. A differing isotopic composition within the massive ice body was found (delta O-18 - 21.3 +/- 1.4 parts per thousand, delta D -165 +/- 11.5 parts per thousand, slope 8.13, deuterium excess 4.9 +/- 3.2 parts per thousand), probably related to the infill of dilation cracks by surface water with quasi-meteoric signature. Currently inactive syngenetic ice wedges formed in the thermokarst basin after lake drainage. The pingo preserves traces of permafrost response to climate variations in terms of ground-ice degradation (thermokarst) during the early and mid-Holocene, and aggradation (wedge-ice and pingo-ice growth) during the late Holocene.
After the United Kingdom has left the European Union it remains unclear whether the two parties can successfully negotiate and sign a trade agreement within the transition period. Ongoing negotiations, practical obstacles and resulting uncertainties make it highly unlikely that economic actors would be fully prepared to a “no-trade-deal” situation. Here we provide an economic shock simulation of the immediate aftermath of such a post-Brexit no-trade-deal scenario by computing the time evolution of more than 1.8 million interactions between more than 6,600 economic actors in the global trade network. We find an abrupt decline in the number of goods produced in the UK and the EU. This sudden output reduction is caused by drops in demand as customers on the respective other side of the Channel incorporate the new trade restriction into their decision-making. As a response, producers reduce prices in order to stimulate demand elsewhere. In the short term consumers benefit from lower prices but production value decreases with potentially severe socio-economic consequences in the longer term.
After the United Kingdom has left the European Union it remains unclear whether the two parties can successfully negotiate and sign a trade agreement within the transition period. Ongoing negotiations, practical obstacles and resulting uncertainties make it highly unlikely that economic actors would be fully prepared to a “no-trade-deal” situation. Here we provide an economic shock simulation of the immediate aftermath of such a post-Brexit no-trade-deal scenario by computing the time evolution of more than 1.8 million interactions between more than 6,600 economic actors in the global trade network. We find an abrupt decline in the number of goods produced in the UK and the EU. This sudden output reduction is caused by drops in demand as customers on the respective other side of the Channel incorporate the new trade restriction into their decision-making. As a response, producers reduce prices in order to stimulate demand elsewhere. In the short term consumers benefit from lower prices but production value decreases with potentially severe socio-economic consequences in the longer term.
We investigate whether the dust content of the circum-galactic medium (CGM) depends on the location of the quasar sightline with respect to the galaxy major-axis using 13 galaxy-Mg II absorber pairs (9-81 kpc distance) from the MusE GAs FLOw and Wind (MEGAFLOW) survey at 0.4 < z < 1.4. The dust content of the CGM is obtained from [Zn/Fe] using ultraviolet and visual echelle spectrograph data. When a direct measurement of [Zn/Fe] is unavailable, we estimate the dust depletion from a method that consists in solving for the depletion from multiple singly ionized ions (e.g. Mn II, Cr II, and Zn II) since each ion depletes on dust grains at different rates. We find a positive correlation between the azimuthal angle and [Zn/Fe] with a Pearson's gamma = 0.70 +/- 0.14. The sightlines along the major axis show [Zn/Fe] < 0.5, whereas the [Zn/Fe] is > 0.8 along the minor axis. These results suggest that the CGM along the minor axis is on average more metal enriched (by approximate to 1 dex) than the gas located along the major axis of galaxies provided that dust depletion is a proxy for metallicity. This anisotropic distribution is consistent with recent results on outflow and accretion in hydro-dynamical simulations.
Cosmic-Ray neutron sensors are widely used to determine soil moisture on the hectare scale. Precise measurements, especially in the case of mobile application, demand for neutron detectors with high counting rates and high signal-to-noise ratios. For a long time Cosmic Ray Neutron Sensing (CRNS) instruments have relied on He-3 as an efficient neutron converter. Its ongoing scarcity demands for technological solutions using alternative converters, which are Li-6 and B-10. Recent developments lead to a modular neutron detector consisting of several B-10-lined proportional counter tubes, which feature high counting rates via its large surface area. The modularity allows for individual shieldings of different segments within the detector featuring the capability of gaining spectral information about the detected neutrons. This opens the possibility for active signal correction, especially useful when applied to mobile measurements, where the influence of constantly changing near-field to the overall signal should be corrected. Furthermore, the signal-to-noise ratio could be increased by combining pulse height and pulse length spectra to discriminate between neutrons and other environmental radiation. This novel detector therefore combines high-selective counting electronics with large-scale instrumentation technology.
The Antennae galaxy (NGC 4038/39) is the closest major interacting galaxy system and is therefore often studied as a merger prototype. We present the first comprehensive integral field spectroscopic dataset of this system, observed with the MUSE instrument at the ESO VLT. We cover the two regions in this system which exhibit recent star formation: the central galaxy interaction and a region near the tip of the southern tidal tail. In these fields, we detect H II regions and diffuse ionized gas to unprecedented depth. About 15% of the ionized gas was undetected by previous observing campaigns. This newly detected faint ionized gas is visible everywhere around the central merger, and shows filamentary structure. We estimate diffuse gas fractions of about 60% in the central field and 10% in the southern region. We are able to show that the southern region contains a significantly different population of H II regions, showing fainter luminosities. By comparing H II region luminosities with the HST catalog of young star clusters in the central field, we estimate that there is enough Lyman-continuum leakage in the merger to explain the amount of diffuse ionized gas that we detect. We compare the Lyman-continuum escape fraction of each H II region against emission line ratios that are sensitive to the ionization parameter. While we find no systematic trend between these properties, the most extreme line ratios seem to be strong indicators of density bounded ionization. Extrapolating the Lyman-continuum escape fractions to the southern region, we conclude that simply from the comparison of the young stellar populations to the ionized gas there is no need to invoke other ionization mechanisms than Lyman-continuum leaking H II regions for the diffuse ionized gas in the Antennae.
Proteins are chain molecules built from amino acids. The precise sequence of the 20 different types of amino acids in a protein chain defines into which structure a protein folds, and the three-dimensional structure in turn specifies the biological function of the protein. The reliable folding of proteins is a prerequisite for their robust function. Misfolding can lead to protein aggregates that cause severe diseases, such as Alzheimer's, Parkinson's, or the variant Creutzfeldt-Jakob disease. Small single-domain proteins often fold without experimentally detectable metastable intermediate states. The folding dynamics of these proteins is thought to be governed by a single transition-state barrier between the unfolded and the folded state. The transition state is highly instable and cannot be observed directly. However, mutations in which a single amino acid of the protein is substituted by another one can provide indirect access. The mutations slightly change the transition-state barrier and, thus, the folding and unfolding times of the protein. The central question is how to reconstruct the transition state from the observed changes in folding times. In this habilitation thesis, a novel method to extract structural information on transition states from mutational data is presented. The method is based on (i) the cooperativity of structural elements such as alpha-helices and beta-hairpins, and (ii) on splitting up mutation-induced free-energy changes into components for these elements. By fitting few parameters, the method reveals the degree of structure formation of alpha-helices and beta-hairpins in the transition state. In addition, it is shown in this thesis that the folding routes of small single-domain proteins are dominated by loop-closure dependencies between the structural elements.
Robotic telescopes & Doppler imaging : measuring differential rotation on long-period active stars
(2004)
The sun shows a wide variety of magnetic-activity related phenomena. The magnetic field responsible for this is generated by a dynamo process which is believed to operate in the tachocline, which is located at the bottom of the convection zone. This dynamo is driven in part by differential rotation and in part by magnetic turbulences in the convection zone. The surface differential rotation, one key ingredient of dynamo theory, can be measured by tracing sunspot positions.To extend the parameter space for dynamo theories, one can extend these measurements to other stars than the sun. The primary obstacle in this endeavor is the lack of resolved surface images on other stars. This can be overcome by the Doppler imaging technique, which uses the rotation-induced Doppler-broadening of spectral lines to compute the surface distribution of a physical parameter like temperature. To obtain the surface image of a star, high-resolution spectroscopic observations, evenly distributed over one stellar rotation period are needed. This turns out to be quite complicated for long period stars. The upcoming robotic observatory STELLA addresses this problem with a dedicated scheduling routine, which is tailored for Doppler imaging targets. This will make observations for Doppler imaging not only easier, but also more efficient.As a preview of what can be done with STELLA, we present results of a Doppler imaging study of seven stars, all of which show evidence for differential rotation, but unfortunately the errors are of the same order of magnitude as the measurements due to unsatisfactory data quality, something that will not happen on STELLA. Both, cross-correlation analysis and the sheared image technique where used to double check the results if possible. For four of these stars, weak anti-solar differential rotation was found in a sense that the pole rotates faster than the equator, for the other three stars weak differential rotation in the same direction as on the sun was found.Finally, these new measurements along with other published measurements of differential rotation using Doppler imaging, were analyzed for correlations with stellar evolution, binarity, and rotation period. The total sample of stars show a significant correlation with rotation period, but if separated into antisolar and solar type behavior, only the subsample showing anti-solar differential rotation shows this correlation. Additionally, there is evidence for binary stars showing less differential rotation as single stars, as is suggested by theory. All other parameter combinations fail to deliver any results due to the still small sample of stars available.
In nature as well as in the context of infection and medical applications, bacteria often have to move in highly complex environments such as soil or tissues. Previous studies have shown that bacteria strongly interact with their surroundings and are often guided by confinements. Here, we investigate theoretically how the dispersal of swimming bacteria can be augmented by microfluidic environments and validate our theoretical predictions experimentally. We consider a system of bacteria performing the prototypical run-and-tumble motion inside a labyrinth with square lattice geometry. Narrow channels between the square obstacles limit the possibility of bacteria to reorient during tumbling events to an area where channels cross. Thus, by varying the geometry of the lattice it might be possible to control the dispersal of cells. We present a theoretical model quantifying diffusive spreading of a run-and-tumble random walker in a square lattice. Numerical simulations validate our theoretical predictions for the dependence of the diffusion coefficient on the lattice geometry. We show that bacteria moving in square labyrinths exhibit enhanced dispersal as compared to unconfined cells. Importantly, confinement significantly extends the duration of the phase with strongly non-Gaussian diffusion, when the geometry of channels is imprinted in the density profiles of spreading cells. Finally, in good agreement with our theoretical findings, we observe the predicted behaviors in experiments with E. coli bacteria swimming in a square lattice labyrinth created in amicrofluidic device. Altogether, our comprehensive understanding of bacterial dispersal in a simple two-dimensional labyrinth makes the first step toward the analysis of more complex geometries relevant for real world applications.
In nature as well as in the context of infection and medical applications, bacteria often have to move in highly complex environments such as soil or tissues. Previous studies have shown that bacteria strongly interact with their surroundings and are often guided by confinements. Here, we investigate theoretically how the dispersal of swimming bacteria can be augmented by microfluidic environments and validate our theoretical predictions experimentally. We consider a system of bacteria performing the prototypical run-and-tumble motion inside a labyrinth with square lattice geometry. Narrow channels between the square obstacles limit the possibility of bacteria to reorient during tumbling events to an area where channels cross. Thus, by varying the geometry of the lattice it might be possible to control the dispersal of cells. We present a theoretical model quantifying diffusive spreading of a run-and-tumble random walker in a square lattice. Numerical simulations validate our theoretical predictions for the dependence of the diffusion coefficient on the lattice geometry. We show that bacteria moving in square labyrinths exhibit enhanced dispersal as compared to unconfined cells. Importantly, confinement significantly extends the duration of the phase with strongly non-Gaussian diffusion, when the geometry of channels is imprinted in the density profiles of spreading cells. Finally, in good agreement with our theoretical findings, we observe the predicted behaviors in experiments with E. coli bacteria swimming in a square lattice labyrinth created in amicrofluidic device. Altogether, our comprehensive understanding of bacterial dispersal in a simple two-dimensional labyrinth makes the first step toward the analysis of more complex geometries relevant for real world applications.
Understanding stars, their magnetic activity phenomena and the underlying dynamo action is the foundation for understanding 'life, the universe and everything' - as stellar magnetic fields play a fundamental role for star and planet formation and for the terrestrial atmosphere and climate. Starspots are the fingerprints of magnetic field lines and thereby the most important sign of activity in a star's photosphere. However, they cannot be observed directly, as it is not (yet) possible to spacially resolve the surfaces of even the nearest neighbouring stars. Therefore, an indirect approach called 'Doppler imaging' is applied, which allows to reconstruct the surface spot distribution on rapidly rotating, active stars. In this work, data from 11 years of continuous spectroscopic observations of the active binary star EI Eridani are reduced and analysed. 34 Doppler maps are obtained and the problem of how to parameterise the information content of Doppler maps is discussed. Three approaches for parameter extraction are introduced and applied to all maps: average temperature, separated for several latitude bands; fractional spottedness; and, for the analysis of structural temperature distribution, longitudinal and latitudinal spot-occurrence functions. The resulting values do not show a distinct correlation with the proposed activity cycle as seen from photometric long-term observations, thereby suggesting that the photometric activity cycle is not accompanied by a spot cycle as seen on the Sun. The general morphology of the spot pattern on EI Eri remains persistent for the whole period of 11 years. In addition, a detailed parameter study is performed. Improved orbital parameters suggest that EI Eri might be complemented by a third star in a wide orbit of about 19 years. Preliminary differential rotation measurements are carried out, indicating an anti-solar orientation.
The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, N(E) proportional to E-p for p greater than or similar to 2. However, in standard Fermi shock acceleration, a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future gamma-ray observatories, such as CTA, very likely.
The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.
The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion.
How does a systematic time-dependence of the diffusion coefficient D(t) affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we answer this question via studying the characteristics of a set of standard statistical quantifiers relevant to single-particle-tracking (SPT) experiments. We examine, for instance, how the behavior of the ensemble- and time-averaged mean-squared displacements-denoted as the standard MSD < x(2)(Delta)> and TAMSD <<(delta(2)(Delta))over bar>> quantifiers-of FBM featuring < x(2) (Delta >> = <<(delta(2)(Delta >)over bar>> proportional to Delta(2H) (where H is the Hurst exponent and Delta is the [lag] time) changes in the presence of a power-law deterministically varying diffusivity D-proportional to(t) proportional to t(alpha-1) -germane to the process of scaled Brownian motion (SBM)-determining the strength of fractional Gaussian noise. The resulting compound "scaled-fractional" Brownian motion or FBM-SBM is found to be nonergodic, with < x(2)(Delta >> proportional to Delta(alpha+)(2H)(-1) and <(delta 2(Delta >) over bar > proportional to Delta(2H). We also detect a stalling behavior of the MSDs for very subdiffusive SBM and FBM, when alpha + 2H - 1 < 0. The distribution of particle displacements for FBM-SBM remains Gaussian, as that for the parent processes of FBM and SBM, in the entire region of scaling exponents (0 < alpha < 2 and 0 < H < 1). The FBM-SBM process is aging in a manner similar to SBM. The velocity autocorrelation function (ACF) of particle increments of FBM-SBM exhibits a dip when the parent FBM process is subdiffusive. Both for sub- and superdiffusive FBM contributions to the FBM-SBM process, the SBM exponent affects the long-time decay exponent of the ACF. Applications of the FBM-SBM-amalgamated process to the analysis of SPT data are discussed. A comparative tabulated overview of recent experimental (mainly SPT) and computational datasets amenable for interpretation in terms of FBM-, SBM-, and FBM-SBM-like models of diffusion culminates the presentation. The statistical aspects of the dynamics of a wide range of biological systems is compared in the table, from nanosized beads in living cells, to chromosomal loci, to water diffusion in the brain, and, finally, to patterns of animal movements.
How do different reset protocols affect ergodicity of a diffusion process in single-particle-tracking experiments? We here address the problem of resetting of an arbitrary stochastic anomalous-diffusion process (ADP) from the general mathematical points of view and assess ergodicity of such reset ADPs for an arbitrary resetting protocol. The process of stochastic resetting describes the events of the instantaneous restart of a particle’s motion via randomly distributed returns to a preset initial position (or a set of those). The waiting times of such resetting events obey the Poissonian, Gamma, or more generic distributions with specified conditions regarding the existence of moments. Within these general approaches, we derive general analytical results and support them by computer simulations for the behavior of the reset mean-squared displacement (MSD), the new reset increment-MSD (iMSD), and the mean reset time-averaged MSD (TAMSD). For parental nonreset ADPs with the MSD(t)∝ tμ we find a generic behavior and a switch of the short-time growth of the reset iMSD and mean reset TAMSDs from ∝ _μ for subdiffusive to ∝ _1 for superdiffusive reset ADPs. The critical condition for a reset ADP that recovers its ergodicity is found to be more general than that for the nonequilibrium stationary state, where obviously the iMSD and the mean TAMSD are equal. The consideration of the new statistical quantifier, the iMSD—as compared to the standard MSD—restores the ergodicity of an arbitrary reset ADP in all situations when the μth moment of the waiting-time distribution of resetting events is finite. Potential applications of these new resetting results are, inter alia, in the area of biophysical and soft-matter systems.
How do different reset protocols affect ergodicity of a diffusion process in single-particle-tracking experiments? We here address the problem of resetting of an arbitrary stochastic anomalous-diffusion process (ADP) from the general mathematical points of view and assess ergodicity of such reset ADPs for an arbitrary resetting protocol. The process of stochastic resetting describes the events of the instantaneous restart of a particle’s motion via randomly distributed returns to a preset initial position (or a set of those). The waiting times of such resetting events obey the Poissonian, Gamma, or more generic distributions with specified conditions regarding the existence of moments. Within these general approaches, we derive general analytical results and support them by computer simulations for the behavior of the reset mean-squared displacement (MSD), the new reset increment-MSD (iMSD), and the mean reset time-averaged MSD (TAMSD). For parental nonreset ADPs with the MSD(t)∝ tμ we find a generic behavior and a switch of the short-time growth of the reset iMSD and mean reset TAMSDs from ∝ _μ for subdiffusive to ∝ _1 for superdiffusive reset ADPs. The critical condition for a reset ADP that recovers its ergodicity is found to be more general than that for the nonequilibrium stationary state, where obviously the iMSD and the mean TAMSD are equal. The consideration of the new statistical quantifier, the iMSD—as compared to the standard MSD—restores the ergodicity of an arbitrary reset ADP in all situations when the μth moment of the waiting-time distribution of resetting events is finite. Potential applications of these new resetting results are, inter alia, in the area of biophysical and soft-matter systems.
Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.
Thanks to the combined effort s of scientist s in several research fields, the preceding decade has witnessed considerable progress in the use of conjugated polymers as emerging thermoelectric materials leading to significant improvements in performance and demonstration of a number of diverse applications. Despite these recent advances, systematic assessments of the impact of molecular design on thermoelectric properties are scarce. Although several reviews marginally highlight the role of chemical structure, the understanding of structure-performance relationships is still fragmented. An in-depth understanding of the relationship between molecular structure and thermoelectric properties will enable the rational design of next-generation thermoelectric polymers. To this end, this review showcases the state-of-the-art thermoelectric polymers, discusses structure-performance relationships, suggests strategies for improving thermoelectric performance that go beyond molecular design, and highlights some of the most impressive applications of thermoelectric polymers.
Inorganic perovskites with cesium (Cs+) as the cation have great potential as photovoltaic materials if their phase purity and stability can be addressed. Herein, a series of inorganic perovskites is studied, and it is found that the power conversion efficiency of solar cells with compositions CsPbI1.8Br1.2, CsPbI2.0Br1.0, and CsPbI2.2Br0.8 exhibits a high dependence on the initial annealing step that is found to significantly affect the crystallization and texture behavior of the final perovskite film. At its optimized annealing temperature, CsPbI1.8Br1.2 exhibits a pure orthorhombic phase and only one crystal orientation of the (110) plane. Consequently, this allows for the best efficiency of up to 14.6% and the longest operational lifetime, T-S80, of approximate to 300 h, averaged of over six solar cells, during the maximum power point tracking measurement under continuous light illumination and nitrogen atmosphere. This work provides essential progress on the enhancement of photovoltaic performance and stability of CsPbI3 - xBrx perovskite solar cells.
Inorganic perovskites with cesium (Cs+) as the cation have great potential as photovoltaic materials if their phase purity and stability can be addressed. Herein, a series of inorganic perovskites is studied, and it is found that the power conversion efficiency of solar cells with compositions CsPbI1.8Br1.2, CsPbI2.0Br1.0, and CsPbI2.2Br0.8 exhibits a high dependence on the initial annealing step that is found to significantly affect the crystallization and texture behavior of the final perovskite film. At its optimized annealing temperature, CsPbI1.8Br1.2 exhibits a pure orthorhombic phase and only one crystal orientation of the (110) plane. Consequently, this allows for the best efficiency of up to 14.6% and the longest operational lifetime, T-S80, of approximate to 300 h, averaged of over six solar cells, during the maximum power point tracking measurement under continuous light illumination and nitrogen atmosphere. This work provides essential progress on the enhancement of photovoltaic performance and stability of CsPbI3 - xBrx perovskite solar cells.
Two new hole selective materials (HSMs) based on dangling methylsulfanyl groups connected to the C-9 position of the fluorene core are synthesized and applied in perovskite solar cells. Being structurally similar to a half of Spiro-OMeTAD molecule, these HSMs (referred as FS and DFS) share similar redox potentials but are endowed with slightly higher hole mobility, due to the planarity and large extension of their structure. Competitive power conversion efficiency (up to 18.6%) is achieved by using the new HSMs in suitable perovskite solar cells. Time-resolved photoluminescence decay measurements and electrochemical impedance spectroscopy show more efficient charge extraction at the HSM/perovskite interface with respect to Spiro-OMeTAD, which is reflected in higher photocurrents exhibited by DFS/FS-integrated perovskite solar cells. Density functional theory simulations reveal that the interactions of methylammonium with methylsulfanyl groups in DFS/FS strengthen their electrostatic attraction with the perovskite surface, providing an additional path for hole extraction compared to the sole presence of methoxy groups in Spiro-OMeTAD. Importantly, the low-cost synthesis of FS makes it significantly attractive for the future commercialization of perovskite solar cells.
Bipolar electrets from polypropylene (PP) are essential, e.g., in electret air filters and in cellular-foam ferroelectrets. Therefore, the mechanism of surface-charge stability enhancement on PP electrets via orthophosphoric-acid surface treatment is investigated in detail. It is shown that the significant charge-stability enhancement can be mainly attributed to deeper surface traps originating from deposited chemicals and topographic features on the modified surfaces. Thermally stimulated discharge of chemically treated and non-treated PP films with different surface-charge densities is used to test the limits of the newly formed deep traps in terms of the capacity for hosting surface charges. When the initial surface-charge density is very high, more charges are forced into shallower original traps on the surface or in the bulk of the treated PP samples, reducing the effect of the deeper surface traps brought by the surface modification. The well-known crossover phenomenon (of the surface-charge decay curves) has been observed between modified PP electrets charged to +/- 2kV and to +/- 3kV. Acoustically probed charge distributions in the thickness direction of PP electrets at different stages of thermal discharging indicate that the deep surface trapping sites may have preference for negative charges, resulting in the observed asymmetric charge stability of the modified PP films.
A double-layer transcrystalline polypropylene (PP) film with a flat central interface layer between its two transcrystalline layers is obtained by recrystallization from the melt between two polytetrafluoroethylene (PTFE) surfaces on both sides of the PP film. Its electret properties are studied and compared with those of a single-layer transcrystalline PP film re-crystallized in contact with only one PTFE surface. Within experimental uncertainty, the two types of transcrystalline films exhibit the same thermal properties and crystallinities. After thermal poling, however, two hetero-charge layers of opposite polarity are found on the internal interfaces of the double-layer transcrystalline films and may together be considered as micrometer-sized dipoles. The unexpected phenomenon does not occur in single-layer transcrystalline samples without a central interface layer, suggesting that the interfaces between the transcrystalline layers and the micrometer-thick central interface layer may be the origin of deeper traps rather than the crystalline structures in the transcrystallites or the spherulites. The origin of the interfacial charges was also studied by means of an injection-blocking charging method, which revealed that intrinsic charge carriers introduced during recrystallization are most likely responsible for the interfacial charges. It is fascinating that a material as familiar as PP can exhibit such intriguing properties with a special bipolar space-charge polarization across the central interface layer after quasi-epitaxial surface moulding into a double-layer transcrystalline form. In addition to applications in electret (micro-)devices for electro-mechanical transduction, the highly ordered structures may also be employed as a new paradigm for studying charge storage and transport in polymer electrets and in dielectrics for DC electrical insulation.
Previous work has shown that surface modification with orthophosphoric acid can significantly enhance the charge stability on polypropylene (PP) surface by generating deeper traps. In the present study, thermally stimulated potential-decay measurements revealed that the chemical treatment may also significantly increase the number of available trapping sites on the surface. Thus, as a consequence, the so-called "cross-over" phenomenon, which is observed on as-received and thermally treated PP electrets, may be overcome in a certain range of initial charge densities. Furthermore, the discharge behavior of chemically modified samples indicates that charges can be injected from the treated surface into the bulk, and/or charges of opposite polarity can be pulled from the rear electrode into the bulk at elevated temperatures and at the high electric fields that are caused by the deposited charges. In the bulk, a lack of deep traps causes rapid charge decay already in the temperature range around 95 degrees C.
Spherulite-related space-charge electret properties of polypropylene (PP) have been widely discussed in the past decades. In the present paper, a less-common crystalline structure in PP-transcrystalline PP-is studied regarding its electret behavior in comparison with the typical spherulitic morphology. Polarized light microscopy and differential scanning calorimetry were employed to characterize the crystallite types and crystallinities of transcrystalline and spherulitic PP. Their electret functionality is investigated by means of thermally stimulated discharge experiments, where the cross-over phenomenon is observed on transcrystalline PP films, whereas surface-potential saturation and undercharging on the surface occur on the spherulitic samples. Besides, an asymmetrical behavior of positive and negative surface-charge stabilities is found on PP with spherulites, the negatively charged spherulitic surfaces show a better charge stability. It is shown that PP electrets are very sensitive to changes in the microscopic crystalline structures and their interfaces as well as in the molecular conformations controlled through adjustments of the respective processing steps. In addition, surface and bulk nanocomposites of PP or low-density polyethylene with inorganic particles are included in the comparison. In view of recent developments in the areas of PP-based electret-fiber filters and cellular-foam ferroelectrets, the observed changes in the charge-storage properties may have particular relevance, as the required film, fiber, or foam processing might significantly modify crystalline morphologies and nano-scale interfaces in PP electrets. Limitations in the charge-storage capabilities of interface structures may also be of interest in the context of high-voltage electrical-insulation materials where reduced space-charge accumulation and slightly increased charge transport can be advantageous.
As one of the most-produced commodity polymers, polypropylene draws considerable scientific and commercial interest as an electret material. In the present thesis, the influence of the surface chemical modification and crystalline reconstruction on the electret properties of the polypropylene thin films will be discussed. The chemical treatment with orthophosphoric acid can significantly improve the surface charge stability of the polypropylene electrets by introducing phosphorus- and oxygen-containing structures onto the modified surface. The thermally stimulated discharge measurement and charge profiling by means of piezoelectrically generated pressure steps are used to investigate the electret behaviour. It is concluded that deep traps of limited number density are created during the treatment with inorganic chemicals. Hence, the improvement dramatically decreases when the surface-charge density is substantially higher than ±1.2×10^(-3) C·m^(-2). The newly formed traps also show a higher trapping energy for negative charges. The energetic distributions of the traps in the non-treated and chemically treated samples offer an insight regarding the surface and foreign-chemical dominance on the charge storage and transport in the polypropylene electrets.
Additionally, different electret properties are observed on the polypropylene films with the spherulitic and transcrystalline structures. It indicates the dependence of the charge storage and transport on the crystallite and molecular orientations in the crystalline phase. In general, a more diverse crystalline growth in the spherulitic samples can result in a more complex energetic trap distribution, in comparison to that in a transcrystalline polypropylene. The double-layer transcrystalline polypropylene film with a crystalline interface in the middle can be obtained by crystallising the film in contact with rough moulding surfaces on both sides. A layer of heterocharges appears on each side of the interface in the double-layer transcrystalline polypropylene electrets after the thermal poling. However, there is no charge captured within the transcrystalline layers. The phenomenon reveals the importance of the crystalline interface in terms of creating traps with the higher activation energy in polypropylene. The present studies highlight the fact that even slight variations in the polypropylene film may lead to dramatic differences in its electret properties.
Tailoring the secondary surface morphology of electro-spun nanofibers has been highly desired, as such delicate structures equip nanofibers with distinct functions. Here, we report a simple strategy to directly reconstruct the surface of polyvinyl alcohol/polyvinylpyrrolidone (PVA/PVP) nanofibers by water evaporation. The roughness and diameter of the nanofibers depend on the temperature during vacuum drying. Surface changes of the nanofibers from smooth to rough were observed at 55 degrees C, with a significant drop in nanofiber diameter. We attribute the formation of the secondary surface morphology to the intermolecular forces in the water vapor, including capillary and the compression forces, on the basis of the results from the Fourier-transform infrared (FTIR) and X-ray photoelectron (XPS) spectroscopy. The strategy is universally effective for various electro-spun polymer nanofibers, thus opening up avenues toward more detailed and sophisticated structure design and implementation for nanofibers.
Chorus waves play an important role in the dynamic evolution of energetic electrons in the Earth's radiation belts and ring current. Using more than 5 years of Van Allen Probe data, we developed a new analytical model for upper‐band chorus (UBC; 0.5fce < f < fce) and lower‐band chorus (LBC; 0.05fce < f < 0.5fce) waves, where fce is the equatorial electron gyrofrequency. By applying polynomial fits to chorus wave root mean square amplitudes, we developed regression models for LBC and UBC as a function of geomagnetic activity (Kp), L, magnetic latitude (λ), and magnetic local time (MLT). Dependence on Kp is separated from the dependence on λ, L, and MLT as Kp‐scaling law to simplify the calculation of diffusion coefficients and inclusion into particle tracing codes. Frequency models for UBC and LBC are also developed, which depends on MLT and magnetic latitude. This empirical model is valid in all MLTs, magnetic latitude up to 20°, Kp ≤ 6, L‐shell range from 3.5 to 6 for LBC and from 4 to 6 for UBC. The dependence of root mean square amplitudes on L are different for different bands, which implies different energy sources for different wave bands. This analytical chorus wave model is convenient for inclusion in quasi‐linear diffusion calculations of electron scattering rates and particle simulations in the inner magnetosphere, especially for the newly developed four‐dimensional codes, which require significantly improved wave parameterizations.
Different modeling methodologies possess different strengths and weakness. For instance, data based models may provide superior accuracy but have a limited spatial coverage while physics based models may provide lower accuracy but provide greater spatial coverage. This study investigates the coupling of a data based model of the electron fluxes at geostationary orbit (GEO) with a numerical model of the radiation belt region to improve the resulting forecasts/pastcasts of electron fluxes over the whole radiation belt region. In particular, two coupling methods are investigated. The first assumes an average value for L* for GEO, namely LGEO* L-GEO* = 6.2. The second uses a value of L* that varies with geomagnetic activity, quantified using the Kp index. As the terrestrial magnetic field responds to variations in geomagnetic activity, the value of L* will vary for a specific location. In this coupling method, the value of L* is calculated using the Kp driven Tsyganenko 89c magnetic field model for field line tracing. It is shown that this addition can result in changes in the initialization of the parameters at the Versatile Electron Radiation Belt model outer boundary. Model outputs are compared to Van Allen Probes MagEIS measurements of the electron fluxes in the inner magnetosphere for the March 2015 geomagnetic storm. It is found that the fixed LGEO* L-GEO* coupling method produces a more realistic forecast.
The maximum entropy method is used to predict flows on water distribution networks. This analysis extends the water distribution network formulation of Waldrip et al. (2016) Journal of Hydraulic Engineering (ASCE), by the use of a continuous relative entropy defined on a reduced parameter set. This reduction in the parameters that the entropy is defined over ensures consistency between different representations of the same network. The performance of the proposed reduced parameter method is demonstrated with a one-loop network case study.
The maximum entropy method is used to derive an alternative gravity model for a transport network. The proposed method builds on previous methods which assign the discrete value of a maximum entropy distribution to equal the traffic flow rate. The proposed method however, uses a distribution to represent each flow rate. The proposed method is shown to be able to handle uncertainty in a more elegant way and give similar results to traditional methods. It is able to incorporate more of the observed data through the entropy function, prior distribution and integration limits potentially allowing better inferences to be made.
A maximum entropy (MaxEnt) method is developed to predict flow rates or pressure gradients in hydraulic pipe networks without sufficient information to give a closed-form (deterministic) solution. This methodology substantially extends existing deterministic flow network analysis methods. It builds on the MaxEnt framework previously developed by the authors. This study uses a continuous relative entropy defined on a reduced parameter set, here based on the external flow rates. This formulation ensures consistency between different representations of the same network. The relative entropy is maximized subject to observable constraints on the mean values of a subset of flow rates or potential differences, the frictional properties of each pipe, and physical constraints arising from Kirchhoff’s first and second laws. The new method is demonstrated by application to a simple one-loop network and a 1,123-node, 1,140-pipe water distribution network in the suburb of Torrens, Australian Capital Territory, Australia.
A task-based parallel elliptic solver for numerical relativity with discontinuous Galerkin methods
(2022)
Elliptic partial differential equations are ubiquitous in physics. In numerical relativity---the study of computational solutions to the Einstein field equations of general relativity---elliptic equations govern the initial data that seed every simulation of merging black holes and neutron stars. In the quest to produce detailed numerical simulations of these most cataclysmic astrophysical events in our Universe, numerical relativists resort to the vast computing power offered by current and future supercomputers. To leverage these computational resources, numerical codes for the time evolution of general-relativistic initial value problems are being developed with a renewed focus on parallelization and computational efficiency. Their capability to solve elliptic problems for accurate initial data must keep pace with the increasing detail of the simulations, but elliptic problems are traditionally hard to parallelize effectively.
In this thesis, I develop new numerical methods to solve elliptic partial differential equations on computing clusters, with a focus on initial data for orbiting black holes and neutron stars. I develop a discontinuous Galerkin scheme for a wide range of elliptic equations, and a stack of task-based parallel algorithms for their iterative solution. The resulting multigrid-Schwarz preconditioned Newton-Krylov elliptic solver proves capable of parallelizing over 200 million degrees of freedom to at least a few thousand cores, and already solves initial data for a black hole binary about ten times faster than the numerical relativity code SpEC. I also demonstrate the applicability of the new elliptic solver across physical disciplines, simulating the thermal noise in thin mirror coatings of interferometric gravitational-wave detectors to unprecedented accuracy. The elliptic solver is implemented in the new open-source SpECTRE numerical relativity code, and set up to support simulations of astrophysical scenarios for the emerging era of gravitational-wave and multimessenger astronomy.
Using a special technique of data analysis, we have found out 34 grand minima of solar activity obtained from a 7,700 years long Δ14C record. The method used rests on a proper filtering of the Δ14C record and the extrapolation of verifiable results for the later history back in time. Additionally, we use a method of nonlinear dynamics, the recurrence rate, to back up the results. Our findings are not contradictory to the record of solar maxima resp. minima by Eddy [5], but constitute a considerable extension. Hence, it has become possible to look closer at the validity of models. This way, we have tested several models for solar activity, esp. the model of Barnes et al. [1]. There are hints for that the grand minima might solely be driven by the 209 year period found in the Δ14C record.
Control over spin and electronic structure of MoS₂ monolayer via interactions with substrates
(2023)
The molybdenum disulfide (MoS2) monolayer is a semiconductor with a direct bandgap while it is a robust and affordable material.
It is a candidate for applications in optoelectronics and field-effect transistors.
MoS2 features a strong spin-orbit coupling which makes its spin structure promising for acquiring the Kane-Mele topological concept with corresponding applications in spintronics and valleytronics.
From the optical point of view, the MoS2 monolayer features two valleys in the regions of K and K' points. These valleys are differentiated by opposite spins and a related valley-selective circular dichroism.
In this study we aim to manipulate the MoS2 monolayer spin structure in the vicinity of the K and K' points to explore the possibility of getting control over the optical and electronic properties.
We focus on two different substrates to demonstrate two distinct routes: a gold substrate to introduce a Rashba effect and a graphene/cobalt substrate to introduce a magnetic proximity effect in MoS2.
The Rashba effect is proportional to the out-of-plane projection of the electric field gradient. Such a strong change of the electric field occurs at the surfaces of a high atomic number materials and effectively influence conduction electrons as an in-plane magnetic field. A molybdenum and a sulfur are relatively light atoms, thus, similar to many other 2D materials, intrinsic Rashba effect in MoS2 monolayer is vanishing small. However, proximity of a high atomic number substrate may enhance Rashba effect in a 2D material as it was demonstrated for graphene previously.
Another way to modify the spin structure is to apply an external magnetic field of high magnitude (several Tesla), and cause a Zeeman splitting, the conduction electrons.
However, a similar effect can be reached via magnetic proximity which allows us to reduce external magnetic fields significantly or even to zero. The graphene on cobalt interface is ferromagnetic and stable for MoS2 monolayer synthesis. Cobalt is not the strongest magnet; therefore, stronger magnets may lead to more significant results.
Nowadays most experimental studies on the dichalcogenides (MoS2 included) are performed on encapsulated heterostructures that are produced by mechanical exfoliation.
While mechanical exfoliation (or scotch-tape method) allows to produce a huge variety of structures, the shape and the size of the samples as well as distance between layers in heterostructures are impossible to control reproducibly.
In our study we used molecular beam epitaxy (MBE) methods to synthesise both MoS2/Au(111) and MoS2/graphene/Co systems.
We chose to use MBE, as it is a scalable and reproducible approach, so later industry may adapt it and take over.
We used graphene/cobalt instead of just a cobalt substrate because direct contact of MoS2\ monolayer and a metallic substrate may lead to photoluminescence (PL) quenching in the metallic substrate. Graphene and hexagonal boron nitride monolayer are considered building blocks of a new generation of electronics also commonly used as encapsulating materials for PL studies. Moreover graphene is proved to be a suitable substrate for the MBE growth of transitional metal dichalcogenides (TMDCs).
In chapter 1,
we start with an introduction to TMDCs. Then we focus on MoS2 monolayer state of the art research in the fields of application scenario; synthesis approaches; electronic, spin, and optical properties; and interactions with magnetic fields and magnetic materials.
We briefly touch the basics of magnetism in solids and move on to discuss various magnetic exchange interactions and magnetic proximity effect.
Then we describe MoS2 optical properties in more detail. We start from basic exciton physics and its manifestation in the MoS2 monolayer. We consider optical selection rules in the MoS2 monolayer and such properties as chirality, spin-valley locking, and coexistence of bright and dark excitons.
Chapter 2 contains an overview of the employed surface science methods: angle-integrated, angle-resolved, and spin-resolved photoemission; low energy electron diffraction and scanning tunneling microscopy.
In chapter 3, we describe MoS2 monolayer synthesis details for two substrates: gold monocrystal with (111) surface and graphene on cobalt thin film with Co(111) surface orientation.
The synthesis descriptions are followed by a detailed characterisation of the obtained structures: fingerprints of MoS2 monolayer formation; MoS2 monolayer symmetry and its relation to the substrate below; characterisation of MoS2 monolayer coverage, domain distribution, sizes and shapes, and moire structures.
In chapter~4, we start our discussion with MoS2/Au(111) electronic and spin structure. Combining density functional theory computations (DFT) and spin-resolved photoemission studies, we demonstrate that the MoS2 monolayer band structure features an in-plane Rashba spin splitting. This confirms the possibility of MoS2 monolayer spin structure manipulation via a substrate.
Then we investigate the influence of a magnetic proximity in the MoS2/graphene/Co system on the MoS2 monolayer spin structure.
We focus our investigation on MoS2 high symmetry points: G and K.
First, using spin-resolved measurements, we confirm that electronic states are spin-split at the G point via a magnetic proximity effect. Second, combining spin-resolved measurements and DFT computations for MoS2 monolayer in the K point region, we demonstrate the appearance of a small in-plane spin polarisation in the valence band top and predict a full in-plane spin polarisation for the conduction band bottom.
We move forward discussing how these findings are related to the MoS2 monolayer optical properties, in particular the possibility of dark exciton observation. Additionally, we speculate on the control of the MoS2 valley energy via magnetic proximity from cobalt.
As graphene is spatially buffering the MoS2 monolayer from the Co thin film, we speculate on the role of graphene in the magnetic proximity transfer by replacing graphene with vacuum and other 2D materials in our computations.
We finish our discussion by investigating the K-doped MoS2/graphene/Co system and the influence of this doping on the electronic and spin structure as well as on the magnetic proximity effect.
In summary, using a scalable MBE approach we synthesised
MoS2/Au(111) and MoS2/graphene/Co systems. We found a Rashba effect taking place in MoS2/Au(111) which proves that the MoS2 monolayer in-plane spin structure can be modified. In MoS2/graphene/Co the in-plane magnetic proximity effect indeed takes place which rises the possibility of fine tuning the MoS2 optical properties via manipulation of the the substrate magnetisation.
Recently, Nocera and co-workers (J. Am. Chem. Soc. 2018, 140, 13711) demonstrated that triaryl borate Lewis acids facilitate the direct electrochemical reduction of triphenylphosphine oxide (TPPO) to triphenylphosphine (TPP). In the present contribution, we report a quantum chemical study unravelling details of the reaction, which also supports the proposed ECrECi mechanism. Alternative electrochemical routes to TPPO reduction facilitated by other Lewis acids (CH3+), or by photocatalysis at semiconductor surfaces, are also briefly discussed.
We employ ultrafast X-ray diffraction to compare the lattice dynamics of laser-excited continuous and granular FePt films on MgO (100) substrates. Contrary to recent results on free-standing granular films, we observe in both cases a pronounced and long-lasting out-of-plane expansion. We attribute this discrepancy to the in-plane expansion, which is suppressed by symmetry in continuous films. Granular films on substrates are less constrained and already show a reduced out-of-plane contraction. Via the Poisson effect, out-of-plane contractions drive in-plane expansion and vice versa. Consistently, the granular film exhibits a short-lived out-of-plane contraction driven by ultrafast demagnetization which is followed by a reduced and delayed expansion. From the acoustic reflections of the observed strain waves at the film-substrate interface, we extract a 13% reduction of the elastic constants in thin 10 nm FePt films compared to bulk-like samples. (C) 2018 Author(s).
We present a temperature and fluence dependent Ultrafast X-Ray Diffraction study of a laser-heated antiferromagnetic dysprosium thin film. The loss of antiferromagnetic order is evidenced by a pronounced lattice contraction. We devise a method to determine the energy flow between the phonon and spin system from calibrated Bragg peak positions in thermal equilibrium. Reestablishing the magnetic order is much slower than the cooling of the lattice, especially around the Néel temperature. Despite the pronounced magnetostriction, the transfer of energy from the spin system to the phonons in Dy is slow after the spin-order is lost.
We present a temperature and fluence dependent Ultrafast X-Ray Diffraction study of a laser-heated antiferromagnetic dysprosium thin film. The loss of antiferromagnetic order is evidenced by a pronounced lattice contraction. We devise a method to determine the energy flow between the phonon and spin system from calibrated Bragg peak positions in thermal equilibrium. Reestablishing the magnetic order is much slower than the cooling of the lattice, especially around the Néel temperature. Despite the pronounced magnetostriction, the transfer of energy from the spin system to the phonons in Dy is slow after the spin-order is lost.
Optical excitation of spin-ordered rare earth metals triggers a complex response of the crystal lattice since expansive stresses from electron and phonon excitations compete with a contractive stress induced by spin disorder. Using ultrafast x-ray diffraction experiments, we study the layer specific strain response of a dysprosium film within a metallic heterostructure upon femtosecond laser-excitation. The elastic and diffusive transport of energy to an adjacent, non-excited detection layer clearly separates the contributions of strain pulses and thermal excitations in the time domain. We find that energy transfer processes to magnetic excitations significantly modify the observed conventional bipolar strain wave into a unipolar pulse. By modeling the spin system as a saturable energy reservoir that generates substantial contractive stress on ultrafast timescales, we can reproduce the observed strain response and estimate the time- and space dependent magnetic stress. The saturation of the magnetic stress contribution yields a non-monotonous total stress within the nanolayer, which leads to unconventional picosecond strain pulses.
Optical excitation of spin-ordered rare earth metals triggers a complex response of the crystal lattice since expansive stresses from electron and phonon excitations compete with a contractive stress induced by spin disorder. Using ultrafast x-ray diffraction experiments, we study the layer specific strain response of a dysprosium film within a metallic heterostructure upon femtosecond laser-excitation. The elastic and diffusive transport of energy to an adjacent, non-excited detection layer clearly separates the contributions of strain pulses and thermal excitations in the time domain. We find that energy transfer processes to magnetic excitations significantly modify the observed conventional bipolar strain wave into a unipolar pulse. By modeling the spin system as a saturable energy reservoir that generates substantial contractive stress on ultrafast timescales, we can reproduce the observed strain response and estimate the time- and space dependent magnetic stress. The saturation of the magnetic stress contribution yields a non-monotonous total stress within the nanolayer, which leads to unconventional picosecond strain pulses.
Magnetic strain contributions in laser-excited metals studied by time-resolved X-ray diffraction
(2021)
In this work I explore the impact of magnetic order on the laser-induced ultrafast strain response of metals. Few experiments with femto- or picosecond time-resolution have so far investigated magnetic stresses. This is contrasted by the industrial usage of magnetic invar materials or magnetostrictive transducers for ultrasound generation, which already utilize magnetostrictive stresses in the low frequency regime.
In the reported experiments I investigate how the energy deposition by the absorption of femtosecond laser pulses in thin metal films leads to an ultrafast stress generation. I utilize that this stress drives an expansion that emits nanoscopic strain pulses, so called hypersound, into adjacent layers. Both the expansion and the strain pulses change the average inter-atomic distance in the sample, which can be tracked with sub-picosecond time resolution using an X-ray diffraction setup at a laser-driven Plasma X-ray source. Ultrafast X-ray diffraction can also be applied to buried layers within heterostructures that cannot be accessed by optical methods, which exhibit a limited penetration into metals. The reconstruction of the initial energy transfer processes from the shape of the strain pulse in buried detection layers represents a contribution of this work to the field of picosecond ultrasonics.
A central point for the analysis of the experiments is the direct link between the deposited energy density in the nano-structures and the resulting stress on the crystal lattice. The underlying thermodynamical concept of a Grüneisen parameter provides the theoretical framework for my work. I demonstrate how the Grüneisen principle can be used for the interpretation of the strain response on ultrafast timescales in various materials and that it can be extended to describe magnetic stresses. The class of heavy rare-earth elements exhibits especially large magnetostriction effects, which can even lead to an unconventional contraction of the laser-excited transducer material. Such a dominant contribution of the magnetic stress to the motion of atoms has not been demonstrated previously. The observed rise time of the magnetic stress contribution in Dysprosium is identical to the decrease in the helical spin-order, that has been found previously using time-resolved resonant X-ray diffraction. This indicates that the strength of the magnetic stress can be used as a proxy of the underlying magnetic order. Such magnetostriction measurements are applicable even in case of antiparallel or non-collinear alignment of the magnetic moments and a vanishing magnetization.
The strain response of metal films is usually determined by the pressure of electrons and lattice vibrations. I have developed a versatile two-pulse excitation routine that can be used to extract the magnetic contribution to the strain response even if systematic measurements above and below the magnetic ordering temperature are not feasible. A first laser pulse leads to a partial ultrafast demagnetization so that the amplitude and shape of the strain response triggered by the second pulse depends on the remaining magnetic order. With this method I could identify a strongly anisotropic magnetic stress contribution in the magnetic data storage material iron-platinum and identify the recovery of the magnetic order by the variation of the pulse-to-pulse delay. The stark contrast of the expansion of iron-platinum nanograins and thin films shows that the different constraints for the in-plane expansion have a strong influence on the out-of-plane expansion, due to the Poisson effect. I show how such transverse strain contributions need to be accounted for when interpreting the ultrafast out-of-plane strain response using thermal expansion coefficients obtained in near equilibrium conditions.
This work contributes an investigation of magnetostriction on ultrafast timescales to the literature of magnetic effects in materials. It develops a method to extract spatial and temporal varying stress contributions based on a model for the amplitude and shape of the emitted strain pulses. Energy transfer processes result in a change of the stress profile with respect to the initial absorption of the laser pulses. One interesting example occurs in nanoscopic gold-nickel heterostructures, where excited electrons rapidly transport energy into a distant nickel layer, that takes up much more energy and expands faster and stronger than the laser-excited gold capping layer. Magnetic excitations in rare earth materials represent a large energy reservoir that delays the energy transfer into adjacent layers. Such magneto-caloric effects are known in thermodynamics but not extensively covered on ultrafast timescales. The combination of ultrafast X-ray diffraction and time-resolved techniques with direct access to the magnetization has a large potential to uncover and quantify such energy transfer processes.
Poly(vinylidene fluoride-trifluoroethylene) (P(VDF-TrFE)) ferroelectric thin films of different molar ratio have been studied with regard to data memory applications. Therefore, films with thicknesses of 200 nm and less have been spin coated from solution. Observations gained from single layers have been extended to multilayer capacitors and three terminal transistor devices.
Besides conventional hysteresis measurements, the measurement of dielectric non-linearities has been used as a main tool of characterisation. Being a very sensitive and non-destructive method, non-linearity measurements are well suited for polarisation readout and property studies. Samples have been excited using a high quality, single-frequency sinusoidal voltage with an amplitude significantly smaller than the coercive field of the samples. The response was then measured at the excitation frequency and its higher harmonics. Using the measurement results, the linear and non-linear dielectric permittivities ɛ₁, ɛ₂ and ɛ₃ have been determined. The permittivities have been used to derive the temperature-dependent polarisation behaviour as well as the polarisation state and the order of the phase transitions.
The coercive field in VDF-TrFE copolymers is high if compared to their ceramic competitors. Therefore, the film thickness had to be reduced significantly. Considering a switching voltage of 5 V and a coercive field of 50 MV/m, the film thickness has to be 100 nm and below. If the thickness becomes substantially smaller than the other dimensions, surface and interface layer effects become more pronounced. For thicker films of P(VDF-TrFE) with a molar fraction of 56/44 a second-order phase transition without a thermal hysteresis for an ɛ₁(T) temperature cycle has been predicted and observed. This however, could not be confirmed by the measurements of thinner films. A shift of transition temperatures as well as a temperature independent, non-switchable polarisation and a thermal hysteresis for P(VDF-TrFE) 56/44 have been observed. The impact of static electric fields on the polarisation and the phase transition has therefore been studied and simulated, showing that all aforementioned phenomena including a linear temperature dependence of the polarisation might originate from intrinsic electric fields.
In further experiments the knowledge gained from single layer capacitors has been extended to bilayer copolymer thin films of different molar composition. Bilayers have been deposited by succeeding cycles of spin coating from solution. Single layers and their bilayer combination have been studied individually in order to prove the layers stability. The individual layers have been found to be physically stable. But while the bilayers reproduced the main ɛ₁(T) properties of the single layers qualitatively, quantitative numbers could not be explained by a simple serial connection of capacitors. Furthermore, a linear behaviour of the polarisation throughout the measured temperature range has been observed. This was found to match the behaviour predicted considering a constant electric field.
Retention time is an important quantity for memory applications. Hence, the retention behaviour of VDF-TrFE copolymer thin films has been determined using dielectric non-linearities. The polarisation loss in P(VDF-TrFE) poled samples has been found to be less than 20% if recorded over several days. The loss increases significantly if the samples have been poled with lower amplitudes, causing an unsaturated polarisation. The main loss was attributed to injected charges. Additionally, measurements of dielectric non-linearities have been proven to be a sensitive and non-destructive tool to measure the retention behaviour.
Finally, a ferroelectric field effect transistor using mainly organic materials (FerrOFET) has been successfully studied. DiNaphtho[2,3-b:2',3'-f]Thieno[3,2-b]Thiophene (DNTT) has proven to be a stable, suitable organic semiconductor to build up ferroelectric memory devices. Furthermore, an oxidised aluminium bottom electrode and additional dielectric layers, i.e. parylene C, have proven to reduce the leakage current and therefore enhance the performance significantly.
The nonlinear interaction of waves excited by the modified two-stream instability (Farley-Buneman instability) is considered. It is found that, during the linear stage of wave growth, the enhanced pressure of the high-frequency part of the waves locally generates a ponderomotive force. This force acts on the plasma particles and redistributes them. Thus an additional electrostatic polarization field occurs, which influences the low-frequency part of the waves. Then, the low-frequency waves also cause a redistribution of the high-frequency waves. In the paper, a self-consistent system of equations is obtained, which describes the nonlinear interaction of the waves. It is shown that the considered mechanism of wave interaction causes a nonlinear stabilization of the high-frequency waves’ growth and a formation of local density structures of the charged particles. The density modifications of the charged particles during the non-linear stage of wave growth and the possible interval of aspect angles of the high-frequency waves are estimated.
Basing on recent solar models, the excitation of ion-acoustic turbulence in the weaklycollisional, fully and partially-ionized regions of the solar atmosphere is investigated. Within the frame of hydrodynamics, conditions are found under which the heating of the plasma by ion-acoustic type waves is more effective than the Joule heating. Taking into account wave and Joule heating effects, a nonlinear differential equation is derived, which describes the evolution of nonlinear ion-acoustic waves in the collisional plasma.
The increase in the performance of organic solar cells observed over the past few years has reinvigorated the search for a deeper understanding of the loss and extraction processes in this class of device. A detailed knowledge of the density of free charge carriers under different operating conditions and illumination intensities is a prerequisite to quantify the recombination and extraction dynamics. Differential charging techniques are a promising approach to experimentally obtain the charge carrier density under the aforementioned conditions. In particular, the combination of transient photovoltage and photocurrent as well as impedance and capacitance spectroscopy have been successfully used in past studies to determine the charge carrier density of organic solar cells. In this Tutorial, these experimental techniques will be discussed in detail, highlighting fundamental principles, practical considerations, necessary corrections, advantages, drawbacks, and ultimately their limitations. Relevant references introducing more advanced concepts will be provided as well. Therefore, the present Tutorial might act as an introduction and guideline aimed at new prospective users of these techniques as well as a point of reference for more experienced researchers. Published under an exclusive license by AIP Publishing.
The determination of the voltage-dependent density of free charge carriers via capacitance spectroscopy is considered an important step in the analysis of emerging photovoltaic technologies, such as organic and perovskite solar cells. In particular, an intimate knowledge of the density of free charge carriers is required for the determination of crucial parameters such as the effective mobility, charge carrier lifetime, nongeminate recombination coefficients, average extraction times, and competition factors. Hence, it is paramount to verify the validity of the commonly employed approaches to obtain the density of free charge carriers. The advantages, drawbacks, and limitations of the most common approaches are investigated in detail and strategies to mitigate misleading values are explored. To this end, two types of nonfullerene organic solar cells based on a PTB7-Th:ITIC-2F blend and a PM6:Y6 blend, respectively, are used as a case study to assess how subsequent analyses of the nongeminate recombination dynamics depend on the chosen approach to calculate the density of free charge carriers via capacitance spectroscopy.
Understanding the recombination dynamics of organic and perovskite solar cells has been a crucial prerequisite in the steadily increasing performance of these promising new types of photovoltaics. Surface recombination in particular has turned out to be one of the last remaining roadblocks, which specifically reduces the open circuit voltage. In this study, the relationship between the rate of surface recombination and the density of charge carriers is analyzed, revealing a cubic dependence between these two parameters. This hypothesis is then tested and verified with the recombination dynamics of an organic solar cell known to exhibit significant surface recombination and a high energy proton irradiated CH3NH3PbI3 pemvskite solar cell during white light illumination. Incidentally, these results can also explain recombination orders exceeding the commonly known threshold for bimolecular recombination that have been observed in some studies without the need for a charge carrier dependent bimolecular recombination coefficient.
We investigate anomalous diffusion processes governed by the fractional Langevin equation and confined to a finite or semi-infinite interval by reflecting potential barriers. As the random and damping forces in the fractional Langevin equation fulfill the appropriate fluctuation-dissipation relation, the probability density on a finite interval converges for long times towards the expected uniform distribution prescribed by thermal equilibrium. In contrast, on a semi-infinite interval with a reflecting wall at the origin, the probability density shows pronounced deviations from the Gaussian behavior observed for normal diffusion. If the correlations of the random force are persistent (positive), particles accumulate at the reflecting wall while antipersistent (negative) correlations lead to a depletion of particles near the wall. We compare and contrast these results with the strong accumulation and depletion effects recently observed for nonthermal fractional Brownian motion with reflecting walls, and we discuss broader implications.
Synchronization of large ensembles of oscillators is an omnipresent phenomenon observed in different fields of science like physics, engineering, life sciences, etc. The most simple setup is that of globally coupled phase oscillators, where all the oscillators contribute to a global field which acts on all oscillators. This formulation of the problem was pioneered by Winfree and Kuramoto. Such a setup gives a possibility for the analysis of these systems in terms of global variables. In this work we describe nontrivial collective dynamics in oscillator populations coupled via mean fields in terms of global variables. We consider problems which cannot be directly reduced to standard Kuramoto and Winfree models.
In the first part of the thesis we adopt a method introduced by Watanabe and Strogatz. The main idea is that the system of identical oscillators of particular type can be described by a low-dimensional system of global equations. This approach enables us to perform a complete analytical analysis for a special but vast set of initial conditions. Furthermore, we show how the approach can be expanded for some nonidentical systems. We apply the Watanabe-Strogatz approach to arrays of Josephson junctions and systems of identical phase oscillators with leader-type coupling.
In the next parts of the thesis we consider the self-consistent mean-field theory method that can be applied to general nonidentical globally coupled systems of oscillators both with or without noise. For considered systems a regime, where the global field rotates uniformly, is the most important one. With the help of this approach such solutions of the self-consistency equation for an arbitrary distribution of frequencies and coupling parameters can be found analytically in the parametric form, both for noise-free and noisy cases.
We apply this method to deterministic Kuramoto-type model with generic coupling and an ensemble of spatially distributed oscillators with leader-type coupling. Furthermore, with the proposed self-consistent approach we fully characterize rotating wave solutions of noisy Kuramoto-type model with generic coupling and an ensemble of noisy oscillators with bi-harmonic coupling.
Whenever possible, a complete analysis of global dynamics is performed and compared with direct numerical simulations of large populations.
The problem of biological motion is a very intriguing and topical issue. Many efforts are being focused on the development of novel modelling approaches for the description of anomalous diffusion in biological systems, such as the very complex and heterogeneous cell environment. Nevertheless, many questions are still open, such as the joint manifestation of statistical features in agreement with different models that can also be somewhat alternative to each other, e.g. continuous time random walk and fractional Brownian motion. To overcome these limitations, we propose a stochastic diffusion model with additive noise and linear friction force (linear Langevin equation), thus involving the explicit modelling of velocity dynamics. The complexity of the medium is parametrized via a population of intensity parameters (relaxation time and diffusivity of velocity), thus introducing an additional randomness, in addition to white noise, in the particle's dynamics. We prove that, for proper distributions of these parameters, we can get both Gaussian anomalous diffusion, fractional diffusion and its generalizations.
We derive. the ensemble-and time-averaged mean-squared displacements (MSD, TAMSD) for Poisson-reset geometric Brownian motion (GBM), in agreement with simulations. We find MSD and TAMSD saturation for frequent resetting, quantify the spread of TAMSDs via the ergodicity-breaking parameter and compute distributions of prices. General MSD-TAMSD nonequivalence proves reset GBM nonergodic.
How do near-bankruptcy events in the past affect the dynamics of stock-market prices in the future? Specifically, what are the long-time properties of a time-local exponential growth of stock-market prices under the influence of stochastically occurring economic crashes? Here, we derive the ensemble- and time-averaged properties of the respective "economic" or geometric Brownian motion (GBM) with a nonzero drift exposed to a Poissonian constant-rate price-restarting process of "resetting." We examine-based both on thorough analytical calculations and on findings from systematic stochastic computer simulations-the general situation of reset GBM with a nonzero [positive] drift and for all special cases emerging for varying parameters of drift, volatility, and reset rate in the model. We derive and summarize all short- and long-time dependencies for the mean-squared displacement (MSD), the variance, and the mean time-averaged MSD (TAMSD) of the process of Poisson-reset GBM under the conditions of both rare and frequent resetting. We consider three main regions of model parameters and categorize the crossovers between different functional behaviors of the statistical quantifiers of this process. The analytical relations are fully supported by the results of computer simulations. In particular, we obtain that Poisson-reset GBM is a nonergodic stochastic process, with generally MSD(Delta) not equal TAMSD(Delta) and Variance(Delta) not equal TAMSD(Delta) at short lag times Delta and for long trajectory lengths T. We investigate the behavior of the ergodicity-breaking parameter in each of the three regions of parameters and examine its dependence on the rate of reset at Delta/T << 1. Applications of these theoretical results to the analysis of prices of reset-containing options are pertinent.
Corona and the climate
(2020)
We study a non-Markovian and nonstationary model of animal mobility incorporating both exploration and memory in the form of preferential returns. Exact results for the probability of visiting a given number of sites are derived and a practical WKB approximation to treat the nonstationary problem is developed. A mean-field version of this model, first suggested by Song et al., [Modelling the scaling properties of human mobility, Nat. Phys. 6, 818 (2010)] was shown to well describe human movement data. We show that our generalized model adequately describes empirical movement data of Egyptian fruit bats (Rousettus aegyptiacus) when accounting for interindividual variation in the population. We also study the probability of visiting any site a given number of times and derive a mean-field equation. Our analysis yields a remarkable phase transition occurring at preferential returns which scale linearly with past visits. Following empirical evidence, we suggest that this phase transition reflects a trade-off between extensive and intensive foraging modes.
Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations (“Joseph effect”), fat-tailed probability density of increments (“Noah effect”), and nonstationarity (“Moses effect”). We show that such a decomposition of real-life data allows us to infer nontrivial behavioral predictions and to resolve open questions in the fields of single-particle tracking in living cells and movement ecology.
Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations (“Joseph effect”), fat-tailed probability density of increments (“Noah effect”), and nonstationarity (“Moses effect”). We show that such a decomposition of real-life data allows us to infer nontrivial behavioral predictions and to resolve open questions in the fields of single-particle tracking in living cells and movement ecology.
We report on a combined theoretical and experimental study of core-excitation spectra of gas and liquid phase methanol as obtained with the use of X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS). The electronic transitions are studied with computational methods that include strict and extended second-order algebraic diagrammatic construction [ADC(2) and ADC(2)-x], restricted active space second-order perturbation theory, and time-dependent density functional theory-providing a complete assignment of the near oxygen K-edge XAS. We show that multimode nuclear dynamics is of crucial importance for explaining the available experimental XAS and RIXS spectra. The multimode nuclear motion was considered in a recently developed "mixed representation" where dissociative states and highly excited vibrational modes are accurately treated with a time-dependent wave packet technique, while the remaining active vibrational modes are described using Franck-Condon amplitudes. Particular attention is paid to the polarization dependence of RIXS and the effects of the isotopic substitution on the RIXS profile in the case of dissociative core-excited states. Our approach predicts the splitting of the 2a RIXS peak to be due to an interplay between molecular and pseudo-atomic features arising in the course of transitions between dissociative core- and valence-excited states. The dynamical nature of the splitting of the 2a peak in RIXS of liquid methanol near pre-edge core excitation is shown. The theoretical results are in good agreement with our liquid phase measurements and gas phase experimental data available from the literature. (C) 2019 Author(s).
Local probes of the electronic ground state are essential for understanding hydrogen bonding in aqueous environments. When tuned to the dissociative core-excited state at the O1s pre-edge of water, resonant inelastic X-ray scattering back to the electronic ground state exhibits a long vibrational progression due to ultrafast nuclear dynamics. We show how the coherent evolution of the OH bonds around the core-excited oxygen provides access to high vibrational levels in liquid water. The OH bonds stretch into the long-range part of the potential energy curve, which makes the X-ray probe more sensitive than infra-red spectroscopy to the local environment. We exploit this property to effectively probe hydrogen bond strength via the distribution of intramolecular OH potentials derived from measurements. In contrast, the dynamical splitting in the spectral feature of the lowest valence-excited state arises from the short-range part of the OH potential curve and is rather insensitive to hydrogen bonding.
Tautomerism is one of the most important forms of isomerism, owing to the facile interconversion between species and the large differences in chemical properties introduced by the proton transfer connecting the tautomers. Spectroscopic techniques are often used for the characterization of tautomers. In this context, separating the overlapping spectral response of coexisting tautomers is a long-standing challenge in chemistry. Here, we demonstrate that by using resonant inelastic X-ray scattering tuned to the core excited states at the site of proton exchange between tautomers one is able to experimentally disentangle the manifold of valence excited states of each tautomer in a mixture. The technique is applied to the prototypical keto-enol equilibrium of 3-hydroxypyridine in aqueous solution. We detect transitions from the occupied orbitals into the LUMO for each tautomer in solution, which report on intrinsic and hydrogen-bond-induced orbital polarization within the pi and sigma manifolds at the proton-transfer site.
We show that, although the equilibrium band dispersion of the Shockley-type surface state of two-dimensional Au(111) quantum films grown on W(110) does not deviate from the expected free-electron-like behavior, its nonequilibrium energy-momentum dispersion probed by time- and angle-resolved photoemission exhibits a remarkable kink above the Fermi level due to a significant enhancement of the effective mass. The kink is pronounced for certain thicknesses of the Au quantum well and vanishes in the very thin limit. We identify the kink as induced by the coupling between the Au(111) surface state and emergent quantum-well states which probe directly the buried gold-tungsten interface. The signatures of the coupling are further revealed by our time-resolved measurements which show that surface state and quantum-well states thermalize together behaving as dynamically locked electron populations. In particular, relaxation of hot carriers following laser excitation is similar for both surface state and quantum-well states and much slower than expected for a bulk metallic system. The influence of quantum confinement on the interplay between elementary scattering processes of the electrons at the surface and ultrafast carrier transport in the direction perpendicular to the surface is shown to be the reason for the slow electron dynamics.
We show that, although the equilibrium band dispersion of the Shockley-type surface state of two-dimensional Au(111) quantum films grown on W(110) does not deviate from the expected free-electron-like behavior, its nonequilibrium energy-momentum dispersion probed by time- and angle-resolved photoemission exhibits a remarkable kink above the Fermi level due to a significant enhancement of the effective mass. The kink is pronounced for certain thicknesses of the Au quantum well and vanishes in the very thin limit. We identify the kink as induced by the coupling between the Au(111) surface state and emergent quantum-well states which probe directly the buried gold-tungsten interface. The signatures of the coupling are further revealed by our time-resolved measurements which show that surface state and quantum-well states thermalize together behaving as dynamically locked electron populations. In particular, relaxation of hot carriers following laser excitation is similar for both surface state and quantum-well states and much slower than expected for a bulk metallic system. The influence of quantum confinement on the interplay between elementary scattering processes of the electrons at the surface and ultrafast carrier transport in the direction perpendicular to the surface is shown to be the reason for the slow electron dynamics.
The Thesis is focused on the properties of self-organized nanostructures. Atomic and electronic properties of different systems have been investigated using methods of electron diffraction, scanning tunneling microscopy and photoelectron spectroscopy. Implementation of the STM technique (including design, construction, and tuning of the UHV experimental set-up) has been done in the framework of present work. This time-consuming work is reported to greater detail in the experimental part of this Thesis. The scientific part starts from the study of quantum-size effects in the electronic structure of a two-dimensional Ag film on the supporting substrate Ni(111). Distinct quantum well states in the sp-band of Ag were observed in photoelectron spectra. Analysis of thickness- and angle-dependent photoemission supplies novel information on the properties of the interface. For the first time the Ni(111) relative band gap was indirectly probed in the ground-state through the electronic structure of quantum well states in the adlayer. This is particularly important for Ni where valence electrons are strongly correlated. Comparison of the experiment with calculations performed in the formalism of the extended phase accumulation model gives the substrate gap which is fully consistent with the one obtained by ab-initio LDA calculations. It is, however, in controversy to the band structure of Ni measured directly by photoemission. These results lend credit to the simplest view of photoemission from Ni, assigning early observed contradictions between theory and experiments to electron correlation effects in the final state of photoemission. Further, nanosystems of lower dimensionality have been studied. Stepped surfaces W(331) and W(551) were used as one-dimensional model systems and as templates for self-organization of Au nanoclusters. Photon energy dependent photoemission revealed a surface resonance which was never observed before on W(110) which is the base plane of the terrace microsurfaces. The dispersion E(k) of this state measured on stepped W(331) and W(551) with angle-resolved photoelectron spectroscopy is modified by a strong umklapp effect. It appears as two parabolas shifted symmetrically relative to the microsurface normal by half of the Brillouin zone of the step superlattice. The reported results are very important for understanding of the electronic properties of low-dimensional nanostructures. It was also established that W(331) and W(551) can serve as templates for self-organization of metallic nanostructures. A combined study of electronic and atomic properties of sub-monolayer amounts of gold deposited on these templates have shown that if the substrate is slightly pre-oxidized and the temperature is elevated, then Au can alloy with the first monolayer of W. As a result, a nanostructure of uniform clusters of a surface alloy is produced all over the steps. Such clusters feature a novel sp-band in the vicinity of the Fermi level, which appears split into constant energy levels due to effects of lateral quantization. The last and main part of this work is devoted to large-scale reconstructions on surfaces and nanostructures self-assembled on top. The two-dimensional surface carbide W(110)/C-R(15x3) has been extensively investigated. Photoemission studies of quantum size effects in the electronic structure of this reconstruction, combined with an investigation of its surface geometry, lead to an advanced structural model of the carbide overlayer. It was discovered that W(110)/C-R(15x3) can control self-organization of adlayers into nanostructures with extremely different electronic and structural properties. Thus, it was established that at elevated temperature the R(15x3) superstructure controls the self-assembly of sub-monolayer amounts of Au into nm-wide nanostripes. Based on the results of core level photoemission, the R(15x3)-induced surface alloying which takes place between Au and W can be claimed as driving force of self-organization. The observed stripes exhibit a characteristic one-dimensional electronic structure with laterally quantized d-bands. Obviously, these are very important for applications, since dimensions of electronic devices have already stepped into the nm-range, where quantum-size phenomena must undoubtedly be considered. Moreover, formation of perfectly uniform molecular clusters of C60 was demonstrated and described in terms of the van der Waals formalism. It is the first experimental observation of two-dimensional fullerene nanoclusters with "magic numbers". Calculations of the cluster potentials using the static approach have revealed characteristic minima in the interaction energy. They are achieved for 4 and 7 molecules per cluster. The obtained "magic numbers" and the corresponding cluster structures are fully consistent with the results of the STM measurements.
The stable operation of a turbulent combustor is not completely silent; instead, there is a background of small amplitude aperiodic acoustic fluctuations known as combustion noise. Pressure fluctuations during this state of combustion noise are multifractal due to the presence of multiple temporal scales that contribute to its dynamics. However, existing models are unable to capture the multifractality in the pressure fluctuations. We conjecture an underlying fractional dynamics for the thermoacoustic system and obtain a fractional-order model for pressure fluctuations. The data from this model has remarkable visual similarity to the experimental data and also has a wide multifractal spectrum during the state of combustion noise. Quantitative similarity with the experimental data in terms of the Hurst exponent and the multifractal spectrum is observed during the state of combustion noise. This model is also able to produce pressure fluctuations that are qualitatively similar to the experimental data acquired during intermittency and thermoacoustic instability. Furthermore, we argue that the fractional dynamics vanish as we approach the state of thermoacoustic instability.
Absorption Tails of Donor
(2017)
In disordered organic semiconductors, the transfer of a rather localized charge carrier from one site to another triggers a deformation of the molecular structure quantified by the intramolecular relaxation energy. A similar structural relaxation occurs upon population of intermolecular charge-transfer (CT) states formed at organic electron donor (D)-acceptor (A) interfaces. Weak CT absorption bands for D A complexes occur at photon energies below the optical gaps of both the donors and the C-60 acceptor as a result of optical transitions from the neutral ground state to the ionic CT state. In this work, we show that temperature-activated intramolecular vibrations of the ground state play a major role in determining the line shape of such CT absorption bands. This allows us to extract values for the relaxation energy related to the geometry change from neutral to ionic CT complexes. Experimental values for the relaxation energies of 20 D:C-60 CT complexes correlate with values calculated within density functional theory. These results provide an experimental method for determining the polaron relaxation energy in solid-state organic D-A blends and show the importance of a reduced relaxation energy, which we introduce to characterize thermally activated CT processes.
Many solar wind observations at 1 au indicate that the proton (as well as electron) temperature anisotropy is limited. The data distribution in the (A(a), beta(a),(parallel to))-plane have a rhombic-shaped form around beta(a),(parallel to) similar to 1. The boundaries of the temperature anisotropy at beta(a),(parallel to) > 1 can be well explained by the threshold conditions of the mirror (whistler) and oblique proton (electron) firehose instabilities in a bi-Maxwellian plasma, whereas the physical mechanism of the similar restriction at beta(a),(parallel to) < 1 is still under debate. One possible option is Coulomb collisions, which we revisit in the current work. We derive the relaxation rate nu(A)(aa) of the temperature anisotropy in a bi-Maxwellian plasma that we then study analytically and by observed proton data from WIND. We found that nu(A)(pp) increases toward small beta(p),(parallel to) < 1. We matched the data distribution in the (A(p), beta(p),(parallel to))-plane with the constant contour nu(A)(pp) = 2.8 . 10(-6) s(-1), corresponding to the minimum value for collisions to play a role. This contour fits rather well the left boundary of the rhombic-shaped data distribution in the (A(p), beta(p),(parallel to))-plane. Thus, Coulomb collisions are an interesting candidate for explaining the limitations of the temperature anisotropy in the solar wind with small beta(a),(parallel to) < 1 at 1 au.
This work revisits the electrostatic instability for blazar-induced pair beams propagating through the intergalactic medium (IGM) using linear analysis and PIC simulations. We study the impact of the realistic distribution function of pairs resulting from the interaction of high-energy gamma-rays with the extragalactic background light. We present analytical and numerical calculations of the linear growth rate of the instability for the arbitrary orientation of wave vectors. Our results explicitly demonstrate that the finite angular spread of the beam dramatically affects the growth rate of the waves, leading to the fastest growth for wave vectors quasi-parallel to the beam direction and a growth rate at oblique directions that is only a factor of 2-4 smaller compared to the maximum. To study the nonlinear beam relaxation, we performed PIC simulations that take into account a realistic wide-energy distribution of beam particles. The parameters of the simulated beam-plasma system provide an adequate physical picture that can be extrapolated to realistic blazar-induced pairs. In our simulations, the beam looses only 1% of its energy, and we analytically estimate that the beam would lose its total energy over about 100 simulation times. An analytical scaling is then used to extrapolate the parameters of realistic blazar-induced pair beams. We find that they can dissipate their energy slightly faster by the electrostatic instability than through inverse-Compton scattering. The uncertainties arising from, e.g., details of the primary gamma-ray spectrum are too large to make firm statements for individual blazars, and an analysis based on their specific properties is required.
Due to the relatively low rate of particle-particle collisions in the solar wind, kinetic instabilities (e.g., the mirror and firehose) play an important role in regulating large deviations from temperature isotropy. These instabilities operate in the high beta(parallel to) > 1 plasmas, and cannot explain the other limits of the temperature anisotropy reported by observations in the low beta beta(parallel to) < 1 regimes. However, the instability conditions are drastically modified in the presence of streaming (or counterstreaming) components, which are ubiquitous in space plasmas. These effects have been analyzed for the solar wind conditions in a large interval of heliospheric distances, 0.3-2.5 AU. It was found that proton counter-streams are much more crucial for plasma stability than electron ones. Moreover, new instability thresholds can potentially explain all observed bounds on the temperature anisotropy, and also the level of differential streaming in the solar wind.
A polymer is a large molecule made up of many elementary chemical units, joined together by covalent bonds (for example, polyethylene). Polyelectrolytes (PELs) are polymer chains containing a certain amount of ionizable monomers. With their specific properties PELs acquire big importance in molecular and cell biology as well as in technology. Compared to neutral polymers the theory of PELs is less understood. In particular, this is valid for PELs in poor solvents. A poor solvent environment causes an effective attraction between monomers. Hence, for PELs in a poor solvent, there occurs a competition between attraction and repulsion. Strong or quenched PELs are completely dissociated at any accessible pH. The position of charges along the chain is fixed by chemical synthesis. On the other hand, in weak or annealed PELs dissociation of charges depends on solution pH. For the first time the simulation results have given direct evidence that at rather poor solvents an annealed PEL indeed undergoes a first-order phase transition when the chemical potential (solution pH) reaches at a certain value. The discontinuous transition occurs between a weakly charged compact globular structure and a strongly charged stretched configuration. At not too poor solvents theory predicts that globule would become unstable with respect to the formation of pearl-necklaces. The results show that pearl-necklaces exist in annealed PELs indeed. Furthermore, as predicted by theory, the simulation results have shown that annealed PELs display a sharp transition from a highly charged stretched state to a weakly charged globule at a critical salt concentration.
Hot localised charge carriers on the Si(111)-7×7 surface are modelled by small charged clusters. Such resonances induce non-local desorption, i.e. more than 10 nm away from the injection site, of chlorobenzene in scanning tunnelling microscope experiments. We used such a cluster model to characterise resonance localisation and vibrational activation for positive and negative resonances recently. In this work, we investigate to which extent the model depends on details of the used cluster or quantum chemistry methods and try to identify the smallest possible cluster suitable for a description of the neutral surface and the ion resonances. Furthermore, a detailed analysis for different chemisorption orientations is performed. While some properties, as estimates of the resonance energy or absolute values for atomic changes, show such a dependency, the main findings are very robust with respect to changes in the model and/or the chemisorption geometry.
Quantum steering describes how local actions on a quantum system can affect another, spacelike separated, quantum state. Lately, quantum steering has been formulated also for timelike scenarios and for quantum channels. We approach all the three scenarios as one using tools from Stinespring dilations of quantum channels. By applying our technique we link all three steering problems one-to-one with the incompatibility of quantum measurements, a result formerly known only for spatial steering. We exploit this connection by showing how measurement uncertainty relations can be used as tight steering inequalities for all three scenarios. Moreover, we show that certain notions of temporal and spatial steering are fully equivalent and prove a hierarchy between temporal steering and macrorealistic hidden variable models.
We report on the adsorption kinetics of azoben-zene-containing surfactants on solid surfaces of different hydrophobicity. The understanding of this processes is of great importance for many interfacial phenomena that can be actuated and triggered by light, since the surfactant molecules contain a photoresponsive azobenzene group in their hydrophobic tail. Three surfactant types are studied, differing in the spacer connecting the headgroup and the azobenzene unit by between 6 and 10 CH2 groups. Under irradiation with light of a suitable wavelength, the azobenzene undergoes reversible photoisomerization between two states, a nonpolar trans-state and a highly polar cis-state. Consequently, the surfactant molecule changes its hydrophobicity and thus affinity to a surface depending on the photoisomerization state of the azobenzene. The adsorption behavior on hydrophilic (glass) and hydrophobic (TeflonAF) surfaces is analyzed using quartz crystal microbalance with dissipation (QCM-D) and zeta-potential measurements. At equilibrium, the adsorbed surfactant amount is almost twice as large on glass compared to TeflonAF for both isomers. The adsorption rate for the trans-isomers on both surfaces is similar, but the desorption rate of the trans-isomers is faster at the glass-water interface than at the Teflon-water interface. This result demonstrates that the trans-isomers have higher affinity for the glass surface, so the trans-to-cis ratios on glass and TeflonAF are 80/1 and 2/1, respectively, with similar trends for all three surfactant types.
My thesis is concerned with several new noise-induced phenomena in excitable neural models, especially those with FitzHugh-Nagumo dynamics. In these effects the fluctuations intrinsically present in any complex neural network play a constructive role and improve functionality. I report the occurrence of Vibrational Resonance in excitable systems. Both in an excitable electronic circuit and in the FitzHugh-Nagumo model, I show that an optimal amplitude of high-frequency driving enhances the response of an excitable system to a low-frequency signal. Additionally, the influence of additive noise and the interplay between Stochastic and Vibrational Resonance is analyzed. Further, I study systems which combine both oscillatory and excitable properties, and hence intrinsically possess two internal frequencies. I show that in such a system the effect of Stochastic Resonance can be amplified by an additional high-frequency signal which is in resonance with the oscillatory frequency. This amplification needs much lower noise intensities than for conventional Stochastic Resonance in excitable systems. I study frequency selectivity in noise-induced subthreshold signal processing in a system with many noise-supported stochastic attractors. I show that the response of the coupled elements at different noise levels can be significantly enhanced or reduced by forcing some elements into resonance with these new frequencies which correspond to appropriate phase-relations. A noise-induced phase transition to excitability is reported in oscillatory media with FitzHugh-Nagumo dynamics. This transition takes place via noise-induced stabilization of a deterministically unstable fixed point of the local dynamics, while the overall phase-space structure of the system is maintained. The joint action of coupling and noise leads to a different type of phase transition and results in a stabilization of the system. The resulting noise-induced regime is shown to display properties characteristic of excitable media, such as Stochastic Resonance and wave propagation. This effect thus allows the transmission of signals through an otherwise globally oscillating medium. In particular, these theoretical findings suggest a possible mechanism for suppressing undesirable global oscillations in neural networks (which are usually characteristic of abnormal medical conditions such as Parkinson′s disease or epilepsy), using the action of noise to restore excitability, which is the normal state of neuronal ensembles.
Intermolecular charge-transfer states at the interface between electron donating (D) and accepting (A) materials are crucial for the operation of organic solar cells but can also be exploited for organic light-emitting diodes(1,2). Non-radiative charge-transfer state decay is dominant in state-of-the-art D-A-based organic solar cells and is responsible for large voltage losses and relatively low power-conversion efficiencies as well as electroluminescence external quantum yields in the 0.01-0.0001% range(3,4). In contrast, the electroluminescence external quantum yield reaches up to 16% in D-A-based organic light-emitting diodes(5-7). Here, we show that proper control of charge-transfer state properties allows simultaneous occurrence of a high photovoltaic and emission quantum yield within a single, visible-light-emitting D-A system. This leads to ultralow-emission turn-on voltages as well as significantly reduced voltage losses upon solar illumination. These results unify the description of the electro-optical properties of charge-transfer states in organic optoelectronic devices and foster the use of organic D-A blends in energy conversion applications involving visible and ultraviolet photons(8-11).
Additive manufacturing (AM) of metals and in particular laser powder bed fusion (LPBF) enables a degree of freedom in design unparalleled by conventional subtractive methods. To ensure that the designed precision is matched by the produced LPBF parts, a full understanding of the interaction between the laser and the feedstock powder is needed. It has been shown that the laser also melts subjacent layers of material underneath. This effect plays a key role when designing small cavities or overhanging structures, because, in these cases, the material underneath is feed-stock powder. In this study, we quantify the extension of the melt pool during laser illumination of powder layers and the defect spatial distribution in a cylindrical specimen. During the LPBF process, several layers were intentionally not exposed to the laser beam at various locations, while the build process was monitored by thermography and optical tomography. The cylinder was finally scanned by X-ray computed tomography (XCT). To correlate the positions of the unmolten layers in the part, a staircase was manufactured around the cylinder for easier registration. The results show that healing among layers occurs if a scan strategy is applied, where the orientation of the hatches is changed for each subsequent layer. They also show that small pores and surface roughness of solidified material below a thick layer of unmolten material (>200 mu m) serve as seeding points for larger voids. The orientation of the first two layers fully exposed after a thick layer of unmolten powder shapes the orientation of these voids, created by a lack of fusion.
Both climate-change damages and climate-change mitigation will incur economic costs. While the risk of severe damages increases with the level of global warming (Dell et al., 2014; IPCC, 2014b, 2018; Lenton et al., 2008), mitigating costs increase steeply with more stringent warming limits (IPCC, 2014a; Luderer et al., 2013; Rogelj et al., 2015). Here, we show that the global warming limit that minimizes this century's total economic costs of climate change lies between 1.9 and 2 ∘C, if temperature changes continue to impact national economic growth rates as observed in the past and if instantaneous growth effects are neither compensated nor amplified by additional growth effects in the following years. The result is robust across a wide range of normative assumptions on the valuation of future welfare and inequality aversion. We combine estimates of climate-change impacts on economic growth for 186 countries (applying an empirical damage function from Burke et al., 2015) with mitigation costs derived from a state-of-the-art energy–economy–climate model with a wide range of highly resolved mitigation options (Kriegler et al., 2017; Luderer et al., 2013, 2015). Our purely economic assessment, even though it omits non-market damages, provides support for the international Paris Agreement on climate change. The political goal of limiting global warming to “well below 2 degrees” is thus also an economically optimal goal given above assumptions on adaptation and damage persistence.
What comes NeXT?
(2019)
Here, we report on a new record in the acquisition time for fast neutron tomography. With an optimized imaging setup, it was possible to acquire single radiographic projection images with 10 ms and full tomographies with 155 projections images and a physical spatial resolution of 200 mu m within 1.5 s. This is about 6.7 times faster than the current record. We used the technique to investigate the water infiltration in the soil with a living lupine root system. The fast imaging setup will be part of the future NeXT instrument at ILL in Grenoble with a great field of possible future applications. (C) 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
Low-dimensional description for ensembles of identical phase oscillators subject to Cauchy noise
(2020)
We study ensembles of globally coupled or forced identical phase oscillators subject to independent white Cauchy noise. We demonstrate that if the oscillators are forced in several harmonics, stationary synchronous regimes can be exactly described with a finite number of complex order parameters. The corresponding distribution of phases is a product of wrapped Cauchy distributions. For sinusoidal forcing, the Ott-Antonsen low-dimensional reduction is recovered.
We formulate a linear phase and frequency response theory for hyperbolic flows, which generalizes phase response theory for autonomous limit cycle oscillators to hyperbolic chaotic dynamics. The theory is based on a shadowing conjecture, stating the existence of a perturbed trajectory shadowing every unperturbed trajectory on the system attractor for any small enough perturbation of arbitrary duration and a corresponding unique time isomorphism, which we identify as phase such that phase shifts between the unperturbed trajectory and its perturbed shadow are well defined. The phase sensitivity function is the solution of an adjoint linear equation and can be used to estimate the average change of phase velocity to small time dependent or independent perturbations. These changes in frequency are experimentally accessible, giving a convenient way to define and measure phase response curves for chaotic oscillators. The shadowing trajectory and the phase can be constructed explicitly in the tangent space of an unperturbed trajectory using co-variant Lyapunov vectors. It can also be used to identify the limits of the regime of linear response.
This work is concerned with the spatio-temporal structures that emerge when non-identical, diffusively coupled oscillators synchronize. It contains analytical results and their confirmation through extensive computer simulations. We use the Kuramoto model which reduces general oscillatory systems to phase dynamics. The symmetry of the coupling plays an important role for the formation of patterns. We have studied the ordering influence of an asymmetry (non-isochronicity) in the phase coupling function on the phase profile in synchronization and the intricate interplay between this asymmetry and the frequency heterogeneity in the system. The thesis is divided into three main parts. Chapter 2 and 3 introduce the basic model of Kuramoto and conditions for stable synchronization. In Chapter 4 we characterize the phase profiles in synchronization for various special cases and in an exponential approximation of the phase coupling function, which allows for an analytical treatment. Finally, in the third part (Chapter 5) we study the influence of non-isochronicity on the synchronization frequency in continuous, reaction diffusion systems and discrete networks of oscillators.
We have developed a method for deriving systems of closed equations for the dynamics of order parameters in the ensembles of phase oscillators. The Ott-Antonsen equation for the complex order parameter is a particular case of such equations. The simplest nontrivial extension of the Ott-Antonsen equation corresponds to two-bunch states of the ensemble. Based on the equations obtained, we study the dynamics of multi-bunch chimera states in coupled Kuramoto-Sakaguchi ensembles. We show an increase in the dimensionality of the system dynamics for two-bunch chimeras in the case of identical phase elements and a transition to one-bunch "Abrams chimeras" for imperfect identity (in the latter case, the one-bunch chimeras become attractive).
In the present dissertation paper an approach which ensures an efficient control of such diverse systems as noisy or chaotic oscillators and neural ensembles is developed. This approach is implemented by a simple linear feedback loop. The dissertation paper consists of two main parts. One part of the work is dedicated to the application of the suggested technique to a population of neurons with a goal to suppress their synchronous collective dynamics. The other part is aimed at investigating linear feedback control of coherence of a noisy or chaotic self-sustained oscillator. First we start with a problem of suppressing synchronization in a large population of interacting neurons. The importance of this task is based on the hypothesis that emergence of pathological brain activity in the case of Parkinson's disease and other neurological disorders is caused by synchrony of many thousands of neurons. The established therapy for the patients with such disorders is a permanent high-frequency electrical stimulation via the depth microelectrodes, called Deep Brain Stimulation (DBS). In spite of efficiency of such stimulation, it has several side effects and mechanisms underlying DBS remain unclear. In the present work an efficient and simple control technique is suggested. It is designed to ensure suppression of synchrony in a neural ensemble by a minimized stimulation that vanishes as soon as the tremor is suppressed. This vanishing-stimulation technique would be a useful tool of experimental neuroscience; on the other hand, control of collective dynamics in a large population of units represents an interesting physical problem. The main idea of suggested approach is related to the classical problem of oscillation theory, namely the interaction between a self-sustained (active) oscillator and a passive load (resonator). It is known that under certain conditions the passive oscillator can suppress the oscillations of an active one. In this thesis a much more complicated case of active medium, which itself consists of thousands of oscillators is considered. Coupling this medium to a specially designed passive oscillator, one can control the collective motion of the ensemble, specifically can enhance or suppress it. Having in mind a possible application in neuroscience, the problem of suppression is concentrated upon. Second, the efficiency of suggested suppression scheme is illustrated by considering more complex case, i.e. when the population of neurons generating the undesired rhythm consists of two non-overlapping subpopulations: the first one is affected by the stimulation, while the collective activity is registered from the second one. Generally speaking, the second population can be by itself both active and passive; both cases are considered here. The possible applications of suggested technique are discussed. Third, the influence of the external linear feedback on coherence of a noisy or chaotic self-sustained oscillator is considered. Coherence is one of the main properties of self-oscillating systems and plays a key role in the construction of clocks, electronic generators, lasers, etc. The coherence of a noisy limit cycle oscillator in the context of phase dynamics is evaluated by the phase diffusion constant, which is in its turn proportional to the width of the spectral peak of oscillations. Many chaotic oscillators can be described within the framework of phase dynamics, and, therefore, their coherence can be also quantified by the way of the phase diffusion constant. The analytical theory for a general linear feedback, considering noisy systems in the linear and Gaussian approximation is developed and validated by numerical results.
Context. On 27 April 2015, when comet 67P/Churyumov-Gerasimenko was at 1.76 au from the Sun and moving toward perihelion, the OSIRIS and VIRTIS-M instruments on board the Rosetta spacecraft simultaneously observed the evolving dust and gas coma during a complete rotation of the comet. Aims. We aim to characterize the spatial distribution of dust, H2O, and CO2 gas in the inner coma. To do this, we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the rate we calculated using a thermophysical model. Methods. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O emission band), and 4200 nm (CO2 emission band). We measured the average signal in a circular annulus to study the spatial variation around the comet, and in a sector of the annulus to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet center. Results. The spatial correlation between dust and water, both coming from the sunlit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0 degrees subsolar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production with respect to the value calculated using a simple thermophysical model is observed when the head lobe and regions of the southern hemisphere with strong seasonal variations are illuminated (subsolar longitude 270 degrees-50 degrees). A drastic decrease in dust production when the water production (both measured and from the model) displays a maximum occurs when typical northern consolidated regions are illuminated and the southern hemisphere regions with strong seasonal variations are instead in shadow (subsolar longitude 50 degrees-90 degrees). Possible explanations of these observations are presented and discussed.
The European potential for renewable electricity is sufficient to enable fully renewable supply on different scales, from self-sufficient, subnational regions to an interconnected continent. We not only show that a continental-scale system is the cheapest, but also that systems on the national scale and below are possible at cost penalties of 20% or less. Transmission is key to low cost, but it is not necessary to vastly expand the transmission system. When electricity is transmitted only to balance fluctuations, the transmission grid size is comparable to today's, albeit with expanded cross-border capacities. The largest differences across scales concern land use and thus social acceptance: in the continental system, generation capacity is concentrated on the European periphery, where the best resources are. Regional systems, in contrast, have more dispersed generation. The key trade-off is therefore not between geographic scale and cost, but between scale and the spatial distribution of required generation and transmission infrastructure.
The European potential for renewable electricity is sufficient to enable fully renewable supply on different scales, from self-sufficient, subnational regions to an interconnected continent. We not only show that a continental-scale system is the cheapest, but also that systems on the national scale and below are possible at cost penalties of 20% or less. Transmission is key to low cost, but it is not necessary to vastly expand the transmission system. When electricity is transmitted only to balance fluctuations, the transmission grid size is comparable to today's, albeit with expanded cross-border capacities. The largest differences across scales concern land use and thus social acceptance: in the continental system, generation capacity is concentrated on the European periphery, where the best resources are. Regional systems, in contrast, have more dispersed generation. The key trade-off is therefore not between geographic scale and cost, but between scale and the spatial distribution of required generation and transmission infrastructure.
In processing and data storage mainly ferromagnetic (FM) materials are being used. Approaching physical limits, new concepts have to be found for faster, smaller switches, for higher data densities and more energy efficiency. Some of the discussed new concepts involve the material classes of correlated oxides and materials with antiferromagnetic coupling. Their applicability depends critically on their switching behavior, i.e., how fast and how energy efficient material properties can be manipulated. This thesis presents investigations of ultrafast non-equilibrium phase transitions on such new materials. In transition metal oxides (TMOs) the coupling of different degrees of freedom and resulting low energy excitation spectrum often result in spectacular changes of macroscopic properties (colossal magneto resistance, superconductivity, metal-to-insulator transitions) often accompanied by nanoscale order of spins, charges, orbital occupation and by lattice distortions, which make these material attractive. Magnetite served as a prototype for functional TMOs showing a metal-to-insulator-transition (MIT) at T = 123 K. By probing the charge and orbital order as well as the structure after an optical excitation we found that the electronic order and the structural distortion, characteristics of the insulating phase in thermal equilibrium, are destroyed within the experimental resolution of 300 fs. The MIT itself occurs on a 1.5 ps timescale. It shows that MITs in functional materials are several thousand times faster than switching processes in semiconductors. Recently ferrimagnetic and antiferromagnetic (AFM) materials have become interesting. It was shown in ferrimagnetic GdFeCo, that the transfer of angular momentum between two opposed FM subsystems with different time constants leads to a switching of the magnetization after laser pulse excitation. In addition it was theoretically predicted that demagnetization dynamics in AFM should occur faster than in FM materials as no net angular momentum has to be transferred out of the spin system. We investigated two different AFM materials in order to learn more about their ultrafast dynamics. In Ho, a metallic AFM below T ≈ 130 K, we found that the AFM Ho can not only be faster but also ten times more energy efficiently destroyed as order in FM comparable metals. In EuTe, an AFM semiconductor below T ≈ 10 K, we compared the loss of magnetization and laser-induced structural distortion in one and the same experiment. Our experiment shows that they are effectively disentangled. An exception is an ultrafast release of lattice dynamics, which we assign to the release of magnetostriction. The results presented here were obtained with time-resolved resonant soft x-ray diffraction at the Femtoslicing source of the Helmholtz-Zentrum Berlin and at the free-electron laser in Stanford (LCLS). In addition the development and setup of a new UHV-diffractometer for these experiments will be reported.