530 Physik
Refine
Year of publication
Document Type
- Article (950) (remove)
Keywords
- diffusion (25)
- anomalous diffusion (22)
- gamma rays: general (20)
- ISM: supernova remnants (12)
- cosmic rays (12)
- stars: massive (11)
- organic solar cells (9)
- stars: atmospheres (9)
- astroparticle physics (8)
- perovskite solar cells (8)
Institute
- Institut für Physik und Astronomie (874)
- Institut für Chemie (36)
- Institut für Mathematik (15)
- Extern (8)
- Institut für Geowissenschaften (6)
- Institut für Umweltwissenschaften und Geographie (6)
- Mathematisch-Naturwissenschaftliche Fakultät (5)
- Fachgruppe Politik- & Verwaltungswissenschaft (4)
- Hasso-Plattner-Institut für Digital Engineering gGmbH (3)
- Institut für Biochemie und Biologie (3)
We report a comparison of two photonic techniques for single-molecule sensing: fluorescence nanoscopy and optoplasmonic sensing. As the test system, oligonucleotides with and without fluorescent labels are transiently hybridized to complementary "docking" strands attached to gold nanorods. Comparing the measured single-molecule kinetics helps to examine the influence of the fluorescent labels as well as factors arising from different sensing geometries. Our results demonstrate that DNA dissociation is not significantly altered by the fluorescent labels and that DNA association is affected by geometric factors in the two techniques. These findings open the door to exploiting plasmonic sensing and fluorescence nanoscopy in a complementary fashion, which will aid in building more powerful sensors and uncovering the intricate effects that influence the behavior of single molecules.
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.
A symmetry-breaking mechanism is investigated that creates bistability between fully and partially synchronized states in oscillator networks. Two populations of oscillators with unimodal frequency distribution and different amplitudes, in the presence of weak global coupling, are shown to simplify to a modular network with asymmetrical coupling. With increasing the coupling strength, a synchronization transition is observed with an isolated fully synchronized state. The results are interpreted theoretically in the thermodynamic limit and confirmed in experiments with chemical oscillators.
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.
The electron-phonon scattering is one of the main microscopic mechanisms responsible for the spin-flip in the transient state of ultrafast demagnetization. Here, we present an experimental determination of the temperature-dependent electron-phonon scattering rate in Gd. Using a static x-ray emission spectroscopy method, where the reduction of the decay peak intensities when increasing the temperature is quantified, we measure independently the electron-phonon scattering rate for the 5d and the 4f electrons. We deduce the temperature dependence of scattering for the 5d electrons, while no effect on the phonon population is observed for the 4f electrons. Our results suggest that the ultrafast magnetization dynamics in Gd is triggered by the spin-flip in the 5d electrons. We also evidence the existence of a temperature threshold, above which spin-flip scattering of the 5d electrons takes place. We deduce that during the transient state of ultrafast demagnetization, the exchange energy between 5d electrons has to be overcome before the microscopic electron-phonon scattering process can occur.
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.
Steric hindrance of hydration and hydrogen bond enhancement by localized charges have been identified as key factors for the massive chemical differences between the hydroxypyridine/pyridone isomers in aqueous solution. While all isomers occur mainly in the hydroxypyridine form in the gas phase, they differ by more than 3 orders of magnitude both in their acidity and tautomeric equilibrium constants upon hydration. By monitoring the electronic and solvation structures as a function of the protonation state and the O- substitution position on the pyridine ring, the amplification of the isomeric differences in aqueous solution has been investigated. Near-edge X-ray absorption fine structure (NEXAFS) measurements at the N K-edge served as the probe of the chemical state. The combination of molecular dynamics simulations, complete active space self-consistent field (CASSCF), and time-dependent density functional theory (TD-DFT) spectral calculations contributes to unraveling the principles of tautomerism and acidity in multiple biochemical systems based on tautomerism.
Laser based powder bed fusion additive manufacturing offers the flexibility to incorporate standard and user-defined scan strategies in a layer or in between the layers for the customized fabrication of metallic components. In the present study, four different scan strategies and their impact on the development of microstructure, texture, and residual stresses in laser powder bed fusion additive manufacturing of a nickel-based superalloy Inconel 718 was investigated. Light microscopy, scanning electron microscopy combined with electron back-scatter diffraction, and neutron diffraction were used as the characterization tools. Strong textures with epitaxially grown columnar grains were observed along the build direction for the two individual scan strategies. Patterns depicting the respective scan strategies were visible in the build plane, which dictated the microstructure development in the other planes. An alternating strategy combining the individual strategies in the successive layers and a 67 degrees rotational strategy weakened the texture by forming finer micro-structural features. Von Mises equivalent stress plots revealed lower stress values and gradients, which translates as lower distortions for the alternating and rotational strategies. Overall results confirmed the scope for manipulating the microstructure, texture, and residual stresses during laser powder bed fusion additive manufacturing by effectively controlling the scan strategies.
Recent experiments on laser-dissociation of aligned homonuclear diatomic molecules show an asymmetric forward-backward (spatial) electron-localization along the laser polarization axis. Most theoretical models attribute this asymmetry to interference effects between gerade and ungerade vibronic states. Presumably due to alignment, these models neglect molecular rotations and hence infer an asymmetric (post-dissociation) charge distribution over the two identical nuclei. In this paper, we question the equivalence that is made between spatial electron-localization, observed in experiments, and atomic electron-localization, alluded by these theoretical models. We show that (seeming) agreement between these models and experiments is due to an unfortunate omission of nuclear permutation symmetry, i.e., quantum statistics. Enforcement of the latter requires mandatory inclusion of the molecular rotational degree of freedom, even for perfectly aligned molecules. Unlike previous interpretations, we ascribe spatial electron-localization to the laser creation of a rovibronic wavepacket that involves field-free molecular eigenstates with opposite space-inversion symmetry i.e., even and odd parity. Space-inversion symmetry breaking would then lead to an asymmetric distribution of the (space-fixed) electronic density over the forward and backward hemisphere. However, owing to the simultaneous coexistence of two indistinguishable molecular orientational isomers, our analytical and computational results show that the post-dissociation electronic density along a specified space-fixed axis is equally shared between the two identical nuclei-a result that is in perfect accordance with the principle of the indistinguishability of identical particles. Published under an exclusive license by AIP Publishing.
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.
Hydrogels with a hierarchical structure were prepared from a new highly water-soluble crosslinker N,N,N',N'-tetramethyl-N,N'-bis(2-ethylmethacrylate)-propyl-1,3-diammonium dibromide and from the sulfobetaine monomer 2-(N-3-sulfopropyl-N,N-dimethyl ammonium)ethyl methacrylate. The free radical polymerization of the two compounds is rapid and yields near-transparent hydrogels with sizes up to 5 cm in diameter. Rheology shows a clear correlation between the monomer-to-crosslinker ratio and the storage and loss moduli of the hydrogels. Cryo-scanning electron microscopy, low-field nuclear magnetic resonance (NMR) spectroscopy, and small-angle X-ray scattering show that the gels have a hierarchical structure with features spanning the nanometer to the sub-millimeter scale. The NMR study is challenged by the marked inhomogeneity of the gels and the complex chemical structure of the sulfobetaine monomer. NMR spectroscopy shows how these complications can be addressed via a novel fitting approach that considers the mobility gradient along the side chain of methacrylate-based monomers.
Self-organized coherence-incoherence patterns, called chimera states, have first been reported in systems of Kuramoto oscillators. For coupled excitable units, similar patterns where coherent units are at rest are called bump states. Here, we study bumps in an array of active rotators coupled by nonlocal attraction and global repulsion. We demonstrate how they can emerge in a supercritical scenario from completely coherent Turing patterns: a single incoherent unit appears in a homoclinic bifurcation, undergoing subsequent transitions to quasiperiodic and chaotic behavior, which eventually transforms into extensive chaos with many incoherent units. We present different types of transitions and explain the formation of coherence-incoherence patterns according to the classical paradigm of short-range activation and long-range inhibition.
Various mathematical Black-Scholes-Merton-like models of option pricing employ the paradigmatic stochastic process of geometric Brownian motion (GBM). The innate property of such models and of real stock-market prices is the roughly exponential growth of prices with time [on average, in crisis-free times]. We here explore the ensemble- and time averages of a multiplicative-noise stochastic process with power-law-like time-dependent volatility, sigma(t) similar to t(alpha), named scaled GBM (SGBM). For SGBM, the mean-squared displacement (MSD) computed for an ensemble of statistically equivalent trajectories can grow faster than exponentially in time, while the time-averaged MSD (TAMSD)-based on a sliding-window averaging along a single trajectory-is always linear at short lag times Delta. The proportionality factor between these the two averages of the time series is Delta/T at short lag times, where T is the trajectory length, similarly to GBM. This discrepancy of the scaling relations and pronounced nonequivalence of the MSD and TAMSD at Delta/T << 1 is a manifestation of weak ergodicity breaking for standard GBM and for SGBM with s (t)-modulation, the main focus of our analysis. The analytical predictions for the MSD and mean TAMSD for SGBM are in quantitative agreement with the results of stochastic computer simulations.
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.
Moving spiral wave chimeras
(2021)
We consider a two-dimensional array of heterogeneous nonlocally coupled phase oscillators on a flat torus and study the bound states of two counter-rotating spiral chimeras, shortly two-core spiral chimeras, observed in this system. In contrast to other known spiral chimeras with motionless incoherent cores, the two-core spiral chimeras typically show a drift motion. Due to this drift, their incoherent cores become spatially modulated and develop specific fingerprint patterns of varying synchrony levels. In the continuum limit of infinitely many oscillators, the two-core spiral chimeras can be studied using the Ott-Antonsen equation. Numerical analysis of this equation allows us to reveal the stability region of different spiral chimeras, which we group into three main classes-symmetric, asymmetric, and meandering spiral chimeras.
Generalized (non-Markovian) diffusion equations with different memory kernels and subordination schemes based on random time change in the Brownian diffusion process are popular mathematical tools for description of a variety of non-Fickian diffusion processes in physics, biology, and earth sciences. Some of such processes (notably, the fluid limits of continuous time random walks) allow for either kind of description, but other ones do not. In the present work we discuss the conditions under which a generalized diffusion equation does correspond to a subordination scheme, and the conditions under which a subordination scheme does possess the corresponding generalized diffusion equation. Moreover, we discuss examples of random processes for which only one, or both kinds of description are applicable.
We investigate whether the dust content of the circum-galactic medium (CGM) depends on the location of the quasar sightline with respect to the galaxy major-axis using 13 galaxy-Mg II absorber pairs (9-81 kpc distance) from the MusE GAs FLOw and Wind (MEGAFLOW) survey at 0.4 < z < 1.4. The dust content of the CGM is obtained from [Zn/Fe] using ultraviolet and visual echelle spectrograph data. When a direct measurement of [Zn/Fe] is unavailable, we estimate the dust depletion from a method that consists in solving for the depletion from multiple singly ionized ions (e.g. Mn II, Cr II, and Zn II) since each ion depletes on dust grains at different rates. We find a positive correlation between the azimuthal angle and [Zn/Fe] with a Pearson's gamma = 0.70 +/- 0.14. The sightlines along the major axis show [Zn/Fe] < 0.5, whereas the [Zn/Fe] is > 0.8 along the minor axis. These results suggest that the CGM along the minor axis is on average more metal enriched (by approximate to 1 dex) than the gas located along the major axis of galaxies provided that dust depletion is a proxy for metallicity. This anisotropic distribution is consistent with recent results on outflow and accretion in hydro-dynamical simulations.
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.
Proximity to heavy sp-elements is considered promising for reaching a band gap in graphene that could host quantum spin Hall states. The recent report of an induced spin-orbit gap of 0.2 eV in Pb-intercalated graphene detectable by spin-resolved photoemission has spurred renewed interest in such systems (Klimovskikh et al 2017 ACS Nano 11, 368). In the case of Bi intercalation an even larger band gap of 0.4 eV has been observed but was assigned to the influence of a dislocation network (Warmuth et al 2016 Phys. Rev. B 93, 165 437). Here, we study Bi intercalation under graphene on Ir(111) and report a nearly ideal graphene dispersion without band replicas and no indication of hybridization with the substrate. The band gap is small (0.19 eV) and can be tuned by +/- 25 meV through the Bi coverage. The Bi atomic density is higher than in the recent report. By spin-resolved photoemission we exclude induced spin-orbit interaction as origin of the gap. Quantitative agreement of a photoemission intensity analysis with the measured band gap suggests sublattice symmetry breaking as one of the possible band gap opening mechanisms. We test several Bi structures by density functional theory. Our results indicate the possibility that Bi intercalates in the phase of bismuthene forming a graphene-bismuthene van der Waals heterostructure.
Simple and robust
(2021)
A spectrum of 7562 publications on Molecularly Imprinted Polymers (MIPs) has been presented in literature within the last ten years (Scopus, September 7, 2020). Around 10 % of the papers published on MIPs describe the recognition of proteins. The straightforward synthesis of MIPs is a significant advantage as compared with the preparation of enzymes or antibodies. MIPs have been synthesized from only one up to six functional monomers while proteins are made up of 20 natural amino acids. Furthermore, they can be synthesized against structures of low immunogenicity and allow multi-analyte measurements via multi-target synthesis. Electrochemical methods allow simple polymer synthesis, removal of the template and readout. Among the different sensor configurations electrochemical MIP-sensors provide the broadest spectrum of protein analytes. The sensitivity of MIP-sensors is sufficiently high for biomarkers in the sub-nanomolar region, nevertheless the cross-reactivity of highly abundant proteins in human serum is still a challenge. MIPs for proteins offer innovative tools not only for clinical and environmental analysis, but also for bioimaging, therapy and protein engineering.
Heat waves are increasingly common in many countries across the globe, and also in Germany, where this study is set. Heat poses severe health risks, especially for vulnerable groups such as the elderly and children. This case study explores visitors' behavior and perceptions during six weekends in the summer of 2018 at a 6-month open-air horticultural show. Data from a face-to-face survey (n = 306) and behavioral observations ( n = 2750) were examined by using correlation analyses, ANOVA, and multiple regression analyses. Differences in weather perception, risk awareness, adaptive behavior, and activity level were observed between rainy days (maximum daily temperature, 25 degrees C), warmsummer days (25 degrees-30 degrees C), and hot days (>30 degrees C). Respondents reported a high level of heat risk awareness, butmost (90%) were unaware of actual heat warnings. During hot days, more adaptive measures were reported and observed. Older respondents reported taking the highest number of adaptive measures. We observed the highest level of adaptation in children, but they also showed the highest activity level. From our results we discuss how to facilitate individual adaptation to heat stress at open-air events by taking the heterogeneity of visitors into account. To mitigate negative health outcomes for citizens in the future, we argue for tailored risk communication aimed at vulnerable groups. <br /> SIGNIFICANCE STATEMENT: People around the world are facing higher average temperatures. While higher temperatures make open-air events a popular leisure time activity in summer, heat waves are a threat to health and life. Since there is not much research on how visitors of such events perceive different weather conditions-especially hot temperatures-we explored this in our case study in southern Germany at an open-air horticultural show in the summer of 2018. We discovered deficits both in people's awareness of current heat risk and the heat adaptation they carry out themselves. Future research should further investigate risk perception and adaptation behavior of private individuals, whereas event organizers and authorities need to continually focus on risk communication and facilitate individual adaptation of their visitors.
Recent experiments and theory suggest that ground state properties and reactivity of molecules can be modified when placed inside a nanoscale cavity, giving rise to strong coupling between vibrational modes and the quantized cavity field. This is commonly thought to be caused either by a cavity-distorted Born-Oppenheimer ground state potential or by the formation of light-matter hybrid states, vibrational polaritons. Here, we systematically study the effect of a cavity on ground state properties and infrared spectra of single molecules, considering vibration-cavity coupling strengths from zero up to the vibrational ultrastrong coupling regime. Using single-mode models for Li-H and O-H stretch modes and for the NH3 inversion mode, respectively, a single cavity mode in resonance with vibrational transitions is coupled to position-dependent molecular dipole functions. We address the influence of the cavity mode on polariton ground state energies, equilibrium bond lengths, dissociation energies, activation energies for isomerization, and on vibro-polaritonic infrared spectra. In agreement with earlier work, we observe all mentioned properties being strongly affected by the cavity, but only if the dipole self-energy contribution in the interaction Hamiltonian is neglected. When this term is included, these properties do not depend significantly on the coupling anymore. Vibro-polaritonic infrared spectra, in contrast, are always affected by the cavity mode due to the formation of excited vibrational polaritons. It is argued that the quantized nature of vibrational polaritons is key to not only interpreting molecular spectra in cavities but also understanding the experimentally observed modification of molecular reactivity in cavities.
We introduce a thermofield-based formulation of the multilayer multiconfigurational time-dependent Hartree (MCTDH) method to study finite temperature effects on non-adiabatic quantum dynamics from a non-stochastic, wave function perspective. Our approach is based on the formal equivalence of bosonic many-body theory at zero temperature with a doubled number of degrees of freedom and the thermal quasi-particle representation of bosonic thermofield dynamics (TFD). This equivalence allows for a transfer of bosonic many-body MCTDH as introduced by Wang and Thoss to the finite temperature framework of thermal quasi-particle TFD. As an application, we study temperature effects on the ultrafast internal conversion dynamics in pyrazine. We show that finite temperature effects can be efficiently accounted for in the construction of multilayer expansions of thermofield states in the framework presented herein. Furthermore, we find our results to agree well with existing studies on the pyrazine model based on the pMCTDH method.
The stable operation of a turbulent combustor is not completely silent; instead, there is a background of small amplitude aperiodic acoustic fluctuations known as combustion noise. Pressure fluctuations during this state of combustion noise are multifractal due to the presence of multiple temporal scales that contribute to its dynamics. However, existing models are unable to capture the multifractality in the pressure fluctuations. We conjecture an underlying fractional dynamics for the thermoacoustic system and obtain a fractional-order model for pressure fluctuations. The data from this model has remarkable visual similarity to the experimental data and also has a wide multifractal spectrum during the state of combustion noise. Quantitative similarity with the experimental data in terms of the Hurst exponent and the multifractal spectrum is observed during the state of combustion noise. This model is also able to produce pressure fluctuations that are qualitatively similar to the experimental data acquired during intermittency and thermoacoustic instability. Furthermore, we argue that the fractional dynamics vanish as we approach the state of thermoacoustic instability.
Sustainable electricity systems need renewable and dispatchable energy sources. Solar energy is an abundant source of renewable energy globally which is, though, by nature only available during the day, and especially in clear weather conditions. We compare three technology configurations able to provide dispatchable solar power at times without sunshine: Photovoltaics (PV) combined with battery (BESS) or thermal energy storage (TES) and concentrating solar power (CSP) with TES. Modeling different periods without sunshine, we find that PV+BESS is competitive for shorter storage durations while CSP+TES gains economic advantages for longer storage periods (also over PV+TES). The corresponding tipping points lie at 2-3 hours (current cost), and 4-10 hours if expectations on future cost developments are taken into consideration. PV+TES becomes only more competitive than CSP+TES with immense additional cost reductions of PV. Hence, there remain distinct niches for two technologies: PV+BESS for short storage durations and CSP+TES for longer ones.
Non-linear dielectric spectroscopy (NLDS) is employed as an effective tool to study relaxation processes and phase transitions of a poly(vinylidenefluoride-trifluoroethylene-chlorofluoroethylene) (P(VDF-TrFE-CFE)) relaxor-ferroelectric (R-F) terpolymer in detail. Measurements of the non-linear dielectric permittivity epsilon 2 ' reveal peaks at 30 and 80 degrees C that cannot be identified in conventional dielectric spectroscopy. By combining the results from NLDS experiments with those from other techniques such as thermally stimulated depolarization and dielectric-hysteresis studies, it is possible to explain the processes behind the additional peaks. The former peak, which is associated with the mid-temperature transition, is found in all other vinylidene fluoride-based polymers and may help to understand the non-zero epsilon 2 ' values that are detected on the paraelectric phase of the terpolymer. The latter peak can also be observed during cooling of P(VDF-TrFE) copolymer samples at 100 degrees C and is due to conduction and space-charge polarization as a result of the accumulation of real charges at the electrode-sample interface.
Spherulite-related space-charge electret properties of polypropylene (PP) have been widely discussed in the past decades. In the present paper, a less-common crystalline structure in PP-transcrystalline PP-is studied regarding its electret behavior in comparison with the typical spherulitic morphology. Polarized light microscopy and differential scanning calorimetry were employed to characterize the crystallite types and crystallinities of transcrystalline and spherulitic PP. Their electret functionality is investigated by means of thermally stimulated discharge experiments, where the cross-over phenomenon is observed on transcrystalline PP films, whereas surface-potential saturation and undercharging on the surface occur on the spherulitic samples. Besides, an asymmetrical behavior of positive and negative surface-charge stabilities is found on PP with spherulites, the negatively charged spherulitic surfaces show a better charge stability. It is shown that PP electrets are very sensitive to changes in the microscopic crystalline structures and their interfaces as well as in the molecular conformations controlled through adjustments of the respective processing steps. In addition, surface and bulk nanocomposites of PP or low-density polyethylene with inorganic particles are included in the comparison. In view of recent developments in the areas of PP-based electret-fiber filters and cellular-foam ferroelectrets, the observed changes in the charge-storage properties may have particular relevance, as the required film, fiber, or foam processing might significantly modify crystalline morphologies and nano-scale interfaces in PP electrets. Limitations in the charge-storage capabilities of interface structures may also be of interest in the context of high-voltage electrical-insulation materials where reduced space-charge accumulation and slightly increased charge transport can be advantageous.
Destabilization of super-rotating Taylor-Couette flows by current-free helical magnetic fields
(2021)
In an earlier paper we showed that the combination of azimuthal magnetic fields and super-rotation in Taylor-Couette flows of conducting fluids can be unstable against non-axisymmetric perturbations if the magnetic Prandtl number of the fluid is Pm not equal 1. Here we demonstrate that the addition of a weak axial field component allows axisymmetric perturbation patterns for Pm of order unity depending on the boundary conditions. The axisymmetric modes only occur for magnetic Mach numbers (of the azimuthal field) of order unity, while higher values are necessary for the non-axisymmetric modes. The typical growth time of the instability and the characteristic time scale of the axial migration of the axisymmetric mode are long compared with the rotation period, but short compared with the magnetic diffusion time. The modes travel in the positive or negative z direction along the rotation axis depending on the sign of B phi Bz. We also demonstrate that the azimuthal components of flow and field perturbations travel in phase if vertical bar B phi vertical bar >> vertical bar B-z vertical bar, independent of the form of the rotation law. Within a short-wave approximation for thin gaps it is also shown (in an appendix) that for ideal fluids the considered helical magnetorotational instability only exists for rotation laws with negative shear.
Symmetry and its breaking crucially define the chemical properties of molecules and their functionality. Resonant inelastic X-ray scattering is a local electronic structure probe reporting on molecular symmetry and its dynamical breaking within the femtosecond scattering duration. Here, we study pyrimidine, a system from the C-2v point group, in an aqueous solution environment, using scattering though its 2a(2) resonance. Despite the absence of clean parity selection rules for decay transitions from in-plane orbitals, scattering channels including decay from the 7b(2) and 11a(1) orbitals with nitrogen lone pair character are a direct probe for molecular symmetry. Computed spectra of explicitly solvated molecules sampled from a molecular dynamics simulation are combined with the results of a quantum dynamical description of the X-ray scattering process. We observe dominant signatures of core-excited Jahn-Teller induced symmetry breaking for resonant excitation. Solvent contributions are separable by shortening of the effective scattering duration through excitation energy detuning.
The DNA interaction with cis-isomers of photosensitive azobenzene-containing surfactants was studied by both experimental methods and computer simulation. It was shown that before the organization of micelles, such surfactants in the cis-conformation form associates of only a single type with a disordered orientation of molecules. In contrast, for trans-isomers, there exist two types of associates with head-to-head or head-to-tail orientations of molecules in dependence on salt concentration in a solution. The comparison of cis- and trans-isomer binding to DNA and the influence of salt concentration on the formation of their complexes with DNA were studied. It was shown that cis-isomers interact with phosphate groups of DNA and that their molecules were also located along the minor groove of DNA.
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.
A theory for diffusivity estimation for spatially extended activator-inhibitor dynamics modeling the evolution of intracellular signaling networks is developed in the mathematical framework of stochastic reaction-diffusion systems. In order to account for model uncertainties, we extend the results for parameter estimation for semilinear stochastic partial differential equations, as developed in Pasemann and Stannat (Electron J Stat 14(1):547-579, 2020), to the problem of joint estimation of diffusivity and parametrized reaction terms. Our theoretical findings are applied to the estimation of effective diffusivity of signaling components contributing to intracellular dynamics of the actin cytoskeleton in the model organism Dictyostelium discoideum.
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.
Fires are a fundamental part of the Earth System. In the last decades, they have been altering ecosystem structure, biogeochemical cycles and atmospheric composition with unprecedented rapidity. In this study, we implement a complex networks-based methodology to track individual fires over space and time. We focus on extreme fires-the 5% most intense fires-in the tropical forests of the Brazilian Legal Amazon over the period 2002-2019. We analyse the interannual variability in the number and spatial patterns of extreme forest fires in years with diverse climatic conditions and anthropogenic pressure to examine potential synergies between climate and anthropogenic drivers. We observe that major droughts, that increase forest flammability, co-occur with high extreme fire years but also that it is fundamental to consider anthropogenic activities to understand the distribution of extreme fires. Deforestation fires, fires escaping from managed lands, and other types of forest degradation and fragmentation provide the ignition sources for fires to ignite in the forests. We find that all extreme forest fires identified are located within a 0.5-km distance from forest edges, and up to 56% of them are within a 1-km distance from roads (which increases to 73% within 5 km), showing a strong correlation that defines spatial patterns of extreme fires.
We report on the multiple response of microgels triggered by a single optical stimulus. Under irradiation, the volume of the microgels is reversibly switched by more than 20 times. The irradiation initiates two different processes: photo-isomerization of the photo-sensitive surfactant, which forms a complex with the anionic microgel, rendering it photo-responsive; and local heating due to a thermo-plasmonic effect within the structured gold layer on which the microgel is deposited. The photo-responsivity is related to the reversible accommodation/release of the photo-sensitive surfactant depending on its photo-isomerization state, while the thermo-sensitivity is intrinsically built in. We show that under exposure to green light, the thermo-plasmonic effect generates a local hot spot in the gold layer, resulting in the shrinkage of the microgel. This process competes with the simultaneous photo-induced swelling. Depending on the position of the laser spot, the spatiotemporal control of reversible particle shrinking/swelling with a predefined extent on a per-second base can be implemented.
We introduce and study a family of lattice equations which may be viewed either as a strongly nonlinear discrete extension of the Gardner equation, or a non-convex variant of the Lotka-Volterra chain. Their deceptively simple form supports a very rich family of complex solitary patterns. Some of these patterns are also found in the quasi-continuum rendition, but the more intriguing ones, like interlaced pairs of solitary waves, or waves which may reverse their direction either spontaneously or due a collision, are an intrinsic feature of the discrete realm.
We investigate a diffusion process with a time-dependent diffusion coefficient, both exponentially increasing and decreasing in time, D(t)=D-0(e +/- 2 alpha t). For this (hypothetical) nonstationary diffusion process we compute-both analytically and from extensive stochastic simulations-the behavior of the ensemble- and time-averaged mean-squared displacements (MSDs) of the particles, both in the over- and underdamped limits. Simple asymptotic relations derived for the short- and long-time behaviors are shown to be in excellent agreement with the results of simulations. The diffusive characteristics in the presence of ageing are also considered, with dramatic differences of the over- versus underdamped regime. Our results for D(t)=D-0(e +/- 2 alpha t) extend and generalize the class of diffusive systems obeying scaled Brownian motion featuring a power-law-like variation of the diffusivity with time, D(t) similar to t(alpha-1). We also examine the logarithmically increasing diffusivity, D(t)=D(0)log[t/tau(0)], as another fundamental functional dependence (in addition to the power-law and exponential) and as an example of diffusivity slowly varying in time. One of the main conclusions is that the behavior of the massive particles is predominantly ergodic, while weak ergodicity breaking is repeatedly found for the time-dependent diffusion of the massless particles at short times. The latter manifests itself in the nonequivalence of the (both nonaged and aged) MSD and the mean time-averaged MSD. The current findings are potentially applicable to a class of physical systems out of thermal equilibrium where a rapid increase or decrease of the particles' diffusivity is inherently realized. One biological system potentially featuring all three types of time-dependent diffusion (power-law-like, exponential, and logarithmic) is water diffusion in the brain tissues, as we thoroughly discuss in the end.
Writing a history of a scientific theory is always difficult because it requires to focus on some key contributors and to "reconstruct" some supposed influences. In the 1970s, a new way of performing science under the name "chaos" emerged, combining the mathematics from the nonlinear dynamical systems theory and numerical simulations. To provide a direct testimony of how contributors can be influenced by other scientists or works, we here collected some writings about the early times of a few contributors to chaos theory. The purpose is to exhibit the diversity in the paths and to bring some elements-which were never published-illustrating the atmosphere of this period. Some peculiarities of chaos theory are also discussed.
Respect the surroundings
(2021)
Fourteen-month-olds' ability to distinguish a just learned word, /bu?k/, from its minimally different word, /du?k/, was assessed under two pre-exposure conditions: one where /b, d/-initial forms occurred in a varying vowel context and another where the vowel was fixed but the final consonant varied. Infants in the experiments benefited from the variable vowel but not from the variable final consonant context, suggesting that vowel variability but not all kinds of variability are beneficial. These results are discussed in the context of time-honored observations on the vowel-dependent nature of place of articulation cues for consonants.
Quenching mechanism of uranyl(VI) by chloride and bromide in aqueous and non-aqueous solutions
(2021)
A major hindrance in utilizing uranyl(VI) luminescence as a standard analytical tool, for example, in environmental monitoring or nuclear industries, is quenching by other ions such as halide ions, which are present in many relevant matrices of uranyl(VI) speciation. Here, we demonstrate through a combination of time-resolved laser-induced fluorescence spectroscopy, transient absorption spectroscopy, and quantum chemistry that coordinating solvent molecules play a crucial role in U(VI) halide luminescence quenching. We show that our previously suggested quenching mechanism based on an internal redox reaction of the 1:2-uranyl-halide-complex holds also true for bromide-induced quenching of uranyl(VI). By adopting specific organic solvents, we were able to suppress the separation of the oxidized halide ligand X-2(center dot-) and the formed uranyl(V) into fully solvated ions, thereby "reigniting" U(VI) luminescence. Time-dependent density functional theory calculations show that quenching occurs through the outer-sphere complex of U(VI) and halide in water, while the ligand-to-metal charge transfer is strongly reduced in acetonitrile.
How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x) = D0|x|gamma and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicitybreaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB similar to(1/r )-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics.
The first detections of black hole-neutron star mergers (GW200105 and GW200115) by the LIGO-Virgo-Kagra Collaboration mark a significant scientific breakthrough. The physical interpretation of pre- and postmerger signals requires careful cross-examination between observational and theoretical modelling results. Here we present the first set of black hole-neutron star simulations that were obtained with the numerical-relativity code BAM. Our initial data are constructed using the public LORENE spectral library, which employs an excision of the black hole interior. BAM, in contrast, uses the moving-puncture gauge for the evolution. Therefore, we need to "stuff" the black hole interior with smooth initial data to evolve the binary system in time. This procedure introduces constraint violations such that the constraint damping properties of the evolution system are essential to increase the accuracy of the simulation and in particular to reduce spurious center-of-mass drifts. Within BAM we evolve the Z4c equations and we compare our gravitational-wave results with those of the SXS collaboration and results obtained with the SACRA code. While we find generally good agreement with the reference solutions and phase differences less than or similar to 0.5 rad at the moment of merger, the absence of a clean convergence order in our simulations does not allow for a proper error quantification. We finally present a set of different initial conditions to explore how the merger of black hole neutron star systems depends on the involved masses, spins, and equations of state.
We investigate the effects of Markovian resetting events on continuous time random walks where the waiting times and the jump lengths are random variables distributed according to power-law probability density functions.
We prove the existence of a nonequilibrium stationary state and finite mean first arrival time.
However, the existence of an optimum reset rate is conditioned to a specific relationship between the exponents of both power-law tails.
We also investigate the search efficiency by finding the optimal random walk which minimizes the mean first arrival time in terms of the reset rate, the distance of the initial position to the target, and the characteristic transport exponents.
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.
State-of-the-art organic solar cells exhibit power conversion efficiencies of 18% and above. These devices benefit from the suppression of free charge recombination with regard to the Langevin limit of charge encounter in a homogeneous medium. It is recognized that the main cause of suppressed free charge recombination is the reformation and resplitting of charge-transfer (CT) states at the interface between donor and acceptor domains. Here, we use kinetic Monte Carlo simulations to understand the interplay between free charge motion and recombination in an energetically disordered phase-separated donor-acceptor blend. We identify conditions for encounter-dominated and resplitting-dominated recombination. In the former regime, recombination is proportional to mobility for all parameters tested and only slightly reduced with respect to the Langevin limit. In contrast, mobility is not the decisive parameter that determines the nongeminate recombination coefficient, k(2), in the latter case, where k2 is a sole function of the morphology, CT and charge-separated (CS) energetics, and CT-state decay properties. Our simulations also show that free charge encounter in the phase-separated disordered blend is determined by the average mobility of all carriers, while CT reformation and resplitting involves mostly states near the transport energy. Therefore, charge encounter is more affected by increased disorder than the resplitting of the CT state. As a consequence, for a given mobility, larger energetic disorder, in combination with a higher hopping rate, is preferred. These findings have implications for the understanding of suppressed recombination in solar cells with nonfullerene acceptors, which are known to exhibit lower energetic disorder than that of fullerenes.
We consider the spatiotemporal states of an ensemble of nonlocally coupled nonidentical phase oscillators, which correspond to different regimes of the long-term evolution of such a system. We have obtained homogeneous, twisted, and nonhomogeneous stationary solutions to the Ott-Antonsen equations corresponding to key variants of the realized collective rotational motion of elements of the medium in question with nonzero mesoscopic characteristics determining the degree of coherence of the dynamics of neighboring particles. We have described the procedures of the search for the class of nonhomogeneous solutions as stationary points of the auxiliary point map and of determining the stability based on analysis of the eigenvalue spectrum of the composite operator. Static and breather cluster regimes have been demonstrated and described, as well as the regimes with an irregular behavior of averaged complex fields including, in particular, the local order parameter.
We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed.
This paper is concerned with correlation functions of stochastic systems with memory, a prominent example being a molecule or colloid moving through a complex (e.g. viscoelastic) fluid environment. Analytical investigations of such systems based on non-Markovian stochastic equations are notoriously difficult. A common approximation is that of a single-exponential memory, corresponding to the introduction of one auxiliary variable coupled to the Markovian dynamics of the main variable. As a generalization, we here investigate a class of 'toy' models with altogether three degrees of freedom, giving rise to more complex forms of memory. Specifically, we consider, mainly on an analytical basis, the under- and overdamped motion of a colloidal particle coupled linearly to two auxiliary variables, where the coupling between variables can be either reciprocal or non-reciprocal. Projecting out the auxiliary variables, we obtain non-Markovian Langevin equations with friction kernels and colored noise, whose structure is similar to that of a generalized Langevin equation. For the present systems, however, the non-Markovian equations may violate the fluctuation-dissipation relation as well as detailed balance, indicating that the systems are out of equilibrium. We then study systematically the connection between the coupling topology of the underlying Markovian system and various autocorrelation functions. We demonstrate that already two auxiliary variables can generate surprisingly complex (e.g. non-monotonic or oscillatory) memory and correlation functions. Finally, we show that a minimal overdamped model with two auxiliary variables and suitable non-reciprocal coupling yields correlation functions resembling those describing hydrodynamic backflow in an optical trap.
How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM.
Development of functional and stable solid polymer electrolytes (SPEs) for battery applications is an important step towards both safer batteries and for the realization of lithium-based or anode-less batteries. The interface between the lithium and the solid polymer electrolyte is one of the bottlenecks, where severe degradation is expected. Here, the stability of three different SPEs - poly(ethylene oxide) (PEO), poly(epsilon-caprolactone) (PCL) and poly(trimethylene carbonate) (PTMC) - together with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt, is investigated after they have been exposed to lithium metal under UHV conditions. Degradation compounds, e.g. Li-O-R, LiF and LixSyOz, are identified for all SPEs using soft X-ray photoelectron spectroscopy. A competing degradation between polymer and salt is identified in the outermost surface region (<7 nm), and is dependent on the polymer host. PTMC:LiTFSI shows the most severe decomposition of both polymer and salt followed by PCL:LiTFSI and PEO:LiTFSI. In addition, the movement of lithium species through the decomposed interface shows large variation depending on the polymer electrolyte system.
We consider an array of nonlocally coupled oscillators on a ring, which for equally spaced units possesses a Kuramoto-Battogtokh chimera regime and a synchronous state. We demonstrate that disorder in oscillators positions leads to a transition from the synchronous to the chimera state. For a static (quenched) disorder we find that the probability of synchrony survival depends on the number of particles, from nearly zero at small populations to one in the thermodynamic limit. Furthermore, we demonstrate how the synchrony gets destroyed for randomly (ballistically or diffusively) moving oscillators. We show that, depending on the number of oscillators, there are different scalings of the transition time with this number and the velocity of the units.
Binary III-V nitrides such as AlN, GaN and InN in the wurtzite-type structure have long been considered as potent semiconducting materials because of their optoelectronic properties, amongst others. With rising concerns over the utilization of scarce elements, a replacement of the trivalent cations by others in ternary and multinary nitrides has led to the development of different variants of nitrides and oxide nitrides crystallizing in lower-symmetry variants of wurtzite. This work presents the symmetry relationships between these structural types specific to nitrides and oxide nitrides and updates some prior work on this matter. The non-existence of compounds crystallizing in Pmc2(1), formally the highest subgroup of the wurtzite type fulfilling Pauling's rules for 1:1:2 stoichiometries, has been puzzling scientists for a while; a rationalization is given, from a crystallographic basis, of why this space group is unlikely to be adopted.
Using magneto-optical Faraday and Kerr measurements, we investigate the magnetic and magnetooptical properties of a thick Bi-substituted gadolinium iron garnet film over a broad range of wavelengths (250-850 nm) and temperatures (150-300 K), including the magnetization compensation point, TM. We observe an exchange-bias-like effect in the vicinity of TM. By slightly changing the sample temperature, we can precisely tune the bias field, which reaches a magnitude 6 times higher than the coercive field. We explain this phenomenon by considering the short-range superexchange interaction and a change in the magnetic behavior when moving from the surface to the bulk of the film. This finding may lead to the development of single-film magneto-optical devices based on the exchange-bias effect.
Most image restoration methods in astronomy rely upon probabilistic tools that infer the best solution for a deconvolution problem. They achieve good performances when the point spread function (PSF) is spatially invariant in the image plane.
However, this condition is not always satisfied in real optical systems. We propose a new method for the restoration of images affected by static and anisotropic aberrations using Deep Neural Networks that can be directly applied to sky images.
The network is trained using simulated sky images corresponding to the T80-S Telescope optical model, a 80-cm survey imager at Cerro Tololo (Chile), which are synthesized using a Zernike polynomial representation of the optical system.
Once trained, the network can be used directly on sky images, outputting a corrected version of the image that has a constant and known PSF across its field of view. The method is to be tested on the T80-S Telescope.
We present the method and results on synthetic data.
Both ground- and satellite-based airglow imaging have significantly contributed to understanding the low-latitude ionosphere, especially the morphology and dynamics of the equatorial ionization anomaly (EIA). The NASA Global-scale Observations of the Limb and Disk (GOLD) mission focuses on far-ultraviolet airglow images from a geostationary orbit at 47.5 degrees W. This region is of particular interest at low magnetic latitudes because of the high magnetic declination (i.e., about -20 degrees) and proximity of the South Atlantic magnetic anomaly. In this study, we characterize an exciting feature of the nighttime EIA using GOLD observations from October 5, 2018 to June 30, 2020. It consists of a wavelike structure of a few thousand kilometers seen as poleward and equatorward displacements of the EIA-crests. Initial analyses show that the synoptic-scale structure is symmetric about the dip equator and appears nearly stationary with time over the night. In quasi-dipole coordinates, maxima poleward displacements of the EIA-crests are seen at about +/- 12 degrees latitude and around 20 and 60 degrees longitude (i.e., in geographic longitude at the dip equator, about 53 degrees W and 14 degrees W). The wavelike structure presents typical zonal wavelengths of about 6.7 x 10(3) km and 3.3 x 10(3) km. The structure's occurrence and wavelength are highly variable on a day-to-day basis with no apparent dependence on geomagnetic activity. In addition, a cluster or quasi-periodic wave train of equatorial plasma depletions (EPDs) is often detected within the synoptic-scale structure. We further outline the difference in observing these EPDs from FUV images and in situ measurements during a GOLD and Swarm mission conjunction.
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.
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.
In contrast to the common conception that the interfacial energy-level alignment is affixed once the interface is formed, we demonstrate that heterojunctions between organic semiconductors and metal-halide perovskites exhibit huge energy-level realignment during photoexcitation. Importantly, the photoinduced level shifts occur in the organic component, including the first molecular layer in direct contact with the perovskite. This is caused by charge-carrier accumulation within the organic semiconductor under illumination and the weak electronic coupling between the junction components.
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.
The spatial magnetic properties, through-space NMR shieldings (TSNMRSs), of stable O, S and Hal analogues of N-heterocyclic carbenes (NHCs) have been calculated using the GIAO perturbation method employing the nucleus-independent chemical shift (NICS) concept and the results visualized as iso-chemical-shielding surfaces (ICSSs) of various sizes and directions. The TSNMRS values (actually the anisotropy effects measurable in H-1 NMR spectroscopy) are employed to qualify and quantify the position of the present mesomeric equilibria (carbenes <-> ylides). The results are confirmed by geometry (bond angles and bond lengths), IR spectra, UV spectra, and C-13 chemical shifts of the electron-deficient carbon centers.
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.
A comet is a highly dynamic object, undergoing a permanent state of change. These changes have to be carefully classified and considered according to their intrinsic temporal and spatial scales. The Rosetta mission has, through its contiguous in-situ and remote sensing coverage of comet 67P/Churyumov-Gerasimenko (hereafter 67P) over the time span of August 2014 to September 2016, monitored the emergence, culmination, and winding down of the gas and dust comae. This provided an unprecedented data set and has spurred a large effort to connect in-situ and remote sensing measurements to the surface. In this review, we address our current understanding of cometary activity and the challenges involved when linking comae data to the surface. We give the current state of research by describing what we know about the physical processes involved from the surface to a few tens of kilometres above it with respect to the gas and dust emission from cometary nuclei. Further, we describe how complex multidimensional cometary gas and dust models have developed from the Halley encounter of 1986 to today. This includes the study of inhomogeneous outgassing and determination of the gas and dust production rates. Additionally, the different approaches used and results obtained to link coma data to the surface will be discussed. We discuss forward and inversion models and we describe the limitations of the respective approaches. The current literature suggests that there does not seem to be a single uniform process behind cometary activity. Rather, activity seems to be the consequence of a variety of erosion processes, including the sublimation of both water ice and more volatile material, but possibly also more exotic processes such as fracture and cliff erosion under thermal and mechanical stress, sub-surface heat storage, and a complex interplay of these processes. Seasons and the nucleus shape are key factors for the distribution and temporal evolution of activity and imply that the heliocentric evolution of activity can be highly individual for every comet, and generalisations can be misleading.
9,10-substituted anthracenes are known for their useful optical properties like fluorescence, which makes them frequently used probes in sensing applications. In this article, we investigate the fundamental photophysical properties of three pyridyl-substituted variants. The nitrogen atoms in the pyridinium six-membered rings are located in the ortho-, meta-, and para-positions in relation to the anthracene core. Absorption, fluorescence, and transient absorption measurements were carried out and were complemented by theoretical calculations. We monitored the photophysics of the anthracene derivatives in chloroform and water investigating the protonated as well as their nonprotonated forms. We found that the optical properties of the nonprotonated forms are strongly determined by the anthracene chromophore, with only small differences to other 9,10-substituted anthracenes, for example diphenyl anthracene. In contrast, protonation leads to a strong decrease in fluorescence intensity and lifetime. Transient absorption measurements and theoretical calculations revealed the formation of a charge-transfer state in the protonated chromophores, where electron density is shifted from the anthracene moiety toward the protonated pyridyl substituents. While the para- and ortho-derivatives' charge transfer is still moderately fluorescent, the meta-derivative is affected much stronger and shows nearly no fluorescence. This nitrogen-atom-position-dependent sensitivity to hydronium activity makes a combination of these fluorophores very attractive for pH-sensing applications covering a broadened pH range.
We consider a perturbation of the de Rham complex on a compact manifold with boundary. This perturbation goes beyond the framework of complexes, and so cohomology does not apply to it. On the other hand, its curvature is "small", hence there is a natural way to introduce an Euler characteristic and develop a Lefschetz theory for the perturbation. This work is intended as an attempt to develop a cohomology theory for arbitrary sequences of linear mappings.
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.
The harmonic oscillator is a powerful model that can appear as a limit case when examining a nonlinear system. A well known fact is that, without driving, the inclusion of a friction term makes the origin of the phase space-which is a fixed point of the system-linearly stable. In this work, we include a telegraph process as perturbation of the oscillator's frequency, for example, to describe the motion of a particle with fluctuating charge gyrating in an external magnetic field. Increasing intensity of this colored noise is capable of changing the quality of the fixed point. To characterize the stability of the system, we use a stability measure that describes the growth of the displacement of the system's phase space position and express it in a closed form. We expand the respective exponent for light friction and low noise intensity and compare both the exact analytic solution and the expansion to numerical values. Our findings allow stability predictions for several physical systems.
We present an efficient technique for control of synchrony in a globally coupled ensemble by pulsatile action. We assume that we can observe the collective oscillation and can stimulate all elements of the ensemble simultaneously. We pay special attention to the minimization of intervention into the system. The key idea is to stimulate only at the most sensitive phase. To find this phase, we implement an adaptive feedback control. Estimating the instantaneous phase of the collective mode on the fly, we achieve efficient suppression using a few pulses per oscillatory cycle. We discuss the possible relevance of the results for neuroscience, namely, for the development of advanced algorithms for deep brain stimulation, a medical technique used to treat Parkinson's disease.
Levy walks (LWs) are spatiotemporally coupled random-walk processes describing superdiffusive heat conduction in solids, propagation of light in disordered optical materials, motion of molecular motors in living cells, or motion of animals, humans, robots, and viruses. We here investigate a key feature of LWs-their response to an external harmonic potential. In this generic setting for confined motion we demonstrate that LWs equilibrate exponentially and may assume a bimodal stationary distribution. We also show that the stationary distribution has a horizontal slope next to a reflecting boundary placed at the origin, in contrast to correlated superdiffusive processes. Our results generalize LWs to confining forces and settle some longstanding puzzles around LWs.
Machine learning control
(2020)
Recently, the term explainable AI came into discussion as an approach to produce models from artificial intelligence which allow interpretation. For a long time, symbolic regression has been used to produce explainable and mathematically tractable models. In this contribution, we extend previous work on symbolic regression methods to infer the optimal control of a dynamical system given one or several optimization criteria, or cost functions. In earlier publications, network control was achieved by automated machine learning control using genetic programming. Here, we focus on the subsequent path continuation analysis of the mathematical expressions which result from the machine learning model. In particular, we use AUTO to analyze the solution properties of the controlled oscillator system which served as our model. As a result, we show that there is a considerable advantage of explainable symbolic regression models over less accessible neural networks. In particular, the roadmap of future works may be to integrate such analyses into the optimization loop itself to filter out robust solutions by construction.
Dielectric elastomer devices operate on the principle of Maxwell stress and their operating performance significantly rely on the elastomer and compliant electrode's electrical and mechanical properties. This paper reports that performing actuation tests at elevated temperatures resulted in an enhanced performance due to the reduction of Young's modulus and the increase of dielectric permittivity. As a result, considerably higher isometric forces and isotonic strains were achieved above the ambient operating temperature. For actuators made of silicone, polyurethane and acrylic elastomers, 166%, 70% and 266% higher isometric forces and 450%, 250% and 54% higher isotonic strains were observed, respectively, when tested at the temperature of 100 degrees C in comparison to ambient temperature values using the same operating voltages. Values of up to 0.4 J kg(-1) and 3.1 W kg(-1) were achieved for the work and power outputs per mass, respectively, on a silicone elastomer driven with a voltage of 1.5 kV at a temperature of 100 degrees C.
Nanoporous carbon materials (NCMs) provide the "function" of high specific surface area and thus have large interface area for interactions with surrounding species, which is of particular importance in applications related to adsorption processes. The strength and mechanism of adsorption depend on the pore architecture of the NCMs. In addition, chemical functionalization can be used to induce changes of electron density and/or electron density distribution in the pore walls, thus further modifying the interactions between carbons and guest species. Typical approaches for functionalization of nanoporous materials with regular atomic construction like porous silica, metal-organic frameworks, or zeolites, cannot be applied to NCMs due to their less defined local atomic construction and abundant defects. Therefore, synthetic strategies that offer a higher degree of control over the process of functionalization are needed. Synthetic approaches for covalent functionalization of NCMs, that is, for the incorporation of heteroatoms into the carbon backbone, are critically reviewed with a special focus on strategies following the concept "from molecules to materials." Approaches for coordinative functionalization with metallic species, and the functionalization by nanocomposite formation between pristine carbon materials and heteroatom-containing carbons, are introduced as well. Particular focus is given to the influences of these functionalizations in adsorption-related applications.
This paper focuses on the ground state phase diagram of a 1D spin-1/2 quantum Ising model with competing first and second nearest neighbour interactions known as the axial next nearest neighbour Ising model in the presence of a transverse magnetic field. Here, using quantum correlations, both numerically and analytically, some evidence is provided to clarify the identification of the ground state phase diagram. Local quantum correlations play a crucial role in detecting the critical lines either revealed or hidden by symmetry-breaking. A non-symmetry-breaking disorder transition line can be identified by the first derivative of both entanglement of formation and quantum discord between nearest neighbour spins. In addition, the quantum correlations between the second neighbour spins can also be used to reveal Kosterlitz-Thouless phase transition when their interaction strength grows and becomes closer to the first nearest neighbour one. The results obtained using the Jordan-Wigner transformation confirm the accuracy of the numerical case.
IMDfence
(2020)
Over the past decade, focus on the security and privacy aspects of implantable medical devices (IMDs) has intensified, driven by the multitude of cybersecurity vulnerabilities found in various existing devices. However, due to their strict computational, energy and physical constraints, conventional security protocols are not directly applicable to IMDs. Custom-tailored schemes have been proposed instead which, however, fail to cover the full spectrum of security features that modern IMDs and their ecosystems so critically require. In this paper we propose IMDfence, a security protocol for IMD ecosystems that provides a comprehensive yet practical security portfolio, which includes availability, non-repudiation, access control, entity authentication, remote monitoring and system scalability. The protocol also allows emergency access that results in the graceful degradation of offered services without compromising security and patient safety. The performance of the security protocol as well as its feasibility and impact on modern IMDs are extensively analyzed and evaluated. We find that IMDfence achieves the above security requirements at a mere less than 7% increase in total IMD energy consumption, and less than 14 ms and 9 kB increase in system delay and memory footprint, respectively.
We consider large networks of theta neurons on a ring, synaptically coupled with an asymmetric kernel. Such networks support stable "bumps" of activity, which move along the ring if the coupling kernel is asymmetric. We investigate the effects of the kernel asymmetry on the existence, stability, and speed of these moving bumps using continuum equations formally describing infinite networks. Depending on the level of heterogeneity within the network, we find complex sequences of bifurcations as the amount of asymmetry is varied, in strong contrast to the behavior of a classical neural field model.
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.
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.
In this study we investigate, using all-atom molecular-dynamics computer simulations, the in-plane diffusion of a doxorubicin drug molecule in a thin film of water confined between two silica surfaces. We find that the molecule diffuses along the channel in the manner of a Gaussian diffusion process, but with parameters that vary according to its varying transversal position. Our analysis identifies that four Gaussians, each describing particle motion in a given transversal region, are needed to adequately describe the data. Each of these processes by itself evolves with time at a rate slower than that associated with classical Brownian motion due to a predominance of anticorrelated displacements. Long adsorption events lead to ageing, a property observed when the diffusion is intermittently hindered for periods of time with an average duration which is theoretically infinite. This study presents a simple system in which many interesting features of anomalous diffusion can be explored. It exposes the complexity of diffusion in nanoconfinement and highlights the need to develop new understanding.
Nonstationary coherence-incoherence patterns in nonlocally coupled heterogeneous phase oscillators
(2020)
We consider a large ring of nonlocally coupled phase oscillators and show that apart from stationary chimera states, this system also supports nonstationary coherence-incoherence patterns (CIPs). For identical oscillators, these CIPs behave as breathing chimera states and are found in a relatively small parameter region only. It turns out that the stability region of these states enlarges dramatically if a certain amount of spatially uniform heterogeneity (e.g., Lorentzian distribution of natural frequencies) is introduced in the system. In this case, nonstationary CIPs can be studied as stable quasiperiodic solutions of a corresponding mean-field equation, formally describing the infinite system limit. Carrying out direct numerical simulations of the mean-field equation, we find different types of nonstationary CIPs with pulsing and/or alternating chimera-like behavior. Moreover, we reveal a complex bifurcation scenario underlying the transformation of these CIPs into each other. These theoretical predictions are confirmed by numerical simulations of the original coupled oscillator system.
Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.
We present a diffusion-based simulation and theoretical models for explanation of the photoluminescence (PL) emission intensity in semiconductor nanoplatelets. It is shown that the shape of the PL intensity curves can be reproduced by the interplay of recombination, diffusion and trapping of excitons. The emission intensity at short times is purely exponential and is defined by recombination. At long times, it is governed by the release of excitons from surface traps and is characterized by a power-law tail. We show that the crossover from one limit to another is controlled by diffusion properties. This intermediate region exhibits a rich behaviour depending on the value of diffusivity. The proposed approach reproduces all the features of experimental curves measured for different nanoplatelet systems.
Context. The High Energy Stereoscopic System Galactic plane survey (HGPS) is to date the most comprehensive census of Galactic gamma -ray sources at very high energies (VHE; 100 GeV <= E <= 100 TeV). As a consequence of the limited sensitivity of this survey, the 78 detected gamma -ray sources comprise only a small and biased subsample of the overall population. The larger part consists of currently unresolved sources, which contribute to large-scale diffuse emission to a still uncertain amount.Aims. We study the VHE gamma -ray source population in the Milky Way. For this purpose population-synthesis models are derived based on the distributions of source positions, extents, and luminosities.Methods. Several azimuth-symmetric and spiral-arm models are compared for spatial source distribution. The luminosity and radius function of the population are derived from the source properties of the HGPS data set and are corrected for the sensitivity bias of the HGPS. Based on these models, VHE source populations are simulated and the subsets of sources detectable according to the HGPS are compared with HGPS sources.Results. The power-law indices of luminosity and radius functions are determined to range between -1.6 and -1.9 for luminosity and -1.1 and -1.6 for radius. A two-arm spiral structure with central bar is discarded as spatial distribution of VHE sources, while azimuth-symmetric distributions and a distribution following a four-arm spiral structure without bar describe the HGPS data reasonably well. The total number of Galactic VHE sources is predicted to be in the range from 800 to 7000 with a total luminosity and flux of (1.6-6.3) x 10(36) ph s(-1) and (3-15) x 10(-10) ph cm(-2) s(-1), respectively.Conclusions. Depending on the model, the HGPS sample accounts for (68-87)% of the emission of the population in the scanned region. This suggests that unresolved sources represent a critical component of the diffuse emission measurable in the HGPS. With the foreseen jump in sensitivity of the Cherenkov Telescope Array, the number of detectable sources is predicted to increase by a factor between 5 and 9.
The Colorado Learning Attitudes about Science Survey (CLASS) is an instrument which is widely used in physics education to characterize students' attitudes toward physics and learning physics and compare them with those of experts. While CLASS has been extensively validated for use in the context of higher education institutions in the United States, there has been less information about its use with European students. We have studied the structural, content, and substantive aspects of validity of CLASS by first doing a confirmatory factor analysis of N = 642 sets of student answers from the University of Helsinki, Finland. The students represented a culturally and demographically different subset of university physics students than in previous studies. The confirmatory factor analysis used a 3-factor, 15-item factor structure as a starting point and the resulting factor structure was similar to the original. Just minor modifications were needed for fit parameters to be in the acceptable range. We explored the differences by student interviews and consultation of experts. With the exception of one item, they supported the new 14-item, 3-factor structure. The results show that the interpretations made from CLASS results are mostly transferable, and CLASS remains a useful instrument for a wide variety of populations.
LetH be a Schrodinger operator defined on a noncompact Riemannianmanifold Omega, and let W is an element of L-infinity (Omega; R). Suppose that the operator H + W is critical in Omega, and let phi be the corresponding Agmon ground state. We prove that if u is a generalized eigenfunction ofH satisfying vertical bar u vertical bar <= C-phi in Omega for some constant C > 0, then the corresponding eigenvalue is in the spectrum of H. The conclusion also holds true if for some K is an element of Omega the operator H admits a positive solution in (Omega) over bar = Omega \ K, and vertical bar u vertical bar <= C psi in (Omega) over bar for some constant C > 0, where psi is a positive solution of minimal growth in a neighborhood of infinity in Omega. Under natural assumptions, this result holds also in the context of infinite graphs, and Dirichlet forms.
This study compares the classification of Azerbaijani fricatives based on two sets of features: (a) spectral moments, spectral peak, amplitude, duration, and (b) cepstral coefficients employing Hidden Markov Models to divide each fricative into three regions such that the variances of the measures within each region are minimized. The cepstral coefficients were found to be more reliable predictors in the classification of all nine Azerbaijani fricatives and the cepstral measures yielded highly successful classification rates (91.21% across both genders) in the identification of the full set of fricatives of Azerbaijani.
We present the use of modern machine learning approaches to suppress self-sustained collective oscillations typically signaled by ensembles of degenerative neurons in the brain. The proposed hybrid model relies on two major components: an environment of oscillators and a policy-based reinforcement learning block. We report a model-agnostic synchrony control based on proximal policy optimization and two artificial neural networks in an Actor-Critic configuration. A class of physically meaningful reward functions enabling the suppression of collective oscillatory mode is proposed. The synchrony suppression is demonstrated for two models of neuronal populations-for the ensembles of globally coupled limit-cycle Bonhoeffer-van der Pol oscillators and for the bursting Hindmarsh-Rose neurons using rectangular and charge-balanced stimuli.
Low-dimensional description for ensembles of identical phase oscillators subject to Cauchy noise
(2020)
We study ensembles of globally coupled or forced identical phase oscillators subject to independent white Cauchy noise. We demonstrate that if the oscillators are forced in several harmonics, stationary synchronous regimes can be exactly described with a finite number of complex order parameters. The corresponding distribution of phases is a product of wrapped Cauchy distributions. For sinusoidal forcing, the Ott-Antonsen low-dimensional reduction is recovered.
The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function.
Motivated by the observation of non-exponential run-time distributions of bacterial swimmers, we propose a minimal phenomenological model for taxis of active particles whose motion is controlled by an internal clock. The ticking of the clock depends on an external concentration field, e.g., a chemical substance. We demonstrate that these particles can detect concentration gradients and respond to them by moving up- or down-gradient depending on the clock design, albeit measurements of these fields are purely local in space and instantaneous in time. Altogether, our results open a new route in the study of directional navigation: we show that the use of a clock to control motility actions represents a generic and versatile toolbox to engineer behavioral responses to external cues, such as light, chemical, or temperature gradients.
We propose Mobius maps as a tool to model synchronization phenomena in coupled phase oscillators. Not only does the map provide fast computation of phase synchronization, it also reflects the underlying group structure of the sinusoidally coupled continuous phase dynamics. We study map versions of various known continuous-time collective dynamics, such as the synchronization transition in the Kuramoto-Sakaguchi model of nonidentical oscillators, chimeras in two coupled populations of identical phase oscillators, and Kuramoto-Battogtokh chimeras on a ring, and demonstrate similarities and differences between the iterated map models and their known continuous-time counterparts.
Alone but not lonely
(2020)
Context. Hot subdwarfs are core-helium burning stars that show lower masses and higher temperatures than canonical horizontal branch stars. They are believed to be formed when a red giant suffers an extreme mass-loss episode. Binary interaction is suggested to be the main formation channel, but the high fraction of apparently single hot subdwarfs (up to 30%) has prompted single star formation scenarios to be proposed.Aims. We investigate the possibility that hot subdwarfs could form without interaction by studying wide binary systems. If single formation scenarios were possible, there should be hot subdwarfs in wide binaries that have undergone no interaction.Methods. Angular momentum accretion during interaction is predicted to cause the hot subdwarf companion to spin up to the critical velocity. The effect of this should still be observable given the timescales of the hot subdwarf phase. To study the rotation rates of companions, we have analysed light curves from the Transiting Exoplanet Survey Satellite for all known hot subdwarfs showing composite spectral energy distributions indicating the presence of a main sequence wide binary companion. If formation without interaction were possible, that would also imply the existence of hot subdwarfs in very wide binaries that are not predicted to interact. To identify such systems, we have searched for common proper motion companions with projected orbital distances of up to 0.1 pc to all known spectroscopically confirmed hot subdwarfs using Gaia DR2 astrometry.Results. We find that the companions in composite hot subdwarfs show short rotation periods when compared to field main sequence stars. They display a triangular-shaped distribution with a peak around 2.5 days, similar to what is observed for young open clusters. We also report a shortage of hot subdwarfs with candidate common proper motion companions. We identify only 16 candidates after probing 2938 hot subdwarfs with good astrometry. Out of those, at least six seem to be hierarchical triple systems, in which the hot subdwarf is part of an inner binary.Conclusions. The observed distribution of rotation rates for the companions in known wide hot subdwarf binaries provides evidence of previous interaction causing spin-up. Additionally, there is a shortage of hot subdwarfs in common proper motion pairs, considering the frequency of such systems among progenitors. These results suggest that binary interaction is always required for the formation of hot subdwarfs.
A conducting Taylor-Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler-Spruit type. The flows are unstable against non-axisymmetric perturbations which form electromotive forces defining a effect and eddy diffusivity. If both degenerated modes with m = +/- 1 are excited with the same power then the global a effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite alpha effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an alpha(2) dynamo. Moreover, an alpha Omega dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of Rm similar or equal to 2000 to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, eta(T) similar or equal to eta. By basically lower or larger eddy diffusivities the dynamo instability is suppressed.
Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges.
Classical option pricing schemes assume that the value of a financial asset follows a geometric Brownian motion (GBM). However, a growing body of studies suggest that a simple GBM trajectory is not an adequate representation for asset dynamics, due to irregularities found when comparing its properties with empirical distributions. As a solution, we investigate a generalisation of GBM where the introduction of a memory kernel critically determines the behaviour of the stochastic process. We find the general expressions for the moments, log-moments, and the expectation of the periodic log returns, and then obtain the corresponding probability density functions using the subordination approach. Particularly, we consider subdiffusive GBM (sGBM), tempered sGBM, a mix of GBM and sGBM, and a mix of sGBMs. We utilise the resulting generalised GBM (gGBM) in order to examine the empirical performance of a selected group of kernels in the pricing of European call options. Our results indicate that the performance of a kernel ultimately depends on the maturity of the option and its moneyness.
As VLSI design moves into ultra-deep-submicron technologies, timing margins added to the clock period are mandatory, to ensure correct circuit behavior under worst-case conditions. Timing resilient architectures emerged as a promising solution to alleviate these worst-case timing margins. These architectures allow improving system performance and reducing energy consumption. Asynchronous systems, on the other hand, have the potential to improve energy efficiency and performance. Blade is an asynchronous timing resilient template that leverages the advantages of both asynchronous and timing resilient techniques. However, Blade still presents challenges regarding its testability, which hinders its commercial or large-scale application. This paper demonstrates that scan chains can be prohibitive for Blade due to their high silicon costs., which can reach more than 100%. Then, it proposes an alternative test approach that allows concurrent testing, stuck-at, and delay testing. The test approach is based on the reuse the Blade features to provide testability, with silicon area overheads between 4 and 7%.
In this paper, we study the problem of formal verification for Answer Set Programming (ASP), namely, obtaining aformal proofshowing that the answer sets of a given (non-ground) logic programPcorrectly correspond to the solutions to the problem encoded byP, regardless of the problem instance. To this aim, we use a formal specification language based on ASP modules, so that each module can be proved to capture some informal aspect of the problem in an isolated way. This specification language relies on a novel definition of (possibly nested, first order)program modulesthat may incorporate local hidden atoms at different levels. Then,verifyingthe logic programPamounts to prove some kind of equivalence betweenPand its modular specification.
The European potential for renewable electricity is sufficient to enable fully renewable supply on different scales, from self-sufficient, subnational regions to an interconnected continent. We not only show that a continental-scale system is the cheapest, but also that systems on the national scale and below are possible at cost penalties of 20% or less. Transmission is key to low cost, but it is not necessary to vastly expand the transmission system. When electricity is transmitted only to balance fluctuations, the transmission grid size is comparable to today's, albeit with expanded cross-border capacities. The largest differences across scales concern land use and thus social acceptance: in the continental system, generation capacity is concentrated on the European periphery, where the best resources are. Regional systems, in contrast, have more dispersed generation. The key trade-off is therefore not between geographic scale and cost, but between scale and the spatial distribution of required generation and transmission infrastructure.
The determination of the voltage-dependent density of free charge carriers via capacitance spectroscopy is considered an important step in the analysis of emerging photovoltaic technologies, such as organic and perovskite solar cells. In particular, an intimate knowledge of the density of free charge carriers is required for the determination of crucial parameters such as the effective mobility, charge carrier lifetime, nongeminate recombination coefficients, average extraction times, and competition factors. Hence, it is paramount to verify the validity of the commonly employed approaches to obtain the density of free charge carriers. The advantages, drawbacks, and limitations of the most common approaches are investigated in detail and strategies to mitigate misleading values are explored. To this end, two types of nonfullerene organic solar cells based on a PTB7-Th:ITIC-2F blend and a PM6:Y6 blend, respectively, are used as a case study to assess how subsequent analyses of the nongeminate recombination dynamics depend on the chosen approach to calculate the density of free charge carriers via capacitance spectroscopy.
Understanding the recombination dynamics of organic and perovskite solar cells has been a crucial prerequisite in the steadily increasing performance of these promising new types of photovoltaics. Surface recombination in particular has turned out to be one of the last remaining roadblocks, which specifically reduces the open circuit voltage. In this study, the relationship between the rate of surface recombination and the density of charge carriers is analyzed, revealing a cubic dependence between these two parameters. This hypothesis is then tested and verified with the recombination dynamics of an organic solar cell known to exhibit significant surface recombination and a high energy proton irradiated CH3NH3PbI3 pemvskite solar cell during white light illumination. Incidentally, these results can also explain recombination orders exceeding the commonly known threshold for bimolecular recombination that have been observed in some studies without the need for a charge carrier dependent bimolecular recombination coefficient.