530 Physik
Refine
Year of publication
- 2024 (11)
- 2023 (33)
- 2022 (190)
- 2021 (95)
- 2020 (183)
- 2019 (237)
- 2018 (238)
- 2017 (125)
- 2016 (36)
- 2015 (31)
- 2014 (26)
- 2013 (8)
- 2012 (18)
- 2011 (15)
- 2010 (13)
- 2009 (11)
- 2008 (7)
- 2007 (22)
- 2006 (49)
- 2005 (26)
- 2004 (30)
- 2003 (12)
- 2002 (8)
- 2001 (8)
- 2000 (4)
- 1999 (5)
- 1998 (6)
- 1997 (7)
- 1996 (7)
- 1995 (9)
- 1994 (8)
- 1989 (1)
- 1972 (1)
Document Type
- Article (886)
- Doctoral Thesis (362)
- Postprint (95)
- Other (49)
- Preprint (36)
- Habilitation Thesis (22)
- Review (12)
- Master's Thesis (10)
- Monograph/Edited Volume (4)
- Course Material (2)
Keywords
- diffusion (43)
- anomalous diffusion (31)
- gamma rays: general (20)
- synchronization (18)
- Synchronisation (16)
- organic solar cells (16)
- stochastic processes (14)
- ISM: supernova remnants (13)
- cosmic rays (13)
- data analysis (12)
Institute
- Institut für Physik und Astronomie (1480) (remove)
We apply the concepts of relative dimensions and mutual singularities to characterize the fractal properties of overlapping attractor and repeller in chaotic dynamical systems. We consider one analytically solvable example (a generalized baker's map); two other examples, the Anosov-Mobius and the Chirikov-Mobius maps, which possess fractal attractor and repeller on a two-dimensional torus, are explored numerically. We demonstrate that although for these maps the stable and unstable directions are not orthogonal to each other, the relative Renyi and Kullback-Leibler dimensions as well as the mutual singularity spectra for the attractor and repeller can be well approximated under orthogonality assumption of two fractals.
The fundamental sensitivity limit of atomic force microscopy is strongly correlated to the thermal noise of cantilever oscillation. A method to suppress this unwanted noise is to reduce the bandwidth of the measurement, but this approach is limited by the speed of the measurement and the width of the cantilever resonance, commonly defined through the quality factor Q. However, it has been shown that optomechanical resonances in interferometers might affect cantilever oscillations resulting in an effective quality factor Q(eff). When the laser power is sufficiently increased cantilever oscillations might even reach the regime of self-oscillation. In this self-oscillation state, the noise of the system is partially determined by the interaction with laser light far from equilibrium. Here, we show and discuss how tuning of laser power leads to nonlinear optomechanical effects that can dramatically increase the effective quality factor of the cantilever leading to out-of-equilibrium noise. We model the effects using a fourth order nonlinearity of the damping coefficient. Published under an exclusive license by AIP Publishing.
Data-driven expectations for electromagnetic counterpart searches based on LIGO/Virgo public alerts
(2022)
Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted.
We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall.
We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up.
We include a data release of our simulated localizations as a public proposal planning resource for astronomers.
Stochastic resetting, a diffusive process whose amplitude is reset to the origin at random times, is a vividly studied strategy to optimize encounter dynamics, e.g., in chemical reactions. Here we generalize the resetting step by introducing a random resetting amplitude such that the diffusing particle may be only partially reset towards the trajectory origin or even overshoot the origin in a resetting step. We introduce different scenarios for the random-amplitude stochastic resetting process and discuss the resulting dynamics. Direct applications are geophysical layering (stratigraphy) and population dynamics or financial markets, as well as generic search processes.
Comb-like geometric constraints leading to emergence of the time-fractional Schrödinger equation
(2021)
This paper presents an overview over several examples, where the comb-like geometric constraints lead to emergence of the time-fractional Schrodinger equation. Motion of a quantum object on a comb structure is modeled by a suitable modification of the kinetic energy operator, obtained by insertion of the Dirac delta function in the Laplacian. First, we consider motion of a free particle on two- and three-dimensional comb structures, and then we extend the study to the interacting cases. A general form of a nonlocal term, which describes the interactions of the particle with the medium, is included in the Hamiltonian, and later on, the cases of constant and Dirac delta potentials are analyzed. At the end, we discuss the case of non-integer dimensions, considering separately the case of fractal dimension between one and two, and the case of fractal dimension between two and three. All these examples show that even though we are starting with the standard time-dependent Schrodinger equation on a comb, the time-fractional equation for the Green's functions appears, due to these specific geometric constraints.
We present the discovery of a new double-detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with a P (orb) = 76.34179(2) minutes orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits, we find a low-mass sdB star, M (sdB) = 0.383 +/- 0.028 M (circle dot) with a massive white dwarf companion, M (WD) = 0.725 +/- 0.026 M (circle dot). From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of approximate to 25 Myr whereas our MESA model predicts an sdB age of approximate to 170 Myr. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion approximate to 25 Myr ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in approximate to 6 Myr and in approximate to 60 Myr the white dwarf will reach a total mass of 0.92 M (circle dot) with a thick helium layer of 0.17 M (circle dot). This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 J2238+7430 is only the second confirmed candidate for a double-detonation thermonuclear supernova. Using both systems we estimate that at least approximate to 1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.
The study addresses the question, if observed changes in terms of Arctic-midlatitude linkages during winter are driven by Arctic Sea ice decline alone or if the increase of global sea surface temperatures plays an additional role. We compare atmosphere-only model experiments with ECHAM6 to ERA-Interim Reanalysis data. The model sensitivity experiment is implemented as a set of four combinations of sea ice and sea surface temperature boundary conditions. Atmospheric circulation regimes are determined and evaluated in terms of their cyclone and blocking characteristics and changes in frequency during winter. As a prerequisite, ECHAM6 reproduces general features of circulation regimes very well. Tropospheric changes induced by the change of boundary conditions are revealed and further impacts on the large-scale circulation up into the stratosphere are investigated. In early winter, the observed increase of atmospheric blocking in the region between Scandinavia and the Urals are primarily related to the changes in sea surface temperatures. During late winter, we f nd a weakened polar stratospheric vortex in the reanalysis that further impacts the troposphere. In the model sensitivity study a climatologically weakened polar vortex occurs only if sea ice is reduced and sea surface temperatures are increased together. This response is delayed compared to the reanalysis. The tropospheric response during late winter is inconclusive in the model, which is potentially related to the weak and delayed response in the stratosphere. The model experiments do not reproduce the connection between early and late winter as interpreted from the reanalysis. Potentially explaining this mismatch, we identify a discrepancy of ECHAM6 to reproduce the weakening of the stratospheric polar vortex through blocking induced upward propagation of planetary waves.
We consider a system of noninteracting particles on a line with initial positions distributed uniformly with density ? on the negative half-line. We consider two different models: (i) Each particle performs independent Brownian motion with stochastic resetting to its initial position with rate r and (ii) each particle performs run -and-tumble motion, and with rate r its position gets reset to its initial value and simultaneously its velocity gets randomized. We study the effects of resetting on the distribution P(Q, t) of the integrated particle current Q up to time t through the origin (from left to right). We study both the annealed and the quenched current distributions and in both cases, we find that resetting induces a stationary limiting distribution of the current at long times. However, we show that the approach to the stationary state of the current distribution in the annealed and the quenched cases are drastically different for both models. In the annealed case, the whole distribution P-an(Q, t) approaches its stationary limit uniformly for all Q. In contrast, the quenched distribution P-qu(Q, t) attains its stationary form for Q < Q(crit)(t), while it remains time dependent for Q > Q(crit)(t). We show that Q(crit)(t) increases linearly with t for large t. On the scale where Q <; Q(crit)(t), we show that P-qu(Q, t) has an unusual large deviation form with a rate function that has a third-order phase transition at the critical point. We have computed the associated rate functions analytically for both models. Using an importance sampling method that allows to probe probabilities as tiny as 10-14000, we were able to compute numerically this nonanalytic rate function for the resetting Brownian dynamics and found excellent agreement with our analytical prediction.
We present a detailed spectroscopic and timing analysis of X-ray observations of the bright pulsar PSR B0656+14. The observations were obtained simultaneously with eROSITA and XMM-Newton during the calibration and performance verification phase of the Spektrum-Roentgen-Gamma mission (SRG). The analysis of the 100 ks deep observation of eROSITA is supported by archival observations of the source, including XMM-Newton, NuSTAR, and NICER. Using XMM-Newton and NICER, we first established an X-ray ephemeris for the time interval 2015 to 2020, which connects all X-ray observations in this period without cycle count alias and phase shifts. The mean eROSITA spectrum clearly reveals an absorption feature originating from the star at 570 eV with a Gaussian sigma of about 70 eV that was tentatively identified in a previous long XMM-Newton observation. A second previously discussed absorption feature occurs at 260-265 eV and is described here as an absorption edge. It could be of atmospheric or of instrumental origin. These absorption features are superposed on various emission components that are phenomenologically described here as the sum of hot (120 eV) and cold (65 eV) blackbody components, both of photospheric origin, and a power law with photon index Gamma = 2 from the magnetosphere. We created energy-dependent light curves and phase-resolved spectra with a high signal-to-noise ratio. The phase-resolved spectroscopy reveals that the Gaussian absorption line at 570 eV is clearly present throughout similar to 60% of the spin cycle, but it is otherwise undetected. Likewise, its parameters were found to be dependent on phase. The visibility of the line strength coincides in phase with the maximum flux of the hot blackbody. If the line originates from the stellar surface, it nevertheless likely originates from a different location than the hot polar cap. We also present three families of model atmospheres: a magnetized atmosphere, a condensed surface, and a mixed model. They were applied to the mean observed spectrum, whose continuum fit the observed data well. The atmosphere model, however, predicts distances that are too short. For the mixed model, the Gaussian absorption may be interpreted as proton cyclotron absorption in a field as high as 10(14) G, which is significantly higher than the field derived from the moderate observed spin-down.
Arctic climate change is marked by intensified warming compared to global trends and a significant reduction in Arctic sea ice which can intricately influence mid-latitude atmospheric circulation through tropo- and stratospheric pathways. Achieving accurate simulations of current and future climate demands a realistic representation of Arctic climate processes in numerical climate models, which remains challenging.
Model deficiencies in replicating observed Arctic climate processes often arise due to inadequacies in representing turbulent boundary layer interactions that determine the interactions between the atmosphere, sea ice, and ocean. Many current climate models rely on parameterizations developed for mid-latitude conditions to handle Arctic turbulent boundary layer processes.
This thesis focuses on modified representation of the Arctic atmospheric processes and understanding their resulting impact on large-scale mid-latitude atmospheric circulation within climate models. The improved turbulence parameterizations, recently developed based on Arctic measurements, were implemented in the global atmospheric circulation model ECHAM6. This involved modifying the stability functions over sea ice and ocean for stable stratification and changing the roughness length over sea ice for all stratification conditions. Comprehensive analyses are conducted to assess the impacts of these modifications on ECHAM6's simulations of the Arctic boundary layer, overall atmospheric circulation, and the dynamical pathways between the Arctic and mid-latitudes.
Through a step-wise implementation of the mentioned parameterizations into ECHAM6, a series of sensitivity experiments revealed that the combined impacts of the reduced roughness length and the modified stability functions are non-linear. Nevertheless, it is evident that both modifications consistently lead to a general decrease in the heat transfer coefficient, being in close agreement with the observations.
Additionally, compared to the reference observations, the ECHAM6 model falls short in accurately representing unstable and strongly stable conditions.
The less frequent occurrence of strong stability restricts the influence of the modified stability functions by reducing the affected sample size. However, when focusing solely on the specific instances of a strongly stable atmosphere, the sensible heat flux approaches near-zero values, which is in line with the observations. Models employing commonly used surface turbulence parameterizations were shown to have difficulties replicating the near-zero sensible heat flux in strongly stable stratification.
I also found that these limited changes in surface layer turbulence parameterizations have a statistically significant impact on the temperature and wind patterns across multiple pressure levels, including the stratosphere, in both the Arctic and mid-latitudes. These significant signals vary in strength, extent, and direction depending on the specific month or year, indicating a strong reliance on the background state.
Furthermore, this research investigates how the modified surface turbulence parameterizations may influence the response of both stratospheric and tropospheric circulation to Arctic sea ice loss.
The most suitable parameterizations for accurately representing Arctic boundary layer turbulence were identified from the sensitivity experiments. Subsequently, the model's response to sea ice loss is evaluated through extended ECHAM6 simulations with different prescribed sea ice conditions.
The simulation with adjusted surface turbulence parameterizations better reproduced the observed Arctic tropospheric warming in vertical extent, demonstrating improved alignment with the reanalysis data. Additionally, unlike the control experiments, this simulation successfully reproduced specific circulation patterns linked to the stratospheric pathway for Arctic-mid-latitude linkages. Specifically, an increased occurrence of the Scandinavian-Ural blocking regime (negative phase of the North Atlantic Oscillation) in early (late) winter is observed. Overall, it can be inferred that improving turbulence parameterizations at the surface layer can improve the ECHAM6's response to sea ice loss.
Perovskite semiconductors are an attractive option to overcome the limitations of established silicon based photovoltaic (PV) technologies due to their exceptional opto-electronic properties and their successful integration into multijunction cells. However, the performance of single- and multijunction cells is largely limited by significant nonradiative recombination at the perovskite/organic electron transport layer junctions. In this work, the cause of interfacial recombination at the perovskite/C-60 interface is revealed via a combination of photoluminescence, photoelectron spectroscopy, and first-principle numerical simulations. It is found that the most significant contribution to the total C-60-induced recombination loss occurs within the first monolayer of C-60, rather than in the bulk of C-60 or at the perovskite surface. The experiments show that the C-60 molecules act as deep trap states when in direct contact with the perovskite. It is further demonstrated that by reducing the surface coverage of C-60, the radiative efficiency of the bare perovskite layer can be retained. The findings of this work pave the way toward overcoming one of the most critical remaining performance losses in perovskite solar cells.
The bimolecular recombination characteristics of conjugated polymer poly[(4,4'-bis(2-ethylhexyl)dithieno[3,2-b:2',3'-d]silole)-2,6-diyl-alt-(2,5-bis 3-tetradecylthiophen-2-y1 thiazolo 5,4-d thiazole)-2,5diy1] (PDTSiTTz) blended with the fullerene series PC60BM, ICMA, ICBA, and ICTA have been investigated using microsecond and femtosecond transient absorption spectroscopy, in conjunction with electroluminescence measurements and ambient photoemission spectroscopy. The non-Langevin polymer PDTSiTTz allows an inspection of intrinsic bimolecular recombination rates uninhibited by diffusion, while the low oscillator strengths of fullerenes allow polymer features to dominate, and we compare our results to those of the well-known polymer Si-PCPDTBT. Using mu s-TAS, we have shown that the trap -limited decay dynamics of the PDTSiTTz polaron becomes progressively slower across the fullerene series, while those of Si-PCPDTBT are invariant. Electroluminescence measurements showed an unusual double peak in pristine PDTSiTTz, attributed to a low energy intragap charge transfer state, likely interchain in nature. Furthermore, while the pristine PDTSiTTz showed a broad, low-intensity density of states, the ICBA and ICTA blends presented a virtually identical DOS to Si-PCPDTBT and its blends. This has been attributed to a shift from a delocalized, interchain highest occupied molecular orbital (HOMO) in the pristine material to a dithienosilole-centered HOMO in the blends, likely a result of the bulky fullerenes increasing interchain separation. This HOMO localization had a side effect of progressively shifting the polymer HOMO to shallower energies, which was correlated with the observed decrease in bimolecular recombination rate and increased "trap" depth. However, since the density of tail states remained the same, this suggests that the traditional viewpoint of "trapping" being dominated by tail states may not encompass the full picture and that the breadth of the DOS may also have a strong influence on bimolecular recombination.
Multijunction solar cells can overcome the fundamental efficiency limits of single-junction devices. The bandgap tunability of metal halide perovskite solar cells renders them attractive for multijunction architectures(1). Combinations with silicon and copper indium gallium selenide (CIGS), as well as all-perovskite tandem cells, have been reported(2-5). Meanwhile, narrow-gap non-fullerene acceptors have unlocked skyrocketing efficiencies for organic solar cells(6,7). Organic and perovskite semiconductors are an attractive combination, sharing similar processing technologies. Currently, perovskite-organic tandems show subpar efficiencies and are limited by the low open-circuit voltage (V-oc) of wide-gap perovskite cells(8) and losses introduced by the interconnect between the subcells(9,10). Here we demonstrate perovskite-organic tandem cells with an efficiency of 24.0 per cent (certified 23.1 per cent) and a high V-oc of 2.15 volts. Optimized charge extraction layers afford perovskite subcells with an outstanding combination of high V-oc and fill factor. The organic subcells provide a high external quantum efficiency in the near-infrared and, in contrast to paradigmatic concerns about limited photostability of non-fullerene cells(11), show an outstanding operational stability if excitons are predominantly generated on the non-fullerene acceptor, which is the case in our tandems. The subcells are connected by an ultrathin (approximately 1.5 nanometres) metal-like indium oxide layer with unprecedented low optical/electrical losses. This work sets a milestone for perovskite-organic tandems, which outperform the best p-i-n perovskite single junctions(12) and are on a par with perovskite-CIGS and all-perovskite multijunctions(13).
When two initially thermal many-body systems start to interact strongly, their transient states quickly become non-Gibbsian, even if the systems eventually equilibrate. To see beyond this apparent lack of structure during the transient regime, we use a refined notion of thermality, which we call g-local. A system is g-locally thermal if the states of all its small subsystems are marginals of global thermal states. We numerically demonstrate for two harmonic lattices that whenever the total system equilibrates in the long run, each lattice remains g-locally thermal at all times, including the transient regime. This is true even when the lattices have long-range interactions within them. In all cases, we find that the equilibrium is described by the generalized Gibbs ensemble, with three-dimensional lattices requiring special treatment due to their extended set of conserved charges. We compare our findings with the well-known two-temperature model. While its standard form is not valid beyond weak coupling, we show that at strong coupling it can be partially salvaged by adopting the concept of a g-local temperature.
Core-collapse supernova remnants are structures of the interstellar medium (ISM) left behind the explosive death of most massive stars ( ?40 M-?). Since they result in the expansion of the supernova shock wave into the gaseous environment shaped by the star's wind history, their morphology constitutes an insight into the past evolution of their progenitor star. Particularly, fast-mo ving massiv e stars can produce asymmetric core-collapse superno va remnants. We inv estigate the mixing of materials in core-collapse supernova remnants generated by a moving massive 35 M-? star, in a magnetized ISM. Stellar rotation and the wind magnetic field are time-dependently included into the models which follow the entire evolution of the stellar surroundings from the zero-age main-sequence to 80 kyr after the supernova explosion. It is found that very little main-sequence material is present in remnants from moving stars, that the Wolf-Rayet wind mixes very efficiently within the 10 kyr after the explosion, while the red supergiant material is still unmixed by 30 per cent within 50 kyr after the supernova. Our results indicate that the faster the stellar motion, the more complex the internal organization of the supernova remnant and the more ef fecti ve the mixing of ejecta therein. In contrast, the mixing of stellar wind material is only weakly affected by progenitor motion, if at all.
Understanding the origin of inefficient photocurrent generation in organic solar cells with low energy offset remains key to realizing high-performance donor-acceptor systems. Here, we probe the origin of field-dependent free-charge generation and photoluminescence in wnon-fullereneacceptor (NFA)-based organic solar cells using the polymer PM6 and the NFA Y5-a non-halogenated sibling to Y6, with a smaller energetic offset to PM6. By performing time-delayed collection field (TDCF) measurements on a variety of samples with different electron transport layers and active layer thickness, we show that the fill factor and photocurrent are limited by field-dependent free charge generation in the bulk of the blend. We also introduce a new method of TDCF called m-TDCF to prove the absence of artifacts from non-geminate recombination of photogenerated and dark charge carriers near the electrodes. We then correlate free charge generation with steady-state photoluminescence intensity and find perfect anticorrelation between these two properties. Through this, we conclude that photocurrent generation in this low-offset system is entirely controlled by the field-dependent dissociation of local excitons into charge-transfer states. (c) 2023 Author(s).
Modern stationary X-ray spectroscopy is unable to resolve rotational structure.
In the present paper, we propose to use time-resolved two color X-ray pump-probe spectroscopy with picosecond resolution for real-time monitoring of the rotational dynamics induced by the recoil effect.
The proposed technique consists of two steps.
The first short pump X-ray pulse ionizes the valence electron, which transfers angular momentum to the molecule.
The second time-delayed short probe X-ray pulse resonantly excites a 1s electron to the created valence hole.
Due to the recoil-induced angular momentum the molecule rotates and changes the orientation of transition dipole moment of core-excitation with respect to the transition dipole moment of the valence ionization, which results in a temporal modulation of the probe X-ray absorption as a function of the delay time between the pulses.
We developed an accurate theory of the X-ray pump-probe spectroscopy of the recoil-induced rotation and study how the energy of the photoelectron and thermal dephasing affect the structure of the time-dependent X-ray absorption using the CO molecule as a case-study.
We also discuss the feasibility of experimental observation of our theoretical findings, opening new perspectives in studies of molecular rotational dynamics.
We characterize finite-time thermodynamic processes of multidimensional quadratic overdamped systems.
Analytic expressions are provided for heat, work, and dissipation for any evolution of the system covariance matrix.
The Bures-Wasserstein metric between covariance matrices naturally emerges as the local quantifier of dissipation.
General principles of how to apply these geometric tools to identify optimal protocols are discussed.
Focusing on the relevant slow-driving limit, we show how these results can be used to analyze cases in which the experimental control over the system is partial.
A magnetic field modifies optical properties and provides valley splitting in a molybdenum disulfide (MoS2) monolayer.
Here we demonstrate a scalable approach to the epitaxial synthesis of MoS2 monolayer on a magnetic graphene/Co system.
Using spin- and angle-resolved photoemission spectroscopy we observe a magnetic proximity effect that causes a 20 meV spin-splitting at the (Gamma) over bar point and canting of spins at the (K) over bar point in the valence band toward the in-plane direction of cobalt magnetization.
Our density functional theory calculations reveal that the in-plane spin component at (K) over bar is localized on Co atoms in the valence band, while in the conduction band it is localized on the MoS2 layer.
The calculations also predict a 16 meV spin-splitting at the (Gamma) over bar point and 8 meV (K) over bar-(K) over bar' valley asymmetry for an out-of-plane magnetization. These findings suggest control over optical transitions in MoS2 via Co magnetization. Our estimations show that the magnetic proximity effect is equivalent to the action of the magnetic field as large as 100 T.
Auger-photoelectron coincidence spectroscopy (APECS) has been used to examine the electron correlation and itinerance effects in transition metals Cu, Ni and Co.
It is shown that the LVV Auger, in coincidence with 2p photoelectrons, spectra can be represented using atomic multiplet positions if the 3d-shell is localized (atomic-like) and with a self-convoluted valence band for band-like (itinerant) materials as explained using the Cini-Sawatzky model.
For transition metals, the 3d band changes from band-like to localized with increasing atomic number, with the possibility of a mixed behavior.
Our result shows that the LVV spectra of Cu can be represented by atomic multiplet calculations, those of Co resemble the self-convolution of the valence band and those of Ni are a mixture of both, consistent with the Cini-Sawatzky model.
Diffusion with stochastic resetting is a paradigm of resetting processes. Standard renewal or master equation approach are typically used to study steady state and other transport properties such as average, mean squared displacement etc.
What remains less explored is the two time point correlation functions whose evaluation is often daunting since it requires the implementation of the exact time dependent probability density functions of the resetting processes which are unknown for most of the problems.
We adopt a different approach that allows us to write a stochastic solution for a single trajectory undergoing resetting.
Moments and the autocorrelation functions between any two times along the trajectory can then be computed directly using the laws of total expectation. Estimation of autocorrelation functions turns out to be pivotal for investigating the ergodic properties of various observables for this canonical model.
In particular, we investigate two observables (i) sample mean which is widely used in economics and (ii) time-averaged-mean-squared-displacement (TAMSD) which is of acute interest in physics.
We find that both diffusion and drift-diffusion processes with resetting are ergodic at the mean level unlike their reset-free counterparts. In contrast, resetting renders ergodicity breaking in the TAMSD while both the stochastic processes are ergodic when resetting is absent. We quantify these behaviors with detailed analytical study and corroborate with extensive numerical simulations.
Our results can be verified in experimental set-ups that can track single particle trajectories and thus have strong implications in understanding the physics of resetting.
With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to 1250 m (2% at 90% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudospectral code ELLIPTICA, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, ELLIPTICA deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of ELLIPTICA and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude ∼0.8.
We perform numerical studies of a thermally driven, overdamped particle in a random quenched force field, known as the Sinai model. We compare the unbounded motion on an infinite 1-dimensional domain to the motion in bounded domains with reflecting boundaries and show that the unbounded motion is at every time close to the equilibrium state of a finite system of growing size. This is due to time scale separation: inside wells of the random potential, there is relatively fast equilibration, while the motion across major potential barriers is ultraslow. Quantities studied by us are the time dependent mean squared displacement, the time dependent mean energy of an ensemble of particles, and the time dependent entropy of the probability distribution. Using a very fast numerical algorithm, we can explore times up top 10(17) steps and thereby also study finite-time crossover phenomena.
Improving permafrost dynamics in land surface models: insights from dual sensitivity experiments
(2024)
The thawing of permafrost and the subsequent release of greenhouse gases constitute one of the most significant and uncertain positive feedback loops in the context of climate change, making predictions regarding changes in permafrost coverage of paramount importance. To address these critical questions, climate scientists have developed Land Surface Models (LSMs) that encompass a multitude of physical soil processes. This thesis is committed to advancing our understanding and refining precise representations of permafrost dynamics within LSMs, with a specific focus on the accurate modeling of heat fluxes, an essential component for simulating permafrost physics.
The first research question overviews fundamental model prerequisites for the representation of permafrost soils within land surface modeling. It includes a first-of-its-kind comparison between LSMs in CMIP6 to reveal their differences and shortcomings in key permafrost physics parameters. Overall, each of these LSMs represents a unique approach to simulating soil processes and their interactions with the climate system. Choosing the most appropriate model for a particular application depends on factors such as the spatial and temporal scale of the simulation, the specific research question, and available computational resources.
The second research question evaluates the performance of the state-of-the-art Community Land Model (CLM5) in simulating Arctic permafrost regions. Our approach overcomes traditional evaluation limitations by individually addressing depth, seasonality, and regional variations, providing a comprehensive assessment of permafrost and soil temperature dynamics. I compare CLM5's results with three extensive datasets: (1) soil temperatures from 295 borehole stations, (2) active layer thickness (ALT) data from the Circumpolar Active Layer Monitoring Network (CALM), and (3) soil temperatures, ALT, and permafrost extent from the ESA Climate Change Initiative (ESA-CCI). The results show that CLM5 aligns well with ESA-CCI and CALM for permafrost extent and ALT but reveals a significant global cold temperature bias, notably over Siberia. These results echo a persistent challenge identified in numerous studies: the existence of a systematic 'cold bias' in soil temperature over permafrost regions. To address this challenge, the following research questions propose dual sensitivity experiments.
The third research question represents the first study to apply a Plant Functional Type (PFT)-based approach to derive soil texture and soil organic matter (SOM), departing from the conventional use of coarse-resolution global data in LSMs. This novel method results in a more uniform distribution of soil organic matter density (OMD) across the domain, characterized by reduced OMD values in most regions. However, changes in soil texture exhibit a more intricate spatial pattern. Comparing the results to observations reveals a significant reduction in the cold bias observed in the control run. This method shows noticeable improvements in permafrost extent, but at the cost of an overestimation in ALT. These findings emphasize the model's high sensitivity to variations in soil texture and SOM content, highlighting the crucial role of soil composition in governing heat transfer processes and shaping the seasonal variation of soil temperatures in permafrost regions.
Expanding upon a site experiment conducted in Trail Valley Creek by \citet{dutch_impact_2022}, the fourth research question extends the application of the snow scheme proposed by \citet{sturm_thermal_1997} to cover the entire Arctic domain. By employing a snow scheme better suited to the snow density profile observed over permafrost regions, this thesis seeks to assess its influence on simulated soil temperatures. Comparing this method to observational datasets reveals a significant reduction in the cold bias that was present in the control run. In most regions, the Sturm run exhibits a substantial decrease in the cold bias. However, there is a distinctive overshoot with a warm bias observed in mountainous areas. The Sturm experiment effectively addressed the overestimation of permafrost extent in the control run, albeit resulting in a substantial reduction in permafrost extent over mountainous areas. ALT results remain relatively consistent compared to the control run. These outcomes align with our initial hypothesis, which anticipated that the reduced snow insulation in the Sturm run would lead to higher winter soil temperatures and a more accurate representation of permafrost physics.
In summary, this thesis demonstrates significant advancements in understanding permafrost dynamics and its integration into LSMs. It has meticulously unraveled the intricacies involved in the interplay between heat transfer, soil properties, and snow dynamics in permafrost regions. These insights offer novel perspectives on model representation and performance.
Non-fullerene acceptors (NFAs) as used in state-of-the-art organic solar cells feature highly crystalline layers that go along with low energetic disorder.
Here, the crucial role of energetic disorder in blends of the donor polymer PM6 with two Y-series NFAs, Y6, and N4 is studied.
By performing temperature-dependent charge transport and recombination studies, a consistent picture of the shape of the density of state distributions for free charges in the two blends is developed, allowing an analytical description of the dependence of the open-circuit voltage V-OC on temperature and illumination intensity.
Disorder is found to influence the value of the V-OC at room temperature, but also its progression with temperature. Here, the PM6:Y6 blend benefits substantially from its narrower state distributions.
The analysis also shows that the energy of the equilibrated free charge population is well below the energy of the NFA singlet excitons for both blends and possibly below the energy of the populated charge transfer manifold, indicating a down-hill driving force for free charge formation.
It is concluded that energetic disorder of charge-separated states has to be considered in the analysis of the photovoltaic properties, even for the more ordered PM6:Y6 blend.
The formation of large polarons has been proposed as reason for the high defect tolerance, low mobility, low charge carrier trapping, and low nonradiative recombination rates of lead halide perovskites. Recently, direct evidence for large-polaron formation has been reported from a 50% effective mass enhancement in angle-resolved photoemission of CsPbBr3 over theory for the orthorhombic structure. We present in-depth band dispersion measurements of CsPbBr3 and GW calculations, which lead to similar effective masses at the valence band maximum of 0.203 1 0.016 m0 in experiment and 0.226 m0 in orthorhombic theory. We argue that the effective mass can be explained solely on the basis of electron-electron correlation and largepolaron formation cannot be concluded from photoemission data.
Only a fast and global transformation towards decarbonization and sustainability can keep the Earth in a civilization-friendly state. As hotspots for (green) innovation and experimentation, cities could play an important role in this transition. They are also known to profit from each other's ideas, with policy and technology innovations spreading to other cities. In this way, cities can be conceptualized as nodes in a globe-spanning learning network. The dynamics of this process are important for society's response to climate change and other challenges, but remain poorly understood on a macroscopic level. In this contribution, we develop an approach to identify whether network-based complex contagion effects are a feature of sustainability policy adoption by cities, based on dose-response contagion and surrogate data models. We apply this methodology to an exemplary data set, comprising empirical data on the spreading of a public transport innovation (Bus Rapid Transit Systems) and a global inter-city connection network based on scheduled flight routes. Although our approach is not able to identify detailed mechanisms, our results point towards a contagious spreading process, and cannot be explained by either the network structure or the increase in global adoption rate alone. Further research on the role of a city's abstract "global neighborhood" regarding its policy and innovation decisions is thus both needed and promising, and may connect with research on social tipping processes. The methodology is generic, and can be used to compare the predictive power for innovation spreading of different kinds of inter-city network connections, e.g. via transport links, trade, or co-membership in political networks.
We investigate the effects of Markovian resetting events on continuous time random walks where the waiting times and the jump lengths are random variables distributed according to power-law probability density functions.
We prove the existence of a nonequilibrium stationary state and finite mean first arrival time.
However, the existence of an optimum reset rate is conditioned to a specific relationship between the exponents of both power-law tails.
We also investigate the search efficiency by finding the optimal random walk which minimizes the mean first arrival time in terms of the reset rate, the distance of the initial position to the target, and the characteristic transport exponents.
Leptonic nonthermal emission from supernova remnants evolving in the circumstellar magnetic field
(2022)
The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their nonthermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting nonthermal emission.
Ultrafast X-ray diffraction is used to quantify the transport of energy in laser-excited nanoscale gold-nickel (Au-Ni) bilayers.
Electron transport and efficient electron-phonon coupling in Ni convert the laser-deposited energy in the conduction electrons within a few picoseconds into a strong non-equilibrium between hot Ni and cold Au phonons at the bilayer interface.
Modeling of the subsequent equilibration dynamics within various two-temperature models confirms that for ultrathin Au films, the thermal transport is dominated by phonons instead of conduction electrons because of the weak electron-phonon coupling in Au.
Magnetic reconnection is a multi-faceted process of energy conversion in astrophysical, space and laboratory plasmas that operates at microscopic scales but has macroscopic drivers and consequences.
Solar flares present a key laboratory for its study, leaving imprints of the microscopic physics in radiation spectra and allowing the macroscopic evolution to be imaged, yet a full observational characterization remains elusive.
Here we combine high resolution imaging and spectral observations of a confined solar flare at multiple wavelengths with data-constrained magnetohydrodynamic modeling to study the dynamics of the flare plasma from the current sheet to the plasmoid scale. The analysis suggests that the flare resulted from the interaction of a twisted magnetic flux rope surrounding a filament with nearby magnetic loops whose feet are anchored in chromospheric fibrils. Bright cusp-shaped structures represent the region around a reconnecting separator or quasi-separator (hyperbolic flux tube).
The fast reconnection, which is relevant for other astrophysical environments, revealed plasmoids in the current sheet and separatrices and associated unresolved turbulent motions.
Solar flares provide wide range of observational details about fundamental processes involved. Here, the authors show evidence for magnetic reconnection in a strong confined solar flare displaying all four reconnection flows with plasmoids in the current sheet and the separatrices.
The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant H-0 that does not rely on a cosmic distance ladder, nor assumes a specific cosmological model.
By using gravitational waves as "standard sirens", this approach holds promise to arbitrate the existing tension between the H-0 value inferred from the cosmic microwave background and those obtained from local measurements.
However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source led to a H-0 value from GW170817 that was not precise enough to resolve the existing tension.
In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on H-0.
We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.
In the last years, electron density profile functions characterized by a linear dependence on the scale height showed good results when approximating the topside ionosphere. The performance above 800 km, however, is not yet well investigated.
This study investigates the capability of the semi-Epstein functions to represent electron density profiles from the peak height up to 20,000 km. Electron density observations recorded by the Van Allen Probes were used to resolve the scale height dependence in the plasmasphere.
It was found that the linear dependence of the scale height in the topside ionosphere cannot be directly used to extrapolate profiles above 800 km.
We find that the dependence of scale heights on altitude is quadratic in the plasmasphere. A statistical model of the scale heights is therefore proposed. After combining the topside ionosphere and plasmasphere by a unified model, we have obtained good estimations not only in the profile shapes, but also in the Total Electron Content magnitude and distributions when compared to actual measurements from 2013, 2014, 2016 and 2017.
Our investigation shows that Van Allen Probes can be merged to radio-occultation data to properly represent the upper ionosphere and plasmasphere by means of a semi-Epstein function.
A large Rashba effect is essential for future applications in spintronics. Particularly attractive is understanding and controlling nonequilibrium properties of ferroelectric Rashba semiconductors. Here, time- and angle-resolved photoemission is utilized to access the ultrafast dynamics of bulk and surface transient Rashba states after femtosecond optical excitation of GeTe. A complex thermalization pathway is observed, wherein three different timescales can be clearly distinguished: intraband thermalization, interband equilibration, and electronic cooling. These dynamics exhibit an unconventional temperature dependence: while the cooling phase speeds up with increasing sample temperature, the opposite happens for interband thermalization. It is demonstrated how, due to the Rashba effect, an interdependence of these timescales on the relative strength of both electron-electron and electron-phonon interactions is responsible for the counterintuitive temperature dependence, with spin-selection constrained interband electron-electron scatterings found both to dominate dynamics away from the Fermi level, and to weaken with increasing temperature. These findings are supported by theoretical calculations within the Boltzmann approach explicitly showing the opposite behavior of all relevant electron-electron and electron-phonon scattering channels with temperature, thus confirming the microscopic mechanism of the experimental findings. The present results are important for future applications of ferroelectric Rashba semiconductors and their excitations in ultrafast spintronics.
Using over-5-year EMFISIS wave measurements from Van Allen Probes, we present a detailed survey to identify the controlling factors among the geomagnetic indices and solar wind parameters for the 1-min root mean square amplitudes of lower band chorus (LBC) and upper band chorus (UBC).
A set of important features are automatically determined by feature selection techniques, namely, Random Forest and Maximum Relevancy Minimum Redundancy. Our analysis results indicate the AE index with zero-time-delay dominates the intensity evolution of LBC and UBC, consistent with the evidence that chorus waves prefer to occur and amplify during enhanced substorm periods. Regarding solar wind parameters, solar wind speed and IMF B-z are identified as the controlling factors for chorus wave intensity. Using the combination of all these important features, a predictive neural network model of chorus wave intensity is established to reconstruct the temporal variations of chorus wave intensity, for which application of Random Forest produces the overall best performance.
Plain Language Summary
Whistler mode chorus waves are electromagnetic waves observed in the low-density region near the geomagnetic equator outside the plasmapause. The dynamics of Earth's radiation belts are largely influenced by chorus waves owing to their dual contributions to both radiation belt electron acceleration and loss. In this study, we use feature selection techniques to identify the controlling geomagnetic and solar wind factors for magnetospheric chorus waves. Feature selection techniques implement the processes which can select the features most influential to the output.
In this study, the inputs are geomagnetic indices and solar wind parameters and the output is the chorus wave intensity. The results indicate that AE index with zerotime delay dominates the chorus wave intensity. Furthermore, solar wind speed and IMF B-z are identified as the most important solar wind drivers for chorus wave intensity.
On basis of the combination of all these important geomagnetic and solar wind controlling factors, we develop a neural network model of chorus wave intensity, and find that the model with the inputs identified using the Random Forest method produces the overall best performance.
The determination of the spin state of iron-bearing compounds at high pressure and temperature is crucial for our understanding of chemical and physical properties of the deep Earth. Studies on the relationship between the coordination of iron and its electronic spin structure in iron-bearing oxides, silicates, carbonates, iron alloys, and other minerals found in the Earth's mantle and core are scarce because of the technical challenges to simultaneously probe the sample at high pressures and temperatures. We used the unique properties of a pulsed and highly brilliant x-ray free electron laser (XFEL) beam at the High Energy Density (HED) instrument of the European XFEL to x-ray heat and probe samples contained in a diamond anvil cell. We heated and probed with the same x-ray pulse train and simultaneously measured x-ray emission and x-ray diffraction of an FeCO3 sample at a pressure of 51 GPa with up to melting temperatures. We collected spin state sensitive Fe K beta(1,3) fluorescence spectra and detected the sample's structural changes via diffraction, observing the inverse volume collapse across the spin transition. During x-ray heating, the carbonate transforms into orthorhombic Fe4C3O12 and iron oxides. Incipient melting was also observed. This approach to collect information about the electronic state and structural changes from samples contained in a diamond anvil cell at melting temperatures and above will considerably improve our understanding of the structure and dynamics of planetary and exoplanetary interiors.
We present a study of the control of electric field induced strain on the magnetic and electrical transport properties in a magnetoelastically coupled artificial multiferroic Fe3O4/BaTiO3 heterostructure.
In this Fe3O4/BaTiO3 heterostructure, the Fe3O4 thin film is epitaxially grown in the form of bilateral domains, analogous to a-c stripe domains of the underlying BaTiO3(001) substrate.
By in situ electric field dependent magnetization measurements, we demonstrate the extrinsic control of the magnetic anisotropy and the characteristic Verwey metal-insulator transition of the epitaxial Fe3O4 thin film in a wide temperature range between 20-300 K, via strain mediated converse magnetoelectric coupling.
In addition, we observe strain induced modulations in the magnetic and electrical transport properties of the Fe3O4 thin film across the thermally driven intrinsic ferroelectric and structural phase transitions of the BaTiO3 substrate.
In situ electric field dependent Raman measurements reveal that the electric field does not significantly modify the antiphase boundary defects in the Fe3O4 thin film once it is thermodynamically stable after deposition and that the modification of the magnetic properties is mainly caused by strain induced lattice distortions and magnetic anisotropy.
These results provide a framework to realize electrical control of the magnetization in a classical highly correlated transition metal oxide.
The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search.
How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM.
We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed.
We consider the spatiotemporal states of an ensemble of nonlocally coupled nonidentical phase oscillators, which correspond to different regimes of the long-term evolution of such a system. We have obtained homogeneous, twisted, and nonhomogeneous stationary solutions to the Ott-Antonsen equations corresponding to key variants of the realized collective rotational motion of elements of the medium in question with nonzero mesoscopic characteristics determining the degree of coherence of the dynamics of neighboring particles. We have described the procedures of the search for the class of nonhomogeneous solutions as stationary points of the auxiliary point map and of determining the stability based on analysis of the eigenvalue spectrum of the composite operator. Static and breather cluster regimes have been demonstrated and described, as well as the regimes with an irregular behavior of averaged complex fields including, in particular, the local order parameter.
The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours.
Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process.
For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations.
Lennard-Jones mixtures represent one of the popular systems for the study of glass-forming liquids.
Spatio/temporal heterogeneity and rare (activated) events are at the heart of the slow dynamics typical of these systems. Such slow dynamics is characterised by the development of a plateau in the mean-squared displacement (MSD) at intermediate times, accompanied by a non-Gaussianity in the displacement distribution identified by exponential tails.
As pointed out by some recent works, the non-Gaussianity persists at times beyond the MSD plateau, leading to a Brownian yet non-Gaussian regime and thus highlighting once again the relevance of rare events in such systems.
Single-particle motion of glass-forming liquids is usually interpreted as an alternation of rattling within the local cage and cage-escape motion and therefore can be described as a sequence of waiting times and jumps. In this work, by using a simple yet robust algorithm, we extract jumps and waiting times from single-particle trajectories obtained via molecular dynamics simulations.
We investigate the presence of correlations between waiting times and find negative correlations, which becomes more and more pronounced when lowering the temperature.
Organic solar cells (OSCs) have progressed rapidly in recent years through the development of novel organic photoactive materials, especially non-fullerene acceptors (NFAs). Consequently, OSCs based on state-of-the-art NFAs have reached significant milestones, such as similar to 19% power conversion efficiencies (PCEs) and small energy losses (less than 0.5 eV). Despite these significant advances, understanding of the interplay between molecular structure and optoelectronic properties lags significantly behind. For example, despite the theoretical framework for describing the energetic disorder being well developed for the case of inorganic semiconductors, the question of the applicability of classical semiconductor theories in analyzing organic semiconductors is still under debate. A general observation in the inorganic field is that inorganic photovoltaic materials possessing a polycrystalline microstructure exhibit suppressed disorder properties and better charge carrier transport compared to their amorphous analogs. Accordingly, this principle extends to the organic semiconductor field as many organic photovoltaic materials are synthesized to pursue polycrystalline-like features. Yet, there appears to be sporadic examples that exhibit an opposite trend. However, full studies decoupling energetic disorder from aggregation effects have largely been left out. Hence, the potential role of the energetic disorder in OSCs has received little attention. Interestingly, recently reported state-of-the-art NFA-based devices could achieve a small energetic disorder and high PCE at the same time; and interest in this investigation related to the disorder properties in OSCs was revived. In this contribution, progress in terms of the correlation between molecular design and energetic disorder is reviewed together with their effects on the optoelectronic mechanism and photovoltaic performance. Finally, the specific challenges and possible solutions in reducing the energetic disorder of OSCs from the viewpoint of materials and devices are proposed.
In organic solar cells, the resulting device efficiency depends strongly on the local morphology and intermolecular interactions of the blend film. Optical spectroscopy was used to identify the spectral signatures of interacting chromophores in blend films of the donor polymer PM6 with two state-of-the-art nonfullerene acceptors, Y6 and N4, which differ merely in the branching point of the side chain. From temperature-dependent absorption and luminescence spectroscopy in solution, it is inferred that both acceptor materials form two types of aggregates that differ in their interaction energy. Y6 forms an aggregate with a predominant J-type character in solution, while for N4 molecules the interaction is predominantly in a H-like manner in solution and freshly spin-cast film, yet the molecules reorient with respect to each other with time or thermal annealing to adopt a more J-type interaction. The different aggregation behavior of the acceptor materials is also reflected in the blend films and accounts for the different solar cell efficiencies reported with the two blends.
The mobile-immobile model (MIM) has been established in geoscience in the context of contaminant transport in groundwater. Here the tracer particles effectively immobilise, e.g., due to diffusion into dead-end pores or sorption. The main idea of the MIM is to split the total particle density into a mobile and an immobile density. Individual tracers switch between the mobile and immobile state following a two-state telegraph process, i.e., the residence times in each state are distributed exponentially. In geoscience the focus lies on the breakthrough curve (BTC), which is the concentration at a fixed location over time. We apply the MIM to biological experiments with a special focus on anomalous scaling regimes of the mean squared displacement (MSD) and non-Gaussian displacement distributions. As an exemplary system, we have analysed the motion of tau proteins, that diffuse freely inside axons of neurons. Their free diffusion thereby corresponds to the mobile state of the MIM. Tau proteins stochastically bind to microtubules, which effectively immobilises the tau proteins until they unbind and continue diffusing. Long immobilisation durations compared to the mobile durations give rise to distinct non-Gaussian Laplace shaped distributions. It is accompanied by a plateau in the MSD for initially mobile tracer particles at relevant intermediate timescales. An equilibrium fraction of initially mobile tracers gives rise to non-Gaussian displacements at intermediate timescales, while the MSD remains linear at all times. In another setting bio molecules diffuse in a biosensor and transiently bind to specific receptors, where advection becomes relevant in the mobile state. The plateau in the MSD observed for the advection-free setting and long immobilisation durations persists also for the case with advection. We find a new clear regime of anomalous diffusion with non-Gaussian distributions and a cubic scaling of the MSD. This regime emerges for initially mobile and for initially immobile tracers. For an equilibrium fraction of initially mobile tracers we observe an intermittent ballistic scaling of the MSD. The long-time effective diffusion coefficient is enhanced by advection, which we physically explain with the variance of mobile durations. Finally, we generalize the MIM to incorporate arbitrary immobilisation time distributions and focus on a Mittag-Leffler immobilisation time distribution with power-law tail ~ t^(-1-mu) with 0<mu<1 and diverging mean immobilisation durations. A fit of our model to the BTC of experimental data from tracer particles in aquifers matches the BTC including the power-law tail. We use the fit parameters for plotting the displacement distributions and the MSD. We find Gaussian normal diffusion at short times and long-time power-law decay of mobile mass accompanied by anomalous diffusion at long times. The long-time diffusion is subdiffusive in the advection-free setting, while it is either subdiffusive for 0<mu<1/2 or superdiffusive for 1/2<mu<1 when advection is present. In the long-time limit we show equivalence of our model to a bi-fractional diffusion equation.
We investigate a class of diffusion-controlled reactions that are initiated at the time instance when a prescribed number K among N particles independently diffusing in a solvent are simultaneously bound to a target region.
In the irreversible target-binding setting, the particles that bind to the target stay there forever, and the reaction time is the Kth fastest first-passage time to the target, whose distribution is well-known. In turn, reversible binding, which is common for most applications, renders theoretical analysis much more challenging and drastically changes the distribution of reaction times.
We develop a renewal-based approach to derive an approximate solution for the probability density of the reaction time.
This approximation turns out to be remarkably accurate for a broad range of parameters.
We also analyze the dependence of the mean reaction time or, equivalently, the inverse reaction rate, on the main parameters such as K, N, and binding/unbinding constants. Some biophysical applications and further perspectives are briefly discussed.
We consider an array of nonlocally coupled oscillators on a ring, which for equally spaced units possesses a Kuramoto-Battogtokh chimera regime and a synchronous state. We demonstrate that disorder in oscillators positions leads to a transition from the synchronous to the chimera state. For a static (quenched) disorder we find that the probability of synchrony survival depends on the number of particles, from nearly zero at small populations to one in the thermodynamic limit. Furthermore, we demonstrate how the synchrony gets destroyed for randomly (ballistically or diffusively) moving oscillators. We show that, depending on the number of oscillators, there are different scalings of the transition time with this number and the velocity of the units.
We report generation of ultra-broadband longitudinal acoustic coherent phonon wavepackets in SrTiO3 (STO) with frequency components extending throughout the first Brillouin zone. The wavepackets are efficiently generated in STO using femtosecond infrared laser excitation of an atomically flat 1.6 nm-thick epitaxial SrRuO3 film. We use femtosecond x-ray diffraction at the European X-Ray Free Electron Laser Facility to study the dispersion and damping of phonon wavepackets. The experimentally determined damping constants for multi-THz frequency phonons compare favorably to the extrapolation of a simple ultrasound damping model over several orders of magnitude.