Researcher profile

Adam S. Jermyn

Adam S. Jermyn contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
14topics
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

25 published item(s)

preprint2022arXiv

A Transparent Window into Early-Type Stellar Variability

Subsurface convection zones are ubiquitous in early-type stars. Driven by narrow opacity peaks, these thin convective regions transport little heat but play an important role in setting the magnetic properties and surface variability of stars. Here we demonstrate that these convection zones are \emph{not} present in as wide a range of stars as previously believed. In particular, there are regions which 1D stellar evolution models report to be convectively unstable but which fall below the critical Rayleigh number for onset of convection. For sub-solar metallicity this opens up a \emph{stability window} in which there are no subsurface convection zones. For LMC metallicity this surface stability region extends roughly between $8M_\odot$ and $16M_\odot$, increasing to $8M_\odot$ -- $35M_\odot$ for SMC metallicity. Such windows are then an excellent target for probing the relative influence of subsurface convection and other sources of photometric variability in massive stars.

preprint2022arXiv

An Atlas of Convection in Main-Sequence Stars

Convection is ubiquitous in stars and occurs under many different conditions. Here we explore convection in main-sequence stars through two lenses: dimensionless parameters arising from stellar structure and parameters which emerge from the application of mixing length theory. We first define each quantity in terms familiar both to the 1D stellar evolution community and the hydrodynamics community. We then explore the variation of these quantities across different convection zones, different masses, and different stages of main-sequence evolution. We find immense diversity across stellar convection zones. Convection occurs in thin shells, deep envelopes, and nearly-spherical cores; it can be efficient of inefficient, rotationally constrained or not, transsonic or deeply subsonic. This atlas serves as a guide for future theoretical and observational investigations by indicating which regimes of convection are active in a given star, and by describing appropriate model assumptions for numerical simulations.

preprint2022arXiv

Convective Penetration in Early-Type Stars

Observations indicate that the convective cores of stars must ingest a substantial amount of material from the overlying radiative zone, but the extent of this mixing and the mechanism which causes it remain uncertain. Recently, Anders et al. (2021) developed a theory of convective penetration and calibrated it with 3D numerical hydrodynamics simulations. Here we employ that theory to predict the extent of convective boundary mixing in early-type main-sequence stars. We find that convective penetration produces enough mixing to explain core masses inferred from asteroseismology and eclipsing binary studies, and matches observed trends in mass and age. While there are remaining uncertainties in the theory, this agreement suggests that most convective boundary mixing in early-type main-sequence stars arises from convective penetration. Finally, we provide a fitting formula for the extent of core convective penetration for main-sequence stars in the mass range from $1.1-60 M_\odot$.

preprint2022arXiv

Effects of an Immortal Stellar Population in AGN Disks

Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial AGN disk. We find numbers of order $10^{2-4}$, confined to the inner $r_{\rm cap} \sim 3000 r_s \sim 0.03\rm pc$. These stars can profoundly alter the chemistry of AGN disks, enriching them in helium and depleting them in hydrogen, both by order-unity amounts. We further consider mergers between these stars and other disk objects, suggesting that star-star mergers result in rapid mass loss from the remnant to restore an equilibrium mass, while star-compact object mergers may result in exotic outcomes and even host binary black hole mergers within themselves. Finally, we examine how these stars react as the disk dissipates towards the end of its life, and find that they may return mass to the disk fast enough to extend its lifetime by a factor of several and/or may drive powerful outflows from the disk. Post-AGN, these stars rapidly lose mass and form a population of stellar mass black holes around $10M_{\odot}$. Due to the complex and uncertain interactions between embedded stars and the disk, their plausible ubiquity, and their order unity impact on disk structure and evolution, they must be included in realistic disk models.

preprint2022arXiv

Implications of a Temperature Dependent IMF I: Photometric Template Fitting

A universal stellar initial mass function (IMF) should not be expected from theoretical models of star formation, but little conclusive observational evidence for a variable IMF has been uncovered. In this paper, a parameterization of the IMF is introduced into photometric template fitting of the COSMOS2015 catalog. The resulting best-fit templates suggest systematic variations in the IMF, with most galaxies exhibiting top-heavier stellar populations than in the Milky Way. At fixed redshift, only a small range of IMFs are found, with the typical IMF becoming progressively top-heavier with increasing redshift. Additionally, subpopulations of ULIRGs, quiescent- and star-forming galaxies are compared with predictions of stellar population feedback and show clear qualitative similarities to the evolution of dust temperatures.

preprint2022arXiv

