BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matteo Cantiello

Matteo Cantiello appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.IMphysics.flu-dynastro-ph.GAastro-ph.EPhep-ph

Catalog footprint

What is connected

33works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

33 published item(s)

preprint2022arXiv

A Measurement of Stellar Surface Gravity Hidden in Radial Velocity Differences of Co-moving Stars

The gravitational redshift induced by stellar surface gravity is notoriously difficult to measure for non-degenerate stars, since its amplitude is small in comparison with the typical Doppler shift induced by stellar radial velocity. In this study, we make use of the large observational data set of the Gaia mission to achieve a significant reduction of noise caused by these random stellar motions. By measuring the differences in velocities between the components of pairs of co-moving stars and wide binaries, we are able to statistically measure the combined effects of gravitational redshift and convective blueshifting of spectral lines, and nullify the effect of the peculiar motions of the stars. For the subset of stars considered in this study, we find a positive correlation between the observed differences in Gaia radial velocities and the differences in surface gravity and convective blueshift inferred from effective temperature and luminosity measurements. The results rule out a null signal at the $5σ$ level for our full data-set. Additionally, we study the sub-dominant effects of binary motion, and possible systematic errors in radial velocity measurements within Gaia. Results from the technique presented in this study are expected to improve significantly with data from the next Gaia data release. Such improvements could be used to constrain the mass-luminosity relation and stellar models which predict the magnitude of convective blueshift.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electromagnetic Signatures of Relativistic Explosions in AGN Disks

The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS-NS and NS-BH mergers are expected. However, their properties in the unique environments of AGN disks have never been studied. Here we show that GRBs in AGNs can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, the external shock develops before the internal shocks, leading to prompt emission powered by a relativistic reverse shock. The transient's time evolution is also compressed, yielding afterglow emission that is much brighter and peaks much earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball is trapped inside the disc photosphere. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of $\sim 10^6 M_\odot$ to years for SMBHs of $\sim 10^8 M_\odot$. Finally, a large fraction of events, especially in AGNs with SMBHs $\lesssim 10^7 M_\odot$, display time-variable absorption in the X-ray band.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'' 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Asteroseismic signatures of evolving internal stellar magnetic fields

Recent asteroseismic analyses have revealed the presence of strong (B $\gtrsim 10^5$ G) magnetic fields in the cores of many red giant stars. Here, we examine the implications of these results for the evolution of stellar magnetic fields, and we make predictions for future observations. Those stars with suppressed dipole modes indicative of strong core fields should exhibit moderate but detectable quadrupole mode suppression. The long magnetic diffusion times within stellar cores ensure that dynamo-generated fields are confined to mass coordinates within the main sequence convective core, and the observed sharp increase in dipole mode suppression rates above $1.5 \, M_\odot$ may be explained by the larger convective core masses and faster rotation of these more massive stars. In clump stars, core fields of $\sim10^5 \, {\rm G}$ can suppress dipole modes, whose visibility should be equal to or less than the visibility of suppressed modes in ascending red giants. High dipole mode suppression rates in low-mass ($M \lesssim 2 \, M_\odot$) clump stars would indicate that magnetic fields generated during the main sequence can withstand subsequent convective phases and survive into the compact remnant phase. Finally, we discuss implications for observed magnetic fields in white dwarfs and neutron stars, as well as the effects of magnetic fields in various types of pulsating stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asteroseismology Can Reveal Strong Internal Magnetic Fields in Red Giant Stars

Internal stellar magnetic fields are inaccessible to direct observations and little is known about their amplitude, geometry and evolution. We demonstrate that strong magnetic fields in the cores of red giant stars can be identified with asteroseismology. The fields can manifest themselves via depressed dipole stellar oscillation modes, which arises from a magnetic greenhouse effect that scatters and traps oscillation mode energy within the core of the star. The Kepler satellite has observed a few dozen red giants with depressed dipole modes which we interpret as stars with strongly magnetized cores. We find field strengths larger than $\sim\! 10^5 \,{\rm G}$ may produce the observed depression, and in one case we infer a minimum core field strength of $\approx \! \! 10^7 \,{\rm G}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Limits On The Spin Up Of Stellar-Mass Black Holes Through A Spiral Stationary Accretion Shock Instability

