Source author record

Paul M. Ricker

Paul M. Ricker appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Flash-X, a multiphysics simulation software instrument

Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Flash-X is meant primarily for solving Eulerian formulations of applications with compressible and/or incompressible reactive flows. It also has a built-in, versatile Lagrangian framework that can be used in many different ways, including implementing tracers, particle-in-cell simulations, and immersed boundary methods.

preprint2020arXiv

Common Envelope Shaping of Planetary Nebulae. II. Magnetic Solutions and Self-Collimated Outflows

Magnetic fields of order $10^1-10^2$ gauss that are present in the envelopes of red giant stars are ejected in common envelope scenarios. These fields could be responsible for the launching of magnetically driven winds in proto-planetary nebulae. Using 2D simulations of magnetized winds interacting with an envelope drawn from a 3D simulation of the common envelope phase, we study the confinement, heating, and magnetic field development of post-common envelope winds. We find that the ejected magnetic field can be enhanced via compression by factors up to $\sim 10^4$ in circumbinary disks during the self-regulated phases. We find values for the kinetic energy of the order of $10^{46}$ erg that explain the large values inferred in proto-planetary nebula outflows. We show that the interaction of the formed circumbinary disk with a spherical, stellar wind produces a "tapered" flow that is almost indistinguishable from an imposed tapered flow. This increases the uncertainty of the origin of proto-planetary nebula winds, which could be either stellar, circumstellar (stellar accretion disk), circumbinary (circumbinary accretion disk), or a combination of all three. Within this framework, a scenario for self-collimation of weakly magnetized winds is discussed, which can explain the two objects where the collimation process is observationally resolved, HD 101584 and Hen 3-1475. An explanation for the equatorial, molecular hydrogen emission in CRL 2688 is also presented.

preprint2018arXiv

Common Envelope Evolution of Massive Stars

The discovery via gravitational waves of binary black hole systems with total masses greater than $60M_\odot$ has raised interesting questions for stellar evolution theory. Among the most promising formation channels for these systems is one involving a common envelope binary containing a low metallicity, core helium burning star with mass $\sim 80-90M_\odot$ and a black hole with mass $\sim 30-40M_\odot$. For this channel to be viable, the common envelope binary must eject more than half the giant star's mass and reduce its orbital separation by as much as a factor of 80. We discuss issues faced in numerically simulating the common envelope evolution of such systems and present a 3D AMR simulation of the dynamical inspiral of a low-metallicity red supergiant with a massive black hole companion.

preprint2016arXiv

A hot big bang theory: magnetic fields and the early evolution of the protolunar disk

The leading theory for the formation of the Earth's moon invokes a collision between a Mars-sized body and the proto-Earth to produce a disk of orbiting material that later condenses to form the Moon. Here we study the early evolution of the protolunar disk. First, we show that the disk opacity is large and cooling is therefore inefficient (t_{cool} Ω>> 1). In this regime angular momentum transport in the disk leads to steady heating unless α< (t_{cool} Ω)^{-1} << 1. Following earlier work by Charnoz and Michaut, and Carballido et al., we show that once the disk is completely vaporized it is well coupled to the magnetic field. We consider a scenario in which turbulence driven by magnetic fields leads to a brief, hot phase where the disk is geometrically thick, wit h strong turbulent mixing. The disk cools by spreading until it decouples from the field. We point out that approximately half the accretion energy is dissipated in the boundary layer where the disk meets the Earth's surface. This creates high entropy material close to the Earth, driving convection and mixing. Finally, a hot, magnetized disk could drive bipolar outflows that remove mass and angular momentum from the Earth-Moon system.

preprint2015arXiv

Ram Pressure Stripping of Hot Coronal Gas from Group and Cluster Galaxies and the Detectability of Surviving X-ray Coronae

Ram pressure stripping can remove hot and cold gas from galaxies in the intracluster medium (ICM), as shown by observations of X-ray and HI galaxy wakes in nearby clusters of galaxies. However, ram pressure stripping, including pre-processing in group environments, does not remove all the hot coronal gas from cluster galaxies. Recent high-resolution Chandra observations have shown that $\sim 1 - 4$ kpc extended, hot galactic coronae are ubiquitous in group and cluster galaxies. To better understand this result, we simulate ram pressure stripping of a cosmologically motivated population of galaxies in isolated group and cluster environments. The galaxies and the host group and cluster are composed of collisionless dark matter and hot gas initially in hydrostatic equilibrium with the galaxy and host potentials. We show that the rate at which gas is lost depends on the galactic and host halo mass. Using synthetic X-ray observations, we evaluate the detectability of stripped galactic coronae in real observations by stacking images on the known galaxy centers. We find that coronal emission should be detected within $\sim 10$ arcsec, or $\sim 5$ kpc up to $\sim 2.3$ Gyr in the lowest (0.1 - 1.2 keV) energy band. Thus the presence of observed coronae in cluster galaxies significantly smaller than the hot X-ray halos of field galaxies indicates that at least some gas removal occurs within cluster environments for recently accreted galaxies. Finally, we evaluate the possibility that existing and future X-ray cluster catalogs can be used in combination with optical galaxy positions to detect galactic coronal emission via stacking analysis. We briefly discuss the effects of additional physical processes on coronal survival, and will address them in detail in future papers in this series.

