Source author record

Pablo Marchant

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

21 published item(s)

preprint2022arXiv

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos.

preprint2022arXiv

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars.

preprint2022arXiv

The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

preprint2021arXiv

Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts

Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances.

preprint2020arXiv

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars. I. The northern WC sequence

It is now well established that the majority of massive stars reside in multiple systems. However, the effect of multiplicity is not sufficiently understood, resulting in a plethora of uncertainties about the end stages of massive-star evolution. In order to investigate these uncertainties, it is useful to study massive stars just before their demise. Classical Wolf-Rayet stars represent the final end stages of stars at the upper-mass end. The multiplicity fraction of these stars was reported to be ${\sim}0.4$ in the Galaxy but no correction for observational biases has been attempted. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited ($V$ $\leq 12$) sample of Galactic Wolf-Rayet stars to derive their bias-corrected multiplicity properties. The present paper focuses on 12 northern Galactic carbon-rich (WC) Wolf-Rayet stars observable with the 1.2m Mercator telescope on the island of La Palma. We homogeneously measured relative radial velocities (RVs) for carbon-rich Wolf-Rayet stars using cross-correlation. Variations in the derived RVs were used to flag binary candidates. We investigated probable orbital configurations and provide a first correction of observational biases through Monte-Carlo simulations. Of the 12 northern Galactic WC stars in our sample, seven show peak-to-peak RV variations larger than 10 km s$^{-1}$, which we adopt as our detection threshold. This results in an observed spectroscopic multiplicity fraction of 0.58 with a binomial error of 0.14. In our campaign, we find a clear lack of short-period (P~$<~\sim$100\,d), indicating that a large number of Galactic WC binaries likely reside in long-period systems. Finally, our simulations show that at the 10% significance level, the intrinsic multiplicity fraction of the Galactic WC population is at least 0.72.

preprint2020arXiv

BAT99 126: A multiple Wolf-Rayet system in the Large Magellanic Cloud with a massive near-contact binary

BAT99 126 is a multiple system in the Large Magellanic Cloud containing a Wolf-Rayet (WR) star, which has a reported spectroscopic (orbital) period of 25.5 days and a photometric (orbital) period of 1.55 days, and hence is potentially one of the shortest WR binaries known to date. Such short-period binary systems containing a WR star in low-metallicity environments are prime candidate progenitors of black-hole (BH) mergers. By thoroughly analysing the spectroscopic and photometric data, we aim to establish the true multiplicity of BAT99 126, characterise the orbit(s) of the system, measure the physical properties of its individual components, and determine the overall evolutionary status of the system. Using newly acquired high resolution spectra, we measured radial velocities and performed a spectral analysis of the individual components. We estimated the age of the system and derived an evolutionary scenario for the 1.55-day system. BAT99 126 comprises at least four components. The 1.55-day photometric signal originates in an eclipsing O+O binary of spectral types O4 V and O6.5 V, and masses of 36+/-5 MSun and 15+/-2 MSun, respectively. For the WR star, we derived a spectral type WN2.5-3 and reject the previously reported 25.5-d period. We find evidence in the composite spectra of a fourth component, a B1 V star, which shows radial velocity variation. The configurations of the B-type star, WR star, and possible additional undetected components remain unknown. We estimated the age of the system to be 4.2 Myr. The O+O binary likely went through a phase of conservative mass transfer and is currently a near-contact system. We show that BAT99 126 is a multiple - quadruple or higher-order - system with a total initial mass > 160 MSun. The 1.55-day O+O binary most likely will not evolve towards a BH+BH merger, but instead will merge before the components collapse to BHs.

preprint2020arXiv

Effects of close binary evolution on the main-sequence morphology of young star clusters

Star clusters are the building blocks of galaxies. They are composed of stars of nearly equal age and chemical composition, allowing us to use them as chronometers and as testbeds for gauging stellar evolution. It has become clear recently that massive stars are formed preferentially in close binaries, in which mass transfer will drastically change the evolution of the stars. This is expected to leave a significant imprint in the distribution of cluster stars in the Hertzsprung-Russell diagram. Our results, based on a dense model grid of more than 50,000 detailed binary-evolution calculations, indeed show several distinct, coeval main-sequence (MS) components, most notably an extended MS turnoff region, and a group of near-critical rotating stars that is spread over a large luminosity range on the red side of the classical MS. We comprehensively demonstrate the time evolution of the features in an animation, and we derive analytic expressions to describe these features. We find quantitative agreement with results based on recent photometric and spectroscopic observations. We conclude that while other factors may also be at play, binary evolution has a major impact on the MS morphology of young star clusters.