The spin of a number of black holes (BHs) in binary systems has been measured. In the case of BHs found in low-mass X-ray binaries (LMXBs) the observed values are in agreement with some theoretical predictions based on binary stellar evolution. However, using the same evolutionary models, the calculated spins of BHs in high-mass X-ray binaries (HMXBs) fall short compared to the observations. A possible solution to this conundrum is the accretion of high-specific-angular-momentum material after the formation of the BH, although this requires accretion above the Eddington limit. Another suggestion is that the observed high values of the BHs spin could be the result of an asymmetry during Core Collapse (CC). The only available energy to spin up the compact object during CC is its binding energy. A way to convert it to rotational kinetic energy is by using a Standing Accretion Shock Instability (SASI), which can develop during CC and push angular momentum into the central compact object through a spiral mode ($m = 1$). Here we study the CC-SASI scenario and discuss, in the case of LMXBs and HMXBs, the limits for the spin of a stellar-mass BHs. Our results predict a strong dichotomy in the maximum spin of low-mass compact objects and massive BHs found in HMXBs. The maximum spin value ($|a_\star|$) for a compact object near the mass boundary between BHs and NSs is found to be somewhere between 0.27 and 0.38, depending on whether secular or dynamical instabilities limit the efficiency of the spin up process. For more massive BHs, such as those found in HMXBs, the natal spin is substantially smaller and for $M_{\rm BH}\!>\!8~M_\odot$ spin is limited to values $|a_\star|\!\lesssim\!0.05$. Therefore we conclude that the observed high spins of BHs in HMXBs cannot be the result of a CC-SASI spin up.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $τ_0$, to $τ_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $τ_0 \gg τ_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $τ_0 \ll τ_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The spin rate of pre-collapse stellar cores: wave-driven angular momentum transport in massive stars

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 30 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 5 \times 10^3 \, {\rm s}$ in typical core-collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 500 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular Momentum Transport via Internal Gravity Waves in Evolving Stars

Recent asteroseismic advances have allowed for direct measurements of the internal rotation rates of many sub-giant and red giant stars. Unlike the nearly rigidly rotating Sun, these evolved stars contain radiative cores that spin faster than their overlying convective envelopes, but slower than they would in the absence of internal angular momentum transport. We investigate the role of internal gravity waves in angular momentum transport in evolving low mass stars. In agreement with previous results, we find that convectively excited gravity waves can prevent the development of strong differential rotation in the radiative cores of Sun-like stars. As stars evolve into sub-giants, however, low frequency gravity waves become strongly attenuated and cannot propagate below the hydrogen burning shell, allowing the spin of the core to decouple from the convective envelope. This decoupling occurs at the base of the sub-giant branch when stars have surface temperatures of roughly 5500 K. However, gravity waves can still spin down the upper radiative region, implying that the observed differential rotation is likely confined to the deep core near the hydrogen burning shell. The torque on the upper radiative region may also prevent the core from accreting high-angular momentum material and slow the rate of core spin-up. The observed spin-down of cores on the red giant branch cannot be totally attributed to gravity waves, but the waves may enhance shear within the radiative region and thus increase the efficacy of viscous/magnetic torques.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular momentum transport within evolved low-mass stars

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Formation and Gravitational-Wave Detection of Massive Stellar Black-Hole Binaries

If binaries consisting of two 100 Msun black holes exist they would serve as extraordinarily powerful gravitational-wave sources, detectable to redshifts of z=2 with the advanced LIGO/Virgo ground-based detectors. Large uncertainties about the evolution of massive stars preclude definitive rate predictions for mergers of these massive black holes. We show that rates as high as hundreds of detections per year, or as low as no detections whatsoever, are both possible. It was thought that the only way to produce these massive binaries was via dynamical interactions in dense stellar systems. This view has been challenged by the recent discovery of several stars with mass above 150 Msun in the R136 region of the Large Magellanic Cloud. Current models predict that when stars of this mass leave the main sequence, their expansion is insufficient to allow common envelope evolution to efficiently reduce the orbital separation. The resulting black-hole--black-hole binary remains too wide to be able to coalesce within a Hubble time. If this assessment is correct, isolated very massive binaries do not evolve to be gravitational-wave sources. However, other formation channels exist. For example, the high multiplicity of massive stars, and their common formation in relatively dense stellar associations, opens up dynamical channels for massive black hole mergers (e.g., via Kozai cycles or repeated binary-single interactions). We identify key physical factors that shape the population of very massive black-hole--black-hole binaries. Advanced gravitational-wave detectors will provide important constraints on the formation and evolution of very massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The fate of fallback matter around newly born compact objects

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Observational Signatures of Convectively Excited Internal Gravity Modes in Main Sequence Stars

