Neutron Stars—Thermal Emitters

Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main properties of the observational data. In particular, we focus on the nearby sources for which a clear thermal component has been detected, without being contaminated by other emission processes (magnetosphere, accretion, nebulae). We also discuss the applications of the modern theoretical models of the formation of spectra of strongly magnetized neutron stars to the observed thermally emitting objects.

Neutron Stars—Cooling and Transport

Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.

On the implication of the recently discovered 5 millisecond binary pulsar PSR 1855+09

e argue that the extraordinary fact that all three known millisecond pulsars are very close to the galactic plane implies that there must be ~100 potentially observable millisecond pulsars within ~4 kpc from the Sun. Our other main conclusion is that the dipole magnetic fields or old neutron stars probably saturate around 5 x 108 gauss.

Predicting spin of compact objects from their QPOs: A global QPO model

We establish a unified model to explain Quasi-Periodic-Oscillation (QPO) observed from black hole and neutron star systems globally. This is based on the accreting systems thought to be damped harmonic oscillators with higher order nonlinearity. The model explains multiple properties parallelly independent of the nature of the compact object. It describes QPOs successfully for several compact sources. Based on it, we predict the spin frequency of the neutron star Sco X-1 and the specific angular momentum of black holes GRO J1655-40, GRS 1915+105.

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.

The importance of spin for observing gravitational waves from coalescing compact binaries with LIGO and Virgo

General Relativity predicts the existence of gravitational waves, which carry information about the physical and dynamical properties of their source. One of the many promising sources of gravitational waves observable by ground-based instruments, such as in LIGO and Virgo, is the coalescence of two compact objects (neutron star or black hole). Black holes and neutron stars sometimes form binaries with short orbital periods, radiating so strongly in gravitational waves that they coalesce on astrophysically short timescales. General Relativity gives precise predictions for the form of the signal emitted by these systems. The most recent searches for theses events used waveform models that neglected the effects of black hole and neutron star spin. However, real astrophysical compact objects, especially black holes, are expected to have large spins. We demonstrate here a data analysis infrastructure which achieves an improved sensitivity to spinning compact binaries by the inclusion of spin effects in the template waveforms. This infrastructure is designed for scalable, low-latency data analysis, ideal for rapid electromagnetic followup of gravitational wave events.

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.

Emissividade de neutrinos pelo processo Urca direto na presença de um campo magnético

As estrelas de nêutrons nascem com altas temperaturas (~ 1011 K) e durante alguns segundos sofrem um rápido resfriamento por emissão de neutrinos. O processo Urca direto é o principal mecanismo para explicar essa perda de energia. O problema do resfriamento das estrelas de nêutrons é um problema de grande interesse porque seu entendimento pode fornecer informações importantes sobre a constituição do interior da estrela. Na literatura existente até o momento, a emissividade de neutrinos é calculada considerando os núcleons como partículas não relativísticas quando considerados todos os níveis de Landau das partículas carregadas. Por outro lado, a emissividade de neutrinos para núcleons relativísticos é calculada considerando somente o primeiro nível de Landau (para campo magnético forte). Para campos magnéticos fracos, onde mais de um nível de Landau é ocupado, é usada a emissividade correspondente à do campo nulo. Neste trabalho aplicamos a teoria de Weinberg-Salan para interações fracas no cálculo da emissividade de neutrinos com e sem campo magnético presente, num cálculo totalmente relativístico para os núcleons e considerando todos os níveis de Landau. Esta é a contribuição original do trabalho. Para descrever a matéria a altas densidades, utilizamos uma teoria relativística de campo médio a temperatura zero que inclui apenas o octeto bariônico e os léptons mais leves. São apresentados os resultados para a emissividade de neutrinos, onde é evidente a ocupação dos diferentes níveis de Landau como função do campo magnético.

Superfluidity in dense nuclear matter

We discuss the onset of superfluidity in neutron stars, where the model of nuclear matter is realized in a high-density and asymmetry state. In particular, we present the study of the effects of microscopic three-body forces on the proton pairing in the 1S0 channel and neutron pairing in 3PF1 channel for β-stable neutron star matter. It is found that the main effects of three-body forces are to shrink the domain of existence of the 1S0 below the threshold of the direct URCA process and to stretch the density range of the 3PF1 pairing in a broad domain so to cover most part of the neutron-star core.

