Refine
Year of publication
- 2019 (308) (remove)
Document Type
- Article (254)
- Doctoral Thesis (19)
- Other (19)
- Postprint (12)
- Review (3)
- Monograph/Edited Volume (1)
Language
- English (308) (remove)
Is part of the Bibliography
- yes (308)
Keywords
- diffusion (9)
- stars: evolution (7)
- stars: winds, outflows (7)
- methods: numerical (6)
- quasars: absorption lines (6)
- stars: massive (6)
- anomalous diffusion (5)
- intergalactic medium (5)
- stars: mass-loss (5)
- subdwarfs (5)
Institute
- Institut für Physik und Astronomie (308) (remove)
Extremely low-mass white dwarf stars (ELMs) are M < 0.3 M-circle dot helium-core white dwarfs born either as a result of a common-envelope phase or after a stable Roche lobe overflow episode in a multiple system. The Universe is not old enough for ELMs to have formed through single-star evolution channels. As remnants of binary evolution, ELMs can shed light onto the poorly understood phase of common-envelope evolution and provide constraints to the physics of mass accretion. Most known ELMs will merge in less than a Hubble time, providing an important contribution to the signal to be detected by upcoming space-based gravitational wave detectors. There are currently less than 150 known ELMs; most were selected by colour, focusing on hot objects, in a magnitude-limited survey of the Northern hemisphere only. Recent theoretical models have predicted a much larger space density for ELMs than estimated observationally based on this limited sample. In order to perform meaningful comparisons with theoretical models and test their predictions, a larger well-defined sample is required. In this work, we present a catalogue of ELM candidates selected from the second data release of Gaia (DR2). We have used predictions from theoretical models and analysed the properties of the known sample to map the space spanned by ELMs in the Gaia Hertzsprung-Russell diagram. Defining a set of colour cuts and quality flags, we have obtained a final sample of 5762 ELM candidates down to T-eff approximate to 5000 K.
It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the Λ cold dark matter (ΛCDM) model. One type of anisotropy recently detected in the Sloan Digital Sky Survey (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called lopsided. Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set-up when comparing pairs of galaxies in ΛCDM simulations indicating that such a set-up is not inconsistent with ΛCDM, the origin has yet to be explained. Here we examine the origin of such lopsided set-ups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at z = 0 display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards z = 0. Secondly, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the z = 0 signal is driven primarily by satellites that are on first approach, who have yet to experience a ‘flyby’. This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.
Galaxies are surrounded by massive gas reservoirs ( i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated MilkyWay-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L-star galaxies is extremely diverse: column densities of commonly observed species span similar to 3-4 dex and their covering fractions range from similar to 5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions ( CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L-star galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.
We numerically and analytically analyze transitions between different synchronous states in a network of globally coupled phase oscillators with attractive and repulsive interactions. The elements within the attractive or repulsive group are identical, but natural frequencies of the groups differ. In addition to a synchronous two-cluster state, the system exhibits a solitary state, when a single oscillator leaves the cluster of repulsive elements, as well as partially synchronous quasiperiodic dynamics. We demonstrate how the transitions between these states occur when the repulsion starts to prevail over attraction.
Aims. In this study, we analyzed a filament system, which expanded between moving magnetic features (MMFs) of a decaying sunspot and opposite flux outside of the active region from the nearby quiet-Sun network. This configuration deviated from a classical arch filament system (AFS), which typically connects two pores in an emerging flux region. Thus, we called this system an extended AFS. We contrasted classical and extended AFSs with an emphasis on the complex magnetic structure of the latter. Furthermore, we examined the physical properties of the extended AFS and described its dynamics and connectivity. Methods. The extended AFS was observed with two instruments at the Dunn Solar Telescope (DST). The Rapid Oscillations in the Solar Atmosphere (ROSA) imager provided images in three different wavelength regions, which covered the dynamics of the extended AFS at different atmospheric heights. The Interferometric Bidimensional Spectropolarimeter (IBIS) provided spectroscopic Ha data and spectropolarimetric data that was obtained in the near-infrared (NIR) Call lambda 8542 angstrom line. We derived the corresponding line-of-sight (LOS) velocities and used He II lambda 304 angstrom extreme ultraviolet (EUV) images of the Atmospheric Imaging Assembly (AIA) and LOS magnetograms of the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) as context data. Results. The NIR Call Stokes-V maps are not suitable to definitively define a clear polarity inversion line and to classify this chromospheric structure. Nevertheless, this unusual AFS connects the MMFs of a decaying sunspot with the network field. At the southern footpoint, we measured that the flux decreases over time. We find strong downflow velocities at the footpoints of the extended AFS, which increase in a time period of 30 min. The velocities are asymmetric at both footpoints with higher velocities at the southern footpoint. An EUV brigthening appears in one of the arch filaments, which migrates from the northern footpoint toward the southern one. This activation likely influences the increasing redshift at the southern footpoint. Conclusions. The extended AFS exhibits a similar morphology as classical AFSs, for example, threaded filaments of comparable length and width. Major differences concern the connection from MMFs around the sunspot with the flux of the neighboring quietSun network, converging footpoint motions, and longer lifetimes of individual arch filaments of about one hour, while the extended AFS is still very dynamic.
The in-depth understanding of charge carrier photogeneration and recombination mechanisms in organic solar cells is still an ongoing effort. In donor:acceptor (bulk) heterojunction organic solar cells, charge photogeneration and recombination are inter-related via the kinetics of charge transfer states-being singlet or triplet states. Although high-charge-photogeneration quantum yields are achieved in many donor:acceptor systems, only very few systems show significantly reduced bimolecular recombination relative to the rate of free carrier encounters, in low-mobility systems. This is a serious limitation for the industrialization of organic solar cells, in particular when aiming at thick active layers. Herein, a meta-analysis of the device performance of numerous bulk heterojunction organic solar cells is presented for which field-dependent photogeneration, charge carrier mobility, and fill factor are determined. Herein, a "spin-related factor" that is dependent on the ratio of back electron transfer of the triplet charge transfer (CT) states to the decay rate of the singlet CT states is introduced. It is shown that this factor links the recombination reduction factor to charge-generation efficiency. As a consequence, it is only in the systems with very efficient charge generation and very fast CT dissociation that free carrier recombination is strongly suppressed, regardless of the spin-related factor.
Supercritical Kuramoto oscillators with distributed frequencies can be separated into two disjoint groups: an ordered one locked to the mean field, and a disordered one consisting of effectively decoupled oscillators-at least so in the thermodynamic limit. In finite ensembles, in contrast, such clear separation fails: The mean field fluctuates due to finite-size effects and thereby induces order in the disordered group. This publication demonstrates this effect, similar to noise-induced synchronization, in a purely deterministic system. We start by modeling the situation as a stationary mean field with additional white noise acting on a pair of unlocked Kuramoto oscillators. An analytical expression shows that the cross-correlation between the two increases with decreasing ratio of natural frequency difference and noise intensity. In a deterministic finite Kuramoto model, the strength of the mean-field fluctuations is inextricably linked to the typical natural frequency difference. Therefore, we let a fluctuating mean field, generated by a finite ensemble of active oscillators, act on pairs of passive oscillators with a microscopic natural frequency difference between which we then measure the cross-correlation, at both super- and subcritical coupling.
The purity of the analysed samples (e.g. quartz) with respect to chemical composition and radionuclide contamination is essential for geomorphologic applications using so-called terrestrial cosmogenic nuclides (TCNs). To guarantee this, numerous cleaning and dissolution procedures have been developed. At the DREsden Accelerator Mass Spectrometry (DREAMS) facility, we also work on enhancing the chemical quartz-enrichment methodology from bulk rock and dissolution of quartz. Repeated exposure of the bulk material to acid mixtures (HCl/H2SiF6) at room temperature for cleaning and its monitoring by optical microscopy works for most quartz-rich samples. The quartz dissolution in HF under rather mild conditions (at room temperature on a shaker-table) has the advantage to leave difficult-to-dissolve minerals (e.g., tourmaline, zircon, rutile, sillimanite, kyanite, chromite, corundum), not separated by other physical methods before, as residue. Our comparison with a high-temperature dissolution method (in a microwave) indicates an additional amount of interfering elements, such as in average about 3 mg of Ti, more than 7 mg of Al, and about 22 mu g of Be (for 50 g SiO2), is added to the sample, hence showing the superiority of our mild method. This way, we reduce problems for chemistry and AMS, but also ensure better comparability to production rates of cleaner stoichiometric quartz from calibration sites.
Long-term river evolution depends partly on crustal deformation, which shapes the topography crossed by rivers. On glacial timescales, ice-sheet growth and decay can produce crustal vertical motion of ∼10 mm/yr resulting from the solid Earth's adjustment to variations in ice and water loads, comparable to tectonically-driven rates in the most rapidly uplifting mountains on Earth. This process of glacial isostatic adjustment (GIA) can influence river courses and drainage basins substantially, particularly near former ice margins. We explore the extent to which GIA influenced the evolution of rivers along the United States east coast during the last glacial cycle. We compute gravitationally self-consistent GIA responses that incorporate recent constraints on the Laurentide Ice Sheet history through the last glacial build-up phase, and we connect the predicted variations in topography to abrupt changes in river dynamics recorded in the Hudson, Delaware, Susquehanna, and Potomac Rivers from 40 ka to present. To the extent that increases in sediment transport capacity imply increases in river incision rate, the GIA-driven changes in slope and drainage area are consistent with episodes of erosion and sedimentation observed in the Hudson, Delaware, and Potomac Rivers, but inconsistent with the observed accelerated river incision in the Susquehanna River at 30-14 ka. These analyses add to a growing body of evidence showing that GIA strongly influences river evolution over millennial timescales.
Simulations of the inner magnetospheric energetic electrons using the IMPTAM-VERB coupled model
(2019)
In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10-100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within similar to 0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.
Holocene glacial chronostratigraphies in glaciated valleys spread throughout the Himalayan-Tibetan orogen, including the Himalaya, Tibet, Pamir, and Tian Shan, are developed using a landsystems approach, detailed geomorphic mapping, and new and published Be-10 surface exposure dating. New studies in the Kulti valley of Lahul and the Parkachik valley of the Nun Kun massif of the Himalaya of northern India define three glacier advances at similar to 14.7, 12.2, 0.5 ka, in addition to one historically dated late 19th Century advance in the Kulti valley, and one Late Holocene advance at similar to 0.2 ka in the Parkachik valley. Three major climatic groups (subdivided into five climatic zones) are defined across the orogen using Cluster Analysis (CA) and Principal Component Analysis (PCA) to identify glaciated regions with comparable climatic characteristics to evaluate the timing, and extent of Holocene glacier advances across these regions. Our regional analyses across the Himalayan-Tibetan orogen suggest at least one Lateglacial (similar to 15.3-11.8 ka) and five Himalayan-Tibetan Holocene glacial stages (HTHS) at similar to 11.5-9.5, similar to 8.8-7.7, similar to 7.0-3.2, similar to 2.3-1.0, and <1 ka. The extent (amplitude) of glacier advances in 77 glaciated valleys is reconstructed and defined using equilibrium-line altitudes (ELAs). Modern glacier hypsometries are also assessed to help explain the intra-regional variations in glacier amplitudes during each regional glacier advance. A linear inverse glacier flow model is used to decipher the net changes in temperature (Delta T) between periods of reconstructed regional glacier advances in 66 glaciated valleys across different climatic regions throughout the orogen. The Be-10, ELAs, and Delta T data suggest enhanced monsoonal and increased precipitation during the Early Holocene, followed by relative cooling and increased aridity during the Mid- and Late Holocene that influenced glaciation. The sublimation-dominated cold-based glaciers in the northern regions of Himalayan-Tibetan orogen are more affected during these shifts in climate than the temperate glaciers in the south. (C) 2019 Elsevier Ltd. All rights reserved.
The intermediate He-enriched hot subdwarf star Feige 46 was recently reported as the second member of the V366 Aqr (or He-sdOBV) pulsating class. Feige 46 is very similar to the prototype of the class, LS IV-14 degrees 116, not only in terms of pulsational properties, but also in terms of atmospheric parameters and kinematic properties. LS IV-14 degrees 116 is additionally characterized by a very peculiar chemical composition, with extreme overabundances of the trans-iron elements Ge, Sr, Y, and Zr. We investigate the possibility that the similarity between the two pulsators extends to their chemical composition. We retrieved archived optical and UV spectroscopic observations of Feige 46 and performed an abundance analysis using model atmospheres and synthetic spectra computed with TLUSTY and SYNSPEC. In total, we derived abundances for 16 elements and provide upper limits for four additional elements. Using absorption lines in the optical spectrum of the star we measure an enrichment of more than 10 000x solar for yttrium and zirconium. The UV spectrum revealed that strontium is equally enriched. Our results confirm that Feige 46 is not only a member of the now growing group of heavy metal subdwarfs, but also has an abundance pattern that is remarkably similar to that of LS IV-14 degrees 116.
Transport in exclusion processes with one-step memory: density dependence and optimal acceleration
(2019)
We study a lattice gas of persistent walkers, in which each site is occupied by at most one particle and the direction each particle attempts to move to depends on its last step. We analyse the mean squared displacement (MSD) of the particles as a function of the particle density and their persistence (the tendency to continue moving in the same direction). For positive persistence the MSD behaves as expected: it increases with the persistence and decreases with the density. However, for strong anti-persistence we find two different regimes, in which the dependence of the MSD on the density is non-monotonic. For very strong anti-persistence there is an optimal density at which the MSD reaches a maximum. In an intermediate regime, the MSD as a function of the density exhibits both a minimum and a maximum, a phenomenon which has not been observed before. We derive a mean-field theory which qualitatively explains this behaviour.
EVR-CB-001: An Evolving, Progenitor, White Dwarf Compact Binary Discovered with the Evryscope
(2019)
We present EVR-CB-001, the discovery of a compact binary with an extremely low-mass (0.21 +/- 0.05M(circle dot)) helium core white dwarf progenitor (pre-He WD) and an unseen low-mass (0.32 +/- 0.06M(circle dot)) helium white dwarf (He WD) companion. He WDs are thought to evolve from the remnant helium-rich core of a main-sequence star stripped during the giant phase by a close companion. Low-mass He WDs are exotic objects (only about 0.2% of WDs are thought to be less than 0.3 M-circle dot), and are expected to be found in compact binaries. Pre-He WDs are even rarer, and occupy the intermediate phase after the core is stripped, but before the star becomes a fully degenerate WD and with a larger radius (approximate to 0.2R(circle dot)) than a typical WD. The primary component of EVR-CB-001 (the pre-He WD) was originally thought to be a hot subdwarf (sdB) star from its blue color and under-luminous magnitude, characteristic of sdBs. The mass, temperature (T-eff = 18,500 +/- 500 K), and surface gravity (log(g) = 4.96 +/- 0.04) solutions from this work are lower than values for typical hot subdwarfs. The primary is likely to be a post-red-giant branch, pre-He WD contracting into a He WD, and at a stage that places it nearest to sdBs on color-magnitude and T-eff-log(g) diagrams. EVR-CB-001 is expected to evolve into a fully double degenerate, compact system that should spin down and potentially evolve into a single hot subdwarf star. Single hot subdwarfs are observed, but progenitor systems have been elusive.
The nitrogen and sulfur K-edge X-ray absorption spectra of aqueous 2-thiopyridone, a model system for excited-state proton transfer in several recent time-resolved measurements, have been simulated from ab initio molecular dynamics. Spectral signatures of the local intra- and inter-molecular structure are identified and rationalized, which facilitates experimental interpretation and optimization. In particular, comparison of aqueous and gas phase spectrum simulations assesses the previously unquantified solvatization effects, where hydrogen bonding is found to yield solvatochromatic shifts up to nearly 1 eV of the main peak positions. Thereby, while each K-edge can still decisively determine the local protonation of its core-excited site, only their combined, complementary fingerprints allow separating all of the three relevant molecular forms, giving a complete picture of the proton transfer.
Over large coastal regions in Greenland and Antarctica the ice sheet calves directly into the ocean. In contrast to ice-shelf calving, an increase in calving from grounded glaciers contributes directly to sea-level rise. Ice cliffs with a glacier freeboard larger than approximate to 100 m are currently not observed, but it has been shown that such ice cliffs are increasingly unstable with increasing ice thickness. This cliff calving can constitute a self-amplifying ice loss mechanism that may significantly alter sea-level projections both of Greenland and Antarctica. Here we seek to derive a minimalist stress-based parametrization for cliff calving from grounded glaciers whose freeboards exceed the 100m stability limit derived in previous studies. This will be an extension of existing calving laws for tidewater glaciers to higher ice cliffs. To this end we compute the stress field for a glacier with a simplified two-dimensional geometry from the two-dimensional Stokes equation. First we assume a constant yield stress to derive the failure region at the glacier front from the stress field within the glacier. Secondly, we assume a constant response time of ice failure due to exceedance of the yield stress. With this strongly constraining but very simple set of assumptions we propose a cliff-calving law where the calving rate follows a power-law dependence on the freeboard of the ice with exponents between 2 and 3, depending on the relative water depth at the calving front. The critical freeboard below which the ice front is stable decreases with increasing relative water depth of the calving front. For a dry water front it is, for example, 75 m. The purpose of this study is not to provide a comprehensive calving law but to derive a particularly simple equation with a transparent and minimalist set of assumptions.
Rotating stellar convection transports angular momentum towards the equator, generating the characteristic equatorial acceleration of the solar rotation while the radial flux of angular momentum is always inwards. New numerical box simulations for the meridional cross-correlation < u(theta)u(phi)>, however, reveal the angular momentum transport towards the poles for slow rotation and towards the equator for fast rotation. The explanation is that for slow rotation a negative radial gradient of the angular velocity always appears, which in combination with a so-far neglected rotation-induced off-diagonal eddy viscosity term nu(perpendicular to) provides "antisolar rotation" laws with a decelerated equator Similarly, the simulations provided positive values for the rotation-induced correlation < u(r)u(theta)>, which is relevant for the resulting latitudinal temperature profiles (cool or warm poles) for slow rotation and negative values for fast rotation. Observations of the differential rotation of slowly rotating stars will therefore lead to a better understanding of the actual stress-strain relation, the heat transport, and the underlying model of the rotating convection.
Degradation of thermal barrier coatings (TBCs) in gas-turbine engines due to calcium-magnesium-aluminosilicate (CMAS) glassy deposits from various sources has been a persistent issue since many years. In this study, state of the art electron microscopy was correlated with X-ray refraction techniques to elucidate the intrusion of CMAS into the porous structure of atmospheric plasma sprayed (APS) TBCs and the formation and growth of cracks under thermal cycling in a burner rig. Results indicate that the sparse nature of the infiltration as well as kinetics in the burner rig are majorly influenced by the wetting behavior of the CMAS. Despite the obvious attack of CMAS on grain boundaries, the interaction of yttria-stabilized zirconia (YSZ) with intruded CMAS has no immediate impact on structure and density of internal surfaces. At a later stage the formation of horizontal cracks is observed in a wider zone of the TBC layer.
Precursor Wave Amplification by Ion-Electron Coupling through Wakefield in Relativistic Shocks
(2019)
We investigated electromagnetic precursor wave emission in relativistic shocks by using two-dimensional particle-in-cell simulations. We found that the wave amplitude is significantly enhanced by a positive feedback process associated with ion-electron coupling through the wakefields for high magnetization. The wakefields collapse during the nonlinear process of the parametric decay instability in the near-upstream region, where nonthermal electrons and ions are generated. The intense coherent emission and the particle acceleration may operate in high-energy astrophysical objects.
We report the discovery of three stars that, along with the prototype LP 40-365, form a distinct class of chemically peculiar runaway stars that are the survivors of thermonuclear explosions. Spectroscopy of the four confirmed LP 40-365 stars finds ONe-dominated atmospheres enriched with remarkably similar amounts of nuclear ashes of partial O- and Si-burning. Kinematic evidence is consistent with ejection from a binary supernova progenitor; at least two stars have rest-frame velocities indicating they are unbound to the Galaxy. With masses and radii ranging between 0.20 and 0.28M(circle dot) and between 0.16 and 0.60 R-circle dot, respectively, we speculate these inflated white dwarfs are the partly burnt remnants of either peculiar Type Iax or electron-capture supernovae. Adopting supernova rates from the literature, we estimate that similar to 20 LP 40-365 stars brighter than 19 mag should be detectable within 2 kpc from the Sun at the end of the Gaia mission. We suggest that as they cool, these stars will evolve in their spectroscopic appearance, and eventually become peculiar O-rich white dwarfs. Finally, we stress that the discovery of new LP 40-365 stars will be useful to further constrain their evolution, supplying key boundary conditions to the modelling of explosion mechanisms, supernova rates, and nucleosynthetic yields of peculiar thermonuclear explosions.
Stellar post asymptotic giant branch (post-AGB) evolution can be completely altered by a final thermal pulse (FTP) which may occur when the star is still leaving the AGB (AFTP), at the departure from the AGB at still constant luminosity (late TP, LTP) or after the entry to the white-dwarf cooling sequence (very late TP, VLTP). Then convection mixes the Herich material with the H-rich envelope. According to stellar evolution models the result is a star with a surface composition of H approximate to 20 per cent by mass (AFTP), approximate to 1 per cent (LTP), or (almost) no H (VLTP). Since FTP stars exhibit intershell material at their surface, spectral analyses establish constraints for AGB nucleosynthesis and stellar evolution. We performed a spectral analysis of the so-called hybrid PG 1159-type central stars (CS) of the planetary nebulae Abell 43 and NGC7094 by means of non-local thermodynamical equilibrium models. We confirm the previously determined effective temperatures of T-eff = 115 000 +/- 5 000K and determine surface gravities of log (g /(cm s(-2))) = 5.6 +/- 0.1 for both. From a comparison with AFTP evolutionary tracks, we derive stellar masses of 0.57(-0.04)(+0.07)M(circle dot) and determine the abundances of H, He, and metals up to Xe. Both CS are likely AFTP stars with a surface H mass fraction of 0.25 +/- 0.03 and 0.15 +/- 0.03, respectively, and an Fe deficiency indicating subsolar initial metallicities. The light metals show typical PG 1159-type abundances and the elemental composition is in good agreement with predictions from AFTP evolutionary models. However, the expansion ages do not agree with evolution time-scales expected from the AFTP scenario and alternatives should be explored.
A comprehensive analysis of the variable temporal and spatial responses of tropical-subtropical high-altitude glaciers to climate change is critical for successful model predictions and environmental risk assessment in the Himalayan-Tibetan orogen. High-frequency Holocene glacier chronostratigraphies are therefore reconstructed in 79 glaciated valleys across the orogen using 519 published and 16 new terrestrial cosmogenic 10Be exposure age dataset. Published 10Be ages are compiled only for moraine boulders (excluding bedrock ages). These ages are recalculated using the latest ICE-D production rate calibration database and the scaling scheme models. Outliers for the individual moraine are detected using the Chauvenet's criterion. In addition, past equilibrium-line altitudes (ELAs) are determined using the area-altitude (AA), area accumulation ratio (AAR), and toe-headwall accumulation ratio (THAR) methods for each glacier advance. The modern maximum elevations of lateral moraines (MELM) are also used to estimate modern ELAs and as an independent check on mean ELAs derived using the above three methods. These data may serve as an essential archive for future studies focusing on the cryospheric and environmental changes in the Himalayan-Tibetan orogen. A more comprehensive analysis of the published and new 10Be ages and ELA results and a list of references are presented in Saha et al. (2019, High-frequency Holocene glacier fluctuations in the Himalayan-Tibetan orogen. Quaternary Science Reviews, 220, 372–400).
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
Photothermal conversion (PTC) nanostructures have great potential for applications in many fields, and therefore, they have attracted tremendous attention. However, the construction of a PTC nanoreactor with multi-compartment structure to achieve the combination of unique chemical properties and structural feature is still challenging due to the synthetic difficulties. Herein, we designed and synthesized a catalytically active, PTC gold (Au)@polydopamine (PDA) nanoreactor driven by infrared irradiation using assembled PS-b-P2VP nanosphere as soft template. The particles exhibit multi-compartment structure which is revealed by 3D electron tomography characterization technique. They feature permeable shells with tunable shell thickness. Full kinetics for the reduction reaction of 4-nitrophenol has been investigated using these particles as nanoreactors and compared with other reported systems. Notably, a remarkable acceleration of the catalytic reaction upon near-infrared irradiation is demonstrated, which reveals for the first time the importance of the synergistic effect of photothermal conversion and complex inner structure to the kinetics of the catalytic reduction. The ease of synthesis and fresh insights into catalysis will promote a new platform for novel nanoreactor studies.
We present new calculations of the mass inflow and outflow rates around the Milky Way (MW), derived from a catalog of ultraviolet metal-line high-velocity clouds (HVCs). These calculations are conducted by transforming the HVC velocities into the Galactic standard of rest (GSR) reference frame, identifying inflowing (vGSR.<.0 km s(-1)) and outflowing (vGSR > 0 km s(-1)) populations, and using observational constraints on the distance, metallicity, dust content, covering fractions, and total silicon column density of each population. After removing HVCs associated with the Magellanic Stream and the Fermi Bubbles, we find inflow and outflow rates in cool (T similar to 10(4) K) ionized gas of dM(in)/dt greater than or similar to.(0.53 +/- 0.23)(d/12 kpc)(Z/0.2Z(circle dot))-1M(circle dot) yr(-1) and dM(out)/dt greater than or similar to (0.16 +/- 0.07)(d/12 kpc)(Z/0.5Z(circle dot))M--1(circle dot) yr(-1). The apparent excess of inflowing over outflowing gas suggests that the MW is currently in an inflow-dominated phase, but the presence of substantial mass flux in both directions supports a Galactic fountain model, in which gas is constantly recycled between the disk and the halo. We also find that the metal flux in both directions (in and out) is indistinguishable. By comparing the outflow rate to the Galactic star formation rate, we present the first estimate of the mass loading factor (eta(HVC)) of the disk-wide MW wind, finding eta(HVC) greater than or similar to (0.10 +/- 0.06)(d/12 kpc)(Z/0.5Z(circle dot))(-1). Including the contributions from low- and intermediatevelocity clouds and from hot gas would increase these inflow and outflow estimates.
Context. Galactic globular clusters (GCs) are now known to host multiple populations displaying particular abundance variations. The different populations within a GC can be well distinguished following their position in the pseudo two-colors diagrams, also referred to as "chromosome maps". These maps are constructed using optical and near-UV photometry available from the Hubble Space Telescope (HST) UV survey of GCs. However, the chemical tagging of the various populations in the chromosome maps is hampered by the fact that HST photometry and elemental abundances are both only available for a limited number of stars. Aims. The spectra collected as part of the MUSE survey of globular clusters provide a spectroscopic counterpart to the HST photometric catalogs covering the central regions of GCs. In this paper, we use the MUSE spectra of 1115 red giant branch (RGB) stars in NGC 2808 to characterize the abundance variations seen in the multiple populations of this cluster. Methods. We used the chromosome map of NGC 2808 to divide the RGB stars into their respective populations. We then combined the spectra of all stars belonging to a given population, resulting in one high signal-to-noise ratio spectrum representative of each population. Results. Variations in the spectral lines of O, Na, Mg, and Al are clearly detected among four of the populations. In order to quantify these variations, we measured equivalent width differences and created synthetic populations spectra that were used to determine abundance variations with respect to the primordial population of the cluster. Our results are in good agreement with the values expected from previous studies based on high-resolution spectroscopy. We do not see any significant variations in the spectral lines of Ca, K, and Ba. We also do not detect abundance variations among the stars belonging to the primordial population of NGC 2808. Conclusions. We demonstrate that in spite of their low resolution, the MUSE spectra can be used to investigate abundance variations in the context of multiple populations.
We utilise multi-epoch MUSE spectroscopy to study binary stars in the core of the Galactic globular cluster NGC 3201. Our sample consists of 3553 stars with 54 883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing radial velocity variations based on the whole inhomogeneous radial velocity sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific cluster we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE field of view and enables us to deduce the underlying true binary frequency of (6.75 +/- 0.72)% in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a large portion of primordial binaries. We can also confirm a radial increase in the binary fraction towards the cluster centre due to mass segregation. We discovered that in the core of NGC 3201 at least (57.5 +/- 7.9)% of blue straggler stars are in a binary system. For the first time in a study of globular clusters, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven blue straggler stars and the connection to SX Phoenicis-type stars. We show evidence that two blue straggler formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole candidates with minimum masses (M sin i) of (7.68 +/- 0.50)M-circle dot, (4.4 +/- 2.8)M-circle dot, and refine the minimum mass estimate on the already published black hole to (4.53 +/- 0.21)M-circle dot, These black holes are consistent with an extensive black hole subsystem hosted by NGC 3201.
New wave frequency and amplitude models for the nightside and dayside chorus waves are built based on measurements from the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrument onboard the Van Allen Probes. The corresponding 3D diffusion coefficients are systematically obtained. Compared with previous commonly-used (typical) parameterizations, the new parameterizations result in differences in diffusion rates that depend on the energy and pitch angle. Furthermore, one-year 3D diffusive simulations are performed using the Versatile Electron Radiation Belt (VERB) code. Both typical and new wave parameterizations simulation results are in a good agreement with observations at 0.9 MeV. However, the new parameterizations for nightside chorus better reproduce the observed electron fluxes. These parameterizations will be incorporated into future modeling efforts.
Most of the matter in the universe consists of hydrogen. The hydrogen in the intergalactic medium (IGM), the matter between the galaxies, underwent a change of its ionisation state at the epoch of reionisation, at a redshift roughly between 6>z>10, or ~10^8 years after the Big Bang. At this time, the mostly neutral hydrogen in the IGM was ionised but the source of the responsible hydrogen ionising emission remains unclear. In this thesis I discuss the most likely candidates for the emission of this ionising radiation, which are a type of galaxy called Lyman alpha emitters (LAEs). As implied by their name, they emit Lyman alpha radiation, produced after a hydrogen atom has been ionised and recombines with a free electron. The ionising radiation itself (also called Lyman continuum emission) which is needed for this process inside the LAEs could also be responsible for ionising the IGM around those galaxies at the epoch of reionisation, given that enough Lyman continuum escapes. Through this mechanism, Lyman alpha and Lyman continuum radiation are closely linked and are both studied to better understand the properties of high redshift galaxies and the reionisation state of the universe.
Before I can analyse their Lyman alpha emission lines and the escape of Lyman continuum emission from them, the first step is the detection and correct classification of LAEs in integral field spectroscopic data, specifically taken with the Multi-Unit Spectroscopic Explorer (MUSE). After detecting emission line objects in the MUSE data, the task of classifying them and determining their redshift is performed with the graphical user interface QtClassify, which I developed during the work on this thesis. It uses the strength of the combination of spectroscopic and photometric information that integral field spectroscopy offers to enable the user to quickly identify the nature of the detected emission lines. The reliable classification of LAEs and determination of their redshifts is a crucial first step towards an analysis of their properties.
Through radiative transfer processes, the properties of the neutral hydrogen clouds in and around LAEs are imprinted on the shape of the Lyman alpha line. Thus after identifying the LAEs in the MUSE data, I analyse the properties of the Lyman alpha emission line, such as the equivalent width (EW) distribution, the asymmetry and width of the line as well as the double peak fraction. I challenge the common method of displaying EW distributions as histograms without taking the limits of the survey into account and construct a more independent EW distribution function that better reflects the properties of the underlying population of galaxies. I illustrate this by comparing the fraction of high EW objects between the two surveys MUSE-Wide and MUSE-Deep, both consisting of MUSE pointings (each with the size of one square arcminute) of different depths. In the 60 MUSE-Wide fields of one hour exposure time I find a fraction of objects with extreme EWs above EW_0>240A of ~20%, while in the MUSE-Deep fields (9 fields with an exposure time of 10 hours and one with an exposure time of 31 hours) I find a fraction of only ~1%, which is due to the differences in the limiting line flux of the surveys. The highest EW I measure is EW_0 = 600.63 +- 110A, which hints at an unusual underlying stellar population, possibly with a very low metallicity.
With the knowledge of the redshifts and positions of the LAEs detected in the MUSE-Wide survey, I also look for Lyman continuum emission coming from these galaxies and analyse the connection between Lyman continuum emission and Lyman alpha emission. I use ancillary Hubble Space Telescope (HST) broadband photometry in the bands that contain the Lyman continuum and find six Lyman continuum leaker candidates. To test whether the Lyman continuum emission of LAEs is coming only from those individual objects or the whole population, I select LAEs that are most promising for the detection of Lyman continuum emission, based on their rest-frame UV continuum and Lyman alpha line shape properties. After this selection, I stack the broadband data of the resulting sample and detect a signal in Lyman continuum with a significance of S/N = 5.5, pointing towards a Lyman continuum escape fraction of ~80%. If the signal is reliable, it strongly favours LAEs as the providers of the hydrogen ionising emission at the epoch of reionisation and beyond.
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.
Diffusion in tilted washboard potentials can paradoxically exceed free normal diffusion. The effect becomes much stronger in the underdamped case due to inertial effects. What happens upon inclusion of usually neglected fractional hydrodynamics memory effects (Basset-Boussinesq frictional force), which result in a heavy algebraic tail of the velocity autocorrelation function of the potential-free diffusion making it transiently superdiffusive? Will a giant enhancement of diffusion become even stronger, and the transient superdiffusion last even longer? These are the questions that we answer in this Letter based on an accurate numerical investigation. We show that a resonancelike enhancement of normal diffusion becomes indeed much stronger and sharper. Moreover, a long-lasting transient regime of superdiffusion, including Richardson-like diffusion, <delta x(2)(t)>proportional to t(3) and ballistic supertransport, <delta x(t)>proportional to t(2), is revealed.
The formation process of zeolite A (Linde Type A) was monitored inline at 1.5 L scale by Photon Density Wave (PDW) spectroscopy as novel process analytical technology for highly turbid liquid suspensions. As a result, the reduced scattering coefficient, being a measure for particle number, size, and morphology, provides distinct process information, including the formation of amorphous particles and their transfer into crystalline zeolite structures. The onset and end of the crystallization process can be detected inline and in real-time. Analyses by powder X-ray diffraction and electron microscopy, based on a sampling approach, support the interpretation of the results obtained by PDW spectroscopy. In addition, the influence of the molar water content was investigated, indicating a linear increase of the time needed to reach the end of the zeolite A crystallization with increasing molar water content. Further experiments indicate a strong influence of the silica source on the course of the crystallization. The applicability of PDW spectroscopy under even more demanding chemical and physical conditions was investigated by monitoring the synthesis of zeolite L (Linde Type L).
We report on the experimental and theoretical interpretation of the diffraction of a probe beam during inscription of a surface relief grating with an interference pattern into a photo-responsive polymer film. For this, we developed a set-up allowing for the simultaneous recording of the diffraction efficiency (DE), the fine structure of the diffraction spot and the topographical changes, in situ and in real time while the film is irradiated. The time dependence of the DE, as the surface relief deepens, follows a Bessel function exhibiting maxima and minima. The size of the probe beam relative to the inscribed grating (i.e., to the size of the writing beams) matters and has to be considered for the interpretation of the DE signal. It is also at the origin of a fine structure within the diffraction spot where ring-shaped features appear once an irradiation time corresponding to the first maximum of the DE has been exceeded.
Brownian motion and beyond: first-passage, power spectrum, non-Gaussianity, and anomalous diffusion
(2019)
Brownian motion is a ubiquitous physical phenomenon across the sciences. After its discovery by Brown and intensive study since the first half of the 20th century, many different aspects of Brownian motion and stochastic processes in general have been addressed in Statistical Physics. In particular, there now exists a very large range of applications of stochastic processes in various disciplines. Here we provide a summary of some of the recent developments in the field of stochastic processes, highlighting both the experimental findings and theoretical frameworks.
We explore electron preacceleration at high-Mach-number nonrelativistic perpendicular shocks at, e.g., young supernova remnants, which are a prerequisite of further acceleration to very high energies via diffusive shock acceleration. Using fully kinetic particle-in-cell simulations of shocks and electron dynamics in them, we investigate the influence of shock-surfing acceleration (SSA) at the shock foot on the nonthermal population of electrons downstream of the shock. The SSA is followed by further energization at the shock ramp where the Weibel instability spawns a type of second-order Fermi acceleration. The combination of these two processes leads to the formation of a nonthermal electron population, but the importance of SSA becomes smaller for larger ion-to-electron mass ratios in the simulation. We discuss the resulting electron spectra and the relevance of our results to the physics of systems with real ion-to-electron mass ratios and fully three-dimensional behavior.
We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded H II regions, active-galactic-nucleus (AGN) narrow-line regions, and radiative shocks computed all using the same physically consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in He II lambda 4686/H beta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed He II lambda 4686/H beta ratios of galaxies with large EW(H beta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias towards massive stars at low metallicity. This would also account for the surprisingly high [O I]/[O III] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.
Polar crown filaments form above the polarity inversion line between the old magnetic flux of the previous cycle and the new magnetic flux of the current cycle. Studying their appearance and their properties can lead to a better understanding of the solar cycle. We use full-disk data of the Chromospheric Telescope (ChroTel) at the Observatorio del Teide, Tenerife, Spain, which were taken in three different chromospheric absorption lines (H alpha lambda 6563 angstrom, Ca II K lambda 3933 angstrom, and He I lambda 10830 angstrom), and we create synoptic maps. In addition, the spectroscopic He I data allow us to compute Doppler velocities and to create synoptic Doppler maps. ChroTel data cover the rising and decaying phase of Solar Cycle 24 on about 1000 days between 2012 and 2018. Based on these data, we automatically extract polar crown filaments with image-processing tools and study their properties. We compare contrast maps of polar crown filaments with those of quiet-Sun filaments. Furthermore, we present a super-synoptic map summarizing the entire ChroTel database. In summary, we provide statistical properties, i.e. number and location of filaments, area, and tilt angle for both the maximum and the declining phase of Solar Cycle 24. This demonstrates that ChroTel provides a promising data set to study the solar cycle.
Storm Time Depletions of Multi-MeV Radiation Belt Electrons Observed at Different Pitch Angles
(2019)
During geomagnetic storms, the rapid depletion of the high-energy (several MeV) outer radiation belt electrons is the result of loss to the interplanetary medium through the magnetopause, outward radial diffusion, and loss to the atmosphere due to wave-particle interactions. We have performed a statistical study of 110 storms using pitch angle resolved electron flux measurements from the Van Allen Probes mission and found that inside of the radiation belt (L* = 3 - 5) the number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 30 degrees is higher than number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 75 degrees. We conclude that this result is consistent with electron scattering by whistler and electromagnetic ion cyclotron waves. At the outer edge of the radiation belt (L* >= 5.2) the number of storms that result in depletion is also large (similar to 40-50%), emphasizing the significance of the magnetopause shadowing effect and outward radial transport.
We consider a simple model for active random walk with general temporal correlations, and investigate the shape of the probability distribution function of the displacement during a short time interval. We find that under certain conditions the distribution exhibits multiple peaks and we show analytically and numerically that the existence of these peaks is governed by the walker?s tendency to move forward, while the correlations between the timing of its active motion control the magnitude and shape of the peaks. In particular, we find that in a homogeneous system such peaks can occur only if the persistence is strong enough.
Digital Volume Correlation (DVC) is a powerful set of techniques used to compute the local shifts of 3D images obtained, for instance, in tomographic experiments. It is utilized to analyze the geometric changes of the investigated object as well as to correct the corresponding image misalignments for further analysis. It can therefore be used to evaluate the local density changes of the same regions of the inspected specimens, which might be shifted between measurements. In recent years, various approaches and corresponding pieces of software were introduced. Accuracies for the computed shift vectors of up to about 1 parts per thousand of a single voxel size have been reported. These results, however, were based either on synthetic datasets or on an unrealistic setup. In this work, we propose two simple methods to evaluate the accuracy of DVC-techniques using more realistic input data and apply them to several DVC programs. We test these methods on three materials (tuff, sandstone, and concrete) that show different contrast and structural features. Published under license by AIP Publishing.
We study the dynamics of the ring of identical phase oscillators with nonlinear nonlocal coupling. Using the Ott - Antonsen approach, the problem is formulated as a system of partial derivative equations for the local complex order parameter. In this framework, we investigate the existence and stability of twisted states. Both fully coherent and partially coherent stable twisted states were found (the latter ones for the first time for identical oscillators). We show that twisted states can be stable starting from a certain critical value of the medium length, or on a length segment. The analytical results are confirmed with direct numerical simulations in finite ensembles.
The actin cytoskeleton and its response to external chemical stimuli is fundamental to the mechano-biology of eukaryotic cells and their functions. One of the key players that governs the dynamics of the actin network is the motor protein myosin II. Based on a phase space embedding we have identified from experiments three phases in the cytoskeletal dynamics of starved Dictyostelium discoideum in response to a precisely controlled chemotactic stimulation. In the first two phases the dynamics of actin and myosin II in the cortex is uncoupled, while in the third phase the time scale for the recovery of cortical actin is determined by the myosin II dynamics. We report a theoretical model that captures the experimental observations quantitatively. The model predicts an increase in the optimal response time of actin with decreasing myosin II-actin coupling strength highlighting the role of myosin II in the robust control of cell contraction.
A stellar census in globular clusters with MUSE: A spectral catalogue of emission-line sources
(2019)
Aims. Globular clusters produce many exotic stars due to a much higher frequency of dynamical interactions in their dense stellar environments. Some of these objects were observed together with several hundred thousand other stars in our MUSE survey of 26 Galactic globular clusters. Assuming that at least a few exotic stars have exotic spectra (i.e. spectra that contain emission lines), we can use this large spectroscopic data set of over a million stellar spectra as a blind survey to detect stellar exotica in globular clusters. Methods. To detect emission lines in each spectrum, we modelled the expected shape of an emission line as a Gaussian curve. This template was used for matched filtering on the di fferences between each observed 1D spectrum and its fitted spectral model. The spectra with the most significant detections of H alpha emission are checked visually and cross-matched with published catalogues. Results. We find 156 stars with H alpha emission, including several known cataclysmic variables (CV) and two new CVs, pulsating variable stars, eclipsing binary stars, the optical counterpart of a known black hole, several probable sub-subgiants and red stragglers, and 21 background emission-line galaxies. We find possible optical counterparts to 39 X-ray sources, as we detected H alpha emission in several spectra of stars that are close to known positions of Chandra X-ray sources. This spectral catalogue can be used to supplement existing or future X-ray or radio observations with spectra of potential optical counterparts to classify the sources.
Lévy flights are paradigmatic generalised random walk processes, in which the independent stationary increments—the 'jump lengths'—are drawn from an -stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering Lévy flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of Lévy flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index and the skewness (asymmetry) parameter . The other approach is based on the stochastic Langevin equation with -stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times.
Precipitation patterns and extremes are significantly influenced by various climatic factors and large-scale atmospheric circulation patterns. This study uses wavelet coherence analysis to detect significant interannual and interdecadal oscillations in monthly precipitation extremes across India and their teleconnections to three prominent climate indices, namely, Nino 3.4, Pacific Decadal Oscillation, and Indian Ocean Dipole (IOD). Further, partial wavelet coherence analysis is used to estimate the standalone relationship between the climate indices and precipitation after removing the effect of interdependency. The wavelet analysis of monthly precipitation extremes at 30 different locations across India reveals that (a) interannual (2-8 years) and interdecadal (8-32 years) oscillations are statistically significant, and (b) the oscillations vary in both time and space. The results from the partial wavelet coherence analysis reveal that Nino 3.4 and IOD are the significant drivers of Indian precipitation at interannual and interdecadal scales. Intriguingly, the study also confirms that the strength of influence of large-scale atmospheric circulation patterns on Indian precipitation extremes varies with spatial physiography of the region.
Understanding and disentangling photophysical properties of long-lived photoexcitations in bulk heterojunction (BHJ) solar cells, which contribute mostly to photocurrent, provide essential guidelines to their improvement. However, to construct improved physical models, their rational design relies on reliable measurement techniques for charge recombination. Here, we combine photocurrent and photoinduced absorption spectroscopy (PCPIA) to directly probe the free carrier concentration and investigate loss mechanisms of long-lived excitations in nearly 10% efficient PPDT2FBT/PC70BM BHJ solar cells under steady-state operational conditions. From the PCPIA data obtained under open- circuit and short-circuit conditions, the absorption cross section and the concentration of photoexcitations are obtained. This material system exhibits an exceptionally low bimolecular recombination rate, about 300 times smaller than the diffusion-controlled electron and hole encounter rate. Furthermore, we observe that the fill factor is limited by losses originating from long-lived photoexcitations undergoing dispersive bimolecular recombination.
We predict the recoil-induced molecular dissociation in hard-x-ray photoionization. The recoil effect is caused by electronic and photon momentum exchange with the molecule. We show the strong role of relativistic effects for the studied molecular fragmentation. The recoil-induced fragmentation of the molecule is caused by elongation of the bond due to the vibrational recoil effect and because of the centrifugal force caused by the rotational recoil. The calculations of the x-ray photoelectron spectra of the H-2 and NO molecules show that the predicted effects can be observed in high-energy synchrotrons like SOLEIL, SPring-8, PETRA, and XFEL SACLA. The relativistic effect enhances the recoil momentum transfer and makes it strongly sensitive to the direction of ejection of the fast photoelectron with respect to the photon momentum.
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.
In this paper, we explore the time-energy domain quantum-classical transition comparing a classical pump-probe experiment on a diatomic molecule to its quantum enhanced counterpart, where the pump and probe pulses are substituted by the signal and idler beams of a spontaneous parametric down conversion (SPDC) source. Absorption of biphotons generated with SPDC exploits quantum time-energy entanglement to enhance the overall yield and selectivity of the process, when compared with a classical pump-probe setup, while maintaining femtosecond time resolution. We systematically study the effects of correlation strength on process efficiency and selectivity, comparing the results to classical pump-probe spectra. An excitation scheme to improve the yield based on spectral narrowing of biphotons is shown. The results indicate that the quantum improvements in yield are caused by a more efficient use of the total power available for the process.
In this Letter, we propose that the x-ray and the TeV observations in the vicinity of Geminga can be understood in the framework of anisotropic diffusion of injected electrons or positrons. This interpretation only requires the turbulence in the vicinity of Geminga to be sub-Alfvenic with the local mean magnetic field direction approximately aligned with our line of sight towards Geminga, without invoking extreme conditions for the environment, such as an extremely small diffusion coefficient and a weak magnetic field of submicrogauss as suggested in previous literature.
Charge extraction methods are popular for measuring the charge carrier density in thin film organic solar cells and to draw conclusions about the order and coefficient of nongeminate charge recombination. However, results from such studies may be falsified by inhomogeneous steady state carrier profiles or surface recombination. Here, we present a detailed drift-diffusion study of two charge extraction methods, bias-assisted charge extraction (BACE) and time-delayed collection field (TDCF). Simulations are performed over a wide range of the relevant parameters. Our simulations reveal that both charge extraction methods provide reliable information about the recombination order and coefficient if the measurements are performed under appropriate conditions. However, results from BACE measurements may be easily affected by surface recombination, in particular for small active layer thicknesses and low illumination densities. TDCF, on the other hand, is more robust against surface recombination due to its transient nature but also because it allows for a homogeneous high carrier density to be inserted into the active layer. Therefore, TDCF is capable to provide meaningful information on the order and coefficient of recombination even if the model conditions are not exactly fulfilled. We demonstrate this for an only 100 nm thick layer of a highly efficient nonfullerene acceptor (NFA) blend, comprising the donor polymer PM6 and the NFA Y6. TDCF measurements were performed as a function of delay time for different laser fluences and bias conditions. The full set of data could be consistently fitted by a strict second order recombination process, with a bias- and fluence-independent bimolecular recombination coefficient k(2) = 1.7 x 10(-17)m(3) s(-1). BACE measurements performed on the very same layer yielded the identical result, despite the very different excitation conditions. This proves that recombination in this blend is mostly through processes in the bulk and that surface recombination is of minor importance despite the small active layer thickness. Published under license by AIP Publishing.
Near ambient pressure - x-ray photoelectron spectroscopy (NAP-XPS) is a less traditional form of XPS that allows samples to be analyzed at relatively high pressures, i.e., at 2500Pa or higher. With NAP-XPS, one can analyze moderately volatile liquids, biological samples, porous materials, and/or polymeric materials that outgas significantly. In this submission we show C 1s, O 1s, and survey NAP-XPS spectra from poly(L-lactic acid). The C 1s and O 1s envelopes were fit with three and two Gaussian-Lorentzian sum functions, respectively. Water vapor (800Pa) was used as the residual gas for charge compensation, which was confirmed by the sharp peak at 535.0 eV in the O 1s narrow scan. The uniqueness plot corresponding to the C 1s fit shows that the fit parameters had statistical significance. C 1s and O 1s spectra of PLLA damaged by exposure to x-rays for ca. 1 hour are also included. Published by the AVS.
We propose an efficient method for demodulation of phase modulated signals via iterated Hilbert transform embeddings. We show that while a usual approach based on one application of the Hilbert transform provides only an approximation to a proper phase, with iterations the accuracy is essentially improved, up to precision limited mainly by discretization effects. We demonstrate that the method is applicable to arbitrarily complex waveforms, and to modulations fast compared to the basic frequency. Furthermore, we develop a perturbative theory applicable to a simple cosine waveform, showing convergence of the technique.
Dielectric Relaxation Spectroscopy (DRS) and Thermally Stimulated Depolarization Current (TSDC) measurements were employed to study dielectric-relaxation processes, structural transitions and electric-polarization phenomena in poly(vinylidenefluoride-trifluoroethylene-chlorofluoroethylene) (P(VDF-TrFE-CFE)) terpolymer films. Results from DRS confirm the existence of two separate dispersion regions related to a para-to-ferroelectric phase transition and to the glass transition. The dipolar TSDC peak correlates with the loss peak of the ? relaxation that represents the glass transition. The electric polarization calculated from the dipolar TSDC peak (glass transition) shows a non-linear electric-field dependence and saturates at high electric poling fields. As the observed behaviour is essentially the same as that of the electric polarization obtained from direct polarization-versus-electric-field hysteresis measurements, TSDC experiments are also suitable for studying the polarization in relaxor-ferroelectric polymers. A saturation polarization of 44 mC m(?2) was found for an electric field of 190 MV m(?1).
The stability of conducting Taylor-Couette flows under the presence of toroidal magnetic background fields is considered. For strong enough magnetic amplitudes such magnetohydrodynamic flows are unstable against non-axisymmetric perturbations which may also transport angular momentum. In accordance with the often used diffusion approximation, one expects the angular momentum transport to be vanishing for rigid rotation. In the sense of a non-diffusive Lambda effect, however, even for rigidly rotating z-pinches, an axisymmetric angular momentum flux appears which is directed outward (inward) for large (small) magnetic Mach numbers. The internal rotation in a magnetized rotating tank can thus never be uniform. Those particular rotation laws are used to estimate the value of the instability-induced eddy viscosity for which the non-diffusive Lambda effect and the diffusive shear-induced transport compensate each other. The results provide the Shakura & Sunyaev viscosity ansatz leading to numerical values linearly growing with the applied magnetic field.
X-ray free-electron lasers have, over the past decade, opened up the possibility of understanding the ultrafast response of matter to intense X-ray pulses. In earlier research on atoms and small molecules, new aspects of this response were uncovered, such as rapid sequences of inner-shell photoionization and Auger ionization. Here, we studied a larger molecule, buckminsterfullerene (C-60), exposed to 640 eV X-rays, and examined the role of chemical effects, such as chemical bonds and charge transfer, on the fragmentation following multiple ionization of the molecule. To provide time resolution, we performed femtosecond-resolved X-ray pump/X-ray probe measurements, which were accompanied by advanced simulations. The simulations and experiment reveal that despite substantial ionization induced by the ultrashort (20 fs) X-ray pump pulse, the fragmentation of C-60 is considerably delayed. This work uncovers the persistence of the molecular structure of C-60, which hinders fragmentation over a timescale of hundreds of femtoseconds. Furthermore, we demonstrate that a substantial fraction of the ejected fragments are neutral carbon atoms. These findings provide insights into X-ray free-electron laser-induced radiation damage in large molecules, including biomolecules.
Optimizing the photoluminescence (PL) yield of a solar cell has long been recognized as a key principle to maximize the power conversion efficiency. While PL measurements are routinely applied to perovskite films and solar cells under open circuit conditions (V-OC), it remains unclear how the emission depends on the applied voltage. Here, we performed PL(V) measurements on perovskite cells with different hole transport layer thicknesses and doping concentrations, resulting in remarkably different fill factors (FFs). The results reveal that PL(V) mirrors the current-voltage (JV) characteristics in the power-generating regime, which highlights an interesting correlation between radiative and nonradiative recombination losses. In particular, high FF devices show a rapid quenching of PL(V) from open-circuit to the maximum power point. We conclude that, while the PL has to be maximized at V-OC at lower biases < V-OC the PL must be rapidly quenched as charges need to be extracted prior to recombination.
We present results from our on-going MusE GAs FLOw and Wind (MEGAFLOW) survey, which consists of 22 quasar lines of sight, each observed with the integral field unit MUSE and the UVES spectrograph at the ESO Very Large Telescopes (VLT). The goals of this survey are to study the properties of the circumgalactic medium around z similar to 1 star-forming galaxies. The absorption-line selected survey consists of 79 strong MgII absorbers (with rest-frame equivalent width greater than or similar to 0.3 angstrom) and, currently, 86 associated galaxies within 100 projected kpc of the quasar with stellar masses (M-star) from 109 to 1011 M-circle dot. We find that the cool halo gas traced by MgII is not isotropically distributed around these galaxies from the strong bi-modal distribution in the azimuthal angle of the apparent location of the quasar with respect to the galaxy major axis. This supports a scenario in which outflows are bi-conical in nature and co-exist with a co-planar gaseous structure extending at least up to 60-80 kpc. Assuming that absorbers near the minor axis probe outflows, the current MEGAFLOW sample allowed us to select 26 galaxy-quasar pairs suitable for studying winds. From this sample, using a simple geometrical model, we find that the outflow velocity only exceeds the escape velocity when M-star less than or similar to 4 x 10(9) M-circle dot, implying the cool material is likely to fall back except in the smallest haloes. Finally, we find that the mass loading factor., the ratio between the ejected mass rate and the star formation rate, appears to be roughly constant with respect to the galaxy mass.
Context. The TESS satellite was launched in 2018 to perform high-precision photometry from space over almost the whole sky in a search for exoplanets orbiting bright stars. This instrument has opened new opportunities to study variable hot subdwarfs, white dwarfs, and related compact objects. Targets of interest include white dwarf and hot subdwarf pulsators, both carrying high potential for asteroseismology. Aims. We present the discovery and detailed asteroseismic analysis of a new g-mode hot B subdwarf (sdB) pulsator, EC 21494-7018 (TIC 278659026), monitored in TESS first sector using 120-s cadence. Methods. The TESS light curve was analyzed with standard prewhitening techniques, followed by forward modeling using our latest generation of sdB models developed for asteroseismic investigations. By simultaneously best-matching all the observed frequencies with those computed from models, we identified the pulsation modes detected and, more importantly, we determined the global parameters and structural configuration of the star. Results. The light curve analysis reveals that EC 21494-7018 is a sdB pulsator counting up to 20 frequencies associated with independent g-modes. The seismic analysis singles out an optimal model solution in full agreement with independent measurements provided by spectroscopy (atmospheric parameters derived from model atmospheres) and astrometry (distance evaluated from Gaia DR2 trigonometric parallax). Several key parameters of the star are derived. Its mass (0.391 +/- 0.009x2006;M-circle dot) is significantly lower than the typical mass of sdB stars and suggests that its progenitor has not undergone the He-core flash; therefore this progenitor could originate from a massive (greater than or similar to 2;M-circle dot) red giant, which is an alternative channel for the formation of sdBs. Other derived parameters include the H-rich envelope mass (0.0037 +/- 0.0010;M-circle dot), radius (0.1694 +/- 0.0081;R-circle dot), and luminosity (8.2 +/- 1.1;L-circle dot). The optimal model fit has a double-layered He+H composition profile, which we interpret as an incomplete but ongoing process of gravitational settling of helium at the bottom of a thick H-rich envelope. Moreover, the derived properties of the core indicate that EC 21494-7018 has burnt similar to 43% (in mass) of its central helium and possesses a relatively large mixed core (M-core;=;0.198 +/- 0.010;M-circle dot), in line with trends already uncovered from other g-mode sdB pulsators analyzed with asteroseismology. Finally, we obtain for the first time an estimate of the amount of oxygen (in mass; X(O)(core) = 0.16(-0.05)(+0.13)X(O)core=0.16-0.05+0.13$ X(mathrm{O})_{mathrm{core}}=0.16_{-0.05}<^>{+0.13} $) produced at this stage of evolution by an helium-burning core. This result, along with the core-size estimate, is an interesting constraint that may help to narrow down the still uncertain C-12(alpha,;gamma)O-16 nuclear reaction rate.
The movement of microswimmers is often described by active Brownian particle models. Here we introduce a variant of these models with several internal states of the swimmer to describe stochastic strategies for directional swimming such as run and tumble or run and reverse that are used by microorganisms for chemotaxis. The model includes a mechanism to generate a directional bias for chemotaxis and interactions with external fields (e.g., gravity, magnetic field, fluid flow) that impose forces or torques on the swimmer. We show how this modified model can be applied to various scenarios: First, the run and tumble motion of E. coli is used to establish a paradigm for chemotaxis and investigate how it is affected by external forces. Then, we study magneto-aerotaxis in magnetotactic bacteria, which is biased not only by an oxygen gradient towards a preferred concentration, but also by magnetic fields, which exert a torque on an intracellular chain of magnets. We study the competition of magnetic alignment with active reorientation and show that the magnetic orientation can improve chemotaxis and thereby provide an advantage to the bacteria, even at rather large inclination angles of the magnetic field relative to the oxygen gradient, a case reminiscent of what is expected for the bacteria at or close to the equator. The highest gain in chemotactic velocity is obtained for run and tumble with a magnetic field parallel to the gradient, but in general a mechanism for reverse motion is necessary to swim against the magnetic field and a run and reverse strategy is more advantageous in the presence of a magnetic torque. This finding is consistent with observations that the dominant mode of directional changes in magnetotactic bacteria is reversal rather than tumbles. Moreover, it provides guidance for the design of future magnetic biohybrid swimmers. Author summary In this paper, we propose a modified Active Brownian particle model to describe bacterial swimming behavior under the influence of external forces and torques, in particular of a magnetic torque. This type of interaction is particularly important for magnetic biohybrids (i.e. motile bacteria coupled to a synthetic magnetic component) and for magnetotactic bacteria (i.e. bacteria with a natural intracellular magnetic chain), which perform chemotaxis to swim along chemical gradients, but are also directed by an external magnetic field. The model allows us to investigate the benefits and disadvantages of such coupling between two different directionality mechanisms. In particular we show that the magnetic torque can speed chemotaxis up in some conditions, while it can hinder it in other cases. In addition to an understanding of the swimming strategies of naturally magnetotactic organisms, the results may guide the design of future biomedical devices.
Near-ambient pressure x-ray photoelectron spectroscopy (NAP-XPS) is a less traditional form of XPS that allows samples to be analyzed at relatively high pressures, i. e., at greater than 2500 Pa. In this study, poly(.- benzyl L- glutamate) (PBLG) with a molar mass of 11.3 kg/mol was analyzed by NAP-XPS; here, we show the survey, C 1s, N 1s, and O 1s narrow scans of PBLG. The C 1s peak envelope was fitted in three different ways, to five, six, or seven synthetic peaks. In each fit, there was also a shake-up signal. The O 1s narrow scan was well fit with three peaks: CZO and CvO in a 1:2 ratio from the polymer, and a higher energy signal from water vapor. Hartree-Fock orbital energies of a model monomer served as a guide to an additional fit of the C 1s envelope.
Germanium oxide is a prototype network-forming oxide with pressure-induced structural changes similar to those found in crystals and amorphous silicate oxides at high pressure. Studying density and coordination changes in amorphous GeO2 allows for insight into structural changes in silicate oxides at very high pressure, with implications for the properties of planetary magmas. Here, we report the density of germanium oxide glass up to 133 GPa using the x-ray absorption technique, with very good agreement with previous experimental data at pressure below 40 GPa and recent calculation up to 140 GPa. Our data highlight four distinct compressibility domains, corresponding to changes of the local structure of GeO2. Above 80 GPa, our density data show a compressibility and bulk modulus similar to the counterpart crystal phase, and we propose that a compact distorted sixfold coordination, similar to the structural motif of the pyritelike crystalline GeO2 polymorph, is likely to be stable in that pressure range. Our density data point to a smooth continuous evolution of the average coordination for pressure above 20 GPa with persistent sixfold coordination, without sharp density or density slope discontinuities. These observations are in very good agreement with theoretical calculations and spectroscopic measurements, and our results indicate that glasses and melts may behave similarly to their high-pressure solid counterparts with comparable densities, compressibility, and possibly average coordination.
Border effect corrections for diagonal line based recurrence quantification analysis measures
(2019)
Recurrence Quantification Analysis (RQA) defines a number of quantifiers, which base upon diagonal line structures in the recurrence plot (RP). Due to the finite size of an RP, these lines can be cut by the borders of the RP and, thus, bias the length distribution of diagonal lines and, consequently, the line based RQA measures. In this letter we investigate the impact of the mentioned border effects and of the thickening of diagonal lines in an RP (caused by tangential motion) on the estimation of the diagonal line length distribution, quantified by its entropy. Although a relation to the Lyapunov spectrum is theoretically expected, the mentioned entropy yields contradictory results in many studies. Here we summarize correction schemes for both, the border effects and the tangential motion and systematically compare them to methods from the literature. We show that these corrections lead to the expected behavior of the diagonal line length entropy, in particular meaning zero values in case of a regular motion and positive values for chaotic motion. Moreover, we test these methods under noisy conditions, in order to supply practical tools for applied statistical research.
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.
In this paper we investigate the magneto-optical (MO) and magnetic properties of bismuth iron garnet Bi3Fe5O12 thin films over a wide range of photon energies (1.6-3.5 eV) and temperatures (5-740 K). Depending on the photon energy range, the Faraday rotation (Theta(F)) and ellipticity (epsilon(F)) vary nonmonotonously with temperature. This behavior cannot be explained by a magnetization variation that can only decrease with increasing temperature. Theta(F) and epsilon(F) spectra have therefore been analyzed using a model based on two optical transitions of a diamagnetic nature, representing the tetrahedral and octahedral iron sites. Thus, the contribution of each magnetic sublattice has been extracted from the global macroscopic MO response and investigated as a function of temperature. The magnetic properties of octahedral and tetrahedral sublattices depend differently on temperature, suggesting a different anisotropy due to oxygen coordination. We have demonstrated that this relatively simple macroscopic measurement with a subsequent analysis can grant access to the information on the properties at a microscopic level. These results can advance the fundamental understanding of MO properties in multisublattice magnetic materials.
Phase reduction is a general tool widely used to describe forced and interacting self-sustained oscillators. Here, we explore the phase coupling functions beyond the usual first-order approximation in the strength of the force. Taking the periodically forced Stuart-Landau oscillator as the paradigmatic model, we determine and numerically analyse the coupling functions up to the fourth order in the force strength. We show that the found nonlinear phase coupling functions can be used for predicting synchronization regions of the forced oscillator.
We introduce a versatile mechanism of light-driven self-propelled motion applied to porous Janus-type particles. The mechanism is based on the generation of local light-driven diffusio-osmotic (l-LDDO) flow around each single porous particle subjected to suitable irradiation conditions. The photosensitivity is introduced by a cationic azobenzene containing surfactant, which undergoes a photoisomerization reaction from a more hydrophobic trans-state to a rather hydrophilic cis-state under illumination with light. The negatively charged porous silica particles are dispersed in a corresponding aqueous solution and absorb molecules in their trans-state but expel them in their cis-state. During illumination with blue light triggering both trans-cis and cis-trans isomerization at the same time, the colloids start to move due to the generation of a steady-state diffusive flow of cis-isomers out of and trans-isomers into the particle. This is because a hemi-spherical metal cap partially sealing the colloid breaks the symmetry of the otherwise radially directed local flow around the particle, leading to self-propelled motion. Janus particles exhibit superdiffusive motion with a velocity of similar to 0.5 mu m/s and a persistence length of ca. 50 mu m, confined to microchannels the direction can be maintained up to 300 mu m before rotational diffusion reverts it. Particles forming dimers of different shapes can be made to travel along circular trajectories. The unique feature of this mechanism is that the strength of self-propulsion can be tuned by convenient external optical stimuli (intensity and irradiation wavelength) such that a broad variety of experimental situations can be realized in a spatiotemporal way and in situ.
Additive Manufacturing (AM) in terms of laser powder-bed fusion (L-PBF) offers new prospects regarding the design of parts and enables therefore the production of lattice structures. These lattice structures shall be implemented in various industrial applications (e.g. gas turbines) for reasons of material savings or cooling channels. However, internal defects, residual stress, and structural deviations from the nominal geometry are unavoidable.
In this work, the structural integrity of lattice structures manufactured by means of L-PBF was non-destructively investigated on a multiscale approach.
A workflow for quantitative 3D powder analysis in terms of particle size, particle shape, particle porosity, inter-particle distance and packing density was established. Synchrotron computed tomography (CT) was used to correlate the packing density with the particle size and particle shape. It was also observed that at least about 50% of the powder porosity was released during production of the struts.
Struts are the component of lattice structures and were investigated by means of laboratory CT. The focus was on the influence of the build angle on part porosity and surface quality. The surface topography analysis was advanced by the quantitative characterisation of re-entrant surface features. This characterisation was compared with conventional surface parameters showing their complementary information, but also the need for AM specific surface parameters.
The mechanical behaviour of the lattice structure was investigated with in-situ CT under compression and successive digital volume correlation (DVC). The deformation was found to be knot-dominated, and therefore the lattice folds unit cell layer wise.
The residual stress was determined experimentally for the first time in such lattice structures. Neutron diffraction was used for the non-destructive 3D stress investigation. The principal stress directions and values were determined in dependence of the number of measured directions. While a significant uni-axial stress state was found in the strut, a more hydrostatic stress state was found in the knot. In both cases, strut and knot, seven directions were at least needed to find reliable principal stress directions.
Theoretical Fluid Dynamics
(2019)
If not oriented perfectly isotropically, the strong dipole moment of polar organic semiconductor materials such as tris-(8-hydroxyquinolate)aluminum (Alq3) will lead to the buildup of a giant surface potential (GSP) and thus to a macroscopic dielectric polarization of the organic film. Despite this having been a known fact for years, the implications of such high potentials within an organic layer stack have only been studied recently. In this work, the influence of the GSP on hole injection into organic layers is investigated. Therefore, we apply a concept called dipolar doping to devices consisting of the prototypical organic materials N,N′-Di(1-naphthyl)-N,N′-diphenyl-(1,1′-biphenyl)-4,4′-diamine (NPB) as nonpolar host and Alq3 as dipolar dopant with different mixing ratios to tune the GSP. The mixtures are investigated in single-layer monopolar devices as well as bilayer metal/insulator/semiconductor structures. Characterization is done electrically using current-voltage (I-V) characteristics, impedance spectroscopy, and charge extraction by linearly increasing voltage and time of flight, as well as with ultraviolet photoelectron spectroscopy. We find a maximum in device performance for moderate to low doping concentrations of the polar species in the host. The observed behavior can be described on the basis of the Schottky effect for image-force barrier lowering, if the changes in the interface dipole, the carrier mobility, and the GSP induced by dipolar doping are taken into account.
The possibility to manufacture perovskite solar cells (PSCs) at low temperatures paves the way to flexible and lightweight photovoltaic (PV) devices manufactured via high-throughput roll-to-roll processes. In order to achieve higher power conversion efficiencies, it is necessary to approach the radiative limit via suppression of non-radiative recombination losses. Herein, we performed a systematic voltage loss analysis for a typical low-temperature processed, flexible PSC in n-i-p configuration using vacuum deposited C-60 as electron transport layer (ETL) and two-step hybrid vacuum-solution deposition for CH3NH3PbI3 perovskite absorber. We identified the ETL/absorber interface as a bottleneck in relation to non-radiative recombination losses, the quasi-Fermi level splitting (QFLS) decreases from similar to 1.23 eV for the bare absorber, just similar to 90 meV below the radiative limit, to similar to 1.10 eV when C-60 is used as ETL. To effectively mitigate these voltage losses, we investigated different interfacial modifications via vacuum deposited interlayers (BCP, B4PyMPM, 3TPYMB, and LiF). An improvement in QFLS of similar to 30-40 meV is observed after interlayer deposition and confirmed by comparable improvements in the open-circuit voltage after implementation of these interfacial modifications in flexible PSCs. Further investigations on absorber/hole transport layer (HTL) interface point out the detrimental role of dopants in Spiro-OMeTAD film (widely employed HTL in the community) as recombination centers upon oxidation and light exposure. [GRAPHICS] .
Electrets are dielectrics with quasi-permanent electric charge and/or dipoles, sometimes can be regarded as an electric analogy to a magnet. Since the discovery of the excellent charge retention capacity of poly(tetrafluoro ethylene) and the invention of the electret microphone, electrets have grown out of a scientific curiosity to an important application both in science and technology. The history of electret research goes hand in hand with the quest for new materials with better capacity at charge and/or dipole retention. To be useful, electrets normally have to be charged/poled to render them electro-active. This process involves electric-charge deposition and/or electric dipole orientation within the dielectrics ` surfaces and bulk. Knowledge of the spatial distribution of electric charge and/or dipole polarization after their deposition and subsequent decay is crucial in the task to improve their stability in the dielectrics.
Likewise, for dielectrics used in electrical insulation applications, there are also needs for accumulated space-charge and polarization spatial profiling. Traditionally, space-charge accumulation and large dipole polarization within insulating dielectrics is considered undesirable and harmful to the insulating dielectrics as they might cause dielectric loss and could lead to internal electric field distortion and local field enhancement. High local electric field could trigger several aging processes and reduce the insulating dielectrics' lifetime. However, with the advent of high-voltage DC transmission and high-voltage capacitor for energy storage, these are no longer the case. There are some overlapped between the two fields of electrets and electric insulation. While quasi-permanently trapped electric-charge and/or large remanent dipole polarization are the requisites for electret operation, stably trapped electric charge in electric insulation helps reduce electric charge transport and overall reduced electric conductivity. Controlled charge trapping can help in preventing further charge injection and accumulation as well as serving as field grading purpose in insulating dielectrics whereas large dipole polarization can be utilized in energy storage applications.
In this thesis, the Piezoelectrically-generated Pressure Steps (PPSs) were employed as a nondestructive method to probe the electric-charge and dipole polarization distribution in a range of thin film (several hundred micron) polymer-based materials, namely polypropylene (PP), low-density polyethylene/magnesium oxide (LDPE/MgO) nanocomposites and poly(vinylidene fluoride-co- trifluoro ethylene) (P(VDF-TrFE)) copolymer. PP film surface-treated with phosphoric acid to introduce surfacial isolated nanostructures serves as example of 2-dimensional nano-composites whereas LDPE/MgO serves as the case of 3-dimensional nano-composites with MgO nano-particles dispersed in LDPE polymer matrix. It is evidenced that the nanoparticles on the surface of acid-treated PP and in the bulk of LDPE/MgO nanocomposites improve charge trapping capacity of the respective material and prevent further charge injection and transport and that the enhanced charge trapping capacity makes PP and LDPE/MgO nanocomposites potential materials for both electret and electrical insulation applications. As for PVDF and VDF-based copolymers, the remanent spatial polarization distribution depends critically on poling method as well as specific parameters used in the respective poling method. In this work, homogeneous polarization poling of P(VDF-TrFE) copolymers with different VDF-contents have been attempted with hysteresis cyclical poling. The behaviour of remanent polarization growth and spatial polarization distribution are reported and discussed. The Piezoelectrically-generated Pressure Steps (PPSs) method has proven as a powerful method for the charge storage and transport characterization of a wide range of polymer material from nonpolar, to polar, to polymer nanocomposites category.
In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.
In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.
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.
Due to advances in science and technology towards smaller and more powerful processing units, the fabrication of micrometer sized machines for different tasks becomes more and more possible. Such micro-robots could revolutionize medical treatment of diseases and shall support to work on other small machines. Nevertheless, scaling down robots and other devices is a challenging task and will probably remain limited in near future. Over the past decade the concept of bio-hybrid systems has proved to be a promising approach in order to advance the further development of micro-robots. Bio-hybrid systems combine biological cells with artificial components, thereby benefiting from the functionality of living biological cells. Cell-driven micro-transport is one of the most prominent applications in the emerging field of these systems. So far, micrometer sized cargo has been successfully transported by means of swimming bacterial cells. The potential of motile adherent cells as transport systems has largely remained unexplored.
This thesis concentrates on the social amoeba Dictyostelium discoideum as a potential candidate for an amoeboid bio-hybrid transport system. The use of this model organism comes with several advantages. Due to the unspecific properties of Dictyostelium adhesion, a wide range of different cargo materials can be used for transport. As amoeboid cells exceed bacterial cells in size by one order of magnitude, also the size of an object carried by a single cell can also be much larger for an amoeba. Finally it is possible to guide the cell-driven transport based on the chemotactic behavior of the amoeba. Since cells undergo a developmentally induced chemotactic aggregation, cargo could be assembled in a self-organized manner into a cluster. It is also possible to impose an external chemical gradient to guide the amoeboid transport system to a desired location.
To establish Dictyostelium discoideum as a possible candidate for bio-hybrid transport systems, this thesis will first investigate the movement of single cells. Secondly, the interaction of cargo and cells will be studied. Eventually, a conceptional proof will be conducted, that the cheomtactic behavior can be exploited either to transport a cargo self-organized or through an external chemical source.
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.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.
Quorum-sensing bacteria in a growing colony of cells send out signalling molecules (so-called “autoinducers”) and themselves sense the autoinducer concentration in their vicinity. Once—due to increased local cell density inside a “cluster” of the growing colony—the concentration of autoinducers exceeds a threshold value, cells in this clusters get “induced” into a communal, multi-cell biofilm-forming mode in a cluster-wide burst event. We analyse quantitatively the influence of spatial disorder, the local heterogeneity of the spatial distribution of cells in the colony, and additional physical parameters such as the autoinducer signal range on the induction dynamics of the cell colony. Spatial inhomogeneity with higher local cell concentrations in clusters leads to earlier but more localised induction events, while homogeneous distributions lead to comparatively delayed but more concerted induction of the cell colony, and, thus, a behaviour close to the mean-field dynamics. We quantify the induction dynamics with quantifiers such as the time series of induction events and burst sizes, the grouping into induction families, and the mean autoinducer concentration levels. Consequences for different scenarios of biofilm growth are discussed, providing possible cues for biofilm control in both health care and biotechnology.
Many studies on biological and soft matter systems report the joint presence of a linear mean-squared displacement and a non-Gaussian probability density exhibiting, for instance, exponential or stretched-Gaussian tails. This phenomenon is ascribed to the heterogeneity of the medium and is captured by random parameter models such as ‘superstatistics’ or ‘diffusing diffusivity’. Independently, scientists working in the area of time series analysis and statistics have studied a class of discrete-time processes with similar properties, namely, random coefficient autoregressive models. In this work we try to reconcile these two approaches and thus provide a bridge between physical stochastic processes and autoregressive models.Westart from the basic Langevin equation of motion with time-varying damping or diffusion coefficients and establish the link to random coefficient autoregressive processes. By exploring that link we gain access to efficient statistical methods which can help to identify data exhibiting Brownian yet non-Gaussian diffusion.
Many studies on biological and soft matter systems report the joint presence of a linear mean-squared displacement and a non-Gaussian probability density exhibiting, for instance, exponential or stretched-Gaussian tails. This phenomenon is ascribed to the heterogeneity of the medium and is captured by random parameter models such as ‘superstatistics’ or ‘diffusing diffusivity’. Independently, scientists working in the area of time series analysis and statistics have studied a class of discrete-time processes with similar properties, namely, random coefficient autoregressive models. In this work we try to reconcile these two approaches and thus provide a bridge between physical stochastic processes and autoregressive models.Westart from the basic Langevin equation of motion with time-varying damping or diffusion coefficients and establish the link to random coefficient autoregressive processes. By exploring that link we gain access to efficient statistical methods which can help to identify data exhibiting Brownian yet non-Gaussian diffusion.
Supermassive black holes reside in the hearts of almost all massive galaxies. Their evolutionary path seems to be strongly linked to the evolution of their host galaxies, as implied by several empirical relations between the black hole mass (M BH ) and different host galaxy properties. The physical driver of this co-evolution is, however, still not understood. More mass measurements over homogeneous samples and a detailed understanding of systematic uncertainties are required to fathom the origin of the scaling relations.
In this thesis, I present the mass estimations of supermassive black holes in the nuclei of one late-type and thirteen early-type galaxies. Our SMASHING sample extends from the intermediate to the massive galaxy mass regime and was selected to fill in gaps in number of galaxies along the scaling relations. All galaxies were observed at high spatial resolution, making use of the adaptive-optics mode of integral field unit (IFU) instruments on state-of-the-art telescopes (SINFONI, NIFS, MUSE). I extracted the stellar kinematics from these observations and constructed dynamical Jeans and Schwarzschild models to estimate the mass of the central black holes robustly. My new mass estimates increase the number of early-type galaxies with measured black hole masses by 15%. The seven measured galaxies with nuclear light deficits (’cores’) augment the sample of cored galaxies with measured black holes by 40%. Next to determining massive black hole masses, evaluating the accuracy of black hole masses is crucial for understanding the intrinsic scatter of the black hole- host galaxy scaling relations. I tested various sources of systematic uncertainty on my derived mass estimates.
The M BH estimate of the single late-type galaxy of the sample yielded an upper limit, which I could constrain very robustly. I tested the effects of dust, mass-to-light ratio (M/L) variation, and dark matter on my measured M BH . Based on these tests, the typically assumed constant M/L ratio can be an adequate assumption to account for the small amounts of dark matter in the center of that galaxy. I also tested the effect of a variable M/L variation on the M BH measurement on a second galaxy. By considering stellar M/L variations in the dynamical modeling, the measured M BH decreased by 30%. In the future, this test should be performed on additional galaxies to learn how an as constant assumed M/L flaws the estimated black hole masses.
Based on our upper limit mass measurement, I confirm previous suggestions that resolving the predicted BH sphere-of-influence is not a strict condition to measure black hole masses. Instead, it is only a rough guide for the detection of the black hole if high-quality, and high signal-to-noise IFU data are used for the measurement. About half of our sample consists of massive early-type galaxies which show nuclear surface brightness cores and signs of triaxiality. While these types of galaxies are typically modeled with axisymmetric modeling methods, the effects on M BH are not well studied yet. The massive galaxies of our presented galaxy sample are well suited to test the effect of different stellar dynamical models on the measured black hole mass in evidently triaxial galaxies. I have compared spherical Jeans and axisymmetric Schwarzschild models and will add triaxial Schwarzschild models to this comparison in the future. The constructed Jeans and Schwarzschild models mostly disagree with each other and cannot reproduce many of the triaxial features of the galaxies (e.g., nuclear sub-components, prolate rotation). The consequence of the axisymmetric-triaxial assumption on the accuracy of M BH and its impact on the black hole - host galaxy relation needs to be carefully examined in the future.
In the sample of galaxies with published M BH , we find measurements based on different dynamical tracers, requiring different observations, assumptions, and methods. Crucially, different tracers do not always give consistent results. I have used two independent tracers (cold molecular gas and stars) to estimate M BH in a regular galaxy of our sample. While the two estimates are consistent within their errors, the stellar-based measurement is twice as high as the gas-based. Similar trends have also been found in the literature. Therefore, a rigorous test of the systematics associated with the different modeling methods is required in the future. I caution to take the effects of different tracers (and methods) into account when discussing the scaling relations.
I conclude this thesis by comparing my galaxy sample with the compilation of galaxies with measured black holes from the literature, also adding six SMASHING galaxies, which were published outside of this thesis. None of the SMASHING galaxies deviates significantly from the literature measurements. Their inclusion to the published early-type galaxies causes a change towards a shallower slope for the M BH - effective velocity dispersion relation, which is mainly driven by the massive galaxies of our sample. More unbiased and homogenous measurements are needed in the future to determine the shape of the relation and understand its physical origin.
Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.
Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.
Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement.
Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement.
In this paper we report on photoswitchable polymer surfaces with dynamically and reversibly fluctuating topographies. It is well known that when azobenzene containing polymer films are irradiated with optical interference patterns the film topography changes to form a surface relief grating. In the simplest case, the film shape mimics the intensity distribution and deforms into a wave like, sinusoidal manner with amplitude that may be as large as the film thickness. This process takes place in the glassy state without photo-induced softening. Here we report on an intriguing discovery regarding the formation of reliefs under special illumination conditions. We have developed a novel setup combining the optical part for creating interference patterns, an AFM for in situ acquisition of topography changes and diffraction efficiency signal measurements. In this way we demonstrate that these gratings can be “set in motion” like water waves or dunes in the desert. We achieve this by applying repetitive polarization changes to the incoming interference pattern. Such light responsive surfaces represent the prerequisite for providing practical applications ranging from conveyer or transport systems for adsorbed liquid objects and colloidal particles to generation of adaptive and dynamic optical devices.
In this paper we report on photoswitchable polymer surfaces with dynamically and reversibly fluctuating topographies. It is well known that when azobenzene containing polymer films are irradiated with optical interference patterns the film topography changes to form a surface relief grating. In the simplest case, the film shape mimics the intensity distribution and deforms into a wave like, sinusoidal manner with amplitude that may be as large as the film thickness. This process takes place in the glassy state without photo-induced softening. Here we report on an intriguing discovery regarding the formation of reliefs under special illumination conditions. We have developed a novel setup combining the optical part for creating interference patterns, an AFM for in situ acquisition of topography changes and diffraction efficiency signal measurements. In this way we demonstrate that these gratings can be “set in motion” like water waves or dunes in the desert. We achieve this by applying repetitive polarization changes to the incoming interference pattern. Such light responsive surfaces represent the prerequisite for providing practical applications ranging from conveyer or transport systems for adsorbed liquid objects and colloidal particles to generation of adaptive and dynamic optical devices.
The goal of this thesis is to broaden the empirical basis for a better, comprehensive understanding
of massive star evolution, star formation and feedback at low metallicity. Low metallicity massive stars are a key to understand the early universe. Quantitative information on metal-poor massive stars was sparse before. The quantitative spectroscopic studies of massive star populations associated with large-scale ISM structures were not performed at low metallicity before, but are important to investigate star-formation histories and feedback in detail. Much of this work relies on spectroscopic observations with VLT-FLAMES of ~500 OB stars in the Magellanic Clouds. When available, the optical spectroscopy was complemented by UV spectra from the HST, IUE, and FUSE archives. The two representative young stellar populations that have been studied are associated with the superbubble N 206 in the Large Magellanic Cloud (LMC) and with the supergiant shell SMC-SGS 1 in the Wing of the Small Magellanic Cloud (SMC), respectively. We performed spectroscopic analyses of the massive stars using the nonLTE Potsdam Wolf-Rayet (PoWR) model atmosphere code. We estimated the stellar, wind, and feedback parameters of the individual massive stars and established their statistical distributions.
The mass-loss rates of N206 OB stars are consistent with theoretical expectations for LMC metallicity. The most massive and youngest stars show nitrogen enrichment at their surface and are found to be slower rotators than the rest of the sample. The N 206 complex has undergone star formation episodes since more than 30 Myr, with a current star formation rate higher than average in the LMC. The spatial age distribution of stars across the complex possibly indicates triggered star formation due to the expansion of the superbubble. Three very massive, young Of stars in the region dominate the ionizing and mechanical feedback among hundreds of other OB stars in the sample. The current stellar wind feedback rate from the two WR stars in the complex is comparable to that released by the whole OB sample. We see only a minor fraction of this stellar wind feedback converted into X-ray emission. In this LMC complex, stellar winds and supernovae equally contribute to the total energy feedback, which eventually powered the central superbubble. However, the total energy input accumulated over the time scale of the superbubble significantly exceeds the observed energy content of the complex. The lack of energy along with the morphology of the complex suggests a leakage of hot gas from the superbubble.
With a detailed spectroscopic study of massive stars in SMC-SGS 1, we provide the stellar and wind parameters of a large sample of OB stars at low metallicity, including those in the lower mass-range. The stellar rotation velocities show a broad, tentatively bimodal distribution, with Be stars being among the fastest. A few very luminous O stars are found close to the main sequence, while all other, slightly evolved stars obey a strict luminosity limit. Considering additional massive stars in evolved stages, with published parameters and located all over the SMC, essentially confirms this picture. The comparison with single-star evolutionary tracks suggests a dichotomy in the fate of massive stars in the SMC. Only stars with an initial mass below 30 solar masses seem to evolve from the main sequence to the cool side of the HRD to become a red supergiant and to explode as type II-P supernova. In contrast, more massive stars appear to stay always hot and might evolve quasi chemically homogeneously, finally collapsing to relatively massive black holes. However, we find no indication that chemical mixing is correlated with rapid rotation. We measured the key parameters of stellar feedback and established the links between the rates of star formation and supernovae. Our study demonstrates that in metal-poor environments stellar feedback is dominated by core-collapse supernovae in combination with winds and ionizing radiation supplied by a few of the most massive stars. We found indications of the stochastic mode of star formation, where the resulting stellar population is fully capable of producing large-scale structures such as the supergiant shell SMC-SGS 1 in the Wing. The low level of feedback in metal-poor stellar populations allows star formation episodes to persist over long timescales.
Our study showcases the importance of quantitative spectroscopy of massive stars with adequate stellar-atmosphere models in order to understand star-formation, evolution, and feedback. The stellar population analyses in the LMC and SMC make us understand that massive stars and their impact can be very different depending on their environment. Obviously, due to their different metallicity, the massive stars in the LMC and the SMC follow different evolutionary paths. Their winds differ significantly, and the key feedback agents are different. As a consequence, the star formation can proceed in different modes.
Binaries play an important role in observational and theoretical astrophysics. Since the mass and the chemical composition are key ingredients for stellar evolution, high-resolution spectroscopy is an important and necessary tool to derive those parameters to high confidence in binaries. This involves carefully measured orbital motion by the determination of radial velocity (RV) shifts and sophisticated techniques to derive the abundances of elements within the stellar atmosphere.
A technique superior to conventional cross-correlation methods to determine RV shifts in known as spectral disentangling. Hence, a major task of this thesis was the design of a sophisticated software package for this approach. In order to investigate secondary effects, such as flux and line-profile variations, imprinting changes on the spectrum the behavior of spectral disentangling on such variability is a key to understand the derived values, to improve them, and to get information about the variability itself. Therefore, the spectral disentangling code presented in this thesis and available to the community combines multiple advantages: separation of the spectra for detailed chemical analysis, derivation of orbital elements, derivation of individual RVs in order to investigate distorted systems (either by third body interaction or relativistic effects), the suppression of telluric contaminations, the derivation of variability, and the possibility to apply the technique to eclipsing binaries (important for orbital inclination) or in general to systems that undergo flux-variations.
This code in combination with the spectral synthesis codes MOOG and SME was used in order to derive the carbon 12C/13C isotope ratio (CIR) of the benchmark binary Capella. The observational result will be set into context with theoretical evolution by the use of MESA models and resolves the discrepancy of theory and observations existing since the first measurement of Capella's CIR in 1976.
The spectral disentangling code has been made available to the community and its applicability to completely different behaving systems, Wolf-Rayet stars, have also been investigated and resulted in a published article.
Additionally, since this technique relies strongly on data quality, continues development of scientific instruments to achieve best observational data is of great importance in observational astrophysics. That is the reason why there has also been effort in astronomical instrumentation during the work on this thesis.