Implications of a Temperature Dependent IMF II: An Updated View of the Star-Forming Main Sequence

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, allows galaxies to be fit with a range of IMFs. Three surprising new features appear: (1) most star-forming galaxies are best fit with a bottom-lighter IMF than the Milky Way; (2) most star-forming galaxies at fixed redshift are fit with a very similar IMF; and (3) the most massive star-forming galaxies at fixed redshift instead exhibit a less bottom-light IMF, similar to that measured in quiescent galaxies. Additionally, since stellar masses and star formation rates both depend on the IMF, these results slightly modify the resulting relationship, while yielding similar qualitative characteristics to previous studies.

preprint2022arXiv

Implications of a Temperature Dependent IMF III: Mass Growth and Quiescence

The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, suggest most galaxies obey an IMF top-heavier than the Galactic IMF. The implications of these revised fits on mass functions, quiescence and turnoff are discussed. At all redshifts the highest mass galaxies become quiescent first with the turnoff mass decreasing towards the present. The synchronous turnoff mass across galaxies suggests quiescence is driven by universal mechanisms rather than by stochastic or environmental processes.

preprint2022arXiv

Inferring the Gas-to-Dust Ratio in the Main Planet-Forming Region of Disks

Measuring the amount of gas and dust in protoplanetary disks is a key challenge in planet formation studies. Here we provide a new set of dust depletion factors and relative mass surface densities of gas and dust for the innermost regions of a sample of protoplanetary disks. We do this by combining stellar theory with observed refractory element abundances in both disk hosts and open cluster stars. Our results are independent of, and complementary to, those obtained from spatially resolved disk observations.

preprint2022arXiv

Moosinesq Convection in the Cores of Moosive Stars

Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure of the convection zones of these stars. These flows are bounded in a manner informed by the majestic terrestrial Alces alces (moose) and could have important consequences for the evolution of these stars. We find that Moosinesq convection results in very interesting flow morphologies and rapid heat transfer, and posit this as a mechanism of biomechanical thermoregulation.

preprint2022arXiv

Probing carrier interactions using electron hydrodynamics

Electron hydrodynamics arises when momentum-relaxing scattering processes are slow compared to momentum-conserving ones. While the microscopic details necessary to satisfy this condition are material-specific, experimentally accessible current densities share remarkable similarities. We study the dependence of electron hydrodynamic flows on the rates of momentum-relaxing and momentum-conserving scattering processes in a microscopics-agnostic way. We develop a framework for generating random collision operators which respect crystal symmetries and conservation laws and which have a tunable ratio between the momentum-conserving and momentum-relaxing lifetimes. Using various random instances of these collision operators, we calculate macroscopic electron viscosity tensors and solve the Boltzmann transport equation (BTE) in a channel geometry over a grid of momentum-conserving and momentum-relaxing lifetimes, and for different crystal symmetry groups. We find that different random collision operators using the same lifetimes produce very similar current density profiles, meaning that the current density is primarily a probe of the overall rates of momentum conservation and relaxation. By contrast, the viscosity tensor varies substantially at fixed lifetimes, meaning that properties like channel resistance provide detailed probes of the underlying scattering processes. This suggests that, while details of the scattering process are imprinted in the electronic viscosity tensor, for many applications theoretical calculations of hydrodynamic electron flows can use experimentally-available lifetimes within a spatially-resolved BTE framework rather than requiring the costly computation of ab initio collision operators.

preprint2022arXiv

Schwarzschild and Ledoux are equivalent on evolutionary timescales

Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild-unstable but Ledoux-stable due to a composition gradient. Over many convective overturn timescales the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists.

preprint2022arXiv

Stellar convective penetration: parameterized theory and dynamical simulations

Most stars host convection zones in which heat is transported directly by fluid motion, but the behavior of convective boundaries is not well understood. Here we present 3D numerical simulations which exhibit penetration zones: regions where the entire luminosity \emph{could} be carried by radiation, but where the temperature gradient is approximately adiabatic and convection is present. To parameterize this effect, we define the "penetration parameter" $\mathcal{P}$ which compares how far the radiative gradient deviates from the adiabatic gradient on either side of the Schwarzschild convective boundary. Following Roxburgh (1989) and Zahn (1991), we construct an energy-based theoretical model in which $\mathcal{P}$ controls the extent of penetration. We test this theory using 3D numerical simulations which employ a simplified Boussinesq model of stellar convection. The convection is driven by internal heating and we use a height-dependent radiative conductivity; this allows us to separately specify $\mathcal{P}$ and the stiffness $\mathcal{S}$ of the radiative-convective boundary. We find significant convective penetration in all simulations. Our simple theory describes the simulations well. Penetration zones can take thousands of overturn times to develop, so long simulations or accelerated evolutionary techniques are required. In stars, we expect $\mathcal{P} \approx 1$ and in this regime our results suggest that convection zones may extend beyond the Schwarzschild boundary by up to $\sim$20-30% of a mixing length. We present a MESA stellar model of the Sun which employs our parameterization of convective penetration as a proof of concept. We discuss prospects for extending these results to more realistic stellar contexts.

