-
The Prevalence of Resonance Among Young, Close-in Planets
Authors:
Fei Dai,
Max Goldberg,
Konstantin Batygin,
Jennifer van Saders,
Eugene Chiang,
Nick Choksi,
Rixin Li,
Erik A. Petigura,
Gregory J. Gilbert,
Sarah C. Millholland,
Yuan-Zhe Dai,
Luke Bouma,
Lauren M. Weiss,
Joshua N. Winn
Abstract:
Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common…
▽ More
Multiple planets undergoing disk migration may be captured into a chain of mean-motion resonances with the innermost planet parked near the disk's inner edge. Subsequent dynamical evolution may disrupt these resonances, leading to the non-resonant configurations typically observed among {\it Kepler} planets that are Gyrs old. In this scenario, resonant configurations are expected to be more common in younger systems. This prediction can now be tested, thanks to recent discoveries of young planets, particularly those in stellar clusters, by NASA's {\it TESS} mission. We divided the known planetary systems into three age groups: young ($<$100-Myr-old), adolescent (0.1-1-Gyr-old), and mature ($>1$-Gyr-old). The fraction of neighboring planet pairs having period ratios within a few percent of a first-order commensurability (e.g.~4:3, 3:2, or 2:1) is 70$\pm$15\% for young pairs, 24$\pm$8\% for adolescent pairs, and 15$\pm$2\% for mature pairs. The fraction of systems with at least one nearly commensurable pair (either first or second-order) is 86$\pm13$\% among young systems, 38$\pm12$\% for adolescent systems, and 23$\pm3$\% for mature systems. First-order commensurabilities prevail across all age groups, with an admixture of second-order commensurabilities. Commensurabilities are more common in systems with high planet multiplicity and low mutual inclinations. Observed period ratios often deviate from perfect commensurability by $\sim$1\% even among young planets, too large to be explained by resonant repulsion with equilibrium eccentricity tides. We also find that super-Earths in the radius gap ($1.5-1.9R_\oplus$) are less likely to be near-resonant (11.9$\pm2.0\%$) compared to Earth-sized planets ($R_p<1R_\oplus$; 25.3$\pm4.4\%$) or mini-Neptunes ($1.9R_\oplus \leq R_p<2.5R_\oplus$; 14.4$\pm1.8\%$).
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Spectral Energy Distributions of Disc-Embedded Accreting Protoplanets
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
Many dozens of protoplanetary discs show signatures of sculpting by planets. To help find these protoplanets by direct imaging, we compute their broadband spectral energy distributions, which overlap with the JWST (James Webb Space Telescope) and ALMA (Atacama Large Millimeter Array) passbands. We consider circumplanetary spherical envelopes, and combinations of circumplanetary discs and envelopes…
▽ More
Many dozens of protoplanetary discs show signatures of sculpting by planets. To help find these protoplanets by direct imaging, we compute their broadband spectral energy distributions, which overlap with the JWST (James Webb Space Telescope) and ALMA (Atacama Large Millimeter Array) passbands. We consider circumplanetary spherical envelopes, and combinations of circumplanetary discs and envelopes, heated by protoplanet accretion and external irradiation. Searches with JWST's NIRCam (Near-Infrared Camera) and the blue portion of MIRI (Mid-Infrared Instrument) are most promising since $\sim$300--1000 K protoplanets outshine their $\sim$20--50 K circumstellar environs at wavelengths of $\sim$2--20 $μ$m. Detection is easier if circumplanetary dust settles into discs (more likely for more massive planets) or is less abundant per unit mass gas (because of grain growth or filtration). At wavelengths longer than 20 $μ$m, circumplanetary material is difficult to see against the circumstellar disc's surface layers that directly absorb starlight and reprocess it to the far-infrared. Such contaminating circumstellar emission can be serious even within the evacuated gaps observed by ALMA. Only in strongly depleted regions, like the cavity of the transitional disc PDS 70 where two protoplanets have been confirmed, may long-wavelength windows open for protoplanet study. We compile a list of candidate protoplanets and identify those with potentially the highest accretion luminosities, all peaking in the near-infrared.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
The maximum accretion rate of a protoplanet: how fast can runaway be?
Authors:
Nick Choksi,
Eugene Chiang,
Jeffrey Fung,
Zhaohuan Zhu
Abstract:
The hunt is on for dozens of protoplanets hypothesised to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations we calculate maximum gas accretion rates for planet masses $M_{\rm p}$ from 1$\,M_{\oplus}$ to $10\,M_{\rm J}$. When the planet is…
▽ More
The hunt is on for dozens of protoplanets hypothesised to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations we calculate maximum gas accretion rates for planet masses $M_{\rm p}$ from 1$\,M_{\oplus}$ to $10\,M_{\rm J}$. When the planet is small enough that its sphere of influence is fully embedded in the disc, with a Bondi radius $r_{\rm Bondi}$ smaller than the disc's scale height $H_{\rm p}$ -- such planets have thermal mass parameters $q_{\rm th} \equiv (M_{\rm p}/M_{\star}) / (H_{\rm p}/R_{\rm p})^3 \lesssim 0.3$, for host stellar mass $M_{\star}$ and orbital radius $R_{\rm p}$ -- the maximum accretion rate follows a Bondi scaling, with $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^2 / (H_{\rm p}/R_{\rm p})^3$ for ambient disc density $ρ_{\rm g}$. For more massive planets with $0.3 \lesssim q_{\rm th} \lesssim 10$, the Hill sphere replaces the Bondi sphere as the gravitational sphere of influence, and $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^1$, with no dependence on $H_{\rm p}/R_{\rm p}$. In the strongly superthermal limit when $q_{\rm th} \gtrsim 10$, the Hill sphere pops well out of the disc, and $\max \dot{M}_{\rm p} \propto ρ_{\rm g} M_{\rm p}^{2/3} (H_{\rm p}/R_{\rm p})^1$. Applied to the two confirmed protoplanets PDS 70b and c, our numerically calibrated maximum accretion rates imply their Jupiter-like masses may increase by up to a factor of $\sim$2 before their parent disc dissipates.
△ Less
Submitted 25 July, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Exciting the TTV Phases of Resonant Sub-Neptunes
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
There are excesses of sub-Neptunes just wide of period commensurabilities like the 3:2 and 2:1, and corresponding deficits narrow of them. Any theory that explains this period ratio structure must also explain the strong transit timing variations (TTVs) observed near resonance. Besides an amplitude and a period, a sinusoidal TTV has a phase. Often overlooked, TTV phases are effectively integration…
▽ More
There are excesses of sub-Neptunes just wide of period commensurabilities like the 3:2 and 2:1, and corresponding deficits narrow of them. Any theory that explains this period ratio structure must also explain the strong transit timing variations (TTVs) observed near resonance. Besides an amplitude and a period, a sinusoidal TTV has a phase. Often overlooked, TTV phases are effectively integration constants, encoding information about initial conditions or the environment. Many TTVs near resonance exhibit non-zero phases. This observation is surprising because dissipative processes that capture planets into resonance also damp TTV phases to zero. We show how both the period ratio structure and the non-zero TTV phases can be reproduced if pairs of sub-Neptunes capture into resonance in a gas disc while accompanied by a third eccentric non-resonant body. Convergent migration and eccentricity damping by the disc drives pairs to orbital period ratios wide of commensurability; then, after the disc clears, secular forcing by the third body phase-shifts the TTVs. The scenario predicts that resonant planets are apsidally aligned and possess eccentricities up to an order of magnitude larger than previously thought.
△ Less
Submitted 19 March, 2023; v1 submitted 28 November, 2022;
originally announced November 2022.
-
An Implementation of Stochastic Forces for the N-body code REBOUND
Authors:
Hanno Rein,
Nick Choksi
Abstract:
We describe the implementation of a new module which can be used to simulate physical systems in which the motion of particles is affected by stochastic forces. Such forces are expected to be present in turbulent circumstellar disks or remnant planetesimal disks. Our implementation offers a convenient way to generate correlated noise with a user-specified amplitude and auto-correlation time for ea…
▽ More
We describe the implementation of a new module which can be used to simulate physical systems in which the motion of particles is affected by stochastic forces. Such forces are expected to be present in turbulent circumstellar disks or remnant planetesimal disks. Our implementation offers a convenient way to generate correlated noise with a user-specified amplitude and auto-correlation time for each particle. The module has minimal memory requirements and is freely available within the REBOUNDx additional effects library.
△ Less
Submitted 13 May, 2022;
originally announced May 2022.
-
The Lick Observatory Supernova Search follow-up program: photometry data release of 70 stripped-envelope supernovae
Authors:
WeiKang Zheng,
Benjamin E. Stahl,
Thomas de Jaeger,
Alexei V. Filippenko,
Shan-Qin Wang,
Wen-Pei Gan,
Thomas G. Brink,
Ivan Altunin,
Raphael Baer-Way,
Andrew Bigley,
Kyle Blanchard,
Peter K. Blanchard,
James Bradley,
Samantha K. Cargill,
Chadwick Casper,
Teagan Chapman,
Vidhi Chander,
Sanyum Channa,
Byung Yun Choi,
Nick Choksi,
Matthew Chu,
Kelsey I. Clubb,
Daniel P. Cohen,
Paul A. Dalba,
Asia deGraw
, et al. (63 additional authors not shown)
Abstract:
We present BVRI and unfiltered Clear light curves of 70 stripped-envelope supernovae (SESNe), observed between 2003 and 2020, from the Lick Observatory Supernova Search (LOSS) follow-up program. Our SESN sample consists of 19 spectroscopically normal SNe~Ib, two peculiar SNe Ib, six SN Ibn, 14 normal SNe Ic, one peculiar SN Ic, ten SNe Ic-BL, 15 SNe IIb, one ambiguous SN IIb/Ib/c, and two superlum…
▽ More
We present BVRI and unfiltered Clear light curves of 70 stripped-envelope supernovae (SESNe), observed between 2003 and 2020, from the Lick Observatory Supernova Search (LOSS) follow-up program. Our SESN sample consists of 19 spectroscopically normal SNe~Ib, two peculiar SNe Ib, six SN Ibn, 14 normal SNe Ic, one peculiar SN Ic, ten SNe Ic-BL, 15 SNe IIb, one ambiguous SN IIb/Ib/c, and two superluminous SNe. Our follow-up photometry has (on a per-SN basis) a mean coverage of 81 photometric points (median of 58 points) and a mean cadence of 3.6d (median of 1.2d). From our full sample, a subset of 38 SNe have pre-maximum coverage in at least one passband, allowing for the peak brightness of each SN in this subset to be quantitatively determined. We describe our data collection and processing techniques, with emphasis toward our automated photometry pipeline, from which we derive publicly available data products to enable and encourage further study by the community. Using these data products, we derive host-galaxy extinction values through the empirical colour evolution relationship and, for the first time, produce accurate rise-time measurements for a large sample of SESNe in both optical and infrared passbands. By modeling multiband light curves, we find that SNe Ic tend to have lower ejecta masses and lower ejecta velocities than SNe~Ib and IIb, but higher $^{56}$Ni masses.
△ Less
Submitted 10 March, 2022;
originally announced March 2022.
-
Testing planet formation from the ultraviolet to the millimeter
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities $Σ_{\rm gas}$ from C$^{18}$O observations, and planet masses $M_{\rm p}$ from simulations fitted to obse…
▽ More
Gaps imaged in protoplanetary discs are suspected to be opened by planets. We compute the present-day mass accretion rates $\dot{M}_{\rm p}$ of seven hypothesized gap-embedded planets, plus the two confirmed planets in the PDS 70 disc. The accretion rates are based on disc gas surface densities $Σ_{\rm gas}$ from C$^{18}$O observations, and planet masses $M_{\rm p}$ from simulations fitted to observed gaps. Assuming accretion is Bondi-like, we find in eight out of nine cases that $\dot{M}_{\rm p}$ is consistent with the time-averaged value given by the current planet mass and system age, $M_{\rm p}/t_{\rm age}$. As system ages are comparable to circumstellar disc lifetimes, these gap-opening planets may be undergoing their last mass doublings, reaching final masses of $M_{\rm p} \sim 10-10^2 \, M_\oplus$ for the non-PDS 70 planets, and $M_{\rm p} \sim 1-10 \, M_{\rm J}$ for the PDS 70 planets. For another fifteen gaps without C$^{18}$O data, we predict $Σ_{\rm gas}$ by assuming their planets are accreting at their time-averaged $\dot{M}_{\rm p}$. Bondi accretion rates for PDS 70b and c are orders of magnitude higher than accretion rates implied by measured U-band and H$α$ fluxes, suggesting most of the accretion shock luminosity emerges in as yet unobserved wavebands, or that the planets are surrounded by dusty, highly extincting, quasi-spherical circumplanetary envelopes. Thermal emission from such envelopes or from circumplanetary discs, on Hill sphere scales, peaks at wavelengths in the mid-to-far-infrared and can reproduce observed mm-wave excesses.
△ Less
Submitted 29 November, 2021; v1 submitted 30 September, 2021;
originally announced October 2021.
-
Chondrules from high-velocity collisions: thermal histories and the agglomeration problem
Authors:
Nick Choksi,
Eugene Chiang,
Harold C. Connolly Jr.,
Zack Gainsforth,
Andrew J. Westphal
Abstract:
We assess whether chondrules, once-molten mm-sized spheres filling the oldest meteorites, could have formed from super-km/s collisions between planetesimals in the solar nebula. High-velocity collisions release hot and dense clouds of silicate vapor which entrain and heat chondrule precursors. Thermal histories of CB chondrules are reproduced for colliding bodies $\sim$10--100 km in radius. The sl…
▽ More
We assess whether chondrules, once-molten mm-sized spheres filling the oldest meteorites, could have formed from super-km/s collisions between planetesimals in the solar nebula. High-velocity collisions release hot and dense clouds of silicate vapor which entrain and heat chondrule precursors. Thermal histories of CB chondrules are reproduced for colliding bodies $\sim$10--100 km in radius. The slower cooling rates of non-CB, porphyritic chondrules point to colliders with radii $\gtrsim$ 500 km. How chondrules, collisionally dispersed into the nebula, agglomerated into meteorite parent bodies remains a mystery. The same orbital eccentricities and inclinations that enable energetic collisions prevent planetesimals from re-accreting chondrules efficiently and without damage; thus the sedimentary laminations of the CB/CH chondrite Isheyevo are hard to explain by direct fallback of collisional ejecta. At the same time, planetesimal surfaces may be littered with the shattered remains of chondrules. The micron-sized igneous particles recovered from comet 81P/Wild-2 may have originated from in-situ collisions and subsequent accretion in the proto-Kuiper belt, obviating the need to transport igneous solids across the nebula. Asteroid sample returns from Hayabusa2 and OSIRIS-REx may similarly contain chondrule fragments.
△ Less
Submitted 28 March, 2021; v1 submitted 21 September, 2020;
originally announced September 2020.
-
Sub-Neptune Formation: The View from Resonant Planets
Authors:
Nick Choksi,
Eugene Chiang
Abstract:
The orbital period ratios of neighbouring sub-Neptunes are distributed asymmetrically near first-order resonances. There are deficits of systems---"troughs" in the period ratio histogram---just short of commensurability, and excesses---"peaks"---just wide of it. We reproduce quantitatively the strongest peak-trough asymmetries, near the 3:2 and 2:1 resonances, using dissipative interactions betwee…
▽ More
The orbital period ratios of neighbouring sub-Neptunes are distributed asymmetrically near first-order resonances. There are deficits of systems---"troughs" in the period ratio histogram---just short of commensurability, and excesses---"peaks"---just wide of it. We reproduce quantitatively the strongest peak-trough asymmetries, near the 3:2 and 2:1 resonances, using dissipative interactions between planets and their natal discs. Disc eccentricity damping captures bodies into resonance and clears the trough, and when combined with disc-driven convergent migration, draws planets initially wide of commensurability into the peak. The migration implied by the magnitude of the peak is modest; reductions in orbital period are $\sim$10\%, supporting the view that sub-Neptunes complete their formation more-or-less in situ. Once captured into resonance, sub-Neptunes of typical mass $\sim$$5$--$15 M_{\oplus}$ stay captured (contrary to an earlier claim), as they are immune to the overstability that afflicts lower mass planets. Driving the limited, short-scale migration is a gas disc depleted in mass relative to a solar-composition disc by 3--5 orders of magnitude. Such gas-poor but not gas-empty environments are quantitatively consistent with sub-Neptune core formation by giant impacts (and not, e.g., pebble accretion). While disc-planet interactions at the close of the planet formation era adequately explain the 3:2 and 2:1 asymmetries at periods $\gtrsim$ $5$--$15$ days, subsequent modification by stellar tides appears necessary at shorter periods, particularly for the 2:1.
△ Less
Submitted 20 May, 2020; v1 submitted 6 March, 2020;
originally announced March 2020.
-
On the initial mass-radius relation of stellar clusters
Authors:
Nick Choksi,
J. M. Diederik Kruijssen
Abstract:
Young stellar clusters across nearly five orders of magnitude in mass appear to follow a power-law mass-radius relationship (MRR), $R_{\star} \propto M_{\star}^α$, with $α\approx 0.2 - 0.33$. We develop a simple analytic model for the cluster mass-radius relation. We consider a galaxy disc in hydrostatic equilibrium, which hosts a population of molecular clouds that fragment into clumps undergoing…
▽ More
Young stellar clusters across nearly five orders of magnitude in mass appear to follow a power-law mass-radius relationship (MRR), $R_{\star} \propto M_{\star}^α$, with $α\approx 0.2 - 0.33$. We develop a simple analytic model for the cluster mass-radius relation. We consider a galaxy disc in hydrostatic equilibrium, which hosts a population of molecular clouds that fragment into clumps undergoing cluster formation and feedback-driven expansion. The model predicts a mass-radius relation of $R_{\star} \propto M_{\star}^{1/2}$ and a dependence on the kpc-scale gas surface density $R_{\star} \propto Σ_{\rm g}^{-1/2}$, which results from the formation of more compact clouds (and cluster-forming clumps within) at higher gas surface densities. This environmental dependence implies that the high-pressure environments in which the most massive clusters can form also induce the formation of clusters with the smallest radii, thereby shallowing the observed MRR at high-masses towards the observed $R_{\star} \propto M_{\star}^{1/3}$. At low cluster masses, relaxation-driven expansion induces a similar shallowing of the MRR. We combine our predicted MRR with a simple population synthesis model and apply it to a variety of star-forming environments, finding good agreement. Our model predicts that the high-pressure formation environments of globular clusters at high redshift naturally led to the formation of clusters that are considerably more compact than those in the local Universe, thereby increasing their resilience to tidal shock-driven disruption and contributing to their survival until the present day.
△ Less
Submitted 2 September, 2021; v1 submitted 11 December, 2019;
originally announced December 2019.
-
Lick Observatory Supernova Search Follow-Up Program: Photometry Data Release of 93 Type Ia Supernovae
Authors:
Benjamin E. Stahl,
WeiKang Zheng,
Thomas de Jaeger,
Alexei V. Filippenko,
Andrew Bigley,
Kyle Blanchard,
Peter K. Blanchard,
Thomas G. Brink,
Samantha K. Cargill,
Chadwick Casper,
Sanyum Channa,
Byung Yun Choi,
Nick Choksi,
Jason Chu,
Kelsey I. Clubb,
Daniel P. Cohen,
Michael Ellison,
Edward Falcon,
Pegah Fazeli,
Kiera Fuller,
Mohan Ganeshalingam,
Elinor L. Gates,
Carolina Gould,
Goni Halevi,
Kevin T. Hayakawa
, et al. (30 additional authors not shown)
Abstract:
We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar event…
▽ More
We present BVRI and unfiltered light curves of 93 Type Ia supernovae (SNe Ia) from the Lick Observatory Supernova Search (LOSS) follow-up program conducted between 2005 and 2018. Our sample consists of 78 spectroscopically normal SNe Ia, with the remainder divided between distinct subclasses (three SN 1991bg-like, three SN 1991T-like, four SNe Iax, two peculiar, and three super-Chandrasekhar events), and has a median redshift of 0.0192. The SNe in our sample have a median coverage of 16 photometric epochs at a cadence of 5.4 days, and the median first observed epoch is ~4.6 days before maximum B-band light. We describe how the SNe in our sample are discovered, observed, and processed, and we compare the results from our newly developed automated photometry pipeline to those from the previous processing pipeline used by LOSS. After investigating potential biases, we derive a final systematic uncertainty of 0.03 mag in BVRI for our dataset. We perform an analysis of our light curves with particular focus on using template fitting to measure the parameters that are useful in standardising SNe Ia as distance indicators. All of the data are available to the community, and we encourage future studies to incorporate our light curves in their analyses.
△ Less
Submitted 24 September, 2019;
originally announced September 2019.
-
Origins of scaling relations of globular cluster systems
Authors:
Nick Choksi,
Oleg Y. Gnedin
Abstract:
Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physica…
▽ More
Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physical origins and evolution of these scaling relations. The ratio of the combined mass in GCs $M_{\rm GC}$ to the host dark matter halo mass $M_h$ is nearly constant at all redshifts, but its normalization evolves by a factor of $\sim$10 from birth to $z=0$. The relation is steeper than linear at halo masses $M_h \lesssim 10^{11.5} M_{\odot}$, primarily due to non-linearity in the stellar mass-halo mass relation. The near constancy of the ratio $M_{\rm GC}/M_h$, combined with the shape of the stellar mass-halo mass relation, sets the characteristic $U-$shape of the GC specific frequency as a function of host galaxy mass. The contribution of accreted satellite galaxies to the buildup of GC systems is a strong function of the host galaxy mass, ranging from $\approx$0% at $M_h \approx 10^{11} M_{\odot}$ to 80% at $M_h \approx 10^{15} M_{\odot}$. The metal-poor clusters are significantly more likely to form ex-situ relative to the metal-rich clusters, but a substantial fraction of metal-poor clusters still form in-situ in lower mass galaxies. Similarly, the fraction of red clusters increases from $\approx 10$% at $M_h = 10^{11} M_{\odot}$ to $\approx 60$% at $M_h \approx 10^{13} M_{\odot}$, and flattens at higher $M_h$. Clusters formation occurs essentially continuously at high redshift, while at low redshift galactic mergers become increasingly important for cluster formation.
△ Less
Submitted 13 May, 2019;
originally announced May 2019.
-
Formation of Globular Cluster Systems II: Impact of the Cutoff of the Cluster Initial Mass Function
Authors:
Nick Choksi,
Oleg Y. Gnedin
Abstract:
Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to…
▽ More
Observations of young star clusters reveal that the high-mass end of the cluster initial mass function (CIMF) deviates from a pure power-law and instead truncates exponentially. We investigate the effects of this truncation on the formation of globular cluster (GC) systems by updating our analytic model for cluster formation and evolution, which is based on dark matter halo merger trees coupled to empirical galactic scaling relations, and has been shown in previous work to match a wide array of observational data. The cutoff masses of $M_c=10^{6.5} M_{\odot}$ or $10^{7}M_{\odot}$ match many scaling relations: between the GC system mass and host halo mass, between the average metallicity of the GC system and host halo mass, and the distribution of cluster masses. This range of $M_c$ agrees with indirect measurements from extragalactic GC systems. Models with $M_c<10^{6.5}M_{\odot}$ cannot reproduce the observed GC metallicity and mass distributions in massive galaxies. The slope of the mass-metallicity relation for metal-poor clusters (blue tilt) for all $M_c$ models is consistent with observations within their errors, when measured using the same method. We introduce an alternative, more robust fitting method, which reveals a trend of increasing tilt slope for lower $M_c$. In our model the blue tilt arises because the metal-poor clusters form in relatively low-mass galaxies which lack sufficient cold gas to sample the CIMF at highest masses. Massive blue clusters form in progressively more massive galaxies and inherit their higher metallicity. The metal-rich clusters do not exhibit such a tilt because they form in significantly more massive galaxies, which have enough cold gas to fully sample the CIMF.
△ Less
Submitted 20 March, 2019; v1 submitted 3 October, 2018;
originally announced October 2018.
-
The star clusters that make black hole binaries across cosmic time
Authors:
Nick Choksi,
Marta Volonteri,
Monica Colpi,
Oleg Y. Gnedin,
Hui Li
Abstract:
We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble…
▽ More
We explore the properties of dense star clusters that are likely to be the nurseries of stellar black holes pairing in close binaries. We combine a cosmological model of globular cluster formation with analytic prescriptions for the dynamical assembly and evolution of black hole binaries to constrain which types of clusters are most likely to form binaries tight enough to coalesce within a Hubble time. We find that black hole binaries which are ejected and later merge ex-situ form in clusters of a characteristic mass $M_{\rm cl} \sim 10^{5.3}M_{\odot}$, whereas binaries which merge in-situ form in more massive clusters, $M_{\rm cl} \sim 10^{5.7}M_{\odot}$. The clusters which dominate the production of black hole binaries are similar in age and metallicity to the entire population. Finally, we estimate an approximate cosmic black hole merger rate of dynamically assembled binaries using the mean black hole mass for each cluster given its metallicity. We find an intrinsic rate of $\sim 6\,\mathrm{Gpc\,yr^{-1}}$ at $z=0$, a weakly increasing merger rate out to $z=1.5$, and then a decrease out to $z=4$. Our results can be used to provide a cosmological context and choose initial conditions in numerical studies of black hole binaries assembled in star clusters.
△ Less
Submitted 4 March, 2019; v1 submitted 4 September, 2018;
originally announced September 2018.
-
The formation and hierarchical assembly of globular cluster populations
Authors:
Kareem El-Badry,
Eliot Quataert,
Daniel R. Weisz,
Nick Choksi,
Michael Boylan-Kolchin
Abstract:
We use a semi-analytic model for globular cluster (GC) formation built on dark matter merger trees to explore the relative role of formation physics and hierarchical assembly in determining the properties of GC populations. Many previous works have argued that the observed linear relation between total GC mass and halo mass points to a fundamental GC -- dark matter connection or indicates that GCs…
▽ More
We use a semi-analytic model for globular cluster (GC) formation built on dark matter merger trees to explore the relative role of formation physics and hierarchical assembly in determining the properties of GC populations. Many previous works have argued that the observed linear relation between total GC mass and halo mass points to a fundamental GC -- dark matter connection or indicates that GCs formed at very high redshift before feedback processes introduced nonlinearity in the baryon-to-dark matter mass relation. We demonstrate that at $M_{\rm vir}(z=0) \gtrsim 10^{11.5} M_{\odot}$, a constant ratio between halo mass and total GC mass is in fact an almost inevitable consequence of hierarchical assembly: by the central limit theorem, it is expected at $z=0$ independent of the GC-to-halo mass relation at the time of GC formation. The GC-to-halo mass relation at $M_{\rm vir}(z=0) < 10^{11.5} M_{\odot}$ is more sensitive to the details of the GC formation process. In our fiducial model, GC formation occurs in galaxies when the gas surface density exceeds a critical value. This model naturally predicts bimodal GC color distributions similar to those observed in nearby galaxies and reproduces the observed relation between GC system metallicity and halo mass. It predicts that the cosmic GC formation rate peaked at $z$ $\sim$ 4, too late for GCs to contribute significantly to the UV luminosity density during reionization.
△ Less
Submitted 2 November, 2018; v1 submitted 9 May, 2018;
originally announced May 2018.
-
Formation of Globular Cluster Systems: From Dwarf Galaxies to Giants
Authors:
Nick Choksi,
Oleg Gnedin,
Hui Li
Abstract:
Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich…
▽ More
Globular cluster (GC) systems around galaxies of a vast mass range show remarkably simple scaling relations. The combined mass of all GCs is a constant fraction of the total galaxy mass and the mean metallicity and metallicity dispersion of the GC system scale up weakly with galaxy mass. The metallicity of massive, metal-poor ("blue") clusters increases with cluster mass, while that of metal-rich ("red") clusters does not. A significant age-metallicity relation emerges from analysis of resolved stellar populations in Galactic GCs and unresolved populations in nearby galaxies. Remarkably, all these trends can be explained by a simple merger-based model developed in previous work and updated here using recent observations of galaxy scaling relations at high redshift. We show that the increasing dispersion of GC metallicity distributions with galaxy mass is a robust prediction of the model. It arises from more massive galaxies having more mergers that combine satellite GC systems. The average metallicity also increases by 0.6~dex over 3~dex in halo mass. The models show a non-linear trend between the GC system mass and host galaxy mass which is consistent with the data. The model does not consider GC self-enrichment, yet predicts a correlation between cluster mass and metallicity for massive blue clusters. The age-metallicity relation is another robust prediction of the model. Half of all clusters are predicted to form within the redshift range $5<z<2.3$, corresponding to ages of $10.8-12.5$~Gyr, in halos of masses $10^{11}-10^{12.5} M_{\odot}$.
△ Less
Submitted 23 August, 2018; v1 submitted 10 January, 2018;
originally announced January 2018.
-
Recoiling Supermassive Black Hole Escape Velocities from Dark Matter Halos
Authors:
Nick Choksi,
Peter Behroozi,
Marta Volonteri,
Raffaela Schneider,
Chung-Pei Ma,
Joseph Silk,
Benjamin Moster
Abstract:
We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $Λ$CDM halo accretion increases escape…
▽ More
We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $Λ$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.
△ Less
Submitted 15 July, 2018; v1 submitted 19 July, 2017;
originally announced July 2017.