Refine
Year of publication
- 2015 (295) (remove)
Document Type
- Article (250)
- Doctoral Thesis (28)
- Postprint (7)
- Other (3)
- Conference Proceeding (2)
- Review (2)
- Habilitation Thesis (1)
- Master's Thesis (1)
- Preprint (1)
Keywords
- anomalous diffusion (8)
- stars: early-type (8)
- Magellanic Clouds (7)
- stars: massive (7)
- stars: atmospheres (6)
- stars: mass-loss (6)
- stars: winds, outflows (6)
- ISM: supernova remnants (5)
- X-rays: binaries (5)
- X-rays: stars (5)
Institute
- Institut für Physik und Astronomie (295) (remove)
We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T−h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.
The electrical conductivity of organic semiconductors can be enhanced by orders of magnitude via doping with strong molecular electron acceptors or donors. Ground-state integer charge transfer and charge-transfer complex formation between organic semiconductors and molecular dopants have been suggested as the microscopic mechanisms causing these profound changes in electrical materials properties. Here, we study charge-transfer interactions between the common molecular p-dopant 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane and a systematic series of thiophene-based copolymers by a combination of spectroscopic techniques and electrical measurements. Subtle variations in chemical structure are seen to significantly impact the nature of the charge-transfer species and the efficiency of the doping process, underlining the need for a more detailed understanding of the microscopic doping mechanism in organic semiconductors to reliably guide targeted chemical design.
We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer-obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer-obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer-crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids.
The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.
Mixing layer manipulation experiment from open-loop forcing to closed-loop machine learning control
(2015)
Context. Blazars are variable sources on various timescales over a broad energy range spanning from radio to very high energy (>100 GeV, hereafter VHE). Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. However, most of the gamma-ray studies performed on Mrk 501 during the past years relate to flaring activity, when the source detection and characterization with the available gamma-ray instrumentation was easier to perform.
Aims. Our goal is to characterize the source gamma-ray emission in detail, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity.
Methods. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. This provided extensive energy and temporal coverage of Mrk 501 throughout the entire campaign.
Results. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy and a tentative correlation between the X-ray and VHE fluxes. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density.
Conclusions. The one-zone SSC model can adequately describe the broadband spectral energy distribution of the source during the two months covered by the MW campaign. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during such a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.
The gamma-ray spectrum of the low-frequency-peaked BL Lac (LBL) object AP Librae is studied, following the discovery of very-high-energy (VHE; E > 100 GeV) gamma-ray emission up to the TeV range by the H.E.S.S. experiment. Thismakes AP Librae one of the few VHE emitters of the LBL type. The measured spectrum yields a flux of (8.8 +/- 1.5(stat) +/- 1.8(sys)) x 10(-12) cm(-2) s(-1) above 130 GeV and a spectral index of Gamma = 2.65 +/- 0.19(stat) +/- 0.20(sys). This study also makes use of Fermi-LAT observations in the high energy (HE, E > 100 MeV) range, providing the longest continuous light curve (5 years) ever published on this source. The source underwent a flaring event between MJD 56 306-56 376 in the HE range, with a flux increase of a factor of 3.5 in the 14 day bin light curve and no significant variation in spectral shape with respect to the low-flux state. While the H.E.S.S. and (low state) Fermi-LAT fluxes are in good agreement where they overlap, a spectral curvature between the steep VHE spectrum and the Fermi-LAT spectrum is observed. The maximum of the gamma-ray emission in the spectral energy distribution is located below the GeV energy range.
We used microfluidic tools and high-speed time-lapse microscopy to record trajectories of the soil bacterium Pseudomonas putida in a confined environment with cells swimming in close proximity to a glass-liquid interface. While the general swimming pattern is preserved, when compared to swimming in the bulk fluid, our results show that cells in the presence of two solid boundaries display more frequent reversals in swimming direction and swim faster. Additionally, we observe that run segments are no longer straight and that cells swim on circular trajectories, which can be attributed to the hydrodynamic wall effect. Using the experimentally observed parameters together with a recently presented analytic model for a run-reverse random walker, we obtained additional insight on how the spreading behavior of a cell population is affected under confinement. While on short time scales, the mean square displacement of confined swimmers grows faster as compared to the bulk fluid case, our model predicts that for large times the situation reverses due to the strong increase in effective rotational diffusion.
The initiation of a marine ice-sheet instability (MISI) is generally discussed from the ocean side of the ice sheet. It has been shown that the reduction in ice-shelf buttressing and softening of the coastal ice can destabilize a marine ice sheet if the bedrock is sloping upward towards the ocean. Using a conceptional flow-line geometry, we investigate the possibility of whether a MISI can be triggered from the direction of the ice divide as opposed to coastal forcing and explore the interaction between connected basins. We find that the initiation of a MISI in one basin can induce a destabilization in the other. The underlying mechanism of basin interaction is based on dynamic thinning and a consecutive motion of the ice divide which induces a thinning in the adjacent basin and a successive initiation of the instability. Our simplified and symmetric topographic setup allows scaling both the geometry and the transition time between both instabilities. We find that the ice profile follows a universal shape that is scaled with the horizontal extent of the ice sheet and that the same exponent of 1/2 applies for the scaling relation between central surface elevation and horizontal extent as in the pure shallow ice approximation (Vialov profile). Altering the central bed elevation, we find that the extent of grounding-line retreat in one basin determines the degree of interaction with the other. Different scenarios of basin interaction are discussed based on our modeling results as well as on a conceptual flux-balance analysis. We conclude that for the three-dimensional case, the possibility of drainage basin interaction on timescales on the order of 1 kyr or larger cannot be excluded and hence needs further investigation.
Paper-based microfluidics provide an inexpensive, easy to use technology for point-of-care diagnostics in developing countries. Here, we combine paper-based microfluidic devices with responsive hydrogels to add an entire new class of functions to these versatile low-cost fluidic systems. The hydrogels serve as fluid reservoirs. In response to an external stimulus, e.g. an increase in temperature, the hydrogels collapse and release fluid into the structured paper substrate. In this way, chemicals that are either stored on the paper substrate or inside the hydrogel pads can be dissolved, premixed, and brought to reaction to fulfill specific analytic tasks. We demonstrate that multi-step sequences of chemical reactions can be implemented in a paper-based system and operated without the need for external precision pumps. We exemplify this technology by integrating an antibody-based E. coli test on a small and easy to use paper device.
What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer-surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density sigma(c)-defining the adsorption-desorption transition-are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of sigma(c) for the concave interfaces versus the Debye screening length 1/kappa and the extent of confinement a for these three interfaces for small kappa a values. For large kappa a the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/kappa. We also rationalize how sigma(c)(kappa) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids.
Brownian motion is ergodic in the Boltzmann-Khinchin sense that long time averages of physical observables such as the mean squared displacement provide the same information as the corresponding ensemble average, even at out-of-equilibrium conditions. This property is the fundamental prerequisite for single particle tracking and its analysis in simple liquids. We study analytically and by event-driven molecular dynamics simulations the dynamics of force-free cooling granular gases and reveal a violation of ergodicity in this Boltzmann-Khinchin sense as well as distinct ageing of the system. Such granular gases comprise materials such as dilute gases of stones, sand, various types of powders, or large molecules, and their mixtures are ubiquitous in Nature and technology, in particular in Space. We treat-depending on the physical-chemical properties of the inter-particle interaction upon their pair collisions-both a constant and a velocity-dependent (viscoelastic) restitution coefficient epsilon. Moreover we compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behaviour of the ensemble mean squared displacement (MSD) and the velocity correlations in the limit of weak dissipation. Qualitatively, the reported non-ergodic behaviour is generic for granular gases with any realistic dependence of epsilon on the impact velocity of particles.
X-ray spectroscopy is a powerful tool to study the local charge distribution of chemical systems. Together with the liquid jet it becomes possible to probe chemical systems in their natural environment, the liquid phase. In this work, we present X-ray absorption (XA), X-ray emission (XE) and resonant inelastic X-ray scattering (RIXS) data of pure water and various salt solutions and show the possibilities these methods offer to elucidate the nature of ion-water interaction.
We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law similar to T-h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.
The plasmon resonance of metal nanoparticles determines their optical response in the visible spectral range. Many details such as the electronic properties of gold near the particle surface and the local environment of the particles influence the spectra. We show how the cheap but highly precise fabrication of composite nanolayers by spin-assisted layer-by-layer deposition of polyelectrolytes can be used to investigate the spectral response of gold nanospheres (GNS) and gold nanorods (GNR) in a self-consistent way, using the established Maxwell-Garnett effective medium (MGEM) theory beyond the limit of homogeneous media. We show that the dielectric function of gold nanoparticles differs from the bulk value and experimentally characterize the shape and the surrounding of the particles thoroughly by SEM, AFM and ellipsometry. Averaging the dielectric functions of the layered surrounding by an appropriate weighting with the electric field intensity yields excellent agreement for the spectra of several nanoparticles and nanorods with various cover-layer thicknesses.
The electrophoretic deposition process was used to produce multi-layered ceramics consisting of alternating layers of fully stabilized cubic zirconia and partially stabilized tetragonal zirconia to make use of their different mechanical behaviour, investigating the possibility to deflect advancing cracks at the interfaces of the different layers. This crack deflection is apparently impacted by a toughening mechanism only found in the tetragonal stabilized zirconia polymorph and is characterized by the stress induced transformation of the metastable tetragonal phase into the monoclinic one, which is accompanied by a volume increase resulting in a closing mechanism for advancing cracks.
While improving the electrophoretic deposition process, we investigated the transformation toughening mechanism at the layer interfaces and their effect on crack propagation. Investigations involved a combination of different imaging methods, including light microscopy, white light interferometry, atomic force microscopy, scanning electron microscopy, energy dispersive X-ray spectroscopy and Raman spectroscopy.
We present the XMM-Newton discovery of X-ray emission from the planetary nebula (PN) A78, the second born-again PN detected in X-rays apart from A30. These two PNe share similar spectral and morphological characteristics: they harbor diffuse soft X-ray emission associated with the interaction between the H-poor ejecta and the current fast stellar wind and a point-like source at the position of the central star (CSPN). We present the spectral analysis of the CSPN, using for the first time an NLTE code for expanding atmospheres that takes line blanketing into account for the UV and optical spectra. The wind abundances are used for the X-ray spectral analysis of the CSPN and the diffuse emission. The X-ray emission from the CSPN in A78 can be modeled by a single C VI emission line, while the X-ray emission from its diffuse component is better described by an optically thin plasma emission model with a temperature of kT = 0.088 keV (T approximate to 1.0 x 10(6) K). We estimate X-ray luminosities in the 0.2-2.0 keV energy band of L-X,L-CSPN =(1.2 +/- 0.3) x 10(31) erg s(-1) and L-X,L-DIFF =(9.2 +/- 2.3) x 10(30) erg s(-1) for the CSPN and diffuse components, respectively.
We present the results of a broadband X-ray study of the enigmatic Be star Gamma Cassiopeiae (herein gamma Cas) based on observations made with both the Suzaku and INTEGRAL observatories.. Cas has long been recognized as the prototypical example of a small subclass of Be stars with moderately strong X-ray emission dominated by a hot thermal component in the 0.5-12 keV energy range (L-x approximate to 10(32)-10(33) erg s(-1)). This places them at the high end of the known luminosity distribution for stellar emission, but several orders of magnitude below typical accretion-powered Be X-ray binaries. The INTEGRAL observations spanned an eight-year baseline and represent the deepest measurement to date at energies above similar to 50 keV. We find that the INTEGRAL data are consistent within statistics to a constant intensity source above 20 keV, with emission extending up to similar to 100 keV, and that searches for all of the previously reported periodicities of the system at lower energies led to null results. We further find that our combined Suzaku and INTEGRAL spectrum, which we suggest is the most accurate broadband X-ray measurement of gamma Cas to date, is fitted extremely well with a thermal plasma emission model with a single absorption component. We found no compelling need for an additional non-thermal high-energy component. We discuss these results in the context of a currently favored models for gamma Cas and its analogs.
Magnetic reconnection in the partially ionized solar chromosphere is studied in 2.5 dimensional magnetohydrodynamic simulations including radiative cooling and ambipolar diffusion. A Harris current sheet with and without a guide field is considered. Characteristic values of the parameters in the middle chromosphere imply a high magnetic Reynolds number of similar to 10(6)-10(7) in the present simulations. Fast magnetic reconnection then develops as a consequence of the plasmoid instability without the need to invoke anomalous resistivity enhancements. Multiple levels of the instability are followed as it cascades to smaller scales, which approach the ion inertial length. The reconnection rate, normalized to the asymptotic values of magnetic field and Alfven velocity in the inflow region, reaches values in the range similar to 0.01-0.03 throughout the cascading plasmoid formation and for zero as well as for strong guide field. The outflow velocity reaches approximate to 40 km s(-1). Slow-mode shocks extend from the X-points, heating the plasmoids up to similar to 8 x 10(4) K. In the case of zero guide field, the inclusion of both ambipolar diffusion and radiative cooling causes a rapid thinning of the current sheet (down to similar to 30 m) and early formation of secondary islands. Both of these processes have very little effect on the plasmoid instability for a strong guide field. The reconnection rates, temperature enhancements, and upward outflow velocities from the vertical current sheet correspond well to their characteristic values in chromospheric jets.
We present results from VERITAS observations of the BL Lac object PG 1553+113 spanning the years 2010, 2011, and 2012. The time-averaged spectrum, measured between 160 and 560 GeV, is well described by a power law with a spectral index of 4.33 +/- 0.09. The time-averaged integral flux above 200 GeV measured for this period was (1.69 +/- 0.06) x 10(-11) photons cm(-2) s(-1), corresponding to 6.9% of the Crab Nebula flux. We also present the combined gamma-ray spectrum from the Fermi Large Area Telescope and VERITAS covering an energy range from 100 MeV to 560 GeV. The data are well fit by a power law with an exponential cutoff at 101.9 +/- 3.2 GeV. The origin of the cutoff could be intrinsic to PG 1553+113 or be due to the gamma-ray opacity of our universe through pair production off the extragalactic background light (EBL). Given lower limits to the redshift of z > 0.395 based on optical/UV observations of PG 1553+113, the cutoff would be dominated by EBL absorption. Conversely, the small statistical uncertainties of the VERITAS energy spectrum have allowed us to provide a robust upper limit on the redshift of PG 1553+113 of z <= 0.62. A strongly elevated mean flux of (2.50 +/- 0.14) x10(-11) photons cm(-2) s(-1) (10.3% of the Crab Nebula flux) was observed during 2012, with the daily flux reaching as high as (4.44 +/- 0.71) x10(-11) photons cm(-2) s(-1) (18.3% of the Crab Nebula flux) on MJD 56048. The light curve measured during the 2012 observing season is marginally inconsistent with a steady flux, giving a chi(2) probability for a steady flux of 0.03%.
We report on multifrequency observations performed during 2012 December-2013 August of the first narrow-line Seyfert 1 galaxy detected in gamma-rays, PMN J0948+0022 (z = 0.5846). A y -ray flare was observed by the Large Area Telescope on board Fermi during 2012 December-2013 January, reaching a daily peak flux in the 0.1-100 GeV energy range of (155 31) x 10 8 ph cm(-2) S-1 on 2013 January 1, corresponding to an apparent isotropic luminosity of similar to 1.5 x 1048 erg s(-1). The y -ray flaring period triggered Swift and Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations in addition to radio and optical monitoring by Owens Valley Radio Observatory, Monitoring Of Jets in Active galactic nuclei with VLBA Experiments, and Catalina Real-time Transient Survey. A strong flare was observed in optical, UV, and X-rays on 2012 December 30, quasi-simultaneously to the y -ray flare, reaching a record flux for this source from optical to y gamma-rays. VERITAS observations at very high energy (E > 100 GeV) during 2013 January 6-17 resulted in an upper limit of F>0.2 Trev < 4.0 x 10(-12) ph cm(-2) s(-1). We compared the spectral energy distribution (SED) of the flaring state in 2013 January with that of an intermediate state observed in 2011. The two SEDs, modelled as synchrotron emission and an external Compton scattering of seed photons from a dust torus, can be modelled by changing both the electron distribution parameters and the magnetic field.
The Large Magellanic Cloud, a satellite galaxy of the Milky Way, has been observed with the High Energy Stereoscopic System (H.E.S.S.) above an energy of 100 billion electron volts for a deep exposure of 210 hours. Three sources of different types were detected: the pulsar wind nebula of the most energetic pulsar known, N 157B; the radio-loud supernova remnant N 132D; and the largest nonthermal x-ray shell, the superbubble 30 Dor C. The unique object SN 1987A is, unexpectedly, not detected, which constrains the theoretical framework of particle acceleration in very young supernova remnants. These detections reveal the most energetic tip of a g-ray source population in an external galaxy and provide via 30 Dor C the unambiguous detection of g-ray emission from a superbubble.
The structure of bulk liquid water was recently probed by x-ray scattering below the temperature limit of homogeneous nucleation (T-H) of similar to 232 K [J. A. Sellberg et al., Nature 510, 381-384 (2014)]. Here, we utilize a similar approach to study the structure of bulk liquid water below T-H using oxygen K-edge x-ray emission spectroscopy (XES). Based on previous XES experiments [T. Tokushima et al., Chem. Phys. Lett. 460, 387-400 (2008)] at higher temperatures, we expected the ratio of the 1b(1)' and 1b(1)" peaks associated with the lone-pair orbital in water to change strongly upon deep supercooling as the coordination of the hydrogen (H-) bonds becomes tetrahedral. In contrast, we observed only minor changes in the lone-pair spectral region, challenging an interpretation in terms of two interconverting species. A number of alternative hypotheses to explain the results are put forward and discussed. Although the spectra can be explained by various contributions from these hypotheses, we here emphasize the interpretation that the line shape of each component changes dramatically when approaching lower temperatures, where, in particular, the peak assigned to the proposed disordered component would become more symmetrical as vibrational interference becomes more important. (C) 2015 AIP Publishing LLC.
We study the linear and nonlinear acoustic response of SrTiO3 across its ferroelastic transition at T-a = 105 K by time domain Brillouin scattering. Above T-a we observe that for a strain amplitude of similar to 0.18% the sound velocity for compressive strain exceeds the tensile strain velocity by 3%. Below T-a we find a giant slowing down of the sound velocity by 12% and attribute this to the coupling of GHz phonons to ferroelastic twin domain walls. We propose a new mechanism for this coupling on the ultrafast time scale, providing an important new test ground for theories used to simulate atomic motion in domain forming crystals.
DNA binding to trans- and cis-isomers of azobenzene containing cationic surfactant in 5 mM NaCl solution was investigated by the methods of dynamic light scattering (DLS), low-gradient viscometry (LGV), atomic force microscopy (AFM), circular dichroism (CD), gel electrophoresis (GE), flow birefringence (FB), UV-Vis spectrophotometry. Light-responsive conformational transitions of DNA in complex with photosensitive surfactant, changes in DNA optical anisotropy and persistent length, phase transition of DNA into nanoparticles induced by high surfactant concentration, as well as transformation of surfactant conformation under its binding to macromolecule were studied. Computer simulations of micelles formation for cis- and trans-isomers of azobenzene containing surfactant, as well as DNA-surfactant interaction, were carried out. Phase diagram for DNA-surfactant solutions was designed. The possibility to reverse the DNA packaging induced by surfactant binding with the dilution and light irradiation was shown. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 109-122, 2015.
Here we generalize our previous model of molecular motors trafficking subdiffusing cargos in viscoelastic cytosol by (i) including mechano-chemical coupling between cyclic conformational fluctuations of the motor protein driven by the reaction of ATP hydrolysis and its translational motion within the simplest two-state model of hand-over-hand motion of kinesin, and also (ii) by taking into account the anharmonicity of the tether between the motor and the cargo (its maximally possible extension length). It is shown that the major earlier results such as occurrence of normal versus anomalous transport depending on the amplitude of binding potential, cargo size and the motor turnover frequency not only survive in this more realistic model, but the results also look very similar for the correspondingly adjusted parameters. However, this more realistic model displays a substantially larger thermodynamic efficiency due to a bidirectional mechano-chemical coupling. For realistic parameters, the maximal thermodynamic efficiency can transiently be about 50% as observed for kinesins, and even larger, surprisingly also in a novel strongly anomalous (sub) transport regime, where the motor enzymatic turnovers become also anomalously slow and cannot be characterized by a turnover rate. Here anomalously slow dynamics of the cargo enforces anomalously slow cyclic kinetics of the motor protein.
Context. radio spectra of many shell-type supernova remnants show deviations from those expected on theoretical grounds.
Aims. In this paper we determine the effect of stochastic reacceleration on the spectra of electrons in the GeV band and at lower energies, and we investigate whether reacceleration can explain the observed variation in radio spectral indices.
Methods. We explicitely calculated the momentum diffusion coefficient for 3 types of turbulence expected downstream of the forward shock: fast-mode waves, small-scale non-resonant modes, and large-scale modes arising from turbulent dynamo activity. After noting that low-energy particles are efficiently coupled to the quasi-thermal plasma, a simplified cosmic-ray transport equation can be formulated and is numerically solved.
Results. Only fast-mode waves can provide momentum diffusion fast enough to significantly modify the spectra of particles. Using a synchrotron emissivity that accurately reflects a highly turbulent magnetic field, we calculated the radio spectral index and find that soft spectra with index a alpha less than or similar to -0.6 can be maintained over more than 2 decades in radio frequency, even if the electrons experience reacceleration for only one acceleration time. A spectral hardening is possible but considerably more frequency-dependent. The spectral modification imposed by stochastic reacceleration downstream of the forward shock depends only weakly on the initial spectrum provided by, e.g., diffusive shock acceleration at the shock itself.
Only a small fraction of massive stars seem to host a measurable structured magnetic field, whose origin is still unknown and whose implications for stellar evolution still need to be assessed. Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD 44743; B1 II/III) and epsilon CMa (HD 52089; B1.5II) in December 2013 and April 2014. For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field, approximately constant with time. We determined the physical parameters of both stars and characterise their X-ray spectrum. For the beta Cep star beta CMa, our mode identification analysis led to determining a rotation period of 13.6 +/- 1.2 days and of an inclination angle of the rotation axis of 57.6 +/- 1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of about 22 degrees and a dipolar magnetic field strength (B-d) of about 100 G (60 < B-d < 230 G within the 1 sigma level), below what is typically found for other magnetic massive stars. This conclusion is strengthened further by considerations of the star's X-ray spectrum. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Our results imply that both stars are expected to have a dynamical magnetosphere, so the magnetic field is not able to support a circumstellar disk. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.
G349.7+0.2 is a young Galactic supernova remnant (SNR) located at the distance of 11.5 kpc and observed across the entire electromagnetic spectrum from radio to high energy (HE; 0.1 GeV < E < 100 GeV) gamma-rays. Radio and infrared observations indicate that the remnant is interacting with a molecular cloud. In this paper, the detection of very high energy (VHE, E > 100 GeV) gamma-ray emission coincident with this SNR with the High Energy Stereoscopic System (HESS.) is reported. This makes it one of the farthest Galactic SNR ever detected in this domain. An integral flux F(E > 400 GeV) = (6.5 +/- 1.1(stat) +/- 1.3(syst)) x 10-11 ph cm(-2) s(-1) corresponding to similar to 0.7% of that of the Crab Nebula and to a luminosity of similar to 10(34) erg s(-1) above the same energy threshold, and a steep photon index Gamma(VHE) = 2.8 +/- 0.27(stat) +/- 0.20(syst) are measured. The analysis of more than 5 yr of Fermi-LAT data towards this source shows a power-law like spectrum with a best-fit photon index Gamma(HE) = 2.2 +/- 0.04.2(stat-0.31sys)(+0.13), The combined gamma-ray spectrum of 0349.7+0.2 can be described by either a broken power law (I3PL) or a power law with exponential (or sub exponential) cutoff (PLC). In the former case, the photon break energy is found at E-br,E-gamma = 551(-30)(+70) GeV, slightly higher than what is usually observed in the HE/VHE gamma-ray emitting middle-aged SNRs known to be interacting with molecular clouds. In the latter case. the exponential (respectively sub-exponential) cutoff energy is measured at E-cat,E-gamma = 1.4(-0.55)(+1.6) (respectively 0.35(-0.21)(+0.75)) TeV. A pion decay process resulting from the interaction of the accelerated protons and nuclei with the dense surrounding medium is clearly the preferred scenario to explain the gamma-ray emission. The BPL with a spectral steepening of 0.5-1 and the PLC provide equally good fits to the data. The product or the average gas density and the total energy content of accelerated protons and nuclei amounts to nu W-p similar to 5 x 10(51) erg cm(-3)
Aims. Previous observations with the High Energy Stereoscopic System (H.E.S.S.) have revealed an extended very-high-energy (VHE; E > 100 GeV) gamma-ray source, HESS J1834-087, coincident with the supernova remnant (SNR) W41. The origin of the gamma-ray emission was investigated in more detail with the H.E.S.S. array and the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope.
Methods. The gamma-ray data provided by 61 h of observations with H.E.S.S., and four years with the Fermi LAT were analyzed, covering over five decades in energy from 1.8 GeV up to 30 TeV. The morphology and spectrum of the TeV and GeV sources were studied and multiwavelength data were used to investigate the origin of the gamma-ray emission toward W41.
Results. The TeV source can be modeled with a sum of two components: one point-like and one significantly extended (sigma(TeV) = 0.17 degrees +/- 0.01 degrees), both centered on SNR W41 and exhibiting spectra described by a power law with index Gamma(TeV) similar or equal to 2.6. The GeV source detected with Fermi LAT is extended (sigma(GeV) = 0.15 degrees +/- 0.03 degrees) and morphologically matches the VHE emission. Its spectrum can be described by a power-law model with an index Gamma(GeV) = 2.15 +/- 0.12 and smoothly joins the spectrum of the whole TeV source. A break appears in the gamma-ray spectra around 100 GeV. No pulsations were found in the GeV range.
Conclusions. Two main scenarios are proposed to explain the observed emission: a pulsar wind nebula (PWN) or the interaction of SNR W41 with an associated molecular cloud. X-ray observations suggest the presence of a point-like source (a pulsar candidate) near the center of the remnant and nonthermal X-ray diffuse emission that could arise from the possibly associated PWN. The PWN scenario is supported by the compatible positions of the TeV and GeV sources with the putative pulsar. However, the spectral energy distribution from radio to gamma-rays is reproduced by a one-zone leptonic model only if an excess of low-energy electrons is injected following a Maxwellian distribution by a pulsar with a high spin-down power (> 10(37) erg s(-1)). This additional low-energy component is not needed if we consider that the point-like TeV source is unrelated to the extended GeV and TeV sources. The interacting SNR scenario is supported by the spatial coincidence between the gamma-ray sources, the detection of OH (1720 MHz) maser lines, and the hadronic modeling.
Most of our knowledge about the Sun's activity cycle arises from sunspot observations over the last centuries since telescopes have been used for astronomy. The German astronomer Gustav Sporer observed almost daily the Sun from 1861 until the beginning of 1894 and assembled a 33-year collection of sunspot data covering a total of 445 solar rotation periods. These sunspot drawings were carefully placed on an equidistant grid of heliographic longitude and latitude for each rotation period, which were then copied to copper plates for a lithographic reproduction of the drawings in astronomical journals. In this article, we describe in detail the process of capturing these data as digital images, correcting for various effects of the aging print materials, and preparing the data for contemporary scientific analysis based on advanced image processing techniques. With the processed data we create a butterfly diagram aggregating sunspot areas, and we present methods to measure the size of sunspots (umbra and penumbra) and to determine tilt angles of active regions. A probability density function of the sunspot area is computed, which conforms to contemporary data after rescaling. (C) 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim
Polymer looping is controlled by macromolecular crowding, spatial confinement, and chain stiffness
(2015)
We study by extensive computer simulations the looping characteristics of linear polymers with varying persistence length inside a spherical cavity in the presence of macromolecular crowding. For stiff chains, the looping probability and looping time reveal wildly oscillating patterns as functions of the chain length. The effects of crowding differ dramatically for flexible versus stiff polymers. While for flexible chains the looping kinetics is slowed down by the crowders, for stiffer chains the kinetics turns out to be either decreased or facilitated, depending on the polymer length. For severe confinement, the looping kinetics may become strongly facilitated by crowding. Our findings are of broad impact for DNA looping in the crowded and compartmentalized interior of living biological cells.
It is well established that spontaneous parametric down-conversion with induced coherence across two coupled interferometers results in high-visibility single-photon interference. We describe experiments in which additional photon channels are introduced such that "which-path" information is made possible and the fringe visibility in single-photon interference is reduced in accordance with basic notions of complementarity. However, these additional pathways result in nearly perfect visibility when photons are counted in coincidence. A simplified theoretical model accounts for these observations and attributes them directly to the vacuum fields at the different crystals.
We investigate hybrid charge transfer states (HCTS) at the planar interface between a-NPD and ZnO by spectrally resolved electroluminescence (EL) and external quantum efficiency (EQE) measurements. Radiative decay of HCTSs is proven by distinct emission peaks in the EL spectra of such bilayer devices in the NIR at energies well below the bulk a-NPD or ZnO emission. The EQE spectra display low energy contributions clearly red-shifted with respect to the a-NPD photocurrent and partially overlapping with the EL emission. Tuning of the energy gap between the ZnO conduction band and a-NPD HOMO level (E-int) was achieved by modifying the ZnO surface with self-assembled monolayers based on phosphonic acids. We find a linear dependence of the peak position of the NIR EL on E-int, which unambiguously attributes the origin of this emission to radiative recombination between an electron on the ZnO and a hole on a-NPD. In accordance with this interpretation, we find a strictly linear relation between the open-circuit voltage and the energy of the charge state for such hybrid organicinorganic interfaces.
The discovery of ultrafast X-ray induced optical reflectivity changes enabled the development of X-ray/optical cross correlation techniques at X-ray free electron lasers worldwide. We have now linked through experiment and theory the fundamental excitation and relaxation steps with the transient optical properties in finite solid samples. Therefore, we gain a thorough interpretation and an optimized detection scheme of X-ray induced changes to the refractive index and the X-ray/optical cross correlation response. (C) 2015 AIP Publishing LLC.
We present the results of 71.6 hr of observations of the Geminga pulsar (PSR J0633+1746) with the VERITAS very-high-energy gamma-ray telescope array. Data taken with VERITAS between 2007 November and 2013 February were phase-folded using a Geminga pulsar timing solution derived from data recorded by the XMM-Newton and Fermi-LAT space telescopes. No significant pulsed emission above 100 GeV is observed, and we report upper limits at the 95% confidence level on the integral flux above 135 GeV (spectral analysis threshold) of 4.0x10(-13) s(-1) cm(-2) and 1.7 x 10(-13) s(-1) cm(-2) for the two principal peaks in the emission profile. These upper limits, placed in context with phase-resolved spectral energy distributions determined from 5 yr of data from the Fermi-Large Area Telescope (LAT), constrain possible hardening of the Geminga pulsar emission spectra above similar to 50 GeV.
Particle diffusion and localized acceleration in inhomogeneous AGN jets - I. Steady-state spectra
(2015)
We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton (SSC) emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axisymmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard. gamma-ray spectra observed in some blazars.
The rhythmic activity observed in nervous systems, in particular in epilepsies and Parkinson's disease, has often been hypothesized to originate from a macroscopic self-sustained neural oscillator. However, this assumption has not been tested experimentally. Here we support this viewpoint with in vivo experiments in a rodent model of absence seizures, by demonstrating frequency locking to external periodic stimuli and finding the characteristic Arnold tongue. This result has important consequences for developing methods for the control of brain activity, such as seizure cancellation.
The PAC2MAN mission
(2015)
An accurate forecast of flare and coronal mass ejection (CME) initiation requires precise measurements of the magnetic energy buildup and release in the active regions of the solar atmosphere. We designed a new space weather mission that performs such measurements using new optical instruments based on the Hanle and Zeeman effects. The mission consists of two satellites, one orbiting the L1 Lagrangian point (Spacecraft Earth, SCE) and the second in heliocentric orbit at 1AU trailing the Earth by 80 degrees (Spacecraft 80, SC80). Optical instruments measure the vector magnetic field in multiple layers of the solar atmosphere. The orbits of the spacecraft allow for a continuous imaging of nearly 73% of the total solar surface. In-situ plasma instruments detect solar wind conditions at 1AU and ahead of our planet. Earth-directed CMEs can be tracked using the stereoscopic view of the spacecraft and the strategic placement of the SC80 satellite. Forecasting of geoeffective space weather events is possible thanks to an accurate surveillance of the magnetic energy buildup in the Sun, an optical tracking through the interplanetary space, and in-situ measurements of the near-Earth environment.
Hybrid multijunction solar cells comprising hydrogenated amorphous silicon and an organic bulk heterojunction are presented, reaching 11.7% power conversion efficiency. The benefits of merging inorganic and organic subcells are pointed out, the optimization of the cells, including optical modeling predictions and tuning of the recombination contact are described, and an outlook of this technique is given.
PCPDTBT, a marginally crystallizable polymer, is crystallized into a new crystal structure using solvent-vapor annealing. Highly ordered areas with three different polymer-chain orientations are identified using TEM/ED, GIWAXS, and polarized Raman spectroscopy. The optical and structural properties differ significantly from films prepared by standard device preparation protocols. Bilayer solar cells, however, show similar performance.
When arranged in a proper nanoaggregate architecture, gold nanoparticles can offer controllable plasmon-related absorption/scattering, yielding distinct color effects that depend critically on the relative orientation and distance between nanoparticle constituents. Herein, we report on the implementation of novel plasmonic nanoarchitectures based on complexes between gold nanoparticles and an azobenzene-modified cationic surfactant that can exhibit a light-tunable plasmonic response. The formation of such complexes becomes possible through the use of strongly negatively charged bare gold nanoparticles (similar to 10-nm diameter) prepared by the method of laser ablation in deionized water. Driven by electrostatic interactions, the cationic surfactant molecules attach and form a shell around the negatively charged nanoparticles, resulting in neutralization of the particle charge or even overcompensation beyond which the nanoparticles become positively charged. At low and high surfactant concentrations, Au nanoparticles are negatively and positively charged, respectively, and are represented by single species due to electric repulsion effects having absorption peaks around 523-527 nm, whereas at intermediate concentrations, the Au nanoparticles become neutral, forming nanoscale 100-nm clusterlike aggregates and exhibiting an additional absorption peak at gimel > 600 nm and a visible change in the color of the solution from red to blue. Because of the presence of the photosensitive azobenzene unit in the surfactant tail that undergoes trans-to-cis isomerization under irradiation with UV light, we then demonstrate a light-controlled nanoclustering of nanoparticles, yielding a switch in the plasmonic absorption band and a related change in the solution color. The formed hybrid architectures with a light-controlled plasmonic response could be important for a variety of tasks, including biomedical, surface-enhanced Raman spectroscopy (SERS), data transmission, and storage applications.
A rather strong mean longitudinal magnetic field of the order of a few hundred gauss was detected a few years ago in the Of?p star CPD -28 degrees 2561 using FORS2 (FOcal Reducer low dispersion Spectrograph 2) low-resolution spectropolarimetric observations. In this work, we present additional low-resolution spectropolarimetric observations obtained during several weeks in 2013 December using FORS 2 mounted at the 8-m Antu telescope of the Very Large Telescope (VLT). These observations cover a little less than half of the stellar rotation period of 73.41 d mentioned in the literature. The behaviour of the mean longitudinal magnetic field is consistent with the assumption of a single-wave variation during the stellar rotation cycle, indicating a dominant dipolar contribution to the magnetic field topology. The estimated polar strength of the surface dipole B-d is larger than 1.15 kG. Further, we compared the behaviour of the line profiles of various elements at different rotation phases associated with different magnetic field strengths. The strongest contribution of the emission component is observed at the phases when the magnetic field shows a negative or positive extremum. The comparison of the spectral behaviour of CPD -28 degrees 2561 with that of another Of?p star, HD 148937 of similar spectral type, reveals remarkable differences in the degree of variability between both stars. Finally, we present new X-ray observations obtained with the Suzaku X-ray Observatory. We report that the star is X-ray bright with log L-X/L-bol approximate to -5.7. The low-resolution X-ray spectra reveal the presence of a plasma heated up to 24 MK. We associate the 24 MK plasma in CPD -28 degrees 2561 with the presence of a kG strong magnetic field capable to confine stellar wind.
Liquid water molecules interact strongly with each other, forming a fluctuating hydrogen bond network and thereby giving rise to the anomalous phase diagram of liquid water. Consequently, symmetric and asymmetric water molecules have been found in the picosecond time average with IR and optical Raman spectroscopy. With subnatural linewidth resonant inelastic x-ray scattering (RIXS) at vibrational resolution, we take sub-femtosecond snapshots of the electronic and structural properties of water molecules in the hydrogen bond network. We derive a strong dominance of nonsymmetric molecules in liquid water in contrast to the gas phase on the sub-femtosecond timescale of RIXS and determine the fraction of highly asymmetrically distorted molecules.
An annihilation signal of dark matter is searched for from the central region of the Milky Way. Data acquired in dedicated on-off observations of the Galactic center region with H.E.S.S. are analyzed for this purpose. No significant signal is found in a total of similar to 9 h of on-off observations. Upper limits on the velocity averaged cross section, <sigma upsilon >, for the annihilation of dark matter particles with masses in the range of similar to 300 GeV to similar to 10 TeV are derived. In contrast to previous constraints derived from observations of the Galactic center region, the constraints that are derived here apply also under the assumption of a central core of constant dark matter density around the center of the Galaxy. Values of <sigma upsilon > that are larger than 3 x 10(-24) cm(3)/s are excluded for dark matter particles with masses between similar to 1 and similar to 4 TeV at 95% C.L. if the radius of the central dark matter density core does not exceed 500 pc. This is the strongest constraint that is derived on <sigma upsilon > for annihilating TeV mass dark matter without the assumption of a centrally cusped dark matter density distribution in the search region.
Femtosecond x-ray laser pulses are used to probe the carbon monoxide (CO) oxidation reaction on ruthenium (Ru) initiated by an optical laser pulse. On a time scale of a few hundred femtoseconds, the optical laser pulse excites motions of CO and oxygen (O) on the surface, allowing the reactants to collide, and, with a transient close to a picosecond (ps), new electronic states appear in the OK-edge x-ray absorption spectrum. Density functional theory calculations indicate that these result from changes in the adsorption site and bond formation between CO and O with a distribution of OC-O bond lengths close to the transition state (TS). After 1 ps, 10% of the CO populate the TS region, which is consistent with predictions based on a quantum oscillator model.
Supernova remnants (SNRs) are among the most important targets for gamma-ray observatories. Being prominent non-thermal sources, they are very likely responsible for the acceleration of the bulk of Galactic cosmic rays (CRS). To firmly establish the SNR paradigm for the origin of cosmic rays, it should be confirmed that protons are indeed accelerated in, and released from, SNRs with the appropriate flux and spectrum. This can be done by detailed theoretical models which account for microphysics of acceleration and various radiation processes of hadrons and leptons. The current generation of Cherenkov telescopes has insufficient sensitivity to constrain theoretical models. A new facility, the Cherenkov Telescope Array (CTA), will have superior capabilities and may finally resolve this long standing issue of high-energy astrophysics. We want to assess the capabilities of CTA to reveal the physics of various types of SNRs in the initial 2000 years of their evolution. During this time, the efficiency to accelerate cosmic rays is highest. We perform time-dependent simulations of the hydrodynamics, the magnetic fields, the cosmic-ray acceleration, and the non-thermal emission for type Ia, Ic and IIP SNRs. We calculate the CTA response to the y-ray emission from these SNRs for various ages and distances, and we perform a realistic analysis of the simulated data. We derive distance limits for the detectability and resolvability of these SNR types at several ages. We test the ability of CTA to reconstruct their morphological and spectral parameters as a function of their distance. Finally, we estimate how well CTA data will constrain the theoretical models. (C) 2014 Elsevier B.V. All rights reserved.
While N-body simulations suggest a cuspy profile in the centra of the dark matter halos of galaxies, the majority of astronomical observations favor a relatively soft cored density distribution of these regions. The routine method of testing the convergence of N-body simulations (in particular, the negligibility of two-body scattering effect) is to find the conditions under which formed structures is insensitive to numerical parameters. The results obtained with this approach suggest a surprisingly minor role of the particle collisions: the central density profile remains untouched and close to the Navarro-Frenk-White shape, even if the simulation time significantly exceeds the collisional relaxation time tau(r). In order to check the influence of the unphysical test body collisions we use the Fokker-Planck equation. It turns out that a profile rho proportional to r(-beta) where beta similar or equal to 1 is an attractor: the Fokker-Planck diffusion transforms any reasonable initial distribution into it in a time shorter than tau(r), and then the cuspy profile should survive much longer than tau(r), since the Fokker-Planck diffusion is self-compensated if beta similar or equal to 1. Thus the purely numerical effect of test body scattering may create a stable NFW-like pseudosolution. Moreover, its stability may be mistaken for the simulation convergence. We present analytical estimations for this potential bias effect and call for numerical tests. For that purpose, we suggest a simple test that can be performed as the simulation progresses and would indicate the magnitude of the collisional influence and the veracity of the simulation results. (C) 2014 Elsevier B.V. All rights reserved.