preprint2020arXiv

Pulsational pair-instability supernovae in very close binaries

Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor ($Z_\odot/10$) single helium stars using the \texttt{MESA} software instrument. For PPISN progenitors with pre-pulse masses $>50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4\;$~yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.

preprint2020arXiv

The impact of stellar rotation on the black hole mass-gap from pair-instability supernovae

Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between $\sim 45M_\odot-120M_\odot$, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics it has become possible to test this prediction, and there is an important associated effort to understand what theoretical uncertainties modify the boundaries of this gap. In this work we study the impact of rotation on the hydrodynamics of PISNe, which leave no compact remnant, as well as the evolution of pulsational-PISNe (PPISNe), which undergo thermonuclear eruptions before forming a compact object. We perform simulations of non-rotating and rapidly-rotating stripped helium stars in a metal poor environment $(Z_\odot/50)$ in order to resolve the lower edge of the upper mass-gap. We find that the outcome of our simulations is dependent on the efficiency of angular momentum transport, with models that include efficient coupling through the Spruit-Tayler dynamo shifting the lower edge of the mass-gap upwards by $\sim 4\%$, while simulations that do not include this effect shift it upwards by $\sim 15\%$. From this, we expect the lower edge of the upper mass-gap to be dependent on BH spin, which can be tested as the number of observed BH mergers increases. Moreover, we show that stars undergoing PPISNe have extended envelopes ($R\sim 10-1000~R_\odot$) at iron-core collapse, making them promising progenitors for ultra-long gamma-ray bursts.

preprint2019arXiv

Eclipses of continuous gravitational waves as a probe of stellar structure

Although gravitational waves only interact weakly with matter, their propagation is affected by a gravitational potential. If a gravitational wave source is eclipsed by a star, measuring these perturbations provides a way to directly measure the distribution of mass throughout the stellar interior. We compute the expected Shapiro time delay, amplification and deflection during an eclipse, and show how this can be used to infer the mass distribution of the eclipsing body. We identify continuous gravitational waves from neutron stars as the best candidates to detect this effect. When the Sun eclipses a far-away source, depending on the depth of the eclipse the time-delay can change by up to $\sim 0.034$ ms, the gravitational-wave strain amplitude can increase by $\sim 4$%, and the apparent position of the source in the sky can vary by $4''$. Accreting neutron stars with Roche-lobe filling companion stars have a high probability of exhibiting eclipses, producing similar time delays but undetectable changes in amplitude and sky location. Even for the most rapidly rotating neutron stars, this time delay only corresponds to a few percent of the phase of the gravitational wave, making it an extremely challenging measurement. However, if sources of continuous gravitational waves exist just below the limit of detection of current observatories, next-generation instruments will be able to observe them with enough precision to measure the signal of an eclipsing star. Detecting this effect would provide a new direct probe to the interior of stars, complementing asteroseismology and the detection of solar neutrinos.

preprint2019arXiv

No signature of the orbital motion of a putative 70 solar mass black hole in LB-1

Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.

preprint2019arXiv

Progenitors of Type IIb Supernovae: I. Evolutionary Pathways and Rates

Type IIb supernovae (SNe) are important candidates to understand mechanisms that drive the stripping of stripped-envelope (SE) supernova (SN) progenitors. While binary interactions and their high incidence are generally cited to favor them as Type IIb SN progenitors, this idea has not been tested using models covering a broad parameter space. In this paper, we use non-rotating single- and binary-star models at solar and low metallicities spanning a wide parameter space in primary mass, mass ratio, orbital period, and mass transfer efficiencies. We find that our single- and binary-star models contribute to roughly equal, however small, numbers of Type IIb SNe at solar metallicity. Binaries only dominate as progenitors at low metallicity. We also find that our models can account for less than half the observationally inferred rate for Type IIb SNe at solar metallicity, with computed rates ~<4% of core-collapse (CC) SNe. On the other hand, our models can account for the rates currently indicated by observations at low metallicity, with computed rates as high as 15% of CC SNe. However, this requires low mass transfer efficiencies (~<0.1) to prevent most progenitors from entering contact. We suggest that the stellar wind mass-loss rates at solar metallicity used in our models are too high. Lower mass-loss rates would widen the parameter space for binary Type IIb SNe at solar metallicity by allowing stars that initiate mass transfer earlier in their evolution to reach CC without getting fully stripped.