preprint2015arXiv

Simulations of the symbiotic recurrent nova V407 Cyg. I. Accretion and shock evolutions

The shock interaction and evolution of nova ejecta with a wind from a red giant star in a symbiotic binary system are investigated via three-dimensional hydrodynamics simulations. We specifically model the March 2010 outburst of the symbiotic recurrent nova V407~Cygni from the quiescent phase to its eruption phase. The circumstellar density enhancement due to wind-white dwarf interaction is studied in detail. It is found that the density-enhancement efficiency depends on the ratio of the orbital speed to the red giant wind speed. Unlike another recurrent nova, RS~Ophiuchi, we do not observe a strong disk-like density enhancement, but instead observe an aspherical density distribution with $\sim 20\%$ higher density in the equatorial plane than at the poles. To model the 2010 outburst, we consider several physical parameters, including the red giant mass loss rate, nova eruption energy, and ejecta mass. A detailed study of the shock interaction and evolution reveals that the interaction of shocks with the red giant wind generates strong Rayleigh-Taylor instabilities. In addition, the presence of the companion and circumstellar density enhancement greatly alter the shock evolution during the nova phase. The ejecta speed after sweeping out most of the circumstellar medium decreases to $\sim 100-300$ km-s$^{-1}$, depending on model, which is consistent with the observed extended redward emission in [N~II] lines in April 2011.

preprint2014arXiv

Search for surviving companions in type Ia supernova remnants

The nature of the progenitor systems of type~Ia supernovae is still unclear. One way to distinguish between the single-degenerate scenario and double-degenerate scenario for their progenitors is to search for the surviving companions. Using a technique that couples the results from multi-dimensional hydrodynamics simulations with calculations of the structure and evolution of main-sequence- and helium-rich surviving companions, the color and magnitude of main-sequence- and helium-rich surviving companions are predicted as functions of time. The surviving companion candidates in Galactic type~Ia supernova remnants and nearby extragalactic type~Ia supernova remnants are discussed. We find that the maximum detectable distance of main-sequence surviving companions (helium-rich surviving companions) is $0.6-4$~Mpc ($0.4-16$~Mpc), if the apparent magnitude limit is 27 in the absence of extinction, suggesting that the Large and Small Magellanic Clouds and the Andromeda Galaxy are excellent environments in which to search for surviving companions. However, only five Ia~SNRs have been searched for surviving companions, showing little support for the standard channels in the singe-degenerate scenario. To better understand the progenitors of type Ia supernovae, we encourage the search for surviving companions in other nearby type Ia supernova remnants.

preprint2014arXiv

The Dynamical Origin of Early-Type Dwarfs in Galaxy Clusters: A Theoretical Investigation

Observations of early-type dwarf galaxies in clusters often show that cluster dwarf members have significantly higher velocities and less symmetric distributions than cluster giant ellipticals, suggesting that these dwarfs are recently accreted galaxies, possibly from an infalling group. We use a series of $N$-body simulations, exploring a parameter space of groups falling into clusters, to study the observed velocity distributions of the infall components along various lines of sight. We show that, as viewed along a line of sight parallel to the group's infall direction, there is a significant peculiar velocity boost during the pericentric passage of the group, and an increase in velocity dispersion that persists for many Gyr after the merger. The remnants of the infalling group, however, do not form a spatially distinct system -- consistent with recent observations of dwarf galaxies in the Virgo and Fornax clusters. This velocity signature is completely absent when viewed along a direction perpendicular to the merger. Additionally, the phase-space distribution of radial velocity along the infall direction versus cluster-centric radius reveals the separate dynamical evolution of the group's central core and outer halo, including the presence of infalling remnants outside the escape velocity envelope of the system. The distinct signature in velocity space of an infalling group's galaxies can therefore prove important in understanding the dynamical history of clusters and their dwarfs. Our results suggest that dwarf galaxies, being insensitive to dynamical friction, are excellent probes of their host clusters' dynamical histories.

preprint2013arXiv

Pre-Processing and Post-Processing in Group-Cluster Mergers

