Pulsar timing irregularities and the imprint of magnetic field evolution

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ~ 3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n < 3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n > 3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (≳ 1014 G) multipolar or toroidal components, even with a weak (~1012 G) dipolar field are consistent with the observed trend of the timing properties.

The outburst decay of the low magnetic field magnetar SGR 0418+5729

We report on the long-term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729 using all the available X-ray data obtained with RXTE, Swift, Chandra, and XMM-Newton observations from the discovery of the source in 2009 June up to 2012 August. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: ˙ P = 4(1) × 10−15 s s−1, significant at a ∼3.5σ confidence level. This leads to a surface dipolar magnetic field of Bdip 6 × 1012 G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ∼1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ∼0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about three orders of magnitude: from about 1.4 × 10−11 to 1.2 × 10−14 erg s−1 cm−2. Deep radio, millimeter, optical, and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ∼550 kyr, and a dipolar magnetic field at birth of ∼1014 G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such a result in several other astrophysical contexts.

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.

Central compact objects and the hidden magnetic field scenario

Central compact objects (CCOs) are X-ray sources lying close to the centre of supernova remnants, with inferred values of the surface magnetic fields significantly lower (≲1011 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and re-emergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 10−4–10−3 M⊙ over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field re-emerges on a typical time-scale of 1–100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as an alternative to the antimagnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

Modelling of the surface emission of the low magnetic field magnetar SGR 0418+5729

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6×1012G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. None the less, the modelling allows us to place constraints on the system geometry (i.e. the angles ψ and ξ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to 65° ≲ ψ + ξ ≲ 125° or 235° ≲ ψ + ξ ≲ 295°; (b) the temperature contrast between the poles and the equator must be at least a factor of ∼6, and (c) the size of the hottest region ranges between 0.2 and 0.7 km (including uncertainties on the source distance). Lastly, we interpret our findings within the context of internal and external heating models.

Isotopic composition of carbonates from the Logachev Hydrothermal Field, MAR 14°45'N

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45' N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a d44/40Ca_HydEnd of -0.95+/-0.07 per mil relative to seawater (SW) and a 87Sr/86Sr isotope ratio of 0.7034(4). We demonstrate that d44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a d44/40Ca of -1.17+/-0.04 per mil (SW) for the host-rocks in the reaction zone and -1.45+/-0.05 per mil (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed d44/40Ca for Bulk Earth of -0.92+/-0.18 per mil (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta 61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta 62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about D44/40Ca = -0.5 per mil relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average d44/40Ca of -1.54+/-0.08 per mil (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.

(Table T1) Permeability, electrical resistivity and density of ODP Hole 193-1188A and Site 193-1189 minicores, PACMANUS field

Magmatic fluids, heat fluxes, and fluid/rock interactions associated with hydrothermal systems along spreading centers and convergent margins have a significant impact on the genesis of major sulfide deposits and biological communities. Circulation of hydrothermal fluids is one of the most fundamental processes associated with localized mineralization and is controlled by inherent porous and permeable properties of the ocean crust. Heat from magmatic intrusions drives circulation of seawater through permeable portions of the oceanic crust and upper mantle, discharging at the seafloor as both focused high-temperature (250°-400°C) fluids and diffuse lower-temperature (<250°C) fluids. This complex interaction between the circulating hydrothermal fluids and the oceanic basement greatly influences the physical properties and the composition of the crust (Thompson, 1983; Jacobson, 1992, doi:10.1029/91RG02811; Johnson and Semyan, 1994, doi:10.1029/93JB00717). During Ocean Drilling Program (ODP) Leg 193, 13 holes were drilled in the PACMANUS hydrothermal system (Binns, Barriga, Miller, et al., 2002, doi:10.2973/odp.proc.ir.193.2002). The hydrothermal system consists of isolated hydrothermal deposits lined along the main crest of the Pual Ridge, a 500- to 700-m-high felsic neovolcanic ridge in the eastern Manus Basin. The principal drilling targets were the Snowcap (Site 1188) and Roman Ruins (Site 1189) active hydrothermal fields. Samples from these two sites were used for a series of permeability, electrical resistivity, and X-ray computed tomography measurements.

Concentrations and isotopic compositions of low molecular weight hydrocarbons in hydrothermal fluids from ODP Field, Middle Valley, northern Juan de Fuca Ridge

The presence of sedimentary organic matter blanketing midocean ridge crests has a potentially strong impact on metal transport in hydrothermal vent fluids. To constrain the role of organic matter in metal mobility during hydrothermal sediment alteration, we reacted organic-rich diatomaceous ooze from Guaymas Basin, Gulf of California, and organic-poor hemipelagic mud from Middle Valley, northern Juan de Fuca Ridge, with seawater and a Na-Ca-K-Cl fluid of seawater chlorinity, at 275° to 400°C, 350 to 500 bars, and initial fluid: sediment mass ratios ranging from 1.6 to 9.8. Reaction of these fluids with both sediment types released CO2 and high concentrations of ore-forming metals (Fe, Mn, Zn, Pb) to solution. Relatively low concentrations of Cu were observed in solution and likely reflect the reducing conditions that resulted from the presence of sedimentary organic matter. Both the concentrations of CO2 and dissolved metals were lower in fluids reacted with Middle Valley sediment compared with aqueous concentrations in fluids reacted with Guaymas Basin sediment. During alteration of both sediment types, metal concentrations varied strongly as a function of temperature, increasing by up to an order of magnitude over the 75°C range of each experiment. Major element fluid chemistry and observed alteration assemblages suggest that during hydrothermal alteration of organic-lean sediment from Middle Valley a feldspar-quartz-illite mineral assemblage buffered in situ pH. In contrast, data from the experimental alteration of organic-rich Guaymas Basin sediment suggest that a calcite-plagioclase-quartz assemblage regulated in situ pH. Fluid speciation calculations suggest that in situ pH during Guaymas Basin sediment alteration was lower than during alteration of Middle Valley sediment and accounts for the substantially greater metal mobility at a given temperature and pressure during the former experiment. Comparison of our results with the results of basalt alteration experiments indicate that except for Cu, hydrothermal sediment alteration results in equal or greater concentrations of ore-forming metals at a given temperature and pressure. Accordingly, the presence of ore-forming metals in fluids currently venting from sediment-covered hydrothermal systems at concentrations substantially lower than in fluids from bare-rock systems may reflect chemical reequilibration during subsurface cooling within the sediment pile.

(Table 1) Sea surface temperature reconstruction for eastern equatorial Pacific surface sediment samples

Significant uncertainties persist in the reconstruction of past sea surface temperatures in the eastern equatorial Pacific, especially regarding the amplitude of the glacial cooling and the details of the post-glacial warming. Here we present the first regional calibration of alkenone unsaturation in surface sediments versus mean annual sea surface temperatures (maSST). Based on 81 new and 48 previously published data points, it is shown that open ocean samples conform to established global regressions of Uk'37 versus maSST and that there is no systematic bias from seasonality in the production or export of alkenones, or from surface ocean nutrient concentrations or salinity. The flattening of the regression at the highest maSSTs is found to be statistically insignificant. For the near-coastal Peru upwelling zone between 11-15°S and 76-79°W, however, we corroborate earlier observations that Uk'37 SST estimates significantly over-estimate maSSTs at many sites. We posit that this is caused either by uncertainties in the determination of maSSTs in this highly dynamic environment, or by biasing of the alkenone paleothermometer toward El Niño events as postulated by Rein et al. (2005).

Temperature gradient and thermal stresses in bodies with uniformly distributed volume heat sources.

Mode of access: Internet.

An isotopic analysis of geothermal brine and calcite scaling from the Blue Mountain geothermal field, Winnemucca, Nevada

Understanding, and controlling, the conditions under which calcite precipitates within geothermal energy production systems is a key step in maintaining production efficiency. In this study, I apply methods of bulk and clumped isotope thermometry to an operating geothermal energy facility in northern Nevada to see how those methods can better inform the facility owner, AltaRock Energy, Inc., about the occurrence of calcite scale in their power plant. I have taken water samples from five production wells, the combined generator effluent, shallow cold-water wells, monitoring wells, and surface water. I also collected calcite scale samples from within the production system. Water samples were analyzed for stable oxygen isotope composition (d18O). Calcite samples were analyzed for stable oxygen and carbon (d13C) composition, and clumped isotope composition (D47). With two exceptions, the water compositions are very similar, likely indicating common origin and a well-mixed hydrothermal system. The calcite samples are likewise similar to one another. Apparent temperatures calculated from d18O values of water and calcite are lower than those recorded for the system. Apparent temperatures calculated from D47 are several degrees higher than the recorded well temperatures. The lower temperatures from the bulk isotope data are consistent with temperatures that could be expected during a de-pressurization of the production system, which would cause boiling in the pipes, a reduction in system temperature, and rapid precipitation of calcite scale. However, the high apparent temperature indicated by the D47 data suggests that the calcite is depleted in clumped isotopes given the known temperature of the system, which is inconsistent with this hypothesis. This depletion could instead result from disequilibrium isotopic fractionation during the aforementioned boil events, which would make both the apparent d18O-based and D47-based temperatures unrepresentative of the actual water temperature. This research can help improve our understanding of how isotopic analyses can better inform us about the movement of water through geothermal systems of the past and how it now moves through modern systems. Increased understanding of water movement in these systems could potentially allow for more efficient utilization of geothermal energy as a renewable resource.

The electric field effect on spin injection in tunneling regime

Using the quantum tunneling theory, we investigate the spin-dependent transport properties of the ferromagnetic metal/Schottky barrier/semiconductor heterojunction under the influence of an external electric field. It is shown that increasing the electric field, similar to increasing the electron density in semiconductor, will result in a slight enhancement of spin injection in tunneling regime, and this enhancement is significantly weakened when the tunneling Schottky barrier becomes stronger. Temperature effect on spin injection is also discussed. (C) 2003 Elsevier B.V. All rights reserved.

[0 0 1] zone-axis bright-field diffraction contrast from coherent Ge(Si) islands on Si(0 0 1): Dedicated to Professor Fang-hua Li on the occasion of her 70th birthday

Coherent Ge(Si)/Si(001) quantum dot islands grown by solid source molecular beam epitaxy at a growth temperature of 700degreesC were investigated using transmission electron microscopy working at 300 kV. The [001] zone-axis bright-field diffraction contrast images of the islands show strong periodicity with the change of the TEM sample substrate thickness and the period is equal to the effective extinction distance of the transmitted beam. Simulated images based on finite element models of the displacement field and using multi-beam dynamical diffraction theory show a high degree of agreement. Studies for a range of electron energies show the power of the technique for investigating composition segregation in quantum dot islands. (C) 2003 Elsevier B.V. All rights reserved.