preprint2016arXiv

A new route towards merging massive black holes

Recent advances in gravitational-wave astronomy make the direct detection of gravitational waves from the merger of two stellar-mass compact objects a realistic prospect. Evolutionary scenarios towards mergers of double compact objects generally invoke common-envelope evolution which is poorly understood, leading to large uncertainties in merger rates. We explore the alternative scenario of massive overcontact binary (MOB) evolution, which involves two very massive stars in a very tight binary which remain fully mixed due to their tidally induced high spin. We use the public stellar-evolution code MESA to systematically study this channel by means of detailed simulations. We find that, at low metallicity, MOBs produce double-black-hole (BH+BH) systems that will merge within a Hubble time with mass ratios close to one, in two mass ranges, ~25...60msun and >~ 130msun, with pair instability supernovae (PISNe) being produced in-between. Our models are also able to reproduce counterparts of various stages in the MOB scenario in the local Universe, providing direct support for it. We map the initial parameter space that produces BH+BH mergers, determine the expected chirp mass distribution, merger times, Kerr parameters and predict event rates. We typically find that for Z~<Z_sun/10, there is one BH+BH merger for ~1000 core-collapse supernovae. The advanced LIGO (aLIGO) detection rate is more uncertain and depends on the metallicity evolution. Deriving upper and lower limits from a local and a global approximation for the metallicity distribution of massive stars, we estimate aLIGO detection rates (at design limit) of ~19-550 yr^(-1) for BH+BH mergers below the PISN gap and of ~2.1-370 yr^(-1) above the PISN gap. Even with conservative assumptions, we find that aLIGO should soon detect BH+BH mergers from the MOB scenario and that these could be the dominant source for aLIGO detections.

preprint2016arXiv

Models of low-mass helium white dwarfs including gravitational settling, thermal and chemical diffusion, and rotational mixing

A large number of extremely low-mass helium white dwarfs (ELM WDs) have been discovered in recent years. The majority of them are found in close binary systems suggesting they are formed either through a common-envelope phase or via stable mass transfer in a low-mass X-ray binary (LMXB) or a cataclysmic variable (CV) system. Here, we investigate the formation of these objects through the LMXB channel with emphasis on the proto-WD evolution in environments with different metallicities. We study, for the first time, the combined effects of rotational mixing and element diffusion (e.g. gravitational settling, thermal and chemical diffusion) on the evolution of proto-WDs and on the cooling properties of the resulting WDs. We present state-of-the-art binary stellar evolution models computed with MESA for metallicities between Z=0.0002 and Z=0.02, producing WDs with masses between 0.16-0.45 M$_{\odot}$. Our results confirm that element diffusion plays a significant role in the evolution of proto-WDs that experience hydrogen shell flashes. The occurrence of these flashes produces a clear dichotomy in the cooling timescales of ELM WDs, which has important consequences e.g. for the age determination of binary millisecond pulsars. Rotational mixing is found to counteract the effect of gravitational settling in the surface layers of young, bloated ELM proto-WDs and therefore plays a key role in determining their surface chemical abundances. We predict that these proto-WDs have helium-rich envelopes through a significant part of their lifetime, a crucial ingredient for understanding the newly observed ELM proto-WD pulsators. The hydrogen envelope at detachment, although small compared to the total mass of the WD, contains enough angular momentum such that the spin frequency of the resulting WD on the cooling track is well above the orbital frequency.

preprint2015arXiv

AM Canum Venaticorum Progenitors with Helium Star Donors and the resultant Explosions

