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.

Outer crust of a cold non-accreting magnetar

The outer-crust structure and composition of a cold, non-accreting magnetar are studied. We model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb crystal embedded in an electron gas. The main effects of the strong magnetic field are those of quantizing the electron motion in Landau levels and of modifying the nuclear single-particle levels producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. The effect of a homogeneous and constant magnetic field on nuclear masses has been predicted by using a covariant density functional in which induced currents and axial deformation due to the presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently in the nucleon and meson equations of motion. Although not yet observed, for Ba 1016 G both effects contribute to produce different compositions - odd-mass nuclei are frequently predicted - and to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a regime, the magnetic-field effects on nuclei favor the appearance of heavier nuclei at low pressures. As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme magnetic field considered, B=1018 G, and for the models studied, almost the whole outer crust is made of 4092Zr52.

On neutron stars in ƒ (R) theories: Small radii, large masses and large energy emitted in a merger

In the context of ƒ (R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star’s edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic R^2 and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable ƒ (R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the ƒ (R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.

Search for chaos in neutron star systems: Is Cyg X-3 a black hole?

The accretion disk around a compact object is a nonlinear general relativistic system involving magnetohydrodynamics. Naturally, the question arises whether such a system is chaotic (deterministic) or stochastic (random) which might be related to the associated transport properties whose origin is still not confirmed. Earlier, the black hole system GRS 1915+105 was shown to be low-dimensional chaos in certain temporal classes. However, so far such nonlinear phenomena have not been studied fairly well for neutron stars which are unique for their magnetosphere and kHz quasi-periodic oscillation (QPO). On the other hand, it was argued that the QPO is a result of nonlinear magnetohydrodynamic effects in accretion disks. If a neutron star exhibits chaotic signature, then what is the chaotic/correlation dimension? We analyze RXTE/PCA data of neutron stars Sco X-1 and Cyg X-2, along with the black hole Cyg X-1 and the unknown source Cyg X-3, and show that while Sco X-1 and Cyg X-2 are low dimensional chaotic systems, Cyg X-1 and Cyg X-3 are stochastic sources. Based on our analysis, we argue that Cyg X-3 may be a black hole.

Diffusion of Pure CH4 and Its Binary Mixture with CO2 in Faujasite NaY: A Combination of Neutron Scattering Experiments and Molecular Dynamics Simulations

The self-diffusion properties of pure CH4 and its binary mixture with CO2 within MY zeolite have been investigated by combining an experimental quasi-elastic neutron scattering (QENS) technique and classical Molecular dynamics simulations. The QENS measurements carried out at 200 K led to an unexpected self-diffusivity profile for Pure CH4 with the presence of a maximum for a loading of 32 CH4/unit cell, which was never observed before for the diffusion of apolar species in azeolite system With large windows. Molecular dynamics simulations were performed using two distinct microscopic models for representing the CH4/NaY interactions. Depending on the model, we are able to fairly reproduce either the magnitude or the profile of the self-diffusivity.Further analysis allowed LIS to provide some molecular insight into the diffusion mechanism in play. The QENS measurements report only a slight decrease of the self-diffusivity of CH4 in the presence of CO2 when the CO2 loading increases. Molecular dynamics simulations successfully capture this experimental trend and suggest a plausible microscopic diffusion mechanism in the case of this binary mixture.

Magnetic and spin evolution of pulsars

We explore the consequences of the model of spin-down-induced flux expulsion for the magnetic field evolution in solitary as well as in binary neutron stars. The spin evolution of pulsars, allowing for their field evolution according to this model, is shown to be consistent with the existing observational constraints in both low- and high-mass X-ray binary systems. The contribution from pulsars recycled in massive binaries to the observed excess in the number of low-field (10(11)-10(12) G) solitary pulsars is argued to be negligible in comparison with that of normal pulsars undergoing a 'restricted' field decay predicted by the adopted field decay model. Magnetic fields of neutron stars born in close binaries with intermediate- or high-mass main-sequence companions are predicted to decay down to values as low as similar to 10(6) G, which would leave them unobservable as pulsars during most of their lifetimes. The post-recycling evolution of some of these systems can, however, account for the observed binary pulsars having neutron star or massive white dwarf companions. Pulsars recycled in the disc population low-mass binaries are expected to have residual fields greater than or similar to 10(8) G, while for those processed in globular clusters larger residual fields are predicted because of the lower field strength of the neutron star at the epoch of binary formation. A value of tau similar to 1-2 x 10(7) yr for the mean value of the Ohmic decay time-scale in the crusts of neutron stars is suggested, based on the consistency of the model predictions with the observed distribution of periods and magnetic fields in the single and binary pulsars.

New Insights from Aperiodic Variability of Young Stars

Nearly all young stars are variable, with the variability traditionally divided into two classes: periodic variables and aperiodic or "irregular" variables. Periodic variables have been studied extensively, typically using periodograms, while aperiodic variables have received much less attention due to a lack of standard statistical tools. However, aperiodic variability can serve as a powerful probe of young star accretion physics and inner circumstellar disk structure. For my dissertation, I analyzed data from a large-scale, long-term survey of the nearby North America Nebula complex, using Palomar Transient Factory photometric time series collected on a nightly or every few night cadence over several years. This survey is the most thorough exploration of variability in a sample of thousands of young stars over time baselines of days to years, revealing a rich array of lightcurve shapes, amplitudes, and timescales.

I have constrained the timescale distribution of all young variables, periodic and aperiodic, on timescales from less than a day to ~100 days. I have shown that the distribution of timescales for aperiodic variables peaks at a few days, with relatively few (~15%) sources dominated by variability on tens of days or longer. My constraints on aperiodic timescale distributions are based on two new tools, magnitude- vs. time-difference (Δm-Δt) plots and peak-finding plots, for describing aperiodic lightcurves; this thesis provides simulations of their performance and presents recommendations on how to apply them to aperiodic signals in other time series data sets. In addition, I have measured the error introduced into colors or SEDs from combining photometry of variable sources taken at different epochs. These are the first quantitative results to be presented on the distributions in amplitude and time scale for young aperiodic variables, particularly those varying on timescales of weeks to months.

The astrophysics of strongly interacting systems

This thesis presents investigations in four areas of theoretical astrophysics: the production of sterile neutrino dark matter in the early Universe, the evolution of small-scale baryon perturbations during the epoch of cosmological recombination, the effect of primordial magnetic fields on the redshifted 21-cm emission from the pre-reionization era, and the nonlinear stability of tidally deformed neutron stars.

In the first part of the thesis, we study the asymmetry-driven resonant production of 7 keV-scale sterile neutrino dark matter in the primordial Universe at temperatures T >~ 100 MeV. We report final DM phase space densities that are robust to uncertainties in the nature of the quark-hadron transition. We give transfer functions for cosmological density fluctuations that are useful for N-body simulations. We also provide a public code for the production calculation.

In the second part of the thesis, we study the instability of small-scale baryon pressure sound waves during cosmological recombination. We show that for relevant wavenumbers, inhomogenous recombination is driven by the transport of ionizing continuum and Lyman-alpha photons. We find a maximum growth factor less than ≈ 1.2 in 10⁷ random realizations of initial conditions. The low growth factors are due to the relatively short duration of the recombination epoch.

In the third part of the thesis, we propose a method of measuring weak magnetic fields, of order 10^-19 G (or 10^-21 G if scaled to the present day), with large coherence lengths in the inter galactic medium prior to and during the epoch of cosmic reionization. The method utilizes the Larmor precession of spin-polarized neutral hydrogen in the triplet state of the hyperfine transition. We perform detailed calculations of the microphysics behind this effect, and take into account all the processes that affect the hyperfine transition, including radiative decays, collisions, and optical pumping by Lyman-alpha photons.

In the final part of the thesis, we study the non-linear effects of tidal deformations of neutron stars (NS) in a compact binary. We compute the largest three- and four-mode couplings among the tidal mode and high-order p- and g-modes of similar radial wavenumber. We demonstrate the near-exact cancellation of their effects, and resolve the question of the stability of the tidally deformed NS to leading order. This result is significant for the extraction of binary parameters from gravitational wave observations.

Collective dynamics of glass-forming polymers at intermediate length scales: a synergetic combination of neutron scattering, atomistic simulations and theoretical modelling

Qens/wins 2014 - 11th International Conference on Quasielastic Neutron Scattering and 6th International Workshop on Inelastic Neutron Spectrometers / editado por:Frick, B; Koza, MM; Boehm, M; Mutka, H

Spectroscopy of neutron-rich Dy168,170: Yrast band evolution close to the NpNn valence maximum

The yrast sequence of the neutron-rich dysprosium isotope Dy168 has been studied using multinucleon transfer reactions following collisions between a 460-MeV Se82 beam and an Er170 target. The reaction products were identified using the PRISMA magnetic spectrometer and the γ rays detected using the CLARA HPGe-detector array. The 2+ and 4+ members of the previously measured ground-state rotational band of Dy168 have been confirmed and the yrast band extended up to 10+. A tentative candidate for the 4+→2+ transition in Dy170 was also identified. The data on these nuclei and on the lighter even-even dysprosium isotopes are interpreted in terms of total Routhian surface calculations and the evolution of collectivity in the vicinity of the proton-neutron valence product maximum is discussed. © 2010 The American Physical Society.

Temperature-stimulated abnormal annealing of neutron-induced damage in high-resistivity silicon detectors

Neutron-irradiated high-resistivity silicon detectors have been subjected to elevated temperature annealing (ETA). It has been found that both detector full depletion voltage and leakage current exhibit abnormal annealing (or ''reverse annealing'') behaviour for highly irradiated detectors: increase with ETA. Laser induced current measurements indicate a net increase of acceptor type space charges associated with the full depletion voltage increase after ETA. Current deep level transient spectroscopy (I-DLTS) and thermally stimulated current (TSC) data show that the dominant effect is the increase of a level at 0.39 eV below the conduction band (E(c) - 0.39 eV) or a level above the valence band (E(v) + 0.39 eV). Candidates tentatively identified for this level are the singly charged double vacancy (V-V-) level at E(c) - 0.39 eV, the carbon interstitial-oxygen interstitial (C-i-O-i) level at E(v) + 0.36 eV, and/or the tri-vacancy-oxygen center (V3O) at E(v) + 0.40 eV.

STUDY OF THE LONG-TERM STABILITY OF THE EFFECTIVE CONCENTRATION OF IONIZED SPACE CHARGES (N-EFF) OF NEUTRON-IRRADIATED SILICON DETECTORS FABRICATED BY VARIOUS THERMAL-OXIDATION PROCESSES

Experimental study of the reverse annealing of the effective concentration of ionized space charges (N-eff, also called effective doping or impurity concentration) of neutron irradiated high resistivity silicon detectors fabricated on wafers with various thermal oxides has been conducted at room temperature (RT) and elevated temperature (ET). Various thermal oxidations with temperatures ranging from 975 degrees C to 1200 degrees C with and without trichlorethane (TCA), which result in different concentrations of oxygen and carbon impurities, have been used. It has been found that, the RT annealing of the N-eff is hindered initially (t < 42 days after the radiation) for detectors made on the oxides with high carbon concentrations, and there was no carbon effect on the long term (t > 42 days after the radiation) N-eff reverse annealing. No apparent effect of oxygen on the stability of N-eff has been observed at RT. At elevated temperature (80 degrees C), no significant difference in annealing behavior has been found for detectors fabricated on silicon wafers with various thermal oxides. It is apparent that for the initial stages (first and/or second) of N-eff reverse annealing, there may tie no dependence on the oxygen and carbon concentrations in the ranges studied.