ABSTRACT In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $$1{{\ \rm per\ cent}}$$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations. We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to $${\gt}40{{\ \rm per\ cent}}$$ more black holes reaching the centre of their host halo when DF is added.
more »
« less
A new discrete dynamical friction estimator based on N -body simulations
ABSTRACT A long-standing problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models based on the traditional Chandrasekhar DF formula have been proposed, yet they suffer from fundamental ambiguities in the definition of some terms in Chandrasekhar’s formula when applied to real galaxies, as well as difficulty in evaluating continuous quantities from (spatially) discrete simulation data. In this work, we present a new sub-grid DF estimator based on the discrete nature of N-body simulations, which also avoids the ambiguously defined quantities in Chandrasekhar’s formula. We test our estimator in the gizmo code and find that it agrees well with high-resolution simulations where DF is fully captured, with negligible additional computational cost. We also compare it with a Chandrasekhar estimator and discuss its applications in real galactic simulations.
more »
« less
- PAR ID:
- 10405744
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 519
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 5543 to 5553
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub-Chandrasekhar mass white dwarfs; double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope; a violent merger model of two white dwarfs; and deflagrations and delayed detonations in Chandrasekhar mass white dwarfs. We focus on the iron peak elements Mn, Zn, and Cu. To this end, we also briefly review the different burning regimes and production sites of these elements, as well as the results of abundance measurements and several galactic chemical evolution studies. We find that super-solar values of [Mn/Fe] are not restricted to Chandrasekhar mass explosion models. Scenarios including a helium detonation can significantly contribute to the production of Mn, in particular the models proposed for calcium-rich transients. Although Type Ia supernovae are often not accounted for as production sites of Zn and Cu, our models involving helium shell detonations can produce these elements in super-solar ratios relative to Fe. Our results suggest a re-consideration of Type Ia supernova yields in galactic chemical evolution models. A detailed comparison with observations can provide new insight into the progenitor and explosion channels of these events.more » « less
-
Abstract On the main sequence, the asteroseismic small frequency separationδν02between radial and quadrupolep-modes is customarily interpreted to be a direct diagnostic of internal structure. Such an interpretation is based on a well-known integral estimator relatingδν02to a radially averaged sound-speed gradient. However, this estimator fails, catastrophically, when evaluated on structural models of red giants: their small separations must therefore be interpreted differently. We derive a single expression that both reduces to the classical estimator when applied to main-sequence stellar models and reproduces the qualitative features of the small separation for stellar models of very evolved red giants. This expression indicates that the small separations of red giants scale primarily with their global seismic properties as , rather than being in any way sensitive to their internal structure. Departures from this asymptotic behavior, during the transition from the main-sequence to red giant regimes, have been recently reported in open-cluster Christensen–Dalsgaard (C-D) diagrams from K2 mission data. Investigating them in detail, we demonstrate that they occur when the convective envelope boundary passes a specific acoustic distance—roughly one-third of a wavelength at —from the center of the star, at which point radial modes become maximally sensitive to the position of the boundary. The shape of the corresponding features onϵpand C-D (orr02) diagrams may be useful in constraining the nature of convective boundary mixing in the context of undershooting beneath a convective envelope.more » « less
-
Walker, Gilbert (Ed.)Microplastics (MPs) and nanoplastics (NPs) in water pose a global threat to human health and the environment. To develop efficient removal strategies, it is crucial to understand how these particles behave as they aggregate. However, our knowledge of the process of aggregate formation from primary particles of different sizes is limited. In this study, we analyzed the growth kinetics and structures of aggregates formed by polystyrene MPs in mono- and bidisperse systems using in situ microscopy and image analysis. Our findings show that the scaling behavior of aggregate growth remains unaffected by the primary particle size distribution, but it does delay the onset of rapid aggregation. We also performed a structural analysis that reveals the power law dependence of aggregate fractal dimension (df) in both mono- and bidisperse systems, with mean df consistent with diffusion-limited cluster aggregation (DLCA) aggregates. Our results also suggest that the df of aggregates is insensitive to the shape anisotropy. We simulated molecular forces driving aggregation of polystyrene NPs of different sizes under high ionic strength conditions. These conditions represent salt concentration in ocean water and wastewater, where the DLVO theory does not apply. Our simulation results show that the aggregation tendency of the NPs increases with the ionic strength. The increase in the aggregation is caused by the depletion of clusters of ions from the NPs surface.more » « less
-
Abstract. Satellite cloud observations have become an indispensable tool for evaluatinggeneral circulation models (GCMs). To facilitate the satellite and GCMcomparisons, the CFMIP (Cloud Feedback Model Inter-comparison Project)Observation Simulator Package (COSP) has been developed and is nowincreasingly used in GCM evaluations. Real-world clouds and precipitation canhave significant sub-grid variations, which, however, are often ignored oroversimplified in the COSP simulation. In this study, we use COSP cloudsimulations from the Super-Parameterized Community Atmosphere Model (SPCAM5)and satellite observations from the Moderate Resolution ImagingSpectroradiometer (MODIS) and CloudSat to demonstrate the importance ofconsidering the sub-grid variability of cloud and precipitation when usingthe COSP to evaluate GCM simulations. We carry out two sensitivity tests:SPCAM5 COSP and SPCAM5-Homogeneous COSP. In the SPCAM5 COSP run, the sub-gridcloud and precipitation properties from the embeddedcloud-resolving model (CRM) of SPCAM5 are used to drive the COSP simulation, while inthe SPCAM5-Homogeneous COSP run only grid-mean cloud and precipitationproperties (i.e., no sub-grid variations) are given to the COSP. We find thatthe warm rain signatures in the SPCAM5 COSP run agree with the MODIS andCloudSat observations quite well. In contrast, the SPCAM5-Homogeneous COSPrun which ignores the sub-grid cloud variations substantially overestimatesthe radar reflectivity and probability of precipitation compared to thesatellite observations, as well as the results from the SPCAM5 COSP run. Thesignificant differences between the two COSP runs demonstrate that it isimportant to take into account the sub-grid variations of cloud andprecipitation when using COSP to evaluate the GCM to avoid confusing andmisleading results.more » « less
An official website of the United States government

