Search

The origin of magnetic fields in white dwarfs (WDs) remains mysterious. Magnetic WDs are traditionally associated with field strengths $\gtrsim1\,\mathrm{MG}$, set by the sensitivity of typical spectroscopic magnetic field measurements. Informed by recent developments in red giant magnetoasteroseismology, we revisit the use of WD pulsations as a seismic magnetometer. WD pulsations primarily probe near-surface magnetic fields, whose effect on oscillation mode frequencies is to asymmetrize rotational multiplets and, if strong enough, suppress gravity-mode propagation altogether. The sensitivity of seismology to magnetic fields increases strongly with mode period and decreases quickly with the depth of the partial ionization-driven surface convective zone. We place upper limits for magnetic fields in $24$ pulsating WDs: $20$ hydrogen-atmosphere (DAV) and three helium-atmosphere (DBV) carbon-oxygen WDs, and one extremely low-mass (helium-core) pulsator. These bounds are typically $\sim1$-$10\,\mathrm{kG}$, although they can reach down to $\sim10$-$100\,\mathrm{G}$ for DAVs and helium-core WDs in which lower-frequency modes are excited. Seismic magnetometry may enable new insights into the formation and evolution of WD magnetism., Comment: 21 pages, 9 figures, submitted to ApJ

Mergers between helium white dwarfs and main-sequence stars are likely common, producing red giant-like remnants making up roughly a few percent of all low-mass ($\lesssim2M_\odot$) red giants. Through detailed modeling, we show that these merger remnants possess distinctive photometric, asteroseismic, and surface abundance signatures through which they may be identified. During hydrogen shell burning, merger remnants reach higher luminosities and possess pulsations which depart from the usual degenerate sequence on the asteroseismic $\Delta\nu$--$\Delta\Pi$ diagram for red giant branch stars. For sufficiently massive helium white dwarfs, merger remnants undergo especially violent helium flashes which can dredge up a large amount of core material (up to $\sim0.1M_\odot$) into the envelope. Such post-dredge-up remnants are more luminous than normal red clump stars, are surface carbon-, helium-, and possibly lithium-rich, and possess a wider range of asteroseismic g-mode period spacings and mixed-mode couplings. Recent asteroseismically determined low-mass ($\lesssim0.8M_\odot$) red clump stars may be core helium-burning remnants of mergers involving lower-mass helium white dwarfs., Comment: 24 pages, 16 figures, accepted to the Open Journal of Astrophysics!

Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively. We provide a compact asymptotic formula for the g-mode period spacing, and universally find that strong magnetism decreases this period spacing substantially more than is predicted by perturbation theory. These results are validated with explicit numerical mode calculations for realistic stellar models. The approach we present is highly versatile: once the eigenvalues $\lambda$ of a certain differential operator are precomputed as a function of the magnetogravity and rotational frequencies (in units of the mode frequency), the non-perturbative impact of the Coriolis and Lorentz forces is understood under a broad domain of validity, and is readily incorporated into asteroseismic modeling., Comment: 18 pages, 8 figures; accepted to MNRAS

When main-sequence stars expand into red giants, they are expected to engulf close-in planets. Until now, the absence of planets with short orbital periods around post-expansion, core-helium-burning red giants has been interpreted as evidence that short-period planets around Sun-like stars do not survive the giant expansion phase of their host stars. Here we present the discovery that the giant planet 8 Ursae Minoris b orbits a core-helium-burning red giant. At a distance of only 0.5 au from its host star, the planet would have been engulfed by its host star, which is predicted by standard single-star evolution to have previously expanded to a radius of 0.7 au. Given the brief lifetime of helium-burning giants, the nearly circular orbit of the planet is challenging to reconcile with scenarios in which the planet survives by having a distant orbit initially. Instead, the planet may have avoided engulfment through a stellar merger that either altered the evolution of the host star or produced 8 Ursae Minoris b as a second-generation planet. This system shows that core-helium-burning red giants can harbour close planets and provides evidence for the role of non-canonical stellar evolution in the extended survival of late-stage exoplanetary systems., Comment: Published in Nature on 28 June 2023. In press

