Investigation on the oscillation modes of neutron stars Department of Physics, South China Uni v. of Tech. 华华华华华华华华华华 () collaborators Bao-An Li, William Newton, Plamen Krastev 文文文 epartment of Physics and astronomy, Texas A&M University-Comm 华华华华华华华华华华华华华华华华华华华华华华华华华 华华 2012. 4. 13
33
Embed
Investigation on the oscillation modes of neutron stars Department of Physics, South China Univ. of Tech. (华南理工大学物理系) collaborators Bao-An Li, William.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Investigation on the oscillation modes of
neutron stars
Department of Physics, South China Univ. of Tech. (华南理工大学物理系)
collaborators
Bao-An Li, William Newton, Plamen Krastev
文德华
Department of Physics and astronomy, Texas A&M University-Commerce
第十四届全国核结构大会暨第十次全国核结构专题讨论会
湖州 2012. 4. 13
Outline
I. W-modes in neutron stars
II. R-modes in neutron stars
Axial mode: under the angular transformation θ→ π − θ, ϕ → π + ϕ, a spherical harmonic function with index ℓ transforms as (−1)ℓ+1 for the expanding metric functions.
The non-radial neutron star oscillations could be triggered by various mechanisms such as gravitational collapse, a pulsar “glitch” or a phase transition of matter in the inner core.
Polar mode: transforms as (−1)ℓOscillating neutron star
• Introduction of axial w-mode
I. W-modes in neutron star
Axial w-mode: not accompanied by any matter motions and only the perturbation of the space-time, exists for all relativistic stars, including neutron star and black holes.
One major characteristic of the axial w-mode is its high frequency accompanied by very rapid damping.
(1) The w-modes are very important for astrophysical applications. The gravitational wave frequency of the axial w-mode depends on the neutron star’s structure and properties, which are determined by the EOS of neutron-rich stellar matter.
(2) It is helpful to the detection of gravitational waves to investigate the imprint of the nuclear symmetry energy constrained by very recent terrestrial nuclear laboratory data on the gravitational waves from the axial w-mode.
• Motivation
Key equation of axial w-mode
Inner the star (l=2)
Outer the star
The equation for oscillation of the axial w-mode is give by1
dr
de
dr
d *
drerr
0*
where
or
]6)(6[ 33
2
mprrr
eV
][6
3
2
Mrr
eV
0)]([ 22
*
2
zrVdr
zd
ii 0
1 S.Chandrasekhar and V. Ferrari, Proc. R. Soc. London A, 432, 247(1991) Nobel prize in 1983
It was shown that only values of x in the range between −1 (MDIx-1) and 0 (MDIx0) are consistent with the isospin-diffusion and isoscaling data at sub-saturation densities.
Here we assume that the EOS can be extrapolated to supra-saturation densities according to the MDI predictions.
It was shown that only values of x in the range between −1 (MDIx-1) and 0 (MDIx0) are consistent with the isospin-diffusion and isoscaling data at sub-saturation densities.
Here we assume that the EOS can be extrapolated to supra-saturation densities according to the MDI predictions.1. L.W.Chen, C. M. Ko, and B. A. Li, Phys. Rev.
Lett. 94, 032701 (2005).
2.B. A. Li, L.W. Chen, and C.M. Ko,
Phys. Rep. 464, 113 (2008).
•EOS constrained by terrestrial laboratory data
M-R relation
Wen D. H., Li B.A. and Krastev P.G., Phys. Rev. C 80, 025801 (2009)
Wen D. H., Li B.A. and Krastev P.G., Phys. Rev. C 80, 025801 (2009)
Frequency damping time
•Numerical Result and Discussion
1
Scaling characteristic
Wen D. H., Li B.A. and Krastev P.G., Phys. Rev. C 80, 025801 (2009)
Based on this linear dependence of the scaled frequency, the wII -mode is found to exist about compactness M/R>0.1078.
Exists linear fit
Wen D. H., Li B.A. and Krastev P.G., Phys. Rev. C 80, 025801 (2009)
Conclusion 1. The density dependence of the nuclear symmetry energy affects significantly both the frequencies and the damping times of axial w-mode.
2. Obtain a better scaling characteristic through scalin
g the eigen-frequency by the gravitational energy.
3. Give a general limit, M/R~0.1078, based on the linea
r scaling characteristic of wII, below this limit, wII-mode
will disappear.
In Newtonian theory, the fundamental dynamical equation (Euler equations) that governs the fluid motion in the co-rotating frame is
AccelerationCoriolis force
centrifugal force
external force
where is the fluid velocity and represents the gravitational potential.u
Φ
•Euler equations in the rotating frame
(I) Background and Motivation
II. R-modes in neutron star
For the rotating stars, the Coriolis force provides a restoring force for the toroidal modes, which leads to the so-called r-modes. Its eigen-frequency is
]1[)1(
2 23
2 M
R
ll
mr
It is shown that the structure parameters (M and R) make sense for the through the second order of .r
•Definition of r-mode
Class. Quantum Grav. 20 (2003) R105P111/p113
)1(
2
ll
mror
•CFS instability and canonical energy APJ,222(1978)281
The function Ec govern the stability to nonaxisymmetric perturbations
as: (1) if , stable; (2) if , unstable.
For the r-mode, The condition Ec < 0 is equivalent to a change of sign in the pattern speed as viewed in the inertial frame, which is always satisfied for r-mode.
gr-qc/0010102v1
canonical energy (conserved in absence of radiation and viscosity):
0)( cE 0)( cE
)1(
2
llr
)1(
)2)(1(2
ll
llri
Seen by a non-rotating observer(star is rotating faster than the r-mode pattern speed)
seen by a co-rotating observer. Looks like it's moving backwards
• The fluid motion has no radial component, and is the same inside the star although smaller by a factor of the square of the distance from the center.
• Fluid elements (red buoys) move in ellipses around their unperturbed locations.
Note: The CFS instability is not only existed in GR, but also existed in Newtonian theory.
Images of the motion of r-modes
•Viscous damping instability
• The r-modes ought to grow fast enough that they are not
completely damped out by viscosity.
•Two kinds of viscosity, bulk and shear viscosity, are normally
considered.
•At low temperatures (below a few times 109 K) the main
viscous dissipation mechanism is the shear viscosity arises
from momentum transport due to particle scattering..
•At high temperature (above a few times 109 K) bulk viscosity
is the dominant dissipation mechanism. Bulk viscosity arises
because the pressure and density variations associated with the
mode oscillation drive the fluid away from beta equilibrium.
•The r-mode instability window
Condition: To have an instability we need tgw to be smaller than both tsv and tbv.
For l = m = 2 r-mode of a canonical neutron star (R = 10 km and M = 1.4M⊙ and Kepler period PK ≈ 0.8 ms (n=1 polytrope)).
Int.J.Mod.Phys. D10 (2001) 381
• Motivations
(a) Old neutron stars (having crust) in LMXBs with rapid rotating fr
equency (such as EXO 0748-676) may have high core temperature (ar
Xiv:1107.5064v1.); which hints that there may exist r-mode instability
in the core.
(b) The discovery of massive neutron star (PRS J1614-2230, Nature 467, 10
81(2010) and EXO 0748-676, Nature 441, 1115(2006)) reminds us restudy the r-mode instability of massive NS, as most of the previous work focused on the 1.4Msun neutron star.
(c) The constraint on the symmetric energy at sub-saturation density
range and the core-crust transition density by the terrestrial nucl
ear laboratory data could provide constraints on the r-mode inst
ability.
PhysRevD.62.084030
Here only considers l=2, I2=0.80411. And the viscosity c is density and temperature dependent:
The subscript c denotes the quantities at the outer edge of the core.
T<109 K:
T>109 K:
The viscous timescale for dissipation in the boundary layer:
(II). Basic equations for r-mode instability window of neutron star with rigid crust
The gravitational radiation timescale:
According to , the critical rotation frequency is obtained:
Based on the Kepler frequency, the critical temperature defined as:
PhysRevD.62.084030
(III). Numerical Results
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
The EOSs are calculated using a model for the energy density of nuclear matter and probe the dependence on the symmetry energy by varying the slope of the symmetry energy at saturation density L from 25 MeV (soft) to 105 MeV (stiff). The crust-core transition density, and thus crustal thickness, is calculated consistently with the core EOS.
Equation of states
W. G. Newton, M. Gearheart, and B.-A. Li, arXiv:1110.4043v1.
The mass-radius relation and the core radius
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
The viscous timescale
Comparing the time scale
The gravitational radiation timescale
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
The lower boundary of the r-mode instability window for a 1.4Msun (a) and a 2.0Msun (b) neutron star over the range of the slope of the symmetry energy L consistent with experiment.
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
The location of the observed short-recurrence-time LMXBs in frequency-temperature space, for a 1.4Msun (a) and a 2.0Msun (b) neutron star.
The temperatures are derived from their observed accretion luminosity and assuming the cooling is dominant by the modified Urca neutrino emission process for normal nucleons or by the modified Urca neutrino emission process for neutrons being super-fluid and protons being super-conduction. Phys. Rev. Lett. 107, 101101(2011)
The critical temperature Tc for the onset of the CFS instability vs the crust-core transition densities over the range of the slope of the symmetry energy L consistent with experiment for 1.4Msun and 2.0Msun stars.
D.H. Wen, W. G. Newton, and B.A. Li , Phys. Rev. C 85, 025801 (2012)
Conclusion
(1)Smaller values of L help stabilize neutron stars against
runaway r-mode oscillations;
(2) A massive neutron star has a wider instability window;
(3)Treating consistently the crust thickness and core EOS,
and concluding that a thicker crust corresponds to a
lower critical temperature.
Thanks!
The standard axial w-mode is categorized as wI . The
high order axial w-modes are marked as the second w-
mode (wI2 -mode), the third mode (wI3 -mode) and so on.
An interesting additionally family of axial w-modes is ca
tegorized as wII .
Constrain by the flow data of relativistic heavy-ion reactions
P. Danielewicz, R. Lacey and W.G. Lynch, Science 298 (2002) 1592
1.M.B. Tsang, et al, Phys. Rev. Lett.
92, 062701 (2004)
2. B. A. Li, L.W. Chen, and C.M. Ko,
Phys. Rep. 464, 113 (2008).
The gravitational energy is calculated from
1S.Weinberg, Gravitation and cosmology, (New York: Wiley,1972)