preprint2022arXiv

Wave heating during the helium flash and lithium-enhanced clump stars

Red Clump stars have been found to be enhanced in lithium relative to stars at the tip of the Red Giant Branch (TRGB), which is unexpected in current stellar models. At the TRGB, stars undergo the helium flash, during which helium burning briefly generates roughly $10^9 \, L_\odot$ of power and drives vigorous convection within the star's core. The helium-burning shell excites large fluxes of internal gravity waves. Here we investigate whether or not these waves can deposit enough heat to destabilize the hydrogen-burning shell, generate a convection zone there, and thereby drive the Cameron-Fowler process to enhance surface $^{7}\mathrm{Li}$. We study this with detailed stellar evolution models, and find that while the waves deposit $\sim 10^6 \, L_\odot$ near the hydrogen-burning shell, and while this generally does produce a convection zone, the resulting convection does not reach high enough to merge with the envelope, and so cannot explain enhancements to surface $^{7}\mathrm{Li}$.

preprint2021arXiv

On the Origin of Stochastic, Low-Frequency Photometric Variability in Massive Stars

High-precision photometric observations have revealed ubiquitous stochastic low-frequency photometric variability in early type stars. It has been suggested that this variability arises due to either subsurface convection or internal gravity waves launched by the convective core. Here we show that relevant properties of convection in subsurface convective layers correlate very well with the timescale and amplitude of stochastic low-frequency photometric variability, as well as with the amplitude of macroturbulence. We suggest that low-frequency, stochastic photometric variability and surface turbulence in massive stars are caused by the the presence of subsurface convection. We show that an explanation for the observed surface photometric variability and macroturbulence relying on convective core driven internal gravity waves encounters a number of difficulties and seems unlikely to be able to explain the observed trends.

preprint2021arXiv

Skye: A Differentiable Equation of State

Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software instrument by computing white dwarf cooling curves.

preprint2021arXiv

The excess of cool supergiants from contemporary stellar evolution models defies the metallicity-independent Humphreys-Davidson limit

The Humphreys-Davidson (HD) limit empirically defines a region of high luminosities (log L > 5.5) and low effective temperatures (T < 20kK) on the Hertzsprung-Russell Diagram in which hardly any supergiant stars are observed. Attempts to explain this limit through instabilities arising in near- or super-Eddington winds have been largely unsuccessful. Using modern stellar evolution we aim to re-examine the HD limit, investigating the impact of enhanced mixing on massive stars. We construct grids of stellar evolution models appropriate for the Small and Large Magellanic Clouds (SMC, LMC), as well as for the Galaxy, spanning various initial rotation rates and convective overshooting parameters. Significantly enhanced mixing apparently steers stellar evolution tracks away from the region of the HD limit. To quantify the excess of over-luminous stars in stellar evolution simulations we generate synthetic populations of massive stars, and make detailed comparisons with catalogues of cool (T < 12.5kK) and luminous (log L > 4.7) stars in the SMC and LMC. We find that adjustments to the mixing parameters can lead to agreement between the observed and simulated red supergiant populations, but for hotter supergiants the simulations always over-predict the number of very luminous (log L > 5.4) stars compared to observations. The excess of luminous supergiants decreases for enhanced mixing, possibly hinting at an important role mixing has in explaining the HD limit. Still, the HD limit remains unexplained for hotter supergiants.

preprint2020arXiv

Algorithms for Tensor Network Contraction Ordering

Contracting tensor networks is often computationally demanding. Well-designed contraction sequences can dramatically reduce the contraction cost. We explore the performance of simulated annealing and genetic algorithms, two common discrete optimization techniques, to this ordering problem. We benchmark their performance as well as that of the commonly-used greedy search on physically relevant tensor networks. Where computationally feasible, we also compare them with the optimal contraction sequence obtained by an exhaustive search. We find that the algorithms we consider consistently outperform a greedy search given equal computational resources, with an advantage that scales with tensor network size. We compare the obtained contraction sequences and identify signs of highly non-local optimization, with the more sophisticated algorithms sacrificing run-time early in the contraction for better overall performance.

preprint2020arXiv

Convective Differential Rotation in Stars and Planets I: Theory

