Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

The imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, magnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here, we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious ‘spectral line’ is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4−4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.

Population synthesis of isolated neutron stars with magneto-rotational evolution

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with P0 ≲ 0.5 s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 0.6–0.7. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and γ-rays).

The many lives of magnetized neutron stars

The magnetic field strength at birth is arguably one of the most important properties to determine the evolutionary path of a neutron star. Objects with very high fields, collectively known as magnetars, are characterized by high X-ray quiescent luminosities, occurrence of outbursts, and, for some of them, sporadic giant flares. While the magnetic field strength is believed to drive their collective behaviour, however, the diversity of their properties, and, especially, the observation of magnetar-like bursts from “low-field” pulsars, has been a theoretical puzzle. In this review, we discuss results of long-term simulations following the coupled evolution of the X-ray luminosity and the timing properties for a large, homogeneous sample of X-ray emitting isolated neutron stars, accounting for a range of initial magnetic field strengths, envelope compositions, and neutron star masses. In addition, by following the evolution of magnetic stresses within the neutron star crust, we can also relate the observed magnetar phenomenology to the physical properties of neutron stars, and in particular to their age and magnetic field strength and topology. The dichotomy of “high-B” field pulsars versus magnetars is naturally explained, and occasional outbursts from old, low B-field neutron stars are predicted. We conclude by speculating on the fate of old magnetars, and by presenting observational diagnostics of the neutron star crustal field topology.

Corneal Biomechanics, Retinal Nerve Fiber Layer, and Optic Disc in Children

Purpose: To evaluate the possible associations between corneal biomechanical parameters, optic disc morphology, and retinal nerve fiber layer (RNFL) thickness in healthy white Spanish children. Methods: This cross-sectional study included 100 myopic children and 99 emmetropic children as a control group, ranging in age from 6 to 17 years. The Ocular Response Analyzer was used to measure corneal hysteresis (CH) and corneal resistance factor. The optic disc morphology and RNFL thickness were assessed using posterior segment optical coherence tomography (Cirrus HD-OCT). The axial length was measured using an IOLMaster, whereas the central corneal thickness was measured by anterior segment optical coherence tomography (Visante OCT). Results: The mean (±SD) age and spherical equivalent were 12.11 (±2.76) years and −3.32 (±2.32) diopters for the myopic group and 11.88 (±2.97) years and +0.34 (±0.41) diopters for the emmetropic group. In a multivariable mixed-model analysis in myopic children, the average RNFL thickness and rim area correlated positively with CH (p = 0.007 and p = 0.001, respectively), whereas the average cup-to-disc area ratio correlated negatively with CH (p = 0.01). We did not observe correlation between RNFL thickness and axial length (p = 0.05). Corneal resistance factor was only positively correlated with the rim area (p = 0.001). The central corneal thickness did not correlate with the optic nerve parameters or with RNFL thickness. These associations were not found in the emmetropic group (p > 0.05 for all). Conclusions: The corneal biomechanics characterized with the Ocular Response Analyzer system are correlated with the optic disc profile and RNFL thickness in myopic children. Low CH values may indicate a reduction in the viscous dampening properties of the cornea and the sclera, especially in myopic children.

Quiescent thermal emission from neutron stars in low-mass X-ray binaries

Context. We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in X-rays (outbursts). Aims. The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. Methods. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659−29, KS 1731−260, XTE J1701−462, EXO 0748−676  and IGR J17480−2446. Results. We find that the evolution of MXB 1659−29, KS 1731−260 and EXO 0748−676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701−462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480−2446 with an additional heat deposition in the outer crust from shallow sources.

The Gaia-ESO Survey: The analysis of high-resolution UVES spectra of FGK-type stars