We predict the flux and surface velocity perturbations produced by convectively excited gravity modes (g-modes) in main sequence stars. Core convection in massive stars can excite g-modes to sufficient amplitudes to be detectable with high precision photometry by Kepler and CoRoT, if the thickness of the convective overshoot region is < 30 per cent of a pressure scale height. The g-modes manifest as excess photometric variability, with amplitudes of ~ 10 micromagnitudes at frequencies ~ 10 microHz (0.8 1/d) near the solar metallicity zero-age main sequence. The flux variations are largest for stars with M > 5 solar masses, but are potentially detectable down to M ~ 2 - 3 solar masses. During the main sequence evolution, radiative damping decreases such that ever lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 100 micromagnitudes at the terminal-age main sequence. Using the same convective excitation model, we confirm previous predictions that solar g-modes produce surface velocity perturbations of < 0.3 mm/s. This implies that stochastically excited g-modes are more easily detectable in the photometry of massive main sequence stars than in the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak magnetic fields in early-type stars: failed fossils

Weak magnetic fields have recently been detected in Vega and Sirius. Here, we explore the possibility that these fields are the remnants of some field inherited or created during or shortly after star formation and, unlike true fossil fields, are still evolving as we observe them. The timescale of this evolution is given in terms of the Alfven timescale and the rotation frequency by tau_evol ~ tau_A^2 Omega, which would be comparable to the age of the star. It is shown that it is likely that all intermediate- and high-mass stars contain fields of at least the order of the strength found so far in Vega and Sirius. Faster rotators are expected to have stronger magnetic fields. Stars may experience an increase in field strength during their early main-sequence, but for most of their lives field strength will decrease slowly. The length scale of the magnetic structure on the surface may be small in very young stars but should quickly increase to at least very approximately a fifth of the stellar radius.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic spots on hot massive stars

Hot luminous stars show a variety of phenomena in their photospheres and winds which still lack clear physical explanation. Among these phenomena are photospheric turbulence, line profile variability (LPV), non-thermal emission, non-radial pulsations, discrete absorption components (DACs) and wind clumping. Cantiello et al. (2009) argued that a convection zone close to the stellar surface could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity associated with iron-group elements and is referred to as the "iron convection zone" (FeCZ). Assuming dynamo action producing magnetic fields at equipartition in the FeCZ, we investigate the occurrence of subsurface magnetism in OB stars. Then we study the surface emergence of these magnetic fields and discuss possible observational signatures of magnetic spots. Simple estimates are made using the subsurface properties of massive stars, as calculated in 1D stellar evolution models. We find that magnetic fields of sufficient amplitude to affect the wind could emerge at the surface via magnetic buoyancy. While at this stage it is difficult to predict the geometry of these features, we show that magnetic spots of size comparable to the local pressure scale height can manifest themselves as hot, bright spots. Localized magnetic fields could be widespread in those early type stars that have subsurface convection. This type of surface magnetism could be responsible for photometric variability and play a role in X-ray emission and wind clumping.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating Massive Main-Sequence Stars I: Grids of Evolutionary Models and Isochrones

We present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rotating Massive Main-Sequence Stars II: Simulating a Population of LMC early B-type Stars as a Test of Rotational Mixing

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stellar Evolution with Pulsation-Driven Mass Loss: The Case of LMC Cepheids

