-
Metal-Poor Stars in the MW Disk: Resonant Cooling of Vertical Oscillations of Halo Stars in Barred Galaxies
Authors:
Xingchen Li,
Isaac Shlosman,
Daniel Pfenniger,
Clayton Heller
Abstract:
Using numerical simulations of a barred disk galaxy embedded in nonspinning and spinning dark matter (DM) halos, we present a novel mechanism of `cooling' the vertical oscillations of halo DM particles, which acquire the disk kinematics. The underlying mechanism consists of resonant interactions between halo particles and the stellar bar. The cooling mechanism acts both on dynamical and secular ti…
▽ More
Using numerical simulations of a barred disk galaxy embedded in nonspinning and spinning dark matter (DM) halos, we present a novel mechanism of `cooling' the vertical oscillations of halo DM particles, which acquire the disk kinematics. The underlying mechanism consists of resonant interactions between halo particles and the stellar bar. The cooling mechanism acts both on dynamical and secular timescales, i.e., from ~ 0.5 Gyr to few Gyr, and the stellar bar acts to absorb the kinetic energy of the vertical motions. Using a Milky Way-type stellar halo, we estimate the population of metal-poor disk stars which have been trapped by the MW disk and analyze its kinematics. We find that the population of metal-poor MW disk stars with $|z|\ltorder 3$\,kpc detected by the Gaia DR3 and other surveys can have their origin in the stellar halo. The cooled population also migrates radially outwards compared by acquiring energy from the spinning bar, and prograde-moving stars have a different distribution from the retrograde ones. Next, we have calculated the ratio of the prograde-to-retrograde orbits of the cooled population and found that this ratio varies radially, with the fast-spinning stellar halo resulting in the shallower radial increase of this ratio outside of the corotation. The nonspinning stellar halo shows a monotonic increase of this ratio with radius outside the corotation. Together with analyzed radial migration of these halo stars, the cooling phenomenon of halo metal-poor stars can explain their current disk population, and has corollaries for the chemical evolution of disk galaxies in general.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Evolution of Stellar Bars in Spinning Dark Matter Halos and Stellar Bulges
Authors:
Xingchen Li,
Isaac Shlosman,
Daniel Pfenniger,
Clayton Heller
Abstract:
We use high-resolution numerical simulations to follow the barred disk evolution in a suite of models with progressively more massive stellar bulges, with bulge-to-total (disk$+$bulge) mass ratios of $B/T\sim 0-0.25$, embedded in dark matter (DM) halos with the spin $λ\sim 0 - 0.09$. We focus on models with a sequence of initial rotational support for bulges, and analyze their spinup and spindown.…
▽ More
We use high-resolution numerical simulations to follow the barred disk evolution in a suite of models with progressively more massive stellar bulges, with bulge-to-total (disk$+$bulge) mass ratios of $B/T\sim 0-0.25$, embedded in dark matter (DM) halos with the spin $λ\sim 0 - 0.09$. We focus on models with a sequence of initial rotational support for bulges, and analyze their spinup and spindown. We find that (1) the presence of a bulge affects the evolution of stellar bars, i.e., the timescale of bar instability, bar pattern speed and its decay, and the vertical buckling instability. The bar strength is nearly independent of $B/T$ in halos with spin $λ= 0$, and is suppressed by a factor $\sim 2$ for halos with $λ= 0.09$; (2) The main effect of the bulge is the destruction of the harmonic core which affects the buckling; (3) The bulge plays a minor role in the exchange of angular momentum between the barred disk and the DM halo, during its spinup and spindown; (4) Most interestingly, the buckling process triggers different response above/below the disk midplane, which anti-correlates with the bulge mass; (5) In spinning halos, the buckling process has a prolonged amplitude tail, extending by few Gyr, as verified by measuring distortions in the Laplace plane; (6) Furthermore, as verified by orbital spectral analysis, the bulge gains its spin from the bar mainly via the inner Lindblad resonance, while losing it via a number of resonances lying between the outer and inner Lindblad resonance.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
Five methods for determining pattern speeds in galaxies
Authors:
Daniel Pfenniger,
Kanak Saha,
Yu-Ting Wu
Abstract:
Aims. New or refined methods for determining instantaneous scalar and vector pattern speeds from a restricted domain are developed, for applications in N-body simulations or in galaxies.
Methods. The general feature used throughout follows from the fact that the time-derivative of a function of the coordinates is linearly proportional to its rotation rate and its particle velocities. Knowing the…
▽ More
Aims. New or refined methods for determining instantaneous scalar and vector pattern speeds from a restricted domain are developed, for applications in N-body simulations or in galaxies.
Methods. The general feature used throughout follows from the fact that the time-derivative of a function of the coordinates is linearly proportional to its rotation rate and its particle velocities. Knowing these allow in general to retrieve the instantaneous pattern speed vector by linear optimization. Similarly, if an invariant function depends both on the position and velocities, then its instantaneous rotation vectors in space can be retrieved. Knowing the accelerations allows to find rotation in velocity space.
Results. The first 3 methods are based on the assumed rotational invariance of functions at each point in space or velocity space:
1) The 6D invariant function method measuring the pattern speed vectors in space and velocity space,
2) The differential/regional 3D Tremaine-Weinberg method evaluated over high signal to noise ratio regions,
3) The 3D Jacobi integral method yielding the potential pattern speed.
Extensions to derive the rotation center position, speed and acceleration are introduced in 1) and 3). The last 2 methods are based on the assumed invariance of average functions of the particle coordinates.
4) The moment of inertia methods by using the derivative of the singular value decomposition,
5) The 2D Fourier method (3D for m = 2 mode) giving the mode rotation speeds. Pattern speed accelerations are also derived in 4) and 5).
Conclusions. Depending on the available data, the different methods provide a choice of approaches.
△ Less
Submitted 5 March, 2023; v1 submitted 6 February, 2023;
originally announced February 2023.
-
Stellar Bars in Spinning Halos: Delayed Buckling and Absence of Slowdown
Authors:
Xingchen Li,
Isaac Shlosman,
Clayton Heller,
Daniel Pfenniger
Abstract:
We use high resolution numerical simulations in order to analyze the stellar bar evolution in spinning dark matter (DM) halos. Previous works have shown that the halo spin has a substantial effect on the bar evolution and can lead to bar dissolution following the vertical buckling instability. Here, we invoke the DM spin sequence, $λ=0-0.09$, and study the effect of DM density along this $λ$-seque…
▽ More
We use high resolution numerical simulations in order to analyze the stellar bar evolution in spinning dark matter (DM) halos. Previous works have shown that the halo spin has a substantial effect on the bar evolution and can lead to bar dissolution following the vertical buckling instability. Here, we invoke the DM spin sequence, $λ=0-0.09$, and study the effect of DM density along this $λ$-sequence by varying the compactness of DM halo. We find that (1) varying the DM density has a profound effect on the stellar bar evolution along the $λ$-sequence, namely, on its amplitude, pattern speed, buckling time, etc.; (2) For $λ\gtrsim 0.04$, the buckling instability has been delayed progressively, and does not occur when the bar has reached its maximal strength; (3) Instead, stellar bars remain near maximal strength, and their amplitude plateau stage extends over $\sim 1-7$ Gyr, terminating with the buckling instability; (4) Although stellar bars remain strong during the plateau, their pattern speed stays nearly constant. The reason for this unusual behavior of stellar bars follows from the highly reduced gravitational torques which they experience due to the DM bar being aligned with the stellar bar. The performed orbital analysis shows that the delayed buckling results from a slow evolution of stellar oscillations along the bar major and vertical axes -- thus postponing the action of the vertical 2:1 resonance which pumps the rotational energy into vertical motions; (5) Peanut/boxy shaped bulges form at the beginning of the plateau and grow with time; (6) Strong stellar bars in spinning halos can avoid fast braking, resolving the long standing discrepancy between observations and $N$-body simulations. This behavior of stellar bars along the $λ$- and DM density-sequences, reveals a wealth of stellar bar properties which require additional study.
△ Less
Submitted 12 September, 2023; v1 submitted 8 November, 2022;
originally announced November 2022.
-
The Origin of Buckling Instability in Galactic Bars: Searching for the Scapegoat
Authors:
Xingchen Li,
Isaac Shlosman,
Daniel Pfenniger,
Clayton Heller
Abstract:
The buckling process in stellar bars is full of unsolved issues. We analyze the origin of the buckling instability in stellar bars using high-resolution N-body simulations. Previous studies have promoted the nonresonant firehose instability to be responsible for the vertical buckling. We have analyzed the buckling process in terms of the resonant excitation of stellar orbits in the bar, which pump…
▽ More
The buckling process in stellar bars is full of unsolved issues. We analyze the origin of the buckling instability in stellar bars using high-resolution N-body simulations. Previous studies have promoted the nonresonant firehose instability to be responsible for the vertical buckling. We have analyzed the buckling process in terms of the resonant excitation of stellar orbits in the bar, which pumps energy into vertical oscillations. We find that (1) the buckling is associated with an abrupt increase in the central mass concentration and triggers velocities along the bar and along its rotation axis. The velocity field projected on one of the main axes forms circulation cells and increases vorticity, which are absent in firehose instability; (2) The bending amplitude is nonlinear when measured by isodensity contours or curvature of the Laplace plane, which has a substantial effect on the stellar motions; (3) In the linear description, the planar and vertical 2:1 resonances appear only with the buckling and quickly reach the overlapping phase, thus supporting the energy transfer; (4) Using nonlinear orbit analysis, we analyze the stellar oscillations along the bar and along the rotation axis and find that stars cross the vertical 2:1 resonance simultaneously with the buckling. The overlapping planar and vertical 2:1 resonances trapping more than 25% of the bar particles provide the 'smoking gun' pointing to a close relationship between the bending of stellar orbits and the resonant action -- these particles provide the necessary ingredient assuring the cohesive response in the growing vertical asymmetry. We conclude that resonant excitation is important in triggering the buckling instability, and the contribution from the firehose instability should be reevaluated. Finally, we discuss some observational implications of buckling.
△ Less
Submitted 2 March, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
Angular momentum bounds in particle systems
Authors:
Daniel Pfenniger
Abstract:
Four expressions involving sums of position and velocity coordinates bounding the total angular momentum of particle systems, and by extension of any continuous or discontinuous material systems, are derived which are tighter for any particle configuration than similar inequalities derived by Sundman (1913), Saari (2005) and Scheeres (2012). Eight distinct inequalities can thus be ordered accordin…
▽ More
Four expressions involving sums of position and velocity coordinates bounding the total angular momentum of particle systems, and by extension of any continuous or discontinuous material systems, are derived which are tighter for any particle configuration than similar inequalities derived by Sundman (1913), Saari (2005) and Scheeres (2012). Eight distinct inequalities can thus be ordered according to their tightness to angular momentum.
△ Less
Submitted 2 December, 2019;
originally announced December 2019.
-
Molecular clouds in the Cosmic Snake normal star-forming galaxy 8 billion years ago
Authors:
Miroslava Dessauges-Zavadsky,
Johan Richard,
Françoise Combes,
Daniel Schaerer,
Wiphu Rujopakarn,
Lucio Mayer,
Antonio Cava,
Frédéric Boone,
Eiichi Egami,
Jean-Paul Kneib,
Pablo G. Pérez-González,
Daniel Pfenniger,
Tim D. Rawle,
Romain Teyssier,
Paul P. van der Werf
Abstract:
The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme s…
▽ More
The cold molecular gas in contemporary galaxies is structured in discrete cloud complexes. These giant molecular clouds (GMCs), with $10^4$-$10^7$ solar masses and radii of 5-100 parsecs, are the seeds of star formation. Highlighting the molecular gas structure at such small scales in distant galaxies is observationally challenging. Only a handful of molecular clouds were reported in two extreme submillimetre galaxies at high redshift. Here we search for GMCs in a typical Milky Way progenitor at z = 1.036. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we mapped the CO(4-3) emission of this gravitationally lensed galaxy at high resolution, reading down to 30 parsecs, which is comparable to the resolution of CO observations of nearby galaxies. We identify 17 molecular clouds, characterized by masses, surface densities and supersonic turbulence all of which are 10-100 times higher than present-day analogues. These properties question the universality of GMCs and suggest that GMCs inherit their properties from ambient interstellar medium. The measured cloud gas masses are similar to the masses of stellar clumps seen in the galaxy in comparable numbers. This corroborates the formation of molecular clouds by fragmentation of distant turbulent galactic gas disks, which then turn into stellar clumps ubiquitously observed in galaxies at cosmic noon.
△ Less
Submitted 17 September, 2019;
originally announced September 2019.
-
Time-dependent Pattern Speeds in Barred Galaxies
Authors:
Yu-Ting Wu,
Daniel Pfenniger,
Ronald E. Taam
Abstract:
Based on a high quality $N$-body simulation of a double bar galaxy model, we investigate the evolution of the bar properties, including their size, strength and instantaneous pattern speed derived by using three distinct methods: the Fourier, Jacobi integral, and moment of inertia methods. The interaction of the two bars, which rotate at distinct speeds, primarily affects the size, strength and pa…
▽ More
Based on a high quality $N$-body simulation of a double bar galaxy model, we investigate the evolution of the bar properties, including their size, strength and instantaneous pattern speed derived by using three distinct methods: the Fourier, Jacobi integral, and moment of inertia methods. The interaction of the two bars, which rotate at distinct speeds, primarily affects the size, strength and pattern speed of the inner bar. When the two bars are perpendicular to each other, the size and the pattern speed of the inner bar decrease and its strength increases. The emergence of a strong Fourier $m=1$ mode increases the oscillation amplitude of the size, strength and pattern speed of the inner bar. On the other hand, the characteristics of the outer bar are substantially influenced by its adjacent spiral structure. When the spiral structure disappears, the size of the outer bar increases and its strength and pattern speed decrease. Consequently, the ratio of the pattern speed of the outer bar with respect to the inner bar is not constant and increases with time. Overall, the double bar and disk system displays substantial high frequency semi-chaotic fluctuations of the pattern strengths and speeds both in space and time, superposed on the slow secular evolution, which invalidates the assumption that the actions of individual stars should be well conserved in barred galaxies, such as the Milky Way.
△ Less
Submitted 18 May, 2018;
originally announced May 2018.
-
Solid H2 in the interstellar medium
Authors:
Andreas Füglistaler,
Daniel Pfenniger
Abstract:
Condensation of H 2 in the interstellar medium (ISM) has long been seen as a possibility, either by deposition on dust grains or thanks to a phase transition combined with self-gravity. H 2 condensation might explain the observed low efficiency of star formation and might help to hide baryons in spiral galaxies.
Our aim is to quantify the solid fraction of H 2 in the ISM due to a phase transitio…
▽ More
Condensation of H 2 in the interstellar medium (ISM) has long been seen as a possibility, either by deposition on dust grains or thanks to a phase transition combined with self-gravity. H 2 condensation might explain the observed low efficiency of star formation and might help to hide baryons in spiral galaxies.
Our aim is to quantify the solid fraction of H 2 in the ISM due to a phase transition including self-gravity for different densities and temperatures in order to use the results in more complex simulations of the ISM as subgrid physics.
We used molecular dynamics simulations of fluids at different temperatures and densities to study the formation of solids. Once the simulations reached a steady state, we calculated the solid mass fraction, energy increase, and timescales. By determining the power laws measured over several orders of magnitude, we extrapolated to lower densities the higher density fluids that can be simulated with current computers.
The solid fraction and energy increase of fluids in a phase transition are above 0.1 and do not follow a power law. Fluids out of a phase transition are still forming a small amount of solids due to chance encounters of molecules. The solid mass fraction and energy increase of these fluids are linearly dependent on density and can easily be extrapolated. The timescale is below one second, the condensation can be considered instantaneous.
The presence of solid H 2 grains has important dynamic implications on the ISM as they may be the building blocks for larger solid bodies when gravity is included. We provide the solid mass fraction, energy increase, and timescales for high density fluids and extrapolation laws for lower densities.
△ Less
Submitted 1 March, 2018; v1 submitted 4 December, 2017;
originally announced December 2017.
-
Time-dependent Corotation Resonance in Barred Galaxies
Authors:
Yu-Ting Wu,
Daniel Pfenniger,
Ronald E. Taam
Abstract:
The effective potential neighboring the corotation resonance region in barred galaxies is shown to be strongly time-dependent in any rotating frame because of the competition of nearby perturbations of similar strengths with differing rotation speeds. Contrary to the generally adopted assumption, that in the bar rotating frame the corotation region should possess four stationary equilibrium points…
▽ More
The effective potential neighboring the corotation resonance region in barred galaxies is shown to be strongly time-dependent in any rotating frame because of the competition of nearby perturbations of similar strengths with differing rotation speeds. Contrary to the generally adopted assumption, that in the bar rotating frame the corotation region should possess four stationary equilibrium points (Lagrange points), with high quality N-body simulations we localize the instantaneous equilibrium points and find that they circulate or oscillate broadly in azimuth with respect to the pattern speeds of the inner or outer perturbations. This implies that at the particle level the Jacobi integral is not well conserved around the corotation radius. That is, angular momentum exchanges decouple from energy exchanges, enhancing the chaotic diffusion of stars through the corotation region.
△ Less
Submitted 3 August, 2016;
originally announced August 2016.
-
Episodic Model For Star Formation History and Chemical Abundances in Giant and Dwarf Galaxies
Authors:
Suma Debsarma,
Tanuka Chattopadhyay,
Sukanta Das,
Daniel Pfenniger
Abstract:
In search for a synthetic understanding, a scenario for the evolution of the star formation rate and the chemical abundances in galaxies is proposed, combining gas infall from galactic halos, outflow of gas by supernova explosions, and an oscillatory star formation process. The oscillatory star formation model is a consequence of the modelling of the fractional masses changes of the hot, warm and…
▽ More
In search for a synthetic understanding, a scenario for the evolution of the star formation rate and the chemical abundances in galaxies is proposed, combining gas infall from galactic halos, outflow of gas by supernova explosions, and an oscillatory star formation process. The oscillatory star formation model is a consequence of the modelling of the fractional masses changes of the hot, warm and cold components of the interstellar medium. The observed periods of oscillation vary in the range $(0.1-3.0)\times10^{7}$\,yr depending on various parameters existing from giant to dwarf galaxies. The evolution of metallicity varies in giant and dwarf galaxies and depends on the outflow process. Observed abundances in dwarf galaxies can be reproduced under fast outflow together with slow evaporation of cold gases into hot gas whereas slow outflow and fast evaporation is preferred for giant galaxies. The variation of metallicities in dwarf galaxies supports the fact that low rate of SNII production in dwarf galaxies is responsible for variation in metallicity in dwarf galaxies of similar masses as suggested by various authors.
△ Less
Submitted 31 July, 2016;
originally announced August 2016.
-
Formation of H2-He substellar bodies in cold conditions: Gravitational stability of binary mixtures in a phase transition
Authors:
Andreas Füglistaler,
Daniel Pfenniger
Abstract:
Molecular clouds typically consist of 3/4 H2, 1/4 He and traces of heavier elements. In an earlier work we showed that at very low temperatures and high densities, H2 can be in a phase transition leading to the formation of ice clumps as large as comets or even planets. However, He has very different chemical properties and no phase transition is expected before H2 in dense interstellar medium (IS…
▽ More
Molecular clouds typically consist of 3/4 H2, 1/4 He and traces of heavier elements. In an earlier work we showed that at very low temperatures and high densities, H2 can be in a phase transition leading to the formation of ice clumps as large as comets or even planets. However, He has very different chemical properties and no phase transition is expected before H2 in dense interstellar medium (ISM) conditions. The gravitational stability of fluid mixtures has been studied before, but these studies did not include a phase transition.
First, we study the gravitational stability of van der Waals fluid mixtures using linearized analysis and examine virial equilibrium conditions using the Lennard-Jones intermolecular potential. Then, combining the Lennard-Jones and gravitational potentials, the non-linear dynamics of fluid mixtures are studied via computer simulations using the molecular dynamics code LAMMPS.
Along with the classical, ideal-gas Jeans instability criterion, a fluid mixture is always gravitationally unstable if it is in a phase transition because compression does not increase pressure. However, the condensed phase fraction increases. In unstable situations the species can separate: in some conditions He precipitates faster than H2, while in other conditions the converse occurs. Also, for an initial gas phase collapse the geometry is essential. Contrary to spherical or filamentary collapses, sheet-like collapses starting below 15K easily reach H2 condensation conditions because then they are fastest and both the increase of heating and opacity are limited.
Depending on density, temperature and mass, either rocky H2 planetoids, or gaseous He planetoids form. H2 planetoids are favoured by high density, low temperature and low mass, while He planetoids need more mass and can form at temperature well above the critical value.
△ Less
Submitted 29 June, 2016; v1 submitted 16 July, 2015;
originally announced July 2015.
-
Structure formation in warm dark matter cosmologies: Top-Bottom Upside-Down
Authors:
Sinziana Paduroiu,
Yves Revaz,
Daniel Pfenniger
Abstract:
The damping on the fluctuation spectrum and the presence of thermal velocities as properties of warm dark matter particles like sterile neutrinos imprint a distinct signature found from the structure formation mechanisms to the internal structures of halos. Using warm dark matter simulations we explore these effects on the structure formation for different particle energies and we find that the fo…
▽ More
The damping on the fluctuation spectrum and the presence of thermal velocities as properties of warm dark matter particles like sterile neutrinos imprint a distinct signature found from the structure formation mechanisms to the internal structures of halos. Using warm dark matter simulations we explore these effects on the structure formation for different particle energies and we find that the formation of structure is more complex than originally assumed, a combination of top-down collapse and hierarchical (bottom-up) clustering on multiple scales. The degree on which one scenario is more prominent with respect to the other depends globally on the energy of the particle and locally on the morphology and architecture of the analyzed region. The presence of shells and caustics in warm dark matter halos is another important effect seen in simulations. Furthermore, we discuss the impact of thermal velocities on the structure formation from theoretical considerations as well as from the analysis of the simulations. We re-examine the assumptions considered when estimating the velocity dispersion for warm dark matter particles that have been adopted in previous works for more than a decade and we give an independent estimation for the velocities. We identify some inconsistencies in previous published results. The relation between the warm dark matter particle mass and its corresponding velocity dispersion is strongly model dependent, hence the constraints on particle mass from simulation results are weak. Finally, we review the technical difficulties that arise in warm dark matter simulations along with possible improvements of the methods.
△ Less
Submitted 23 June, 2015; v1 submitted 11 June, 2015;
originally announced June 2015.
-
Substellar fragmentation in self-gravitating fluids with a major phase transition
Authors:
Andreas Füglistaler,
Daniel Pfenniger
Abstract:
The existence of substellar cold H2 globules in planetary nebulae and the mere existence of comets suggest that the physics of cold interstellar gas might be much richer than usually envisioned.
We study the case of a cold gaseous medium in ISM conditions which is subject to a gas-liquid/solid phase transition.
First the equilibrium of general non-ideal fluids is studied using the virial the…
▽ More
The existence of substellar cold H2 globules in planetary nebulae and the mere existence of comets suggest that the physics of cold interstellar gas might be much richer than usually envisioned.
We study the case of a cold gaseous medium in ISM conditions which is subject to a gas-liquid/solid phase transition.
First the equilibrium of general non-ideal fluids is studied using the virial theorem and linear stability analysis. Then the non-linear dynamics is studied by using simulations to characterize the expected formation of solid bodies analogous to comets. The simulations are run with a state of the art molecular dynamics code (LAMMPS). The long-range gravitational forces can be taken into account with short-range molecular forces with finite limited computational resources by using super-molecules, provided the right scaling is followed.
The concept of super-molecule is tested with simulations, allowing us to correctly satisfy the Jeans instability criterion for one-phase fluids. The simulations show that fluids presenting a phase transition are gravitationally unstable as well, independent of the strength of the gravitational potential, producing two distinct kinds of sub-stellar bodies, those dominated by gravity ("planetoids") and those dominated by molecular attractive force ("comets").
Observations, formal analysis and computer simulations suggest the possibility of the formation of substellar H2 clumps in cold molecular clouds due to the combination of phase transition and gravity. Fluids presenting a phase transition are gravitationally unstable, independent of the strength of the gravitational potential. Arbitrarily small H2 clumps may form even at relatively high temperatures up to 400 - 600K, according to virial analysis. The combination of phase transition and gravity may be relevant for a wider range of astrophysical situations, such as proto-planetary disks.
△ Less
Submitted 16 April, 2015; v1 submitted 16 March, 2015;
originally announced March 2015.
-
Numerical Studies of Hamiltonian Systems and Application to Galactic Potentials
Authors:
Daniel Pfenniger
Abstract:
The talk consisted mainly in commenting in a linear way the seminal paper in 1964 by Michel Henon and graduate student Carl Heiles at Princeton University: "The applicability of the third integral of motion: Some numerical experiments" in the field of galactic dynamics. Instead of repeating here the lecture of the paper, we advise the reader interested in dynamical systems to study this "must" rea…
▽ More
The talk consisted mainly in commenting in a linear way the seminal paper in 1964 by Michel Henon and graduate student Carl Heiles at Princeton University: "The applicability of the third integral of motion: Some numerical experiments" in the field of galactic dynamics. Instead of repeating here the lecture of the paper, we advise the reader interested in dynamical systems to study this "must" reading.
△ Less
Submitted 17 November, 2014;
originally announced November 2014.
-
CO map and steep Kennicutt-Schmidt relation in the extended UV disk of M63
Authors:
M. Dessauges-Zavadsky,
C. Verdugo,
F. Combes,
D. Pfenniger
Abstract:
Results from the UV satellite GALEX revealed large extensions of disks in some nearby spiral galaxies, extending out to 3 to 4 times the isophotal radius, r25. M63 is a remarkable example of a spiral galaxy with one of the most extended UV disks, so it offers the opportunity to search for the molecular gas and characterize the star formation in outer disk regions as revealed by the UV emission. We…
▽ More
Results from the UV satellite GALEX revealed large extensions of disks in some nearby spiral galaxies, extending out to 3 to 4 times the isophotal radius, r25. M63 is a remarkable example of a spiral galaxy with one of the most extended UV disks, so it offers the opportunity to search for the molecular gas and characterize the star formation in outer disk regions as revealed by the UV emission. We obtained deep CO(1-0) and CO(2-1) observations on the IRAM 30 m telescope along the major axis of the M63 disk from the center out to the galactocentric radius rgal = 1.6 r25 and over a bright UV region at rgal = 1.36 r25. CO(1-0) is detected all along the M63 major axis out to r25, and CO(2-1) is confined to rgal = 0.68 r25, which may betray lower excitation temperatures in the outer disk. CO(1-0) is also detected in the external bright UV region of M63. The radial profiles of the CO emission and of the Halpha, 24 micron, NUV and FUV star formation tracers and HI taken from the literature show a severe drop with the galactocentric radius, such that beyond r25 they are all absent with the exception of a faint UV emission and HI. The CO emission detection in the external UV region, where the UV flux is higher than the UV flux observed beyond r25, highlights a tight correlation between the CO and UV fluxes, namely the amount of molecular gas and the intensity of star formation. This external UV region is dominated by the atomic gas, suggesting that HI is more likely the precursor of H2 rather than the product of UV photodissociation. A broken power law needs to be invoked to describe the Kennicutt-Schmidt (K-S) relation of M63 from the center of the galaxy out to rgal = 1.36 r25. While all along the major axis out to r25 the K-S relation is almost linear, in the external UV region the SFR regime is highly nonlinear and characterized by a steep K-S relation and very low star formation efficiency.
△ Less
Submitted 2 June, 2014;
originally announced June 2014.
-
Meridional tilt of the stellar velocity ellipsoid during bar buckling instability
Authors:
Kanak Saha,
Daniel Pfenniger,
Ronald E. Taam
Abstract:
The structure and evolution of the stellar velocity ellipsoid plays an important role in shaping galaxies undergoing bar driven secular evolution and the eventual formation of a boxy/peanut bulge such as present in the Milky Way. Using collisionless N-body simulations, we show that during the formation of such a boxy/peanut bulge, the meridional shear stress of stars, which can be measured by the…
▽ More
The structure and evolution of the stellar velocity ellipsoid plays an important role in shaping galaxies undergoing bar driven secular evolution and the eventual formation of a boxy/peanut bulge such as present in the Milky Way. Using collisionless N-body simulations, we show that during the formation of such a boxy/peanut bulge, the meridional shear stress of stars, which can be measured by the meridional tilt of the velocity ellipsoid, reaches a characteristic peak in its time evolution. It is shown that the onset of a bar buckling instability is closely connected to the maximum meridional tilt of the stellar velocity ellipsoid. Our findings bring new insight to this complex gravitational instability of the bar which complements the buckling instability studies based on orbital models. We briefly discuss the observed diagnostics of the stellar velocity ellipsoid during such a phenomenon.
△ Less
Submitted 24 December, 2012; v1 submitted 6 November, 2012;
originally announced November 2012.
-
On the detection of point sources in CMB maps
Authors:
H. G. Khachatryan,
G. Nurbaeva,
D. Pfenniger,
G. Meylan
Abstract:
We use the Wilkinson Microwave Anisotropy Probe 7-year data (WMAP7) to further probe point source detection technique in the sky maps of the cosmic microwave background (CMB) radiation. The method by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of WMAP satellite CMB temperature data. Part of the detected points coincide with…
▽ More
We use the Wilkinson Microwave Anisotropy Probe 7-year data (WMAP7) to further probe point source detection technique in the sky maps of the cosmic microwave background (CMB) radiation. The method by Tegmark et al. for foreground reduced maps and the Kolmogorov parameter as the descriptor are adopted for the analysis of WMAP satellite CMB temperature data. Part of the detected points coincide with point sources already revealed by other methods. However, we have also found 2 source candidates for which still no counterparts are known, and identified 7 point sources listed in Planck Early Release Compact Source Catalogue as high reliability sources.
△ Less
Submitted 29 July, 2013; v1 submitted 29 June, 2012;
originally announced June 2012.
-
Stellar diffusion in barred spiral galaxies
Authors:
Maura Brunetti,
Cristina Chiappini,
Daniel Pfenniger
Abstract:
We characterize empirically the radial diffusion of stars in the plane of a typical barred disk galaxy by calculating the local spatial diffusion coefficient and diffusion time-scale for bulge-disk-halo N-body self-consistent systems which initially differ in the Safronov-Toomre-Q_T parameter. We find different diffusion scenarios that depend on the bar strength and on the degree of instability of…
▽ More
We characterize empirically the radial diffusion of stars in the plane of a typical barred disk galaxy by calculating the local spatial diffusion coefficient and diffusion time-scale for bulge-disk-halo N-body self-consistent systems which initially differ in the Safronov-Toomre-Q_T parameter. We find different diffusion scenarios that depend on the bar strength and on the degree of instability of the disk. Marginally stable disks, with Q_T \sim 1, have two families of bar orbits with different values of angular momentum and energy, which determine a large diffusion in the corotation region. In hot disks, Q_T> 1, stellar diffusion is reduced with respect to the case of marginally stable disks. In cold models, we find that spatial diffusion is not constant in time and strongly depends on the activity of the bar, which can move stars all over the disk recurrently. We conclude that to realistically study the impact of radial migration on the chemical evolution modeling of the Milky Way the role of the bar has to be taken into account.
△ Less
Submitted 29 August, 2011;
originally announced August 2011.
-
Diffusion in barred-spiral galaxies
Authors:
Maura Brunetti,
Cristina Chiappini,
Daniel Pfenniger
Abstract:
We characterize the radial migration of stars in the disk plane by calculating the diffusion coefficient and the diffusion time-scale for a bulge-disk N-body self-consistent system with a marginally-stable Toomre-Q parameter. We find that diffusion is not constant in time, but follows the evolution of the bar, and becomes maximum near the corotation region and in the external disk region, where as…
▽ More
We characterize the radial migration of stars in the disk plane by calculating the diffusion coefficient and the diffusion time-scale for a bulge-disk N-body self-consistent system with a marginally-stable Toomre-Q parameter. We find that diffusion is not constant in time, but follows the evolution of the bar, and becomes maximum near the corotation region and in the external disk region, where asymmetric patterns develop.
△ Less
Submitted 3 December, 2010;
originally announced December 2010.
-
Modeling proper motions beyond the Galactic bulge
Authors:
Maura Brunetti,
Daniel Pfenniger
Abstract:
We analyse the radial and tangential velocity fields in the Galaxy as seen from the Sun by using as a first approximation a simple axisymmetric model, which we then compare with the corresponding fields in a barred N-body model of the Milky Way. This provides a global description of these quantities missing in the literature, showing where they take large values susceptible to be used in future…
▽ More
We analyse the radial and tangential velocity fields in the Galaxy as seen from the Sun by using as a first approximation a simple axisymmetric model, which we then compare with the corresponding fields in a barred N-body model of the Milky Way. This provides a global description of these quantities missing in the literature, showing where they take large values susceptible to be used in future observations even for sources well beyond the Galactic center. Absolute largest proper motions occur at a distance slightly behind the Galactic Center, which are there 1.5 times larger than the highest local proper motions due to the Galactic differential rotation. Large proper motions well beyond the Galactic center are well within the current astrometric accuracy.
△ Less
Submitted 16 December, 2009; v1 submitted 12 November, 2009;
originally announced November 2009.
-
Simulations of galactic disks including an additional dark baryonic component
Authors:
Y. Revaz,
D. Pfenniger,
F. Combes,
F. Bournaud
Abstract:
The near proportionality between HI and dark matter in outer galactic disks prompted us to run N-body simulations of galactic disks in which the observed gas content is supplemented by a dark gas component representing between zero and five times the visible gas content. While adding baryons in the disk of galaxies may solve some issues, it poses the problem of disk stability. We show that the g…
▽ More
The near proportionality between HI and dark matter in outer galactic disks prompted us to run N-body simulations of galactic disks in which the observed gas content is supplemented by a dark gas component representing between zero and five times the visible gas content. While adding baryons in the disk of galaxies may solve some issues, it poses the problem of disk stability. We show that the global stability is ensured if the ISM is multiphased, composed of two partially coupled phases, a visible warm gas phase and a weakly collisionless cold dark phase corresponding to a fraction of the unseen baryons. The phases are subject to stellar and UV background heating and gas cooling, and their transformation into each other is studied as a function of the coupling strength. This new model, which still possesses a dark matter halo, fits the rotation curves as well as the classical CDM halos, but is the only one to explain the existence of an open and contrasting spiral structure, as observed in the outer HI disks
△ Less
Submitted 29 April, 2009;
originally announced April 2009.
-
Chaos in Galaxies
Authors:
Daniel Pfenniger
Abstract:
After general considerations about limits of theories and models, where small changes may imply large effects, we discuss three cases in galactic astrophysics illustrating how galactic dynamics models may become insufficient when previously neglected effects are taken into account: 1) Like in 3D hydrodynamics, the non-linearity of the Poisson-Boltzmann system may imply dissipation through the gr…
▽ More
After general considerations about limits of theories and models, where small changes may imply large effects, we discuss three cases in galactic astrophysics illustrating how galactic dynamics models may become insufficient when previously neglected effects are taken into account: 1) Like in 3D hydrodynamics, the non-linearity of the Poisson-Boltzmann system may imply dissipation through the growth of discontinuous solutions. 2) The relationship between the microscopic exponential sensitivity of N-body systems and the stability of mean field galaxy models. 3) The role of quantum physics in the dynamics of structure formation, considering that cosmological neutrinos are massive and semi-degenerate fermions.
△ Less
Submitted 22 February, 2008;
originally announced February 2008.
-
Simulations of galactic disks including a dark baryonic component
Authors:
Y. Revaz,
D. Pfenniger,
F. Combes,
F. Bournaud
Abstract:
$Λ$CDM numerical simulations predict that the "missing baryons" reside in a Warm-Hot gas phase in the over-dense cosmic filaments. However, there are now several theoretical and observational arguments that support the fact that galactic disks may be more massive than usually thought, containing a substantial fraction of the "missing baryons". Hereafter, we present new N-body simulations of gala…
▽ More
$Λ$CDM numerical simulations predict that the "missing baryons" reside in a Warm-Hot gas phase in the over-dense cosmic filaments. However, there are now several theoretical and observational arguments that support the fact that galactic disks may be more massive than usually thought, containing a substantial fraction of the "missing baryons". Hereafter, we present new N-body simulations of galactic disks, where the gas content has been multiplied by a factor 5. The stability of the disk is ensured by assuming that the ISM is composed out of two partially coupled phases, a warm phase, corresponding the observed CO and HI gas and a cold collisionless phase corresponding to the unseen baryons.
△ Less
Submitted 8 January, 2008;
originally announced January 2008.
-
Molecular gas in high-velocity clouds: revisited scenario
Authors:
M. Dessauges-Zavadsky,
F. Combes,
D. Pfenniger
Abstract:
We report a new search for 12CO(1-0) emission in high-velocity clouds (HVCs) performed with the IRAM 30 m telescope. This search was motivated by the recent detection of cold dust emission in the HVCs of Complex C. Despite a spatial resolution which is three times better and sensitivity twice as good compared to previous studies, no CO emission is detected in the HVCs of Complex C down to a best…
▽ More
We report a new search for 12CO(1-0) emission in high-velocity clouds (HVCs) performed with the IRAM 30 m telescope. This search was motivated by the recent detection of cold dust emission in the HVCs of Complex C. Despite a spatial resolution which is three times better and sensitivity twice as good compared to previous studies, no CO emission is detected in the HVCs of Complex C down to a best 5 sigma limit of 0.16 K km/s at a 22'' resolution. The CO emission non-detection does not provide any evidence in favor of large amounts of molecular gas in these HVCs and hence in favor of the infrared findings. We discuss different configurations which, however, allow us to reconcile the negative CO result with the presence of molecular gas and cold dust emission. H2 column densities higher than our detection limit, N(H2) = 3x10^{19} cm^{-2}, are expected to be confined in very small and dense clumps with 20 times smaller sizes than the 0.5 pc clumps resolved in our observations according to the results obtained in cirrus clouds, and might thus still be highly diluted. As a consequence, the inter-clump gas at the 1 pc scale has a volume density lower than 20 cm^{-3} and already appears as too diffuse to excite the CO molecules. The observed physical conditions in the HVCs of Complex C also play an important role against CO emission detection. It has been shown that the CO-to-H2 conversion factor in low metallicity media is 60 times higher than at the solar metallicity, leading for a given H2 column density to a 60 times weaker integrated CO intensity. And the very low dust temperature estimated in these HVCs implies the possible presence of gas cold enough (< 20 K) to cause CO condensation onto dust grains under interstellar medium pressure conditions and thus CO depletion in gas-phase observations.
△ Less
Submitted 8 August, 2007;
originally announced August 2007.
-
Virial statistical description of non-extensive hierarchical systems
Authors:
Daniel Pfenniger
Abstract:
In a first part the scope of classical thermodynamics and statistical mechanics is discussed in the broader context of formal dynamical systems, including computer programmes. In this context classical thermodynamics appears as a particular theory suited to a subset of all dynamical systems. A statistical mechanics similar to the one derived with the microcanonical ensemble emerges from dynamica…
▽ More
In a first part the scope of classical thermodynamics and statistical mechanics is discussed in the broader context of formal dynamical systems, including computer programmes. In this context classical thermodynamics appears as a particular theory suited to a subset of all dynamical systems. A statistical mechanics similar to the one derived with the microcanonical ensemble emerges from dynamical systems provided it contains, 1) a finite non-integrable part of its phase space which is, 2) ergodic at a satisfactory degree after a finite time. The integrable part of phase space provides the constraints that shape the particular system macroscopical properties, and the chaotic part provides well behaved statistical properties over a relevant finite time. More generic semi-ergodic systems lead to intermittent behaviour, thus may be unsuited for a statistical description of steady states. Following these lines of thought, in a second part non-extensive hierarchical systems with statistical scale-invariance and power law interactions are explored. Only the virial constraint, consistent with their microdynamics, is included. No assumptions of classical thermodynamics are used, in particular extensivity and local homogeneity. In the limit of a large hierarchical range new constraints emerge in some conditions that depend on the interaction law range. In particular for the gravitational case, a velocity-site scaling relation is derived which is consistant with the ones empirically observed in the fractal interstellar medium.
△ Less
Submitted 26 May, 2006;
originally announced May 2006.
-
Cosmological Neutrino Entanglement and Quantum Pressure
Authors:
Daniel Pfenniger,
Veruska Muccione
Abstract:
Context: The widespread view that cosmological neutrinos, even if massive, are well described since the decoupling redshift z~10^10 down to the present epoch by an almost perfectly collisionless fluid of classical point particles is re-examined. Aims: In view of the likely sub-eV rest mass of neutrinos, the main effects due to their fermionic nature are studied. Methods: By numerical means we ca…
▽ More
Context: The widespread view that cosmological neutrinos, even if massive, are well described since the decoupling redshift z~10^10 down to the present epoch by an almost perfectly collisionless fluid of classical point particles is re-examined. Aims: In view of the likely sub-eV rest mass of neutrinos, the main effects due to their fermionic nature are studied. Methods: By numerical means we calculate the accurate entropy, fugacity and pressure of cosmological neutrinos in the Universe expansion. By solving the Schroedinger equation we derive how and how fast semi-degenerate identical free fermions become entangled. Results: We find that for sub-eV neutrinos the exchange degeneracy has significantly increased during the relativistic to non-relativistic transition epoch at z~10^4-10^5. At all times neutrinos become entangled in less than 10^-6 s, much faster than any plausible decoherence time. The total pressure is increased by quantum effect from 5% at high redshifts to 68% at low redshifts with respect to a collisionless classical fluid. Conclusions: The quantum overpressure has no dynamical consequences in the homogeneous regime at high redshifts, but must be significant for neutrino clustering during the non-linear structure formation epoch at low redshifts.
△ Less
Submitted 19 June, 2006; v1 submitted 15 May, 2006;
originally announced May 2006.
-
A new scenario for the origin of galactic warps
Authors:
Yves Revaz,
Daniel Pfenniger
Abstract:
Galactic warps represent an old unresolved problem, since the discovery, at the end of the fifties, of the HI warp of the Milky Way. In this paper, we propose a new scenario explaining a large fraction of the observed optical warps. Based on N-body simulations, we show that realistic galactic disks, where the dark matter is essentially distributed in a disk, are subject to bending instabilities.…
▽ More
Galactic warps represent an old unresolved problem, since the discovery, at the end of the fifties, of the HI warp of the Milky Way. In this paper, we propose a new scenario explaining a large fraction of the observed optical warps. Based on N-body simulations, we show that realistic galactic disks, where the dark matter is essentially distributed in a disk, are subject to bending instabilities. S, U-shaped, as well as asymmetric warps are spontaneously generated and in some cases are long-lived. While this scenario presents the advantage of explaining the three observed types of warps, it also brings new constraints on the dark matter distribution in spiral galaxies. Finally, it gives us a unified picture of galaxies where galactic asymmetries, like bars, spirals and warps result from gravitational instabilities.
△ Less
Submitted 12 September, 2005;
originally announced September 2005.
-
Bending instabilities at the origin of persistent warps
Authors:
Y. Revaz,
D. Pfenniger
Abstract:
Based on N-body simulations, we show that realistic galactic disks are subject to bending instabilities of fire-hose type when the disks are substantially self-gravitating, that is, if they contain dark matter distributed in the disk. Depending on the degree of instability, S and U-shaped, as well as asymmetric warps are generated. In some cases, the warp may last several galactic rotations, par…
▽ More
Based on N-body simulations, we show that realistic galactic disks are subject to bending instabilities of fire-hose type when the disks are substantially self-gravitating, that is, if they contain dark matter distributed in the disk. Depending on the degree of instability, S and U-shaped, as well as asymmetric warps are generated. In some cases, the warp may last several galactic rotations, particularly when the instability is marginal. Since the bending instability is very sensitive to the disk flattening, the fractions of dark matter distributed in the disk and in the dark halo are constrained. For a Milky Way like galaxy the extended dark halo can not exceed 30-40% of the total mass within 35 kpc if the Milky Way warp results from a bending instability. This mode of warping provides a unified picture of spiral galaxies, where bars, spiral arms and warps result all from disk gravitational instabilities, radial or transverse, which are constantly regenerated by the dissipative gas component.
△ Less
Submitted 16 December, 2004;
originally announced December 2004.
-
The Baryonic Tully-Fisher relation revisited
Authors:
D. Pfenniger,
Y. Revaz
Abstract:
The Baryonic Tully-Fisher relation (BTF) can be substantially improved when considering that the galactic baryonic mass is likely to consist not only from the detected baryons, stars and gas, but also from a dark baryonic component proportional to the HI gas. The BTF relation is optimally improved when the HI mass is multiplied by a factor of about 3, but larger factors up to 11-16 still improve…
▽ More
The Baryonic Tully-Fisher relation (BTF) can be substantially improved when considering that the galactic baryonic mass is likely to consist not only from the detected baryons, stars and gas, but also from a dark baryonic component proportional to the HI gas. The BTF relation is optimally improved when the HI mass is multiplied by a factor of about 3, but larger factors up to 11-16 still improve the fit over the original one using only the detected baryons. The strength of this improved relation is quantified with up-to-date statistical tests such as the Akaike Information Criterion or the Bayesian Information Criterion. In particular they allow to show that supposing a variable $M_\star/L$ ratio instead is much less significant. This result reinforces the suggestion made in several recent works that mass within galactic disks must be a multiple of the HI mass, and that galactic disks are substantially, but not necessarily fully, self-gravitating.
△ Less
Submitted 27 September, 2004;
originally announced September 2004.
-
A Unified Picture of Disk Galaxies where Bars, Spirals and Warps Result from the Same Fundamental Causes
Authors:
Daniel Pfenniger,
Yves Revaz
Abstract:
Bars and spiral arms have played an important role as constraints on the dynamics and on the distribution of dark matter in the optical parts of disk galaxies. Dynamics linked to the dissipative nature of gas, and its transformation into stars provide clues that spiral galaxies are driven by dissipation close to a state of \textit{marginal stability} with respect to the dynamics in the galaxy pl…
▽ More
Bars and spiral arms have played an important role as constraints on the dynamics and on the distribution of dark matter in the optical parts of disk galaxies. Dynamics linked to the dissipative nature of gas, and its transformation into stars provide clues that spiral galaxies are driven by dissipation close to a state of \textit{marginal stability} with respect to the dynamics in the galaxy plane. Here we present numerical evidences that warps play a similar role but in the transverse direction. N-body simulations show that typical galactic disks are also marginally stable with respect to a bending instability leading to typical observed warps. The frequent occurrence of warps and asymmetries in the outer galactic disks give therefore, like bars in the inner disks, new constraints on the dark matter, but this time in the outer disks. If disks are marginally stable with respect to bending instabilities, our models suggest that the mass within the HI disks must be a multiple of the detected HI and stars, i.e., disks must be heavier than seen. But the models do not rule out a traditional thick halo with a mass within the HI disk radius similar to the total disk mass.
△ Less
Submitted 25 June, 2004;
originally announced June 2004.
-
Bending instabilities at the origin of persistent warps: a new constraint on dark matter halos
Authors:
Yves Revaz,
Daniel Pfenniger
Abstract:
A substantial fraction of the warps in spiral galaxies may result from bending instabilities if the disks are essentially self-gravitating. With N-body simulations, we show that galaxies with self-gravitating disks as thick as HI disks are subject to bending instabilities generating S-shaped, U-shaped or asymmetric warps. S-shaped warps persist during several rotations and keep the line of node…
▽ More
A substantial fraction of the warps in spiral galaxies may result from bending instabilities if the disks are essentially self-gravitating. With N-body simulations, we show that galaxies with self-gravitating disks as thick as HI disks are subject to bending instabilities generating S-shaped, U-shaped or asymmetric warps. S-shaped warps persist during several rotations and keep the line of node straight. The warp amplitudes generated by bending instabilities remain however modest. Other factors must be invoked for extreme warped disks. However, bending instabilities can account for most of the cases reported in optical surveys, where the warp angle is generally less than $5^\circ$. This mode of warping is very sensitive to the disk flattening. It also constrains the fraction of dark matter distributed in the disk and in the dark halo.
△ Less
Submitted 15 June, 2004;
originally announced June 2004.
-
Scale Interactions and Galaxy Evolution
Authors:
Daniel Pfenniger
Abstract:
To understand galaxies and their evolution, it is necessary to describe how the different scales interact: how the microscopic physics, such as star formation, or the large scale physics, such as galaxy interactions may modify the galaxy global shapes. The purpose of this review is to point out some general or recent topics related to such scale interactions, both observational and theoretical,…
▽ More
To understand galaxies and their evolution, it is necessary to describe how the different scales interact: how the microscopic physics, such as star formation, or the large scale physics, such as galaxy interactions may modify the galaxy global shapes. The purpose of this review is to point out some general or recent topics related to such scale interactions, both observational and theoretical, which are relevant in the present understanding of galaxies.
△ Less
Submitted 13 December, 2002;
originally announced December 2002.
-
Possible flakes of molecular hydrogen in the early Universe
Authors:
D. Pfenniger,
D. Puy
Abstract:
The thermochemistry of H2 and HD in non-collapsed, non-reionized primordial gas up to the end of the dark age is investigated with recent radiation-matter and chemical reaction rates taking into account the efficient coolant HD, and the possibility of a gas-solid phase transition of H2. In the standard big-bang model we find that these molecules can freeze out and lead to the growth of flakes of…
▽ More
The thermochemistry of H2 and HD in non-collapsed, non-reionized primordial gas up to the end of the dark age is investigated with recent radiation-matter and chemical reaction rates taking into account the efficient coolant HD, and the possibility of a gas-solid phase transition of H2. In the standard big-bang model we find that these molecules can freeze out and lead to the growth of flakes of solid molecular hydrogen at redshifts z ~ 6-12 in the unperturbed medium and under-dense regions. While this freezing caused by the mere adiabatic cooling of the expanding matter is less likely to occur in collapsed regions due to their higher than radiation background temperature, on the other hand the super-adiabatic expansion in voids strongly favors it. Later reionization (at z ~ 5-6) eventually destroys all these H2 flakes. The possible occurrence of H2 flakes is important for the degree of coupling between matter and radiation, as well as for the existence of a gas-grain chemistry at the end of the dark age.
△ Less
Submitted 18 November, 2002;
originally announced November 2002.
-
Long-Range Correlations in Self-Gravitating N-Body Systems
Authors:
Daniel Huber,
Daniel Pfenniger
Abstract:
Observed self-gravitating systems reveal often fragmented non-equilibrium structures that feature characteristic long-range correlations. However, models accounting for non-linear structure growth are not always consistent with observations and a better understanding of self-gravitating $N$-body systems appears necessary. Because unstable gravitating systems are sensitive to non-gravitational pe…
▽ More
Observed self-gravitating systems reveal often fragmented non-equilibrium structures that feature characteristic long-range correlations. However, models accounting for non-linear structure growth are not always consistent with observations and a better understanding of self-gravitating $N$-body systems appears necessary. Because unstable gravitating systems are sensitive to non-gravitational perturbations we study the effect of different dissipative factors as well as different small and large scale boundary conditions on idealized $N$-body systems. We find, in the interval of negative specific heat, equilibrium properties differing from theoretical predictions made for gravo-thermal systems, substantiating the importance of microscopic physics and the lack of consistent theoretical tools to describe self-gravitating gas. Also, in the interval of negative specific heat, yet outside of equilibrium, unforced systems fragment and establish transient long-range correlations. The strength of these correlations depends on the degree of granularity, suggesting to make the resolution of mass and force coherent. Finally, persistent correlations appear in model systems subject to an energy flow.
△ Less
Submitted 29 January, 2002;
originally announced January 2002.
-
Lumpy Structures in Self-Gravitating Disks
Authors:
Daniel Huber,
Daniel Pfenniger
Abstract:
Following Toomre & Kalnajs (1991), local models of slightly dissipative self-gravitating disks show how inhomogeneous structures can be maintained over several galaxy rotations. Their basic physical ingredients are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes on the kpc scale, local simulation of self-gravitating disks are…
▽ More
Following Toomre & Kalnajs (1991), local models of slightly dissipative self-gravitating disks show how inhomogeneous structures can be maintained over several galaxy rotations. Their basic physical ingredients are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes on the kpc scale, local simulation of self-gravitating disks are performed in this paper in 2D as well as in 3D. The third dimension becomes a priori important as soon as matter clumping causes a tight coupling of the 3D equations of motion. The physically simple and general framework of the model permits to make conclusions beyond the here considered scales. A time dependent affine coordinate system is used, allowing to calculate the gravitational forces via a particle-mesh FFT-method, increasing the performance with respect to previous direct force calculations. Persistent patterns, formed by transient structures, whose intensity and morphological characteristic depend on the dissipation rate are obtained and described. Some of our simulations reveal first signs of mass-size and velocity dispersion-size power-law relations, but a clear scale invariant behavior will require more powerful computer techniques.
△ Less
Submitted 29 May, 2001;
originally announced May 2001.
-
Periodic orbits in warped disk
Authors:
Y. Revaz,
D. Pfenniger
Abstract:
It is often assumed that a warped galaxy can be modeled by a set of rings. This paper verifies numerically the validity of this assumption by the study of periodic orbits populating a heavy self-gravitating warped disk. The phase space structure of a warped model reveals that the circular periodic orbits of a flat disk are transformed in quasi annular periodic orbits which conserve their stabili…
▽ More
It is often assumed that a warped galaxy can be modeled by a set of rings. This paper verifies numerically the validity of this assumption by the study of periodic orbits populating a heavy self-gravitating warped disk. The phase space structure of a warped model reveals that the circular periodic orbits of a flat disk are transformed in quasi annular periodic orbits which conserve their stability. This lets us also explore the problem of the persistence of a large outer warp. In particular, the consistency of its orbits with the density distribution is checked as a function of the pattern speed.
△ Less
Submitted 10 May, 2001;
originally announced May 2001.
-
Fragmentation in Kinematically Cold Disks
Authors:
Daniel Huber,
Daniel Pfenniger
Abstract:
Gravity is scale free. Thus gravity may form similar structures in self-gravitating systems on different scales. Indeed, observations of the interstellar medium, spiral disks and cosmic structures, reveal similar characteristics. The structures in these systems are very lumpy and inhomogeneous. Moreover some of these structures do not seem to be of random nature, but obey certain power laws.
M…
▽ More
Gravity is scale free. Thus gravity may form similar structures in self-gravitating systems on different scales. Indeed, observations of the interstellar medium, spiral disks and cosmic structures, reveal similar characteristics. The structures in these systems are very lumpy and inhomogeneous. Moreover some of these structures do not seem to be of random nature, but obey certain power laws.
Models of slightly dissipative self-gravitating disks show how such inhomogeneous structures can be maintained on the kpc-scale. The basic physical processes in these models are self-gravity, dissipation and differential rotation. In order to explore the structures resulting from these processes, local simulations of self-gravitating disks are performed in 2D and 3D. We observe persistent patterns, formed by transient structures, whose intensity and morphological characteristic depend on the dissipation rate.
△ Less
Submitted 16 January, 2001;
originally announced January 2001.
-
N-body simulations of warped galaxies
Authors:
Yves Revaz,
Daniel Pfenniger
Abstract:
Two methods generating warped galaxies with N-body simulations are presented. One uses an external potential as a disturber while the other is based on material accretion. The results of both methods are compared. A particular attention is given on the shape of the line of node (LON).
Two methods generating warped galaxies with N-body simulations are presented. One uses an external potential as a disturber while the other is based on material accretion. The results of both methods are compared. A particular attention is given on the shape of the line of node (LON).
△ Less
Submitted 5 September, 2000;
originally announced September 2000.
-
Maximum Galactic Disks vs. Hot Dark Halos
Authors:
Daniel Pfenniger
Abstract:
A series of arguments is presented for heavy galaxy disks not only in the optical regions, but also in the dark matter dominated regions of spirals. We are testing this possibility with extreme maximum disk N-body models without any conventional spheroidal dark halo.
A series of arguments is presented for heavy galaxy disks not only in the optical regions, but also in the dark matter dominated regions of spirals. We are testing this possibility with extreme maximum disk N-body models without any conventional spheroidal dark halo.
△ Less
Submitted 4 September, 2000;
originally announced September 2000.
-
Scaling Laws in Self-Gravitating Disks
Authors:
Daniel Huber,
Daniel Pfenniger
Abstract:
The interstellar medium (ISM) reveals strongly inhomogeneous structures at every scale. These structures do not seem completely random since they obey certain power laws. Larson's law (\citeyear{Larson81}) $σ\propto R^δ$ and the plausible assumption of virial equilibrium justify to consider fractals as a possible description of the ISM. In the following we investigate how self-gravitation, diffe…
▽ More
The interstellar medium (ISM) reveals strongly inhomogeneous structures at every scale. These structures do not seem completely random since they obey certain power laws. Larson's law (\citeyear{Larson81}) $σ\propto R^δ$ and the plausible assumption of virial equilibrium justify to consider fractals as a possible description of the ISM. In the following we investigate how self-gravitation, differential rotation and dissipation affect the matter distribution in galaxies. To this end we have performed 3D-simulations for self-gravitating local boxes embedded in a larger disk, extending the 2D-method of Toomre & Kalnajs (\citeyear{Toomre91}) and Wisdom & Tremaine (\citeyear{Wisdom88}). Our simulations lead to the conclusion that gravitation, shearing and dissipation can be dominantly responsible for maintaining an inhomogeneous and eventually a fractal distribution of the matter.
△ Less
Submitted 16 April, 1999;
originally announced April 1999.
-
Formation of massive counter-rotating disks: an alternative scenario
Authors:
I. Puerari,
D. Pfenniger
Abstract:
In this contribution, we investigate the scenario of major mergers between galaxies of comparable mass to the formation of massive counter-rotating disks. Such a scenario can be ruled out in general because they produce rather ellipticals. Yet for a narrow range of initial conditions, major mergers can nevertheless produce anyway rare galaxies such as NGC 4550, and so they are viable alternative…
▽ More
In this contribution, we investigate the scenario of major mergers between galaxies of comparable mass to the formation of massive counter-rotating disks. Such a scenario can be ruled out in general because they produce rather ellipticals. Yet for a narrow range of initial conditions, major mergers can nevertheless produce anyway rare galaxies such as NGC 4550, and so they are viable alternatives for the formation of these counter-rotating systems. The scenario is successful in producing a remarkably axisymmetric disk possessing strongly counter-rotating populations.
△ Less
Submitted 5 March, 1999;
originally announced March 1999.
-
Which thermal physics for gravitationally unstable media?
Authors:
Daniel Pfenniger
Abstract:
We remind that the assumptions almost universally adopted among astronomers concerning the physics to use to describe rarefied cosmic gases remain often without justifications, mainly because the long range of gravitation invalidates the use of classical thermal physics. In turn, without sufficiently good local thermal equilibrium, macroscopic quantities, such as temperature and pressure, are no…
▽ More
We remind that the assumptions almost universally adopted among astronomers concerning the physics to use to describe rarefied cosmic gases remain often without justifications, mainly because the long range of gravitation invalidates the use of classical thermal physics. In turn, without sufficiently good local thermal equilibrium, macroscopic quantities, such as temperature and pressure, are not defined and the fundamental assumption that locally the medium is smoothed by ``molecular chaos'' to justify the use of differential equations is not granted. The highly inhomogeneous fractal state of the interstellar gas is probably a plain symptom of the large discrepancy between the available theoretical tools, predicting local homogeneity after a few sound crossing times, and reality. Such fundamental problems begin to occur in optically thin media such as stellar atmospheres, but become exacerbated in the interstellar medium, in cooling flows, and in the post-recombination gas, particularly when gravitation becomes energetically dominant, i.e., when the medium is Jeans unstable.
△ Less
Submitted 10 June, 1998;
originally announced June 1998.
-
Clumpuscule Formation at High Redshift
Authors:
F. Combes,
D. Pfenniger
Abstract:
Over the last two decades, realistic studies have often concluded that the first bound objects to form can be very small, much smaller than a solar mass. After recombination, the Jeans mass drops rapidly to the order of a Giant Molecular Cloud (GMC) mass ($\sim 10^5$ M$_\odot$), and the H$_2$ cooling can make the collapse quasi-isothermal; this leads to recursive fragmentation, and formation of…
▽ More
Over the last two decades, realistic studies have often concluded that the first bound objects to form can be very small, much smaller than a solar mass. After recombination, the Jeans mass drops rapidly to the order of a Giant Molecular Cloud (GMC) mass ($\sim 10^5$ M$_\odot$), and the H$_2$ cooling can make the collapse quasi-isothermal; this leads to recursive fragmentation, and formation of clumps so dense that 3-body reactions transform the gas almost entirely to the molecular phase. This could lead to star formation in some places, but since star formation is very inefficient, most of the molecular gas could consist of a fractal built on clumpuscules thermalized with the background radiation, and filling a tiny fraction of the volume.
The bulk of the gas mass can therefore be trapped in this phase, well before the first stars re-heat and re-ionize the diffuse gaseous medium}. This results in a very contrasted multi-phase baryonic medium, that has partly remained until the present time.
△ Less
Submitted 30 January, 1998;
originally announced January 1998.
-
Making Spirals with Counter-Rotating Disks, the case of NGC 4550
Authors:
D. Pfenniger
Abstract:
A single merger scenario for making galaxies such as NGC 4550 possessing equal coplanar counter-rotating stellar disks is investigated by collisionless N-body technique. The scenario is successful in producing an axisymmetric disk made of two almost equal counter-rotating populations. The final disk shows a clear bimodal line profile in the outer part, which demonstrates that disk-disk mergers d…
▽ More
A single merger scenario for making galaxies such as NGC 4550 possessing equal coplanar counter-rotating stellar disks is investigated by collisionless N-body technique. The scenario is successful in producing an axisymmetric disk made of two almost equal counter-rotating populations. The final disk shows a clear bimodal line profile in the outer part, which demonstrates that disk-disk mergers do not always produce ellipticals.
△ Less
Submitted 2 December, 1997;
originally announced December 1997.
-
ISOCAM images at 6.75 and 15 microns: the circumnuclear region of NGC 4321 (M 100)
Authors:
H. Wozniak,
D. Friedli,
L. Martinet,
D. Pfenniger
Abstract:
We present ISOCAM images (1.5"/px) in the LW2 (6.75 microns) and LW3 (15 microns) broad band filters of the circumnuclear starburst of NGC 4321 (M 100). A detailed comparison with the emission at other wavelengths (U, B, Halpha, and CO) is also presented. The nature and intensity of the mid-infrared emission is shown to be remarkably different from one star forming region to the next. For instan…
▽ More
We present ISOCAM images (1.5"/px) in the LW2 (6.75 microns) and LW3 (15 microns) broad band filters of the circumnuclear starburst of NGC 4321 (M 100). A detailed comparison with the emission at other wavelengths (U, B, Halpha, and CO) is also presented. The nature and intensity of the mid-infrared emission is shown to be remarkably different from one star forming region to the next. For instance, star forming regions at the ends of the nuclear bar may contain enshrouded Wolf-Rayet stars.
△ Less
Submitted 13 November, 1997;
originally announced November 1997.
-
First Results from ISOCAM Images of Circumnuclear Starbursts in Barred Galaxies
Authors:
H. Wozniak,
D. Friedli,
L. Martinet,
D. Pfenniger
Abstract:
For a sample of ten nearby barred galaxies, images at 6.75\um and 15\um of circumnuclear starbursts are being taken by ISOCAM on board of the ISO satellite. We first recall the main goals of our project, give its current status, and finally, as an example, preliminary results for NGC 4321 (M100) are presented.
For a sample of ten nearby barred galaxies, images at 6.75\um and 15\um of circumnuclear starbursts are being taken by ISOCAM on board of the ISO satellite. We first recall the main goals of our project, give its current status, and finally, as an example, preliminary results for NGC 4321 (M100) are presented.
△ Less
Submitted 14 August, 1997;
originally announced August 1997.
-
Perspectives for detecting cold H2 in outer galactic disks
Authors:
F. Combes,
D. Pfenniger
Abstract:
We review here the main direct or indirect ways to detect the possible presence of large amounts of cold molecular hydrogen in the outer parts of disk galaxies, an hypothesis that we have recently developed. Direct ways range from H2 absorption in the UV domain to detection of the radio hyperfine structure: the ortho-H2 molecule has an hyperfine, or ultrafine, structure in its fundamental state,…
▽ More
We review here the main direct or indirect ways to detect the possible presence of large amounts of cold molecular hydrogen in the outer parts of disk galaxies, an hypothesis that we have recently developed. Direct ways range from H2 absorption in the UV domain to detection of the radio hyperfine structure: the ortho-H2 molecule has an hyperfine, or ultrafine, structure in its fundamental state, due to the coupling between the rotation-induced magnetic moment, and the nuclear spin. This gives rise to 2 magnetic dipole transitions, at the wavelengths of 0.5 and 5.5 km. Indirect ways are essentially the detection of the HD and LiH transitions, and in some environments like clusters of galaxies, more heavy trace molecules such as CO. We discuss from this point of view the recent discovery by COBE/FIRAS of a very cold Galactic dust component (4-7 K) which could correspond to a dominating gas mass component of the ISM, if interpreted as standard dust emission. Some of the proposed means could be applied to the well-known molecular clouds, to bring some new light to the problem of the H2/CO conversion ratio.
△ Less
Submitted 18 July, 1997;
originally announced July 1997.
-
Very Cold Gas and Dark Matter
Authors:
F. Combes,
D. Pfenniger
Abstract:
We have recently proposed a new candidate for baryonic dark matter: very cold molecular gas, in near-isothermal equilibrium with the cosmic background radiation at 2.73 K. The cold gas, of quasi-primordial abundances, is condensed in a fractal structure, resembling the hierarchical structure of the detected interstellar medium.
We present some perspectives of detecting this very cold gas, eith…
▽ More
We have recently proposed a new candidate for baryonic dark matter: very cold molecular gas, in near-isothermal equilibrium with the cosmic background radiation at 2.73 K. The cold gas, of quasi-primordial abundances, is condensed in a fractal structure, resembling the hierarchical structure of the detected interstellar medium.
We present some perspectives of detecting this very cold gas, either directly or indirectly. The H$_2$ molecule has an "ultrafine" structure, due to the interaction between the rotation-induced magnetic moment and the nuclear spins. But the lines fall in the km domain, and are very weak. The best opportunity might be the UV absorption of H$_2$ in front of quasars. The unexpected cold dust component, revealed by the COBE/FIRAS submillimetric results, could also be due to this very cold H$_2$ gas, through collision-induced radiation, or solid H$_2$ grains or snowflakes. The $γ$-ray distribution, much more radially extended than the supernovae at the origin of cosmic rays acceleration, also points towards and extended gas distribution.
△ Less
Submitted 27 February, 1996;
originally announced February 1996.
-
Fundamental Aspects of the ISM Fractality
Authors:
D. Pfenniger
Abstract:
The ubiquitous clumpy state of the ISM raises a fundamental and open problem of physics, which is the correct statistical treatment of systems dominated by long range interactions. A simple solvable hierarchical model is presented which explains why systems dominated by gravity prefer to adopt a fractal dimension around 2 or less, like the cold ISM and large scale structures. This has direct rel…
▽ More
The ubiquitous clumpy state of the ISM raises a fundamental and open problem of physics, which is the correct statistical treatment of systems dominated by long range interactions. A simple solvable hierarchical model is presented which explains why systems dominated by gravity prefer to adopt a fractal dimension around 2 or less, like the cold ISM and large scale structures. This has direct relation with the general transparency, or blackness, of the Universe.
△ Less
Submitted 27 February, 1996;
originally announced February 1996.