530 Physik
Refine
Has Fulltext
- no (93)
Year of publication
- 2021 (93) (remove)
Document Type
- Article (93) (remove)
Language
- English (93)
Is part of the Bibliography
- yes (93)
Keywords
- diffusion (5)
- anomalous diffusion (4)
- synchronization (3)
- ARPES (1)
- Actin cytoskeleton dynamics (1)
- Adaptation (1)
- Additive manufacturing (1)
- Approximation algorithms (1)
- Atlantic meridional overturning circulation (1)
- Auger–Meitner (1)
Institute
- Institut für Physik und Astronomie (71)
- Institut für Chemie (10)
- Institut für Mathematik (7)
- Institut für Umweltwissenschaften und Geographie (2)
- Department Linguistik (1)
- Extern (1)
- Hasso-Plattner-Institut für Digital Engineering gGmbH (1)
- Institut für Biochemie und Biologie (1)
- Sozialwissenschaften (1)
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.
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.
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.
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.
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.
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.
Bifurcations in rotating spherical shell convection under the influence of differential rotation
(2021)
The bifurcations of thermal convection in a rotating spherical shell heated from the inner sphere and driven by the buoyancy of a central gravity field are studied numerically. This model of spherical Rayleigh-Benard convection describes large-scale convection in planets and in the outer zones of celestial bodies. In this work, the influence of an additionally imposed differential rotation of the inner sphere with respect to the outer one on the heat transfer and, more generally, on the whole bifurcation structure is investigated. In addition to numerical simulations, path-following techniques are applied in order to compute both stable and unstable solution branches. The dynamics and the heat transfer are essentially determined by a global bifurcation, which we have identified as a homoclinic bifurcation that consists of a collision of a stable modulated rotating with an unstable rotating wave.
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 provide an overview of the tools and techniques of resurgence theory used in the Borel-ecalle resummation method, which we then apply to the massless Wess-Zumino model. Starting from already known results on the anomalous dimension of the Wess-Zumino model, we solve its renormalisation group equation for the two-point function in a space of formal series. We show that this solution is 1-Gevrey and that its Borel transform is resurgent. The Schwinger-Dyson equation of the model is then used to prove an asymptotic exponential bound for the Borel transformed two-point function on a star-shaped domain of a suitable ramified complex plane. This proves that the two-point function of the Wess-Zumino model is Borel-ecalle summable.
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.
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.
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.
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.
Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise
(2021)
We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments.
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.
Comb-like geometric constraints leading to emergence of the time-fractional Schrödinger equation
(2021)
This paper presents an overview over several examples, where the comb-like geometric constraints lead to emergence of the time-fractional Schrodinger equation. Motion of a quantum object on a comb structure is modeled by a suitable modification of the kinetic energy operator, obtained by insertion of the Dirac delta function in the Laplacian. First, we consider motion of a free particle on two- and three-dimensional comb structures, and then we extend the study to the interacting cases. A general form of a nonlocal term, which describes the interactions of the particle with the medium, is included in the Hamiltonian, and later on, the cases of constant and Dirac delta potentials are analyzed. At the end, we discuss the case of non-integer dimensions, considering separately the case of fractal dimension between one and two, and the case of fractal dimension between two and three. All these examples show that even though we are starting with the standard time-dependent Schrodinger equation on a comb, the time-fractional equation for the Green's functions appears, due to these specific geometric constraints.
We present a supervised learning method to learn the propagator map of a dynamical system from partial and noisy observations. In our computationally cheap and easy-to-implement framework, a neural network consisting of random feature maps is trained sequentially by incoming observations within a data assimilation procedure. By employing Takens's embedding theorem, the network is trained on delay coordinates. We show that the combination of random feature maps and data assimilation, called RAFDA, outperforms standard random feature maps for which the dynamics is learned using batch data.
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.
'Complex systems are information processors' is a statement that is frequently made. Here we argue for the distinction between information processing-in the sense of encoding and transmitting a symbolic representation-and the formation of correlations (pattern formation/self-organisation). The study of both uses tools from information theory, but the purpose is very different in each case: explaining the mechanisms and understanding the purpose or function in the first case, versus data analysis and correlation extraction in the latter. We give examples of both and discuss some open questions. The distinction helps focus research efforts on the relevant questions in each case.
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.
In this paper, we report X-ray absorption and core-level electron spectra of the nucleobase derivative 2-thiouracil at the sulfur L1- and L2,3-edges. We used soft X-rays from the free-electron laser FLASH2 for the excitation of isolated molecules and dispersed the outgoing electrons with a magnetic bottle spectrometer. We identified photoelectrons from the 2p core orbital, accompanied by an electron correlation satellite, as well as resonant and non-resonant Coster–Kronig and Auger–Meitner emission at the L1- and L2,3-edges, respectively. We used the electron yield to construct X-ray absorption spectra at the two edges. The experimental data obtained are put in the context of the literature currently available on sulfur core-level and 2-thiouracil spectroscopy.
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.
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.
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.
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.
With the advent of big data and data lakes, data are often integrated from multiple sources. Such integrated data are often of poor quality, due to inconsistencies, errors, and so forth. One way to check the quality of data is to infer functional dependencies (fds). However, in many modern applications it might be necessary to extract properties and relationships that are not captured through fds, due to the necessity to admit exceptions, or to consider similarity rather than equality of data values. Relaxed fds (rfds) have been introduced to meet these needs, but their discovery from data adds further complexity to an already complex problem, also due to the necessity of specifying similarity and validity thresholds. We propose Domino, a new discovery algorithm for rfds that exploits the concept of dominance in order to derive similarity thresholds of attribute values while inferring rfds. An experimental evaluation on real datasets demonstrates the discovery performance and the effectiveness of the proposed algorithm.
We consider an array of nonlocally coupled oscillators on a ring, which for equally spaced units possesses a Kuramoto-Battogtokh chimera regime and a synchronous state. We demonstrate that disorder in oscillators positions leads to a transition from the synchronous to the chimera state. For a static (quenched) disorder we find that the probability of synchrony survival depends on the number of particles, from nearly zero at small populations to one in the thermodynamic limit. Furthermore, we demonstrate how the synchrony gets destroyed for randomly (ballistically or diffusively) moving oscillators. We show that, depending on the number of oscillators, there are different scalings of the transition time with this number and the velocity of the units.
We 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 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.
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.
A characterization of the essential spectrum of Schrodinger operators on infinite graphs is derived involving the concept of R-limits. This concept, which was introduced previously for operators on N and Z(d) as "right-limits," captures the behaviour of the operator at infinity. For graphs with sub-exponential growth rate, we show that each point in sigma(ss)(H) corresponds to a bounded generalized eigenfunction of a corresponding R-limit of H. If, additionally, the graph is of uniform sub-exponential growth, also the converse inclusion holds.
Human migration is often studied using gravity models. These models, however, have known limitations, including analytic inconsistencies and a dependence on empirical data to calibrate multiple parameters for the region of interest. Overcoming these limitations, the radiation model has been proposed as an alternative, universal approach to predicting different forms of human mobility, but has not been adopted for studying migration. Here we show, using data on within-country migration from the USA and Mexico, that the radiation model systematically underpredicts long-range moves, while the traditional gravity model performs well for large distances. The universal opportunity model, an extension of the radiation model, shows an improved fit of long-range moves compared to the original radiation model, but at the cost of introducing two additional parameters. We propose a more parsimonious extension of the radiation model that introduces a single parameter. We demonstrate that it fits the data over the full distance spectrum and also-unlike the universal opportunity model-preserves the analytical property of the original radiation model of being equivalent to a gravity model in the limit of a uniform population distribution.
High-resolution observations of polar crown and high-latitude filaments are scarce. We present a unique sample of such filaments observed in high-resolution H alpha narrow-band filtergrams and broad-band images, which were obtained with a new fast camera system at the Vacuum Tower Telescope (VTT), Tenerife, Spain. The Chromospheric Telescope (ChroTel) provided full-disk context observations in H alpha, CaiiK, and Hei 10830 angstrom. The Helioseismic and Magnetic Imager (HMI) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) provided line-of-sight magnetograms and ultraviolet (UV) 1700 angstrom filtergrams, respectively. We study filigree in the vicinity of polar crown and high-latitude filaments and relate their locations to magnetic concentrations at the filaments' footpoints. Bright points are a well studied phenomenon in the photosphere at low latitudes, but they were not yet studied in the quiet network close to the poles. We examine size, area, and eccentricity of bright points and find that their morphology is very similar to their counterparts at lower latitudes, but their sizes and areas are larger. Bright points at the footpoints of polar crown filaments are preferentially located at stronger magnetic flux concentrations, which are related to bright regions at the border of supergranules as observed in UV filtergrams. Examining the evolution of bright points on three consecutive days reveals that their amount increases while the filament decays, which indicates they impact the equilibrium of the cool plasma contained in filaments.
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.
Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions
(2021)
We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.
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.
Global quantum thermometry
(2021)
A paradigm shift in quantum thermometry is proposed. To date, thermometry has relied on local estimation, which is useful to reduce statistical fluctuations once the temperature is very well known. In order to estimate temperatures in cases where few measurement data or no substantial prior knowledge are available, we build instead a method for global quantum thermometry. Based on scaling arguments, a mean logarithmic error is shown here to be the correct figure of merit for thermometry. Its full minimization provides an operational and optimal rule to postprocess measurements into a temperature reading, and it establishes a global precision limit. We apply these results to the simulated outcomes of measurements on a spin gas, finding that the local approach can lead to biased temperature estimates in cases where the global estimator converges to the true temperature. The global framework thus enables a reliable approach to data analysis in thermometry experiments.
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.
The non-equilibrium state of two oscillators with a mutual interaction and coupled to separate heat baths is discussed. Bosonic baths are considered, and an exact spectral representation for the elements of the covariance matrix is provided analytically. A wide class of spectral densities for the relevant bath modes is allowed for. The validity of the fluctuation-dissipation relation is established for global equilibrium (both baths at the same temperature) in the stationary state. Spectral measures of entanglement are suggested by comparing to the equilibrium spectrum of zero-point fluctuations. No rotating-wave approximation is applied, and anomalous heat transport from cold to hot bath, as reported in earlier work, is demonstrated not to occur.
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.
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.
The stratosphere is one of the main potential sources for subseasonal to seasonal predictability in midlatitudes in winter. The ability of an atmospheric model to realistically simulate the stratospheric dynamics is essential in order to move forward in the field of seasonal predictions in midlatitudes. Earlier studies with the ICOsahedral Nonhydrostatic atmospheric model (ICON) point out that stratospheric westerlies in ICON are underestimated. This is the first extensive study on the evaluation of Northern Hemisphere stratospheric winter circulation with ICON in numerical weather prediction (NWP) mode. Seasonal experiments with the default setup are able to reproduce the basic climatology of the stratospheric polar vortex. However, westerlies are too weak and major stratospheric warmings too frequent in ICON. Both a reduction of the nonorographic, and a reduction of the orographic gravity wave and wake drag lead to a strengthening of the stratospheric vortex and a bias reduction, in particular in January. However, the effect of the nonorographic gravity wave drag scheme on the stratosphere is stronger. Stratosphere-troposphere coupling is intensified and more realistic due to a reduced gravity wave drag. Furthermore, an adjustment of the subgrid-scale orographic drag parameterization leads to a significant error reduction in the mean sea level pressure. As a result of these findings, we present our current suggested improved setup for seasonal experiments with ICON-NWP. <br /> Plain Language Summary Although seasonal forecasts for midlatitudes have the potential to be highly beneficial to the public sector, they are still characterized by a large amount of uncertainty. Exact simulations of the circulation in the stratosphere can help to improve tropospheric predictability on seasonal time scales. For this reason, we investigate how well the new German atmospheric model is able to simulate the stratospheric circulation. The model reproduces the basic behavior of the Northern Hemisphere stratospheric polar vortex, but the westerly circulation in winter is underestimated. The stratospheric circulation is influenced by gravity waves that exert drag on the flow. These processes are only partly physically represented in the model, but are very important and are hence parameterized. By adjusting the parameterizations for the gravity wave drag, the stratospheric polar vortex is strengthened, thereby yielding a more realistic stratospheric circulation. In addition, the altered parameterizations improve the simulated surface pressure pattern. Based upon this, we present our current suggested improved model setup for seasonal experiments.
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.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.
Master equations are a vital tool to model heat flow through nanoscale thermodynamic systems. Most practical devices are made up of interacting subsystems and are often modelled using either local master equations (LMEs) or global master equations (GMEs). While the limiting cases in which either the LME or the GME breaks down are well understood, there exists a 'grey area' in which both equations capture steady-state heat currents reliably but predict very different transient heat flows. In such cases, which one should we trust? Here we show that, when it comes to dynamics, the local approach can be more reliable than the global one for weakly interacting open quantum systems. This is due to the fact that the secular approximation, which underpins the GME, can destroy key dynamical features. To illustrate this, we consider a minimal transport setup and show that its LME displays exceptional points (EPs). These singularities have been observed in a superconducting-circuit realisation of the model [1]. However, in stark contrast to experimental evidence, no EPs appear within the global approach. We then show that the EPs are a feature built into the Redfield equation, which is more accurate than the LME and the GME. Finally, we show that the local approach emerges as the weak-interaction limit of the Redfield equation, and that it entirely avoids the secular approximation.
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.
The Rarita-Schwinger operator is the twisted Dirac operator restricted to 3/2-spinors. Rarita-Schwinger fields are solutions of this operator which are in addition divergence-free. This is an overdetermined problem and solutions are rare; it is even more unexpected for there to be large dimensional spaces of solutions. In this paper we prove the existence of a sequence of compact manifolds in any given dimension greater than or equal to 4 for which the dimension of the space of Rarita-Schwinger fields tends to infinity. These manifolds are either simply connected Kahler-Einstein spin with negative Einstein constant, or products of such spaces with flat tori. Moreover, we construct Calabi-Yau manifolds of even complex dimension with more linearly independent Rarita-Schwinger fields than flat tori of the same dimension.
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.