Modeling Neutron Star Oscillations in a Fixed General Relativistic Background Including Solid Crust Dynamics
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
κ Andromedae, an early-type star that hosts a directly imaged low-mass companion, is expected to be oblate due to its rapid rotational velocity (v sin i = ∼162 km s{sup −1}). We observed the star with the CHARA Array’s optical beam combiner, PAVO, measuring its size at multiple orientations and determining its oblateness. The interferometric measurements, combined with photometry and this v sin i value are used to constrain an oblate star model that yields the fundamental properties of the star and finds a rotation speed that is ∼85% of the critical rate and a low inclination of ∼30°. Three modeled properties (the average radius, bolometric luminosity, and equatorial velocity) are compared to MESA evolution models to determine an age and mass for the star. In doing so, we determine an age for the system of 47{sub −40}{sup +27} Myr. Based on this age and previous measurements of the companion’s temperature, the BHAC15 evolution models imply a mass for the companion of 22{sub −9}{sup +8} M {sub J}.
New models of star-forming filamentary clouds are presented in order to quantify their properties and to predict their evolution. These 2D axisymmetric models describe filaments that have no core, one low-mass core, and one cluster-forming core. They are based on Plummer-like cylinders and spheroids that are bounded by a constant-density surface of finite extent. In contrast to 1D Plummer-like models, they have specific values of length and mass, they approximate observed column density maps, and their distributions of column density ( N -pdfs) are pole-free. Each model can estimate the star-forming potential of a core-filament system by identifying the zone of gas dense enough to form low-mass stars and by counting the number of enclosed thermal Jeans masses. This analysis suggests that the Musca central filament may be near the start of its star-forming life, with enough dense gas to make its first ∼3 protostars, while the Coronet filament is near the midpoint of its star formation, with enough dense gas to add ∼8 protostars to its ∼20 known stars. In contrast, L43 appears to be near the end of its star-forming life, since it lacks enough dense gas to add any new protostars to the two young stellar objectsalready known.
The Qatar-2 transiting exoplanet system was recently observed in short-cadence mode by Kepler as part of K2 Campaign 6. We identify dozens of starspot-crossing events, when the planet eclipses a relatively dark region of the stellar photosphere. The observed patterns of these events demonstrate that the planet always transits over the same range of stellar latitudes and, therefore, that the stellar obliquity is less than about 10°. We support this conclusion with two different modeling approaches: one based on explicit identification and timing of the events and the other based on fitting the light curves with a spotted-star model. We refine the transit parameters and measure the stellar rotation period (18.5 ± 1.9 days), which corresponds to a “gyrochronological” age of 1.4 ± 0.3 Gyr. Coherent flux variations with the same period as the transits are well modeled as the combined effects of ellipsoidal light variations (15.4 ± 4.8 ppm) and Doppler boosting (14.6 ± 5.1 ppm). The magnitudes of these effects correspond to a planetary mass of 2.6±0.9 M{sub Jup} and 3.9±1.5 M{sub Jup}, respectively. Both of these independent mass estimates agree with the mass determined by the spectroscopic Doppler technique (2.487±0.086 M{sub Jup}). No occultations are detected, giving a 2σ upper limit of 0.06 on the planet’s visual geometric albedo. We find no evidence for orbital decay, although we are only able to place a weak lower bound on the relevant tidal quality factor: Q{sub ⋆}{sup ′}>1.5×10{sup 4} (95% confidence).
The hypotheses underlying theoretical studies of the evolution of massive model stars with and without mass loss are summarized. The evolutionary tracks followed by the models across theoretical Hertzsprung-Russell (HR) diagrams are compared with the observed distribution of B stars in an HR diagram. The pulsational properties of models of massive star are also described.
We present contemporaneous, broadband, near-infrared spectroscopy (0.9-2.45 {mu}m) and H-band photometry of the black hole X-ray binary, XTE J1118+480. We determined the fractional dilution of the NIR ellipsoidal light curves of the donor star from other emission sources in the system by comparing the absorption features in the spectrum with field stars of known spectral type. We constrained the donor star spectral type to K7 V-M1 V and determined that the donor star contributed 54% {+-} 27% of the H-band flux at the epoch of our observations. This result underscores the conclusion that the donor star cannot be assumed to be the only NIR emission source in quiescent X-ray binaries. The H-band light curve shows a double-humped asymmetric modulation with extra flux at orbital phase 0.75. The light curve was fitted with a donor star model light curve, taking into account a constant second flux component based on the dilution analysis. We also fitted models that included emission from the donor star, a constant component from the accretion disk, and a phase-variable component from the bright spot where the mass accretion stream impacts the disk. These simple models with reasonable estimates for the component physical parameters can fully account for the observed light curve, including the extra emission at phase 0.75. From our fits, we constrained the binary inclination to 68 Degree-Sign {<=} i {<=} 79 Degree-Sign . This leads to a black hole mass of 6.9 M{sub Sun} {<=} M{sub BH} {<=} 8.2 M{sub Sun }. Long-term variations in the NIR light curve shape in XTE J1118+480 are similar to those seen in other X-ray binaries and demonstrate the presence of continued activity and variability in these systems even when in full quiescence.
Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.
We report C, Si, and S isotope measurements on 34 presolar silicon carbide grains of Type AB, characterized by {sup 12}C/{sup 13}C < 10. Nitrogen, Mg-Al-, and Ca-Ti-isotopic compositions were measured on a subset of these grains. Three grains show large {sup 32}S excesses, a signature that has been previously observed for grains from supernovae (SNe). Enrichments in {sup 32}S may be due to contributions from the Si/S zone and the result of S molecule chemistry in still unmixed SN ejecta or due to incorporation of radioactive {sup 32}Si from C-rich explosive He shell ejecta. However, a SN origin remains unlikely for the three AB grains considered here, because of missing evidence for {sup 44}Ti, relatively low {sup 26}Al/{sup 27}Al ratios (a few times 10{sup –3}), and radiogenic {sup 32}S along with low {sup 12}C/{sup 13}C ratios. Instead, we show that born-again asymptotic giant branch (AGB) stars that have undergone a very-late thermal pulse (VLTP), known to have low {sup 12}C/{sup 13}C ratios and enhanced abundances of the light s-process elements, can produce {sup 32}Si, which makes such stars attractive sources for AB grains with {sup 32}S excesses. This lends support to the proposal that at least some AB grains originate from born-again AGB stars, although uncertainties in the born-again AGB star models and possible variations of initial S-isotopic compositions in the parent stars of AB grains make it difficult to draw a definitive conclusion.
Presolar SiC grains fall into several groups based on C, N, and Si isotopic compositions. Approximately 93% are defined as mainstream, having 10 less than C-12/C-13 less than 100 and N-14/N-15 ranging from 50 to 20,000. A number of studies have shown that the most likely sources of mainstream grains are low mass asymptotic giant branch stars. Models of nucleosynthesis in AGB stars reproduce the s-process enhancements seen in the heavy elements in mainstream SiC grains. Among the less common grains, A+B grains, which comprise approximately 3-4% of presolar SiC, are perhaps the least well understood. Recent studies by Amari et al. show that A+B grains can be divided into at least 4 groups based on their trace element concentration patterns. Of 20 grains studied, 7 showed trace element patterns consistent with condensation from a gas of solar system composition, while the rest had varying degrees of process enhancements. Our previous measurements on 3 A+B grains showed Mo of solar isotopic composition, but Zr with a strong enhancement in 96Zr, which is an r-process isotope but can be made in an sprocess if the neutron density is high enough to bridge the unstable Zr-95 (T(sub 1/2)= 64 d). The observation of Mo with solar system isotopic composition in the same grains is puzzling however. Meyer et al. have recently shown that a neutron burst mechanism can produce a high Zr-96/Zr-94 without enhancing Mo-100, however this model leads to enhancements in Mo-95 and Mo-97 not observed in A+B grains. We report here results of Mo measurements on 7 additional A+B grains, and Ba measurements on 2 A+B grains, and compare these to the previous studies.
Relations between neutron star properties that do not depend on the nuclear equation of state offer insights on neutron star physics and have practical applications in data analysis. Such relations are obtained by fitting to a range of phenomenological or nuclear physics equation of state models, each of which may have varying degrees of accuracy. In this study we revisit commonly used relations and reassess them with a very flexible set of phenomenological nonparametric equation of state models that are based on Gaussian processes. Our models correspond to two sets: equations of state which mimic hadronic models, and equations of state with rapidly changing behavior that resemble phase transitions. Here we quantify the accuracy of relations under both sets and discuss their applicability with respect to expected upcoming statistical uncertainties of astrophysical observations. We further propose a goodness-of-fit metric which provides an estimate for the systematic error introduced by using the relation to model a certain equation-of-state set. Overall, the nonparametric distribution is more poorly fit with existing relations, with the I–Love–Q relations retaining the highest degree of universality. Fits degrade for relations involving the tidal deformability, such as the binary-Love and compactness-Love relations, and when introducing phase transition phenomenology. For most relations, systematic errors are comparable to current statistical uncertainties under the nonparametric equation of state distributions.
The internal rotation of post-main sequence stars is investigated, in response to the convective pumping of angular momentum toward the stellar core, combined with a tight magnetic coupling between core and envelope. The spin evolution is calculated using model stars of initial mass 1, 1.5, and 5 M{sub ⊙}, taking into account mass loss on the giant branches. We also include the deposition of orbital angular momentum from a sub-stellar companion, as influenced by tidal drag along with the excitation of orbital eccentricity by a fluctuating gravitational quadrupole moment. A range of angular velocity profiles Ω(r) is considered in the envelope, extending from solid rotation to constant specific angular momentum. We focus on the backreaction of the Coriolis force, and the threshold for dynamo action in the inner envelope. Quantitative agreement with measurements of core rotation in subgiants and post-He core flash stars by Kepler is obtained with a two-layer angular velocity profile: uniform specific angular momentum where the Coriolis parameter Co≡Ωτ{sub con}≲1 (here τ{sub con} is the convective time), and Ω(r)∝r{sup −1} where Co≳1. The inner profile is interpreted in terms of a balance between the Coriolis force and angular pressure gradients driven by radially extended convective plumes. Inward angular momentum pumping reduces the surface rotation of subgiants, and the need for a rejuvenated magnetic wind torque. The co-evolution of internal magnetic fields and rotation is considered in Kissin and Thompson, along with the breaking of the rotational coupling between core and envelope due to heavy mass loss.
Isospin-violating dark matter (IVDM), which couples differently with protons and neutrons, provides a promising mechanism to ameliorate the tension among recent direct detection experiments. Assuming dark matter (DM) is non-interacting bosonic asymmetric IVDM, we investigate how the existence of old neutron stars limits the DM–proton scattering cross-section σ{sub p}, especially the effects of the isospin-violating DM–nucleon interactions and the symmetry energy in the equation of state (EOS) of isospin asymmetric nuclear matter. Our calculations are completely based on general relativity and the structure of neutron stars is obtained by solving the Tolman–Oppenheimer–Volkoff equations with nuclear matter EOS constrained by terrestrial experiments. We find that, by considering the more realistic neutron star model rather than a simple uniform neutron sphere as usual, the σ{sub p} bounds from old neutron stars can be varied by more than an order of magnitude depending on the specific values of the DM neutron-to-proton coupling ratio f{sub n}/f{sub p}, and they can be further varied by more than a factor of two depending on the density dependence of the symmetry energy. In particular, we demonstrate that the observed nearby isolated old neutron star PSR B1257+12 can set a very strong limit on σ{sub p} for low-mass DM particles (⩽20 GeV) that reaches a sensitivity beyond the current best limits from direct detection experiments and disfavors the DM interpretation of previously reported positive experimental results, including the IVDM.
We study the mass-transfer rates and disk stability conditions of ultracompact X-ray binaries (UCXBs) using empirical time-averaged X-ray luminosities from Paper I and compiled information from the literature. The majority of UCXBs are consistent with evolutionary tracks for white dwarf donors. Three UCXBs with orbital periods longer than 40 minutes have mass-transfer rates above 10{sup -10} M{sub Sun} yr{sup -1}, inconsistent with white dwarf donor tracks. We show that if helium star donors can retain their initial high entropy, they can explain the observed mass-transfer rates of these UCXBs. Several UCXBs show persistent luminosities apparently below the disk instability limit for irradiated He accretion disks. We point out that a predominantly C and/or O disk (as observed in the optical spectra of several) lowers the disk instability limit, explaining this disagreement. The orbital period and low time-averaged mass-transfer rate of 2S 0918-549 provide evidence that the donor star is a low-entropy C/O white dwarf, consistent with optical spectra. We combine existing information to constrain the masses of the donors in 4U 1916-053 (0.064 {+-} 0.010 M{sub Sun }) and 4U 1626-67 (<0.036 M{sub Sun} for a 1.4 M{sub Sun} neutron star). We show that 4U 1626-67 is indeed persistent, and not undergoing a transient outburst, leaving He star models as the best explanation for the donor.
Analytical model of eclipsing binary star systems, testing validity through numerical integration error analysis
By using updated stellar low-mass stars models, we systematically investigate the nucleosynthesis processes occurring in asymptotic giant branch (AGB) stars. In this paper, we present a database dedicated to the nucleosynthesis of AGB stars: FRANEC Repository of Updated Isotopic Tables and Yields (FRUITY). An interactive Web-based interface allows users to freely download the full (from H to Bi) isotopic composition, as it changes after each third dredge-up (TDU) episode and the stellar yields the models produce. A first set of AGB models, having masses in the range 1.5 {<=}M/M{sub Sun} {<=} 3.0 and metallicities 1 Multiplication-Sign 10{sup -3} {<=} Z {<=} 2 Multiplication-Sign 10{sup -2}, is discussed. For each model, a detailed description of the physical and the chemical evolution is provided. In particular, we illustrate the details of the s-process and we evaluate the theoretical uncertainties due to the parameterization adopted to model convection and mass loss. The resulting nucleosynthesis scenario is checked by comparing the theoretical [hs/ls] and [Pb/hs] ratios to those obtained from the available abundance analysis of s-enhanced stars. On the average, the variation with the metallicity of these spectroscopic indexes is well reproduced by theoretical models, although the predicted spread at a given metallicity is substantially smaller than the observed one. Possible explanations for such a difference are briefly discussed. An independent check of the TDU efficiency is provided by the C-stars luminosity function. Consequently, theoretical C-stars luminosity functions for the Galactic disk and the Magellanic Clouds have been derived. We generally find good agreement with observations.
Tables of evolutionary sequences of models for massive stars have been prepared for a variety of physical input parameters that are normally treated as free. These parameters include the interior convective mixing scheme, the mixing length in the outer convective envelope, the rate of stellar-wind mass loss, the initial stellar mass, and the initial chemical composition. Ranges of specified initial mass and initial chemical composition are M = 10-120 solar masses, Xe = 0.602-0.739, and Ze = 0.021-0.044. The tables cover evolution of the star from the ZAMS to either the end of core H burning or the end of core He burning. Differences among the evolutionary tracks are illustrated primarily in terms of the interior mixing scheme, since the amount and timing of stellar wind mass loss are still very uncertain for initial masses above about 30 solar masses.
We performed detailed chemical abundance analysis of the extremely metal-poor ([Ar/H] ∼ –2) halo planetary nebula (PN) H4-1 based on the multi-wavelength spectra from Subaru/HDS, GALEX, SDSS, and Spitzer/IRS and determined the abundances of 10 elements. The C and O abundances were derived from collisionally excited lines (CELs) and are almost consistent with abundances from recombination lines (RLs). We demonstrated that the large discrepancy in the C abundance between CEL and RL in H4-1 can be solved using the temperature fluctuation model. We reported the first detection of the [Xe III] λ5846 line in H4-1 and determination of its elemental abundance ([Xe/H] > +0.48). H4-1 is the most Xe-rich PN among the Xe-detected PNe. The observed abundances are close to the theoretical prediction by a 2.0 M {sub ☉} single star model with an initially element rich ([r/Fe] = +2.0 dex) rapid neutron-capture process (r-process). The observed Xe abundance would be a product of the r-process in primordial supernovae. The [C/O]-[Ba/(Eu or Xe)] diagram suggests that the progenitor of H4-1 shares the evolution with carbon-enhanced metal-poor (CEMP)-r/s and CEMP-no stars. The progenitor of H4-1 is presumably a binary formed in an r-process-rich environment.
The measurement of neutron star mass and radius is one of the most direct ways to distinguish between various dense matter equations of state. The mass and radius of accreting neutron stars hosted in low-mass X-ray binaries can be constrained by several methods, including photospheric radius expansion from type I X-ray bursts and from quiescent spectra. In this paper, we apply for the first time these two methods simultaneously to constrain the mass and radius of Aql X–1. The quiescent spectra from Chandra and XMM-Newton , and photospheric radius expansion bursts from RXTE are used. The determination of the mass and radius of Aql X–1 is also used to verify the consistency between the two methods and to narrow down the uncertainties of the neutron star mass and radius. It is found that the distance to Aql X–1 should be in the range of 4.0–5.75 kpc, based on the overlapping confidence regions between photospheric radius expansion burst and quiescent spectra methods. In addition, we show that the mass and radius determined for the compact star in Aql X–1 are compatible with strange star equations of state and conventional neutron star models.