Refine
Has Fulltext
- no (283) (remove)
Year of publication
- 2019 (283) (remove)
Document Type
- Article (256)
- Other (19)
- Doctoral Thesis (4)
- Review (3)
- Monograph/Edited Volume (1)
Is part of the Bibliography
- yes (283)
Keywords
- stars: evolution (7)
- stars: winds, outflows (7)
- diffusion (6)
- methods: numerical (6)
- quasars: absorption lines (6)
- stars: massive (6)
- intergalactic medium (5)
- stars: mass-loss (5)
- subdwarfs (5)
- ISM: supernova remnants (4)
Institute
- Institut für Physik und Astronomie (283) (remove)
One of the biggest successes of the Cassini mission is the detection of small moons (moonlets) embedded in Saturns rings that cause S-shaped density structures in their close vicinity, called propellers. Here, we present isothermal hydrodynamic simulations of moonlet-induced propellers in Saturn's A ring that denote a further development of the original model. We find excellent agreement between these new hydrodynamic and corresponding N-body simulations. Furthermore, the hydrodynamic simulations confirm the predicted scaling laws and the analytical solution for the density in the propeller gaps. Finally, this mean field approach allows us to simulate the pattern of the giant propeller Blériot, which is too large to be modeled by direct N-body simulations. Our results are compared to two stellar occultation observations by the Cassini Ultraviolet Imaging Spectrometer (UVIS), which intersect the propeller Blériot. Best fits to the UVIS optical depth profiles are achieved for a Hill radius of 590 m, which implies a moonlet diameter of about 860 m. Furthermore, the model favors a kinematic shear viscosity of the surrounding ring material of ν0 = 340 cm2 s−1, a dispersion velocity in the range of 0.3 cm s−1 < c0 < 1.5 cm s−1, and a fairly high bulk viscosity 7 < ξ0/ν0 < 17. These large transport values might be overestimated by our isothermal ring model and should be reviewed by an extended model including thermal fluctuations.
The oceans and atmosphere interact via a multiplicity of feedback mechanisms, shaping to a large extent the global climate and its variability. To deepen our knowledge of the global climate system, characterizing and investigating this interdependence is an important task of contemporary research. However, our present understanding of the underlying large-scale processes is greatly limited due to the manifold interactions between essential climatic variables at different temporal scales. To address this problem, we here propose to extend the application of complex network techniques to capture the interdependence between global fields of sea-surface temperature (SST) and precipitation (P) at multiple temporal scales. For this purpose, we combine time-scale decomposition by means of a discrete wavelet transform with the concept of coupled climate network analysis. Our results demonstrate the potential of the proposed approach to unravel the scale-specific interdependences between atmosphere and ocean and, thus, shed light on the emerging multiscale processes inherent to the climate system, which traditionally remain undiscovered when investigating the system only at the native resolution of existing climate data sets. Moreover, we show how the relevant spatial interdependence structures between SST and P evolve across time-scales. Most notably, the strongest mutual correlations between SST and P at annual scale (8-16 months) concentrate mainly over the Pacific Ocean, while the corresponding spatial patterns progressively disappear when moving toward longer time-scales. Published under license by AIP Publishing.
The conversion of light energy into other molecular energetic degrees of freedom is often dominated by ultrafast, non-adiabatic processes. Femtosecond spectroscopy with optical pulses has helped in shaping our understanding of crucial processes in molecular energy-conversion. The advent of new, ultrashort and bright X-ray free electron laser sources opens the possibility to use X-ray-typical element and site sensitivity for ultrafast molecular research. We present two types of spectroscopy, ultrafast Auger and ultrafast X-ray absorption spectroscopy, and discuss their sensitivity to molecular processes. While Auger spectroscopy is able to monitor bond distance changes in the vicinity of an X-ray created core hole, near-edge absorption spectroscopy can deliver high-fidelity information on non-adiabatic transitions involving lone-pair orbitals. We demonstrate these features on the example of the UV-excited nucleobase thymine, investigated at the oxygen K-edge. We find a C-O bond elongation in the Auger data in addition to pi pi*/n pi* non-adiabatic transition in X-ray near-edge absorption. We compare the results from both methods and draw a conclusive scenario of non-adiabatic molecular relaxation after UV excitation.
The Chromospheric Telescope (ChroTel) is a small 10-cm robotic telescope at Observatorio del Teide on Tenerife (Spain), which observes the entire sun in Hα, Ca ii K, and He i 10 830 Å. We present a new calibration method that includes limb-darkening correction, removal of nonuniform filter transmission, and determination of He i Doppler velocities. Chromospheric full-disk filtergrams are often obtained with Lyot filters, which may display nonuniform transmission causing large-scale intensity variations across the solar disk. Removal of a 2D symmetric limb-darkening function from full-disk images results in a flat background. However, transmission artifacts remain and are even more distinct in these contrast-enhanced images. Zernike polynomials are uniquely appropriate to fit these large-scale intensity variations of the background. The Zernike coefficients show a distinct temporal evolution for ChroTel data, which is likely related to the telescope's alt-azimuth mount that introduces image rotation. In addition, applying this calibration to sets of seven filtergrams that cover the He i triplet facilitates the determination of chromospheric Doppler velocities. To validate the method, we use three datasets with varying levels of solar activity. The Doppler velocities are benchmarked with respect to cotemporal high-resolution spectroscopic data of the GREGOR Infrared Spectrograph (GRIS). Furthermore, this technique can be applied to ChroTel Hα and Ca ii K data. The calibration method for ChroTel filtergrams can be easily adapted to other full-disk data exhibiting unwanted large-scale variations. The spectral region of the He i triplet is a primary choice for high-resolution near-infrared spectropolarimetry. Here, the improved calibration of ChroTel data will provide valuable context data.
Electromagnetic ion cyclotron waves have long been recognized to play a crucial role in the dynamic loss of ring current protons. While the field-aligned propagation approximation of electromagnetic ion cyclotron waves was widely used to quantify the scattering loss of ring current protons, in this study, we find that the wave normal distribution strongly affects the pitch angle scattering efficiency of protons. Increase of peak normal angle or angular width can considerably reduce the scattering rates of <= 10 keV protons. For >10 keV protons, the field-aligned propagation approximation results in a pronounced underestimate of the scattering of intermediate equatorial pitch angle protons and overestimates the scattering of high equatorial pitch angle protons by orders of magnitude. Our results suggest that the wave normal distribution of electromagnetic ion cyclotron waves plays an important role in the pitch angle evolution and scattering loss of ring current protons and should be incorporated in future global modeling of ring current dynamics.
For the layered transition metal dichalcogenide 1T-TaS2, we establish through a unique experimental approach and density functional theory, how ultrafast charge transfer in 1T-TaS2 takes on isotropic three-dimensional character or anisotropic two-dimensional character, depending on the commensurability of the charge density wave phases of 1T-TaS2. The X-ray spectroscopic core-hole-clock method prepares selectively in-and out-of-plane polarized sulfur 3p orbital occupation with respect to the 1T-TaS2 planes and monitors sub-femtosecond wave packet delocalization. Despite being a prototypical two-dimensional material, isotropic three-dimensional charge transfer is found in the commensurate charge density wave phase (CCDW), indicating strong coupling between layers. In contrast, anisotropic two-dimensional charge transfer occurs for the nearly commensurate phase (NCDW). In direct comparison, theory shows that interlayer interaction in the CCDW phase - not layer stacking variations - causes isotropic three-dimensional charge transfer. This is presumably a general mechanism for phase transitions and tailored properties of dichalcogenides with charge density waves.
Based on data from the ESA Gaia Data Release 2 (DR2) and several ground-based, multi-band photometry surveys we have compiled an all-sky catalogue of 39 800 hot subluminous star candidates selected in Gaia DR2 by means of colour, absolute magnitude, and reduced proper motion cuts. We expect the majority of the candidates to be hot subdwarf stars of spectral type B and O, followed by blue horizontal branch stars of late B-type (HBB), hot post-AGB stars, and central stars of planetary nebulae. The contamination by cooler stars should be about 10%. The catalogue is magnitude limited to Gaia G < 19 mag and covers the whole sky. Except within the Galactic plane and LMC/SMC regions, we expect the catalogue to be almost complete up to about 1.5 kpc. The main purpose of this catalogue is to serve as input target list for the large-scale photometric and spectroscopic surveys which are ongoing or scheduled to start in the coming years. In the long run, securing a statistically significant sample of spectroscopically confirmed hot subluminous stars is key to advance towards a more detailed understanding of the latest stages of stellar evolution for single and binary stars.
We study the requirement of the jet power in the conventional p-gamma models (photopion production and Bethe-Heitler pair production) for TeV BL Lac objects. We select a sample of TeV BL Lac objects whose spectral energy distributions are difficult to explain by the one-zone leptonic model. Based on the relation between the p-gamma interaction efficiency and the opacity of gamma gamma absorption, we find that the detection of TeV emission poses upper limits on the p-gamma interaction efficiencies in these sources and hence minimum jet powers can be derived accordingly. We find that the obtained minimum jet powers exceed the Eddington luminosity of the supermassive black holes (SMBHs). Implications for the accretion mode of the SMBHs in these BL Lac objects and the origin of their TeV emissions are discussed.
Various studies have implied the existence of a gaseous halo around the Galaxy extending out to similar to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the CR distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, and CR spectral index. We find that the current measurement of isotropic gamma-ray background (IGRB) at less than or similar to TeV with the Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic CR (e.g., with CR luminosity L-CR <= 1041 erg s(-1)). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between CRs and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restrictions on the relevant quantities (such as Galactic CR luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.
Theoretical Fluid Dynamics
(2019)
We present a survey for metal absorption systems traced by neutral oxygen over 3.2 < z < 6.5. Our survey uses Keck/ESI and VLT/X-Shooter spectra of 199 QSOs with redshifts up to 6.6. In total, we detect 74 OI absorbers, of which 57 are separated from the background QSO by more than 5000 km s(-1). We use a maximum likelihood approach to fit the distribution of OI lambda 1302 equivalent widths in bins of redshift and from this determine the evolution in number density of absorbers with W-1302 > 0.05 angstrom, of which there are 49 nonproximate systems in our sample. We find that the number density does not monotonically increase with decreasing redshift, as would naively be expected from the buildup of metal-enriched circumgalactic gas with time. The number density over 4.9 < z < 5.7 is a factor of 1.7-4.1 lower (68% confidence) than that over 5.7 < z < 6.5, with a lower value at z < 5.7 favored with 99% confidence. This decrease suggests that the fraction of metals in a low-ionization phase is larger at z similar to 6 than at lower redshifts. Absorption from highly ionized metals traced by CIV is also weaker in higher-redshift OI systems, supporting this picture. The evolution of OI absorbers implies that metal-enriched circumgalactic gas at z similar to 6 is undergoing an ionization transition driven by a strengthening ultraviolet background. This in turn suggests that the reionization of the diffuse intergalactic medium may still be ongoing at or only recently ended by this epoch.
We have investigated the structural dynamics in photoexcited 1,2-diiodotetrafluoroethane molecules (C2F4I2) in the gas phase experimentally using ultrafast electron diffraction and theoretically using FOMO-CASCI excited-state dynamics simulations. The molecules are excited by an ultraviolet femtosecond laser pulse to a state characterized by a transition from the iodine 5p perpendicular to orbital to a mixed 5p parallel to sigma hole and CF2 center dot antibonding orbital, which results in the cleavage of one of the carbon-iodine bonds. We have observed, with sub-Angstrom resolution, the motion of the nuclear wave packet of the dissociating iodine atom followed by coherent vibrations in the electronic ground state of the C2F4I radical. The radical reaches a stable classical (nonbridged) structure in less than 200 fs.
Current algorithms for the real-time prediction of the Kp index use a combination of models empirically driven by solar wind measurements at the L1 Lagrange point and historical values of the index. In this study, we explore the limitations of this approach, examining the forecast for short and long lead times using measurements at L1 and Kp time series as input to artificial neural networks. We explore the relative efficiency of the solar wind-based predictions, predictions based on recurrence, and predictions based on persistence. Our modeling results show that for short-term forecasts of approximately half a day, the addition of the historical values of Kp to the measured solar wind values provides a barely noticeable improvement. For a longer-term forecast of more than 2 days, predictions can be made using recurrence only, while solar wind measurements provide very little improvement for a forecast with long horizon times. We also examine predictions for disturbed and quiet geomagnetic activity conditions. Our results show that the paucity of historical measurements of the solar wind for high Kp results in a lower accuracy of predictions during disturbed conditions. Rebalancing of input data can help tailor the predictions for more disturbed conditions.
Strong light-matter coupling can re-arrange the exciton energies in organic semiconductors. Here, we exploit strong coupling by embedding a fullerene-free organic solar cell (OSC) photo-active layer into an optical microcavity, leading to the formation of polariton peaks and a red-shift of the optical gap. At the same time, the open-circuit voltage of the device remains unaffected. This leads to reduced photon energy losses for the low-energy polaritons and a steepening of the absorption edge. While strong coupling reduces the optical gap, the energy of the charge-transfer state is not affected for large driving force donor-acceptor systems. Interestingly, this implies that strong coupling can be exploited in OSCs to reduce the driving force for electron transfer, without chemical or microstructural modifications of the photoactive layer. Our work demonstrates that the processes determining voltage losses in OSCs can now be tuned, and reduced to unprecedented values, simply by manipulating the device architecture.
We studied the photoinduced ultrafast relaxation dynamics of the nucleobase thymine using gas-phase time-resolved photoelectron spectroscopy. By employing extreme ultraviolet pulses from high harmonic generation for photoionization, we substantially extend our spectral observation window with respect to previous studies. This enables us to follow relaxation of the excited state population all the way to low-lying electronic states including the ground state. In thymine, we observe relaxation from the optically bright (1)pi pi* state of thymine to a dark (1)n pi* state within 80 +/- 30 fs. The (1)n pi* state relaxes further within 3.5 +/- 0.3 ps to a low-lying electronic state. By comparison with quantum chemical simulations, we can unambiguously assign its spectroscopic signature to the (3)pi pi* state. Hence, our study draws a comprehensive picture of the relaxation mechanism of thymine including ultrafast intersystem crossing to the triplet manifold.
On 2015 July 18, near perihelion at a heliocentric distance of 1.28 au, the Visible InfraRed Thermal Imaging Spectrometer (VIRTIS-M) on board the Rosetta spacecraft had the opportunity of observing dust activity in the inner coma with a view of the night side (shadowed side) of comet 67P/Churyumov-Gerasimenko. At the time of the measurements we present here, we observe a dust plume that originates on the far side of the nucleus. We are able to identify the approximate location of its source at the boundary between the Hapi and Anuket regions, and we find that it has been in darkness for some hours before the observation. Assuming that this time span is equal to the conductive time scale, we obtain a thermal inertia in the range 25-36 W K-1 m(-2) s(-1/2). These thermal inertia values can be used to verify with a 3D finite-element method (REM) numerical code whether the surface and subsurface temperatures agree with the values found in the literature. We explored three different configurations: (1) a layer of water ice mixed with dust beneath a dust mantle of 5 mm with thermal inertia of 36 J m(-2) K-1 S-0.5 ; (2) the same structure, but with thermal inertia of 100 J m(-2) K-1 S-0.5; (3) an ice-dust mixture that is directly exposed. Of these three configurations, the first seems to be the most reasonable, both for the low thermal inertia and for the agreement with the surface and subsurface temperatures that have been found for the comet 67P/Churyumov-Gerasimenko. The spectral properties of the plume show that the visible dust color ranged from 16 +/- 4.8%/100 nm to 13 +/- 2.6%/100 nm, indicating that this plume has no detectable color gradient. The morphology of the plume can be classified as a narrow jet that has an estimated total ejected mass of between 6 and 19 tons when we assume size distribution indices between -2.5 and -3.
Context. On 27 April 2015, when comet 67P/Churyumov-Gerasimenko was at 1.76 au from the Sun and moving toward perihelion, the OSIRIS and VIRTIS-M instruments on board the Rosetta spacecraft simultaneously observed the evolving dust and gas coma during a complete rotation of the comet. Aims. We aim to characterize the spatial distribution of dust, H2O, and CO2 gas in the inner coma. To do this, we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the rate we calculated using a thermophysical model. Methods. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O emission band), and 4200 nm (CO2 emission band). We measured the average signal in a circular annulus to study the spatial variation around the comet, and in a sector of the annulus to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet center. Results. The spatial correlation between dust and water, both coming from the sunlit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0 degrees subsolar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production with respect to the value calculated using a simple thermophysical model is observed when the head lobe and regions of the southern hemisphere with strong seasonal variations are illuminated (subsolar longitude 270 degrees-50 degrees). A drastic decrease in dust production when the water production (both measured and from the model) displays a maximum occurs when typical northern consolidated regions are illuminated and the southern hemisphere regions with strong seasonal variations are instead in shadow (subsolar longitude 50 degrees-90 degrees). Possible explanations of these observations are presented and discussed.
We present an analysis of two pulsating subdwarf B stars PHL 457 and EQ Psc observed during the K2 mission. The K2 light curves of both stars show variation consistent with irradiation of a cooler companion by the hot subdwarf. They also show higher frequency oscillations consistent with pulsation. Using new spectroscopic data, we measured the radial velocity, effective temperature, surface gravity, and helium abundance of both hot subdwarfs as a function of orbital phase. We confirm the previously published spectroscopic orbit of PHL 457, and present the first spectroscopic orbit of EQ Psc. The orbital periods are 0.313 and 0.801 d, respectively. For EQPsc, we find a strong correlation between T-eff and orbital phase, due to contribution of light from the irradiated companion. We calculated amplitude spectra, identified significant pulsation frequencies, and searched for multiplets and asymptotic period spacings. By means of multiplets and period spacing, we identified the degrees of several pulsation modes in each star. The g-mode multiplets indicate subsynchronous core rotation with periods of 4.6 d (PHL 457) and 9.4 d (EQ Psc). We made spectral energy disctribution (SED) fits of PHL 457 and EQ Psc using available broad-band photometry and Gaia data. While the SED of PHL 457 shows no evidence of a cool companion, the SED for EQPsc clearly shows an infrared (IR) access consistent with a secondary with a temperature of about 6800K and a radius of 0.23 R-circle dot. This is the first detection of an IR access in any sdB + dM binary.
Diffusion of cosmic rays (CRs) is the key process for understanding their propagation and acceleration. We employ the description of spatial separation of magnetic field lines in magnetohydrodynamic turbulence in Lazarian & Vishniac to quantify the divergence of the magnetic field on scales less than the injection scale of turbulence and show that this divergence induces superdiffusion of CR in the direction perpendicular to the mean magnetic field. The perpendicular displacement squared increases, not as the distance x along the magnetic field, which is the case for a regular diffusion, but as the x 3 for freely streaming CRs. The dependence changes to x 3/2 for the CRs propagating diffusively along the magnetic field. In the latter case, we show that it is important to distinguish the perpendicular displacement with respect to the mean field and to the local magnetic field. We consider how superdiffusion changes the acceleration of CRs in shocks and show how it decreases efficiency of the CRs acceleration in perpendicular shocks. We also demonstrate that in the case when the small-scale magnetic field is generated in the pre-shock region, an efficient acceleration can take place for the CRs streaming without collisions along the magnetic loops.
“A chain is only as strong as its weakest link” says the proverb. But what about a collection of statistically identical chains: How long till all chains fail? The answer to this question is given by the max-min of a matrix whose (i,j)entry is the failure time of link j of chain i: take the minimum of each row, and then the maximum of the rows' minima. The corresponding min-max is obtained by taking the maximum of each column, and then the minimum of the columns' maxima. The min-max applies to the storage of critical data. Indeed, consider multiple backup copies of a set of critical data items, and consider the (i,j) matrix entry to be the time at which item j on copy i is lost; then, the min-max is the time at which the first critical data item is lost. In this paper we address random matrices whose entries are independent and identically distributed random variables. We establish Poisson-process limit laws for the row's minima and for the columns' maxima. Then, we further establish Gumbel limit laws for the max-min and for the min-max. The limit laws hold whenever the entries' distribution has a density, and yield highly applicable approximation tools and design tools for the max-min and min-max of large random matrices. A brief of the results presented herein is given in: Gumbel central limit theorem for max-min and min-max
The max-min and min-max of matrices arise prevalently in science and engineering. However, in many real-world situations the computation of the max-min and min-max is challenging as matrices are large and full information about their entries is lacking. Here we take a statistical-physics approach and establish limit laws—akin to the central limit theorem—for the max-min and min-max of large random matrices. The limit laws intertwine random-matrix theory and extreme-value theory, couple the matrix dimensions geometrically, and assert that Gumbel statistics emerge irrespective of the matrix entries' distribution. Due to their generality and universality, as well as their practicality, these results are expected to have a host of applications in the physical sciences and beyond.
In this paper we extract the aerosol microphysical properties for a collection of mineral dust cases measured by multi-wavelength depolarization Raman lidar systems located at the National Technical University of Athens (NTUA, Athens, Greece) and the Andalusian Institute for Earth System Research (IISTA-CEAMA, Granada, Spain). The lidar-based retrievals were carried out with the Spheroidal Inversion eXperiments software tool (SphInX) developed at the University of Potsdam (Germany). The software uses regularized inversion of a two-dimensional enhancement of the Mie model based on the spheroid-particle approximation with the aspect ratio determining the particle shape. The selection of the cases was based on the transport time from the source regions to the measuring sites. The aerosol optical depth as measured by AERONET ranged from 0.27 to 0.54 (at 500 nm) depending on the intensity of each event. Our analysis showed the hourly mean particle linear depolarization ratio and particle lidar ratio values at 532 nm ranging from 11 to 34% and from 42 to 79 sr respectively, depending on the mixing status, the corresponding air mass pathways and their transport time. Cases with shorter transport time showed good agreement in terms of the optical and SphInX-retrieved microphysical properties between Athens and Granada providing a complex refractive index value equal to 1.4 + 0.004i. On the other hand, the results for cases with higher transport time deviated from the aforementioned ones as well as from each other, providing, in particular, an imaginary part of the refractive index ranging from 0.002 to 0.005. Reconstructions of two-dimensional shape-size distributions for each selected layer showed that the dominant effective particle shape was prolate with diverse spherical contributions. The retrieved volume concentrations reflect overall the intensity of the episodes.
Melting experiments attempting to reproduce some of the processes affecting asteroidal and cometary material during atmospheric entry have been performed in a high enthalpy facility. For the first time with the specific experimental setup, the resulting material has been recovered, studied, and compared with natural analogues, focusing on the thermal and redox reactions triggered by interaction between the melt and the atmospheric gases under high temperature and low pressure conditions. Experimental conditions were tested across a range of parameters, such as heat flux, experiment duration, and pressure, using two types of sample holders materials, namely cork and graphite. A basalt served as asteroidal analog and to calibrate the experiments, before melting a H5 ordinary chondrite meteorite. The quenched melt recovered after the experiments has been analyzed by mu-XRF, EDS-SEM, EMPA, LA-ICP-MS, and XANES spectroscopy. The glass formed from the basalt is fairly homogeneous, depleted in highly volatile elements (e.g., Na, K), relatively enriched in moderately siderophile elements (e.g., Co, Ni), and has reached an equilibrium redox state with a lower Fe3+/Fe-tot ratio than that in the starting material. Spherical objects, enriched in SiO2, Na2O and K2O, were observed, inferring condensation from the vaporized material. Despite instantaneous quenching, the melt formed from the ordinary chondrite shows extensive crystallization of mostly olivine and magnetite, the latter indicative of oxygen fugacity compatible with presence of both Fe2+ and Fe3+. Similar features have been observed in natural meteorite fusion crusts and in micrometeorites, implying that, at least in terms of maximum temperature reached and chemical reactions, the experiments have successfully reproduced the conditions likely encountered by extraterrestrial material following atmospheric entry.
New wave frequency and amplitude models for the nightside and dayside chorus waves are built based on measurements from the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrument onboard the Van Allen Probes. The corresponding 3D diffusion coefficients are systematically obtained. Compared with previous commonly-used (typical) parameterizations, the new parameterizations result in differences in diffusion rates that depend on the energy and pitch angle. Furthermore, one-year 3D diffusive simulations are performed using the Versatile Electron Radiation Belt (VERB) code. Both typical and new wave parameterizations simulation results are in a good agreement with observations at 0.9 MeV. However, the new parameterizations for nightside chorus better reproduce the observed electron fluxes. These parameterizations will be incorporated into future modeling efforts.
Long-term river evolution depends partly on crustal deformation, which shapes the topography crossed by rivers. On glacial timescales, ice-sheet growth and decay can produce crustal vertical motion of ∼10 mm/yr resulting from the solid Earth's adjustment to variations in ice and water loads, comparable to tectonically-driven rates in the most rapidly uplifting mountains on Earth. This process of glacial isostatic adjustment (GIA) can influence river courses and drainage basins substantially, particularly near former ice margins. We explore the extent to which GIA influenced the evolution of rivers along the United States east coast during the last glacial cycle. We compute gravitationally self-consistent GIA responses that incorporate recent constraints on the Laurentide Ice Sheet history through the last glacial build-up phase, and we connect the predicted variations in topography to abrupt changes in river dynamics recorded in the Hudson, Delaware, Susquehanna, and Potomac Rivers from 40 ka to present. To the extent that increases in sediment transport capacity imply increases in river incision rate, the GIA-driven changes in slope and drainage area are consistent with episodes of erosion and sedimentation observed in the Hudson, Delaware, and Potomac Rivers, but inconsistent with the observed accelerated river incision in the Susquehanna River at 30-14 ka. These analyses add to a growing body of evidence showing that GIA strongly influences river evolution over millennial timescales.
Aims. In this study, we analyzed a filament system, which expanded between moving magnetic features (MMFs) of a decaying sunspot and opposite flux outside of the active region from the nearby quiet-Sun network. This configuration deviated from a classical arch filament system (AFS), which typically connects two pores in an emerging flux region. Thus, we called this system an extended AFS. We contrasted classical and extended AFSs with an emphasis on the complex magnetic structure of the latter. Furthermore, we examined the physical properties of the extended AFS and described its dynamics and connectivity. Methods. The extended AFS was observed with two instruments at the Dunn Solar Telescope (DST). The Rapid Oscillations in the Solar Atmosphere (ROSA) imager provided images in three different wavelength regions, which covered the dynamics of the extended AFS at different atmospheric heights. The Interferometric Bidimensional Spectropolarimeter (IBIS) provided spectroscopic Ha data and spectropolarimetric data that was obtained in the near-infrared (NIR) Call lambda 8542 angstrom line. We derived the corresponding line-of-sight (LOS) velocities and used He II lambda 304 angstrom extreme ultraviolet (EUV) images of the Atmospheric Imaging Assembly (AIA) and LOS magnetograms of the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) as context data. Results. The NIR Call Stokes-V maps are not suitable to definitively define a clear polarity inversion line and to classify this chromospheric structure. Nevertheless, this unusual AFS connects the MMFs of a decaying sunspot with the network field. At the southern footpoint, we measured that the flux decreases over time. We find strong downflow velocities at the footpoints of the extended AFS, which increase in a time period of 30 min. The velocities are asymmetric at both footpoints with higher velocities at the southern footpoint. An EUV brigthening appears in one of the arch filaments, which migrates from the northern footpoint toward the southern one. This activation likely influences the increasing redshift at the southern footpoint. Conclusions. The extended AFS exhibits a similar morphology as classical AFSs, for example, threaded filaments of comparable length and width. Major differences concern the connection from MMFs around the sunspot with the flux of the neighboring quietSun network, converging footpoint motions, and longer lifetimes of individual arch filaments of about one hour, while the extended AFS is still very dynamic.
The intermediate He-enriched hot subdwarf star Feige 46 was recently reported as the second member of the V366 Aqr (or He-sdOBV) pulsating class. Feige 46 is very similar to the prototype of the class, LS IV-14 degrees 116, not only in terms of pulsational properties, but also in terms of atmospheric parameters and kinematic properties. LS IV-14 degrees 116 is additionally characterized by a very peculiar chemical composition, with extreme overabundances of the trans-iron elements Ge, Sr, Y, and Zr. We investigate the possibility that the similarity between the two pulsators extends to their chemical composition. We retrieved archived optical and UV spectroscopic observations of Feige 46 and performed an abundance analysis using model atmospheres and synthetic spectra computed with TLUSTY and SYNSPEC. In total, we derived abundances for 16 elements and provide upper limits for four additional elements. Using absorption lines in the optical spectrum of the star we measure an enrichment of more than 10 000x solar for yttrium and zirconium. The UV spectrum revealed that strontium is equally enriched. Our results confirm that Feige 46 is not only a member of the now growing group of heavy metal subdwarfs, but also has an abundance pattern that is remarkably similar to that of LS IV-14 degrees 116.
Holocene glacial chronostratigraphies in glaciated valleys spread throughout the Himalayan-Tibetan orogen, including the Himalaya, Tibet, Pamir, and Tian Shan, are developed using a landsystems approach, detailed geomorphic mapping, and new and published Be-10 surface exposure dating. New studies in the Kulti valley of Lahul and the Parkachik valley of the Nun Kun massif of the Himalaya of northern India define three glacier advances at similar to 14.7, 12.2, 0.5 ka, in addition to one historically dated late 19th Century advance in the Kulti valley, and one Late Holocene advance at similar to 0.2 ka in the Parkachik valley. Three major climatic groups (subdivided into five climatic zones) are defined across the orogen using Cluster Analysis (CA) and Principal Component Analysis (PCA) to identify glaciated regions with comparable climatic characteristics to evaluate the timing, and extent of Holocene glacier advances across these regions. Our regional analyses across the Himalayan-Tibetan orogen suggest at least one Lateglacial (similar to 15.3-11.8 ka) and five Himalayan-Tibetan Holocene glacial stages (HTHS) at similar to 11.5-9.5, similar to 8.8-7.7, similar to 7.0-3.2, similar to 2.3-1.0, and <1 ka. The extent (amplitude) of glacier advances in 77 glaciated valleys is reconstructed and defined using equilibrium-line altitudes (ELAs). Modern glacier hypsometries are also assessed to help explain the intra-regional variations in glacier amplitudes during each regional glacier advance. A linear inverse glacier flow model is used to decipher the net changes in temperature (Delta T) between periods of reconstructed regional glacier advances in 66 glaciated valleys across different climatic regions throughout the orogen. The Be-10, ELAs, and Delta T data suggest enhanced monsoonal and increased precipitation during the Early Holocene, followed by relative cooling and increased aridity during the Mid- and Late Holocene that influenced glaciation. The sublimation-dominated cold-based glaciers in the northern regions of Himalayan-Tibetan orogen are more affected during these shifts in climate than the temperate glaciers in the south. (C) 2019 Elsevier Ltd. All rights reserved.
Simulations of the inner magnetospheric energetic electrons using the IMPTAM-VERB coupled model
(2019)
In this study, we present initial results of the coupling between the Inner Magnetospheric Particle Transport and Acceleration Model (IMPTAM) and the Versatile Electron Radiation Belt (VERB-3D) code. IMPTAM traces electrons of 10-100 keV energies from the plasma sheet (L = 9 Re) to inner L-shell regions. The flux evolution modeled by IMPTAM is used at the low energy and outer L* computational boundaries of the VERB code (assuming a dipole approximation) to perform radiation belt simulations of energetic electrons. The model was tested on the March 17th, 2013 storm, for a six-day period. Four different simulations were performed and their results compared to satellites observations from Van Allen probes and GOES. The coupled IMPTAM-VERB model reproduces evolution and storm-time features of electron fluxes throughout the studied storm in agreement with the satellite data (within similar to 0.5 orders of magnitude). Including dynamics of the low energy population at L* = 6.6 increases fluxes closer to the heart of the belt and has a strong impact in the VERB simulations at all energies. However, inclusion of magnetopause losses leads to drastic flux decreases even below L* = 3. The dynamics of low energy electrons (max. 10s of keV) do not affect electron fluxes at energies >= 900 keV. Since the IMPTAM-VERB coupled model is only driven by solar wind parameters and the Dst and Kp indexes, it is suitable as a forecasting tool. In this study, we demonstrate that the estimation of electron dynamics with satellite-data-independent models is possible and very accurate.
In this Letter, we propose that the x-ray and the TeV observations in the vicinity of Geminga can be understood in the framework of anisotropic diffusion of injected electrons or positrons. This interpretation only requires the turbulence in the vicinity of Geminga to be sub-Alfvenic with the local mean magnetic field direction approximately aligned with our line of sight towards Geminga, without invoking extreme conditions for the environment, such as an extremely small diffusion coefficient and a weak magnetic field of submicrogauss as suggested in previous literature.
In this paper, we explore the time-energy domain quantum-classical transition comparing a classical pump-probe experiment on a diatomic molecule to its quantum enhanced counterpart, where the pump and probe pulses are substituted by the signal and idler beams of a spontaneous parametric down conversion (SPDC) source. Absorption of biphotons generated with SPDC exploits quantum time-energy entanglement to enhance the overall yield and selectivity of the process, when compared with a classical pump-probe setup, while maintaining femtosecond time resolution. We systematically study the effects of correlation strength on process efficiency and selectivity, comparing the results to classical pump-probe spectra. An excitation scheme to improve the yield based on spectral narrowing of biphotons is shown. The results indicate that the quantum improvements in yield are caused by a more efficient use of the total power available for the process.
X-ray free-electron lasers have, over the past decade, opened up the possibility of understanding the ultrafast response of matter to intense X-ray pulses. In earlier research on atoms and small molecules, new aspects of this response were uncovered, such as rapid sequences of inner-shell photoionization and Auger ionization. Here, we studied a larger molecule, buckminsterfullerene (C-60), exposed to 640 eV X-rays, and examined the role of chemical effects, such as chemical bonds and charge transfer, on the fragmentation following multiple ionization of the molecule. To provide time resolution, we performed femtosecond-resolved X-ray pump/X-ray probe measurements, which were accompanied by advanced simulations. The simulations and experiment reveal that despite substantial ionization induced by the ultrashort (20 fs) X-ray pump pulse, the fragmentation of C-60 is considerably delayed. This work uncovers the persistence of the molecular structure of C-60, which hinders fragmentation over a timescale of hundreds of femtoseconds. Furthermore, we demonstrate that a substantial fraction of the ejected fragments are neutral carbon atoms. These findings provide insights into X-ray free-electron laser-induced radiation damage in large molecules, including biomolecules.
The stability of conducting Taylor-Couette flows under the presence of toroidal magnetic background fields is considered. For strong enough magnetic amplitudes such magnetohydrodynamic flows are unstable against non-axisymmetric perturbations which may also transport angular momentum. In accordance with the often used diffusion approximation, one expects the angular momentum transport to be vanishing for rigid rotation. In the sense of a non-diffusive Lambda effect, however, even for rigidly rotating z-pinches, an axisymmetric angular momentum flux appears which is directed outward (inward) for large (small) magnetic Mach numbers. The internal rotation in a magnetized rotating tank can thus never be uniform. Those particular rotation laws are used to estimate the value of the instability-induced eddy viscosity for which the non-diffusive Lambda effect and the diffusive shear-induced transport compensate each other. The results provide the Shakura & Sunyaev viscosity ansatz leading to numerical values linearly growing with the applied magnetic field.
We present results from our on-going MusE GAs FLOw and Wind (MEGAFLOW) survey, which consists of 22 quasar lines of sight, each observed with the integral field unit MUSE and the UVES spectrograph at the ESO Very Large Telescopes (VLT). The goals of this survey are to study the properties of the circumgalactic medium around z similar to 1 star-forming galaxies. The absorption-line selected survey consists of 79 strong MgII absorbers (with rest-frame equivalent width greater than or similar to 0.3 angstrom) and, currently, 86 associated galaxies within 100 projected kpc of the quasar with stellar masses (M-star) from 109 to 1011 M-circle dot. We find that the cool halo gas traced by MgII is not isotropically distributed around these galaxies from the strong bi-modal distribution in the azimuthal angle of the apparent location of the quasar with respect to the galaxy major axis. This supports a scenario in which outflows are bi-conical in nature and co-exist with a co-planar gaseous structure extending at least up to 60-80 kpc. Assuming that absorbers near the minor axis probe outflows, the current MEGAFLOW sample allowed us to select 26 galaxy-quasar pairs suitable for studying winds. From this sample, using a simple geometrical model, we find that the outflow velocity only exceeds the escape velocity when M-star less than or similar to 4 x 10(9) M-circle dot, implying the cool material is likely to fall back except in the smallest haloes. Finally, we find that the mass loading factor., the ratio between the ejected mass rate and the star formation rate, appears to be roughly constant with respect to the galaxy mass.
We investigate the ultrafast magnetization dynamics of FePt in the L1(0) phase after an optical heating pulse, as used in heat-assisted magnetic recording. We compare continuous and nano-granular thin films and emphasize the impact of the finite size on the remagnetization dynamics. The remagnetization speeds up significantly with increasing external magnetic field only for the continuous film, where domain-wall motion governs the dynamics. The ultrafast remagnetization dynamics in the continuous film are only dominated by heat transport in the regime of high magnetic fields, whereas the timescale required for cooling is prevalent in the granular film for all magnetic field strengths. These findings highlight the necessary conditions for studying the intrinsic heat transport properties in magnetic materials.
Border effect corrections for diagonal line based recurrence quantification analysis measures
(2019)
Recurrence Quantification Analysis (RQA) defines a number of quantifiers, which base upon diagonal line structures in the recurrence plot (RP). Due to the finite size of an RP, these lines can be cut by the borders of the RP and, thus, bias the length distribution of diagonal lines and, consequently, the line based RQA measures. In this letter we investigate the impact of the mentioned border effects and of the thickening of diagonal lines in an RP (caused by tangential motion) on the estimation of the diagonal line length distribution, quantified by its entropy. Although a relation to the Lyapunov spectrum is theoretically expected, the mentioned entropy yields contradictory results in many studies. Here we summarize correction schemes for both, the border effects and the tangential motion and systematically compare them to methods from the literature. We show that these corrections lead to the expected behavior of the diagonal line length entropy, in particular meaning zero values in case of a regular motion and positive values for chaotic motion. Moreover, we test these methods under noisy conditions, in order to supply practical tools for applied statistical research.
Germanium oxide is a prototype network-forming oxide with pressure-induced structural changes similar to those found in crystals and amorphous silicate oxides at high pressure. Studying density and coordination changes in amorphous GeO2 allows for insight into structural changes in silicate oxides at very high pressure, with implications for the properties of planetary magmas. Here, we report the density of germanium oxide glass up to 133 GPa using the x-ray absorption technique, with very good agreement with previous experimental data at pressure below 40 GPa and recent calculation up to 140 GPa. Our data highlight four distinct compressibility domains, corresponding to changes of the local structure of GeO2. Above 80 GPa, our density data show a compressibility and bulk modulus similar to the counterpart crystal phase, and we propose that a compact distorted sixfold coordination, similar to the structural motif of the pyritelike crystalline GeO2 polymorph, is likely to be stable in that pressure range. Our density data point to a smooth continuous evolution of the average coordination for pressure above 20 GPa with persistent sixfold coordination, without sharp density or density slope discontinuities. These observations are in very good agreement with theoretical calculations and spectroscopic measurements, and our results indicate that glasses and melts may behave similarly to their high-pressure solid counterparts with comparable densities, compressibility, and possibly average coordination.
Galaxies are surrounded by massive gas reservoirs ( i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated MilkyWay-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L-star galaxies is extremely diverse: column densities of commonly observed species span similar to 3-4 dex and their covering fractions range from similar to 5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions ( CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L-star galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.
It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the Λ cold dark matter (ΛCDM) model. One type of anisotropy recently detected in the Sloan Digital Sky Survey (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called lopsided. Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set-up when comparing pairs of galaxies in ΛCDM simulations indicating that such a set-up is not inconsistent with ΛCDM, the origin has yet to be explained. Here we examine the origin of such lopsided set-ups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at z = 0 display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards z = 0. Secondly, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the z = 0 signal is driven primarily by satellites that are on first approach, who have yet to experience a ‘flyby’. This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.
We numerically and analytically analyze transitions between different synchronous states in a network of globally coupled phase oscillators with attractive and repulsive interactions. The elements within the attractive or repulsive group are identical, but natural frequencies of the groups differ. In addition to a synchronous two-cluster state, the system exhibits a solitary state, when a single oscillator leaves the cluster of repulsive elements, as well as partially synchronous quasiperiodic dynamics. We demonstrate how the transitions between these states occur when the repulsion starts to prevail over attraction.
On the ALMA observability of nascent massive multiple systems formed by gravitational instability
(2019)
Massive young stellar objects (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (alma) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at similar or equal to a few from the protostar can be resolved by interferometric alma Cycle 7 C43-8 and C43-10 observations at band 6 (), using a maximal 0.015 aracsec beam angular resolution and at least exposure time for sources at distances of . Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.
The concept of a "flow network"-a set of nodes and links which carries one or more flows-unites many different disciplines, including pipe flow, fluid flow, electrical, chemical reaction, ecological, epidemiological, neurological, communications, transportation, financial, economic and human social networks. This Feature Paper presents a generalized maximum entropy framework to infer the state of a flow network, including its flow rates and other properties, in probabilistic form. In this method, the network uncertainty is represented by a joint probability function over its unknowns, subject to all that is known. This gives a relative entropy function which is maximized, subject to the constraints, to determine the most probable or most representative state of the network. The constraints can include "observable" constraints on various parameters, "physical" constraints such as conservation laws and frictional properties, and "graphical" constraints arising from uncertainty in the network structure itself. Since the method is probabilistic, it enables the prediction of network properties when there is insufficient information to obtain a deterministic solution. The derived framework can incorporate nonlinear constraints or nonlinear interdependencies between variables, at the cost of requiring numerical solution. The theoretical foundations of the method are first presented, followed by its application to a variety of flow networks.
The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.
The movement of microswimmers is often described by active Brownian particle models. Here we introduce a variant of these models with several internal states of the swimmer to describe stochastic strategies for directional swimming such as run and tumble or run and reverse that are used by microorganisms for chemotaxis. The model includes a mechanism to generate a directional bias for chemotaxis and interactions with external fields (e.g., gravity, magnetic field, fluid flow) that impose forces or torques on the swimmer. We show how this modified model can be applied to various scenarios: First, the run and tumble motion of E. coli is used to establish a paradigm for chemotaxis and investigate how it is affected by external forces. Then, we study magneto-aerotaxis in magnetotactic bacteria, which is biased not only by an oxygen gradient towards a preferred concentration, but also by magnetic fields, which exert a torque on an intracellular chain of magnets. We study the competition of magnetic alignment with active reorientation and show that the magnetic orientation can improve chemotaxis and thereby provide an advantage to the bacteria, even at rather large inclination angles of the magnetic field relative to the oxygen gradient, a case reminiscent of what is expected for the bacteria at or close to the equator. The highest gain in chemotactic velocity is obtained for run and tumble with a magnetic field parallel to the gradient, but in general a mechanism for reverse motion is necessary to swim against the magnetic field and a run and reverse strategy is more advantageous in the presence of a magnetic torque. This finding is consistent with observations that the dominant mode of directional changes in magnetotactic bacteria is reversal rather than tumbles. Moreover, it provides guidance for the design of future magnetic biohybrid swimmers. Author summary In this paper, we propose a modified Active Brownian particle model to describe bacterial swimming behavior under the influence of external forces and torques, in particular of a magnetic torque. This type of interaction is particularly important for magnetic biohybrids (i.e. motile bacteria coupled to a synthetic magnetic component) and for magnetotactic bacteria (i.e. bacteria with a natural intracellular magnetic chain), which perform chemotaxis to swim along chemical gradients, but are also directed by an external magnetic field. The model allows us to investigate the benefits and disadvantages of such coupling between two different directionality mechanisms. In particular we show that the magnetic torque can speed chemotaxis up in some conditions, while it can hinder it in other cases. In addition to an understanding of the swimming strategies of naturally magnetotactic organisms, the results may guide the design of future biomedical devices.
We utilise multi-epoch MUSE spectroscopy to study binary stars in the core of the Galactic globular cluster NGC 3201. Our sample consists of 3553 stars with 54 883 spectra in total comprising 3200 main-sequence stars up to 4 magnitudes below the turn-off. Each star in our sample has between 3 and 63 (with a median of 14) reliable radial velocity measurements within five years of observations. We introduce a statistical method to determine the probability of a star showing radial velocity variations based on the whole inhomogeneous radial velocity sample. Using HST photometry and an advanced dynamical MOCCA simulation of this specific cluster we overcome observational biases that previous spectroscopic studies had to deal with. This allows us to infer a binary frequency in the MUSE field of view and enables us to deduce the underlying true binary frequency of (6.75 +/- 0.72)% in NGC 3201. The comparison of the MUSE observations with the MOCCA simulation suggests a large portion of primordial binaries. We can also confirm a radial increase in the binary fraction towards the cluster centre due to mass segregation. We discovered that in the core of NGC 3201 at least (57.5 +/- 7.9)% of blue straggler stars are in a binary system. For the first time in a study of globular clusters, we were able to fit Keplerian orbits to a significant sample of 95 binaries. We present the binary system properties of eleven blue straggler stars and the connection to SX Phoenicis-type stars. We show evidence that two blue straggler formation scenarios, the mass transfer in binary (or triple) star systems and the coalescence due to binary-binary interactions, are present in our data. We also describe the binary and spectroscopic properties of four sub-subgiant (or red straggler) stars. Furthermore, we discovered two new black hole candidates with minimum masses (M sin i) of (7.68 +/- 0.50)M-circle dot, (4.4 +/- 2.8)M-circle dot, and refine the minimum mass estimate on the already published black hole to (4.53 +/- 0.21)M-circle dot, These black holes are consistent with an extensive black hole subsystem hosted by NGC 3201.
Context. Galactic globular clusters (GCs) are now known to host multiple populations displaying particular abundance variations. The different populations within a GC can be well distinguished following their position in the pseudo two-colors diagrams, also referred to as "chromosome maps". These maps are constructed using optical and near-UV photometry available from the Hubble Space Telescope (HST) UV survey of GCs. However, the chemical tagging of the various populations in the chromosome maps is hampered by the fact that HST photometry and elemental abundances are both only available for a limited number of stars. Aims. The spectra collected as part of the MUSE survey of globular clusters provide a spectroscopic counterpart to the HST photometric catalogs covering the central regions of GCs. In this paper, we use the MUSE spectra of 1115 red giant branch (RGB) stars in NGC 2808 to characterize the abundance variations seen in the multiple populations of this cluster. Methods. We used the chromosome map of NGC 2808 to divide the RGB stars into their respective populations. We then combined the spectra of all stars belonging to a given population, resulting in one high signal-to-noise ratio spectrum representative of each population. Results. Variations in the spectral lines of O, Na, Mg, and Al are clearly detected among four of the populations. In order to quantify these variations, we measured equivalent width differences and created synthetic populations spectra that were used to determine abundance variations with respect to the primordial population of the cluster. Our results are in good agreement with the values expected from previous studies based on high-resolution spectroscopy. We do not see any significant variations in the spectral lines of Ca, K, and Ba. We also do not detect abundance variations among the stars belonging to the primordial population of NGC 2808. Conclusions. We demonstrate that in spite of their low resolution, the MUSE spectra can be used to investigate abundance variations in the context of multiple populations.
We present new calculations of the mass inflow and outflow rates around the Milky Way (MW), derived from a catalog of ultraviolet metal-line high-velocity clouds (HVCs). These calculations are conducted by transforming the HVC velocities into the Galactic standard of rest (GSR) reference frame, identifying inflowing (vGSR.<.0 km s(-1)) and outflowing (vGSR > 0 km s(-1)) populations, and using observational constraints on the distance, metallicity, dust content, covering fractions, and total silicon column density of each population. After removing HVCs associated with the Magellanic Stream and the Fermi Bubbles, we find inflow and outflow rates in cool (T similar to 10(4) K) ionized gas of dM(in)/dt greater than or similar to.(0.53 +/- 0.23)(d/12 kpc)(Z/0.2Z(circle dot))-1M(circle dot) yr(-1) and dM(out)/dt greater than or similar to (0.16 +/- 0.07)(d/12 kpc)(Z/0.5Z(circle dot))M--1(circle dot) yr(-1). The apparent excess of inflowing over outflowing gas suggests that the MW is currently in an inflow-dominated phase, but the presence of substantial mass flux in both directions supports a Galactic fountain model, in which gas is constantly recycled between the disk and the halo. We also find that the metal flux in both directions (in and out) is indistinguishable. By comparing the outflow rate to the Galactic star formation rate, we present the first estimate of the mass loading factor (eta(HVC)) of the disk-wide MW wind, finding eta(HVC) greater than or similar to (0.10 +/- 0.06)(d/12 kpc)(Z/0.5Z(circle dot))(-1). Including the contributions from low- and intermediatevelocity clouds and from hot gas would increase these inflow and outflow estimates.
Photothermal conversion (PTC) nanostructures have great potential for applications in many fields, and therefore, they have attracted tremendous attention. However, the construction of a PTC nanoreactor with multi-compartment structure to achieve the combination of unique chemical properties and structural feature is still challenging due to the synthetic difficulties. Herein, we designed and synthesized a catalytically active, PTC gold (Au)@polydopamine (PDA) nanoreactor driven by infrared irradiation using assembled PS-b-P2VP nanosphere as soft template. The particles exhibit multi-compartment structure which is revealed by 3D electron tomography characterization technique. They feature permeable shells with tunable shell thickness. Full kinetics for the reduction reaction of 4-nitrophenol has been investigated using these particles as nanoreactors and compared with other reported systems. Notably, a remarkable acceleration of the catalytic reaction upon near-infrared irradiation is demonstrated, which reveals for the first time the importance of the synergistic effect of photothermal conversion and complex inner structure to the kinetics of the catalytic reduction. The ease of synthesis and fresh insights into catalysis will promote a new platform for novel nanoreactor studies.
A comprehensive analysis of the variable temporal and spatial responses of tropical-subtropical high-altitude glaciers to climate change is critical for successful model predictions and environmental risk assessment in the Himalayan-Tibetan orogen. High-frequency Holocene glacier chronostratigraphies are therefore reconstructed in 79 glaciated valleys across the orogen using 519 published and 16 new terrestrial cosmogenic 10Be exposure age dataset. Published 10Be ages are compiled only for moraine boulders (excluding bedrock ages). These ages are recalculated using the latest ICE-D production rate calibration database and the scaling scheme models. Outliers for the individual moraine are detected using the Chauvenet's criterion. In addition, past equilibrium-line altitudes (ELAs) are determined using the area-altitude (AA), area accumulation ratio (AAR), and toe-headwall accumulation ratio (THAR) methods for each glacier advance. The modern maximum elevations of lateral moraines (MELM) are also used to estimate modern ELAs and as an independent check on mean ELAs derived using the above three methods. These data may serve as an essential archive for future studies focusing on the cryospheric and environmental changes in the Himalayan-Tibetan orogen. A more comprehensive analysis of the published and new 10Be ages and ELA results and a list of references are presented in Saha et al. (2019, High-frequency Holocene glacier fluctuations in the Himalayan-Tibetan orogen. Quaternary Science Reviews, 220, 372–400).
EVR-CB-001: An Evolving, Progenitor, White Dwarf Compact Binary Discovered with the Evryscope
(2019)
We present EVR-CB-001, the discovery of a compact binary with an extremely low-mass (0.21 +/- 0.05M(circle dot)) helium core white dwarf progenitor (pre-He WD) and an unseen low-mass (0.32 +/- 0.06M(circle dot)) helium white dwarf (He WD) companion. He WDs are thought to evolve from the remnant helium-rich core of a main-sequence star stripped during the giant phase by a close companion. Low-mass He WDs are exotic objects (only about 0.2% of WDs are thought to be less than 0.3 M-circle dot), and are expected to be found in compact binaries. Pre-He WDs are even rarer, and occupy the intermediate phase after the core is stripped, but before the star becomes a fully degenerate WD and with a larger radius (approximate to 0.2R(circle dot)) than a typical WD. The primary component of EVR-CB-001 (the pre-He WD) was originally thought to be a hot subdwarf (sdB) star from its blue color and under-luminous magnitude, characteristic of sdBs. The mass, temperature (T-eff = 18,500 +/- 500 K), and surface gravity (log(g) = 4.96 +/- 0.04) solutions from this work are lower than values for typical hot subdwarfs. The primary is likely to be a post-red-giant branch, pre-He WD contracting into a He WD, and at a stage that places it nearest to sdBs on color-magnitude and T-eff-log(g) diagrams. EVR-CB-001 is expected to evolve into a fully double degenerate, compact system that should spin down and potentially evolve into a single hot subdwarf star. Single hot subdwarfs are observed, but progenitor systems have been elusive.
Transport in exclusion processes with one-step memory: density dependence and optimal acceleration
(2019)
We study a lattice gas of persistent walkers, in which each site is occupied by at most one particle and the direction each particle attempts to move to depends on its last step. We analyse the mean squared displacement (MSD) of the particles as a function of the particle density and their persistence (the tendency to continue moving in the same direction). For positive persistence the MSD behaves as expected: it increases with the persistence and decreases with the density. However, for strong anti-persistence we find two different regimes, in which the dependence of the MSD on the density is non-monotonic. For very strong anti-persistence there is an optimal density at which the MSD reaches a maximum. In an intermediate regime, the MSD as a function of the density exhibits both a minimum and a maximum, a phenomenon which has not been observed before. We derive a mean-field theory which qualitatively explains this behaviour.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
Plain Language Summary Radiation belts of the Earth, which are the zones of charged energetic particles trapped by the geomagnetic field, comprise enormous and dynamic systems. While the inner radiation belt, composed mainly of high-energy protons, is relatively stable, the outer belt, filled with energetic electrons, is highly variable and depends substantially on solar activity. Hence, extended reliable observations and the improved models of the electron intensities in the outer belt depending on solar wind parameters are necessary for prediction of their dynamics. The Cluster mission has been measuring electron flux intensities in the radiation belts since its launch in 2000, thus providing a huge dataset that can be used for radiation belts analysis. Using 16 years of electron measurements by the Cluster mission corrected for background contamination, we derived a uniform linear-logarithmic dependence of electron fluxes in the outer belt on the solar wind dynamic pressure.
Space radiation is one of the main concerns for human space flights. The prediction of the radiation dose for the actual spacecraft geometry is very important for the planning of long-duration missions. We present a numerical method for the fast calculation of the radiation dose rate during a space flight. We demonstrate its application for dose calculations during the first and the second sessions of the MATROSHKA-R space experiment with a spherical tissue-equivalent phantom. The main advantage of the method is the short simulation time, so it can be applied for urgent radiation dose calculations for low-Earth orbit space missions. The method uses depth-dose curve and shield-and-composition distribution functions to calculate a radiation dose at the point of interest. The spacecraft geometry is processed into a shield-and-composition distribution function using a ray-tracing method. Depth-dose curves are calculated using the GEANT4 Monte-Carlo code (version 10.00.P02) for a double-layer aluminum-water shielding. Aluminum-water shielding is a good approximation of the real geometry, as water is a good equivalent for biological tissues, and aluminum is the major material of spacecraft bodies.
A nova is a cataclysmic event on the surface of a white dwarf in a binary system that increases the overall brightness by several orders of magnitude. Although binary systems with a white dwarf are expected to be overabundant in globular clusters compared with in the Galaxy, only two novae from Galactic globular clusters have been observed. We present the discovery of an emission nebula in the Galactic globular cluster M 22 (NGC 6656) in observations made with the integral-field spectrograph MUSE. We extracted the spectrum of the nebula and used the radial velocity determined from the emission lines to confirm that the nebula is part of NGC 6656. Emission-line ratios were used to determine the electron temperature and density. It is estimated to have a mass of 1-17 x 10(-5) M-circle dot. This mass and the emission-line ratios indicate that the nebula is a nova remnant. Its position coincides with the reported location of a "guest star", an ancient Chinese term for transients, observed in May 48 BCE. With this discovery, this nova may be one of the oldest confirmed extra-solar events recorded in human history.
Short period double degenerate white dwarf (WD) binaries with periods of less than similar to 1 day are considered to be one of the likely progenitors of type Ia supernovae. These binaries have undergone a period of common envelope evolution. If the core ignites helium before the envelope is ejected, then a hot subdwarf remains prior to contracting into a WD. Here we present a comparison of two very rare systems that contain two hot subdwarfs in short period orbits. We provide a quantitative spectroscopic analysis of the systems using synthetic spectra from state-of-the-art non-LTE models to constrain the atmospheric parameters of the stars. We also use these models to determine the radial velocities, and thus calculate dynamical masses for the stars in each system.
Extremely low-mass white dwarf stars (ELMs) are M < 0.3 M-circle dot helium-core white dwarfs born either as a result of a common-envelope phase or after a stable Roche lobe overflow episode in a multiple system. The Universe is not old enough for ELMs to have formed through single-star evolution channels. As remnants of binary evolution, ELMs can shed light onto the poorly understood phase of common-envelope evolution and provide constraints to the physics of mass accretion. Most known ELMs will merge in less than a Hubble time, providing an important contribution to the signal to be detected by upcoming space-based gravitational wave detectors. There are currently less than 150 known ELMs; most were selected by colour, focusing on hot objects, in a magnitude-limited survey of the Northern hemisphere only. Recent theoretical models have predicted a much larger space density for ELMs than estimated observationally based on this limited sample. In order to perform meaningful comparisons with theoretical models and test their predictions, a larger well-defined sample is required. In this work, we present a catalogue of ELM candidates selected from the second data release of Gaia (DR2). We have used predictions from theoretical models and analysed the properties of the known sample to map the space spanned by ELMs in the Gaia Hertzsprung-Russell diagram. Defining a set of colour cuts and quality flags, we have obtained a final sample of 5762 ELM candidates down to T-eff approximate to 5000 K.
Electrospray ionization-ion mobility spectrometry was employed for the determination of collision cross sections (CCS) of 25 synthetically produced peptides in the mass range between 540-3310 Da. The experimental measurement of the CCS is complemented by their calculation applying two different methods. One prediction method is the intrinsic size parameter (ISP) method developed by the Clemmer group. The second new method is based on the evaluation of molecular dynamics (MD) simulation trajectories as a whole, resulting in a single, averaged collision cross-section value for a given peptide in the gas phase. A high temperature MD simulation is run in order to scan through the whole conformational space. The lower temperature conformational distribution is obtained through thermodynamic reweighting. In the first part, various correlations, e.g. CCS vs. mass and inverse mobility vs. m/z correlations, are presented. Differences in CCS between peptides are also discussed in terms of their respective mass and m/z differences, as well as their respective structures. In the second part, measured and calculated CCS are compared. The agreement between the prediction results and the experimental values is in the same range for both calculation methods. While the calculation effort of the ISP method is much lower, the MD method comprises several tools providing deeper insights into the conformations of peptides. Advantages and limitations of both methods are discussed. Based on the separation of two pairs of linear and cyclic peptides of virtually the same mass, the influence of the structure on the cross sections is discussed. The shift in cross section differences and peak shape after transition from the linear to the cyclic peptide can be well understood by applying different MD tools, e.g. the root-mean-square deviation (RMSD) and the root mean square fluctuation (RMSF). (C) 2018 Elsevier B.V. All rights reserved.
Charge extraction methods are popular for measuring the charge carrier density in thin film organic solar cells and to draw conclusions about the order and coefficient of nongeminate charge recombination. However, results from such studies may be falsified by inhomogeneous steady state carrier profiles or surface recombination. Here, we present a detailed drift-diffusion study of two charge extraction methods, bias-assisted charge extraction (BACE) and time-delayed collection field (TDCF). Simulations are performed over a wide range of the relevant parameters. Our simulations reveal that both charge extraction methods provide reliable information about the recombination order and coefficient if the measurements are performed under appropriate conditions. However, results from BACE measurements may be easily affected by surface recombination, in particular for small active layer thicknesses and low illumination densities. TDCF, on the other hand, is more robust against surface recombination due to its transient nature but also because it allows for a homogeneous high carrier density to be inserted into the active layer. Therefore, TDCF is capable to provide meaningful information on the order and coefficient of recombination even if the model conditions are not exactly fulfilled. We demonstrate this for an only 100 nm thick layer of a highly efficient nonfullerene acceptor (NFA) blend, comprising the donor polymer PM6 and the NFA Y6. TDCF measurements were performed as a function of delay time for different laser fluences and bias conditions. The full set of data could be consistently fitted by a strict second order recombination process, with a bias- and fluence-independent bimolecular recombination coefficient k(2) = 1.7 x 10(-17)m(3) s(-1). BACE measurements performed on the very same layer yielded the identical result, despite the very different excitation conditions. This proves that recombination in this blend is mostly through processes in the bulk and that surface recombination is of minor importance despite the small active layer thickness. Published under license by AIP Publishing.
Determinations of the ultraviolet (UV) luminosity function of active galactic nuclei (AGN) at high redshifts are important for constraining the AGN contribution to reionization and understanding the growth of supermassive black holes. Recent inferences of the luminosity function suffer from inconsistencies arising from inhomogeneous selection and analysis of data. We address this problem by constructing a sample of more than 80 000 colour-selected AGN from redshift z= 0 to 7.5 using multiple data sets homogenized to identical cosmologies, intrinsic AGN spectra, and magnitude systems. Using this sample, we derive the AGN UV luminosity function from redshift z= 0 to 7.5. The luminosity function has a double power-law form at all redshifts. The break magnitude M-* shows a steep brightening from M-* similar to -24 at z = 0.7 to M-* similar to -29 at z = 6. The faint-end slope beta significantly steepens from -1.9 at z < 2.2 to -2.4 at z similar or equal to 6. In spite of this steepening, the contribution of AGN to the hydrogen photoionization rate at z similar to 6 is subdominant (< 3 per cent), although it can be non-negligible (similar to 10 per cent) if these luminosity functions hold down to M-1450 = -18. Under reasonable assumptions, AGN can reionize He II by redshift z = 2.9. At low redshifts (z < 0.5), AGN can produce about half of the hydrogen photoionization rate inferred from the statistics of HI absorption lines in the intergalactic medium. Our analysis also reveals important systematic errors in the data, which need to be addressed and incorporated in the AGN selection function in future in order to improve our results. We make various fitting functions, codes, and data publicly available.
Models of ring current electron dynamics unavoidably contain uncertainties in boundary conditions, electric and magnetic fields, electron scattering rates, and plasmapause location. Model errors can accumulate with time and result in significant deviations of model predictions from observations. Data assimilation offers useful tools which can combine physics-based models and measurements to improve model predictions. In this study, we systematically analyze performance of the Kalman filter applied to a log-transformed convection model of ring current electrons and Van Allen Probe data. We consider long-term dynamics of mu = 2.3 MeV/G and K = 0.3 G(1/2) R-E electrons from 1 February 2013 to 16 June 2013. By using synthetic data, we show that the Kalman filter is capable of correcting errors in model predictions associated with uncertainties in electron lifetimes, boundary conditions, and convection electric fields. We demonstrate that reanalysis retains features which cannot be fully reproduced by the convection model such as storm-time earthward propagation of the electrons down to 2.5 R-E. The Kalman filter can adjust model predictions to satellite measurements even in regions where data are not available. We show that the Kalman filter can adjust model predictions in accordance with observations for mu = 0.1, 2.3, and 9.9 MeV/G and constant K = 0.3 G(1/2) R-E electrons. The results of this study demonstrate that data assimilation can improve performance of ring current models, better quantify model uncertainties, and help deeper understand the physics of the ring current particles.
Nanoscale heating by optical excitation of plasmonic nanoparticles offers a new perspective of controlling chemical reactions, where heat is not spatially uniform as in conventional macroscopic heating but strong temperature gradients exist around microscopic hot spots. In nanoplasmonics, metal particles act as a nanosource of light, heat, and energetic electrons driven by resonant excitation of their localized surface plasmon resonance. As an example of the coupling reaction of 4-nitrothiophenol into 4,4′-dimercaptoazobenzene, we show that besides the nanoscopic heat distribution at hot spots, the microscopic distribution of heat dictated by the spot size of the light focus also plays a crucial role in the design of plasmonic nanoreactors. Small sizes of laser spots enable high intensities to drive plasmon-assisted catalysis. This facilitates the observation of such reactions by surface-enhanced Raman scattering, but it challenges attempts to scale nanoplasmonic chemistry up to large areas, where the excess heat must be dissipated by one-dimensional heat transport.
The duration of quasar accretion episodes is a key quantity for distinguishing between models for the formation and growth of supermassive black holes, the evolution of quasars, and their potential feedback effects on their host galaxies. However, this critical time-scale, often referred to as the quasar lifetime, is still uncertain by orders of magnitude (tQ≃0.01Myr−1Gyr). Absorption spectra of quasars exhibiting transmission in the He II Ly α forest provide a unique opportunity to make precise measurements of the quasar lifetime. Indeed, the size of a quasar’s He II proximity zone, the region near the quasar where its own radiation dramatically alters the ionization state of the surrounding intergalactic medium (IGM), depends sensitively on its lifetime for tQ≲30Myr, comparable to the expected e-folding time-scale for SMBH growth tS=45Myr. In this study, we compare the sizes of He II proximity zones in the Hubble Space Telescope (HST) spectra of six z ∼ 4 quasars to theoretical models generated by post-processing cosmological hydrodynamical simulations with a 1D radiative transfer algorithm. We introduce a Bayesian statistical method to infer the lifetimes of individual quasars which allows us to fully marginalize over the unknown ionization state of the surrounding IGM. We measure lifetimes 0.63+0.82−0.40 Myr and 5.75+4.72−2.74 Myr for two objects. For the other four quasars, large redshift uncertainties undermine our sensitivity allowing us to only place upper or lower limits. However, a joint analysis of these four systems yields a measurement of their average lifetime of ⟨tQ⟩=1.17+1.77−0.84 Myr. We discuss our short ∼1Myr inferred lifetimes in the context of other quasar lifetime constraints and the growth of SMBHs.
We review the evidence for a putative early 21st-century divergence between global mean surface temperature (GMST) and Coupled Model Intercomparison Project Phase 5 (CMIP5) projections. We provide a systematic comparison between temperatures and projections using historical versions of GMST products and historical versions of model projections that existed at the times when claims about a divergence were made. The comparisons are conducted with a variety of statistical techniques that correct for problems in previous work, including using continuous trends and a Monte Carlo approach to simulate internal variability. The results show that there is no robust statistical evidence for a divergence between models and observations. The impression of a divergence early in the 21st century was caused by various biases in model interpretation and in the observations, and was unsupported by robust statistics.
We use high-resolution hydrodynamical simulation to test the difference of halo properties in cold dark matter (CDM) and a self-interacting dark matter (SIDM) scenario with a constant cross-section of sigma(T)/m(x) = 1 cm(2) g(-1). We find that the interplay between dark matter self-interaction and baryonic physics induces a complex evolution of the halo properties, which depends on the halo mass and morphological type, as well as on the halo mass accretion history. While high-mass haloes, selected as analogues of early-type galaxies, show cored profiles in the SIDM run, systems of intermediate mass and with a significant disc component can develop a profile that is similar or cuspier than in CDM. The final properties of SIDM haloes - measured at z = 0.2 - correlate with the halo concentration and formation time, suggesting that the differences between different systems are due to the fact that we are observing the impact of self-interaction. We also search for signatures of SIDM in the lensing signal of the main haloes and find hints of potential differences in the distribution of Einstein radii, which suggests that future wide-field survey might be able to distinguish between CDM and SIDM models on this basis. Finally, we find that the subhalo abundances are not altered in the adopted SIDM model with respect to CDM.
Aims. We investigated the ionisation conditions and distances of Galactic high-velocity clouds (HVCs) in the Galactic halo and beyond in the direction of the Local Group (LG) barycentre and anti-barycentre, by studying spectral data of 29 extragalactic background sources obtained with the Cosmic Origins Spectropgraph (COS) installed on the Hubble Space Telescope (HST). Methods. We model column-densities of low, intermediate, and high ions such as Si ii, C ii, Si iii, Si vi, and C iv, and use these data to construct a set of Cloudy ionisation models. Results. In total, we found 69 high-velocity absorption components along the 29 lines of sight. The components in the direction of the LG barycentre span the entire range of studied velocities, 100 less than or similar to vertical bar nu(LSR)vertical bar less than or similar to 400 km s(-1), while those in the anti-barycentre sample have velocities up to about 300 km s(-1). For 49 components, we infer the gas densities. In the direction of the LG barycentre, the gas densities exhibit a wide range from log nH = -3.96 to -2.55, while in the anti-barycentre direction the densities are systematically higher, log nH > -3.25. The barycentre absorbers can be split into two groups based on their density: a high-density group with log nH > -3.54, which can be affected by the Milky Way radiation field, and a low-density group (log nH <= -3.54). The latter has very low thermal pressures of P/k < 7.3 Kcm(-3). Conclusions. Our study shows that part of the absorbers in the LG barycentre direction trace gas at very low gas densities and thermal pressures. These properties indicate that the absorbers are located beyond the virial radius of the Milky Way. Our study also confirms results from earlier, single-sightline studies, suggesting the presence of a metal-enriched intragroup medium filling the LG near its barycentre.
Two new hole selective materials (HSMs) based on dangling methylsulfanyl groups connected to the C-9 position of the fluorene core are synthesized and applied in perovskite solar cells. Being structurally similar to a half of Spiro-OMeTAD molecule, these HSMs (referred as FS and DFS) share similar redox potentials but are endowed with slightly higher hole mobility, due to the planarity and large extension of their structure. Competitive power conversion efficiency (up to 18.6%) is achieved by using the new HSMs in suitable perovskite solar cells. Time-resolved photoluminescence decay measurements and electrochemical impedance spectroscopy show more efficient charge extraction at the HSM/perovskite interface with respect to Spiro-OMeTAD, which is reflected in higher photocurrents exhibited by DFS/FS-integrated perovskite solar cells. Density functional theory simulations reveal that the interactions of methylammonium with methylsulfanyl groups in DFS/FS strengthen their electrostatic attraction with the perovskite surface, providing an additional path for hole extraction compared to the sole presence of methoxy groups in Spiro-OMeTAD. Importantly, the low-cost synthesis of FS makes it significantly attractive for the future commercialization of perovskite solar cells.
We investigate the predictive power of recurrent neural networks for oscillatory systems not only on the attractor but in its vicinity as well. For this, we consider systems perturbed by an external force. This allows us to not merely predict the time evolution of the system but also study its dynamical properties, such as bifurcations, dynamical response curves, characteristic exponents, etc. It is shown that they can be effectively estimated even in some regions of the state space where no input data were given. We consider several different oscillatory examples, including self-sustained, excitatory, time-delay, and chaotic systems. Furthermore, with a statistical analysis, we assess the amount of training data required for effective inference for two common recurrent neural network cells, the long short-term memory and the gated recurrent unit. Published under license by AIP Publishing.
The functionality of organic semiconductor devices crucially depends on molecular energies, namely the ionisation energy and the electron affinity. Ionisation energy and electron affinity values of thin films are, however, sensitive to film morphology and composition, making their prediction challenging. In a combined experimental and simulation study on zinc-phthalocyanine and its fluorinated derivatives, we show that changes in ionisation energy as a function of molecular orientation in neat films or mixing ratio in blends are proportional to the molecular quadrupole component along the p-p-stacking direction. We apply these findings to organic solar cells and demonstrate how the electrostatic interactions can be tuned to optimise the energy of the charge-transfer state at the donor-acceptor interface and the dissociation barrier for free charge carrier generation. The confirmation of the correlation between interfacial energies and quadrupole moments for other materials indicates its relevance for small molecules and polymers.
In this study, rapid loss of relativistic radiation belt electrons at low L* values (2.4-3.2) during a strong geomagnetic storm on 22 June 2015 is investigated along with five possible loss mechanisms. Both the particle and wave data are obtained from the Van Allen Probes. Duskside H+ band electromagnetic ion cyclotron (EMIC) waves were observed during a rapid decrease of relativistic electrons with energy above 5.2 MeV occurring outside the plasma sphere during extreme magnetopause compression. Lower He+ composition and enriched O+ composition are found compared to typical values assumed in other studies of cyclotron resonant scattering of relativistic electrons by EMIC waves. Quantitative analysis demonstrates that even with the existence of He+ band EMIC waves, it is the H+ band EMIC waves that are likely to cause the depletion at small pitch angles and strong gradients in pitch angle distributions of relativistic electrons with energy above 5.2 MeV at low L values for this event. Very low frequency wave activity at other magnetic local time can be favorable for the loss of relativistic electrons at higher pitch angles. An illustrative calculation that combines the nominal pitch angle scattering rate due to whistler mode chorus at high pitch angles with the H+ band EMIC wave loss rate at low pitch angles produces loss on time scale observed at L = 2.4-3.2. At high L values and lower energies, radial loss to the magnetopause is a viable explanation.
Storm Time Depletions of Multi-MeV Radiation Belt Electrons Observed at Different Pitch Angles
(2019)
During geomagnetic storms, the rapid depletion of the high-energy (several MeV) outer radiation belt electrons is the result of loss to the interplanetary medium through the magnetopause, outward radial diffusion, and loss to the atmosphere due to wave-particle interactions. We have performed a statistical study of 110 storms using pitch angle resolved electron flux measurements from the Van Allen Probes mission and found that inside of the radiation belt (L* = 3 - 5) the number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 30 degrees is higher than number of storms that result in depletion of electrons with equatorial pitch angle alpha(eq) = 75 degrees. We conclude that this result is consistent with electron scattering by whistler and electromagnetic ion cyclotron waves. At the outer edge of the radiation belt (L* >= 5.2) the number of storms that result in depletion is also large (similar to 40-50%), emphasizing the significance of the magnetopause shadowing effect and outward radial transport.
We consider a simple model for active random walk with general temporal correlations, and investigate the shape of the probability distribution function of the displacement during a short time interval. We find that under certain conditions the distribution exhibits multiple peaks and we show analytically and numerically that the existence of these peaks is governed by the walker?s tendency to move forward, while the correlations between the timing of its active motion control the magnitude and shape of the peaks. In particular, we find that in a homogeneous system such peaks can occur only if the persistence is strong enough.
Polar crown filaments form above the polarity inversion line between the old magnetic flux of the previous cycle and the new magnetic flux of the current cycle. Studying their appearance and their properties can lead to a better understanding of the solar cycle. We use full-disk data of the Chromospheric Telescope (ChroTel) at the Observatorio del Teide, Tenerife, Spain, which were taken in three different chromospheric absorption lines (H alpha lambda 6563 angstrom, Ca II K lambda 3933 angstrom, and He I lambda 10830 angstrom), and we create synoptic maps. In addition, the spectroscopic He I data allow us to compute Doppler velocities and to create synoptic Doppler maps. ChroTel data cover the rising and decaying phase of Solar Cycle 24 on about 1000 days between 2012 and 2018. Based on these data, we automatically extract polar crown filaments with image-processing tools and study their properties. We compare contrast maps of polar crown filaments with those of quiet-Sun filaments. Furthermore, we present a super-synoptic map summarizing the entire ChroTel database. In summary, we provide statistical properties, i.e. number and location of filaments, area, and tilt angle for both the maximum and the declining phase of Solar Cycle 24. This demonstrates that ChroTel provides a promising data set to study the solar cycle.
A stellar census in globular clusters with MUSE: A spectral catalogue of emission-line sources
(2019)
Aims. Globular clusters produce many exotic stars due to a much higher frequency of dynamical interactions in their dense stellar environments. Some of these objects were observed together with several hundred thousand other stars in our MUSE survey of 26 Galactic globular clusters. Assuming that at least a few exotic stars have exotic spectra (i.e. spectra that contain emission lines), we can use this large spectroscopic data set of over a million stellar spectra as a blind survey to detect stellar exotica in globular clusters. Methods. To detect emission lines in each spectrum, we modelled the expected shape of an emission line as a Gaussian curve. This template was used for matched filtering on the di fferences between each observed 1D spectrum and its fitted spectral model. The spectra with the most significant detections of H alpha emission are checked visually and cross-matched with published catalogues. Results. We find 156 stars with H alpha emission, including several known cataclysmic variables (CV) and two new CVs, pulsating variable stars, eclipsing binary stars, the optical counterpart of a known black hole, several probable sub-subgiants and red stragglers, and 21 background emission-line galaxies. We find possible optical counterparts to 39 X-ray sources, as we detected H alpha emission in several spectra of stars that are close to known positions of Chandra X-ray sources. This spectral catalogue can be used to supplement existing or future X-ray or radio observations with spectra of potential optical counterparts to classify the sources.
The actin cytoskeleton and its response to external chemical stimuli is fundamental to the mechano-biology of eukaryotic cells and their functions. One of the key players that governs the dynamics of the actin network is the motor protein myosin II. Based on a phase space embedding we have identified from experiments three phases in the cytoskeletal dynamics of starved Dictyostelium discoideum in response to a precisely controlled chemotactic stimulation. In the first two phases the dynamics of actin and myosin II in the cortex is uncoupled, while in the third phase the time scale for the recovery of cortical actin is determined by the myosin II dynamics. We report a theoretical model that captures the experimental observations quantitatively. The model predicts an increase in the optimal response time of actin with decreasing myosin II-actin coupling strength highlighting the role of myosin II in the robust control of cell contraction.
Electron acceleration at Saturn due to whistler mode chorus waves has previously been assumed to be ineffective; new data closer to the planet show it can be very rapid (factor of 104 flux increase at 1 MeV in 10 days compared to factor of 2). A full survey of chorus waves at Saturn is combined with an improved plasma density model to show that where the plasma frequency falls below the gyrofrequency additional strong resonances are observed favoring electron acceleration. This results in strong chorus acceleration between approximately 2.5 R-S and 5.5 R-S outside which adiabatic transport may dominate. Strong pitch angle dependence results in butterfly pitch angle distributions that flatten over a few days at 100s keV, tens of days at MeV energies which may explain observations of butterfly distributions of MeV electrons near L = 3. Including cross terms in the simulations increases the tendency toward butterfly distributions. Plain Language Summary Radiation belts are hazardous regions found around several of the planets in our Solar System. They consist of very hot, electrically charged particles trapped in the magnetic field of the planet. At Saturn the most important way to heat these particles has for many years been thought to involve the particles drifting closer toward the planet. This paper adds to the emerging idea at Saturn that a different way to heat the particles is also possible where the heating is done by waves, in a similar way to what we find at the Earth. We use recent information from the Cassini spacecraft on the number and location of particles and also of the waves strength and location combined with computer simulations to show that a particular wave called chorus is excellent at heating the particles where the surrounding number of cold particles is low.
The electric field-dependence of structural dynamics in a tetragonal ferroelectric lead zirconate titanate thin film is investigated under subcoercive and above-coercive fields using time-resolved X-ray diffraction. The domain nucleation and growth are monitored in real time during the application of an external field to the prepoled thin film capacitor. We propose the observed broadening of the in-plane peak width of the symmetric 002 Bragg reflection as an indicator of the domain disorder and discuss the processes that change the measured peak intensity. Subcoercive field switching results in remnant disordered domain configurations. Published under license by AIP Publishing.
We present measurements of the large-scale (≈40 comoving Mpc) effective optical depth of He ii Lyα absorption, ${\tau }_{\mathrm{eff}}$, at 2.54 < z < 3.86 toward 16 He ii-transparent quasars observed with the Cosmic Origins Spectrograph on the Hubble Space Telescope, to characterize the ionization state of helium in the intergalactic medium (IGM). We provide the first statistical sample of ${\tau }_{\mathrm{eff}}$ measurements in six signal-to-noise ratio gsim3 He ii sightlines at z > 3.5, and study the redshift evolution and sightline-to-sightline variance of ${\tau }_{\mathrm{eff}}$ in 24 He ii sightlines. We confirm an increase of the median ${\tau }_{\mathrm{eff}}$ from sime2 at z = 2.7 to ${\tau }_{\mathrm{eff}}\gtrsim 5$ at z > 3, and a scatter in ${\tau }_{\mathrm{eff}}$ that increases with redshift. The z > 3.5 He ii absorption is predominantly saturated, but isolated narrow (Δv < 650 km s−1) transmission spikes indicate patches of reionized helium. We compare our measurements to predictions for a range of UV background models applied to outputs of a large-volume (146 comoving Mpc)3 hydrodynamical simulation by forward-modeling our sample's quality and size. At z > 2.74, the variance in ${\tau }_{\mathrm{eff}}$ significantly exceeds expectations for a spatially uniform UV background, but is consistent with a fluctuating radiation field sourced by variations in the quasar number density and the mean free path in the post-reionization IGM. We develop a method to infer the approximate median He ii photoionization rate ${{\rm{\Gamma }}}_{\mathrm{He}{\rm{II}}}$ of a fluctuating UV background from the median ${\tau }_{\mathrm{eff}}$, finding a factor sime5 decrease in ${{\rm{\Gamma }}}_{\mathrm{He}{\rm{II}}}$ between z sime 2.6 and z sime 3.1. At z sime 3.1, ${{\rm{\Gamma }}}_{\mathrm{He}{\rm{II}}}=\left[{9.1}_{-1.2}^{+1.1}\,(\mathrm{stat}.){\,}_{-3.4}^{+2.4}\,(\mathrm{sys}.)\right]\times {10}^{-16}$ s−1 corresponds to a median He ii fraction of sime2.5%, indicating that our data probe the tail end of He ii reionization.
Intermolecular charge-transfer states at the interface between electron donating (D) and accepting (A) materials are crucial for the operation of organic solar cells but can also be exploited for organic light-emitting diodes(1,2). Non-radiative charge-transfer state decay is dominant in state-of-the-art D-A-based organic solar cells and is responsible for large voltage losses and relatively low power-conversion efficiencies as well as electroluminescence external quantum yields in the 0.01-0.0001% range(3,4). In contrast, the electroluminescence external quantum yield reaches up to 16% in D-A-based organic light-emitting diodes(5-7). Here, we show that proper control of charge-transfer state properties allows simultaneous occurrence of a high photovoltaic and emission quantum yield within a single, visible-light-emitting D-A system. This leads to ultralow-emission turn-on voltages as well as significantly reduced voltage losses upon solar illumination. These results unify the description of the electro-optical properties of charge-transfer states in organic optoelectronic devices and foster the use of organic D-A blends in energy conversion applications involving visible and ultraviolet photons(8-11).
Wealth and income distributions are known to feature country-specific Pareto exponents for their long power-law tails. To propose a rationale for this, we introduce an agent-based dynamic model and use Monte Carlo simulations to unveil the wealth distributions in closed and open economical systems. The standard money-exchange scenario is supplemented with the position-exchange agent dynamics that vitally affects the Pareto law. Specifically, in closed systems with position-exchange dynamics the power law changes to an exponential shape, while for open systems with traps the Pareto law remains valid.
Pitch Angle Scattering of Sub-MeV Relativistic Electrons by Electromagnetic Ion Cyclotron Waves
(2019)
Electromagnetic ion cyclotron (EMIC) waves have long been considered to be a significant loss mechanism for relativistic electrons. This has most often been attributed to resonant interactions with the highest amplitude waves. But recent observations have suggested that the dominant energy of electrons precipitated to the atmosphere may often be relatively low, less than 1 MeV, whereas the minimum resonant energy of the highest amplitude waves is often greater than 2 MeV. Here we use relativistic electron test particle simulations in the wavefields of a hybrid code simulation of EMIC waves in dipole geometry in order to show that significant pitch angle scattering can occur due to interaction with low-amplitude short-wavelength EMIC waves. In the case we examined, these waves are in the H band (at frequencies above the He+ gyrofrequency), even though the highest amplitude waves were in the He band frequency range (below the He+ gyrofrequency). We also present wave power distributions for 29 EMIC simulations in straight magnetic field line geometry that show that the high wave number portion of the spectrum is in every case mostly due to the H band waves. Though He band waves are often associated with relativistic electron precipitation, it is possible that the He band waves do not directly scatter the sub-megaelectron volts (sub-MeV) electrons, but that the presence of He band waves is associated with high plasma density which lowers the minimum resonant energy so that these electrons can more easily resonate with the H band waves.
In magnetized plasmas such as the ionosphere, electric currents develop in regions of strong density gradients to balance the resulting plasma pressure gradients. These currents, usually known as diamagnetic currents decrease the magnetic pressure where the plasma pressure increases, and vice versa. In the low‐latitude ionosphere, equatorial plasma depletions (EPDs) are well known for their steep plasma density gradients and adverse effect on radio wave propagation. In this paper, we use continuous measurements of the magnetic field and electron density from the European Space Agency's Swarm constellation mission to assess the balance between plasma and magnetic pressure across large‐scale EPDs. The analysis is based on the magnetic fluctuations related to diamagnetic currents flowing at the edges of EPDs. This study shows that most of the EPDs detected by Swarm present a decrease of the plasma pressure relative to the ambient plasma. However, EPDs with high plasma pressure are also identified mainly in the vicinity of the South Atlantic magnetic anomaly. From the electron density measurements, we deduce that such an increase in plasma pressure within EPDs might be possible by temperatures inside the EPD as high as twice the temperature of the ambient plasma. Due to the distinct location of the high‐pressure EPDs, we suggest that a possible heating mechanism might be due to precipitation of particle from the radiation belts. This finding corresponds to the first observational evidence of plasma pressure enhancements in regions of depleted plasma density in the ionosphere.
Modeling and observations have shown that energy diffusion by chorus waves is an important source of acceleration of electrons to relativistic energies. By performing long-term simulations using the three-dimensional Versatile Electron Radiation Belt code, in this study, we test how the latitudinal dependence of chorus waves can affect the dynamics of the radiation belt electrons. Results show that the variability of chorus waves at high latitudes is critical for modeling of megaelectron volt (MeV) electrons. We show that, depending on the latitudinal distribution of chorus waves under different geomagnetic conditions, they cannot only produce a net acceleration but also a net loss of MeV electrons. Decrease in high-latitude chorus waves can tip the balance between acceleration and loss toward acceleration, or alternatively, the increase in high-latitude waves can result in a net loss of MeV electrons. Variations in high-latitude chorus may account for some of the variability of MeV electrons.
We demonstrate full control of acoustic and thermal periodic deformations at solid surfaces down to subnanosecond time scales and few-micrometer length scales via independent variation of the temporal and spatial phase of two optical transient grating (TG) excitations. For this purpose, we introduce an experimental setup that exerts control of the spatial phase of subsequent time-delayed TG excitations depending on their polarization state. Specific exemplary coherent control cases are discussed theoretically and corresponding experimental data are presented in which time-resolved x-ray reflectivity measures the spatiotemporal surface distortion of nanolayered heterostructures. Finally, we discuss examples where the application of our method may enable the control of functional material properties via tailored spatiotemporal strain fields.
The changing temperature of the nucleus of comet 67P induced by morphological and seasonal effects
(2019)
Knowledge of the surface temperature distribution on a comet’s nucleus and its temporal evolution at different timescales is key to constraining its thermophysical properties and understanding the physical processes that take place at and below the surface. Here we report on time-resolved maps of comet 67P/Churyumov–Gerasimenko retrieved on the basis of infrared data acquired by the Visible InfraRed and Thermal Imaging Spectrometer (VIRTIS) onboard the Rosetta orbiter in 2014, over a roughly two-month period in the pre-perihelion phase at heliocentric distances between 3.62 and 3.31 au from the Sun. We find that at a spatial resolution ≤15 m per pixel, the measured temperatures point out the major effect that self-heating, due to the complex shape of the nucleus, has on the diurnal temperature variation. The bilobate nucleus of comet 67P also induces daytime shadowing effects, which result in large thermal gradients. Over longer periods, VIRTIS-derived temperature values reveal seasonal changes driven by decreasing heliocentric distance combined with an increasing abundance of ice within the uppermost centimetre-thick layer, which implies the possibility of having a largely pristine nucleus interior already in the shallow subsurface
Pump-probe photoelectron spectroscopy (PES) is a versatile tool to investigate the dynamics of transient states of excited matter. Vacuum space-charge effects can mask these dynamics and complicate the interpretation of electron spectra. Here we report on space-charge effects in Au 4f photoemission from a polycrystalline gold surface, excited with moderately intense 90 ps (FWHM) soft x-ray probe pulses, under the influence of the Coulomb forces exerted by a pump electron cloud, which was produced by intense 40 fs laser pulses. The experimentally observed kinetic energy shift and spectral broadening of the Au 4f lines, measured with highly-efficient time-of-flight spectroscopy, are in good agreement with simulations utilizing a mean-field model of the electrostatic pump electron potential. This confirms that the line broadening is predominantly caused by variations in the take-off time of the probe electrons without appreciable influence of local scattering events. Our findings might be of general interest for pump-probe PES with picosecond-pulse-length sources.
By conducting helicity-dependent ultrafast magnetization dynamics in a CoTb ferrimagnetic alloy, we are able to quantitatively determine the magnetic circular dichroism (MCD) and resolve its role in the helicity-dependent all-optical switching (AOS). Unequivocal interpretation of the sign of the dichroism is provided by performing AOS and femtosecond laser-induced domain wall motion experiments. We demonstrate that AOS occurs when the magnetization is initially in the most absorbent state, according to the light helicity. Moreover, we evidence that the MCD creates a thermal gradient that drives a domain wall toward hotter regions. Our experimental results are in agreement with the purely thermal models of AOS.
Saturn’s main ring system is associated with a set of small moons that either are embedded within it or interact with the rings to alter their shape and composition. Five close flybys of the moons Pan, Daphnis, Atlas, Pandora, and Epimetheus were performed between December 2016 and April 2017 during the ring-grazing orbits of the Cassini mission. Data on the moons’ morphology, structure, particle environment, and composition were returned, along with images in the ultraviolet and thermal infrared. We find that the optical properties of the moons’ surfaces are determined by two competing processes: contamination by a red material formed in Saturn’s main ring system and accretion of bright icy particles or water vapor from volcanic plumes originating on the moon Enceladus.
The optimized design of multilayer-coated blazed gratings (MLBG) for high-flux tender X-ray monochromators was systematically studied by numerical simulations. The resulting correlation between the multilayer d-spacing and grating blaze angle significantly deviated from the one predicted by conventional equations. Three high line density gratings with different blaze angles were fabricated and coated by the same Cr/C multilayer. The MLBG with an optimal blaze angle of 1.0 degrees showed a record efficiency reaching 60% at 3.1 keV and 4.1 keV. The measured efficiencies of all three gratings were consistent with calculated results proving the validity of the numerical simulation and indicating a more rigorous way to design the optimal MLBG structure. Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License.
An understanding of the factors limiting the open-circuit voltage (V-oc) and related photon energy loss mechanisms is critical to increase the power conversion efficiency (PCE) of small-molecule organic solar cells (OSCs), especially those with near-infrared (NIR) absorbers. In this work, two NIR boron dipyrromethene (BODIPY) molecules are characterized for application in planar (PHJ) and bulk (BHJ) heterojunction OSCs. When two H atoms are substituted by F atoms on the peripheral phenyl rings of the molecules, the molecular aggregation type in the thin film changes from the H-type to J-type. For PHJ devices, the nonradiative voltage loss of 0.35 V in the J-aggregated BODIPY is lower than that of 0.49 V in the H-aggregated device. In BHJ devices with a nonradiative voltage loss of 0.35 V, a PCE of 5.5% is achieved with an external quantum efficiency (EQE) maximum of 68% at 700 nm.
The tremendous success of metal-halide perovskites, especially in the field of photovoltaics, has triggered a substantial number of studies in understanding their optoelectronic properties. However, consensus regarding the electronic properties of these perovskites is lacking due to a huge scatter in the reported key parameters, such as work function (Φ) and valence band maximum (VBM) values. Here, we demonstrate that the surface photovoltage (SPV) is a key phenomenon occurring at the perovskite surfaces that feature a non-negligible density of surface states, which is more the rule than an exception for most materials under study. With ultraviolet photoelectron spectroscopy (UPS) and Kelvin probe, we evidence that even minute UV photon fluxes (500 times lower than that used in typical UPS experiments) are sufficient to induce SPV and shift the perovskite Φ and VBM by several 100 meV compared to dark. By combining UV and visible light, we establish flat band conditions (i.e., compensate the surface-state-induced surface band bending) at the surface of four important perovskites, and find that all are p-type in the bulk, despite a pronounced n-type surface character in the dark. The present findings highlight that SPV effects must be considered in all surface studies to fully understand perovskites’ photophysical properties.