Researcher profile

Morgan MacLeod

Morgan MacLeod contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2026arXiv

Betelgeuse: Detection of the Expanding Wake of the Companion Star

Recent analyses conclude that Betelgeuse, a red supergiant star (HD 39801), likely has a companion object with a period of about 2000 days orbiting at only 2.3 stellar radii, deep in the chromosphere of the supergiant. A probable detection of such a companion, named Siwarha, has just occurred from speckle imaging. This study finds that Betelgeuse spectra in the optical region and ultraviolet exhibit signatures of variable circumstellar absorption and chromospheric outflows. These variations are consistent with the ~ 2000-day period of the companion object. Circumstellar absorption evident in optical Mn I lines, and mass outflow marked by ultraviolet Fe II, Si I, and Mg I lines increase after the transit of the companion across the disk of Betelgeuse. Following the eclipse of the companion, the absorption and outflow slowly decrease in advance of the next transit. The occurrence and variation of this plasma appear consistent with the presence of a trailing and expanding wake caused by a companion star orbiting within the atmosphere of Betelgeuse.

preprint2025arXiv

On the coincidence between the close passage of HD7977 and the Pliocene-Pleistocene transition

The Oort Cloud's dynamical evolution is significantly influenced by both the galactic tide and stellar flybys. This study investigates the particular case of HD7977's close encounter 2.47 Myr ago, which likely repopulated the Inner Oort Cloud and potentially triggered a significant comet shower on Earth. Our results demonstrate that the shower's intensity strongly depends on HD7977's impact parameter ($b$), with possible flyby distances ranging from 2,300 AU to $\sim$ 13,000 AU. For the closest approach ($b \sim 2,300$ AU), the terrestrial impact probability of 1 km comets increases by an order of magnitude compared to the steady state, slightly exceeding the asteroid impact probability at this size scale. We propose an analytical method to compute the probability of comet showers impacting Earth, which saves considerable computation time compared to N-body simulations. We identify a threshold diameter $D_0 = 2.25$ km for which yields $P = 1$ in our model, with $D_0$ following a logarithmic dependence on $b$. These findings suggest that HD7977's flyby may have caused an enhanced comet flux during the Pliocene-Pleistocene transition, which could plausibly be related to the environmental changes at this era.

preprint2023arXiv

Formation of black holes in the pair-instability mass gap: Hydrodynamical simulations of a head-on massive star collision

The detection of the binary black hole merger GW190521, with primary black hole mass $85^{+21}_{-14}$ ${\rm M}_{\odot}$, proved the existence of black holes in the theoretically predicted pair-instability gap ($\sim60-120 \, {\rm M}_{\odot}$) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon-oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code {\sc StarSmasher}, finding that a head-on collision can remove up to 12\% of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product.

preprint2022arXiv

Eccentricity Evolution in Gaseous Dynamical Friction

We analyse how drag forces modify the orbits of objects moving through extended gaseous distributions. We consider how hydrodynamic (surface area) drag forces and dynamical friction (gravitational) drag forces drive the evolution of orbital eccentricity. While hydrodynamic drag forces cause eccentric orbits to become more circular, dynamical friction drag can cause orbits to become more eccentric. We develop a semi-analytic model that accurately predicts these changes by comparing the total work and torque applied to the orbit at periapse and apoapse. We use a toy model of a radial power-law density profile, $ρ\propto r^{-γ}$, to determine that there is a critical $γ= 3$ power index which separates the eccentricity evolution in dynamical friction: orbits become more eccentric for $γ< 3$ and circularize for $γ> 3$. We apply these findings to the infall of a Jupiter-like planet into the envelope of its host star. The hydrostatic envelopes of stars are defined by steep density gradients near the limb and shallower gradients in the interior. Under the influence of gaseous dynamical friction, an infalling object&#39;s orbit will first decrease in eccentricity, then increase. The critical separation that delineates these regimes is predicted by the local density slope and is linearly dependent on polytropic index. More broadly, our findings indicate that binary systems may routinely emerge from common envelope phases with non-zero eccentricities that were excited by the dynamical friction forces that drove their orbital tightening.