Context. The ongoing Gaia-ESO Public Spectroscopic Survey is using FLAMES at the VLT to obtain high-quality medium-resolution Giraffe spectra for about 105 stars and high-resolution UVES spectra for about 5000 stars. With UVES, the Survey has already observed 1447 FGK-type stars. Aims. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO second internal release and will be part of its first public release of advanced data products. Methods. The final parameter scale is tied to the scale defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. In addition, a set of open and globular clusters is used to evaluate the physical soundness of the results. Each of the implemented methodologies is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted medians of those from the individual methods. Results. The recommended results successfully reproduce the atmospheric parameters of the benchmark stars and the expected Teff-log  g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log  g and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50−100 K for Teff, 0.10−0.25 dex for log  g and 0.05−0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex. Conclusions. The Gaia-ESO sample of high-resolution spectra of FGK-type stars will be among the largest of its kind analyzed in a homogeneous way. The extensive list of elemental abundances derived in these stars will enable significant advances in the areas of stellar evolution and Milky Way formation and evolution.

The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and α-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpcto 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages >9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more α-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.

Investigating the potential magnetic origin of wind variability in OB stars

In this thesis, the origin of large-scale structures in hot star winds, believed to be responsible for the presence of discrete absorption components (DACs) in the absorption troughs of ultraviolet resonance lines, is constrained using both observations and numerical simulations. These structures are understood as arising from bright regions on the stellar surface, although their physical cause remains unknown. First, we use high quality circular spectropolarimetric observations of 13 well-studied OB stars to evaluate the potential role of dipolar magnetic fields in producing DACs. We perform longitudinal field measurements and place limits on the field strength using Bayesian inference, assuming that it is dipolar. No magnetic field was detected within this sample. The derived constraints statistically refute any significant dynamical influence from a magnetic dipole on the wind for all of these stars, ruling out such fields as a cause for DACs. Second, we perform numerical simulations using bright spots constrained by broadband optical photometric observations. We calculate hydrodynamical wind models using three sets of spot sizes and strengths. Co-rotating interaction regions are yielded in each model, and radiative transfer shows that the properties of the variations in the UV resonance lines synthesized from these models are consistent with those found in observed UV spectra, establishing the first consistent link between UV spectroscopic line profile variability and photometric variations and thus supporting the bright spot paradigm (BSP). Finally, we develop and apply a phenomenological model to quantify the measurable effects co-rotating bright spots would have on broadband optical photometry and on the profiles of photopheric lines in optical spectra. This model can be used to evaluate the existence of these spots, and, in the event of their detection, characterize them. Furthermore, a tentative spot evolution model is presented. A preliminary analysis of its output, compared to the observed photometric variations of xi Persei, suggests the possible existence of “active longitudes” on the surface of this star. Future work will expand the range of observational diagnostics that can be interpreted within the BSP, and link phenomenology (bright spots) to physical processes (magnetic spots or non-radial pulsations).

Drusen do disco óptico e tilted disc

Trabalho Final do Curso de Mestrado Integrado em Medicina, Faculdade de Medicina, Universidade de Lisboa, 2014

Non-Viral Gene Therapy to Cells of the bovine Intervertebral Disc

