Numerical simulations of the possible origin of the two sub-parsec scale and counter-rotating stellar disks around SgrA*

We present a high resolution simulation of an idealized model to explain the origin of the two young, counter-rotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the Center of the Milky Way. In our model, the collision of a single molecular cloud with a circum-nuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing towards the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circum-nuclear disk material. The infalling gas creates two inclined, counter-rotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second, counter-rotating disk forms. Fragmentation of the second disk would lead to the two inclined, counter-rotating stellar disks which are observed at the Galactic Center. A similar event might be happening again right now at the Milky Way Galactic Center. Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic Center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so called mini-spiral that has been detected in the Galactic Center.Comment: 14 pages, 12 figures, submitted to Ap

Time-resolved infrared emission from radiation-driven central obscuring structures in Active Galactic Nuclei

The central engines of Seyfert galaxies are thought to be enshrouded by geometrically thick gas and dust structures. In this article, we derive observable properties for a self-consistent model of such toroidal gas and dust distributions, where the geometrical thickness is achieved and maintained with the help of X-ray heating and radiation pressure due to the central engine. Spectral energy distributions (SEDs) and images are obtained with the help of dust continuum radiative transfer calculations with RADMC-3D. For the first time, we are able to present time-resolved SEDs and images for a physical model of the central obscurer. Temporal changes are mostly visible at shorter wavelengths, close to the combined peak of the dust opacity as well as the central source spectrum and are caused by variations in the column densities of the generated outflow. Due to the three-component morphology of the hydrodynamical models -- a thin disc with high density filaments, a surrounding fluffy component (the obscurer) and a low density outflow along the rotation axis -- we find dramatic differences depending on wavelength: whereas the mid-infrared images are dominated by the elongated appearance of the outflow cone, the long wavelength emission is mainly given by the cold and dense disc component. Overall, we find good agreement with observed characteristics, especially for those models, which show clear outflow cones in combination with a geometrically thick distribution of gas and dust, as well as a geometrically thin, but high column density disc in the equatorial plane.Comment: 16 pages, 12 figures, accepted for publication in MNRA

Comment on The Evidence for a Pentaquark and Kinematic Reflections

