skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Understanding the Origin and Dynamical Evolution of the Unique Open Star Cluster Berkeley 20 Using FIRE Simulations
Open clusters (OCs) act as key probes that can be leveraged to constrain the formation and evolution of the Milky Way (MW)’s disk, as each has a unique chemical fingerprint and well-constrained age. Significant Galactic dynamic interactions can leave imprints on the orbital properties of OCs, allowing us to use the present-day properties of long-lived OCs to reconstruct the MW’s dynamic history. To explore these changes, we identify OC analogs in FIRE-2 simulations of MW-mass galaxies. For this work, we focus on one particular FIRE-2 OC, which we identify as an analog to the old, subsolar, distant, and high-Galactic-latitude MW OC, Berkeley 20. Our simulated OC resides ∼6 kpc from the galactic center and ultimately reaches a height Z max > 2 kpc from the galactic disk, similar to Berkeley 20. We trace the simulated cluster’s orbital and environmental history, identifying key perturbative episodes, including (1) an interaction with a gas overdensity in a spiral arm that prompts an outward migration event and (2) a substantial interaction with a Sagittarius Dwarf Spheroidal Galaxy–mass satellite that causes significant orbital modification. Our simulated OC shows significant resilience to disruption during both its outward migration and the satellite-driven heating event that causes subsequent inward migration. Ultimately, we find these two key processes—migration and satellite heating—are essential to include when assessing OC orbital dynamics in the era of Gaia.  more » « less
Award ID(s):
2045928
PAR ID:
10655946
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal Letters
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
995
Issue:
1
ISSN:
2041-8205
Page Range / eLocation ID:
L25
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present the discovery of Aquarius III, an ultra-faint Milky Way satellite galaxy identified in the second data release of the DECam Local Volume Exploration survey. Based on deeper follow-up imaging with DECam, we find that Aquarius III is a low-luminosity ( M V = 2.5 0.5 + 0.3 ; L V = 850 260 + 380 L ), extended ( r 1 / 2 = 41 8 + 9 pc) stellar system located in the outer halo (D= 85 ± 4 kpc). From medium-resolution Keck/DEIMOS spectroscopy, we identify 11 member stars and measure a mean heliocentric radial velocity of v sys = 13.1 0.9 + 1.0 km s 1 for the system and place an upper limit ofσv< 3.5 km s−1v< 1.6 km s−1) on its velocity dispersion at the 95% (68%) credible level. Based on calcium-triplet metallicities of the six brightest red giant members, we find that Aquarius III is very metal-poor ([Fe/H]= − 2.61 ± 0.21) with a statistically significant metallicity spread ( σ [ Fe / H ] = 0.46 0.14 + 0.26 dex). We interpret this metallicity spread as strong evidence that the system is a dwarf galaxy as opposed to a star cluster. Combining our velocity measurement with Gaia proper motions, we find that Aquarius III is currently situated near its orbital pericenter in the outer halo (rperi= 78 ± 7 kpc) and that it is plausibly on first infall onto the Milky Way. This orbital history likely precludes significant tidal disruption from the Galactic disk, notably unlike other satellites with comparably low velocity dispersion limits in the literature. Thus, if further velocity measurements confirm that its velocity dispersion is truly belowσv≲ 2 km s−1, Aquarius III may serve as a useful laboratory for probing galaxy formation physics in low-mass halos. 
    more » « less
  2. Abstract We present the analysis of a microlensing event KMT-2022-BLG-0086 of which the overall light curve is not described by a binary-lens single-source (2L1S) model, which suggests the existence of an extra lens or an extra source. We found that the event is best explained by the binary-lens binary-source (2L2S) model, but the 2L2S model is only favored over the triple-lens single-source (3L1S) model by Δχ2 ≃ 9. Although the event has noticeable anomalies around the peak of the light curve, they are not enough covered to constrain the angular Einstein radiusθE, thus we only measure the minimum angular Einstein radius θ E , min . From the Bayesian analysis, it is found that that the binary lens system is a binary star with masses of ( m 1 , m 2 ) = ( 0.4 6 0.25 + 0.35 M , 0.7 5 0.55 + 0.67 M ) at a distance of D L = 5.8 7 1.79 + 1.21 kpc, while the triple lens system is a brown dwarf or a massive giant planet in a low-mass binary-star system with masses of ( m 1 , m 2 , m 3 ) = ( 0.4 3 0.35 + 0.41 M , 0.05 6 0.047 + 0.055 M , 20.8 4 17.04 + 20.20 M J ) at a distance of D L = 4.0 6 3.28 + 1.39 kpc, indicating a disk lens system. The 2L2S model yields the relative lens-source proper motion ofμrel ≥ 4.6 mas yr−1that is consistent with the Bayesian result, whereas the 3L1S model yieldsμrel ≥ 18.9 mas yr−1, which is more than three times larger than that of a typical disk object of ∼6 mas yr−1and thus is not consistent with the Bayesian result. This suggests that the event is likely caused by the binary-lens binary-source model. 
    more » « less
  3. Abstract A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum f 0 ( , p ) with the Green’s function G ( r , p ; p ) , which describes the monoenergetic spectrum solution in which f 0 δ ( p p ) asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for G ( r , p ; p ) . In this paper, we explore for the first time, solutions for more general and realistic forms for f 0 ( , p ) . The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time τ ( r , p ) = τ 0 ( p / p 0 ) α in the shear flow region 0 <r<r2, and τ ( r , p ) = τ 0 ( p / p 0 ) α ( r / r 2 ) s , wheres> 0 in the regionr>r2is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution ψ p ( r , p ; p ) that particles observed at (r,p) originated fromr→ ∞ with momentum p . The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described. 
    more » « less
  4. Abstract We show that the affine vertex superalgebra V k ( o s p 1 | 2 n ) V^{k}(\mathfrak{osp}_{1|2n})at generic level 𝑘 embeds in the equivariant 𝒲-algebra of s p 2 n \mathfrak{sp}_{2n}times 4 n 4nfree fermions.This has two corollaries:(1) it provides a new proof that, for generic 𝑘, the coset Com ( V k ( s p 2 n ) , V k ( o s p 1 | 2 n ) ) \operatorname{Com}(V^{k}(\mathfrak{sp}_{2n}),V^{k}(\mathfrak{osp}_{1|2n}))is isomorphic to W ( s p 2 n ) \mathcal{W}^{\ell}(\mathfrak{sp}_{2n})for = ( n + 1 ) + ( k + n + 1 ) / ( 2 k + 2 n + 1 ) \ell=-(n+1)+(k+n+1)/(2k+2n+1), and(2) we obtain the decomposition of ordinary V k ( o s p 1 | 2 n ) V^{k}(\mathfrak{osp}_{1|2n})-modules into V k ( s p 2 n ) W ( s p 2 n ) V^{k}(\mathfrak{sp}_{2n})\otimes\mathcal{W}^{\ell}(\mathfrak{sp}_{2n})-modules.Next, if 𝑘 is an admissible level and ℓ is a non-degenerate admissible level for s p 2 n \mathfrak{sp}_{2n}, we show that the simple algebra L k ( o s p 1 | 2 n ) L_{k}(\mathfrak{osp}_{1|2n})is an extension of the simple subalgebra L k ( s p 2 n ) W ( s p 2 n ) L_{k}(\mathfrak{sp}_{2n})\otimes{\mathcal{W}}_{\ell}(\mathfrak{sp}_{2n}).Using the theory of vertex superalgebra extensions, we prove that the category of ordinary L k ( o s p 1 | 2 n ) L_{k}(\mathfrak{osp}_{1|2n})-modules is a semisimple, rigid vertex tensor supercategory with only finitely many inequivalent simple objects.It is equivalent to a certain subcategory of W ( s p 2 n ) \mathcal{W}_{\ell}(\mathfrak{sp}_{2n})-modules.A similar result also holds for the category of Ramond twisted modules.Due to a recent theorem of Robert McRae, we get as a corollary that categories of ordinary L k ( s p 2 n ) L_{k}(\mathfrak{sp}_{2n})-modules are rigid. 
    more » « less
  5. Abstract We introduce a distributional Jacobian determinant det D V β ( D v ) \det DV_{\beta}(Dv)in dimension two for the nonlinear complex gradient V β ( D v ) = | D v | β ( v x 1 , v x 2 ) V_{\beta}(Dv)=\lvert Dv\rvert^{\beta}(v_{x_{1}},-v_{x_{2}})for any β > 1 \beta>-1, whenever v W loc 1 , 2 v\in W^{1\smash{,}2}_{\mathrm{loc}}and β | D v | 1 + β W loc 1 , 2 \beta\lvert Dv\rvert^{1+\beta}\in W^{1\smash{,}2}_{\mathrm{loc}}.This is new when β 0 \beta\neq 0.Given any planar ∞-harmonic function 𝑢, we show that such distributional Jacobian determinant det D V β ( D u ) \det DV_{\beta}(Du)is a nonnegative Radon measure with some quantitative local lower and upper bounds.We also give the following two applications. Applying this result with β = 0 \beta=0, we develop an approach to build up a Liouville theorem, which improves that of Savin.Precisely, if 𝑢 is an ∞-harmonic function in the whole R 2 \mathbb{R}^{2}with lim inf R inf c R 1 R B ( 0 , R ) | u ( x ) c | d x < , \liminf_{R\to\infty}\inf_{c\in\mathbb{R}}\frac{1}{R}\barint_{B(0,R)}\lvert u(x)-c\rvert\,dx<\infty,then u = b + a x u=b+a\cdot xfor some b R b\in\mathbb{R}and a R 2 a\in\mathbb{R}^{2}.Denoting by u p u_{p}the 𝑝-harmonic function having the same nonconstant boundary condition as 𝑢, we show that det D V β ( D u p ) det D V β ( D u ) \det DV_{\beta}(Du_{p})\to\det DV_{\beta}(Du)as p p\to\inftyin the weak-⋆ sense in the space of Radon measure.Recall that V β ( D u p ) V_{\beta}(Du_{p})is always quasiregular mappings, but V β ( D u ) V_{\beta}(Du)is not in general. 
    more » « less