Strong magnetic fields in the cores of stars are expected to significantly modify the behavior of gravity waves: this is likely the origin of suppressed dipole modes observed in many red giants. However, a detailed understanding of how such fields alter the spectrum and spatial structure of magnetogravity waves has been elusive. For a dipole field, we analytically characterize the horizontal eigenfunctions of magnetogravity modes, assuming that the wavevector is primarily radial. For axisymmetric modes ($m=0$), the magnetogravity wave eigenfunctions become Hough functions, and they have a radial turning point for sufficiently strong magnetic fields. For non-axisymmetric modes ($m\neq0$), the interaction between the discrete $g$ mode spectrum and a continuum of Alfv\'en waves produces nearly discontinuous features in the fluid displacements at critical latitudes associated with a singularity in the fluid equations. effects. We find that magnetogravity modes cannot propagate in regions with sufficiently strong magnetic fields, instead becoming evanescent. When encountering strong magnetic fields, ingoing gravity waves are likely refracted into outgoing slow magnetic waves. These outgoing waves approach infinite radial wavenumbers, which are likely to be damped efficiently. However, it may be possible for a small fraction of the wave power to escape the stellar core as pure Alfv\'en waves or magnetogravity waves confined to a very narrow equatorial band. The artificially sharp features in the WKB-separated solutions suggest the need for global mode solutions which include small terms neglected in our analysis., Comment: 22 pages, 10 figures, 1 table, comments welcome! (submitted to MNRAS)

The globular cluster 47~Tucanae (47~Tuc) is one of the most massive star clusters in the Milky Way and is exceptionally rich in exotic stellar populations. For several decades it has been a favorite target of observers, and yet it is computationally very challenging to model because of its large number of stars ($N\gtrsim 10^6$) and high density. Here we present detailed and self-consistent 47~Tuc models computed with the \texttt{Cluster Monte Carlo} code (\texttt{CMC}). The models include all relevant dynamical interactions coupled to stellar and binary evolution, and reproduce various observations, including the surface brightness and velocity dispersion profiles, pulsar accelerations, and numbers of compact objects. We show that the present properties of 47~Tuc are best reproduced by adopting an initial stellar mass function that is both bottom-heavy and top-light relative to standard assumptions \citep[as in, e.g.,][]{Kroupa2001}, and an initial Elson profile \citep{Elson1987} that is overfilling the cluster's tidal radius. We include new prescriptions in \texttt{CMC} for the formation of binaries through giant star collisions and tidal captures, and we show that these mechanisms play a crucial role in the formation of neutron star binaries and millisecond pulsars in 47~Tuc; our best-fit model contains $\sim 50$ millisecond pulsars, $70\%$ of which are formed through giant collisions and tidal captures. Our models also suggest that 47~Tuc presently contains up to $\sim 200$ stellar-mass black holes, $\sim 5$ binary black holes, $\sim 15$ low-mass X-ray binaries, and $\sim 300$ cataclysmic variables., Comment: 23 pages, 15 figures, 4 tables, published at ApJ

Stellar mergers are important processes in stellar evolution, dynamics, and transient science. However, it is difficult to identify merger remnant stars because they cannot easily be distinguished from single stars based on their surface properties. We demonstrate that merger remnants can potentially be identified through asteroseismology of red giant stars using measurements of the gravity mode period spacing together with the asteroseismic mass. For mergers that occur after the formation of a degenerate core, remnant stars have over-massive envelopes relative to their cores, which is manifested asteroseismically by a g~mode period spacing smaller than expected for the star's mass. Remnants of mergers which occur when the primary is still on the main sequence or whose total mass is less than $\approx\! 2 \, M_\odot$ are much harder to distinguish from single stars. Using the red giant asteroseismic catalogs of Vrard et al. 2016 and Yu et al. 2018, we identify $24$ promising candidates for merger remnant stars. In some cases, merger remnants could also be detectable using only their temperature, luminosity, and asteroseismic mass, a technique that could be applied to a larger population of red giants without a reliable period spacing measurement., Comment: 13 pages, 6 figures, 1 table, submitted to MNRAS

We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using H\'enon's method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release of CMC is pinned directly to the COSMIC population synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest ($N = 10^8$) star-by-star $N$-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, a recently reported merger with component masses in the pulsational pair-instability mass gap)., Comment: Code is available at https://clustermontecarlo.github.io/ 25 pages, 8 Figures, Matches version accepted by ApJS