We derive the scaling of differential rotation in both slowly- and rapidly-rotating convection zones using order of magnitude methods. Our calculations apply across stars and fluid planets and all rotation rates, as well as to both magnetized and purely hydrodynamic systems. We find shear $|R\nablaΩ|$ of order the angular frequency $Ω$ for slowly-rotating systems with $Ω\ll |N|$, where $N$ is the \brvs\ frequency, and find that it declines as a power-law in $Ω$ for rapidly-rotating systems with $Ω\gg |N|$. We further calculate the meridional circulation rate and baroclinicity and examine the magnetic field strength in the rapidly rotating limit. Our results are in general agreement with simulations and observations and we perform a detailed comparison with those in a companion paper.

preprint2020arXiv

Convective Differential Rotation in Stars and Planets II: Observational and Numerical Tests

Differential rotation is central to a great many mysteries in stars and planets. In Part I we predicted the order of magnitude and scaling of the differential rotation in both hydrodynamic and magnetohydrodynamic convection zones. Our results apply to both slowly- and rapidly-rotating systems, and provide a general picture of differential rotation in stars and fluid planets. We further calculated the scalings of the meridional circulation, entropy gradient and baroclinicity. In this companion paper we compare these predictions with a variety of observations and numerical simulations. With a few exceptions we find that these are consistent in both the slowly-rotating and rapidly-rotating limits. Our results help to localize core-envelope shear in red~giant stars, suggest a rotation-dependent frequency shift in the internal gravity waves of massive stars and potentially explain observed deviations from von Zeipel&#39;s gravity darkening in late-type stars.

preprint2020arXiv

Generalized Electron Hydrodynamics, Vorticity Coupling, and Hall Viscosity in Crystals

Theoretical and experimental studies have revealed that electrons in condensed matter can behave hydrodynamically, exhibiting fluid phenomena such as Stokes flow and vortices. Unlike classical fluids, preferred directions inside crystals lift isotropic restrictions, necessitating a generalized treatment of electron hydrodynamics. We explore electron fluid behaviors arising from the most general viscosity tensors in two and three dimensions, constrained only by thermodynamics and crystal symmetries. Hexagonal 2D materials such as graphene support flows indistinguishable from those of an isotropic fluid. By contrast 3D materials including Weyl semimetals, exhibit significant deviations from isotropy. Breaking time-reversal symmetry, for example in magnetic topological materials, introduces a non-dissipative Hall component to the viscosity tensor. While this vanishes by isotropy in 3D, anisotropic materials can exhibit nonzero Hall viscosity components. We show that in 3D anisotropic materials the electronic fluid stress can couple to the vorticity without breaking time-reversal symmetry. Our work demonstrates the anomalous landscape for electron hydrodynamics in systems beyond graphene, and presents experimental geometries to quantify the effects of electronic viscosity.

preprint2020arXiv

Multiphase Gas and the Fractal Nature of Radiative Turbulent Mixing Layers

A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $ρ_{\rm cold}/ρ_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $ρ_{\rm cold}/ρ_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling.

preprint2020arXiv

The Origin of the Bimodal Distribution of Magnetic Fields in Early-type Stars

In early-type stars a fossil magnetic field may be generated during the star formation process or be the result of a stellar merger event. Surface magnetic fields are thought to be erased by (sub)surface convection layers, which typically leave behind weak disordered fields. However, if the fossil field is strong enough it can prevent the onset of (sub)surface convection and so be preserved onto the main sequence. We calculate the critical field strength at which this occurs, and find that it corresponds well with the lower limit amplitude of observed fields in strongly magnetised Ap/Bp stars ($\approx$ 300 G). The critical field strength is predicted to increase slightly during the main sequence evolution, which could also explain the observed decline in the fraction of magnetic stars. This supports the conclusion that the bimodal distribution of observed magnetic fields in early-type stars reflects two different field origin stories: strongly magnetic fields are fossils fields inherited from star formation or a merger event, and weak fields are the product of on-going dynamo action.

preprint2019arXiv

Automatic Contraction of Unstructured Tensor Networks

The evaluation of partition functions is a central problem in statistical physics. For lattice systems and other discrete models the partition function may be expressed as the contraction of a tensor network. Unfortunately computing such contractions is difficult, and many methods to make this tractable require periodic or otherwise structured networks. Here I present a new algorithm for contracting unstructured tensor networks. This method makes no assumptions about the structure of the network and performs well in both structured and unstructured cases so long as the correlation structure is local.

preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise&#39;&#39; profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

preprint2019arXiv

Thermal Regulation and the Star-Forming Main Sequence

We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming &#34;main sequence&#34;. To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.