preprint2022arXiv

Keck/NIRSPEC studies of He I in the atmospheres of two inflated hot gas giants orbiting K dwarfs: WASP-52b and WASP-177b

We present the detection of neutral helium at 10833A in the atmosphere of WASP-52b and tentative evidence of helium in the atmosphere of the grazing WASP-177b, using high-resolution observations acquired with the NIRSPEC instrument on the Keck II telescope. We detect excess absorption by helium in WASP-52b&#39;s atmosphere of $3.44 \pm 0.31$% ($11σ$), or equivalently $66 \pm 5$ atmospheric scale heights. This absorption is centered on the planet&#39;s rest frame ($Δv = 0.00 \pm 1.19$km s$^{-1}$). We model the planet&#39;s escape using a 1D Parker wind model and calculate its mass-loss rate to be $\sim 1.4 \times 10^{11}$g s$^{-1}$, or equivalently 0.5% of its mass per Gyr. For WASP-177b, we see evidence for red-shifted ($Δv = 6.02 +/- 1.88$km s$^{-1}$) helium-like absorption of $1.28 \pm 0.29$% (equal to $23 \pm 5$ atmospheric scale heights). However, due to residual systematics in the transmission spectrum of similar amplitude, we do not interpret this as significant evidence for He absorption in the planet&#39;s atmosphere. Using a 1D Parker wind model, we set a $3σ$ upper limit on WASP-177b&#39;s escape rate of $7.9 \times 10^{10}$ g s$^{-1}$. Our results, taken together with recent literature detections, suggest the tentative relation between XUV irradiation and He I absorption amplitude may be shallower than previously suggested. Our results highlight how metastable helium can advance our understanding of atmospheric loss and its role in shaping the exoplanet population.

preprint2022arXiv

Stellar Wind Confinement of Evaporating Exoplanet Atmospheres and its Signatures in 1083 nm Observations

Atmospheric escape from close-in exoplanets is thought to be crucial in shaping observed planetary populations. Recently, significant progress has been made in observing this process in action through excess absorption in transit spectra and narrowband light curves. We model the escape of initially-homogeneous planetary winds interacting with a stellar wind. The ram pressure balance of the two winds governs this interaction. When the impingement of the stellar wind on the planetary outflow is mild or moderate, the planetary outflow expands nearly spherically through its sonic surface before forming a shocked boundary layer. When the confinement is strong, the planetary outflow is redirected into a cometary tail before it expands to its sonic radius. The resultant transmission spectra at the He 1083 nm line are accurately represented by a 1D spherical wind solution in cases of mild to moderate stellar wind interaction. In cases of strong stellar wind interaction, the degree of absorption is enhanced and the cometary tail leads to an extended egress from transit. The crucial features of the wind--wind interaction are, therefore, encapsulated in the light curve of He 1083 nm equivalent width as a function of time. The possibility of extended He 1083 nm absorption well beyond the optical transit carries important implications for planning &#34;out-of-transit&#34; observations that serve as a baseline for in-transit data.

preprint2020arXiv

Common Envelope Wind Tunnel: Range of Applicability and Self-Similarity in Realistic Stellar Envelopes

Common envelope evolution, the key orbital tightening phase of the traditional formation channel for close binaries, is a multistage process that presents many challenges to the establishment of a fully descriptive, predictive theoretical framework. In an approach complementary to global 3D hydrodynamical modeling, we explore the range of applicability for a simplified drag formalism that incorporates the results of local hydrodynamic &#34;wind tunnel&#34; simulations into a semi-analytical framework in the treatment of the common envelope dynamical inspiral phase using a library of realistic giant branch stellar models across the low, intermediate, and high mass regimes. In terms of a small number of key dimensionless parameters, we characterize a wide range of common envelope events, revealing the broad range of applicability of the drag formalism as well its self-similar nature across mass regimes and ages. Limitations arising from global binary properties and local structural quantities are discussed together with the opportunity for a general prescriptive application for this formalism.

