Gold Open-Access
Refine
Document Type
- Article (197)
- Other (2)
- Review (2)
- Monograph/Edited Volume (1)
Is part of the Bibliography
- yes (202)
Keywords
- diffusion (13)
- Reflexion (4)
- chemotaxis (4)
- anomalous diffusion (3)
- electrons (3)
- magnetosphere (3)
- model (3)
- perovskite solar cells (3)
- synchronization (3)
- Atlantic meridional overturning circulation (2)
Institute
- Institut für Physik und Astronomie (202) (remove)
Using time-resolved x-ray diffraction, we demonstrate the manipulation of the picosecond strain response of a metallic heterostructure consisting of a dysprosium (Dy) transducer and a niobium (Nb) detection layer by an external magnetic field. We utilize the first-order ferromagnetic–antiferromagnetic phase transition of the Dy layer, which provides an additional large contractive stress upon laser excitation compared to its zerofield response. This enhances the laser-induced contraction of the transducer and changes the shape of the picosecond strain pulses driven in Dy and detected within the buried Nb layer. Based on our experiment with rare-earth metals we discuss required properties for functional transducers, which may allow for novel field-control of the emitted picosecond strain pulses.
Reflexion ist eine Schlüsselkategorie für die professionelle Entwicklung von Lehrkräften, welche als Ausbildungsziel in den Bildungsstandards für die Lehrkräftebildung verankert ist. Eine Verstetigung universitär geprägter Forschung und Modellierung in der praxisnahen Anwendung im schulischen Kontext bietet Potentiale nachhaltiger Professionalisierung. Die Stärkung reflexionsbezogener Kompetenzen durch Empirie und Anwendung scheint eine phasenübergreifende Herausforderung der Lehrkräftebildung zu sein, die es zu bewältigen gilt. Ziele des Tagungsbandes Reflexion in der Lehrkräftebildung sind eine theoretische Schärfung des Konzeptes „Reflexive Professionalisierung“ und der Austausch über Fragen der Einbettung wirksamer reflexionsbezogener Lerngelegenheiten in die Lehrkräftebildung. Forschende und Lehrende der‚ drei Phasen (Studium, Referendariat sowie Fort- und Weiterbildung) der Lehrkräftebildung stellen Lehrkonzepte und Forschungsprojekte zum Thema Reflexion in der Lehrkräftebildung vor und diskutieren diese. Gemeinsam mit Teilnehmenden aller Phasen und von verschiedenen Standorten der Lehrkräftebildung werden zukünftige Herausforderungen identifiziert und Lösungsansätze herausgearbeitet.
Schulpraktische Phasen stellen eine bedeutende praxisnahe Lerngelegenheit im Lehramtsstudium dar, da sie Raum für umfangreiche Reflexionen der eigenen Lernerfahrung bieten. Das im Studium erworbene theoretisch-formale Wissen steht hierbei dem praktischen Wissen und Können gegenüber. Mit der professionellen Entwicklung im Referendariat, besonders im Kompetenzbereich des Unterrichtens, kann geschlussfolgert werden, dass sich eine Reflexion über eher fachliche Aspekte unter den Studierenden im Referendariat auf eine Reflexion über eher überfachliche und pädagogische Aspekte weitet. Infolge der Analyse von N = 55 schriftlichen Fremdreflexionen von angehenden Physiklehrkräften aus Studium und Referendariat konnte diese Hypothese für den Bereich der Unterrichtsanalyse und -reflexion unterstützt werden. Weiter wurde aus der Videovignette ein Workshopangebot für Lehrkräfte der zweiten und dritten Phase der Lehrkräftebildung entwickelt, erprobt und evaluiert.
Für die Entwicklung professioneller Handlungskompetenzen angehender Lehrkräfte stellt die Unterrichtsreflexion ein wichtiges Instrument dar, um Theoriewissen und Praxiserfahrungen in Beziehung zu setzen. Die Auswertung von Unterrichtsreflexionen und eine entsprechende Rückmeldung stellt Forschende und Dozierende allerdings vor praktische wie theoretische Herausforderungen. Im Kontext der Forschung zu Künstlicher Intelligenz (KI) entwickelte Methoden bieten hier neue Potenziale. Der Beitrag stellt überblicksartig zwei Teilstudien vor, die mit Hilfe von KI-Methoden wie dem maschinellen Lernen untersuchen, inwieweit eine Auswertung von Unterrichtsreflexionen angehender Physiklehrkräfte auf Basis eines theoretisch abgeleiteten Reflexionsmodells und die automatisierte Rückmeldung hierzu möglich sind. Dabei wurden unterschiedliche Ansätze des maschinellen Lernens verwendet, um modellbasierte Klassifikation und Exploration von Themen in Unterrichtsreflexionen umzusetzen. Die Genauigkeit der Ergebnisse wurde vor allem durch sog. Große Sprachmodelle gesteigert, die auch den Transfer auf andere Standorte und Fächer ermöglichen. Für die fachdidaktische Forschung bedeuten sie jedoch wiederum neue Herausforderungen, wie etwa systematische Verzerrungen und Intransparenz von Entscheidungen. Dennoch empfehlen wir, die Potenziale der KI-basierten Methoden gründlicher zu erforschen und konsequent in der Praxis (etwa in Form von Webanwendungen) zu implementieren.
Reflexion – unhinterfragt eines der wichtigsten Worte im Kontext der Lehrkräftebildung. Fest verankert in den bundesdeutschen Bildungsstandards sind in Forschung und Lehre die Suche nach Evidenz und die Unterstützung (angehender) Lehrkräfte ständiger Antrieb unzähliger Akteur:innen aller Phasen der Lehrkräftebildung. Wenngleich begriff liche Unklarheiten die Kommunikation von Forschungsergebnissen nicht immer intuitiv und die Unterstützung in der Lehre nicht immer praktikabel werden lassen, besteht Einigkeit darüber, dass ein Diskurs zur reflexiven Professionalisierung von Lehrkräften geführt werden muss. Aus diesem Grund veranstalteten die beiden QLB-Projekte PSI-Potsdam der Universität Potsdam und K2teach der Freien Universität Berlin vom 5. bis 7. Oktober 2022 die Onlinetagung „Reflexion in der Lehrkräftebildung. Empirisch – Phasenübergreifend – Interdisziplinär“. Ausgehend von den verschiedensten Fachdisziplinen diskutierten Akteur:innen aller Phasen der Lehrkräftebildung unterschiedlicher Standorte Ergebnisse empirischer Studien und Erfahrungen aus der Arbeit mit (angehenden) Lehrkräften. Beiträge der Tagung sind in diesem Buch festgehalten und sind als Momentaufnahme eines sich ständig entwickelnden Themenfelds zu verstehen. Forschende und Lehrende haben mit dieser Momentaufnahme die Möglichkeit, Eindrücke für die eigene Arbeit aufzunehmen und weiterzuentwickeln.
Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large-scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which comove with the upstream flow and as nonlinear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex nonlinear structures bound by a strong magnetic field around four times the upstream value. Although these nonlinear structures do not in general interact with cospatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around 0.8% of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately three times more efficient than, stochastic shock drift acceleration.
In recent years, the search for more efficient and environmentally friendly materials to be employed in the next generation of thin film solar cell devices has seen a shift towards hybrid halide perovskites and chalcogenide materials crystallising in the kesterite crystal structure. Prime examples for the latter are Cu2ZnSnS4, Cu2ZnSnSe4, and their solid solution Cu2ZnSn(SxSe1-x)(4), where actual devices already demonstrated power conversion efficiencies of about 13 %. However, in their naturally occurring kesterite crystal structure, the so-called Cu-Zn disorder plays an important role and impacts the structural, electronic, and optical properties. To understand the influence of Cu-Zn disorder, we perform first-principles calculations based on density functional theory combined with special quasirandom structures to accurately model the cation disorder. Since the electronic band gaps and derived optical properties are severely underestimated by (semi)local exchange and correlation functionals, supplementary hybrid functional calculations have been performed. Concerning the latter, we additionally employ a recently devised technique to speed up structural relaxations for hybrid functional calculations. Our calculations show that the Cu-Zn disorder leads to a slight increase in the unit cell volume compared to the conventional kesterite structure showing full cation order, and that the band gap gets reduced by about 0.2 eV, which is in very good agreement with earlier experimental and theoretical findings. Our detailed results on structural, electronic, and optical properties will be discussed with respect to available experimental data, and will provide further insights into the atomistic origin of the disorder-induced band gap lowering in these promising kesterite type materials.
There is a large variety of goals instructors have for laboratory courses, with different courses focusing on different subsets of goals. An often implicit, but crucial, goal is to develop students’ attitudes, views, and expectations about experimental physics to align with practicing experimental physicists. The assessment of laboratory courses upon this one dimension of learning has been intensively studied in U.S. institutions using the Colorado Learning Attitudes about Science Survey for Experimental Physics (E-CLASS). However, there is no such an instrument available to use in Germany, and the influence of laboratory courses on students views about the nature of experimental physics is still unexplored at German-speaking institutions. Motivated by the lack of an assessment tool to investigate this goal in laboratory courses at German-speaking institutions, we present a translated version of the E-CLASS adapted to the context at German-speaking institutions. We call the German version of the E-CLASS, the GE-CLASS. We describe the translation process and the creation of an automated web-based system for instructors to assess their laboratory courses. We also present first results using GE-CLASS obtained at the University of Potsdam. A first comparison between E-CLASS and GE-CLASS results shows clear differences between University of Potsdam and U.S. students’ views and beliefs about experimental physics.
Isoflux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f βc for the translocation time, where β = −1 is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics and the force exponent β = −1 regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant and the resulting entropic force causes the force exponent to drop below unity.
Anomalous-diffusion, the departure of the spreading dynamics of diffusing particles from the traditional law of Brownian-motion, is a signature feature of a large number of complex soft-matter and biological systems. Anomalous-diffusion emerges due to a variety of physical mechanisms, e.g., trapping interactions or the viscoelasticity of the environment. However, sometimes systems dynamics are erroneously claimed to be anomalous, despite the fact that the true motion is Brownian—or vice versa. This ambiguity in establishing whether the dynamics as normal or anomalous can have far-reaching consequences, e.g., in predictions for reaction- or relaxation-laws. Demonstrating that a system exhibits normal- or anomalous-diffusion is highly desirable for a vast host of applications. Here, we present a criterion for anomalous-diffusion based on the method of power-spectral analysis of single trajectories. The robustness of this criterion is studied for trajectories of fractional-Brownian-motion, a ubiquitous stochastic process for the description of anomalous-diffusion, in the presence of two types of measurement errors. In particular, we find that our criterion is very robust for subdiffusion. Various tests on surrogate data in absence or presence of additional positional noise demonstrate the efficacy of this method in practical contexts. Finally, we provide a proof-of-concept based on diverse experiments exhibiting both normal and anomalous-diffusion.
Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations (“Joseph effect”), fat-tailed probability density of increments (“Noah effect”), and nonstationarity (“Moses effect”). We show that such a decomposition of real-life data allows us to infer nontrivial behavioral predictions and to resolve open questions in the fields of single-particle tracking in living cells and movement ecology.
Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusionmodel and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a wellcalibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output.
Spin precession in magnetic materials is commonly modelled with the classical phenomenological Landau-Lifshitz-Gilbert (LLG) equation. Based on a quantized three-dimensional spin + environment Hamiltonian, we here derive a spin operator equation of motion that describes precession and includes a general form of damping that consistently accounts for memory, coloured noise and quantum statistics. The LLG equation is recovered as its classical, Ohmic approximation. We further introduce resonant Lorentzian system-reservoir couplings that allow a systematic comparison of dynamics between Ohmic and non-Ohmic regimes. Finally, we simulate the full non-Markovian dynamics of a spin in the semi-classical limit. At low temperatures, our numerical results demonstrate a characteristic reduction and flattening of the steady state spin alignment with an external field, caused by the quantum statistics of the environment. The results provide a powerful framework to explore general three-dimensional dissipation in quantum thermodynamics.
Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusion model and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a well-calibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output. <br /> Diffusive motions in complex environments such as living biological cells or soft matter systems can be analyzed with single-particle-tracking approaches, where accuracy of output may vary. The authors involve a machine-learning technique for decoding anomalous-diffusion data and provide an uncertainty estimate together with predicted output.
The coupling of the internal mechanisms of cell polarization to cell shape deformations and subsequent cell crawling poses many interdisciplinary scientific challenges. Several mathematical approaches have been proposed to model the coupling of both processes, where one of the most successful methods relies on a phase field that encodes the morphology of the cell, together with the integration of partial differential equations that account for the polarization mechanism inside the cell domain as defined by the phase field. This approach has been previously employed to model the motion of single cells of the social amoeba Dictyostelium discoideum, a widely used model organism to study actin-driven motility and chemotaxis of eukaryotic cells. Besides single cell motility, Dictyostelium discoideum is also well-known for its collective behavior. Here, we extend the previously introduced model for single cell motility to describe the collective motion of large populations of interacting amoebae by including repulsive interactions between the cells. We performed numerical simulations of this model, first characterizing the motion of single cells in terms of their polarity and velocity vectors. We then systematically studied the collisions between two cells that provided the basic interaction scenarios also observed in larger ensembles of interacting amoebae. Finally, the relevance of the cell density was analyzed, revealing a systematic decrease of the motility with density, associated with the formation of transient cell clusters that emerge in this system even though our model does not include any attractive interactions between cells. This model is a prototypical active matter system for the investigation of the emergent collective dynamics of deformable, self-driven cells with a highly complex, nonlinear coupling of cell shape deformations, self-propulsion and repulsive cell-cell interactions. Understanding these self-organization processes of cells like their autonomous aggregation is of high relevance as collective amoeboid motility is part of wound healing, embryonic morphogenesis or pathological processes like the spreading of metastatic cancer cells.
We quantitatively address the conjecture that magnetic helicity must be shed from the Sun by eruptions launching coronal mass ejections in order to limit its accumulation in each hemisphere. By varying the ratio of guide and strapping field and the flux rope twist in a parametric simulation study of flux rope ejection from approximately marginally stable force-free equilibria, different ratios of self- and mutual helicity are set and the onset of the torus or helical kink instability is obtained. The helicity shed is found to vary over a broad range from a minor to a major part of the initial helicity, with self helicity being largely or completely shed and mutual helicity, which makes up the larger part of the initial helicity, being shed only partly. Torus-unstable configurations with subcritical twist and without a guide field shed up to about two-thirds of the initial helicity, while a highly twisted, kink-unstable configuration sheds only about one-quarter. The parametric study also yields stable force-free flux rope equilibria up to a total flux-normalized helicity of 0.25, with a ratio of self- to total helicity of 0.32 and a ratio of flux rope to external poloidal flux of 0.94. These results numerically demonstrate the conjecture of helicity shedding by coronal mass ejections and provide a first account of its parametric dependence. Both self- and mutual helicity are shed significantly; this reduces the total initial helicity by a fraction of ∼0.4--0.65 for typical source region parameters.
The detection of internal irregularities is crucial for quality assessment in metal-based additive manufacturing (AM) technologies such as laser powder bed fusion (L-PBF). The utilization of in-process thermography as an in situ monitoring tool in combination with post-process X-ray micro computed tomography (XCT) as a reference technique has shown great potential for this aim. Due to the small irregularity dimensions, a precise registration of the datasets is necessary as a requirement for correlation. In this study, the registration of thermography and XCT reference datasets of a cylindric specimen containing keyhole pores is carried out for the development of a porosity prediction model. The considered datasets show variations in shape, data type and dimensionality, especially due to shrinkage and material elevation effects present in the manufactured part. Since the resulting deformations are challenging for registration, a novel preprocessing methodology is introduced that involves an adaptive volume adjustment algorithm which is based on the porosity distribution in the specimen. Thus, the implementation of a simple three-dimensional image-to-image registration is enabled. The results demonstrate the influence of the part deformation on the resulting porosity location and the importance of registration in terms of irregularity prediction.
The Dirac point of a topological surface state (TSS) is protected against gapping by time-reversal symmetry. Conventional wisdom stipulates, therefore, that only through magnetisation may a TSS become gapped. However, non-magnetic gaps have now been demonstrated in Bi2Se3 systems doped with Mn or In, explained by hybridisation of the Dirac cone with induced impurity resonances. Recent photoemission experiments suggest that an analogous mechanism applies even when Bi2Se3 is surface dosed with Au. Here, we perform a systematic spin- and angle-resolved photoemission study of Au-dosed Bi2Se3. Although there are experimental conditions wherein the TSS appears gapped due to unfavourable photoemission matrix elements, our photon-energy-dependent spectra unambiguously demonstrate the robustness of the Dirac cone against high Au coverage. We further show how the spin textures of the TSS and its accompanying surface resonances remain qualitatively unchanged following Au deposition, and discuss the mechanism underlying the suppression of the spectral weight.
In this work, the fabrication and characterization of a simple, inexpensive, and effective microfluidic paper analytic device (mu PAD) for monitoring DNA samples is reported. The glass microfiber-based chip has been fabricated by a new wax-based transfer-printing technique and an electrode printing process. It is capable of moving DNA effectively in a time-dependent fashion. The nucleic acid sample is not damaged by this process and is accumulated in front of the anode, but not directly on the electrode. Thus, further DNA processing is feasible. The system allows the DNA to be purified by separating it from other components in sample mixtures such as proteins. Furthermore, it is demonstrated that DNA can be moved through several layers of the glass fiber material. This proof of concept will provide the basis for the development of rapid test systems, e.g., for the detection of pathogens in water samples.
Neutrophil granulocytes are essential for the first host defense. After leaving the blood circulation they migrate efficiently towards sites of inflammation. They are guided by chemoattractants released from cells within the inflammatory foci. On a cellular level, directional migration is a consequence of cellular front-rear asymmetry which is induced by the concentration gradient of the chemoattractants. The generation and maintenance of this asymmetry, however, is not yet fully understood. Here we analyzed the paths of chemotacting neutrophils with different stochastic models to gain further insight into the underlying mechanisms. Wildtype chemotacting neutrophils show an anomalous superdiffusive behavior. CXCR2 blockade and TRPC6-knockout cause the tempering of temporal correlations and a reduction of chemotaxis. Importantly, such tempering is found both in vitro and in vivo. These findings indicate that the maintenance of anomalous dynamics is crucial for chemotactic behavior and the search efficiency of neutrophils.
The motility of neutrophils and their ability to sense and to react to chemoattractants in their environment are of central importance for the innate immunity. Neutrophils are guided towards sites of inflammation following the activation of G-protein coupled chemoattractant receptors such as CXCR2 whose signaling strongly depends on the activity of Ca2+ permeable TRPC6 channels. It is the aim of this study to analyze data sets obtained in vitro (murine neutrophils) and in vivo (zebrafish neutrophils) with a stochastic mathematical model to gain deeper insight into the underlying mechanisms. The model is based on the analysis of trajectories of individual neutrophils. Bayesian data analysis, including the covariances of positions for fractional Brownian motion as well as for exponentially and power-law tempered model variants, allows the estimation of parameters and model selection. Our model-based analysis reveals that wildtype neutrophils show pure superdiffusive fractional Brownian motion. This so-called anomalous dynamics is characterized by temporal long-range correlations for the movement into the direction of the chemotactic CXCL1 gradient. Pure superdiffusion is absent vertically to this gradient. This points to an asymmetric 'memory' of the migratory machinery, which is found both in vitro and in vivo. CXCR2 blockade and TRPC6-knockout cause tempering of temporal correlations in the chemotactic gradient. This can be interpreted as a progressive loss of memory, which leads to a marked reduction of chemotaxis and search efficiency of neutrophils. In summary, our findings indicate that spatially differential regulation of anomalous dynamics appears to play a central role in guiding efficient chemotactic behavior.
The motility of adherent eukaryotic cells is driven by the dynamics of the actin cytoskeleton. Despite the common force-generating actin machinery, different cell types often show diverse modes of locomotion that differ in their shape dynamics, speed, and persistence of motion. Recently, experiments in Dictyostelium discoideum have revealed that different motility modes can be induced in this model organism, depending on genetic modifications, developmental conditions, and synthetic changes of intracellular signaling. Here, we report experimental evidence that in a mutated D. discoideum cell line with increased Ras activity, switches between two distinct migratory modes, the amoeboid and fan-shaped type of locomotion, can even spontaneously occur within the same cell. We observed and characterized repeated and reversible switchings between the two modes of locomotion, suggesting that they are distinct behavioral traits that coexist within the same cell. We adapted an established phenomenological motility model that combines a reaction-diffusion system for the intracellular dynamics with a dynamic phase field to account for our experimental findings.
How do different reset protocols affect ergodicity of a diffusion process in single-particle-tracking experiments? We here address the problem of resetting of an arbitrary stochastic anomalous-diffusion process (ADP) from the general mathematical points of view and assess ergodicity of such reset ADPs for an arbitrary resetting protocol. The process of stochastic resetting describes the events of the instantaneous restart of a particle’s motion via randomly distributed returns to a preset initial position (or a set of those). The waiting times of such resetting events obey the Poissonian, Gamma, or more generic distributions with specified conditions regarding the existence of moments. Within these general approaches, we derive general analytical results and support them by computer simulations for the behavior of the reset mean-squared displacement (MSD), the new reset increment-MSD (iMSD), and the mean reset time-averaged MSD (TAMSD). For parental nonreset ADPs with the MSD(t)∝ tμ we find a generic behavior and a switch of the short-time growth of the reset iMSD and mean reset TAMSDs from ∝ _μ for subdiffusive to ∝ _1 for superdiffusive reset ADPs. The critical condition for a reset ADP that recovers its ergodicity is found to be more general than that for the nonequilibrium stationary state, where obviously the iMSD and the mean TAMSD are equal. The consideration of the new statistical quantifier, the iMSD—as compared to the standard MSD—restores the ergodicity of an arbitrary reset ADP in all situations when the μth moment of the waiting-time distribution of resetting events is finite. Potential applications of these new resetting results are, inter alia, in the area of biophysical and soft-matter systems.
Focusing on transient chaos
(2022)
Recent advances in the field of complex, transiently chaotic dynamics are reviewed, based on the results published in the focus issue of J. Phys. Complex. on this topic. One group of achievements concerns network dynamics where transient features are intimately related to the degree and stability of synchronization, as well as to the network topology. A plethora of various applications of transient chaos are described, ranging from the collective motion of active particles, through the operation of power grids, cardiac arrhythmias, and magnetohydrodynamical dynamos, to the use of machine learning to predict time evolutions. Nontraditional forms of transient chaos are also explored, such as the temporal change of the chaoticity in the transients (called doubly transient chaos), as well as transients in systems subjected to parameter drift, the paradigm of which is climate change.
We introduce and study a Lévy walk (LW) model of particle spreading with a finite propagation speed combined with soft resets, stochastically occurring periods in which an harmonic external potential is switched on and forces the particle towards a specific position. Soft resets avoid instantaneous relocation of particles that in certain physical settings may be considered unphysical. Moreover, soft resets do not have a specific resetting point but lead the particle towards a resetting point by a restoring Hookean force. Depending on the exact choice for the LW waiting time density and the probability density of the periods when the harmonic potential is switched on, we demonstrate a rich emerging response behaviour including ballistic motion and superdiffusion. When the confinement periods of the soft-reset events are dominant, we observe a particle localisation with an associated non-equilibrium steady state. In this case the stationary particle probability density function turns out to acquire multimodal states. Our derivations are based on Markov chain ideas and LWs with multiple internal states, an approach that may be useful and flexible for the investigation of other generalised random walks with soft and hard resets. The spreading efficiency of soft-rest LWs is characterised by the first-passage time statistic.
We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials.
Perovskite solar cells represent one of the recent success stories in photovoltaics. The device efficiency has been steadily increasing over the past years, but further work is needed to enhance the performance, for example, through the reduction of defects to prevent carrier recombination. SCAPS-1D simulations were performed to assess efficiency limits and identify approaches to decrease the impact of defects, through the selection of an optimal hole-transport material and a hole-collecting electrode. Particular attention was given to evaluation of the influence of bulk defects within light-absorbing CH3NH3SnI3 layers. In addition, the study demonstrates the influence of interface defects at the TiO2/CH3NH3SnI3 (IL1) and CH3NH3SnI3/HTL (IL2) interfaces across the similar range of defect densities. Finally, the optimal device architecture TiO2/CH3NH3SnI3/Cu2O is proposed for the given absorber layer using the readily available Cu2O hole-transporting material with PCE = 27.95%, FF = 84.05%, V-OC = 1.02 V and J(SC) = 32.60 mA/cm(2), providing optimal performance and enhanced resistance to defects.
The application of the fractional calculus in the mathematical modelling of relaxation processes in complex heterogeneous media has attracted a considerable amount of interest lately.
The reason for this is the successful implementation of fractional stochastic and kinetic equations in the studies of non-Debye relaxation.
In this work, we consider the rotational diffusion equation with a generalised memory kernel in the context of dielectric relaxation processes in a medium composed of polar molecules. We give an overview of existing models on non-exponential relaxation and introduce an exponential resetting dynamic in the corresponding process.
The autocorrelation function and complex susceptibility are analysed in detail.
We show that stochastic resetting leads to a saturation of the autocorrelation function to a constant value, in contrast to the case without resetting, for which it decays to zero. The behaviour of the autocorrelation function, as well as the complex susceptibility in the presence of resetting, confirms that the dielectric relaxation dynamics can be tuned by an appropriate choice of the resetting rate.
The presented results are general and flexible, and they will be of interest for the theoretical description of non-trivial relaxation dynamics in heterogeneous systems composed of polar molecules.
Science education researchers have developed a refined understanding of the structure of science teachers’ pedagogical content knowledge (PCK), but how to develop applicable and situation-adequate PCK remains largely unclear. A potential problem lies in the diverse conceptualisations of the PCK used in PCK research. This study sought to systematize existing science education research on PCK through the lens of the recently proposed refined consensus model (RCM) of PCK. In this review, the studies’ approaches to investigating PCK and selected findings were characterised and synthesised as an overview comparing research before and after the publication of the RCM. We found that the studies largely employed a qualitative case-study methodology that included specific PCK models and tools. However, in recent years, the studies focused increasingly on quantitative aspects. Furthermore, results of the reviewed studies can mostly be integrated into the RCM. We argue that the RCM can function as a meaningful theoretical lens for conceptualizing links between teaching practice and PCK development by proposing pedagogical reasoning as a mechanism and/or explanation for PCK development in the context of teaching practice.
Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean.
Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower.
The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction.
Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.
The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth.
Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model.
Our model results show that for warming levels above 2 degrees C, Greenland could become essentially ice-free within several millennia, mainly as a result of surface melting and acceleration of ice flow. These ice losses are mitigated, however, in some cases with strong GIA feedback even promoting an incomplete recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming - from incomplete recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse.
In the incomplete recovery regime, the initial ice loss due to warming is essentially reversed within 50 000 years, and the ice volume stabilizes at 61 %-93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods between 74 000 and over 300 000 years and oscillation amplitudes between 15 %-70 % of present-day ice volume.
This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on timescales on the order of 100 000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet.
Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy
(2022)
The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.
In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of approximate to 11% initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.
We present observations of three-dimensional magnetic power spectra in wavevector space to investigate the anisotropy and scalings of sub-Alfvenic solar wind turbulence at magnetohydrodynamic (MHD) scale using the Magnetospheric Multiscale spacecraft. The magnetic power distributions are organized in a new coordinate determined by wavevectors ((kappa) over cap) and background magnetic field ((b) over cap (0)) in Fourier space. This study utilizes two approaches to determine wavevectors: the singular value decomposition method and multispacecraft timing analysis. The combination of the two methods allows an examination of the properties of magnetic field fluctuations in terms of mode compositions without any spatiotemporal hypothesis. Observations show that fluctuations (delta B-perpendicular to 1) in the direction perpendicular to (kappa) over cap and (b) over cap (0) prominently cascade perpendicular to (b) over cap (0), and such anisotropy increases with wavenumbers. The reduced power spectra of 6.8 11 follow Goldreich-Sridhar scalings: (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-5/3) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-2). In contrast, fluctuations within the (k) over cap(b) over cap (0) plane show isotropic behaviors: perpendicular power distributions are approximately the same as parallel distributions. The reduced power spectra of fluctuations within the (k) over cap(b) over cap (0) plane follow the scalings (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-3/2) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-3/2). Comparing frequency-wavevector spectra with theoretical dispersion relations of MHD modes, we find that delta B-perpendicular to 1 are probably associated with Alfven modes. On the other hand, magnetic field fluctuations within the (k) over cap(b) over cap (0) plane more likely originate from fast modes based on their isotropic behaviors. The observations of anisotropy and scalings of different magnetic field components are consistent with the predictions of current compressible MHD theory. Moreover, for the Alfvenic component, the ratio of cascading time to the wave period is found to be a factor of a few, consistent with critical balance in the strong turbulence regime. These results are valuable for further studies of energy compositions of plasma turbulence and their effects on energetic particle transport.
Starting from the observation that the reduced state of a system strongly coupled to a bath is, in general, an athermal state, we introduce and study a cyclic battery-charger quantum device that is in thermal equilibrium, or in a ground state, during the charge storing stage. The cycle has four stages: the equilibrium storage stage is interrupted by disconnecting the battery from the charger, then work is extracted from the battery, and then the battery is reconnected with the charger; finally, the system is brought back to equilibrium. At no point during the cycle are the battery-charger correlations artificially erased. We study the case where the battery and charger together comprise a spin-1/2 Ising chain, and show that the main characteristics-the extracted energy and the thermodynamic efficiency-can be enhanced by operating the cycle close to the quantum phase transition point. When the battery is just a single spin, we find that the output work and efficiency show a scaling behavior at criticality and derive the corresponding critical exponents. Due to always present correlations between the battery and the charger, operations that are equivalent from the perspective of the battery can entail different energetic costs for switching the battery-charger coupling. This happens only when the coupling term does not commute with the battery's bare Hamiltonian, and we use this purely quantum leverage to further optimize the performance of the device.
This paper studies cosmic-ray (CR) transport in magnetohydrodynamic (MHD) turbulence. CR transport is strongly dependent on the properties of the magnetic turbulence.
We perform test particle simulations to study the interactions of CR with both total MHD turbulence and decomposed MHD modes.
The spatial diffusion coefficients and the pitch angle scattering diffusion coefficients are calculated from the test particle trajectories in turbulence.
Our results confirm that the fast modes dominate the CR propagation, whereas Alfven and slow modes are much less efficient and have shown similar pitch-angle scattering rates.
We investigate the cross field transport on large and small scales. On large/global scales, normal diffusion is observed and the diffusion coefficient is suppressed by M-A(zeta) compared to the parallel diffusion coefficients, with zeta closer to 4 in Alfven modes than that in total turbulence, as theoretically expected.
For the CR transport on scales smaller than the turbulence injection scale, both the local and global magnetic reference frames are adopted. Superdiffusion is observed on such small scales in all the cases. Particularly, CR transport in Alfven modes show clear Richardson diffusion in the local reference frame. The diffusion transitions smoothly from the Richardson's one with index 1.5 to normal diffusion as the particle mean free path decreases from lambda(parallel to) >> L to lambda(parallel to) << L, where L is the injection/coherence length of turbulence.
Our results have broad applications to CRs in various astrophysical environments.
Tailoring the secondary surface morphology of electro-spun nanofibers has been highly desired, as such delicate structures equip nanofibers with distinct functions. Here, we report a simple strategy to directly reconstruct the surface of polyvinyl alcohol/polyvinylpyrrolidone (PVA/PVP) nanofibers by water evaporation. The roughness and diameter of the nanofibers depend on the temperature during vacuum drying. Surface changes of the nanofibers from smooth to rough were observed at 55 degrees C, with a significant drop in nanofiber diameter. We attribute the formation of the secondary surface morphology to the intermolecular forces in the water vapor, including capillary and the compression forces, on the basis of the results from the Fourier-transform infrared (FTIR) and X-ray photoelectron (XPS) spectroscopy. The strategy is universally effective for various electro-spun polymer nanofibers, thus opening up avenues toward more detailed and sophisticated structure design and implementation for nanofibers.
Point-of-care and in-vivo bio-diagnostic tools are the current need for the present critical scenarios in the healthcare industry. The past few decades have seen a surge in research activities related to solving the challenges associated with precise on-site bio-sensing. Cutting-edge fiber optic technology enables the interaction of light with functionalized fiber surfaces at remote locations to develop a novel, miniaturized and cost-effective lab on fiber technology for bio-sensing applications. The recent remarkable developments in the field of nanotechnology provide innumerable functionalization methodologies to develop selective bio-recognition elements for label free biosensors. These exceptional methods may be easily integrated with fiber surfaces to provide highly selective light-matter interaction depending on various transduction mechanisms. In the present review, an overview of optical fiber-based biosensors has been provided with focus on physical principles used, along with the functionalization protocols for the detection of various biological analytes to diagnose the disease. The design and performance of these biosensors in terms of operating range, selectivity, response time and limit of detection have been discussed. In the concluding remarks, the challenges associated with these biosensors and the improvement required to develop handheld devices to enable direct target detection have been highlighted.
Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice melange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice melange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-melange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.
We consider a one-dimensional oscillatory medium with a coupling through a diffusive linear field. In the limit of fast diffusion this setup reduces to the classical Kuramoto–Battogtokh model. We demonstrate that for a finite diffusion stable chimera solitons, namely localized synchronous domain in an infinite asynchronous environment, are possible. The solitons are stable also for finite density of oscillators, but in this case they sway with a nearly constant speed. This finite-density-induced motility disappears in the continuum limit, as the velocity of the solitons is inverse proportional to the density. A long-wave instability of the homogeneous asynchronous state causes soliton turbulence, which appears as a sequence of soliton mergings and creations. As the instability of the asynchronous state becomes stronger, this turbulence develops into a spatio-temporal intermittency.
Diffraction enhanced imaging (DEI) is an advanced digital radiographic imaging technique employing the refraction of X-rays to contrast internal interfaces. This study aims to qualitatively and quantitatively evaluate images acquired using this technique and to assess how different fitting functions to the typical rocking curves (RCs) influence the quality of the images. RCs are obtained for every image pixel. This allows the separate determination of the absorption and the refraction properties of the material in a position-sensitive manner. Comparison of various types of fitting functions reveals that the Pseudo-Voigt (PsdV) function is best suited to fit typical RCs. A robust algorithm was developed in the Python programming language, which reliably extracts the physically meaningful information from each pixel of the image. We demonstrate the potential of the algorithm with two specimens: a silicone gel specimen that has well-defined interfaces, and an additively manufactured polycarbonate specimen.
In this study we analyze the storm-time evolution of equatorial electron pitch angle distributions (PADs) in the outer radiation belt region using observations from the Magnetic Electron Ion Spectrometer (MagEIS) instrument aboard the Van Allen Probes in 2012-2019. The PADs are approximated using a sum of the first, third and fifth sine harmonics. Different combinations of the respective coefficients refer to the main PAD shapes within the outer radiation belt, namely the pancake, flat-top, butterfly and cap PADs. We conduct a superposed epoch analysis of 129 geomagnetic storms and analyze the PAD evolution for day and night MLT sectors. PAD shapes exhibit a strong energy-dependent response. At energies of tens of keV, the PADs exhibit little variation throughout geomagnetic storms. Cap PADs are mainly observed at energies < 300 keV, and their extent in L shrinks with increasing energy. The cap distributions transform into the pancake PADs around the main phase of the storm on the nightside, and then come back to their original shapes during the recovery phase. At higher energies on the dayside, the PADs are mainly pancake during pre-storm conditions and become more anisotropic during the main phase. The quiet-time butterfly PADs can be observed on the nightside at L> 5.6. During the main phase, butterfly PADs have stronger 90 degrees-minima and can be observed at lower L-shells (down to L = 5), then transitioning into flat-top PADs at L similar to 4.5 - 5 and pancake PADs at L < 4.5. The resulting PAD coefficients for different energies, locations and storm epochs can be used to test the wave models and physics-based radiation belt codes in terms of pitch angle distributions.
Predicting the electron population of Earth's ring current during geomagnetic storms still remains a challenging task.
In this work, we investigate the sensitivity of 10 keV ring current electrons to different driving processes, parameterised by the Kp index, during several moderate and intense storms.
Results are validated against measurements from the Van Allen Probes satellites. Perturbing the Kp index allows us to identify the most dominant processes for moderate and intense storms respectively.
We find that during moderate storms (Kp < 6) the drift velocities mostly control the behaviour of low energy electrons, while loss from wave-particle interactions is the most critical parameter for quantifying the evolution of intense storms (Kp > 6). Perturbations of the Kp index used to drive the boundary conditions at GEO and set the plasmapause location only show a minimal effect on simulation results over a limited L range.
It is further shown that the flux at L & SIM; 3 is more sensitive to changes in the Kp index compared to higher L shells, making it a good proxy for validating the source-loss balance of a ring current model.
Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.
Assessing the impact of hydrogen absorption on the characteristics of the Galactic center excess
(2022)
We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation transport modeling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic center.
The main benefits of the new gas model are (a) the ability to reproduce the negative line signals seen with the HI4PI survey and (b) the accounting for gas that primarily manifests itself through absorption.
We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter.
We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H i model is used to estimate the cosmic-ray induced diffuse gamma-ray emission.
The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence of a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favored by the data.
This finding is robust under the variation of various parameters, for example, the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.
To undergo diffusive shock acceleration, electrons need to be preaccelerated to increase their energies by several orders of magnitude, else their gyroradii will be smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel nor perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the preacceleration in this region, we undertake particle-in-cell simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for theta (Bn) = 74.degrees 3, it increases to R similar to 5% for theta (Bn) = 30 degrees, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to R (0.6) and a wavelength approximate to 2 lambda (se), where lambda (se) is the electron skin length. While the initial reflection mechanism is typically a combination of shock-surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream, they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In less than or similar to 1% of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction is rescattered back toward the shock after reflection, reinteracts with the shock with energies much greater than thermal, and crosses into the downstream.
Suppression of the TeV Pair-beam-Plasma Instability by a Tangled Weak Intergalactic Magnetic Field
(2022)
We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kiloparsec, we find that weak magnetic fields increase the transverse momentum of the pair-beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission.
The current paradigm of cosmic-ray (CR) origin states that the greater part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after a supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. With this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely to be triggered at high Mach number shocks.
We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta.
We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times.
We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient.
When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity.
Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the first binary neutron star simulations of this kind. We find that the presence of dark matter reduces the lifetime of the merger remnant and favors a prompt collapse to a black hole. Furthermore, we find differences in the merger time for systems with the same total mass and mass ratio, but different amounts of dark matter. Finally, we find that electromagnetic signals produced by the merger of binary neutron stars admixed with dark matter are very unlikely to be as bright as their dark matter-free counterparts. Given the increased sensitivity of multi-messenger facilities, our analysis gives a new perspective on how to probe the presence of dark matter.
To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB Lambda phi. Overall, our test simulations have good agreement with those reported in the literature.
Organic photovoltaics (PV) is an energy-harvesting technology that offers many advantages, such as flexibility, low weight and cost, as well as environmentally benign materials and manufacturing techniques. Despite growth of power conversion efficiencies to around 19 % in the last years, organic PVs still lag behind inorganic PV technologies, mainly due to high losses in open-circuit voltage. Understanding and improving open circuit voltage in organic solar cells is challenging, as it is controlled by the properties of a donor-acceptor interface where the optical excitations are separated into charge carriers. Here, we provide an electrostatic model of a rough donor-acceptor interface and test it experimentally on small molecule PV materials systems. The model provides concise relationships between the open-circuit voltage, photovoltaic gap, charge-transfer state energy, and interfacial morphology. In particular, we show that the electrostatic bias generated across the interface reduces the photovoltaic gap. This negative influence on open-circuit voltage can, however, be circumvented by adjusting the morphology of the donor-acceptor interface.
Organic solar cells, despite their high power conversion efficiencies, suffer from open circuit voltage losses making them less appealing in terms of applications. Here, the authors, supported with experimental data on small molecule photovoltaic cells, relate open circuit voltage to photovoltaic gap, charge-transfer state energy, and donor-acceptor interfacial morphology.
Inverted perovskite solar cells still suffer from significant non-radiative recombination losses at the perovskite surface and across the perovskite/C-60 interface, limiting the future development of perovskite-based single- and multi-junction photovoltaics. Therefore, more effective inter- or transport layers are urgently required. To tackle these recombination losses, we introduce ortho-carborane as an interlayer material that has a spherical molecular structure and a three-dimensional aromaticity. Based on a variety of experimental techniques, we show that ortho-carborane decorated with phenylamino groups effectively passivates the perovskite surface and essentially eliminates the non-radiative recombination loss across the perovskite/C-60 interface with high thermal stability. We further demonstrate the potential of carborane as an electron transport material, facilitating electron extraction while blocking holes from the interface. The resulting inverted perovskite solar cells deliver a power conversion efficiency of over 23% with a low non-radiative voltage loss of 110mV, and retain >97% of the initial efficiency after 400h of maximum power point tracking. Overall, the designed carborane based interlayer simultaneously enables passivation, electron-transport and hole-blocking and paves the way toward more efficient and stable perovskite solar cells. Effective transport layers are essential to suppress non-radiative recombination losses. Here, the authors introduce phenylamino-functionalized ortho-carborane as an interfacial layer, and realise inverted perovskite solar cells with efficiency of over 23% and operational stability of T97=400h.
Inverted perovskite solar cells still suffer from significant non-radiative recombination losses at the perovskite surface and across the perovskite/C₆₀ interface, limiting the future development of perovskite-based single- and multi-junction photovoltaics. Therefore, more effective inter- or transport layers are urgently required. To tackle these recombination losses, we introduce ortho-carborane as an interlayer material that has a spherical molecular structure and a three-dimensional aromaticity. Based on a variety of experimental techniques, we show that ortho-carborane decorated with phenylamino groups effectively passivates the perovskite surface and essentially eliminates the non-radiative recombination loss across the perovskite/C₆₀ interface with high thermal stability. We further demonstrate the potential of carborane as an electron transport material, facilitating electron extraction while blocking holes from the interface. The resulting inverted perovskite solar cells deliver a power conversion efficiency of over 23% with a low non-radiative voltage loss of 110 mV, and retain >97% of the initial efficiency after 400 h of maximum power point tracking. Overall, the designed carborane based interlayer simultaneously enables passivation, electron-transport and hole-blocking and paves the way toward more efficient and stable perovskite solar cells.
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.
We use ultrafast x-ray diffraction to investigate the effect of expansive phononic and contractive magnetic stress driving the picosecond strain response of a metallic perovskite SrRuO3 thin film upon femtosecond laser excitation. We exemplify how the anisotropic bulk equilibrium thermal expansion can be used to predict the response of the thin film to ultrafast deposition of energy. It is key to consider that the laterally homogeneous laser excitation changes the strain response compared to the near-equilibrium thermal expansion because the balanced in-plane stresses suppress the Poisson stress on the picosecond timescale. We find a very large negative Grüneisen constant describing the large contractive stress imposed by a small amount of energy in the spin system. The temperature and fluence dependence of the strain response for a double-pulse excitation scheme demonstrates the saturation of the magnetic stress in the high-fluence regime.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.
Reciprocal space slicing
(2021)
An experimental technique that allows faster assessment of out-of-plane strain dynamics of thin film heterostructures via x-ray diffraction is presented. In contrast to conventional high-speed reciprocal space-mapping setups, our approach reduces the measurement time drastically due to a fixed measurement geometry with a position-sensitive detector. This means that neither the incident (ω) nor the exit (2θ) diffraction angle is scanned during the strain assessment via x-ray diffraction. Shifts of diffraction peaks on the fixed x-ray area detector originate from an out-of-plane strain within the sample. Quantitative strain assessment requires the determination of a factor relating the observed shift to the change in the reciprocal lattice vector. The factor depends only on the widths of the peak along certain directions in reciprocal space, the diffraction angle of the studied reflection, and the resolution of the instrumental setup. We provide a full theoretical explanation and exemplify the concept with picosecond strain dynamics of a thin layer of NbO2.
'Complex systems are information processors' is a statement that is frequently made. Here we argue for the distinction between information processing-in the sense of encoding and transmitting a symbolic representation-and the formation of correlations (pattern formation/self-organisation). The study of both uses tools from information theory, but the purpose is very different in each case: explaining the mechanisms and understanding the purpose or function in the first case, versus data analysis and correlation extraction in the latter. We give examples of both and discuss some open questions. The distinction helps focus research efforts on the relevant questions in each case.
Non-geminate recombination, as one of the most relevant loss mechanisms in organic and perovskite solar cells, deserves special attention in research efforts to further increase device performance. It can be subdivided into first, second, and third order processes, which can be elucidated by the effects that they have on the time-dependent open-circuit voltage decay. In this study, analytical expressions for the open-circuit voltage decay exhibiting one of the aforementioned recombination mechanisms were derived. It was possible to support the analytical models with experimental examples of three different solar cells, each of them dominated either by first (PBDBT:CETIC-4F), second (PM6:Y6), or third (irradiated CH3NH3PbI3) order recombination. Furthermore, a simple approach to estimate the dominant recombination process was also introduced and tested on these examples. Moreover, limitations of the analytical models and the measurement technique itself were discussed.
BACKGROUND
Anticipating changes in international migration patterns is useful for demographic studies and for designing policies that support the well-being of those involved. Existing forecasting methods do not account for a number of stylized facts that emerge from large-scale migration observations and theories: existing migrant communities - diasporas - act to lower migration costs and thereby provide a mechanism of self-amplification; return migration and transit migration are important components of global migration flows; and poverty constrains emigration.
OBJECTIVE
Here we present hindcasts and future projections of international migration that explicitly account for these nonlinear features.
METHODS
We develop a dynamic model that simulates migration flows by origin, destination, and place of birth. We calibrate the model using recently constructed global datasets of bilateral migration.
RESULTS
We show that the model reproduces past patterns and trends well based only on initial migrant stocks and changes in national incomes. We then project migration flows under future scenarios of global socioeconomic development.
CONCLUSIONS
Different assumptions about income levels and between-country inequality lead to markedly different migration trajectories, with migration flows either converging towards net zero if incomes in presently poor countries catch up with the rest of the world; or remaining high or even rising throughout the 21st century if economic development is slower and more unequal. Importantly, diasporas induce significant inertia and sizable return migration flows.
The layered dichalcogenide MoS2 is relevant for electrochemical Li adsorption/intercalation, in the course of which the material undergoes a concomitant structural phase transition from semiconducting 2H-MoS2 to metallic 1T-LixMoS2. With the core hole clock approach at the S L1 X-ray absorption edge we quantify the ultrafast directional charge transfer of excited S3p electrons in-plane () and out-of-plane (perpendicular to) for 2H-MoS2 as tau 2H,=0.38 +/- 0.08 fs and tau 2H,perpendicular to =0.33 +/- 0.06 fs and for 1T-LixMoS2 as tau 1T,=0.32 +/- 0.12 fs and tau 1T,perpendicular to =0.09 +/- 0.07 fs. The isotropic charge delocalization of S3p electrons in the semiconducting 2H phase within the S-Mo-S sheets is assigned to the specific symmetry of the Mo-S bonding arrangement. Formation of 1T-LixMoS2 by lithiation accelerates the in-plane charge transfer by a factor of similar to 1.2 due to electron injection to the Mo-S covalent bonds and concomitant structural repositioning of S atoms within the S-Mo-S sheets. For excitation into out-of-plane orbitals, an accelerated charge transfer by a factor of similar to 3.7 upon lithiation occurs due to S-Li coupling.
The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Pots-dam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3 degrees respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice-ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.
The reconstruction of cone-beam computed tomography data using filtered back-projection algorithms unavoidably results in severe artefacts. We describe how the Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) algorithm can be combined with a model of the artefacts for the reconstruction of such data. The implementation of DIRECTT results in reconstructed volumes of superior quality compared to the conventional algorithms.
Organic solar cells have the potential to become the cheapest form of electricity. Rapid increase in the power conversion efficiency of organic solar cells (OSCs) has been achieved with the development of non-fullerene small-molecule acceptors. Next generation photovoltaics based upon environmentally benign "green solvent" processing of organic semiconductors promise a step-change in the adaptability and versatility of solar technologies and promote sustainable development. However, high-performing OSCs are still processed by halogenated (non-environmentally friendly) solvents, so hindering their large-scale manufacture. In this perspective, we discuss the recent progress in developing highly efficient OSCs processed from eco-compatible solvents, and highlight research challenges that should be addressed for the future development of high power conversion efficiencies devices.
We derive modified reflection coefficients for electromagnetic waves in the THz and far infrared range. The idea is based on hydrodynamic boundary conditions for metallic conduction electrons. The temperature-dependent part of the Casimir pressure between metal plates is evaluated. The results should shed light on the "thermal anomaly," where measurements deviate from the standard fluctuation electrodynamics for conducting metals.
In this paper, we report X-ray absorption and core-level electron spectra of the nucleobase derivative 2-thiouracil at the sulfur L1- and L2,3-edges. We used soft X-rays from the free-electron laser FLASH2 for the excitation of isolated molecules and dispersed the outgoing electrons with a magnetic bottle spectrometer. We identified photoelectrons from the 2p core orbital, accompanied by an electron correlation satellite, as well as resonant and non-resonant Coster–Kronig and Auger–Meitner emission at the L1- and L2,3-edges, respectively. We used the electron yield to construct X-ray absorption spectra at the two edges. The experimental data obtained are put in the context of the literature currently available on sulfur core-level and 2-thiouracil spectroscopy.
The stratosphere is one of the main potential sources for subseasonal to seasonal predictability in midlatitudes in winter. The ability of an atmospheric model to realistically simulate the stratospheric dynamics is essential in order to move forward in the field of seasonal predictions in midlatitudes. Earlier studies with the ICOsahedral Nonhydrostatic atmospheric model (ICON) point out that stratospheric westerlies in ICON are underestimated. This is the first extensive study on the evaluation of Northern Hemisphere stratospheric winter circulation with ICON in numerical weather prediction (NWP) mode. Seasonal experiments with the default setup are able to reproduce the basic climatology of the stratospheric polar vortex. However, westerlies are too weak and major stratospheric warmings too frequent in ICON. Both a reduction of the nonorographic, and a reduction of the orographic gravity wave and wake drag lead to a strengthening of the stratospheric vortex and a bias reduction, in particular in January. However, the effect of the nonorographic gravity wave drag scheme on the stratosphere is stronger. Stratosphere-troposphere coupling is intensified and more realistic due to a reduced gravity wave drag. Furthermore, an adjustment of the subgrid-scale orographic drag parameterization leads to a significant error reduction in the mean sea level pressure. As a result of these findings, we present our current suggested improved setup for seasonal experiments with ICON-NWP. <br /> Plain Language Summary Although seasonal forecasts for midlatitudes have the potential to be highly beneficial to the public sector, they are still characterized by a large amount of uncertainty. Exact simulations of the circulation in the stratosphere can help to improve tropospheric predictability on seasonal time scales. For this reason, we investigate how well the new German atmospheric model is able to simulate the stratospheric circulation. The model reproduces the basic behavior of the Northern Hemisphere stratospheric polar vortex, but the westerly circulation in winter is underestimated. The stratospheric circulation is influenced by gravity waves that exert drag on the flow. These processes are only partly physically represented in the model, but are very important and are hence parameterized. By adjusting the parameterizations for the gravity wave drag, the stratospheric polar vortex is strengthened, thereby yielding a more realistic stratospheric circulation. In addition, the altered parameterizations improve the simulated surface pressure pattern. Based upon this, we present our current suggested improved model setup for seasonal experiments.
I study deterministic dynamics of chiral active particles in two dimensions. Particles are considered as discs interacting with elastic repulsive forces. An ensemble of particles, started from random initial conditions, demonstrates chaotic collisions resulting in their normal diffusion. This chaos is transient, as rather abruptly a synchronous collisionless state establishes. The life time of chaos grows exponentially with the number of particles. External forcing (periodic or chaotic) is shown to facilitate the synchronization transition.
We consider a social-type network of coupled phase oscillators. Such a network consists of an active core of mutually interacting elements, and of a flock of passive units, which follow the driving from the active elements, but otherwise are not interacting. We consider a ring geometry with a long-range coupling, where active oscillators form a fluctuating chimera pattern. We show that the passive elements are strongly correlated. This is explained by negative transversal Lyapunov exponents.
Computation of the instantaneous phase and amplitude via the Hilbert Transform is a powerful tool of data analysis. This approach finds many applications in various science and engineering branches but is not proper for causal estimation because it requires knowledge of the signal’s past and future. However, several problems require real-time estimation of phase and amplitude; an illustrative example is phase-locked or amplitude-dependent stimulation in neuroscience. In this paper, we discuss and compare three causal algorithms that do not rely on the Hilbert Transform but exploit well-known physical phenomena, the synchronization and the resonance. After testing the algorithms on a synthetic data set, we illustrate their performance computing phase and amplitude for the accelerometer tremor measurements and a Parkinsonian patient’s beta-band brain activity.
Master equations are a vital tool to model heat flow through nanoscale thermodynamic systems. Most practical devices are made up of interacting subsystems and are often modelled using either local master equations (LMEs) or global master equations (GMEs). While the limiting cases in which either the LME or the GME breaks down are well understood, there exists a 'grey area' in which both equations capture steady-state heat currents reliably but predict very different transient heat flows. In such cases, which one should we trust? Here we show that, when it comes to dynamics, the local approach can be more reliable than the global one for weakly interacting open quantum systems. This is due to the fact that the secular approximation, which underpins the GME, can destroy key dynamical features. To illustrate this, we consider a minimal transport setup and show that its LME displays exceptional points (EPs). These singularities have been observed in a superconducting-circuit realisation of the model [1]. However, in stark contrast to experimental evidence, no EPs appear within the global approach. We then show that the EPs are a feature built into the Redfield equation, which is more accurate than the LME and the GME. Finally, we show that the local approach emerges as the weak-interaction limit of the Redfield equation, and that it entirely avoids the secular approximation.
We analyze historical data of stock-market prices for multiple financial indices using the concept of delay-time averaging for the financial time series (FTS). The region of validity of our recent theoretical predictions [Cherstvy A G et al 2017 New J. Phys. 19 063045] for the standard and delayed time-averaged mean-squared 'displacements' (TAMSDs) of the historical FTS is extended to all lag times. As the first novel element, we perform extensive computer simulations of the stochastic differential equation describing geometric Brownian motion (GBM) which demonstrate a quantitative agreement with the analytical long-term price-evolution predictions in terms of the delayed TAMSD (for all stock-market indices in crisis-free times). Secondly, we present a robust procedure of determination of the model parameters of GBM via fitting the features of the price-evolution dynamics in the FTS for stocks and cryptocurrencies. The employed concept of single-trajectory-based time averaging can serve as a predictive tool (proxy) for a mathematically based assessment and rationalization of probabilistic trends in the evolution of stock-market prices.
We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.
We study the first passage dynamics for a diffusing particle experiencing a spatially varying diffusion coefficient while driven by correlated additive Gaussian white noise and multiplicative coloured non-Gaussian noise. We consider three functional forms for position dependence of the diffusion coefficient: power-law, exponential, and logarithmic. The coloured non-Gaussian noise is distributed according to Tsallis' q-distribution. Tracks of the non-Markovian systems are numerically simulated by using the fourth-order Runge-Kutta algorithm and the first passage times (FPTs) are recorded. The FPT density is determined along with the mean FPT (MFPT). Effects of the noise intensity and self-correlation of the multiplicative noise, the intensity of the additive noise, the cross-correlation strength, and the non-extensivity parameter on the MFPT are discussed.
How different microscopic mechanisms of ultrafast spin dynamics coexist and interplay is not only relevant for the development of spintronics but also for the thorough description of physical systems out-of-equilibrium. In pure crystalline ferromagnets, one of the main microscopic mechanism of spin relaxation is the electron-phonon (el-ph) driven spin-flip, or Elliott-Yafet, scattering. Unexpectedly, recent experiments with ferro- and ferrimagnetic alloys have shown different dynamics for the different sublattices. These distinct sublattice dynamics are contradictory to the Elliott-Yafet scenario. In order to rationalize this discrepancy, it has been proposed that the intra- and intersublattice exchange interaction energies must be considered in the microscopic demagnetization mechanism, too. Here, using a temperature-dependent x-ray emission spectroscopy (XES) method, we address experimentally the element specific el-ph angular momentum transfer rates, responsible for the spin-flips in the respective (sub)lattices of Fe20Ni80, Fe50Ni50 and pure nickel single crystals. We establish how the deduced rate evolution with the temperature is linked to the exchange coupling constants reported for different alloy stoichiometries and how sublattice exchange energies threshold the related el-ph spin-flip channels. Thus, these results evidence that the Elliott-Yafet spin-flip scattering, thresholded by sublattice exchange energies, is the relevant microscopic process to describe sublattice dynamics in alloys and elemental magnetic systems.
We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary.
In their comment on our paper (Caesar et al 2020 Environ. Res. Lett. 15 024003), Chen and Tung (hereafter C&T) argue that our analysis, showing that over the last decades Atlantic meridional overturning circulation (AMOC) strength and global mean surface temperature (GMST) were positively correlated, is incorrect. Their claim is mainly based on two arguments, neither of which is justified: first, C&T claim that our analysis is based on 'established evidence' that was only true for preindustrial conditions-this is not the case. Using data from the modern period (1947-2012), we show that the established understanding (i.e. deep-water formation in the North Atlantic cools the deep ocean and warms the surface) is correct, but our analysis is not based on this fact. Secondly, C&T claim that our results are based on a statistical analysis of only one cycle of data which was furthermore incorrectly detrended. This, too, is not true. Our conclusion that a weaker AMOC delays the current surface warming rather than enhances it, is based on several independent lines of evidence. The data we show to support this covers more than one cycle and the detrending (which was performed to avoid spurious correlations due to a common trend) does not affect our conclusion: the correlation between AMOC strength and GMST is positive. We do not claim that this is strong evidence that the two time series are in phase, but rather that this means that the two time series are not anti-correlated.
Analysis of protrusion dynamics in amoeboid cell motility by means of regularized contour flows
(2021)
Amoeboid cell motility is essential for a wide range of biological processes including wound healing, embryonic morphogenesis, and cancer metastasis. It relies on complex dynamical patterns of cell shape changes that pose long-standing challenges to mathematical modeling and raise a need for automated and reproducible approaches to extract quantitative morphological features from image sequences. Here, we introduce a theoretical framework and a computational method for obtaining smooth representations of the spatiotemporal contour dynamics from stacks of segmented microscopy images. Based on a Gaussian process regression we propose a one-parameter family of regularized contour flows that allows us to continuously track reference points (virtual markers) between successive cell contours. We use this approach to define a coordinate system on the moving cell boundary and to represent different local geometric quantities in this frame of reference. In particular, we introduce the local marker dispersion as a measure to identify localized membrane expansions and provide a fully automated way to extract the properties of such expansions, including their area and growth time. The methods are available as an open-source software package called AmoePy, a Python-based toolbox for analyzing amoeboid cell motility (based on time-lapse microscopy data), including a graphical user interface and detailed documentation. Due to the mathematical rigor of our framework, we envision it to be of use for the development of novel cell motility models. We mainly use experimental data of the social amoeba Dictyostelium discoideum to illustrate and validate our approach. <br /> Author summary Amoeboid motion is a crawling-like cell migration that plays an important key role in multiple biological processes such as wound healing and cancer metastasis. This type of cell motility results from expanding and simultaneously contracting parts of the cell membrane. From fluorescence images, we obtain a sequence of points, representing the cell membrane, for each time step. By using regression analysis on these sequences, we derive smooth representations, so-called contours, of the membrane. Since the number of measurements is discrete and often limited, the question is raised of how to link consecutive contours with each other. In this work, we present a novel mathematical framework in which these links are described by regularized flows allowing a certain degree of concentration or stretching of neighboring reference points on the same contour. This stretching rate, the so-called local dispersion, is used to identify expansions and contractions of the cell membrane providing a fully automated way of extracting properties of these cell shape changes. We applied our methods to time-lapse microscopy data of the social amoeba Dictyostelium discoideum.
We performed numerical simulations with the Kuramoto model and experiments with oscillatory nickel electrodissolution to explore the dynamical features of the transients from random initial conditions to a fully synchronized (one-cluster) state. The numerical simulations revealed that certain networks (e.g., globally coupled or dense Erdos-Renyi random networks) showed relatively simple behavior with monotonic increase of the Kuramoto order parameter from the random initial condition to the fully synchronized state and that the transient times exhibited a unimodal distribution. However, some modular networks with bridge elements were identified which exhibited non-monotonic variation of the order parameter with local maximum and/or minimum. In these networks, the histogram of the transients times became bimodal and the mean transient time scaled well with inverse of the magnitude of the second largest eigenvalue of the network Laplacian matrix. The non-monotonic transients increase the relative standard deviations from about 0.3 to 0.5, i.e., the transient times became more diverse. The non-monotonic transients are related to generation of phase patterns where the modules are synchronized but approximately anti-phase to each other. The predictions of the numerical simulations were demonstrated in a population of coupled oscillatory electrochemical reactions in global, modular, and irregular tree networks. The findings clarify the role of network structure in generation of complex transients that can, for example, play a role in intermittent desynchronization of the circadian clock due to external cues or in deep brain stimulations where long transients are required after a desynchronization stimulus.
Additive manufacturing (AM) of metals and in particular laser powder bed fusion (LPBF) enables a degree of freedom in design unparalleled by conventional subtractive methods. To ensure that the designed precision is matched by the produced LPBF parts, a full understanding of the interaction between the laser and the feedstock powder is needed. It has been shown that the laser also melts subjacent layers of material underneath. This effect plays a key role when designing small cavities or overhanging structures, because, in these cases, the material underneath is feed-stock powder. In this study, we quantify the extension of the melt pool during laser illumination of powder layers and the defect spatial distribution in a cylindrical specimen. During the LPBF process, several layers were intentionally not exposed to the laser beam at various locations, while the build process was monitored by thermography and optical tomography. The cylinder was finally scanned by X-ray computed tomography (XCT). To correlate the positions of the unmolten layers in the part, a staircase was manufactured around the cylinder for easier registration. The results show that healing among layers occurs if a scan strategy is applied, where the orientation of the hatches is changed for each subsequent layer. They also show that small pores and surface roughness of solidified material below a thick layer of unmolten material (>200 mu m) serve as seeding points for larger voids. The orientation of the first two layers fully exposed after a thick layer of unmolten powder shapes the orientation of these voids, created by a lack of fusion.
Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea level rise. Increased surface melt may lower the reflectivity of the ice sheet surface and thereby increase melt rates: the so-called melt-albedo feedback describes this self-sustaining increase in surface melting. In order to test the effect of the melt-albedo feedback in a prognostic ice sheet model, we implement dEBM-simple, a simplified version of the diurnal Energy Balance Model dEBM, in the Parallel Ice Sheet Model (PISM). The implementation includes a simple representation of the melt-albedo feedback and can thereby replace the positive-degree-day melt scheme. Using PISM-dEBM-simple, we find that this feedback increases ice loss through surface warming by 60 % until 2300 for the high-emission scenario RCP8.5 when compared to a scenario in which the albedo remains constant at its present-day values. With an increase of 90 % compared to a fixed-albedo scenario, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be 70 % in the RCP8.5 scenario and a more than 4-fold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt-albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming.
Based on suggested interactions of potential tipping elements in the Earth's climate and in ecological systems, tipping cascades as possible dynamics are increasingly discussed and studied. The activation of such tipping cascades would impose a considerable risk for human societies and biosphere integrity. However, there are ambiguities in the description of tipping cascades within the literature so far. Here we illustrate how different patterns of multiple tipping dynamics emerge from a very simple coupling of two previously studied idealized tipping elements. In particular, we distinguish between a two phase cascade, a domino cascade and a joint cascade. A mitigation of an unfolding two phase cascade may be possible and common early warning indicators are sensitive to upcoming critical transitions to a certain degree. In contrast, a domino cascade may hardly be stopped once initiated and critical slowing down-based indicators fail to indicate tipping of the following element. These different potentials for intervention and anticipation across the distinct patterns of multiple tipping dynamics should be seen as a call to be more precise in future analyses of cascading dynamics arising from tipping element interactions in the Earth system.
The most complex but potentially most severe impacts of climate change are caused by extreme weather events. In a globally connected economy, damages can cause remote perturbations and cascading consequences-a ripple effect along supply chains. Here we show an economic ripple resonance that amplifies losses when consecutive or overlapping weather extremes and their repercussions interact. This amounts to an average amplification of 21% for climate-induced heat stress, river floods, and tropical cyclones. Modeling the temporal evolution of 1.8 million trade relations between >7000 regional economic sectors, we find that the regional responses to future extremes are strongly heterogeneous also in their resonance behavior. The induced effect on welfare varies between gains due to increased demand in some regions and losses due to demand or supply shortages in others. Within the current global supply network, the ripple resonance effect of extreme weather is strongest in high-income economies-an important effect to consider when evaluating past and future economic climate impacts.
Complex networks are abundant in nature and many share an important structural property: they contain a few nodes that are abnormally highly connected (hubs). Some of these hubs are called influencers because they couple strongly to the network and play fundamental dynamical and structural roles. Strikingly, despite the abundance of networks with influencers, little is known about their response to stochastic forcing. Here, for oscillatory dynamics on influencer networks, we show that subjecting influencers to an optimal intensity of noise can result in enhanced network synchronization. This new network dynamical effect, which we call coherence resonance in influencer networks, emerges from a synergy between network structure and stochasticity and is highly nonlinear, vanishing when the noise is too weak or too strong. Our results reveal that the influencer backbone can sharply increase the dynamical response in complex systems of coupled oscillators. Influencer networks include a small set of highly-connected nodes and can reach synchrony only via strong node interaction. Tonjes et al. show that introducing an optimal amount of noise enhances synchronization of such networks, which may be relevant for neuroscience or opinion dynamics applications.
Tropical cyclones range among the costliest disasters on Earth. Their economic repercussions along the supply and trade network also affect remote economies that are not directly affected. We here simulate possible global repercussions on consumption for the example case of Hurricane Sandy in the US (2012) using the shock-propagation model Acclimate. The modeled shock yields a global three-phase ripple: an initial production demand reduction and associated consumption price decrease, followed by a supply shortage with increasing prices, and finally a recovery phase. Regions with strong trade relations to the US experience strong magnitudes of the ripple. A dominating demand reduction or supply shortage leads to overall consumption gains or losses of a region, respectively. While finding these repercussions in historic data is challenging due to strong volatility of economic interactions, numerical models like ours can help to identify them by approaching the problem from an exploratory angle, isolating the effect of interest. For this, our model simulates the economic interactions of over 7000 regional economic sectors, interlinked through about 1.8 million trade relations. Under global warming, the wave-like structures of the economic response to major hurricanes like the one simulated here are likely to intensify and potentially overlap with other weather extremes.
Due to climate change the frequency and character of precipitation are changing as the hydrological cycle intensifies. With regards to snowfall, global warming has two opposing influences; increasing humidity enables intense snowfall, whereas higher temperatures decrease the likelihood of snowfall. Here we show an intensification of extreme snowfall across large areas of the Northern Hemisphere under future warming. This is robust across an ensemble of global climate models when they are bias-corrected with observational data. While mean daily snowfall decreases, both the 99th and the 99.9th percentiles of daily snowfall increase in many regions in the next decades, especially for Northern America and Asia. Additionally, the average intensity of snowfall events exceeding these percentiles as experienced historically increases in many regions. This is likely to pose a challenge to municipalities in mid to high latitudes. Overall, extreme snowfall events are likely to become an increasingly important impact of climate change in the next decades, even if they will become rarer, but not necessarily less intense, in the second half of the century.
Laser-based additive manufacturing methods allow the production of complex metal structures within a single manufacturing step. However, the localized heat input and the layer-wise manufacturing manner give rise to large thermal gradients. Therefore, large internal stress (IS) during the process (and consequently residual stress (RS) at the end of production) is generated within the parts. This IS or RS can either lead to distortion or cracking during fabrication or in-service part failure, respectively. With this in view, the knowledge on the magnitude and spatial distribution of RS is important to develop strategies for its mitigation. Specifically, diffraction-based methods allow the spatial resolved determination of RS in a non-destructive fashion. In this review, common diffraction-based methods to determine RS in laser-based additive manufactured parts are presented. In fact, the unique microstructures and textures associated to laser-based additive manufacturing processes pose metrological challenges. Based on the literature review, it is recommended to (a) use mechanically relaxed samples measured in several orientations as appropriate strain-free lattice spacing, instead of powder, (b) consider that an appropriate grain-interaction model to calculate diffraction-elastic constants is both material- and texture-dependent and may differ from the conventionally manufactured variant. Further metrological challenges are critically reviewed and future demands in this research field are discussed.
Increasing greenhouse gas emissions are likely to impact not only natural systems but economies worldwide. If these impacts alter future economic development, the financial losses will be significantly higher than the mere direct damages. So far, potentially aggravating investment responses were considered negligible. Here we consistently incorporate an empirically derived temperature-growth relation into the simple integrated assessment model DICE. In this framework we show that, if in the next eight decades varying temperatures impact economic growth as has been observed in the past three decades, income is reduced by similar to 20% compared to an economy unaffected by climate change. Hereof similar to 40% are losses due to growth effects of which similar to 50% result from reduced incentive to invest. This additional income loss arises from a reduced incentive for future investment in anticipation of a reduced return and not from an explicit climate protection policy. Under economically optimal climate-change mitigation, however, optimal investment would only be reduced marginally as mitigation efforts keep returns high.
Both ice sheets in Greenland and Antarctica are discharging ice into the ocean. In many regions along the coast of the ice sheets, the icebergs calve into a bay. If the addition of icebergs through calving is faster than their transport out of the embayment, the icebergs will be frozen into a melange with surrounding sea ice in winter. In this case, the buttressing effect of the ice melange can be considerably stronger than any buttressing by mere sea ice would be. This in turn stabilizes the glacier terminus and leads to a reduction in calving rates. Here we propose a simple parametrization of ice melange buttressing which leads to an upper bound on calving rates and can be used in numerical and analytical modelling.
Topic and aim. Synchronization in populations of coupled oscillators can be characterized with order parameters that describe collective order in ensembles. A dependence of the order parameter on the coupling constants is well-known for coupled periodic oscillators. The goal of the study is to extend this analysis to ensembles of oscillators with chaotic phases, moreover with phases possessing hyperbolic chaos. Models and methods. Two models are studied in the paper. One is an abstract discrete-time map, composed with a hyperbolic Bernoulli transformation and with Kuramoto dynamics. Another model is a system of coupled continuous-time chaotic oscillators, where each individual oscillator has a hyperbolic attractor of Smale-Williams type. Results. The discrete-time model is studied with the Ott-Antonsen ansatz, which is shown to be invariant under the application of the Bernoulli map. The analysis of the resulting map for the order parameter shows, that the asynchronouis state is always stable, but the synchronous one becomes stable above a certain coupling strength. Numerical analysis of the continuous-time model reveals a complex sequence of transitions from an asynchronous state to a completely synchronous hyperbolic chaos, with intermediate stages that include regimes with periodic in time mean field, as well as with weakly and strongly irregular mean field variations. Discussion. Results demonstrate that synchronization of systems with hyperbolic chaos of phases is possible, although a rather strong coupling is required. The approach can be applied to other systems of interacting units with hyperbolic chaotic dynamics.
A detailed investigation of the energy levels of perylene-3,4,9,10-tetracarboxylic tetraethylester as a representative compound for the whole family of perylene esters was performed. It was revealed via electrochemical measurements that one oxidation and two reductions take place. The bandgaps determined via the electrochemical approach are in good agreement with the optical bandgap obtained from the absorption spectra via a Tauc plot. In addition, absorption spectra in dependence of the electrochemical potential were the basis for extensive quantum-chemical calculations of the neutral, monoanionic, and dianionic molecules. For this purpose, calculations based on density functional theory were compared with post-Hartree-Fock methods and the CAM-B3LYP functional proved to be the most reliable choice for the calculation of absorption spectra. Furthermore, spectral features found experimentally could be reproduced with vibronic calculations and allowed to understand their origins. In particular, the two lowest energy absorption bands of the anion are not caused by absorption of two distinct electronic states, which might have been expected from vertical excitation calculations, but both states exhibit a strong vibronic progression resulting in contributions to both bands.
Accurate and precise characterization of cirrus cloud geometrical and optical properties is essential for better constraining their radiative footprint. A lidar-based retrieval scheme is proposed here, with its performance assessed on fine spatio-temporal observations over the Arctic site of Ny-Alesund, Svalbard. Two contributions related to cirrus geometrical (dynamic Wavelet Covariance Transform (WCT)) and optical properties (constrained Klett) are reported. The dynamic WCT rendered cirrus detection more robust, especially for thin cirrus layers that frequently remained undetected by the classical WCT method. Regarding optical characterization, we developed an iterative scheme for determining the cirrus lidar ratio (LRci) that is a crucial parameter for aerosol - cloud discrimination. Building upon the Klett-Fernald method, the LRci was constrained by an additional reference value. In established methods, such as the double-ended Klett, an aerosol-free reference value is applied. In the proposed constrained Klett, however, the reference value was approximated from cloud-free or low cloud optical depth (COD up to 0.2) profiles and proved to agree with independent Raman estimates. For optically thin cirrus, the constrained Klett inherent uncertainties reached 50% (60-74%) in terms of COD (LRci). However, for opaque cirrus COD (LRci) uncertainties were lower than 10% (15%). The detection method discrepancies (dynamic versus static WCT) had a higher impact on the optical properties of low COD layers (up to 90%) compared to optically thicker ones (less than 10%). The constrained Klett presented high agreement with two established retrievals. For an exemplary cirrus cloud, the constrained Klett estimated the COD355 (LRci355) at 0.28 +/- 0.17 (29 +/- 4 sr), the double-ended Klett at 0.27 +/- 0.15 (32 +/- 4 sr) and the Raman retrievals at 0.22 +/- 0.12 (26 +/- 11 sr). Our approach to determine the necessary reference value can also be applied in established methods and increase their accuracy. In contrast, the classical aerosol-free assumption led to 44 sr LRci overestimation in optically thin layers and 2-8 sr in thicker ones. The multiple scattering effect was corrected using Eloranta (1998) and accounted for 50-60% extinction underestimation near the cloud base and 20-30% within the cirrus layers.
Cirrus is the only cloud type capable of inducing daytime cooling or heating at the top of the atmosphere (TOA) and the sign of its radiative effect highly depends on its optical depth. However, the investigation of its geometrical and optical properties over the Arctic is limited. In this work the long-term properties of cirrus clouds are explored for the first time over an Arctic site (Ny-Alesund, Svalbard) using lidar and radiosonde measurements from 2011 to 2020. The optical properties were quality assured, taking into account the effects of specular reflections and multiple-scattering. Cirrus clouds were generally associated with colder and calmer wind conditions compared to the 2011-2020 climatology. However, the dependence of cirrus properties on temperature and wind speed was not strong. Even though the seasonal cycle was not pronounced, the winter-time cirrus appeared under lower temperatures and stronger wind conditions. Moreover, in winter, geometrically- and optically-thicker cirrus were found and their ice particles tended to be more spherical. The majority of cirrus was associated with westerly flow and westerly cirrus tended to be geometrically-thicker. Overall, optically-thinner layers tended to comprise smaller and less spherical ice crystals, most likely due to reduced water vapor deposition on the particle surface. Compared to lower latitudes, the cirrus layers over Ny-Alesund were more absorbing in the visible spectral region and they consisted of more spherical ice particles.
The Ornstein–Uhlenbeck process is a stationary and ergodic Gaussian process, that is fully determined by its covariance function and mean. We show here that the generic definitions of the ensemble- and time-averaged mean squared displacements fail to capture these properties consistently, leading to a spurious ergodicity breaking. We propose to remedy this failure by redefining the mean squared displacements such that they reflect unambiguously the statistical properties of any stochastic process. In particular we study the effect of the initial condition in the Ornstein–Uhlenbeck process and its fractional extension. For the fractional Ornstein–Uhlenbeck process representing typical experimental situations in crowded environments such as living biological cells, we show that the stationarity of the process delicately depends on the initial condition.
Electric currents flowing in the terrestrial ionosphere have conventionally been diagnosed by low-earth-orbit (LEO) satellites equipped with science-grade magnetometers and long booms on magnetically clean satellites. In recent years, there are a variety of endeavors to incorporate platform magnetometers, which are initially designed for navigation purposes, to study ionospheric currents. Because of the suboptimal resolution and significant noise of the platform magnetometers, however, most of the studies were confined to high-latitude auroral regions, where magnetic field deflections from ionospheric currents easily exceed 100 nT. This study aims to demonstrate the possibility of diagnosing weak low-/mid-latitude ionospheric currents based on platform magnetometers. We use navigation magnetometer data from two satellites, CryoSat-2 and the Gravity Recovery and Climate Experiment Follow-On (GRACE-FO), both of which have been intensively calibrated based on housekeeping data and a high-precision geomagnetic field model. Analyses based on 8 years of CryoSat-2 data as well as similar to 1.5 years of GRACE-FO data reproduce well-known climatology of inter-hemispheric field-aligned currents (IHFACs), as reported by previous satellite missions dedicated to precise magnetic observations. Also, our results show that C-shaped structures appearing in noontime IHFAC distributions conform to the shape of the South Atlantic Anomaly. The F-region dynamo currents are only partially identified in the platform magnetometer data, possibly because the currents are weaker than IHFACs in general and depend significantly on altitude and solar activity. Still, this study evidences noontime F-region dynamo currents at the highest altitude (717 km) ever reported. We expect that further data accumulation from continuously operating missions may reveal the dynamo currents more clearly during the next solar maximum.
We report on the detection of very high energy (VHE; E > 100 GeV) gamma-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multiwavelength observations with Fermi/LAT, XRT and UVOT onboard the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE gamma-ray regime, we deduce a 95 per cent confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98 and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results, the redshift of KUV 00311-1938 is constrained to 0.51 <= z < 0.98 and of PKS 1440-389 to 0.14 (sic) z < 0.53.
We study the experimentally measured ciprofloxacin antibiotic diffusion through a gel-like artificial sputum medium (ASM) mimicking physiological conditions typical for a cystic fibrosis layer, in which regions occupied by Pseudomonas aeruginosa bacteria are present. To quantify the antibiotic diffusion dynamics we employ a phenomenological model using a subdiffusion-absorption equation with a fractional time derivative. This effective equation describes molecular diffusion in a medium structured akin Thompson’s plumpudding model; here the ‘pudding’ background represents the ASM and the ‘plums’ represent the bacterial biofilm. The pudding is a subdiffusion barrier for antibiotic molecules that can affect bacteria found in plums. For the experimental study we use an interferometric method to determine the time evolution of the amount of antibiotic that has diffused through the biofilm. The theoretical model shows that this function is qualitatively different depending on whether or not absorption of the antibiotic in the biofilm occurs. We show that the process can be divided into three successive stages: (1) only antibiotic subdiffusion with constant biofilm parameters, (2) subdiffusion and absorption of antibiotic molecules with variable biofilm transport parameters, (3) subdiffusion and absorption in the medium but the biofilm parameters are constant again. Stage 2 is interpreted as the appearance of an intensive defence build–up of bacteria against the action of the antibiotic, and in the stage 3 it is likely that the bacteria have been inactivated. Times at which stages change are determined from the experimentally obtained temporal evolution of the amount of antibiotic that has diffused through the ASM with bacteria. Our analysis shows good agreement between experimental and theoretical results and is consistent with the biologically expected biofilm response. We show that an experimental method to study the temporal evolution of the amount of a substance that has diffused through a biofilm is useful in studying the processes occurring in a biofilm. We also show that the complicated biological process of antibiotic diffusion in a biofilm can be described by a fractional subdiffusion-absorption equation with subdiffusion and absorption parameters that change over time.
We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed.
During the last decade, intracellular actin waves have attracted much attention due to their essential role in various cellular functions, ranging from motility to cytokinesis. Experimental methods have advanced significantly and can capture the dynamics of actin waves over a large range of spatio-temporal scales. However, the corresponding coarse-grained theory mostly avoids the full complexity of this multi-scale phenomenon. In this perspective, we focus on a minimal continuum model of activator–inhibitor type and highlight the qualitative role of mass conservation, which is typically overlooked. Specifically, our interest is to connect between the mathematical mechanisms of pattern formation in the presence of a large-scale mode, due to mass conservation, and distinct behaviors of actin waves.
Whereas self-propelled hard discs undergo motility-induced phase separation, self-propelled rods exhibit a variety of nonequilibrium phenomena, including clustering, collective motion, and spatio-temporal chaos. In this work, we present a theoretical framework representing active particles by continuum fields. This concept combines the simplicity of alignment-based models, enabling analytical studies, and realistic models that incorporate the shape of self-propelled objects explicitly. By varying particle shape from circular to ellipsoidal, we show how nonequilibrium stresses acting among self-propelled rods destabilize motility-induced phase separation and facilitate orientational ordering, thereby connecting the realms of scalar and vectorial active matter. Though the interaction potential is strictly apolar, both, polar and nematic order may emerge and even coexist. Accordingly, the symmetry of ordered states is a dynamical property in active matter. The presented framework may represent various systems including bacterial colonies, cytoskeletal extracts, or shaken granular media. Interacting self-propelled particles exhibit phase separation or collective motion depending on particle shape. A unified theory connecting these paradigms represents a major challenge in active matter, which the authors address here by modeling active particles as continuum fields.