The Regge exchange model used by Dzierba et al. is shown to be questionable, since the pion pole term is not allowed. Hence the Regge amplitudes in their calculation are exaggerated. The amount of kinematic reflection in the mass spectrum of the (nK+) system, which is one decay channel of a possible pentaquark, is not well justified in the fitting procedure used by Dzierba et al., as shown by comparison with the (K+K-) invariant mass spectrum, which is one decay channel of the a_2 and f_2 tensor mesons. While kinematic reflections are still a concern in some papers that have presented evidence for the pentaquark, better quantitative calculations are needed to demonstrate the significance of this effect.Comment: Comment submitted to Phys. Rev. D (no figures

Kinematics of massive star ejecta in the Milky Way as traced by 26^26Al

Context. Massive stars form in groups and their winds and supernova explosions create superbubbles up to kpc in size. The fate of their ejecta is of vital importance for the dynamics of the interstellar medium, for chemical evolution models, and the chemical enrichment of galactic halos and the intergalactic medium. However, ejecta kinematics and the characteristic scales in space and time have not been explored in great detail beyond ~10 Ka. Aims: Through measurement of radioactive 26Al with its decay time constant at ~106 years, we aim to trace the kinematics of cumulative massive-star and supernova ejecta independent of the uncertain gas parameters over million-year time scales. Our goal is to identify the mixing time scale and the spatio-kinematics of such ejecta from the pc to kpc scale in our Milky Way. Methods: We use the SPI spectrometer on the INTEGRAL observatory and its observations along the Galactic ridge to trace the detailed line shape systematics of the 1808.63 keV gamma-ray line from 26Al decay. We determine line centroids and compare these to Doppler shift expectations from large-scale systematic rotation around the Galaxy centre, as observed in other Galactic objects. Results: We measure the radial velocities of gas traced by 26Al, averaged over the line of sight, as a function of Galactic longitude. We find substantially higher velocities than expected from Galactic rotation, the average bulk velocity being ~200 km s-1 larger than predicted from Galactic rotation. The observed radial velocity spread implies a Doppler broadening of the gamma-ray line that is consistent with our measurements of the overall line width. We can reproduce the observed characteristics with 26Al sources located along the inner spiral arms, when we add a global blow-out preference into the forward direction away from arms into the inter-arm region, as is expected when massive stars are offset towards the spiral-arm leading edge. With the known connection of superbubbles to the gaseous halo, this implies angular-momentum transfer in the disk-halo system and consequently also radial gas flows. The structure of the interstellar gas above the disk affects how ionizing radiation may escape and ionize intergalactic gas.Peer reviewe

The Structure of Dark Matter Haloes in Dwarf Galaxies

Recent observations indicate that dark matter haloes have flat central density profiles. Cosmological simulations with non-baryonic dark matter predict however self similar haloes with central density cusps. This contradiction has lead to the conclusion that dark matter must be baryonic. Here it is shown that the dark matter haloes of dwarf spiral galaxies represent a one parameter family with self similar density profiles. The observed global halo parameters are coupled with each other through simple scaling relations which can be explained by the standard cold dark matter model if one assumes that all the haloes formed from density fluctuations with the same primordial amplitude. We find that the finite central halo densities correlate with the other global parameters. This result rules out scenarios where the flat halo cores formed subsequently through violent dynamical processes in the baryonic component. These cores instead provide important information on the origin and nature of dark matter in dwarf galaxies.Comment: uuencoded Z-compressed postscript file, 10 pages, 3 figures included, to appear in ApJ Letter

The life cycle of starbursting circumnuclear gas discs

High-resolution observations from the sub-mm to the optical wavelength regime resolve the central few 100pc region of nearby galaxies in great detail. They reveal a large diversity of features: thick gas and stellar discs, nuclear starbursts, in- and outflows, central activity, jet interaction, etc. Concentrating on the role circumnuclear discs play in the life cycles of galactic nuclei, we employ 3D adaptive mesh refinement hydrodynamical simulations with the RAMSES code to self-consistently trace the evolution from a quasi-stable gas disc, undergoing gravitational (Toomre) instability, the formation of clumps and stars and the disc's subsequent, partial dispersal via stellar feedback. Our approach builds upon the observational finding that many nearby Seyfert galaxies have undergone intense nuclear starbursts in their recent past and in many nearby sources star formation is concentrated in a handful of clumps on a few 100pc distant from the galactic centre. We show that such observations can be understood as the result of gravitational instabilities in dense circumnuclear discs. By comparing these simulations to available integral field unit observations of a sample of nearby galactic nuclei, we find consistent gas and stellar masses, kinematics, star formation and outflow properties. Important ingredients in the simulations are the self-consistent treatment of star formation and the dynamical evolution of the stellar distribution as well as the modelling of a delay time distribution for the supernova feedback. The knowledge of the resulting simulated density structure and kinematics on pc scale is vital for understanding inflow and feedback processes towards galactic scales.Comment: accepted by MNRA

Model Analysis of the ep→ep′π+π−ep \to ep'\pi^+\pi^- Electroproduction Reaction on the Proton

Recent CLAS data on the $p\pi^+\pi^-$ electroproduction off protons at 1.3$<$W$<$1.57 GeV and 0.25$<$$Q^{2}$$<$0.6 GeV$^{2}$ have been analyzed using a meson-baryon phenomenological model. By fitting nine 1-fold differential cross section data for each $W$ and $Q^{2}$ bin, the charged double pion electroproduction mechanisms are identified from their manifestations in the observables. We have extracted the cross sections from amplitudes of each of the considered isobar channels as well as from their coherent sum. We also obtained non-resonant partial wave amplitudes of all contributing isobar channels which could be useful for advancing a complete coupled-channel analysis of all meson electroproduction data.Comment: Experiment Numbers: E93-006, E94-005 Group: Hall

Connecting Angular Momentum and Galactic Dynamics: The complex Interplay between Spin, Mass, and Morphology

The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use the Magneticum simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that (1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.Comment: 25 pages, accepted by ApJ, www.magneticum.or

Can multistate dark matter annihilation explain the high-energy cosmic ray lepton anomalies?

Multistate dark matter (DM) models with small mass splittings and couplings to light hidden sector bosons have been proposed as an explanation for the PAMELA/Fermi/H.E.S.S. high-energy lepton excesses. We investigate this proposal over a wide range of DM density profiles, in the framework of concrete models with doublet or triplet dark matter and a hidden SU(2) gauge sector that mixes with standard model hypercharge. The gauge coupling is bounded from below by the DM relic density, and the Sommerfeld enhancement factor is explicitly computable for given values of the DM and gauge boson masses M, mu and the (largest) dark matter mass splitting delta M_{12}. Sommerfeld enhancement is stronger at the galactic center than near the Sun because of the radial dependence of the DM velocity profile, which strengthens the inverse Compton (IC) gamma ray constraints relative to usual assumptions. We find that the PAMELA/Fermi/H.E.S.S. lepton excesses are marginally compatible with the model predictions, and with CMB and Fermi gamma ray constraints, for M ~ 800 GeV, mu ~ 200 MeV, and a dark matter profile with noncuspy Einasto parameters alpha > 0.20, r_s ~ 30 kpc. We also find that the annihilating DM must provide only a subdominant (< 0.4) component of the total DM mass density, since otherwise the boost factor due to Sommerfeld enhancement is too large.Comment: 20 pages, 12 figures; v2: Corrected branching ratio for ground state DM annihilations into leptons, leading to boost factors that are larger than allowed. Added explicit results for doublet DM model. Some conclusions changed; main conclusion of tension between inverse Compton constraints and N-body simulations of halo profiles is unchange

Self-Consistent Data Analysis of the Proton Structure Function g1 and Extraction of its Moments

The reanalysis of all available world data on the longitudinal asymmetry A|| is presented. The proton structure function g1 was extracted within a unique framework of data inputs and assumptions. These data allowed for a reliable evaluation of moments of the structure function g1 in the Q2 range from 0.2 up to 30 GeV2. The Q2 evolution of the moments was studied in QCD by means of Operator Product Expansion (OPE).Comment: Proceeding of 3rd International Symposium on the Gerasimov-Drell-Hearn Sum Rule and its extensions, Old Dominion University, Norfolk, Virginia June 2-5, 200