Isospin dependence of nucleon (PF2)-P-3 superfluidity in asymmetric nuclear matter

We have investigated the isospin dependence of the neutron and proton (PF2)-P-3 superfluidity in isospin-asymmetric nuclear matter within the framework of the Brueckner-Hartree-Fock approach and the BCS theory. We show that the (PF2)-P-3 neutron and proton pairing gaps depend sensitively on isospin asymmetry of asymmetric nuclear matter. As the isospin asymmetry increases, the neutron (PF2)-P-3 superfluidity becomes stronger and the peak value of the neutron (PF2)-P-3 pairing gap increases rapidly. The isospin dependence of the proton (PF2)-P-3 superfluidity is shown to be opposite to the neutron one. The proton (PF2)-P-3 superfluidity becomes weaker at a higher asymmetry and it even vanishes at high enough asymmetries. At high asymmetries, the neutron (PF2)-P-3 superfluidity turns out to be much stronger than the proton one, implying that the neutron (PF2)-P-3 superfluidity is dominated in the highly asymmetric dense interior of neutron stars.

Hot Nuclear Matter Equation of State and Finite Temperature Kaon Condensation

We perform a systematic calculation of the equation of state of asymmetric nuclear matter at finite temperature within the framework of the Brueckner-Hartree-Fock approach with a microscopic three-body force. When applying it to the study of hotka on condensed matter, we find that the thermal effect is more profound in comparison with normal matter, in particular around the threshold density. Also, the increase of temperature makes the equation of state slightly stiffer through suppression of kaon condensation.

System size and beam energy effects on probing the high-density behavior of nuclear symmetry energy with pion ratio

Based on the isospin-and momentum-dependent hadronic transport model IBUU04, we have investigated the pi(-)/pi(+) ratio in the following three reactions: Ca-48+Ca-48, Sn-124 +Sn-124 and Au-197+Au-197 with nearly the same isospin asymmetry but different masses, at the bombarding energies from 0.25 to 0.6 A GeV. It is shown that the sensitivity of probing the E-sym (rho) with pi(-)/pi(+) increases with increasing the system size or decreasing the beam energy, showing a correlation to the degree of isospin fractionation. Therefore, with a given isospin asymmetry, heavier system at energies near the pion threshold is preferential to study the behavior Of nuclear symmetry energy at supra-saturation densities.

Effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) ratios in heavy-ion collisions at intermediate energies

Based on the isospin-dependent Boltzmann-Uehling-Uhlenbeck transport model and the scaling model according to nucleon effective mass, effects of elastic and inelastic NN scattering cross sections on pi(-)/pi(+) in the neutron-rich reaction of Ca-48 + Ca-48 at a beam energy of 400 MeV/nucleon are studied. It is found that cross-section effects of both NN elastic and inelastic scatterings affect Delta(1232), pi(-) and pi(+) production, as well as the value of pi(-)/pi(+).

Dynamics of strangeness production in heavy-ion collisions near threshold energies

Within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model, the dynamics of strangeness (K-0,K-+, Lambda, and Sigma(-,0,+)) production in heavy-ion collisions near threshold energies is investigated systematically, with the strange particles considered to be produced mainly by inelastic collisions of baryon-baryon and pion-baryon. Collisions in the region of suprasaturation densities of the dense baryonic matter formed in heavy-ion collisions dominate the yields of strangeness production. Total multiplicities as functions of incident energies and collision centralities are calculated with the Skyrme parameter SLy6. The excitation function of strangeness production is analyzed and also compared with the KaoS data for K+ production in the reactions C-12 + C-12 and Au-197 + Au-197.

Can the nuclear symmetry potential at supra-saturation densities be negative?

In the framework of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, for the central Au-197 + Au-197 reaction at an incident beam energy of 400 MeV/nucleon, the effect of nuclear symmetry potential at supra-saturation densities on the preequilibrium clusters emission is studied. It is found that for the positive symmetry potential at supra-saturation densities the neutron-to-proton ratio of lighter clusters with mass number A less than or similar to 3 [(n/p)(A less than or similar to 3)] is larger than that of the heavier clusters with mass number A > 3 [(n/p)(A>3)], whereas for the negative symmetry potential at supra-saturation densities the (n/p)(A less than or similar to 3) is smaller than the (n/p)(A>3). This may be considered as a probe of the negative symmetry potential at supra-saturation densities.