We explore the outcome of mass transfer via Roche lobe overflow (RLOF) of $M_{\rm He}\lesssim0.51 M_\odot$ pure helium burning stars in close binaries with white dwarfs (WDs). The evolution is driven by the loss of angular momentum through gravitational wave radiation (GWR), and both stars are modeled using Modules for Experiments in Stellar Astrophysics (MESA). The donors have masses of $M_{\rm He}=0.35, 0.4, \&\ 0.51M_\odot$ and accrete onto WDs of mass $M_{\rm WD}$ from $0.6M_\odot$ to $1.26M_\odot$. The initial orbital periods ($P_{\rm{orb}}$) span 20 to 80 minutes. For all cases, the accretion rate onto the WD is below the stable helium burning range, leading to accumulation of helium followed by unstable ignition. The mass of the convective core in the donors is small enough so that the WD accretes enough helium-rich matter to undergo a thermonuclear runaway in the helium shell before any carbon-oxygen enriched matter is transferred. The mass of the accumulated helium shell depends on $M_{\rm WD}$ and the accretion rate. We show that for $M_{\rm He}\gtrsim0.4 M_\odot$ and $M_{\rm WD}\gtrsim0.8 M_\odot$, the first flash is likely vigorous enough to trigger a detonation in the helium layer. These thermonuclear runaways may be observed as either faint and fast .Ia SNe, or, if the carbon in the core is also detonated, Type Ia SNe. Those that survive the first flash and eject mass will have a temporary increase in orbital separation, but GWR drives the donor back into contact, resuming mass transfer and triggering several subsequent weaker flashes.

preprint2015arXiv

Testing eccentricity pumping mechanisms to model eccentric long period sdB binaries with MESA

Hot subdwarf-B stars in long-period binaries are found to be on eccentric orbits, even though current binary-evolution theory predicts those objects to be circularised before the onset of Roche-lobe overflow (RLOF). We aim to find binary-evolution mechanisms that can explain these eccentric long-period orbits, and reproduce the currently observed period-eccentricity diagram. Three different processes are considered; tidally-enhanced wind mass-loss, phase-dependent RLOF on eccentric orbits and the interaction between a circumbinary disk and the binary. The binary module of the stellar-evolution code MESA (Modules for Experiments in Stellar Astrophysics) is extended to include the eccentricity-pumping processes. The effects of different input parameters on the final period and eccentricity of a binary-evolution model are tested with MESA. The end products of models with only tidally-enhanced wind mass-loss can indeed be eccentric, but these models need to lose too much mass, and invariably end up with a helium white dwarf that is too light to ignite helium. Within the tested parameter space, no sdBs in eccentric systems are formed. Phase-dependent RLOF can reintroduce eccentricity during RLOF, and could help to populate the short-period part of the period-eccentricity diagram. When phase-dependent RLOF is combined with eccentricity pumping via a circumbinary disk, the higher eccentricities can be reached as well. A remaining problem is that these models favour a distribution of higher eccentricities at lower periods, while the observed systems show the opposite. The models presented here are potentially capable of explaining the period-eccentricity distribution of long-period sdB binaries, but further theoretical work on the physical mechanisms is necessary.

preprint2015arXiv

Tidally-driven Roche-Lobe Overflow of Hot Jupiters with MESA

Many exoplanets have now been detected in orbits with ultra-short periods, very close to the Roche limit. Building upon our previous work, we study the possibility that mass loss through Roche lobe overflow (RLO) may affect the evolution of these planets, and could possibly transform a hot Jupiter into a lower-mass planet (hot Neptune or super-Earth). We focus here on systems in which the mass loss occurs slowly ("stable mass transfer" in the language of binary star evolution) and we compute their evolution in detail with the binary evolution code MESA. We include the effects of tides, RLO, irradiation and photo-evaporation of the planet, as well as the stellar wind and magnetic braking. Our calculations all start with a hot Jupiter close to its Roche limit, in orbit around a sun-like star. The initial orbital decay and onset of RLO are driven by tidal dissipation in the star. We confirm that such a system can indeed evolve to produce lower-mass planets in orbits of a few days. The RLO phase eventually ends and, depending on the details of the mass transfer and on the planetary core mass, the orbital period can remain around a few days for several Gyr. The remnant planets have a rocky core and some amount of envelope material, which is slowly removed via photo-evaporation at nearly constant orbital period; these have properties resembling many of the observed super-Earths and sub-Neptunes. For these remnant planets we also predict an anti-correlation between mass and orbital period; very low-mass planets ($M_{\rm pl}\,\lesssim\,5\,M_{\oplus}$) in ultra-short periods ($P_{\rm orb}$<1d) cannot be produced through this type of evolution.

preprint2014arXiv