preprint2020arXiv

Pre-Common-Envelope Mass Loss from Coalescing Binary Systems

Binary systems undergoing unstable Roche Lobe overflow spill gas into their circumbinary environment as their orbits decay toward coalescence. In this paper, we use a suite of hydrodynamic models of coalescing binaries involving an extended donor and a more compact accretor. We focus on the period of unstable Roche Lobe overflow that ends as the accretor plunges within the envelope of the donor at the onset of a common envelope phase. During this stage, mass is removed from the donor and flung into the circumbinary environment. Across a wide range of binary mass ratios, we find that the mass expelled as the separation decreases from the Roche limit to the donor&#39;s original radius is of the order of 25% of the accretor&#39;s mass. We study the kinematics of this ejecta and its dependencies on binary properties and find that it assembles into a toroidal circumbinary distribution. These circumbinary tori have approximately constant specific angular momentum due to momentum transport by spiral shocks launched from the orbiting binary. We show that an analytic model with these torus properties captures many of the main features of the azimuthally-averaged profiles of our hydrodynamic simulations. Our results, in particular the simple relationship between accretor mass and expelled mass and its spatial distribution, may be useful in interpreting stellar coalescence transients like luminous red novae, and in initializing hydrodynamic simulations of the subsequent common envelope phase.

preprint2020arXiv

Runaway Coalescence of Pre-Common-Envelope Stellar Binaries

We study the process of runaway, unstable Roche lobe overflow in coalescing binary systems and its dependence on the properties of the binary involved. We create three-dimensional hydrodynamic models of binary coalescences, and follow them through a phase of increasing Roche lobe overflow until the accretor is engulfed by the donor at the onset of a common envelope phase. In these models, we vary binary properties of mass ratio, donor structure and spin, and equation of state through the gas adiabatic index. We compare the numerical results to semi-analytic models of binary orbit evolution based on mass and angular momentum exchange between two point masses. Using our hydrodynamic simulations, we measure the key parameters: the donor mass loss rate and the angular momentum exchanged per unit mass loss from the donor. Using these calibrations, the semi-analytic model closely reproduces the escalating mass loss and runaway orbital decay observed in the hydrodynamic models. The semi-analytic model accurately reproduces the major differences in orbit evolution that arise with varying mass ratio and donor structure. We encapsulate the semi-analytic model in a publicly-released python package RLOF. We apply this model to the observed period decay and subsequent merger of the binary V1309 Sco, and find that it can simultaneously reproduce the observed orbital decay and time of outburst. We further demonstrate that there is a relationship between period derivative and second derivative that can be a useful metric for evaluating candidate merging binaries.

preprint2020arXiv

The Effect of Star--Disk Interactions on Highly Eccentric Stellar Orbits in Active Galactic Nuclei: A Disk Loss Cone and Implications for Stellar Tidal Disruption Events

Galactic center black holes appear to be nearly universally surrounded by dense stellar clusters. When these black holes go through an active accretion phase, the multiple components of the accretion disk, stellar cluster, and black hole system all coexist. We analyze the effect of drag forces on highly eccentric stellar orbits incurred as stars puncture through the disk plane. Disk crossings dissipate orbital energy, drawing eccentric stars into more circular orbits. For high surface density disks, such as those found around black holes accreting near the Eddington mass accretion limit, the magnitude of this energy dissipation can be larger than the mean scatterings that stars receive by two body relaxation. One implication of this is the presence of a disk &#34;loss cone&#34; for highly eccentric stellar orbits where the dissipation from disk interaction outweighs scatter via two body relaxation. The disk loss cone is larger than the tidal disruption loss cone for near-Eddington black hole accretion rates. Stellar orbits within the disk loss cone are lost from the overall cluster as stellar orbits are circularized and stars are potentially ablated by their high-velocity impacts with the disk. We find, however, that the presence of the disk loss cone has a minimal effect on the overall rate of stellar tidal disruptions. Stars are still efficiently fed to the black hole from more-distant stellar orbits that receive large-enough per orbit scatter to jump over the disk loss cone and end up tidally disrupted.

