DOAJ gelistet
Refine
Language
- English (150)
Is part of the Bibliography
- yes (150)
Keywords
- chemotaxis (4)
- diffusion (4)
- radiation belts (3)
- synchronization (3)
- Atlantic meridional overturning circulation (2)
- Ti-6Al-4V (2)
- X-ray computed tomography (CT) (2)
- additive manufacturing (2)
- cell migration (2)
- electrons (2)
Institute
- Institut für Physik und Astronomie (150) (remove)
Leptonic nonthermal emission from supernova remnants evolving in the circumstellar magnetic field
(2022)
The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their nonthermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting nonthermal emission.
Magnetic reconnection is a multi-faceted process of energy conversion in astrophysical, space and laboratory plasmas that operates at microscopic scales but has macroscopic drivers and consequences.
Solar flares present a key laboratory for its study, leaving imprints of the microscopic physics in radiation spectra and allowing the macroscopic evolution to be imaged, yet a full observational characterization remains elusive.
Here we combine high resolution imaging and spectral observations of a confined solar flare at multiple wavelengths with data-constrained magnetohydrodynamic modeling to study the dynamics of the flare plasma from the current sheet to the plasmoid scale. The analysis suggests that the flare resulted from the interaction of a twisted magnetic flux rope surrounding a filament with nearby magnetic loops whose feet are anchored in chromospheric fibrils. Bright cusp-shaped structures represent the region around a reconnecting separator or quasi-separator (hyperbolic flux tube).
The fast reconnection, which is relevant for other astrophysical environments, revealed plasmoids in the current sheet and separatrices and associated unresolved turbulent motions.
Solar flares provide wide range of observational details about fundamental processes involved. Here, the authors show evidence for magnetic reconnection in a strong confined solar flare displaying all four reconnection flows with plasmoids in the current sheet and the separatrices.
The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant H-0 that does not rely on a cosmic distance ladder, nor assumes a specific cosmological model.
By using gravitational waves as "standard sirens", this approach holds promise to arbitrate the existing tension between the H-0 value inferred from the cosmic microwave background and those obtained from local measurements.
However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source led to a H-0 value from GW170817 that was not precise enough to resolve the existing tension.
In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on H-0.
We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.
In the last years, electron density profile functions characterized by a linear dependence on the scale height showed good results when approximating the topside ionosphere. The performance above 800 km, however, is not yet well investigated.
This study investigates the capability of the semi-Epstein functions to represent electron density profiles from the peak height up to 20,000 km. Electron density observations recorded by the Van Allen Probes were used to resolve the scale height dependence in the plasmasphere.
It was found that the linear dependence of the scale height in the topside ionosphere cannot be directly used to extrapolate profiles above 800 km.
We find that the dependence of scale heights on altitude is quadratic in the plasmasphere. A statistical model of the scale heights is therefore proposed. After combining the topside ionosphere and plasmasphere by a unified model, we have obtained good estimations not only in the profile shapes, but also in the Total Electron Content magnitude and distributions when compared to actual measurements from 2013, 2014, 2016 and 2017.
Our investigation shows that Van Allen Probes can be merged to radio-occultation data to properly represent the upper ionosphere and plasmasphere by means of a semi-Epstein function.
The determination of the spin state of iron-bearing compounds at high pressure and temperature is crucial for our understanding of chemical and physical properties of the deep Earth. Studies on the relationship between the coordination of iron and its electronic spin structure in iron-bearing oxides, silicates, carbonates, iron alloys, and other minerals found in the Earth's mantle and core are scarce because of the technical challenges to simultaneously probe the sample at high pressures and temperatures. We used the unique properties of a pulsed and highly brilliant x-ray free electron laser (XFEL) beam at the High Energy Density (HED) instrument of the European XFEL to x-ray heat and probe samples contained in a diamond anvil cell. We heated and probed with the same x-ray pulse train and simultaneously measured x-ray emission and x-ray diffraction of an FeCO3 sample at a pressure of 51 GPa with up to melting temperatures. We collected spin state sensitive Fe K beta(1,3) fluorescence spectra and detected the sample's structural changes via diffraction, observing the inverse volume collapse across the spin transition. During x-ray heating, the carbonate transforms into orthorhombic Fe4C3O12 and iron oxides. Incipient melting was also observed. This approach to collect information about the electronic state and structural changes from samples contained in a diamond anvil cell at melting temperatures and above will considerably improve our understanding of the structure and dynamics of planetary and exoplanetary interiors.
Characterization of binding interactions of SARS-CoV-2 spike protein and DNA-peptide nanostructures
(2022)
Binding interactions of the spike proteins of the severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) to a peptide fragment derived from the human angiotensin converting enzyme 2 (hACE2) receptor are investigated.
The peptide is employed as capture moiety in enzyme linked immunosorbent assays (ELISA) and quantitative binding interaction measurements that are based on fluorescence proximity sensing (switchSENSE).
In both techniques, the peptide is presented on an oligovalent DNA nanostructure, in order to assess the impact of mono- versus trivalent binding modes.
As the analyte, the spike protein and several of its subunits are tested as well as inactivated SARS-CoV-2 and pseudo viruses. While binding of the peptide to the full-length spike protein can be observed, the subunits RBD and S1 do not exhibit binding in the employed concentrations.
Variations of the amino acid sequence of the recombinant full-length spike proteins furthermore influence binding behavior. The peptide was coupled to DNA nanostructures that form a geometric complement to the trimeric structure of the spike protein binding sites.
An increase in binding strength for trimeric peptide presentation compared to single peptide presentation could be generally observed in ELISA and was quantified in switchSENSE measurements. Binding to inactivated wild type viruses could be shown as well as qualitatively different binding behavior of the Alpha and Beta variants compared to the wild type virus strain in pseudo virus models.
Fast-localized electron loss, resulting from interactions with electromagnetic ion cyclotron (EMIC) waves, can produce deepening minima in phase space density (PSD) radial profiles. Here, we perform a statistical analysis of local PSD minima to quantify how readily these are associated with radiation belt depletions. The statistics of PSD minima observed over a year are compared to the Versatile Electron Radiation Belts (VERB) simulations, both including and excluding EMIC waves. The observed minima distribution can only be achieved in the simulation including EMIC waves, indicating their importance in the dynamics of the radiation belts. By analyzing electron flux depletions in conjunction with the observed PSD minima, we show that, in the heart of the outer radiation belt (L* < 5), on average, 53% of multi-MeV electron depletions are associated with PSD minima, demonstrating that fast localized loss by interactions with EMIC waves are a common and crucial process for ultra-relativistic electron populations.
The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search.
We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta.
We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times.
We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient.
When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity.
Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large-scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which comove with the upstream flow and as nonlinear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex nonlinear structures bound by a strong magnetic field around four times the upstream value. Although these nonlinear structures do not in general interact with cospatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around 0.8% of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately three times more efficient than, stochastic shock drift acceleration.
Assessing the impact of hydrogen absorption on the characteristics of the Galactic center excess
(2022)
We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation transport modeling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic center.
The main benefits of the new gas model are (a) the ability to reproduce the negative line signals seen with the HI4PI survey and (b) the accounting for gas that primarily manifests itself through absorption.
We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter.
We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H i model is used to estimate the cosmic-ray induced diffuse gamma-ray emission.
The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence of a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favored by the data.
This finding is robust under the variation of various parameters, for example, the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.
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.
The current paradigm of cosmic-ray (CR) origin states that the greater part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after a supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. With this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely to be triggered at high Mach number shocks.
A detailed investigation of the energy levels of perylene-3,4,9,10-tetracarboxylic tetraethylester as a representative compound for the whole family of perylene esters was performed. It was revealed via electrochemical measurements that one oxidation and two reductions take place. The bandgaps determined via the electrochemical approach are in good agreement with the optical bandgap obtained from the absorption spectra via a Tauc plot. In addition, absorption spectra in dependence of the electrochemical potential were the basis for extensive quantum-chemical calculations of the neutral, monoanionic, and dianionic molecules. For this purpose, calculations based on density functional theory were compared with post-Hartree-Fock methods and the CAM-B3LYP functional proved to be the most reliable choice for the calculation of absorption spectra. Furthermore, spectral features found experimentally could be reproduced with vibronic calculations and allowed to understand their origins. In particular, the two lowest energy absorption bands of the anion are not caused by absorption of two distinct electronic states, which might have been expected from vertical excitation calculations, but both states exhibit a strong vibronic progression resulting in contributions to both bands.
The most complex but potentially most severe impacts of climate change are caused by extreme weather events. In a globally connected economy, damages can cause remote perturbations and cascading consequences-a ripple effect along supply chains. Here we show an economic ripple resonance that amplifies losses when consecutive or overlapping weather extremes and their repercussions interact. This amounts to an average amplification of 21% for climate-induced heat stress, river floods, and tropical cyclones. Modeling the temporal evolution of 1.8 million trade relations between >7000 regional economic sectors, we find that the regional responses to future extremes are strongly heterogeneous also in their resonance behavior. The induced effect on welfare varies between gains due to increased demand in some regions and losses due to demand or supply shortages in others. Within the current global supply network, the ripple resonance effect of extreme weather is strongest in high-income economies-an important effect to consider when evaluating past and future economic climate impacts.
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.
This paper studies cosmic-ray (CR) transport in magnetohydrodynamic (MHD) turbulence. CR transport is strongly dependent on the properties of the magnetic turbulence.
We perform test particle simulations to study the interactions of CR with both total MHD turbulence and decomposed MHD modes.
The spatial diffusion coefficients and the pitch angle scattering diffusion coefficients are calculated from the test particle trajectories in turbulence.
Our results confirm that the fast modes dominate the CR propagation, whereas Alfven and slow modes are much less efficient and have shown similar pitch-angle scattering rates.
We investigate the cross field transport on large and small scales. On large/global scales, normal diffusion is observed and the diffusion coefficient is suppressed by M-A(zeta) compared to the parallel diffusion coefficients, with zeta closer to 4 in Alfven modes than that in total turbulence, as theoretically expected.
For the CR transport on scales smaller than the turbulence injection scale, both the local and global magnetic reference frames are adopted. Superdiffusion is observed on such small scales in all the cases. Particularly, CR transport in Alfven modes show clear Richardson diffusion in the local reference frame. The diffusion transitions smoothly from the Richardson's one with index 1.5 to normal diffusion as the particle mean free path decreases from lambda(parallel to) >> L to lambda(parallel to) << L, where L is the injection/coherence length of turbulence.
Our results have broad applications to CRs in various astrophysical environments.
Cirrus is the only cloud type capable of inducing daytime cooling or heating at the top of the atmosphere (TOA) and the sign of its radiative effect highly depends on its optical depth. However, the investigation of its geometrical and optical properties over the Arctic is limited. In this work the long-term properties of cirrus clouds are explored for the first time over an Arctic site (Ny-Alesund, Svalbard) using lidar and radiosonde measurements from 2011 to 2020. The optical properties were quality assured, taking into account the effects of specular reflections and multiple-scattering. Cirrus clouds were generally associated with colder and calmer wind conditions compared to the 2011-2020 climatology. However, the dependence of cirrus properties on temperature and wind speed was not strong. Even though the seasonal cycle was not pronounced, the winter-time cirrus appeared under lower temperatures and stronger wind conditions. Moreover, in winter, geometrically- and optically-thicker cirrus were found and their ice particles tended to be more spherical. The majority of cirrus was associated with westerly flow and westerly cirrus tended to be geometrically-thicker. Overall, optically-thinner layers tended to comprise smaller and less spherical ice crystals, most likely due to reduced water vapor deposition on the particle surface. Compared to lower latitudes, the cirrus layers over Ny-Alesund were more absorbing in the visible spectral region and they consisted of more spherical ice particles.
Complex networks are abundant in nature and many share an important structural property: they contain a few nodes that are abnormally highly connected (hubs). Some of these hubs are called influencers because they couple strongly to the network and play fundamental dynamical and structural roles. Strikingly, despite the abundance of networks with influencers, little is known about their response to stochastic forcing. Here, for oscillatory dynamics on influencer networks, we show that subjecting influencers to an optimal intensity of noise can result in enhanced network synchronization. This new network dynamical effect, which we call coherence resonance in influencer networks, emerges from a synergy between network structure and stochasticity and is highly nonlinear, vanishing when the noise is too weak or too strong. Our results reveal that the influencer backbone can sharply increase the dynamical response in complex systems of coupled oscillators. Influencer networks include a small set of highly-connected nodes and can reach synchrony only via strong node interaction. Tonjes et al. show that introducing an optimal amount of noise enhances synchronization of such networks, which may be relevant for neuroscience or opinion dynamics applications.
Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea level rise. Increased surface melt may lower the reflectivity of the ice sheet surface and thereby increase melt rates: the so-called melt-albedo feedback describes this self-sustaining increase in surface melting. In order to test the effect of the melt-albedo feedback in a prognostic ice sheet model, we implement dEBM-simple, a simplified version of the diurnal Energy Balance Model dEBM, in the Parallel Ice Sheet Model (PISM). The implementation includes a simple representation of the melt-albedo feedback and can thereby replace the positive-degree-day melt scheme. Using PISM-dEBM-simple, we find that this feedback increases ice loss through surface warming by 60 % until 2300 for the high-emission scenario RCP8.5 when compared to a scenario in which the albedo remains constant at its present-day values. With an increase of 90 % compared to a fixed-albedo scenario, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be 70 % in the RCP8.5 scenario and a more than 4-fold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt-albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming.
Starting from the observation that the reduced state of a system strongly coupled to a bath is, in general, an athermal state, we introduce and study a cyclic battery-charger quantum device that is in thermal equilibrium, or in a ground state, during the charge storing stage. The cycle has four stages: the equilibrium storage stage is interrupted by disconnecting the battery from the charger, then work is extracted from the battery, and then the battery is reconnected with the charger; finally, the system is brought back to equilibrium. At no point during the cycle are the battery-charger correlations artificially erased. We study the case where the battery and charger together comprise a spin-1/2 Ising chain, and show that the main characteristics-the extracted energy and the thermodynamic efficiency-can be enhanced by operating the cycle close to the quantum phase transition point. When the battery is just a single spin, we find that the output work and efficiency show a scaling behavior at criticality and derive the corresponding critical exponents. Due to always present correlations between the battery and the charger, operations that are equivalent from the perspective of the battery can entail different energetic costs for switching the battery-charger coupling. This happens only when the coupling term does not commute with the battery's bare Hamiltonian, and we use this purely quantum leverage to further optimize the performance of the device.
Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet, the melting increases the ice discharge into the ocean.
Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an isothermal, inherently buttressed ice-sheet-shelf system. In the idealized numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's (1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and (2) shear margins, where the ice flow is slower.
The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction.
Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.
The application of the fractional calculus in the mathematical modelling of relaxation processes in complex heterogeneous media has attracted a considerable amount of interest lately.
The reason for this is the successful implementation of fractional stochastic and kinetic equations in the studies of non-Debye relaxation.
In this work, we consider the rotational diffusion equation with a generalised memory kernel in the context of dielectric relaxation processes in a medium composed of polar molecules. We give an overview of existing models on non-exponential relaxation and introduce an exponential resetting dynamic in the corresponding process.
The autocorrelation function and complex susceptibility are analysed in detail.
We show that stochastic resetting leads to a saturation of the autocorrelation function to a constant value, in contrast to the case without resetting, for which it decays to zero. The behaviour of the autocorrelation function, as well as the complex susceptibility in the presence of resetting, confirms that the dielectric relaxation dynamics can be tuned by an appropriate choice of the resetting rate.
The presented results are general and flexible, and they will be of interest for the theoretical description of non-trivial relaxation dynamics in heterogeneous systems composed of polar molecules.
The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth.
Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model.
Our model results show that for warming levels above 2 degrees C, Greenland could become essentially ice-free within several millennia, mainly as a result of surface melting and acceleration of ice flow. These ice losses are mitigated, however, in some cases with strong GIA feedback even promoting an incomplete recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming - from incomplete recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse.
In the incomplete recovery regime, the initial ice loss due to warming is essentially reversed within 50 000 years, and the ice volume stabilizes at 61 %-93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods between 74 000 and over 300 000 years and oscillation amplitudes between 15 %-70 % of present-day ice volume.
This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on timescales on the order of 100 000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet.
Point-of-care and in-vivo bio-diagnostic tools are the current need for the present critical scenarios in the healthcare industry. The past few decades have seen a surge in research activities related to solving the challenges associated with precise on-site bio-sensing. Cutting-edge fiber optic technology enables the interaction of light with functionalized fiber surfaces at remote locations to develop a novel, miniaturized and cost-effective lab on fiber technology for bio-sensing applications. The recent remarkable developments in the field of nanotechnology provide innumerable functionalization methodologies to develop selective bio-recognition elements for label free biosensors. These exceptional methods may be easily integrated with fiber surfaces to provide highly selective light-matter interaction depending on various transduction mechanisms. In the present review, an overview of optical fiber-based biosensors has been provided with focus on physical principles used, along with the functionalization protocols for the detection of various biological analytes to diagnose the disease. The design and performance of these biosensors in terms of operating range, selectivity, response time and limit of detection have been discussed. In the concluding remarks, the challenges associated with these biosensors and the improvement required to develop handheld devices to enable direct target detection have been highlighted.
Fluorination of organic spacer impacts on the structural and optical response of 2D perovskites
(2020)
Low-dimensional hybrid perovskites have triggered significant research interest due to their intrinsically tunable optoelectronic properties and technologically relevant material stability. In particular, the role of the organic spacer on the inherent structural and optical features in two-dimensional (2D) perovskites is paramount for material optimization. To obtain a deeper understanding of the relationship between spacers and the corresponding 2D perovskite film properties, we explore the influence of the partial substitution of hydrogen atoms by fluorine in an alkylammonium organic cation, resulting in (Lc)(2)PbI4 and (Lf)(2)PbI4 2D perovskites, respectively. Consequently, optical analysis reveals a clear 0.2 eV blue-shift in the excitonic position at room temperature. This result can be mainly attributed to a band gap opening, with negligible effects on the exciton binding energy. According to Density Functional Theory (DFT) calculations, the band gap increases due to a larger distortion of the structure that decreases the atomic overlap of the wavefunctions and correspondingly bandwidth of the valence and conduction bands. In addition, fluorination impacts the structural rigidity of the 2D perovskite, resulting in a stable structure at room temperature and the absence of phase transitions at a low temperature, in contrast to the widely reported polymorphism in some non-fluorinated materials that exhibit such a phase transition. This indicates that a small perturbation in the material structure can strongly influence the overall structural stability and related phase transition of 2D perovskites, making them more robust to any phase change. This work provides key information on how the fluorine content in organic spacer influence the structural distortion of 2D perovskites and their optical properties which possess remarkable importance for future optoelectronic applications, for instance in the field of light-emitting devices or sensors.
Both ice sheets in Greenland and Antarctica are discharging ice into the ocean. In many regions along the coast of the ice sheets, the icebergs calve into a bay. If the addition of icebergs through calving is faster than their transport out of the embayment, the icebergs will be frozen into a melange with surrounding sea ice in winter. In this case, the buttressing effect of the ice melange can be considerably stronger than any buttressing by mere sea ice would be. This in turn stabilizes the glacier terminus and leads to a reduction in calving rates. Here we propose a simple parametrization of ice melange buttressing which leads to an upper bound on calving rates and can be used in numerical and analytical modelling.
In classical thermodynamic processes the unavoidable presence of irreversibility, quantified by the entropy production, carries two energetic footprints: the reduction of extractable work from the optimal, reversible case, and the generation of a surplus of heat that is irreversibly dissipated to the environment. Recently it has been shown that in the quantum regime an additional quantum irreversibility occurs that is linked to decoherence into the energy basis. Here we employ quantum trajectories to construct distributions for classical heat and quantum heat exchanges, and show that the heat footprint of quantum irreversibility differs markedly from the classical case. We also quantify how quantum irreversibility reduces the amount of work that can be extracted from a state with coherences. Our results show that decoherence leads to both entropic and energetic footprints which both play an important role in the optimization of controlled quantum operations at low temperature.
In classical thermodynamics irreversibility occurs whenever a non-thermal system is brought into contact with a thermal environment. Using quantum trajectories the authors here establish two energetic footprints of quantum irreversible processes, and find that while quantum irreversibility leads to the occurrence of a quantum heat and a reduction of work production, the two are not linked in the same manner as the classical laws of thermodynamics would dictate.
Estimating global mean sea-level rise and its uncertainties by 2100 and 2300 from an expert survey
(2020)
Sea-level rise projections and knowledge of their uncertainties are vital to make informed mitigation and adaptation decisions. To elicit projections from members of the scientific community regarding future global mean sea-level (GMSL) rise, we repeated a survey originally conducted five years ago. Under Representative Concentration Pathway (RCP) 2.6, 106 experts projected a likely (central 66% probability) GMSL rise of 0.30-0.65 m by 2100, and 0.54-2.15 m by 2300, relative to 1986-2005. Under RCP 8.5, the same experts projected a likely GMSL rise of 0.63-1.32 m by 2100, and 1.67-5.61 m by 2300. Expert projections for 2100 are similar to those from the original survey, although the projection for 2300 has extended tails and is higher than the original survey. Experts give a likelihood of 42% (original survey) and 45% (current survey) that under the high-emissions scenario GMSL rise will exceed the upper bound (0.98 m) of the likely range estimated by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, which is considered to have an exceedance likelihood of 17%. Responses to open-ended questions suggest that the increases in upper-end estimates and uncertainties arose from recent influential studies about the impact of marine ice cliff instability on the meltwater contribution to GMSL rise from the Antarctic Ice Sheet.
To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB Lambda phi. Overall, our test simulations have good agreement with those reported in the literature.
Diffraction enhanced imaging (DEI) is an advanced digital radiographic imaging technique employing the refraction of X-rays to contrast internal interfaces. This study aims to qualitatively and quantitatively evaluate images acquired using this technique and to assess how different fitting functions to the typical rocking curves (RCs) influence the quality of the images. RCs are obtained for every image pixel. This allows the separate determination of the absorption and the refraction properties of the material in a position-sensitive manner. Comparison of various types of fitting functions reveals that the Pseudo-Voigt (PsdV) function is best suited to fit typical RCs. A robust algorithm was developed in the Python programming language, which reliably extracts the physically meaningful information from each pixel of the image. We demonstrate the potential of the algorithm with two specimens: a silicone gel specimen that has well-defined interfaces, and an additively manufactured polycarbonate specimen.
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.
To undergo diffusive shock acceleration, electrons need to be preaccelerated to increase their energies by several orders of magnitude, else their gyroradii will be smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel nor perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the preacceleration in this region, we undertake particle-in-cell simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for theta (Bn) = 74.degrees 3, it increases to R similar to 5% for theta (Bn) = 30 degrees, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to R (0.6) and a wavelength approximate to 2 lambda (se), where lambda (se) is the electron skin length. While the initial reflection mechanism is typically a combination of shock-surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream, they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In less than or similar to 1% of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction is rescattered back toward the shock after reflection, reinteracts with the shock with energies much greater than thermal, and crosses into the downstream.
Increasing greenhouse gas emissions are likely to impact not only natural systems but economies worldwide. If these impacts alter future economic development, the financial losses will be significantly higher than the mere direct damages. So far, potentially aggravating investment responses were considered negligible. Here we consistently incorporate an empirically derived temperature-growth relation into the simple integrated assessment model DICE. In this framework we show that, if in the next eight decades varying temperatures impact economic growth as has been observed in the past three decades, income is reduced by similar to 20% compared to an economy unaffected by climate change. Hereof similar to 40% are losses due to growth effects of which similar to 50% result from reduced incentive to invest. This additional income loss arises from a reduced incentive for future investment in anticipation of a reduced return and not from an explicit climate protection policy. Under economically optimal climate-change mitigation, however, optimal investment would only be reduced marginally as mitigation efforts keep returns high.
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.
Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice melange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice melange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-melange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.
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.
In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of approximate to 11% initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.
In this study we analyze the storm-time evolution of equatorial electron pitch angle distributions (PADs) in the outer radiation belt region using observations from the Magnetic Electron Ion Spectrometer (MagEIS) instrument aboard the Van Allen Probes in 2012-2019. The PADs are approximated using a sum of the first, third and fifth sine harmonics. Different combinations of the respective coefficients refer to the main PAD shapes within the outer radiation belt, namely the pancake, flat-top, butterfly and cap PADs. We conduct a superposed epoch analysis of 129 geomagnetic storms and analyze the PAD evolution for day and night MLT sectors. PAD shapes exhibit a strong energy-dependent response. At energies of tens of keV, the PADs exhibit little variation throughout geomagnetic storms. Cap PADs are mainly observed at energies < 300 keV, and their extent in L shrinks with increasing energy. The cap distributions transform into the pancake PADs around the main phase of the storm on the nightside, and then come back to their original shapes during the recovery phase. At higher energies on the dayside, the PADs are mainly pancake during pre-storm conditions and become more anisotropic during the main phase. The quiet-time butterfly PADs can be observed on the nightside at L> 5.6. During the main phase, butterfly PADs have stronger 90 degrees-minima and can be observed at lower L-shells (down to L = 5), then transitioning into flat-top PADs at L similar to 4.5 - 5 and pancake PADs at L < 4.5. The resulting PAD coefficients for different energies, locations and storm epochs can be used to test the wave models and physics-based radiation belt codes in terms of pitch angle distributions.
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.
The remarkable progress of metal halide perovskites in photovoltaics has led to the power conversion efficiency approaching 26%. However, practical applications of perovskite-based solar cells are challenged by the stability issues, of which the most critical one is photo-induced degradation. Bare CH3NH3PbI3 perovskite films are known to decompose rapidly, with methylammonium and iodine as volatile species and residual solid PbI2 and metallic Pb, under vacuum under white light illumination, on the timescale of minutes. We find, in agreement with previous work, that the degradation is non-uniform and proceeds predominantly from the surface, and that illumination under N-2 and ambient air (relative humidity 20%) does not induce substantial degradation even after several hours. Yet, in all cases the release of iodine from the perovskite surface is directly identified by X-ray photoelectron spectroscopy. This goes in hand with a loss of organic cations and the formation of metallic Pb. When CH3NH3PbI3 films are covered with a few nm thick organic capping layer, either charge selective or non-selective, the rapid photodecomposition process under ultrahigh vacuum is reduced by more than one order of magnitude, and becomes similar in timescale to that under N-2 or air. We conclude that the light-induced decomposition reaction of CH3NH3PbI3, leading to volatile methylammonium and iodine, is largely reversible as long as these products are restrained from leaving the surface. This is readily achieved by ambient atmospheric pressure, as well as a thin organic capping layer even under ultrahigh vacuum. In addition to explaining the impact of gas pressure on the stability of this perovskite, our results indicate that covalently "locking" the position of perovskite components at the surface or an interface should enhance the overall photostability.
Tropical cyclones range among the costliest disasters on Earth. Their economic repercussions along the supply and trade network also affect remote economies that are not directly affected. We here simulate possible global repercussions on consumption for the example case of Hurricane Sandy in the US (2012) using the shock-propagation model Acclimate. The modeled shock yields a global three-phase ripple: an initial production demand reduction and associated consumption price decrease, followed by a supply shortage with increasing prices, and finally a recovery phase. Regions with strong trade relations to the US experience strong magnitudes of the ripple. A dominating demand reduction or supply shortage leads to overall consumption gains or losses of a region, respectively. While finding these repercussions in historic data is challenging due to strong volatility of economic interactions, numerical models like ours can help to identify them by approaching the problem from an exploratory angle, isolating the effect of interest. For this, our model simulates the economic interactions of over 7000 regional economic sectors, interlinked through about 1.8 million trade relations. Under global warming, the wave-like structures of the economic response to major hurricanes like the one simulated here are likely to intensify and potentially overlap with other weather extremes.
Due to climate change the frequency and character of precipitation are changing as the hydrological cycle intensifies. With regards to snowfall, global warming has two opposing influences; increasing humidity enables intense snowfall, whereas higher temperatures decrease the likelihood of snowfall. Here we show an intensification of extreme snowfall across large areas of the Northern Hemisphere under future warming. This is robust across an ensemble of global climate models when they are bias-corrected with observational data. While mean daily snowfall decreases, both the 99th and the 99.9th percentiles of daily snowfall increase in many regions in the next decades, especially for Northern America and Asia. Additionally, the average intensity of snowfall events exceeding these percentiles as experienced historically increases in many regions. This is likely to pose a challenge to municipalities in mid to high latitudes. Overall, extreme snowfall events are likely to become an increasingly important impact of climate change in the next decades, even if they will become rarer, but not necessarily less intense, in the second half of the century.
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.
Master equations are a vital tool to model heat flow through nanoscale thermodynamic systems. Most practical devices are made up of interacting subsystems and are often modelled using either local master equations (LMEs) or global master equations (GMEs). While the limiting cases in which either the LME or the GME breaks down are well understood, there exists a 'grey area' in which both equations capture steady-state heat currents reliably but predict very different transient heat flows. In such cases, which one should we trust? Here we show that, when it comes to dynamics, the local approach can be more reliable than the global one for weakly interacting open quantum systems. This is due to the fact that the secular approximation, which underpins the GME, can destroy key dynamical features. To illustrate this, we consider a minimal transport setup and show that its LME displays exceptional points (EPs). These singularities have been observed in a superconducting-circuit realisation of the model [1]. However, in stark contrast to experimental evidence, no EPs appear within the global approach. We then show that the EPs are a feature built into the Redfield equation, which is more accurate than the LME and the GME. Finally, we show that the local approach emerges as the weak-interaction limit of the Redfield equation, and that it entirely avoids the secular approximation.
We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
The Parallel Ice Sheet Model (PISM) is applied to the Antarctic Ice Sheet over the last two glacial cycles (approximate to 210 000 years) with a resolution of 16 km. An ensemble of 256 model runs is analyzed in which four relevant model parameters have been systematically varied using full-factorial parameter sampling. Parameters and plausible parameter ranges have been identified in a companion paper (Albrecht et al., 2020) and are associated with ice dynamics, climatic forcing, basal sliding and bed deformation and represent distinct classes of model uncertainties. The model is scored against both modern and geologic data, including reconstructed grounding-line locations, elevation-age data, ice thickness, surface velocities and uplift rates. An aggregated score is computed for each ensemble member that measures the overall model-data misfit, including measurement uncertainty in terms of a Gaussian error model (Briggs and Tarasov, 2013). The statistical method used to analyze the ensemble simulation results follows closely the simple averaging method described in Pollard et al. (2016).
This analysis reveals clusters of best-fit parameter combinations, and hence a likely range of relevant model and boundary parameters, rather than individual best-fit parameters. The ensemble of reconstructed histories of Antarctic Ice Sheet volumes provides a score-weighted likely range of sea-level contributions since the Last Glacial Maximum (LGM) of 9.4 +/- 4.1m (or 6.5 +/- 2.0 x 10(6) km(3)), which is at the upper range of most previous studies. The last deglaciation occurs in all ensemble simulations after around 12 000 years before present and hence after the meltwater pulse 1A (MWP1a). Our ensemble analysis also provides an estimate of parametric uncertainty bounds for the present-day state that can be used for PISM projections of future sea-level contributions from the Antarctic Ice Sheet.
Neutrophil granulocytes are essential for the first host defense. After leaving the blood circulation they migrate efficiently towards sites of inflammation. They are guided by chemoattractants released from cells within the inflammatory foci. On a cellular level, directional migration is a consequence of cellular front-rear asymmetry which is induced by the concentration gradient of the chemoattractants. The generation and maintenance of this asymmetry, however, is not yet fully understood. Here we analyzed the paths of chemotacting neutrophils with different stochastic models to gain further insight into the underlying mechanisms. Wildtype chemotacting neutrophils show an anomalous superdiffusive behavior. CXCR2 blockade and TRPC6-knockout cause the tempering of temporal correlations and a reduction of chemotaxis. Importantly, such tempering is found both in vitro and in vivo. These findings indicate that the maintenance of anomalous dynamics is crucial for chemotactic behavior and the search efficiency of neutrophils.
The motility of neutrophils and their ability to sense and to react to chemoattractants in their environment are of central importance for the innate immunity. Neutrophils are guided towards sites of inflammation following the activation of G-protein coupled chemoattractant receptors such as CXCR2 whose signaling strongly depends on the activity of Ca2+ permeable TRPC6 channels. It is the aim of this study to analyze data sets obtained in vitro (murine neutrophils) and in vivo (zebrafish neutrophils) with a stochastic mathematical model to gain deeper insight into the underlying mechanisms. The model is based on the analysis of trajectories of individual neutrophils. Bayesian data analysis, including the covariances of positions for fractional Brownian motion as well as for exponentially and power-law tempered model variants, allows the estimation of parameters and model selection. Our model-based analysis reveals that wildtype neutrophils show pure superdiffusive fractional Brownian motion. This so-called anomalous dynamics is characterized by temporal long-range correlations for the movement into the direction of the chemotactic CXCL1 gradient. Pure superdiffusion is absent vertically to this gradient. This points to an asymmetric 'memory' of the migratory machinery, which is found both in vitro and in vivo. CXCR2 blockade and TRPC6-knockout cause tempering of temporal correlations in the chemotactic gradient. This can be interpreted as a progressive loss of memory, which leads to a marked reduction of chemotaxis and search efficiency of neutrophils. In summary, our findings indicate that spatially differential regulation of anomalous dynamics appears to play a central role in guiding efficient chemotactic behavior.
Perovskite solar cells represent one of the recent success stories in photovoltaics. The device efficiency has been steadily increasing over the past years, but further work is needed to enhance the performance, for example, through the reduction of defects to prevent carrier recombination. SCAPS-1D simulations were performed to assess efficiency limits and identify approaches to decrease the impact of defects, through the selection of an optimal hole-transport material and a hole-collecting electrode. Particular attention was given to evaluation of the influence of bulk defects within light-absorbing CH3NH3SnI3 layers. In addition, the study demonstrates the influence of interface defects at the TiO2/CH3NH3SnI3 (IL1) and CH3NH3SnI3/HTL (IL2) interfaces across the similar range of defect densities. Finally, the optimal device architecture TiO2/CH3NH3SnI3/Cu2O is proposed for the given absorber layer using the readily available Cu2O hole-transporting material with PCE = 27.95%, FF = 84.05%, V-OC = 1.02 V and J(SC) = 32.60 mA/cm(2), providing optimal performance and enhanced resistance to defects.
Focusing on transient chaos
(2022)
Recent advances in the field of complex, transiently chaotic dynamics are reviewed, based on the results published in the focus issue of J. Phys. Complex. on this topic. One group of achievements concerns network dynamics where transient features are intimately related to the degree and stability of synchronization, as well as to the network topology. A plethora of various applications of transient chaos are described, ranging from the collective motion of active particles, through the operation of power grids, cardiac arrhythmias, and magnetohydrodynamical dynamos, to the use of machine learning to predict time evolutions. Nontraditional forms of transient chaos are also explored, such as the temporal change of the chaoticity in the transients (called doubly transient chaos), as well as transients in systems subjected to parameter drift, the paradigm of which is climate change.
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.
The Dirac point of a topological surface state (TSS) is protected against gapping by time-reversal symmetry. Conventional wisdom stipulates, therefore, that only through magnetisation may a TSS become gapped. However, non-magnetic gaps have now been demonstrated in Bi2Se3 systems doped with Mn or In, explained by hybridisation of the Dirac cone with induced impurity resonances. Recent photoemission experiments suggest that an analogous mechanism applies even when Bi2Se3 is surface dosed with Au. Here, we perform a systematic spin- and angle-resolved photoemission study of Au-dosed Bi2Se3. Although there are experimental conditions wherein the TSS appears gapped due to unfavourable photoemission matrix elements, our photon-energy-dependent spectra unambiguously demonstrate the robustness of the Dirac cone against high Au coverage. We further show how the spin textures of the TSS and its accompanying surface resonances remain qualitatively unchanged following Au deposition, and discuss the mechanism underlying the suppression of the spectral weight.
Suppression of the TeV Pair-beam-Plasma Instability by a Tangled Weak Intergalactic Magnetic Field
(2022)
We study the effect of a tangled sub-fG level intergalactic magnetic field (IGMF) on the electrostatic instability of a blazar-induced pair beam. Sufficiently strong IGMF may significantly deflect the TeV pair beams, which would reduce the flux of secondary cascade emission below the observational limits. A similar flux reduction may result from the electrostatic beam-plasma instability, which operates the best in the absence of IGMF. Considering IGMF with correlation lengths smaller than a kiloparsec, we find that weak magnetic fields increase the transverse momentum of the pair-beam particles, which dramatically reduces the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. We show that the beam-plasma instability is eliminated as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission.
In this work, the fabrication and characterization of a simple, inexpensive, and effective microfluidic paper analytic device (mu PAD) for monitoring DNA samples is reported. The glass microfiber-based chip has been fabricated by a new wax-based transfer-printing technique and an electrode printing process. It is capable of moving DNA effectively in a time-dependent fashion. The nucleic acid sample is not damaged by this process and is accumulated in front of the anode, but not directly on the electrode. Thus, further DNA processing is feasible. The system allows the DNA to be purified by separating it from other components in sample mixtures such as proteins. Furthermore, it is demonstrated that DNA can be moved through several layers of the glass fiber material. This proof of concept will provide the basis for the development of rapid test systems, e.g., for the detection of pathogens in water samples.
Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the first binary neutron star simulations of this kind. We find that the presence of dark matter reduces the lifetime of the merger remnant and favors a prompt collapse to a black hole. Furthermore, we find differences in the merger time for systems with the same total mass and mass ratio, but different amounts of dark matter. Finally, we find that electromagnetic signals produced by the merger of binary neutron stars admixed with dark matter are very unlikely to be as bright as their dark matter-free counterparts. Given the increased sensitivity of multi-messenger facilities, our analysis gives a new perspective on how to probe the presence of dark matter.
Corona and the climate
(2020)
The detection of internal irregularities is crucial for quality assessment in metal-based additive manufacturing (AM) technologies such as laser powder bed fusion (L-PBF). The utilization of in-process thermography as an in situ monitoring tool in combination with post-process X-ray micro computed tomography (XCT) as a reference technique has shown great potential for this aim. Due to the small irregularity dimensions, a precise registration of the datasets is necessary as a requirement for correlation. In this study, the registration of thermography and XCT reference datasets of a cylindric specimen containing keyhole pores is carried out for the development of a porosity prediction model. The considered datasets show variations in shape, data type and dimensionality, especially due to shrinkage and material elevation effects present in the manufactured part. Since the resulting deformations are challenging for registration, a novel preprocessing methodology is introduced that involves an adaptive volume adjustment algorithm which is based on the porosity distribution in the specimen. Thus, the implementation of a simple three-dimensional image-to-image registration is enabled. The results demonstrate the influence of the part deformation on the resulting porosity location and the importance of registration in terms of irregularity prediction.
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.
Tailoring the secondary surface morphology of electro-spun nanofibers has been highly desired, as such delicate structures equip nanofibers with distinct functions. Here, we report a simple strategy to directly reconstruct the surface of polyvinyl alcohol/polyvinylpyrrolidone (PVA/PVP) nanofibers by water evaporation. The roughness and diameter of the nanofibers depend on the temperature during vacuum drying. Surface changes of the nanofibers from smooth to rough were observed at 55 degrees C, with a significant drop in nanofiber diameter. We attribute the formation of the secondary surface morphology to the intermolecular forces in the water vapor, including capillary and the compression forces, on the basis of the results from the Fourier-transform infrared (FTIR) and X-ray photoelectron (XPS) spectroscopy. The strategy is universally effective for various electro-spun polymer nanofibers, thus opening up avenues toward more detailed and sophisticated structure design and implementation for nanofibers.
We present observations of three-dimensional magnetic power spectra in wavevector space to investigate the anisotropy and scalings of sub-Alfvenic solar wind turbulence at magnetohydrodynamic (MHD) scale using the Magnetospheric Multiscale spacecraft. The magnetic power distributions are organized in a new coordinate determined by wavevectors ((kappa) over cap) and background magnetic field ((b) over cap (0)) in Fourier space. This study utilizes two approaches to determine wavevectors: the singular value decomposition method and multispacecraft timing analysis. The combination of the two methods allows an examination of the properties of magnetic field fluctuations in terms of mode compositions without any spatiotemporal hypothesis. Observations show that fluctuations (delta B-perpendicular to 1) in the direction perpendicular to (kappa) over cap and (b) over cap (0) prominently cascade perpendicular to (b) over cap (0), and such anisotropy increases with wavenumbers. The reduced power spectra of 6.8 11 follow Goldreich-Sridhar scalings: (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-5/3) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-2). In contrast, fluctuations within the (k) over cap(b) over cap (0) plane show isotropic behaviors: perpendicular power distributions are approximately the same as parallel distributions. The reduced power spectra of fluctuations within the (k) over cap(b) over cap (0) plane follow the scalings (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-3/2) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-3/2). Comparing frequency-wavevector spectra with theoretical dispersion relations of MHD modes, we find that delta B-perpendicular to 1 are probably associated with Alfven modes. On the other hand, magnetic field fluctuations within the (k) over cap(b) over cap (0) plane more likely originate from fast modes based on their isotropic behaviors. The observations of anisotropy and scalings of different magnetic field components are consistent with the predictions of current compressible MHD theory. Moreover, for the Alfvenic component, the ratio of cascading time to the wave period is found to be a factor of a few, consistent with critical balance in the strong turbulence regime. These results are valuable for further studies of energy compositions of plasma turbulence and their effects on energetic particle transport.
We derive modified reflection coefficients for electromagnetic waves in the THz and far infrared range. The idea is based on hydrodynamic boundary conditions for metallic conduction electrons. The temperature-dependent part of the Casimir pressure between metal plates is evaluated. The results should shed light on the "thermal anomaly," where measurements deviate from the standard fluctuation electrodynamics for conducting metals.
The Paris Climate Agreement aims to keep temperature rise well below 2 degrees C. This implies mitigation costs as well as avoided climate damages. Here we show that independent of the normative assumptions of inequality aversion and time preferences, the agreement constitutes the economically optimal policy pathway for the century. To this end we consistently incorporate a damage-cost curve reproducing the observed relation between temperature and economic growth into the integrated assessment model DICE. We thus provide an intertemporally optimizing cost-benefit analysis of this century's climate problem. We account for uncertainties regarding the damage curve, climate sensitivity, socioeconomic future, and mitigation costs. The resulting optimal temperature is robust as can be understood from the generic temperature-dependence of the mitigation costs and the level of damages inferred from the observed temperature-growth relationship. Our results show that the politically motivated Paris Climate Agreement also represents the economically favourable pathway, if carried out properly.
Several large-scale cryosphere elements such as the Arctic summer sea ice, the mountain glaciers, the Greenland and West Antarctic Ice Sheet have changed substantially during the last century due to anthropogenic global warming. However, the impacts of their possible future disintegration on global mean temperature (GMT) and climate feedbacks have not yet been comprehensively evaluated. Here, we quantify this response using an Earth system model of intermediate complexity. Overall, we find a median additional global warming of 0.43 degrees C (interquartile range: 0.39-0.46 degrees C) at a CO2 concentration of 400 ppm. Most of this response (55%) is caused by albedo changes, but lapse rate together with water vapour (30%) and cloud feedbacks (15%) also contribute significantly. While a decay of the ice sheets would occur on centennial to millennial time scales, the Arctic might become ice-free during summer within the 21st century. Our findings imply an additional increase of the GMT on intermediate to long time scales. The disintegration of cryosphere elements such as the Arctic summer sea ice, mountain glaciers, Greenland and West Antarctica is associated with temperature and radiative feedbacks. In this work, the authors quantify these feedbacks and find an additional global warming of 0.43 degrees C.
The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Pots-dam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3 degrees respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice-ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.
Organic photovoltaics (PV) is an energy-harvesting technology that offers many advantages, such as flexibility, low weight and cost, as well as environmentally benign materials and manufacturing techniques. Despite growth of power conversion efficiencies to around 19 % in the last years, organic PVs still lag behind inorganic PV technologies, mainly due to high losses in open-circuit voltage. Understanding and improving open circuit voltage in organic solar cells is challenging, as it is controlled by the properties of a donor-acceptor interface where the optical excitations are separated into charge carriers. Here, we provide an electrostatic model of a rough donor-acceptor interface and test it experimentally on small molecule PV materials systems. The model provides concise relationships between the open-circuit voltage, photovoltaic gap, charge-transfer state energy, and interfacial morphology. In particular, we show that the electrostatic bias generated across the interface reduces the photovoltaic gap. This negative influence on open-circuit voltage can, however, be circumvented by adjusting the morphology of the donor-acceptor interface.
Organic solar cells, despite their high power conversion efficiencies, suffer from open circuit voltage losses making them less appealing in terms of applications. Here, the authors, supported with experimental data on small molecule photovoltaic cells, relate open circuit voltage to photovoltaic gap, charge-transfer state energy, and donor-acceptor interfacial morphology.
Accurately predicting total electron content (TEC) during geomagnetic storms is still a challenging task for ionospheric models. In this work, a neural-network (NN)-based model is proposed which predicts relative TEC with respect to the preceding 27-day median TEC, during storm time for the European region (with longitudes 30 degrees W-50 degrees E and latitudes 32.5 degrees N-70 degrees N). The 27-day median TEC (referred to as median TEC), latitude, longitude, universal time, storm time, solar radio flux index F10.7, global storm index SYM-H and geomagnetic activity index Hp30 are used as inputs and the output of the network is the relative TEC. The relative TEC can be converted to the actual TEC knowing the median TEC. The median TEC is calculated at each grid point over the European region considering data from the last 27 days before the storm using global ionosphere maps (GIMs) from international GNSS service (IGS) sources. A storm event is defined when the storm time disturbance index Dst drops below 50 nanotesla. The model was trained with storm-time relative TEC data from the time period of 1998 until 2019 (2015 is excluded) and contains 365 storms. Unseen storm data from 33 storm events during 2015 and 2020 were used to test the model. The UQRG GIMs were used because of their high temporal resolution (15 min) compared to other products from different analysis centers. The NN-based model predictions show the seasonal behavior of the storms including positive and negative storm phases during winter and summer, respectively, and show a mixture of both phases during equinoxes. The model's performance was also compared with the Neustrelitz TEC model (NTCM) and the NN-based quiet-time TEC model, both developed at the German Aerospace Agency (DLR). The storm model has a root mean squared error (RMSE) of 3.38 TEC units (TECU), which is an improvement by 1.87 TECU compared to the NTCM, where an RMSE of 5.25 TECU was found. This improvement corresponds to a performance increase by 35.6%. The storm-time model outperforms the quiet-time model by 1.34 TECU, which corresponds to a performance increase by 28.4% from 4.72 to 3.38 TECU. The quiet-time model was trained with Carrington averaged TEC and, therefore, is ideal to be used as an input instead of the GIM derived 27-day median. We found an improvement by 0.8 TECU which corresponds to a performance increase by 17% from 4.72 to 3.92 TECU for the storm-time model using the quiet-time-model predicted TEC as an input compared to solely using the quiet-time model.
Salinization is a well-known problem in agricultural areas worldwide. In the last 20-30 yr, rising salinity in the upper, unconfined aquifer has been observed in the Freepsumer Meer, a grassland near the German North Sea coast. For investigating long-term development of salinity and water balance during 1961-2099, the one-dimensional Soil-Water-Atmosphere-Plant (SWAP) model was set up and calibrated for a soil column in the area. The model setup involves a deep aquifer as the source of salt through upward seepage. In the vertical salt transport equation, dispersion and advection are included. Six different regional outputs of statistical downscaling methods were used as climate scenarios. These comprise different rates of increasing surface temperature and different trends in seasonal rainfall. The simulation results exhibit opposing salinity trends for topsoil and deeper layers. Although projections of some scenarios entail decreasing salinities near the surface, most of them project a rise in subsoil salinity, with the strongest trends of up to +0.9 mg cm(-3) 100 yr(-1) at -65 cm. The results suggest that topsoil salinity trends in the study area are affected by the magnitude of winter rainfall trends, whereas high subsoil salinities correspond to low winter rainfall and high summer temperature. How these projected trends affect the vegetation and thereby future land use will depend on the future management of groundwater levels in the area.
Inverted perovskite solar cells still suffer from significant non-radiative recombination losses at the perovskite surface and across the perovskite/C-60 interface, limiting the future development of perovskite-based single- and multi-junction photovoltaics. Therefore, more effective inter- or transport layers are urgently required. To tackle these recombination losses, we introduce ortho-carborane as an interlayer material that has a spherical molecular structure and a three-dimensional aromaticity. Based on a variety of experimental techniques, we show that ortho-carborane decorated with phenylamino groups effectively passivates the perovskite surface and essentially eliminates the non-radiative recombination loss across the perovskite/C-60 interface with high thermal stability. We further demonstrate the potential of carborane as an electron transport material, facilitating electron extraction while blocking holes from the interface. The resulting inverted perovskite solar cells deliver a power conversion efficiency of over 23% with a low non-radiative voltage loss of 110mV, and retain >97% of the initial efficiency after 400h of maximum power point tracking. Overall, the designed carborane based interlayer simultaneously enables passivation, electron-transport and hole-blocking and paves the way toward more efficient and stable perovskite solar cells. Effective transport layers are essential to suppress non-radiative recombination losses. Here, the authors introduce phenylamino-functionalized ortho-carborane as an interfacial layer, and realise inverted perovskite solar cells with efficiency of over 23% and operational stability of T97=400h.
We study populations of globally coupled noisy rotators (oscillators with inertia) allowing a nonequilibrium transition from a desynchronized state to a synchronous one (with the nonvanishing order parameter). The newly developed analytical approaches resulted in solutions describing the synchronous state with constant order parameter for weakly inertial rotators, including the case of zero inertia, when the model is reduced to the Kuramoto model of coupled noise oscillators. These approaches provide also analytical criteria distinguishing supercritical and subcritical transitions to the desynchronized state and indicate the universality of such transitions in rotator ensembles. All the obtained analytical results are confirmed by the numerical ones, both by direct simulations of the large ensembles and by solution of the associated Fokker-Planck equation. We also propose generalizations of the developed approaches for setups where different rotators parameters (natural frequencies, masses, noise intensities, strengths and phase shifts in coupling) are dispersed.
The motility of adherent eukaryotic cells is driven by the dynamics of the actin cytoskeleton. Despite the common force-generating actin machinery, different cell types often show diverse modes of locomotion that differ in their shape dynamics, speed, and persistence of motion. Recently, experiments in Dictyostelium discoideum have revealed that different motility modes can be induced in this model organism, depending on genetic modifications, developmental conditions, and synthetic changes of intracellular signaling. Here, we report experimental evidence that in a mutated D. discoideum cell line with increased Ras activity, switches between two distinct migratory modes, the amoeboid and fan-shaped type of locomotion, can even spontaneously occur within the same cell. We observed and characterized repeated and reversible switchings between the two modes of locomotion, suggesting that they are distinct behavioral traits that coexist within the same cell. We adapted an established phenomenological motility model that combines a reaction-diffusion system for the intracellular dynamics with a dynamic phase field to account for our experimental findings.
I study deterministic dynamics of chiral active particles in two dimensions. Particles are considered as discs interacting with elastic repulsive forces. An ensemble of particles, started from random initial conditions, demonstrates chaotic collisions resulting in their normal diffusion. This chaos is transient, as rather abruptly a synchronous collisionless state establishes. The life time of chaos grows exponentially with the number of particles. External forcing (periodic or chaotic) is shown to facilitate the synchronization transition.
'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.
According to established understanding, deep-water formation in the North Atlantic and Southern Ocean keeps the deep ocean cold, counter-acting the downward mixing of heat from the warmer surface waters in the bulk of the world ocean. Therefore, periods of strong Atlantic meridional overturning circulation (AMOC) are expected to coincide with cooling of the deep ocean and warming of the surface waters. It has recently been proposed that this relation may have reversed due to global warming, and that during the past decades a strong AMOC coincides with warming of the deep ocean and relative cooling of the surface, by transporting increasingly warmer waters downward. Here we present multiple lines of evidence, including a statistical evaluation of the observed global mean temperature, ocean heat content, and different AMOC proxies, that lead to the opposite conclusion: even during the current ongoing global temperature rise a strong AMOC warms the surface. The observed weakening of the AMOC has therefore delayed global surface warming rather than enhancing it.
Social Media Abstract:
The overturning circulation in the Atlantic Ocean has weakened in response to global warming, as predicted by climate models. Since it plays an important role in transporting heat, nutrients and carbon, a slowdown will affect global climate processes and the global mean temperature. Scientists have questioned whether this slowdown has worked to cool or warm global surface temperatures. This study analyses the overturning strength and global mean temperature evolution of the past decades and shows that a slowdown acts to reduce the global mean temperature. This is because a slower overturning means less water sinks into the deep ocean in the subpolar North Atlantic. As the surface waters are cold there, the sinking normally cools the deep ocean and thereby indirectly warms the surface, thus less sinking implies less surface warming and has a cooling effect. For the foreseeable future, this means that the slowing of the overturning will likely continue to slightly reduce the effect of the general warming due to increasing greenhouse gas concentrations.
Wire Arc Additive Manufacturing (WAAM) features high deposition rates and, thus, allows production of large components that are relevant for aerospace applications. However, a lot of aerospace parts are currently produced by forging or machining alone to ensure fast production and to obtain good mechanical properties; the use of these conventional process routes causes high tooling and material costs. A hybrid approach (a combination of forging and WAAM) allows making production more efficient. In this fashion, further structural or functional features can be built in any direction without using additional tools for every part. By using a combination of forging basic geometries with one tool set and adding the functional features by means of WAAM, the tool costs and material waste can be reduced compared to either completely forged or machined parts. One of the factors influencing the structural integrity of additively manufactured parts are (high) residual stresses, generated during the build process. In this study, the triaxial residual stress profiles in a hybrid WAAM part are reported, as determined by neutron diffraction. The analysis is complemented by microstructural investigations, showing a gradient of microstructure (shape and size of grains) along the part height. The highest residual stresses were found in the transition zone (between WAAM and forged part). The total stress range showed to be lower than expected for WAAM components. This could be explained by the thermal history of the component.
Organic solar cells have the potential to become the cheapest form of electricity. Rapid increase in the power conversion efficiency of organic solar cells (OSCs) has been achieved with the development of non-fullerene small-molecule acceptors. Next generation photovoltaics based upon environmentally benign "green solvent" processing of organic semiconductors promise a step-change in the adaptability and versatility of solar technologies and promote sustainable development. However, high-performing OSCs are still processed by halogenated (non-environmentally friendly) solvents, so hindering their large-scale manufacture. In this perspective, we discuss the recent progress in developing highly efficient OSCs processed from eco-compatible solvents, and highlight research challenges that should be addressed for the future development of high power conversion efficiencies devices.
We 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.
Non-geminate recombination, as one of the most relevant loss mechanisms in organic and perovskite solar cells, deserves special attention in research efforts to further increase device performance. It can be subdivided into first, second, and third order processes, which can be elucidated by the effects that they have on the time-dependent open-circuit voltage decay. In this study, analytical expressions for the open-circuit voltage decay exhibiting one of the aforementioned recombination mechanisms were derived. It was possible to support the analytical models with experimental examples of three different solar cells, each of them dominated either by first (PBDBT:CETIC-4F), second (PM6:Y6), or third (irradiated CH3NH3PbI3) order recombination. Furthermore, a simple approach to estimate the dominant recombination process was also introduced and tested on these examples. Moreover, limitations of the analytical models and the measurement technique itself were discussed.
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.
The reconstruction of cone-beam computed tomography data using filtered back-projection algorithms unavoidably results in severe artefacts. We describe how the Direct Iterative Reconstruction of Computed Tomography Trajectories (DIRECTT) algorithm can be combined with a model of the artefacts for the reconstruction of such data. The implementation of DIRECTT results in reconstructed volumes of superior quality compared to the conventional algorithms.
Three-component molecular brushes with a polyimide backbone and amphiphilic block copolymer side chains with different contents of the "inner" hydrophilic (poly(methacrylic acid)) and "outer" hydrophobic (poly(methyl methacrylate)) blocks were synthesized and characterized by molecular hydrodynamics and optics methods in solutions of chloroform, dimethylformamide, tetrahydrofuran and ethanol. The peculiarity of the studied polymers is the amphiphilic structure of the grafted chains. The molar masses of the molecular brushes were determined by static and dynamic light scattering in chloroform in which polymers form molecularly disperse solutions. Spontaneous self-assembly of macromolecules was detected in dimethylformamide, tetrahydrofuran and ethanol. The aggregates size depended on the thermodynamic quality of the solvent as well as on the macromolecular architectural parameters. In dimethylformamide and tetrahydrofuran, the distribution of hydrodynamic radii of aggregates was bimodal, while in ethanol, it was unimodal. Moreover, in ethanol, an increase in the poly(methyl methacrylate) content caused a decrease in the hydrodynamic radius of aggregates. A significant difference in the nature of the blocks included in the brushes determines the selectivity of the used solvents, since their thermodynamic quality with respect to the blocks is different. The macromolecules of the studied graft copolymers tend to self-organization in selective solvents with formation of a core-shell structure with an insoluble solvophobic core surrounded by the solvophilic shell of side chains.
Spin precession in magnetic materials is commonly modelled with the classical phenomenological Landau-Lifshitz-Gilbert (LLG) equation. Based on a quantized three-dimensional spin + environment Hamiltonian, we here derive a spin operator equation of motion that describes precession and includes a general form of damping that consistently accounts for memory, coloured noise and quantum statistics. The LLG equation is recovered as its classical, Ohmic approximation. We further introduce resonant Lorentzian system-reservoir couplings that allow a systematic comparison of dynamics between Ohmic and non-Ohmic regimes. Finally, we simulate the full non-Markovian dynamics of a spin in the semi-classical limit. At low temperatures, our numerical results demonstrate a characteristic reduction and flattening of the steady state spin alignment with an external field, caused by the quantum statistics of the environment. The results provide a powerful framework to explore general three-dimensional dissipation in quantum thermodynamics.
We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes-namely, the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving while varying the Alfven Mach number (M-A), plasma beta, and the sonic Mach number from subsonic to transsonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven-mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high M-A. The slow-mode properties are similar to the Alfven mode. On the other hand, the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast-mode behavior does not show any transition in going from low to high M-A. We find indications that there is some interaction between the different modes, and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic-ray transport and magnetic reconnection are discussed.
In this work, which is part of a larger research program, a framework called "virtual data fusion" was developed to provide an automated and consistent crack detection method that allows for the cross-comparison of results from large quantities of X-ray computed tomography (CT) data. A partial implementation of this method in a custom program was developed for use in research focused on crack quantification in alkali-silica reaction (ASR)-sensitive concrete aggregates. During the CT image processing, a series of image analyses tailored for detecting specific, individual crack-like characteristics were completed. The results of these analyses were then "fused" in order to identify crack-like objects within the images with much higher accuracy than that yielded by any individual image analysis procedure. The results of this strategy demonstrated the success of the program in effectively identifying crack-like structures and quantifying characteristics, such as surface area and volume. The results demonstrated that the source of aggregate has a very significant impact on the amount of internal cracking, even when the mineralogical characteristics remain very similar. River gravels, for instance, were found to contain significantly higher levels of internal cracking than quarried stone aggregates of the same mineralogical type.
In recent years, the search for more efficient and environmentally friendly materials to be employed in the next generation of thin film solar cell devices has seen a shift towards hybrid halide perovskites and chalcogenide materials crystallising in the kesterite crystal structure. Prime examples for the latter are Cu2ZnSnS4, Cu2ZnSnSe4, and their solid solution Cu2ZnSn(SxSe1-x)(4), where actual devices already demonstrated power conversion efficiencies of about 13 %. However, in their naturally occurring kesterite crystal structure, the so-called Cu-Zn disorder plays an important role and impacts the structural, electronic, and optical properties. To understand the influence of Cu-Zn disorder, we perform first-principles calculations based on density functional theory combined with special quasirandom structures to accurately model the cation disorder. Since the electronic band gaps and derived optical properties are severely underestimated by (semi)local exchange and correlation functionals, supplementary hybrid functional calculations have been performed. Concerning the latter, we additionally employ a recently devised technique to speed up structural relaxations for hybrid functional calculations. Our calculations show that the Cu-Zn disorder leads to a slight increase in the unit cell volume compared to the conventional kesterite structure showing full cation order, and that the band gap gets reduced by about 0.2 eV, which is in very good agreement with earlier experimental and theoretical findings. Our detailed results on structural, electronic, and optical properties will be discussed with respect to available experimental data, and will provide further insights into the atomistic origin of the disorder-induced band gap lowering in these promising kesterite type materials.
The mixture of ammonium nitrate (AN) prills and fuel oil (FO), usually referred to as ANFO, is extensively used in the mining industry as a bulk explosive. One of the major performance predictors of ANFO mixtures is the fuel oil retention, which is itself governed by the complex pore structure of the AN prills. In this study, we present how X-ray computed tomography (XCT), and the associated advanced data processing workflow, can be used to fully characterise the structure and morphology of AN prills. We show that structural parameters such as volume fraction of the different phases and morphological parameters such as specific surface area and shape factor can be reliably extracted from the XCT data, and that there is a good agreement with the measured oil retention values. Importantly, oil retention measurements (qualifying the efficiency of ANFO as explosives) correlate well with the specific surface area determined by XCT. XCT can therefore be employed non-destructively; it can accurately evaluate and characterise porosity in ammonium nitrate prills, and even predict their efficiency.
The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120-600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.
At supranuclear densities, explored in the core of neutron stars, a strong phase transition from hadronic matter to more exotic forms of matter might be present. To test this hypothesis, binary neutron-star mergers offer a unique possibility to probe matter at densities that we cannot create in any existing terrestrial experiment. In this work, we show that, if present, strong phase transitions can have a measurable imprint on the binary neutron-star coalescence and the emitted gravitational-wave signal. We construct a new parametrization of the supranuclear equation of state that allows us to test for the existence of a strong phase transition and extract its characteristic properties purely from the gravitational-wave signal of the inspiraling neutron stars. We test our approach using a Bayesian inference study simulating 600 signals with three different equations of state and find that for current gravitational-wave detector networks already 12 events might be sufficient to verify the presence of a strong phase transition. Finally, we use our methodology to analyze GW170817 and GW190425 but do not find any indication that a strong phase transition is present at densities probed during the inspiral.
Electric currents flowing in the terrestrial ionosphere have conventionally been diagnosed by low-earth-orbit (LEO) satellites equipped with science-grade magnetometers and long booms on magnetically clean satellites. In recent years, there are a variety of endeavors to incorporate platform magnetometers, which are initially designed for navigation purposes, to study ionospheric currents. Because of the suboptimal resolution and significant noise of the platform magnetometers, however, most of the studies were confined to high-latitude auroral regions, where magnetic field deflections from ionospheric currents easily exceed 100 nT. This study aims to demonstrate the possibility of diagnosing weak low-/mid-latitude ionospheric currents based on platform magnetometers. We use navigation magnetometer data from two satellites, CryoSat-2 and the Gravity Recovery and Climate Experiment Follow-On (GRACE-FO), both of which have been intensively calibrated based on housekeeping data and a high-precision geomagnetic field model. Analyses based on 8 years of CryoSat-2 data as well as similar to 1.5 years of GRACE-FO data reproduce well-known climatology of inter-hemispheric field-aligned currents (IHFACs), as reported by previous satellite missions dedicated to precise magnetic observations. Also, our results show that C-shaped structures appearing in noontime IHFAC distributions conform to the shape of the South Atlantic Anomaly. The F-region dynamo currents are only partially identified in the platform magnetometer data, possibly because the currents are weaker than IHFACs in general and depend significantly on altitude and solar activity. Still, this study evidences noontime F-region dynamo currents at the highest altitude (717 km) ever reported. We expect that further data accumulation from continuously operating missions may reveal the dynamo currents more clearly during the next solar maximum.
Strong events of long-range transported biomass burning aerosol were detected during July 2013 at three EARLINET (European Aerosol Research Lidar Network) stations, namely Granada (Spain), Leipzig (Germany) and Warsaw (Poland). Satellite observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) instruments, as well as modeling tools such as HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) and NAAPS (Navy Aerosol Analysis and Prediction System), have been used to estimate the sources and transport paths of those North American forest fire smoke particles. A multiwavelength Raman lidar technique was applied to obtain vertically resolved particle optical properties, and further inversion of those properties with a regularization algorithm allowed for retrieving microphysical information on the studied particles. The results highlight the presence of smoke layers of 1-2 km thickness, located at about 5 km a.s.l. altitude over Granada and Leipzig and around 2.5 km a.s.l. at Warsaw. These layers were intense, as they accounted for more than 30% of the total AOD (aerosol optical depth) in all cases, and presented optical and microphysical features typical for different aging degrees: color ratio of lidar ratios (LR532/LR355) around 2, alpha-related angstrom exponents of less than 1, effective radii of 0.3 mu m and large values of single scattering albedos (SSA), nearly spectrally independent. The intensive microphysical properties were compared with columnar retrievals form co-located AERONET (Aerosol Robotic Network) stations. The intensity of the layers was also characterized in terms of particle volume concentration, and then an experimental relationship between this magnitude and the particle extinction coefficient was established.
Earth’s surface temperature will continue to rise for another 20 to 30 years even with the strongest carbon emission reduction currently considered. The associated changes in rainfall patterns can result in an increased flood risk worldwide. We compute the required increase in flood protection to keep high-end fluvial flood risk at present levels. The analysis is carried out worldwide for subnational administrative units. Most of the United States, Central Europe, and Northeast and West Africa, as well as large parts of India and Indonesia, require the strongest adaptation effort. More than half of the United States needs to at least double their protection within the next two decades. Thus, the need for adaptation to increased river flood is a global problem affecting industrialized regions as much as developing countries.
We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes-namely, the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving while varying the Alfven Mach number (M-A), plasma beta, and the sonic Mach number from subsonic to transsonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven-mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high M-A. The slow-mode properties are similar to the Alfven mode. On the other hand, the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast-mode behavior does not show any transition in going from low to high M-A. We find indications that there is some interaction between the different modes, and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic-ray transport and magnetic reconnection are discussed.
The stratospheric polar vortex can influence the tropospheric circulation and thereby winter weather in the mid-latitudes. Weak vortex states, often associated with sudden stratospheric warmings (SSW), have been shown to increase the risk of cold-spells especially over Eurasia, but its role for North American winters is less clear. Using cluster analysis, we show that there are two dominant patterns of increased polar cap heights in the lower stratosphere. Both patterns represent a weak polar vortex but they are associated with different wave mechanisms and different regional tropospheric impacts. The first pattern is zonally symmetric and associated with absorbed upward-propagating wave activity, leading to a negative phase of the North Atlantic Oscillation (NAO) and cold-air outbreaks over northern Eurasia. This coupling mechanism is well-documented in the literature and is consistent with the downward migration of the northern annular mode (NAM). The second pattern is zonally asymmetric and linked to downward reflected planetary waves over Canada followed by a negative phase of the Western Pacific Oscillation (WPO) and cold-spells in Central Canada and the Great Lakes region. Causal effect network (CEN) analyses confirm the atmospheric pathways associated with this asymmetric pattern. Moreover, our findings suggest the reflective mechanism to be sensitive to the exact region of upward wave-activity fluxes and to be state-dependent on the strength of the vortex. Identifying the causal pathways that operate on weekly to monthly timescales can pave the way for improved sub-seasonal to seasonal forecasting of cold spells in the mid-latitudes.
The power spectral density (PSD) of any time-dependent stochastic processX (t) is ameaningful feature of its spectral content. In its text-book definition, the PSD is the Fourier transform of the covariance function of X-t over an infinitely large observation timeT, that is, it is defined as an ensemble-averaged property taken in the limitT -> infinity. Alegitimate question is what information on the PSD can be reliably obtained from single-trajectory experiments, if one goes beyond the standard definition and analyzes the PSD of a single trajectory recorded for a finite observation timeT. In quest for this answer, for a d-dimensional Brownian motion (BM) we calculate the probability density function of a single-trajectory PSD for arbitrary frequency f, finite observation time T and arbitrary number k of projections of the trajectory on different axes. We show analytically that the scaling exponent for the frequency-dependence of the PSD specific to an ensemble of BM trajectories can be already obtained from a single trajectory, while the numerical amplitude in the relation between the ensemble-averaged and single-trajectory PSDs is afluctuating property which varies from realization to realization. The distribution of this amplitude is calculated exactly and is discussed in detail. Our results are confirmed by numerical simulations and single-particle tracking experiments, with remarkably good agreement. In addition we consider a truncated Wiener representation of BM, and the case of a discrete-time lattice random walk. We highlight some differences in the behavior of a single-trajectory PSD for BM and for the two latter situations. The framework developed herein will allow for meaningful physical analysis of experimental stochastic trajectories.
Here, we present novel equations for the large-scale zonal-mean wind as well as those for planetary waves. Together with synoptic parameterization (as presented by Coumou et al., 2011), these form the mathematical description of the dynamical core of Aeolus 1.0. The regions of high azonal wind velocities (planetary waves) are accurately captured for all validation experiments. The zonal-mean zonal wind and the integrated lower troposphere mass flux show good results in particular in the Northern Hemisphere. In the Southern Hemisphere, the model tends to produce too-weak zonal-mean zonal winds and a too-narrow Hadley circulation. We discuss possible reasons for these model biases as well as planned future model improvements and applications.
Magnetic nanoparticles are met across many biological species ranging from magnetosensitive bacteria, fishes, bees, bats, rats, birds, to humans. They can be both of biogenetic origin and due to environmental contamination, being either in paramagnetic or ferromagnetic state. The energy of such naturally occurring single-domain magnetic nanoparticles can reach up to 10-20 room k(B)T in the magnetic field of the Earth, which naturally led to supposition that they can serve as sensory elements in various animals. This work explores within a stochastic modeling framework a fascinating hypothesis of magnetosensitive ion channels with magnetic nanoparticles serving as sensory elements, especially, how realistic it is given a highly dissipative viscoelastic interior of living cells and typical sizes of nanoparticles possibly involved.
In two-dimensional reaction-diffusion systems, local curvature perturbations on traveling waves are typically damped out and vanish. However, if the inhibitor diffuses much faster than the activator, transversal instabilities can arise, leading from flat to folded, spatio-temporally modulated waves and to spreading spiral turbulence. Here, we propose a scheme to induce or inhibit these instabilities via a spatio-temporal feedback loop. In a piecewise-linear version of the FitzHugh-Nagumo model, transversal instabilities and spiral turbulence in the uncontrolled system are shown to be suppressed in the presence of control, thereby stabilizing plane wave propagation. Conversely, in numerical simulations with the modified Oregonator model for the photosensitive Belousov-Zhabotinsky reaction, which does not exhibit transversal instabilities on its own, we demonstrate the feasibility of inducing transversal instabilities and study the emerging wave patterns in a well-controlled manner.
The fact that organic solar cells perform efficiently despite the low dielectric constant of most photoactive blends initiated a long-standing debate regarding the dominant pathways of free charge formation. Here, we address this issue through the accurate measurement of the activation energy for free charge photogeneration over a wide range of photon energy, using the method of time-delayed collection field. For our prototypical low bandgap polymer:fullerene blends, we find that neither the temperature nor the field dependence of free charge generation depend on the excitation energy, ruling out an appreciable contribution to free charge generation though hot carrier pathways. On the other hand, activation energies are on the order of the room temperature thermal energy for all studied blends. We conclude that charge generation in such devices proceeds through thermalized charge transfer states, and that thermal energy is sufficient to separate most of these states into free charges.
Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of sub-shelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice-ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. We identify a set of parameters that yield two-dimensional melt rate fields that qualitatively reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 ma(-1) for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 ma(-1) for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.
Understanding of wave environments is critical for the understanding of how particles are accelerated and lost in space. This study shows that in the vicinity of Europa and Ganymede, that respectively have induced and internal magnetic fields, chorus wave power is significantly increased. The observed enhancements are persistent and exceed median values of wave activity by up to 6 orders of magnitude for Ganymede. Produced waves may have a pronounced effect on the acceleration and loss of particles in the Jovian magnetosphere and other astrophysical objects. The generated waves are capable of significantly modifying the energetic particle environment, accelerating particles to very high energies, or producing depletions in phase space density. Observations of Jupiter’s magnetosphere provide a unique opportunity to observe how objects with an internal magnetic field can interact with particles trapped in magnetic fields of larger scale objects.