Refine
Year of publication
- 2016 (258) (remove)
Document Type
- Article (203)
- Doctoral Thesis (21)
- Postprint (18)
- Other (10)
- Review (4)
- Habilitation Thesis (2)
Is part of the Bibliography
- yes (258)
Keywords
- Magellanic Clouds (8)
- Sun: magnetic fields (6)
- intergalactic medium (6)
- methods: data analysis (6)
- stars: Wolf-Rayet (6)
- techniques: spectroscopic (6)
- galaxies: formation (5)
- stars: atmospheres (5)
- turbulence (5)
- X-rays: binaries (4)
Institute
- Institut für Physik und Astronomie (258) (remove)
Galactic cosmic rays reach energies of at least a few petaelectronvolts (of the order of 1015 electronvolts). This implies that our Galaxy contains petaelectronvolt accelerators (‘PeVatrons’), but all proposed models of Galactic cosmic-ray accelerators encounter difficulties at exactly these energies. Dozens of Galactic accelerators capable of accelerating particles to energies of tens of teraelectronvolts (of the order of 1013 electronvolts) were inferred from recent γ-ray observations3. However, none of the currently known accelerators—not even the handful of shell-type supernova remnants commonly believed to supply most Galactic cosmic rays—has shown the characteristic tracers of petaelectronvolt particles, namely, power-law spectra of γ-rays extending without a cut-off or a spectral break to tens of teraelectronvolts4. Here we report deep γ-ray observations with arcminute angular resolution of the region surrounding the Galactic Centre, which show the expected tracer of the presence of petaelectronvolt protons within the central 10 parsecs of the Galaxy. We propose that the supermassive black hole Sagittarius A* is linked to this PeVatron. Sagittarius A* went through active phases in the past, as demonstrated by X-ray outbursts5and an outflow from the Galactic Centre6. Although its current rate of particle acceleration is not sufficient to provide a substantial contribution to Galactic cosmic rays, Sagittarius A* could have plausibly been more active over the last 106–107 years, and therefore should be considered as a viable alternative to supernova remnants as a source of petaelectronvolt Galactic cosmic rays.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
The very high energy (VHE; E > 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54 938-54 956), 2009, as part of an extensive multiwavelength campaign from radio to VHE. Strong VHE yray activity was detected on May 1st with Whipple and VERITAS, when the flux (E > 400 GeV) increased to 10 times the preflare baseline flux (3.9 x 10(-11) ph cm(-2) s(-1)), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15. This VHE flare showed a fast flux variation with an increase of a factor similar to 4 in 25 min, and a falling time of similar to 50 min. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.
A SEARCH FOR VERY HIGH ENERGY GAMMA RAYS FROM THE MISSING LINK BINARY PULSAR J1023+0038 WITH VERITAS
(2016)
The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than similar to 2 G before the disappearance of the radio pulsar and greater than similar to 10 G afterward.
Gravity dictates the structure of the whole Universe and, although it is triumphantly described by the theory of General Relativity, it is the force that we least understand in nature. One of the cardinal predictions of this theory are black holes. Massive, dark objects are found in the majority of galaxies. Our own galactic center very contains such an object with a mass of about four million solar masses. Are these objects supermassive black holes (SMBHs), or do we need alternatives? The answer lies in the event horizon, the characteristic that defines a black hole. The key to probe the horizon is to model the movement of stars around a SMBH, and the interactions between them, and look for deviations from real observations. Nuclear star clusters harboring a massive, dark object with a mass of up to ~ ten million solar masses are good testbeds to probe the event horizon of the potential SMBH with stars. The channel for interactions between stars and the central MBH are the fact that (a) compact stars and stellar-mass black holes can gradually inspiral into the SMBH due to the emission of gravitational radiation, which is known as an “Extreme Mass Ratio Inspiral” (EMRI), and (b) stars can produce gases which will be accreted by the SMBH through normal stellar evolution, or by collisions and disruptions brought about by the strong central tidal field. Such processes can contribute significantly to the mass of the SMBH. These two processes involve different disciplines, which combined will provide us with detailed information about the fabric of space and time. In this habilitation I present nine articles of my recent work directly related with these topics.
The TeV binary system LS I +61 degrees 303 is known for its regular, non-thermal emission pattern that traces the orbital period of the compact object in its 26.5 day orbit around its B0 Ve star companion. The system typically presents elevated TeV emission around apastron passage with flux levels between 5% and 15% of the steady flux from the Crab Nebula (> 300 GeV). In this article, VERITAS observations of LS I + 61 degrees. 303 taken in late 2014 are presented, during which bright TeV flares around apastron at flux levels peaking above 30% of the Crab Nebula flux were detected. This is the brightest such activity from this source ever seen in the TeV regime. The strong outbursts have rise and fall times of less than a day. The short timescale of the flares, in conjunction with the observation of 10 TeV photons from LS I + 61 degrees 303 during the flares, provides constraints on the properties of the accelerator in the source.
Between the beginning of its full-scale scientific operations in 2007 and 2012, the VERITAS Cherenkov telescope array observed more than 130 blazars; of these, 26 were detected as very-high-energy (VHE; E > 100 GeV) gamma-ray sources. In this work, we present the analysis results of a sample of 114 undetected objects. The observations constitute a total live-time of similar to 570 hr. The sample includes several unidentified Fermi-Large Area Telescope (LAT) sources (located at high Galactic latitude) as well as all the sources from the second Fermi-LAT catalog that are contained within the field of view of the VERITAS observations. We have also performed optical spectroscopy measurements in order to estimate the redshift of some of these blazars that do not have spectroscopic distance estimates. We present new optical spectra from the Kast instrument on the Shane telescope at the Lick observatory for 18 blazars included in this work, which allowed for the successful measurement or constraint on the redshift of four of them. For each of the blazars included in our sample, we provide the flux upper limit in the VERITAS energy band. We also study the properties of the significance distributions and we present the result of a stacked analysis of the data set, which shows a 4s excess.
The Galactic Center ridge has been observed extensively in the past by both GeV and TeV gamma-ray instruments revealing a wealth of structure, including a diffuse component and the point sources G0.9+0.1 (a composite supernova remnant) and Sgr A* (believed to be associated with the supermassive black hole located at the center of our Galaxy). Previous very high energy (VHE) gamma-ray observations with the H.E.S.S.. experiment have also detected an extended TeV gamma-ray component along the Galactic plane in the >300 GeV gamma-ray regime. Here we report on observations of the Galactic Center ridge from 2010 to 2014 by the VERITAS telescope array in the >2 TeV energy range. From these observations we (1) provide improved measurements of the differential energy spectrum for Sgr A* in the >2 TeV gamma-ray regime, (2) provide a detection in the >2 TeV gamma-ray emission from the composite SNR G0.9+0.1 and an improved determination of its multi-TeV gamma-ray energy spectrum, and. (3) report on the detection of VER J1746-289, a localized enhancement of >2 TeV gamma-ray emission along the Galactic plane.
Methods. In most cases, the given orientation of the ecliptic is used to set up the heliographic coordinate system for the drawings. Positions and sizes are measured manually on screen. Very early drawings have no indication of their orientation. A rotational matching using common spots of adjacent days is used in some cases, while in other cases, the assumption that images were aligned with a zenith-horizon coordinate system appeared to be the most probable.
In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.
In the current paradigm of cosmology, the formation of large-scale structures is mainly driven by non-radiating dark matter, making up the dominant part of the matter budget of the Universe. Cosmological observations however, rely on the detection of luminous galaxies, which are biased tracers of the underlying dark matter. In this thesis I present cosmological reconstructions of both, the dark matter density field that forms the cosmic web, and cosmic velocities, for which both aspects of my work are delved into, the theoretical formalism and the results of its applications to cosmological simulations and also to a galaxy redshift survey.The foundation of our method is relying on a statistical approach, in which a given galaxy catalogue is interpreted as a biased realization of the underlying dark matter density field. The inference is computationally performed on a mesh grid by sampling from a probability density function, which describes the joint posterior distribution of matter density and the three dimensional velocity field. The statistical background of our method is described in Chapter ”Implementation of argo”, where the introduction in sampling methods is given, paying special attention to Markov Chain Monte-Carlo techniques. In Chapter ”Phase-Space Reconstructions with N-body Simulations”, I introduce and implement a novel biasing scheme to relate the galaxy number density to the underlying dark matter, which I decompose into a deterministic part, described by a non-linear and scale-dependent analytic expression, and a stochastic part, by presenting a negative binomial (NB) likelihood function that models deviations from Poissonity. Both bias components had already been studied theoretically, but were so far never tested in a reconstruction algorithm. I test these new contributions againstN-body simulations to quantify improvements and show that, compared to state-of-the-art methods, the stochastic bias is inevitable at wave numbers of k≥0.15h Mpc^−1 in the power spectrum in order to obtain unbiased results from the reconstructions. In the second part of Chapter ”Phase-Space Reconstructions with N-body Simulations” I describe and validate our approach to infer the three dimensional cosmic velocity field jointly with the dark matter density. I use linear perturbation theory for the large-scale bulk flows and a dispersion term to model virialized galaxy motions, showing that our method is accurately recovering the real-space positions of the redshift-space distorted galaxies. I analyze the results with the isotropic and also the two-dimensional power spectrum.Finally, in Chapter ”Phase-space Reconstructions with Galaxy Redshift Surveys”, I show how I combine all findings and results and apply the method to the CMASS (for Constant (stellar) Mass) galaxy catalogue of the Baryon Oscillation Spectroscopic Survey (BOSS). I describe how our method is accounting for the observational selection effects inside our reconstruction algorithm. Also, I demonstrate that the renormalization of the prior distribution function is mandatory to account for higher order contributions in the structure formation model, and finally a redshift-dependent bias factor is theoretically motivated and implemented into our method. The various refinements yield unbiased results of the dark matter until scales of k≤0.2 h Mpc^−1in the power spectrum and isotropize the galaxy catalogue down to distances of r∼20h^−1 Mpc in the correlation function. We further test the results of our cosmic velocity field reconstruction by comparing them to a synthetic mock galaxy catalogue, finding a strong correlation between the mock and the reconstructed velocities. The applications of both, the density field without redshift-space distortions, and the velocity reconstructions, are very broad and can be used for improved analyses of the baryonic acoustic oscillations, environmental studies of the cosmic web, the kinematic Sunyaev-Zel’dovic or integrated Sachs-Wolfe effect.
Arch filament systems occur in active sunspot groups, where a fibril structure connects areas of opposite magnetic polarity, in contrast to active region filaments that follow the polarity inversion line. We used the GREGOR Infrared Spectrograph (GRIS) to obtain the full Stokes vector in the spectral lines SiI lambda 1082.7 nm, He I lambda 1083.0 nm, and Ca I lambda 1083.9 nm. We focus on the near-infrared calcium line to investigate the photospheric magnetic field and velocities, and use the line core intensities and velocities of the helium line to study the chromospheric plasma. The individual fibrils of the arch filament system connect the sunspot with patches of magnetic polarity opposite to that of the spot. These patches do not necessarily coincide with pores, where the magnetic field is strongest. Instead, areas are preferred not far from the polarity inversion line. These areas exhibit photospheric downflows of moderate velocity, but significantly higher downflows of up to 30 km s(-1) in the chromospheric helium line. Our findings can be explained with new emerging flux where the matter flows downward along the field lines of rising flux tubes, in agreement with earlier results. (C) 2016 WILEY-VCH Verlag GmbH& Co. KGaA, Weinheim
Traditional economic theory could not explain, much less predict, the near collapse of the financial system and its long-lasting effects on the global economy. Since the 2008 crisis, there has been increasing interest in using ideas from complexity theory to make sense of economic and financial markets. Concepts, such as tipping points, networks, contagion, feedback, and resilience have entered the financial and regulatory lexicon, but actual use of complexity models and results remains at an early stage. Recent insights and techniques offer potential for better monitoring and management of highly interconnected economic and financial systems and, thus, may help anticipate and manage future crises.
A very general cosmological consideration suggests that, along with galactic dark matter halos, much smaller dark matter structures may exist. These structures are usually called `clumps', and their mass extends to 10−6 M ⊙ or even lower. The clumps should give the main contribution into the signal of dark matter annihilation, provided that they have survived until the present time. Recent observations favor a cored profile for low-mass astrophysical halos. We consider cored clumps and show that they are significantly less firm than the standard NFW ones. In contrast to the standard scenario, the cored clumps should have been completely destroyed inside ~ 20 kpc from the Milky Way center. The dwarf spheroidals should not contain any dark matter clumps. On the other hand, even under the most pessimistic assumption about the clump structure, the clumps should have survived in the Milky Way at a distance exceeding 50 kpc from the center, as well as in low-density cosmic structures. There they significantly boost the dark matter annihilation. We show that at least 70% of the clumps endured the primordial structure formation should still exist untouched in the present-day Universe.
We do magnetohydrodynamic (MHD) simulations of local box models of turbulent Interstellar Medium (ISM) and analyse the process of amplification and saturation of mean magnetic fields with methods of mean field dynamo theory. It is shown that the process of saturation of mean fields can be partially described by the prolonged diffusion time scales in presence of the dynamically significant magnetic fields. However, the outward wind also plays an essential role in the saturation in higher SN rate case. Algebraic expressions for the back reaction of the magnetic field onto the turbulent transport coefficients are derived, which allow a complete description of the nonlinear dynamo. We also present the effects of dynamically significant mean fields on the ISM configuration and pressure distribution. We further add the cosmic ray component in the simulations and investigate the kinematic growth of mean fields with a dynamo perspective.
In this Comment, we review the results of pattern formation in a reaction-diffusion-advection system following the kinetics of the Gray-Scott model. A recent paper by Das [Phys. Rev. E 92, 052914 (2015)] shows that spatiotemporal chaos of the intermittency type can disappear as the advective flow is increased. This study, however, refers to a single point in the space of kinetic parameters of the original Gray-Scott model. Here we show that the wealth of patterns increases substantially as some of these parameters are changed. In addition to spatiotemporal intermittency, defect-mediated turbulence can also be found. In all cases, however, the chaotic behavior is seen to disappear as the advective flow is increased, following a scenario similar to what was reported in our earlier work [I. Berenstein and C. Beta, Phys. Rev. E 86, 056205 (2012)] as well as by Das. We also point out that a similar phenomenon can be found in other reaction-diffusion-advection models, such as the Oregonator model for the Belousov-Zhabotinsky reaction under flow conditions.
To turn or not to turn?
(2016)
Bacteria typically swim in straight runs, interruped by sudden turning events. In particular, some species are limited to a reversal in the swimming direction as the only turning maneuver at their disposal. In a recent article, Grossmann et al (2016 New J. Phys. 18 043009) introduce a theoretical framework to analyze the diffusive properties of active particles following this type of run-and-reverse pattern. Based on a stochastic clock model to mimic the regulatory pathway that triggers reversal events, they show that a run-and-reverse swimmer can optimize its diffusive spreading by tuning the reversal rate according to the level of rotational noise. With their approach, they open up promising new perspectives of how to incorporate the dynamics of intracellular signaling into coarse-grained active particle descriptions.
The concept of bonding and antibonding orbitals is fundamental in chemistry. The population of those orbitals and the energetic difference between the two reflect the strength of the bonding interaction. Weakening the bond is expected to reduce this energetic splitting, but the transient character of bond-activation has so far prohibited direct experimental access. Here we apply time-resolved soft X-ray spectroscopy at a free electron laser to directly observe the decreased bonding antibonding splitting following bond-activation using an ultrashort optical laser pulse.
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion
(2016)
It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases.
This publications-based thesis summarizes my contribution to the scientific field of ultrafast structural dynamics. It consists of 16 publications, about the generation, detection and coupling of coherent gigahertz longitudinal acoustic phonons, also called hypersonic waves. To generate such high frequency phonons, femtosecond near infrared laser pulses were used to heat nanostructures composed of perovskite oxides on an ultrashort timescale. As a consequence the heated regions of such a nanostructure expand and a high frequency acoustic phonon pulse is generated. To detect such coherent acoustic sound pulses I use ultrafast variants of optical Brillouin and x-ray scattering. Here an incident optical or x-ray photon is scattered by the excited sound wave in the sample. The scattered light intensity measures the occupation of the phonon modes.
The central part of this work is the investigation of coherent high amplitude phonon wave packets which can behave nonlinearly, quite similar to shallow water waves which show a steepening of wave fronts or solitons well known as tsunamis. Due to the high amplitude of the acoustic wave packets in the solid, the acoustic properties can change significantly in the vicinity of the sound pulse. This may lead to a shape change of the pulse. I have observed by time-resolved Brillouin scattering, that a single cycle hypersound pulse shows a wavefront steepening. I excited hypersound pulses with strain amplitudes until 1% which I have calibrated by ultrafast x-ray diffraction (UXRD).
On the basis of this first experiment we developed the idea of the nonlinear mixing of narrowband phonon wave packets which we call "nonlinear phononics" in analogy with the nonlinear optics, which summarizes a kaleidoscope of surprising optical phenomena showing up at very high electric fields. Such phenomena are for instance Second Harmonic Generation, four-wave-mixing or solitons. But in case of excited coherent phonons the wave packets have usually very broad spectra which make it nearly impossible to look at elementary scattering processes between phonons with certain momentum and energy.
For that purpose I tested different techniques to excite narrowband phonon wave packets which mainly consist of phonons with a certain momentum and frequency. To this end epitaxially grown metal films on a dielectric substrate were excited with a train of laser pulses. These excitation pulses drive the metal film to oscillate with the frequency given by their inverse temporal displacement and send a hypersonic wave of this frequency into the substrate. The monochromaticity of these wave packets was proven by ultrafast optical Brillouin and x-ray scattering.
Using the excitation of such narrowband phonon wave packets I was able to observe the Second Harmonic Generation (SHG) of coherent phonons as a first example of nonlinear wave mixing of nanometric phonon wave packets.
We consider two coupled populations of leaky integrate-and-fire neurons. Depending on the coupling strength, mean fields generated by these populations can have incommensurate frequencies or become frequency locked. In the observed 2:1 locking state of the mean fields, individual neurons in one population are asynchronous with the mean fields, while in another population they have the same frequency as the mean field. These synchronous neurons form a chimera state, where part of them build a fully synchronized cluster, while other remain scattered. We explain this chimera as a marginal one, caused by a self-organized neutral dynamics of the effective circle map.
We use a background quasar to detect the presence of circumgalactic gas around a z = 0.91 low-mass star-forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope show that the galaxy has a dust-corrected star formation rate (SFR) of 4.7 +/- 2.0. M-circle dot yr(-1), with no companion down to 0.22 M-circle dot yr(-1) (5 sigma) within 240 h(-1) kpc ("30"). Using a high-resolution spectrum of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle alpha of only 15 degrees), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold-flow disk" extending at least 12 kpc (3 x R-1/2). We estimate the mass accretion rate M-in to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the H (I) column density of log N-H (I)/cm(-2) similar or equal to 20.4 obtained from a Hubble Space Telescope/COS near-UV spectrum. From a detailed analysis of the low-ionization lines (e.g., Zn II, Cr II, Ti II, MnII, Si II), the accreting material appears to be enriched to about 0.4 Z(circle dot) (albeit with large uncertainties: log Z/Z(circle dot) = -0.4 +/- 0.4), which is comparable to the galaxy metallicity (12 + log O/H = 8.7 +/- 0.2), implying a large recycling fraction from past outflows. Blueshifted Mg II and Fe II absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The Mg II and Fe II absorption line ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent Fe II* emission.
We present a statistical analysis of phase space density data from the first 26 months of the Van Allen Probes mission. In particular, we investigate the relationship between the tens and hundreds of keV seed electrons and >1 MeV core radiation belt electron population. Using a cross-correlation analysis, we find that the seed and core populations are well correlated with a coefficient of approximate to 0.73 with a time lag of 10-15 h. We present evidence of a seed population threshold that is necessary for subsequent acceleration. The depth of penetration of the seed population determines the inner boundary of the acceleration process. However, we show that an enhanced seed population alone is not enough to produce acceleration in the higher energies, implying that the seed population of hundreds of keV electrons is only one of several conditions required for MeV electron radiation belt acceleration.
In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting approximately 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high, we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detections of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousand seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.
The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the nonstationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total timescale of several hours), the transitions of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the nonstationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favored in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.
The Milky Way is only one out of billions of galaxies in the universe. However, it is a special galaxy because it allows to explore the main mechanisms involved in its evolution and formation history by unpicking the system star-by-star. Especially, the chemical fingerprints of its stars provide clues and evidence of past events in the Galaxy’s lifetime. These information help not only to decipher the current structure and building blocks of the Milky Way, but to learn more about the general formation process of galaxies.
In the past decade a multitude of stellar spectroscopic Galactic surveys have scanned millions of stars far beyond the rim of the solar neighbourhood. The obtained spectroscopic information provide unprecedented insights to the chemo-dynamics of the Milky Way. In addition analytic models and numerical simulations of the Milky Way provide necessary descriptions and predictions suited for comparison with observations in order to decode the physical properties that underlie the complex system of the Galaxy.
In the thesis various approaches are taken to connect modern theoretical modelling of galaxy formation and evolution with observations from Galactic stellar surveys. With its focus on the chemo-kinematics of the Galactic disk this work aims to determine new observational constraints on the formation of the Milky Way providing also proper comparisons with two different models. These are the population synthesis model TRILEGAL based on analytical distribution functions, which aims to simulate the number and distribution of stars in the Milky Way and its different components, and a hybrid model (MCM) that combines an N-body simulation of a Milky Way like galaxy in the cosmological framework with a semi-analytic chemical evolution model for the Milky Way. The major observational data sets in use come from two surveys, namely the “Radial Velocity Experiment” (RAVE) and the “Sloan Extension for Galactic Understanding and Exploration” (SEGUE).
In the first approach the chemo-kinematic properties of the thin and thick disk of the Galaxy as traced by a selection of about 20000 SEGUE G-dwarf stars are directly compared to the predictions by the MCM model. As a necessary condition for this, SEGUE's selection function and its survey volume are evaluated in detail to correct the spectroscopic observations for their survey specific selection biases. Also, based on a Bayesian method spectro-photometric distances with uncertainties below 15% are computed for the selection of SEGUE G-dwarfs that are studied up to a distance of 3 kpc from the Sun.
For the second approach two synthetic versions of the SEGUE survey are generated based on the above models. The obtained synthetic stellar catalogues are then used to create mock samples best resembling the compiled sample of observed SEGUE G-dwarfs. Generally, mock samples are not only ideal to compare predictions from various models. They also allow validation of the models' quality and improvement as with this work could be especially achieved for TRILEGAL. While TRILEGAL reproduces the statistical properties of the thin and thick disk as seen in the observations, the MCM model has shown to be more suitable in reproducing many chemo-kinematic correlations as revealed by the SEGUE stars. However, evidence has been found that the MCM model may be missing a stellar component with the properties of the thick disk that the observations clearly show. While the SEGUE stars do indicate a thin-thick dichotomy of the stellar Galactic disk in agreement with other spectroscopic stellar studies, no sign for a distinct metal-poor disk is seen in the MCM model.
Usually stellar spectroscopic surveys are limited to a certain volume around the Sun covering different regions of the Galaxy’s disk. This often prevents to obtain a global view on the chemo-dynamics of the Galactic disk. Hence, a suitable combination of stellar samples from independent surveys is not only useful for the verification of results but it also helps to complete the picture of the Milky Way. Therefore, the thesis closes with a comparison of the SEGUE G-dwarfs and a sample of RAVE giants. The comparison reveals that the chemo-kinematic relations agree in disk regions where the samples of both surveys show a similar number of stars. For those parts of the survey volumes where one of the surveys lacks statistics they beautifully complement each other. This demonstrates that the comparison of theoretical models on the one side, and the combined observational data gathered by multiple surveys on the other side, are key ingredients to understand and disentangle the structure and formation history of the Milky Way.
This work reports about new high-resolution imaging and spectroscopic observations of solar type III radio bursts at low radio frequencies in the range from 30 to 80 MHz. Solar type III radio bursts are understood as result of the beam-plasma interaction of electron beams in the corona. The Sun provides a unique opportunity to study these plasma processes of an active star. Its activity appears in eruptive events like flares, coronal mass ejections and radio bursts which are all accompanied by enhanced radio emission. Therefore solar radio emission carries important information about plasma processes associated with the Sun’s activity. Moreover, the Sun’s atmosphere is a unique plasma laboratory with plasma processes under conditions not found in terrestrial laboratories. Because of the Sun’s proximity to Earth, it can be studied in greater detail than any other star but new knowledge about the Sun can be transfer to them. This “solar stellar connection” is important for the understanding of processes on other stars.
The novel radio interferometer LOFAR provides imaging and spectroscopic capabilities to study these processes at low frequencies. Here it was used for solar observations.
LOFAR, the characteristics of its solar data and the processing and analysis of the latter with the Solar Imaging Pipeline and Solar Data Center are described. The Solar Imaging Pipeline is the central software that allows using LOFAR for solar observations. So its development was necessary for the analysis of solar LOFAR data and realized here. Moreover a new density model with heat conduction and Alfvén waves was developed that provides the distance of radio bursts to the Sun from dynamic radio spectra.
Its application to the dynamic spectrum of a type III burst observed on March 16, 2016 by LOFAR shows a nonuniform radial propagation velocity of the radio emission. The analysis of an imaging observation of type III bursts on June 23, 2012 resolves a burst as bright, compact region localized in the corona propagating in radial direction along magnetic field lines with an average velocity of 0.23c. A nonuniform propagation velocity is revealed. A new beam model is presented that explains the nonuniform motion of the radio source as a propagation effect of an electron ensemble with a spread velocity distribution and rules out a monoenergetic electron distribution. The coronal electron number density is derived in the region from 1.5 to 2.5 R☉ and fitted with the newly developed density model. It determines the plasma density for the interplanetary space between Sun and Earth. The values correspond to a 1.25- and 5-fold Newkirk model for harmonic and fundamental emission, respectively. In comparison to data from other radio instruments the LOFAR data shows a high sensitivity and resolution in space, time and frequency.
The new results from LOFAR’s high resolution imaging spectroscopy are consistent with current theories of solar type III radio bursts and demonstrate its capability to track fast moving radio sources in the corona. LOFAR solar data is found to be a valuable source for solar radio physics and opens a new window for studying plasma processes associated with highly energetic electrons in the solar corona.
Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback.
Adsorption of amino acids on the magnetite-(111)-surface: a force field study (vol 19, 851, 2013)
(2016)
Particle diffusion and localized acceleration in inhomogeneous AGN jets - II. Stochastic variation
(2016)
We study the stochastic variation of blazar emission under a 2D spatially resolved leptonic jet model we previously developed. Random events of particle acceleration and injection in small zones within the emission region are assumed to be responsible for flux variations. In addition to producing spectral energy distributions that describe the observed flux of Mrk 421, we further analyse the timing properties of the simulated light curves, such as the power spectral density (PSD) at different bands, flux-flux correlations, aswell as the cross-correlation function between X-rays and TeV gamma-rays. We find spectral breaks in the PSD at a time-scale comparable to the dominant characteristic time-scale in the system, which is usually the predefined decay time-scale of an acceleration event. Cooling imposes a delay, and so PSDs taken at lower energy bands in each emission component (synchrotron or inverse Compton) generally break at longer time-scales. The flux-flux correlation between X-rays and TeV gamma-rays can be either quadratic or linear, depending on whether or not there are large variation of the injection into the particle acceleration process. When the relationship is quadratic, the TeV flares lag the X-ray flares, and the optical and GeV flares are large enough to be comparable to the ones in X-ray. When the relationship is linear, the lags are insignificant, and the optical and GeV flares are small.
We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.
We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.
We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion.
XVII. Proper motions of the Small Magellanic Cloud and the Milky Way globular cluster 47 Tucanae
(2016)
Aims. In this study we use multi-epoch near-infrared observations from the VISTA survey of the Magellanic Cloud system (VMC) to measure the proper motions of different stellar populations in a tile of 1.5 deg2 in size in the direction of the Galactic globular cluster 47 Tuc. We obtain the proper motion of the cluster itself, of the Small Magellanic Cloud (SMC), and of the field Milky Way stars.
Methods. Stars of the three main stellar components are selected according to their spatial distributions and their distributions in colour−magnitude diagrams. Their average coordinate displacement is computed from the difference between multiple Ks-band observations for stars as faint as Ks = 19 mag. Proper motions are derived from the slope of the best-fitting line among ten VMC epochs over a time baseline of ~1 yr. Background galaxies are used to calibrate the absolute astrometric reference frame.
Results. The resulting absolute proper motion of 47 Tuc is (μαcos(δ), μδ) = (+7.26 ± 0.03, −1.25 ± 0.03) mas yr-1. This measurement refers to about 35 000 sources distributed between 10′ and 60′ from the cluster centre. For the SMC we obtain (μαcos(δ), μδ) = (+1.16 ± 0.07, −0.81 ± 0.07) mas yr-1 from about 5250 red clump and red giant branch stars. The absolute proper motion of the Milky Way population in the line of sight (l = 305.9, b = −44.9) of this VISTA tile is (μαcos(δ), μδ) = (+10.22 ± 0.14, −1.27 ± 0.12) mas yr-1 and has been calculated from about 4000 sources. Systematic uncertainties associated with the astrometric reference system are 0.18 mas yr-1. Thanks to the proper motion we detect 47 Tuc stars beyond its tidal radius.
Consequences of twenty-first-century policy for multi-millennial climate
and sea-level change
(2016)
Most of the policy debate surrounding the actions needed to mitigate and adapt to anthropogenic climate change has been framed by observations of the past 150 years as well as climate and sea-level projections for the twenty-first century. The focus on this 250-year window, however, obscures some of the most profound problems associated with climate change. Here, we argue that the twentieth and twenty-first centuries, a period during which the overwhelming majority of human-caused carbon emissions are likely to occur, need to be placed into a long-term context that includes the past 20 millennia, when the last Ice Age ended and human civilization developed, and the next ten millennia, over which time the projected impacts of anthropogenic climate change will grow and persist. This long-term perspective illustrates that policy decisions made in the next few years to decades will have profound impacts on global climate, ecosystems and human societies - not just for this century, but for the next ten millennia and beyond.
We show that self-consistent partial synchrony in globally coupled oscillatory ensembles is a general phenomenon. We analyze in detail appearance and stability properties of this state in possibly the simplest setup of a biharmonic Kuramoto-Daido phase model as well as demonstrate the effect in limit-cycle relaxational Rayleigh oscillators. Such a regime extends the notion of splay state from a uniform distribution of phases to an oscillating one. Suitable collective observables such as the Kuramoto order parameter allow detecting the presence of an inhomogeneous distribution. The characteristic and most peculiar property of self-consistent partial synchrony is the difference between the frequency of single units and that of the macroscopic field.
Thermophony in real gases
(2016)
A thermophone is an electrical device for sound generation. The advantages of thermophones over conventional sound transducers such as electromagnetic, electrostatic or piezoelectric transducers are their operational principle which does not require any moving parts, their resonance-free behavior, their simple construction and their low production costs.
In this PhD thesis, a novel theoretical model of thermophonic sound generation in real gases has been developed. The model is experimentally validated in a frequency range from 2 kHz to 1 MHz by testing more then fifty thermophones of different materials, including Carbon nano-wires, Titanium, Indium-Tin-Oxide, different sizes and shapes for sound generation in gases such as air, argon, helium, oxygen, nitrogen and sulfur hexafluoride.
Unlike previous approaches, the presented model can be applied to different kinds of thermophones and various gases, taking into account the thermodynamic properties of thermophone materials and of adjacent gases, degrees of freedom and the volume occupied by the gas atoms and molecules, as well as sound attenuation effects, the shape and size of the thermophone surface and the reduction of the generated acoustic power due to photonic emission. As a result, the model features better prediction accuracy than the existing models by a factor up to 100. Moreover, the new model explains previous experimental findings on thermophones which can not be explained with the existing models.
The acoustic properties of the thermophones have been tested in several gases using unique, highly precise experimental setups comprising a Laser-Doppler-Vibrometer combined with a thin polyethylene film which acts as a broadband and resonance-free sound-pressure detector. Several outstanding properties of the thermophones have been demonstrated for the first time, including the ability to generate arbitrarily shaped acoustic signals, a greater acoustic efficiency compared to conventional piezoelectric and electrostatic airborne ultrasound transducers, and applicability as powerful and tunable sound sources with a bandwidth up to the megahertz range and beyond.
Additionally, new applications of thermophones such as the study of physical properties of gases, the thermo-acoustic gas spectroscopy, broad-band characterization of transfer functions of sound and ultrasound detection systems, and applications in non-destructive materials testing are discussed and experimentally demonstrated.
We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems.
We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems.
We study the adsorption-desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces-are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye-Huckel approximation is often not feasible and the nonlinear Poisson-Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson-Boltzmann equation is smaller than the Debye-Huckel result, such that the required critical surface charge density for polyelectrolyte adsorption sigma(c) increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical-chemical and biophysical systems.
Context. Most theoretical investigations of dust charging processes in space have treated the current balance condition as independent of grain size. However, for small grains, since they are often observed in space environments, a dependence on grain size is expected owing to secondary electron emission (SEE). Here, by the term "small" we mean a particle size comparable to the typical penetration depth for given primary electron energy. The results are relevant for the dynamics of small, charged dust particles emitted by the volcanic moon Io, which forms the Jovian dust streams. Aims. We revise the theory of charging of small (submicron sized) micrometeoroids to take into account a high production of secondary electrons for small grains immersed in an isotropic flux of electrons. We apply our model to obtain an improved estimate for the charge of the dust streams leaving the Jovian system, detected by several spacecraft. Methods. We apply a continuum model to describe the penetration of primary electrons in a grain and the emission of secondary electrons along the path. Averaging over an isotropic flux of primaries, we derive a new expression for the secondary electron yield, which can be used to express the secondary electron current on a grain. Results. For the Jupiter plasma environment we derive the surface potential of grains composed of NaCl (believed to be the major constituent of Jovian dust stream particles) or silicates. For small particles, the potential depends on grain size and the secondary electron current induces a sensitivity to material properties. As a result of the small particle effect, the estimates for the charging times and for the fractional charge fluctuations of NaCl grains obtained using our general approach to SEE give results qualitatively different from the analogous estimates derived from the traditional approach to SEE. We find that for the charging environment considered in this paper field emission does not limit the charging of NaCl grains.
The speciation of 2-Mercaptopyridine in aqueous solution has been investigated with nitrogen 1s Near Edge X-ray Absorption Fine Structure spectroscopy and time dependent Density Functional Theory. The prevalence of distinct species as a function of the solvent basicity is established. No indications of dimerization towards high concentrations are found. The determination of different molecular structures of 2-Mercaptopyridine in aqueous solution is put into the context of proton-transfer in keto-enol and thione-thiol tautomerisms. (C) 2016 The Authors. Published by Elsevier B.V.