Galaxies in clusters are more likely to be of early type and to have lower star formation rates than galaxies in the field. Recent observations and simulations suggest that cluster galaxies may be `pre-processed' by group or filament environments and that galaxies that fall into a cluster as part of a larger group can stay coherent within the cluster for up to one orbital period (`post-processing'). We investigate these ideas by means of a cosmological $N$-body simulation and idealized $N$-body plus hydrodynamics simulations of a group-cluster merger. We find that group environments can contribute significantly to galaxy pre-processing by means of enhanced galaxy-galaxy merger rates, removal of galaxies' hot halo gas by ram pressure stripping, and tidal truncation of their galaxies. Tidal distortion of the group during infall does not contribute to pre-processing. Post-processing is also shown to be effective: galaxy-galaxy collisions are enhanced during a group's pericentric passage within a cluster, the merger shock enhances the ram pressure on group and cluster galaxies, and an increase in local density during the merger leads to greater galactic tidal truncation.

preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

preprint2012arXiv

An AMR Study of the Common Envelope Phase of Binary Evolution

The hydrodynamic evolution of the common envelope phase of a low mass binary composed of a 1.05 Msun red giant and a 0.6 Msun companion has been followed for five orbits of the system using a high resolution method in three spatial dimensions. During the rapid inspiral phase, the interaction of the companion with the red giant's extended atmosphere causes about 25% of the common envelope to be ejected from the system, with mass continuing to be lost at the end of the simulation at a rate ~ 2 Msun/yr. In the process the resulting loss of angular momentum and energy reduces the orbital separation by a factor of seven. After this inspiral phase the eccentricity of the orbit rapidly decreases with time. The gravitational drag dominates hydrodynamic drag at all times in the evolution, and the commonly-used Bondi-Hoyle-Lyttleton prescription for estimating the accretion rate onto the companion significantly overestimates the true rate. On scales comparable to the orbital separation, the gas flow in the orbital plane in the vicinity of the two cores is subsonic with the gas nearly corotating with the red giant core and circulating about the red giant companion. On larger scales, 90% of the outflow is contained within 30 degrees of the orbital plane, and the spiral shocks in this material leave an imprint on the density and velocity structure. Of the energy released by the inspiral of the cores, only about 25% goes toward ejection of the envelope.

preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

preprint2010arXiv

Impact of Type Ia Supernova Ejecta on a Helium-star Binary Companion

The impact of Type Ia supernova ejecta on a helium-star companion is investigated via high-resolution, two-dimensional hydrodynamic simulations. For a range of helium-star models and initial binary separations it is found that the mass unbound in the interaction, $δM_{\rm ub}$, is related to the initial binary separation, $a$, by a power law of the form $δM_{\rm ub} \propto a^{m}$. This power-law index is found to vary from -3.1 to -4.0, depending on the mass of the helium star. The small range of this index brackets values found previously for hydrogen-rich companions, suggesting that the dependence of the unbound mass on orbital separation is not strongly sensitive to the nature of the binary companion. The kick velocity is also related to the initial binary separation by a power law with an index in a range from -2.7 to -3.3, but the power-law index differs from those found in previous studies for hydrogen-rich companions. The space motion of the companion after the supernova is dominated by its orbital velocity in the pre-supernova binary system. The level of Ni/Fe contamination of the companion resulting from the passage of the supernova ejecta is difficult to estimate, but an upper limit on the mass of bound nickel is found to be $\sim 5\times 10^{-4}\ M_\odot$.

preprint2010arXiv

The Impact of Cluster Structure and Dynamical State on Scatter in the Sunyaev-Zel'dovich Flux-Mass Relation

Cosmological constraints from cluster surveys rely on accurate mass estimates from the mass-observable relations. In order to avoid systematic biases and reduce uncertainties, we study the form and physical origin of the intrinsic scatter about the mean Sunyaev-Zel'dovich (SZ) flux-mass relation using a hydrodynamical simulation of galaxy cluster formation. We examine the assumption of lognormal scatter and detect non-negligible positive skewness and kurtosis (> 0.5) for a wide range of limiting masses and redshifts. These higher-order moments should be included in the parametrization of scatter in order not to bias cosmological constraints. We investigate the sources of the scatter by correlating it with measures of cluster morphology, halo concentration, and dynamical state, and we quantify the individual contribution from each source. We find that statistically the impact of dynamical state is weak, so the selection bias due to mergers is negligible. On the other hand, there is a strong correlation between the scatter and halo concentration, which can be used to reduce the scatter significantly (from 12.07% to 7.34% or by ~40% for clusters at z = 0). We also show that a cross-calibration by combining information from X-ray followups can be used to reduce the scatter in the flux-mass relation and also identify outliers in both X-ray and SZ cluster surveys.