The Cepheid mass discrepancy, the difference between mass estimates from stellar pulsation and stellar evolution models, is a long standing challenge for the understanding of stellar astrophysics. We discuss the current state of the mass discrepancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass loss is a significant contributor to the mass discrepancy, but it is not clear if the metallicity dependence of Cepheid mass loss is consistent with the measured metallicity dependence of the mass discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Cepheid mass discrepancy and pulsation-driven mass loss

Context. A longstanding challenge for understanding classical Cepheids is the Cepheid mass discrepancy, where theoretical mass estimates using stellar evolution and stellar pulsation calculations have been found to differ by approximately 10 - 20%. Aims. We study the role of pulsation-driven mass loss during the Cepheid stage of evolution as a possible solution to this mass discrepancy. Methods. We computed stellar evolution models with a Cepheid mass-loss prescription and various amounts of convective core overshooting. The contribution of mass loss towards the mass discrepancy is determined using these models, Results. Pulsation-driven mass loss is found to trap Cepheid evolution on the instability strip, allowing them to lose about 5 - 10% of their total mass when moderate convective core overshooting, an amount consistent with observations of other stars, is included in the stellar models. Conclusions. We find that the combination of moderate convective core overshooting and pulsation-driven mass loss can solve the Cepheid mass discrepancy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convective Core Overshoot & Mass Loss in Classical Cepheids: A Solution to the Mass Discrepancy?

We explore the role of mass loss and convective core overshoot in the evolution of Classical Cepheids. Stellar evolution models are computed with a recipe for pulsation-driven mass loss and it is found that mass loss alone is unable to account for the long-standing Cepheid mass discrepancy. However, the combination of mass loss and moderate convective core overshooting does provide a solution, bringing the amount of convective core overshooting in Cepheids closer to that found for other stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolution of massive stars with pulsation-driven superwinds during the RSG phase

Pulsations driven by partial ionization of hydrogen in the envelope are often considered important for driving winds from red supergiants (RSGs). In particular, it has been suggested by some authors that the pulsation growth rate in a RSG can be high enough to trigger an unusually strong wind (or a super-wind), when the luminosity to mass ratio becomes sufficiently large. Using both hydrostatic and hydrodynamic stellar evolution models with initial masses ranging from 15 to 40 \Msun, we investigate 1) how the pulsation growth rate depends on the global parameters of supergiant stars, and 2) what would be the consequences of a pulsation-driven super-wind, if it occurred, for the late stages of massive star evolution. We suggest that such a super-wind history would be marked by a runaway increase, followed by a sudden decrease, of the winds mass loss rate. The impact on the late evolution of massive stars would be substantial, with stars losing a huge fraction of their H-envelope even with a significantly lower initial mass than previously predicted. This might explain the observed lack of Type II-P supernova progenitors having initial mass higher than about 17 \Msun. We also discuss possible implications for a subset of Type IIn supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermohaline mixing in evolved low-mass stars

Thermohaline mixing has recently been proposed to occur in low-mass red giants, with large consequence for the chemical yields of low-mass stars. We investigate the role of thermohaline mixing during the evolution of stars between 1Msun and 3Msun, in comparison to other mixing processes acting in these stars. We use a stellar evolution code which includes rotational mixing, internal magnetic fields and thermohaline mixing. We confirm that during the red giant stage, thermohaline mixing has the potential to decrease the abundance of ^3He which is produced earlier on the main sequence. In our models we find that this process is working on the RGB only in stars with initial mass M \simle 1.5Msun. Moreover we report that thermohaline mixing is present also during core helium burning and beyond, and has the potential to change the surface abundances of AGB stars. While we find rotational and magnetic mixing to be negligible compared to the thermohaline mixing in the relevant layers, the interaction of thermohaline motions with the differential rotation may be essential to establish the time scale of thermohaline mixing in red giants. To explain the surface abundances observed at the bump in the luminosity function, the speed of the mixing process needs to be more than two orders of magnitude higher than in our models. However it is not clear if thermohaline mixing is the only physical process responsible for these surface abundance anomalies. Therefore, at this stage, it is not possible to calibrate the efficiency of thermohaline mixing against the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulence and magnetic spots at the surface of hot massive stars

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

0 citations0 reviewsNo code linkedOpen access