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.

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.

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.

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.

We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical equations of state and neutrino cooling. Our studies show such an instability develops generically and within a short dynamical time to strengths that are comparable or stronger than the $m=2$ mode which is the strongest during the early post-merger stage. We estimate the signal to noise ratio that might be obtained for the $m=1$ mode and discuss the prospects for observing this signal with available Earth-based detectors. Because the $m=1$ occurs at roughly half the frequency of the more powerful $m=2$ signal and because it can potentially be long-lived, targeted searches could be devised to observe it. We estimate that with constant amplitude direct detection of the mode could occur up to a distance of roughly $14\,\mathrm{Mpc}$ whereas a search triggered by the inspiral signal could extend this distance to roughly $100\,\mathrm{Mpc}$.

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

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.

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.

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.

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.

This white paper briefly describes the astrophysics of ultra-compact binaries, with emphasis of the challenges and opportunities in the next decade.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

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.

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.

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.

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.

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.

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.

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.

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.

Withdrawn at the current time due to errors.

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.

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.

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.

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.

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.

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.

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.

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.

We investigate the predicted present-day temperature profiles of the hot, X-ray emitting gas in galaxy clusters for two cosmological models - a current best-guess LCDM model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" of clusters are derived from high-resolution (15/h kpc) simulations which make use of a sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully captures the shocks which are essential for correctly modelling cluster temperatures. We show that the temperature structure on Mpc-scales is highly complex and non-isothermal. However, the temperature profiles of the simulated LCDM and SCDM clusters are remarkably similar and drop-off as $T +AFw-propto (1+-r/a_x)^{-+AFw-delta}$ where $a_x +AFw-sim r_{vir}/1.5$ and $+AFw-delta +AFw-sim 1.6$. This decrease is in good agreement with the observational results of Markevitch et al.(1998) but diverges, primarily in the innermost regions, from their fit which assumes a polytropic equation of state. Our result is also in good agreement with a recent sample of clusters observed by BeppoSAX though there is some indication of missing physics at small radii ($r<0.2 r_{vir}$). We discuss the interpretation of our results and make predictions for new x-ray observations that will extend to larger radii than previously possible. Finally, we show that, for $r>0.2 r_{vir}$, our universal temperature profile is consistent with our most recent simulations which include both radiative cooling and supernovae feedback.

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.