Stability of Hall equilibria in neutron star crusts

In the solid crusts of neutron stars, the advection of the magnetic field by the current-carrying electrons, an effect known as Hall drift, should play a very important role as the ions remain essentially fixed (as long as the solid does not break). Although Hall drift preserves the magnetic field energy, it has been argued that it may drive a turbulent cascade to scales at which Ohmic dissipation becomes effective, allowing a much faster decay in objects with very strong fields. On the other hand, it has been found that there are "Hall equilibria", i.e., field configurations that are unaffected by Hall drift. Here, we address the crucial question of the stability of these equilibria through axially symmetric (2D) numerical simulations of Hall drift and Ohmic diffusion, with the simplifying assumption of uniform electron density and conductivity. We demonstrate the 2D-stability of a purely poloidal equilibrium, for which Ohmic dissipation makes the field evolve towards an attractor state through adjacent stable configurations, around which damped oscillations occur. For this field, the decay scales with the Ohmic timescale. We also study the case of an unstable equilibrium consisting of both poloidal and toroidal field components that are confined within the crust. This field evolves into a stable configuration, which undergoes damped oscillations superimposed on a slow evolution towards an attractor, just as the purely poloidal one.

preprint2013arXiv

Magnetohydrodynamic equilibria in barotropic stars

Although barotropic matter does not constitute a realistic model for magnetic stars, it would be interesting to confirm a recent conjecture that states that magnetized stars with a barotropic equation of state would be dynamically unstable (Reisenegger 2009). In this work we construct a set of barotropic equilibria, which can eventually be tested using a stability criterion. A general description of the ideal MHD equations governing these equilibria is summarized, allowing for both poloidal and toroidal magnetic field components. A new finite-difference numerical code is developed in order to solve the so-called Grad-Shafranov equation describing the equilibrium of these configurations, and some properties of the equilibria obtained are briefly discussed.

preprint2013arXiv

Stability of magnetic fields in non-barotropic stars: an analytic treatment

Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are known to persist for timescales comparable to their lifetimes. From a theoretical perspective this is problematic, as it can be shown that simple magnetic field configurations are always unstable. In non-barotropic stars, stable stratification allows for a much wider range of magnetic field structures than in barotropic stars, and helps stabilize them by making it harder to induce radial displacements. Recent simulations by Braithwaite and collaborators have shown that, in stably stratified stars, random initial magnetic fields evolve into nearly axisymmetric configurations with both poloidal and toroidal components, which then remain stable for some time. It is desirable to provide an analytic study of the stability of such fields. We write an explicit expression for a plausible equilibrium structure of an axially symmetric magnetic field with both poloidal and toroidal components of adjustable strengths, in a non-barotropic static fluid star, and study its stability using the energy principle. We construct a displacement field that should be a reasonable approximation to the most unstable mode of a toroidal field, and confirm Braithwaite's result that a given toroidal field can be stabilized by a poloidal field containing much less energy than the former. This is consistent with the speculation that the toroidal field is the main reservoir powering magnetar activity. The deformation of a neutron star caused by the hidden toroidal field can also cause emission of gravitational waves.

preprint2011arXiv

Revisiting the Flowers-Ruderman instability of magnetic stars

In 1977, Flowers and Ruderman described a perturbation that destabilises a purely dipolar magnetic field in a fluid star. They considered the effect of cutting the star in half along a plane containing the symmetry axis and rotating each half by $90\degr$ in opposite directions, which would cause the energy of the magnetic field in the exterior of the star to be greatly reduced, just as it happens with a pair of aligned magnets. We formally solve for the energy of the external magnetic field and check that it decreases monotonously along the entire rotation. We also describe the instability using perturbation theory, and see that it happens due to the work done by the interaction of the magnetic field with surface currents. Finally, we consider the stabilising effect of adding a toroidal field by studying the potential energy perturbation when the rotation is not done along a sharp cut, but with a continuous displacement field that switches the direction of rotation across a region of small but finite width. Using these results, we estimate the relative strengths of the toroidal and poloidal field needed to make the star stable to this displacement and see that the energy of the toroidal field required for stabilisation is much smaller than the energy of the poloidal field. We also show that, contrary to a common argument, the Flowers-Ruderman instability cannot be applied many times in a row to reduce the external magnetic energy indefinitely.