Numerical and observational evidence suggests that massive white dwarfs dominate the innermost regions of core-collapsed globular clusters by both number and total mass. Using NGC 6397 as a test case, we constrain the features of white dwarf populations in core-collapsed clusters, both at present day and throughout their lifetimes. The dynamics of these white dwarf subsystems have a number of astrophysical implications. We demonstrate that the collapse of globular cluster cores is ultimately halted by the dynamical burning of white dwarf binaries. We predict core-collapsed clusters in the local universe yield a white dwarf merger rate of $\mathcal{O}(10\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$, roughly $0.1-1\%$ of the observed Type Ia supernova rate. We show that prior to merger, inspiraling white dwarf binaries will be observable as gravitational wave sources at milli- and decihertz frequencies. Over $90\%$ of these mergers have a total mass greater than the Chandrasekhar limit. If the merger/collision remnants are not destroyed completely in an explosive transient, we argue the remnants may be observed in core-collapsed clusters as either young neutron stars/pulsars/magnetars (in the event of accretion-induced collapse) or as young massive white dwarfs offset from the standard white dwarf cooling sequence. Finally, we show collisions between white dwarfs and main sequence stars, which may be detectable as bright transients, occur at a rate of $\mathcal{O}(100\rm{)\,Gpc}^{-3}\,\rm{yr}^{-1}$ in the local universe. We find that these collisions lead to depletion of blue straggler stars and main sequence star binaries in the centers of core-collapsed clusters., Comment: Submitted to ApJ, 26 pages, 12 figures, 3 tables. Comments welcome

Recently, \citet{vitral2021does} detected a central concentration of dark objects in the core-collapsed globular cluster NGC 6397, which could be interpreted as a subcluster of stellar-mass black holes. However, it is well established theoretically that any significant number of black holes in the cluster would provide strong dynamical heating and is fundamentally inconsistent with this cluster's core-collapsed profile. Claims of intermediate-mass black holes in core-collapsed clusters should similarly be treated with suspicion, for reasons that have been understood theoretically for many decades. Instead, the central dark population in NGC 6397 is exactly accounted for by a compact subsystem of white dwarfs, as we demonstrate here by inspection of a previously published model that provides a good fit to this cluster. These central subclusters of heavy white dwarfs are in fact a generic feature of core-collapsed clusters, while central black hole subclusters are present in all {\em non\/}-collapsed clusters., Comment: 3 pages, 1 figure, to submit to AAS journals

As ancient, gravitationally bound stellar populations, globular clusters are abundant, vibrant laboratories characterized by high frequencies of dynamical interactions coupled to complex stellar evolution. Using surface brightness and velocity dispersion profiles from the literature, we fit $59$ Milky Way globular clusters to dynamical models from the \texttt{CMC Cluster Catalog}. Without doing any interpolation, and without any directed effort to fit any particular cluster, $26$ globular clusters are well-matched by at least one of our models. We discuss in particular the core-collapsed clusters NGC 6293, NGC 6397, NGC 6681, and NGC 6624, and the non-core-collapsed clusters NGC 288, NGC 4372, and NGC 5897. As NGC 6624 lacks well-fitting snapshots on the main \texttt{CMC Cluster Catalog}, we run six additional models in order to refine the fit. We calculate metrics for mass segregation, explore the production of compact object sources such as millisecond pulsars, cataclysmic variables, low-mass X-ray binaries, and stellar-mass black holes, finding reasonable agreement with observations. Additionally, closely mimicking observational cuts, we extract the binary fraction from our models, finding good agreement except in the dense core regions of core-collapsed clusters. Accompanying this paper are a number of \textsf{python} methods for examining the publicly accessible \texttt{CMC Cluster Catalog}, as well as any other models generated using \texttt{CMC}., Comment: 22 pages, 11 figures, 5 tables; accepted to ApJ

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters., Comment: 50 pages, 15 figures, 10 tables. Accepted for publication in ApJS

The Quintuplet star cluster is one of only three known young ($<10$ Myr) massive (M $>10^4$ M$_\odot$) clusters within $\sim100$ pc of the Galactic Center. In order to explore star cluster formation and evolution in this extreme environment, we analyze the Quintuplet's dynamical structure. Using the HST WFC3-IR instrument, we take astrometric and photometric observations of the Quintuplet covering a $120''\times120''$ field-of-view, which is $19$ times larger than those of previous proper motion studies of the Quintuplet. We generate a catalog of the Quintuplet region with multi-band, near-infrared photometry, proper motions, and cluster membership probabilities for $10,543$ stars. We present the radial density profile of $715$ candidate Quintuplet cluster members with $M\gtrsim4.7$ M$_\odot$ out to $3.2$ pc from the cluster center. A $3\sigma$ lower limit of $3$ pc is placed on the tidal radius, indicating the lack of a tidal truncation within this radius range. Only weak evidence for mass segregation is found, in contrast to the strong mass segregation found in the Arches cluster, a second and slightly younger massive cluster near the Galactic Center. It is possible that tidal stripping hampers a mass segregation signature, though we find no evidence of spatial asymmetry. Assuming that the Arches and Quintuplet formed with comparable extent, our measurement of the Quintuplet's comparatively large core radius of $0.62^{+0.10}_{-0.10}$ pc provides strong empirical evidence that young massive clusters in the Galactic Center dissolve on a several Myr timescale., Comment: 25 pages (21-page main text, 4-page appendix), 18 figures, submitted to ApJ