preprint2020arXiv

Tidal Disruptions of White Dwarfs: Theoretical Models and Observational Prospects

White dwarf stars that enter the tidal radius of black holes with masses $\lesssim$ $10^5$ M$_{\odot}$ are doomed to be ripped apart by tidal forces. Black holes in this mass range between stellar black holes and supermassive black holes have not been conclusively identified so the detection of a tidal disruption of a white dwarf would provide clear evidence for the existence of intermediate-mass black holes. In this review, we present a theoretical and observational overview of the transient events that result from the tidal disruptions of white dwarfs by intermediate-mass black holes. This includes discussion of the latest simulations and predicted properties, the results of observational searches, as well as a summary of the potential for gravitational wave emission to be detected with upcoming missions.

preprint2019arXiv

Common Envelope Wind Tunnel: The Effects of Binary Mass Ratio and Implications for the Accretion-Driven Growth of LIGO Binary Black Holes

We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a &#34;wind tunnel&#34; formalism. Our simulations model the common envelope dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object&#39;s motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object&#39;s masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially non-rotating black holes change by of order $1\%$ and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.

preprint2019arXiv

Explosions Driven by the Coalescence of a Compact Object with the Core of a Massive-Star Companion Inside a Common Envelope: Circumstellar Properties, Light Curves, and Population Statistics

We model explosions driven by the coalescence of a black hole or neutron star with the core of its massive-star companion. Upon entering a common envelope phase, a compact object may spiral all the way to the core. The concurrent release of energy is likely to be deposited into the surrounding common envelope, powering a merger-driven explosion. We use hydrodynamic models of binary coalescence to model the common envelope density distribution at the time of coalescence. We find toroidal profiles of material, concentrated in the binary&#39;s equatorial plane and extending to many times the massive star&#39;s original radius. We use the spherically-averaged properties of this circumstellar material (CSM) to estimate the emergent light curves that result from the interaction between the blast wave and the CSM. We find that typical merger-driven explosions are brightened by up to three magnitudes by CSM interaction. From population synthesis models we discover that the brightest merger-driven explosions, $M_V \sim -18$ to $-19$, are those involving black holes because they have the most massive and extended CSM. Black hole coalescence events are also common; they represent about 50% of all merger-driven explosions and approximately 0.3% of the core-collapse rate. Merger-driven explosions offer a window into the highly-uncertain physics of common envelope interactions in binary systems by probing the properties of systems that merge rather than eject their envelopes.

preprint2017arXiv

Low-mass White Dwarfs with Hydrogen Envelopes as a Missing Link in the Tidal Disruption Menu

We construct a menu of objects that can give rise to bright flares when disrupted by massive black holes (BHs), ranging from planets to evolved stars. Through their tidal disruption, main sequence and evolved stars can effectively probe the existence of otherwise quiescent supermassive BHs and white dwarfs can probe intermediate mass BHs. Many low-mass white dwarfs possess extended hydrogen envelopes, which allow for the production of prompt flares in disruptive encounters with moderately massive BHs of $10^5$ to $10^7~M_\odot$--masses that may constitute the majority of massive BHs by number. These objects are a missing link in two ways: (1) for probing moderately massive BHs and (2) for understanding the hydrodynamics of the disruption of objects with tenuous envelopes. A flare arising from the tidal disruption of a $0.17~M_\odot$ white dwarf by a $10^5~M_\odot$ BH reaches a maximum between 0.6 and 11 days, with a peak fallback rate that is usually super-Eddington and results in a flare that is likely brighter than a typical tidal disruption event. Encounters stripping only the envelope can provide hydrogen-only fallback, while encounters disrupting the core evolve from H- to He-rich fallback. While most tidal disruption candidates observed thus far are consistent with the disruptions of main sequence stars, the rapid timescales of nuclear transients such as Dougie and PTF10iya are naturally explained by the disruption of low-mass white dwarfs. As the number of observed flares continues to increase, the menu presented here will be essential for characterizing nuclear BHs and their environments through tidal disruptions.