Background. Low back pain is an increasing global health problem, which is associated with intervertebral disc (IVD) damage and degeneration. Major changes occur in the nucleus pulposus (NP), with the degradation of the extracellular matrix (ECM).1 Further studies showed that growth factors from transforming growth factor β (TGFβ) and bone morphogenic proteins (BMP) family may induce chondrogenic differentiation of mesenchymal stem cells (MSC).2 Focusing on non-viral gene therapies and their possible translation into the clinics, we investigated if GDF6 (syn. BMP13 or CDMP2) can induce regeneration of degraded NP. We hypothesized that IVD transfected with plasmid over-expressing GDF6 also up-regulates other NP- and chondrogenic cell markers and enhances ECM deposition. Methods. Bovine nucleus pulposus (bNPC) and annulus fibrosus cells (bAFC) were harvested from bovine coccygeal IVD. Primary cells were then electroporized with plasmid GDF6 (Origene, vector RG211366) by optimizing parameters using the Neon Transfection system (Life Technologies, Basel). After transfection, cells were cultured in 2D monolayer or 3D alginate beads for 7, 14 or 21 days. Transfection efficiency of pGDF6 was analyzed by immunohistochemistry and fluorescent microscopy. Cell phenotype was quantified by real-time RT-PCR. To test a non-viral gene therapy applied directly to 3D whole organ culture, coccygeal bovine IVDs were harvested as previously described. Bovine IVDs were transfected by injection of plasmid GDF6 into the center. Electroporation was performed with ECM830 Square Wave Electroporation System (Harvard Apparatus, MA) using 2-needle array electrode or tweezertrodes. 72 h after tranfection discs were fixed and cryosectioned and analyzed by immunofluorescence against GDF6. Results. RT-PCR and immunohistochemistry confirmed up-regulation of GFP and GDF6 in the primary bNPC/bAFC culture. The GFP-tagged GDF6 protein, however, was not visible, possibly due to failure of dimer formation as a result of fusion structure. Organ IVD culture transfection revealed GDF6 positive staining in the center of the disc using 2-needle array electrode. Results from tweezertrodes did not show any GDF6 positive cells. Conclusion. Non-viral transfection is an appealing approach for gene therapy as it fulfills the translational safety aspects of transiency and lacks the toxic effects of viral transduction. We identified novel parameters to successfully transfect primary bovine IVD cells. For transfection of whole IVD explants electroporation parameters need to be further optimized. Acknowledgements. This project was funded by the Lindenhof Foundation (Funds “Research & Teaching”) Project no. 13-02-F. The imaging part of this study was performed with the facility of the Microscopy Imaging Center (MIC), University of Bern. References. Roughly PJ (2004): Spine (Phila), 29:2691-2699 Clarke LE, McConell JC, Sherratt MJ, Derby B, Richardson SM, Hoyland JA (2014), Arthritis Research & Therapy, 16:R67

Compression Loading Induced Cellular Stress Response of Intervertebral Disc Cells in Organ Culture

Introduction: Mechanical stress is often associated to interverterbal disc (IVD) degeneration and the effect of mechanical loading on IVD has been studied and reviewed.1,2 Previously, expression of heat shock proteins, HSP70 and HSP27 has been found in pathological discs.3 However, there is no direct evidence on whether IVD cells respond to the mechanical loading by expression of HSPs. The objective of this study is to investigate the stress response of IVD cells during compressive loading in an organ culture. Materials and Methods: Fresh adult bovine caudal discs were cultured with compressive loading applied at physiological range. Effect of loading type (static and dynamic) and repeated loading (2 hours per day for 2 days) were studied. Nucleus pulposus (NP) and annulus fibrosus (AF) of the IVD were retrieved at different time points: right after loading and right after resting. Positive control discs were heat shocked (43°C). Cell activity was assessed and expression of stress response genes (HSP70 and HSF1) and matrix remodeling genes (ACAN, COL2, COL1, ADAMTS4, MMP3 and MMP13) were studied. Results: Cell activity was maintained in all groups. Both NP and AF expressed high level of HSP70 in heat shock groups, confirming their expression in response to stress. In NP, expression of HSP70 was up-regulated after static loading and dynamic loading with higher fold change was observed after static loading. During repeated loading, HSP70 appeared to be upregulated right after loading and decreased after resting. Such trend was not observed in AF and HSF1 levels. Expressions of matrix remodeling genes did not change significantly with loading except ADAMTS4 decreased in AF during static loading. Conclusion: This study demonstrated that NP cells upregulate expression of HSP70 in response to loading induced stress without changing cell activity and matrix remodeling significantly. Acknowledgments: This project was funded by AO Spine (AOSPN) (grant number: SRN_2011_14) and a fellowship exchange award by AO Spine Scientific Research Network (SRN).