-
Hydrodynamic simulations of WD-WD mergers and the origin of RCB stars
Authors:
Sagiv Shiber,
Orsola De Marco,
Patrick M. Motl,
Bradley Munson,
Dominic C. Marcello,
Juhan Frank,
Patrick Diehl,
Geoffrey C. Clayton,
Bennett N. Skinner,
Hartmut Kaiser,
Gregor Daiss,
Dirk Pfluger,
Jan E. Staff
Abstract:
We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code Octo-Tiger. We follow the orbital evolution of DWD systems of mass ratio q=0.7 for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, findin…
▽ More
We study the properties of double white dwarf (DWD) mergers by performing hydrodynamic simulations using the new and improved adaptive mesh refinement code Octo-Tiger. We follow the orbital evolution of DWD systems of mass ratio q=0.7 for tens of orbits until and after the merger to investigate them as a possible origin for R Coronae Borealis (RCB) type stars. We reproduce previous results, finding that during the merger, the Helium WD donor star is tidally disrupted within 20-80 minutes since the beginning of the simulation onto the accretor Carbon-Oxygen WD, creating a high temperature shell around the accretor. We investigate the possible Helium burning in this shell and the merged object's general structure. Specifically, we are interested in the amount of Oxygen-16 dredged-up to the hot shell and the amount of Oxygen-18 produced. This is critical as the discovery of very low Oxygen-16 to Oxygen-18 ratios in RCB stars pointed out the merger scenario as a favorable explanation for their origin. A small amount of hydrogen in the donor may help keep the Oxygen-16 to Oxygen-18 ratios within observational bounds, even if moderate dredge-up from the accretor occurs. In addition, we perform a resolution study to reconcile the difference found in the amount of Oxygen-16 dredge-up between smoothed-particle hydrodynamics and grid-based simulations.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Octo-Tiger: A New, 3D Hydrodynamic Code for Stellar Mergers that uses HPX Parallelisation
Authors:
Dominic C. Marcello,
Sagiv Shiber,
Orsola De Marco,
Juhan Frank,
Geoffrey C. Clayton,
Patrick M. Motl,
Patrick Diehl,
Hartmut Kaiser
Abstract:
OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs or main sequence stars. The gravity solver conserves an…
▽ More
OCTO-TIGER is an astrophysics code to simulate the evolution of self-gravitating and rotat-ing systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. OCTO-TIGER is currently optimised to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs or main sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a correction algorithm. This code uses HPX parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, white-dwarf binary simulations. Results are compared to analytic solutions, when known, and to other grid based codes such as FLASH. We also compute the interaction between two white dwarfs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure OCTO-TIGERs scaling properties up to a core count of 80,000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients.
△ Less
Submitted 10 August, 2021; v1 submitted 20 January, 2021;
originally announced January 2021.
-
Numerical Simulations of Mass Transfer in Binaries with Bipolytropic Components
Authors:
Kundan Kadam,
Patrick M. Motl,
Dominic C. Marcello,
Juhan Frank,
Geoffrey C. Clayton
Abstract:
We present the first self-consistent, three dimensional study of hydrodynamic simulations of mass transfer in binary systems with bipolytropic (composite polytropic) components. In certain systems, such as contact binaries or during the common envelope phase, the core-envelope structure of the stars plays an important role in binary interactions. In this paper, we compare mass transfer simulations…
▽ More
We present the first self-consistent, three dimensional study of hydrodynamic simulations of mass transfer in binary systems with bipolytropic (composite polytropic) components. In certain systems, such as contact binaries or during the common envelope phase, the core-envelope structure of the stars plays an important role in binary interactions. In this paper, we compare mass transfer simulations of bipolytropic binary systems in order to test the suitability of our numerical tools for investigating the dynamical behaviour of such systems. The initial, equilibrium binary models possess a core-envelope structure and are obtained using the bipolytropic self-consistent field technique. We conduct mass transfer simulations using two independent, fully three-dimensional, Eulerian codes - Flow-ER and Octo-tiger. These hydrodynamic codes are compared across binary systems undergoing unstable as well as stable mass transfer, and the former at two resolutions. The initial conditions for each simulation and for each code are chosen to match closely so that the simulations can be used as benchmarks. Although there are some key differences, the detailed comparison of the simulations suggests that there is remarkable agreement between the results obtained using the two codes. This study puts our numerical tools on a secure footing, and enables us to reliably simulate specific mass transfer scenarios of binary systems involving components with a core-envelope structure.
△ Less
Submitted 16 September, 2018; v1 submitted 13 September, 2018;
originally announced September 2018.
-
The role of dredge-up in double white dwarf mergers
Authors:
Jan E. Staff,
Brandon Wiggins,
Dominic Marcello,
Patrick M. Motl,
Wesley Even,
Christopher L. Fryer,
Cody Raskin,
Geoffrey C. Clayton,
Juhan Frank
Abstract:
We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested…
▽ More
We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested as a possible formation channel for R Coronae Borealis type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of $^{18}{\mathrm O}$, the merger must avoid dredging up large amounts of $^{16}{\mathrm O}$, or else it will be difficult to produce sufficient $^{18}{\mathrm O}$ to explain the oxygen ratio observed to be of order unity. We performed a total of 9 simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes, and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, $>10^{-2}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO white dwarf with a $\sim 0.1~{\rm M_\odot}$ layer of helium on top, we find that no $^{16}{\mathrm O}$ is being dredged up, while in the $q=0.8$ simulation $<10^{-4}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.
△ Less
Submitted 1 June, 2018;
originally announced June 2018.
-
A Numerical Method for Generating Rapidly Rotating Bipolytropic Structures in Equilibrium
Authors:
Kundan Kadam,
Patrick M. Motl,
Juhan Frank,
Geoffrey C. Clayton,
Dominic C. Marcello
Abstract:
We demonstrate that rapidly rotating bipolytropic (composite polytropic) stars and toroidal disks can be obtained using Hachisu's self consistent field technique. The core and the envelope in such a structure can have different polytropic indices and also different average molecular weights. The models converge for high $T/|W|$ cases, where T is the kinetic energy and W is the gravitational energy…
▽ More
We demonstrate that rapidly rotating bipolytropic (composite polytropic) stars and toroidal disks can be obtained using Hachisu's self consistent field technique. The core and the envelope in such a structure can have different polytropic indices and also different average molecular weights. The models converge for high $T/|W|$ cases, where T is the kinetic energy and W is the gravitational energy of the system. The agreement between our numerical solutions with known analytical as well as previously calculated numerical results is excellent. We show that the uniform rotation lowers the maximum core mass fraction or the Sch$\ddot{\rm{o}}$nberg-Chandrasekhar limit for a bipolytropic sequence. We also discuss the applications of this method to magnetic braking in low mass stars with convective envelopes.
△ Less
Submitted 31 August, 2017;
originally announced September 2017.
-
A Comparison of Grid-Based and SPH Binary Mass-Transfer and Merger Simulations
Authors:
Patrick M Motl,
Juhan Frank,
Jan Staff,
Geoffrey C Clayton,
Christopher L Fryer,
Wesley Even,
Steven Diehl,
Joel E Tohline
Abstract:
Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different nu…
▽ More
Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.
△ Less
Submitted 12 February, 2017;
originally announced February 2017.
-
Do R Coronae Borealis Stars Form from Double White Dwarf Mergers?
Authors:
Jan. E. Staff,
Athira Menon,
Falk Herwig,
Wesley Even,
Chris L. Fryer,
Patrick M. Motl,
Tom Geballe,
Marco Pignatari,
Geoffrey C. Clayton,
Joel E. Tohline
Abstract:
A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the resu…
▽ More
A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.
△ Less
Submitted 2 August, 2012;
originally announced August 2012.
-
Boosting jet power in black hole spacetimes
Authors:
David Neilsen,
Luis Lehner,
Carlos Palenzuela,
Eric W. Hirschmann,
Steven L. Liebling,
Patrick M. Motl,
T. Garrett
Abstract:
The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the i…
▽ More
The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux.
△ Less
Submitted 11 January, 2011; v1 submitted 27 December, 2010;
originally announced December 2010.
-
Mergers of Magnetized Neutron Stars with Spinning Black Holes: Disruption, Accretion and Fallback
Authors:
Sarvnipun Chawla,
Matthew Anderson,
Michael Besselman,
Luis Lehner,
Steven L. Liebling,
Patrick M. Motl,
David Neilsen
Abstract:
We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects…
▽ More
We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.
△ Less
Submitted 5 October, 2010; v1 submitted 14 June, 2010;
originally announced June 2010.
-
Perturbed disks get shocked. Binary black hole merger effects on accretion disks
Authors:
Miguel Megevand,
Matthew Anderson,
Juhan Frank,
Eric W. Hirschmann,
Luis Lehner,
Steven L. Liebling,
Patrick M. Motl,
David Neilsen
Abstract:
The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, osc…
▽ More
The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.
△ Less
Submitted 7 July, 2009; v1 submitted 20 May, 2009;
originally announced May 2009.
-
Magnetized Neutron Star Mergers and Gravitational Wave Signals
Authors:
Matthew Anderson,
Eric W. Hirschmann,
Luis Lehner,
Steven L. Liebling,
Patrick M. Motl,
David Neilsen,
Carlos Palenzuela,
Joel E. Tohline
Abstract:
We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially…
▽ More
We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay the final merger of the two stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal fields. These effects not only produce electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.
△ Less
Submitted 9 June, 2008; v1 submitted 28 January, 2008;
originally announced January 2008.
-
Simulating binary neutron stars: dynamics and gravitational waves
Authors:
Matthew Anderson,
Eric W. Hirschmann,
Luis Lehner,
Steven L. Liebling,
Patrick M. Motl,
David Neilsen,
Carlos Palenzuela,
Joel E. Tohline
Abstract:
We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial pe…
▽ More
We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.
△ Less
Submitted 4 January, 2008; v1 submitted 20 August, 2007;
originally announced August 2007.
-
Why Do Only Some Galaxy Clusters Have Cool Cores?
Authors:
Jack O. Burns,
Eric J. Hallman,
Brennan Gantner,
Patrick M. Motl,
Michael L. Norman
Abstract:
Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian hydrodynamic simulations, including radiative cooling and heating, are used to address this question as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. T…
▽ More
Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian hydrodynamic simulations, including radiative cooling and heating, are used to address this question as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. These adaptive mesh refinement simulations produce both CC and NCC clusters in the same volume. They have a peak resolution of 15.6 h^{-1} kpc within a (256 h^{-1} Mpc)^3 box. Our simulations suggest that there are important evolutionary differences between CC clusters and their NCC counterparts. Many of the numerical CC clusters accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CC's survived the collisions. By contrast, NCC clusters experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevented CC re-formation. As a result, our simulations predict observationally testable distinctions in the properties of CC and NCC beyond the core regions in clusters. In particular, we find differences between CC versus NCC clusters in the shapes of X-ray surface brightness profiles, between the temperatures and hardness ratios beyond the cores, between the distribution of masses, and between their supercluster environs. It also appears that CC clusters are no closer to hydrostatic equilibrium than NCC clusters, an issue important for precision cosmology measurements.
△ Less
Submitted 15 November, 2007; v1 submitted 14 August, 2007;
originally announced August 2007.
-
Cluster Structure in Cosmological Simulations I: Correlation to Observables, Mass Estimates, and Evolution
Authors:
Tesla E. Jeltema,
Eric J. Hallman,
Jack O. Burns,
Patrick M. Motl
Abstract:
We use Enzo, a hybrid Eulerian AMR/N-body code including non-gravitational heating and cooling, to explore the morphology of the X-ray gas in clusters of galaxies and its evolution in current generation cosmological simulations. We employ and compare two observationally motivated structure measures: power ratios and centroid shift. Overall, the structure of our simulated clusters compares remark…
▽ More
We use Enzo, a hybrid Eulerian AMR/N-body code including non-gravitational heating and cooling, to explore the morphology of the X-ray gas in clusters of galaxies and its evolution in current generation cosmological simulations. We employ and compare two observationally motivated structure measures: power ratios and centroid shift. Overall, the structure of our simulated clusters compares remarkably well to low-redshift observations, although some differences remain that may point to incomplete gas physics. We find no dependence on cluster structure in the mass-observable scaling relations, T_X-M and Y_X-M, when using the true cluster masses. However, estimates of the total mass based on the assumption of hydrostatic equilibrium, as assumed in observational studies, are systematically low. We show that the hydrostatic mass bias strongly correlates with cluster structure and, more weakly, with cluster mass. When the hydrostatic masses are used, the mass-observable scaling relations and gas mass fractions depend significantly on cluster morphology, and the true relations are not recovered even if the most relaxed clusters are used. We show that cluster structure, via the power ratios, can be used to effectively correct the hydrostatic mass estimates and mass-scaling relations, suggesting that we can calibrate for this systematic effect in cosmological studies. Similar to observational studies, we find that cluster structure, particularly centroid shift, evolves with redshift. This evolution is mild but will lead to additional errors at high redshift. Projection along the line of sight leads to significant uncertainty in the structure of individual clusters: less than 50% of clusters which appear relaxed in projection based on our structure measures are truly relaxed.
△ Less
Submitted 23 January, 2008; v1 submitted 10 August, 2007;
originally announced August 2007.
-
The Beta-Model Problem: The Incompatibility of X-ray and Sunyaev-Zeldovich Effect Model Fitting for Galaxy Clusters
Authors:
Eric J. Hallman,
Jack O. Burns,
Patrick M. Motl,
Michael L. Norman
Abstract:
We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse effects resulting from use of X-ray fitted beta-model parameters with Sunyaev-Zeldovich effect (SZE) data. There is a fundamental incompatibility between model fits to X-ray surface brightness profiles and those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X-ra…
▽ More
We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse effects resulting from use of X-ray fitted beta-model parameters with Sunyaev-Zeldovich effect (SZE) data. There is a fundamental incompatibility between model fits to X-ray surface brightness profiles and those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X-ray parameters are often used in SZE analysis. We show that this leads to biased estimates of the integrated Compton y-parameter inside r_{500} calculated from clusters. We suggest a simple correction of the method, using a non-isothermal beta-model modified by a universal temperature profile, which brings these calculated quantities into closer agreement with the true values.
△ Less
Submitted 3 May, 2007;
originally announced May 2007.
-
Further Evidence for an Elliptical Instability in Rotating Fluid Bars and Ellipsoidal Stars
Authors:
Shangli Ou,
Joel E. Tohline,
Patrick M. Motl
Abstract:
Using a three-dimensional nonlinear hydrodynamic code, we examine the dynamical stability of more than twenty self-gravitating, compressible, ellipsoidal fluid configurations that initially have the same velocity structure as Riemann S-type ellipsoids. Our focus is on ``adjoint'' configurations, in which internal fluid motions dominate over the collective spin of the ellipsoidal figure; Dedekind…
▽ More
Using a three-dimensional nonlinear hydrodynamic code, we examine the dynamical stability of more than twenty self-gravitating, compressible, ellipsoidal fluid configurations that initially have the same velocity structure as Riemann S-type ellipsoids. Our focus is on ``adjoint'' configurations, in which internal fluid motions dominate over the collective spin of the ellipsoidal figure; Dedekind-like configurations are among this group. We find that, although some models are stable and some are moderately unstable, the majority are violently unstable toward the development of $m=1$, $m=3$, and higher-order azimuthal distortions that destroy the coherent, $m=2$ bar-like structure of the initial ellipsoidal configuration on a dynamical time scale.
The parameter regime over which our models are found to be unstable generally corresponds with the regime over which incompressible Riemann S-type ellipsoids have been found to be susceptible to an elliptical strain instability \citep{LL96}. We therefore suspect that an elliptical instability is responsible for the destruction of our compressible analogs of Riemann ellipsoids. The existence of the elliptical instability raises concerns regarding the final fate of neutron stars that encounter the secular bar-mode instability and regarding the spectrum of gravitational waves that will be radiated from such systems.
△ Less
Submitted 28 February, 2007;
originally announced February 2007.
-
The Stability of Double White Dwarf Binaries Undergoing Direct Impact Accretion
Authors:
Patrick M. Motl,
Juhan Frank,
Joel E. Tohline,
Mario C. R. D'Souza
Abstract:
We present numerical simulations of dynamically unstable mass transfer in a double white dwarf binary with initial mass ratio, q = 0.4. The binary components are approximated as polytropes of index n = 3/2 and the initially synchronously rotating, semi-detached equilibrium binary is evolved hydrodynamically with the gravitational potential being computed through the solution of Poisson's equatio…
▽ More
We present numerical simulations of dynamically unstable mass transfer in a double white dwarf binary with initial mass ratio, q = 0.4. The binary components are approximated as polytropes of index n = 3/2 and the initially synchronously rotating, semi-detached equilibrium binary is evolved hydrodynamically with the gravitational potential being computed through the solution of Poisson's equation. Upon initiating deep contact in our baseline simulation, the mass transfer rate grows by more than an order of magnitude over approximately ten orbits, as would be expected for dynamically unstable mass transfer. However, the mass transfer rate then reaches a peak value, the binary expands and the mass transfer event subsides. The binary must therefore have crossed the critical mass ratio for stability against dynamical mass transfer. Despite the initial loss of orbital angular momentum into the spin of the accreting star, we find that the accretor's spin saturates and angular momentum is returned to the orbit more efficiently than has been previously suspected for binaries in the direct impact accretion mode. To explore this surprising result, we directly measure the critical mass ratio for stability by imposing artificial angular momentum loss at various rates to drive the binary to an equilibrium mass transfer rate. For one of these driven evolutions, we attain equilibrium mass transfer and deduce that effectively q_crit has evolved to approximately 2/3. Despite the absence of a fully developed disk, tidal interactions appear effective in returning excess spin angular momentum to the orbit.
△ Less
Submitted 19 August, 2007; v1 submitted 14 February, 2007;
originally announced February 2007.
-
Beyond the Cool Core: The Formation of Cool Core Galaxy Clusters
Authors:
Jack O. Burns,
Eric J. Hallman,
Brennan Gantner,
Patrick M. Motl,
Michael L. Norman
Abstract:
Why do some clusters have cool cores while others do not? In this paper, cosmological simulations, including radiative cooling and heating, are used to examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. Numerical CC clusters at z=0 accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the…
▽ More
Why do some clusters have cool cores while others do not? In this paper, cosmological simulations, including radiative cooling and heating, are used to examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. Numerical CC clusters at z=0 accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CCs survived the collisions. By contrast, NCC clusters of similar mass experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevent CC re-formation. We discuss observational consequences.
△ Less
Submitted 27 October, 2006;
originally announced October 2006.
-
Morphology and Evolution of Simulated and Optical Clusters: A Comparative Analysis
Authors:
Nurur Rahman,
Janusz Krywult,
Patrick M. Motl,
Piotr Flin,
Sergei F. Shandarin
Abstract:
We have made a comparative study of morphological evolution in simulated DM halos and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters include star formation with supernovae feedback, radiative cooling, and simulation in the adiabatic limit at three different redshifts, z = 0.0, 0.10, and 0.25. The optical sample contains 208 ACO clusters within redshift,…
▽ More
We have made a comparative study of morphological evolution in simulated DM halos and X-ray brightness distribution, and in optical clusters. Samples of simulated clusters include star formation with supernovae feedback, radiative cooling, and simulation in the adiabatic limit at three different redshifts, z = 0.0, 0.10, and 0.25. The optical sample contains 208 ACO clusters within redshift, $z \leq 0.25$. Cluster morphology, within 0.5 and 1.0 h$^{-1}$ Mpc from cluster center, is quantified by multiplicity and ellipticity.
We find that the distribution of the dark matter halos in the adiabatic simulation appear to be more elongated than the galaxy clusters. Radiative cooling brings halo shapes in excellent agreement with observed clusters, however, cooling along with feedback mechanism make the halos more flattened.
Our results indicate relatively stronger structural evolution and more clumpy distributions in observed clusters than in the structure of simulated clusters, and slower increase in simulated cluster shapes compared to those in the observed one.
Within $z \leq 0.1$, we notice an interesting agreement in the shapes of clusters obtained from the cooling simulations and observation. We also notice that the different samples of observed clusters differ significantly in morphological evolution with redshift. We highlight a few possibilities responsible for the discrepancy in morphological evolution of simulated and observed clusters.
△ Less
Submitted 29 December, 2005;
originally announced December 2005.
-
Numerical Simulations of the Onset and Stability of Dynamical Mass Transfer in Binaries
Authors:
Mario C. R. D'Souza,
Patrick M. Motl,
Joel E. Tohline,
Juhan Frank
Abstract:
Hydrodynamical simulations of semi-detached, polytropic binary stars are presented in an effort to study the onset and stability of dynamical mass transfer events. Initial, synchronously rotating equilibrium models are constructed using a self-consistent-field technique and then evolved with an Eulerian hydrodynamics code in a fully self-consistent manner. We describe code improvements introduce…
▽ More
Hydrodynamical simulations of semi-detached, polytropic binary stars are presented in an effort to study the onset and stability of dynamical mass transfer events. Initial, synchronously rotating equilibrium models are constructed using a self-consistent-field technique and then evolved with an Eulerian hydrodynamics code in a fully self-consistent manner. We describe code improvements introduced over the past few years that permit us to follow dynamical mass-transfer events through more than 30 orbits. Mass-transfer evolutions are presented for two different initial configurations: A dynamically unstable binary with initial mass ratio (donor/accretor) $q_0 = 1.3$ that leads to a complete merger in $\sim 10$ orbits; and a double-degenerate binary with initial mass ratio $q_0 = 0.5$ that, after some initial unstable growth of mass transfer, tends to separate as the mass-transfer rate levels off.
△ Less
Submitted 5 December, 2005;
originally announced December 2005.
-
Ruling the Universe: An Improved Method for Measuring H_0 with Galaxy Clusters
Authors:
Eric J. Hallman,
Jack O. Burns,
Patrick M. Motl,
Michael L. Norman
Abstract:
Withdrawn at the current time due to errors.
Withdrawn at the current time due to errors.
△ Less
Submitted 7 July, 2006; v1 submitted 26 October, 2005;
originally announced October 2005.
-
Challenges for Precision Cosmology with X-ray and Sunyaev-Zeldovich Effect Gas Mass Measurements of Galaxy Clusters
Authors:
Eric J. Hallman,
Patrick M. Motl,
Jack O. Burns,
Michael L. Norman
Abstract:
We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numerically simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect (SZE), we reduce the observations to obtain measurements of the cluster gas mass. We are thus able to…
▽ More
We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numerically simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect (SZE), we reduce the observations to obtain measurements of the cluster gas mass. We are thus able to quantify the possible sources of uncertainty and systematic bias associated with the common simplifying assumptions used in reducing real cluster observations including isothermality and hydrostatic equilibrium. We find that intrinsic variations in clusters limit the precision of observational gas mass estimation to ~10% to 1 sigma confidence excluding instrumental effects. Gas mass estimates performed via all methods surprisingly show little or no trending in the scatter as a function of cluster redshift. For the full cluster sample, methods that use SZE profiles out to roughly the virial radius are the simplest, most accurate, and unbiased way to estimate cluster mass. X-ray methods are systematically more precise mass estimators than are SZE methods if merger and cool core systems are removed, but X-ray methods slightly overestimate (5-10%) the cluster gas mass on average. We find that cool core clusters in our samples are particularly poor candidates for observational mass estimation, even when excluding emission from the core region. The effects of cooling in the cluster gas alter the radial profile of the X-ray and SZE surface brightness outside the cool core region, leading to poor gas mass estimates in cool core clusters. Finally, we find that methods using a universal temperature profile estimate cluster masses to higher precision than those assuming isothermality.
△ Less
Submitted 24 April, 2006; v1 submitted 15 September, 2005;
originally announced September 2005.
-
The Integrated Sunyaev-Zeldovich Effect as the Superior Method for Measuring the Mass of Clusters of Galaxies
Authors:
Patrick M. Motl,
Eric J. Hallman,
Jack O. Burns,
Michael L. Norman
Abstract:
We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four different samples that differ only in their assumed baryonic physics. We show that the strength of the thermal SZE integrated over a significant fraction of the virialized region of the clusters is re…
▽ More
We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four different samples that differ only in their assumed baryonic physics. We show that the strength of the thermal SZE integrated over a significant fraction of the virialized region of the clusters is relatively insensitive to the detailed heating and cooling processes in the cores of clusters by demonstrating that the derived scaling relations are nearly identical between the four cluster samples considered. For our synthetic images, the central Comptonization parameter shows significant boosting during transient merging events, but the integrated SZE appears to be relatively insensitive to these events. Most importantly, the integrated SZE closely tracks the underlying cluster mass. Observations through the thermal SZE allow a strikingly accurate mass estimation from relatively simple measurements that do not require either parametric modeling or geometric deprojection and thus avoid assumptions regarding the physics of the ICM or the symmetry of the cluster. This result offers significant promise for precision cosmology using clusters of galaxies.
△ Less
Submitted 22 March, 2005; v1 submitted 10 February, 2005;
originally announced February 2005.
-
Morphology and Evolution in Galaxy Clusters I: Simulated Clusters in the Adiabatic limit and with Radiative Cooling
Authors:
Nurur Rahman,
Sergei F. Shandarin,
Patrick M. Motl,
Adrian L. Melott,
;
Abstract:
We have studied morphological evolution in clusters simulated in the adiabatic limit and with radiative cooling. Cluster morphology in the redshift range, $0 < z < 0.5$, is quantified by multiplicity and ellipticity. In terms of ellipticity, our result indicates slow evolution in cluster shapes compared to those observed in the X-ray and optical wavelengths. The result is consistent with Floor,…
▽ More
We have studied morphological evolution in clusters simulated in the adiabatic limit and with radiative cooling. Cluster morphology in the redshift range, $0 < z < 0.5$, is quantified by multiplicity and ellipticity. In terms of ellipticity, our result indicates slow evolution in cluster shapes compared to those observed in the X-ray and optical wavelengths. The result is consistent with Floor, Melott & Motl (2003). In terms of multiplicity, however, the result indicate relatively stronger evolution (compared to ellipticity but still weaker than observation) in the structure of simulated clusters suggesting that for comparative studies of simulation and observation, sub-structure measures are more sensitive than the shape measures. We highlight a few possibilities responsible for the discrepancy in the shape evolution of simulated and real clusters.
△ Less
Submitted 5 May, 2004;
originally announced May 2004.
-
On the Formation of Cool, Non-Flowing Cores in Galaxy Clusters via Hierarchical Mergers
Authors:
J. O. Burns,
P. M. Motl,
M. L. Norman,
G. L. Bryan
Abstract:
We present a new model for the creation of cool cores in rich galaxy clusters within a LambdaCDM cosmological framework using the results from high spatial dynamic range, adaptive mesh hydro/N-body simulations. It is proposed that cores of cool gas first form in subclusters and these subclusters merge to create rich clusters with cool, central X-Ray excesses. The rich cool clusters do not posses…
▽ More
We present a new model for the creation of cool cores in rich galaxy clusters within a LambdaCDM cosmological framework using the results from high spatial dynamic range, adaptive mesh hydro/N-body simulations. It is proposed that cores of cool gas first form in subclusters and these subclusters merge to create rich clusters with cool, central X-Ray excesses. The rich cool clusters do not possess ``cooling flows'' due to the presence of bulk velocities in the intracluster medium in excess of 1000 km/sec produced by on-going accretion of gas from supercluster filaments. This new model has several attractive features including the presence of substantial core substructure within the cool cores, and it predicts the appearance of cool bullets, cool fronts, and cool filaments all of which have been recently observed with X-Ray satellites. This hierarchical formation model is also consistent with the observation that cool cores in Abell clusters occur preferentially in dense supercluster environments. On the other hand, our simulations overproduce cool cores in virtually all of our numerical clusters, the central densities are high, and physical core temperatures are often below 1 keV (in contrast to recent observations). We will discuss additional preliminary simulations to ``soften'' the cool cores involving star formation and supernova feedback.
△ Less
Submitted 30 September, 2003;
originally announced September 2003.
-
The Impact of Star Formation on Cool Core Galaxy Clusters
Authors:
P. M. Motl,
J. O. Burns,
M. L. Norman,
G. L. Bryan
Abstract:
We present results from recent simulations of the formation and evolution of clusters of galaxies in a LambdaCDM cosmology. These simulations contain our most physically complete input physics to date including radiative cooling, star formation that transforms rapidly cooling material into aggregate star particles and we also model the thermal feedback from resulting supernovae in the star parti…
▽ More
We present results from recent simulations of the formation and evolution of clusters of galaxies in a LambdaCDM cosmology. These simulations contain our most physically complete input physics to date including radiative cooling, star formation that transforms rapidly cooling material into aggregate star particles and we also model the thermal feedback from resulting supernovae in the star particles. We use an adaptive mesh refinement (AMR) Eulerian hydrodynamics scheme to obtain very high spatial resolution (~ 2 kpc) in a computational volume 256 Mpc on a side with mass resolution for dark matter and star particles of ~ 10^8 M_solar. We examine in detail the appearance and evolution of the core region of our simulated clusters.
△ Less
Submitted 30 September, 2003;
originally announced September 2003.
-
Formation of Cool Cores in Galaxy Clusters via Hierarchical Mergers
Authors:
Patrick M. Motl,
Jack O. Burns,
Chris Loken,
Michael L. Norman,
Greg Bryan
Abstract:
We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the…
▽ More
We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the accretion of discrete, stable subclusters. Any ``cooling flow'' present is overwhelmed by the velocity field within the cluster - the bulk flow of gas through the cluster typically has speeds up to about 2,000 km s^-1 and significant rotation is frequently present in the cluster core. The inclusion of consistent initial cosmological conditions for the cluster within its surrounding supercluster environment is crucial when simulating the evolution of cool cores in rich galaxy clusters. This new model for the hierarchical assembly of cool gas naturally explains the high frequency of cool cores in rich galaxy clusters despite the fact that a majority of these clusters show evidence of substructure which is believed to arise from recent merger activity. Furthermore, our simulations generate complex cluster cores in concordance with recent X-ray observations of cool fronts, cool ``bullets'', and filaments in a number of galaxy clusters. Our simulations were computed with a coupled N-body, Eulerian, adaptive mesh refinement, hydrodynamics cosmology code that properly treats the effects of shocks and radiative cooling by the gas. We employ up to seven levels of refinement to attain a peak resolution of 15.6 h^-1 kpc within a volume 256 h-1 Mpc on a side and assume a standard LambdaCDM cosmology.
△ Less
Submitted 20 February, 2003;
originally announced February 2003.
-
Numerical Methods for the Simulation of Dynamical Mass Transfer in Binaries
Authors:
Patrick M Motl,
Joel E Tohline,
Juhan Frank
Abstract:
We describe computational tools that have been developed to simulate dynamical mass transfer in semi-detached, polytropic binaries that are initially executing synchronous rotation upon circular orbits. Initial equilibrium models are generated with a self-consistent field algorithm; models are then evolved in time with a parallel, explicit, Eulerian hydrodynamics code with no assumptions made ab…
▽ More
We describe computational tools that have been developed to simulate dynamical mass transfer in semi-detached, polytropic binaries that are initially executing synchronous rotation upon circular orbits. Initial equilibrium models are generated with a self-consistent field algorithm; models are then evolved in time with a parallel, explicit, Eulerian hydrodynamics code with no assumptions made about the symmetry of the system. Poisson's equation is solved along with the equations of ideal fluid mechanics to allow us to treat the nonlinear tidal distortion of the components in a fully self-consistent manner. We present results from several standard numerical experiments that have been conducted to assess the general viability and validity of our tools, and from benchmark simulations that follow the evolution of two detached systems through five full orbits (up to approximately 90 stellar dynamical times). These benchmark runs allow us to gauge the level of quantitative accuracy with which simulations of semi-detached systems can be performed using presently available computing resources. We find that we should be able to resolve mass transfer at levels $\dot{M} / M > few x 10^-5$ per orbit through approximately 20 orbits with each orbit taking about 30 hours of computing time on parallel computing platforms.
△ Less
Submitted 17 July, 2001;
originally announced July 2001.
