Search NASA⌕ Search

SEARCH · Search NASA

Results for “STAR MODELS”

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.

At least 91 records · Page 5

NS 1987A in SN 1987A

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object’s predicted position. The infrared excess could be due to the decay of isotopes like 44 Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spin down of a pulsar, or to thermal emission from an embedded, cooling neutron star. In this work it is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications that the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 yr old neutron star. Furthermore, models of cooling neutron stars within the minimal cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n- 1 S 0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations do not overlap, then pulsar spin down or accretion might be more likely, but the pulsar’s period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling UV Radiation Feedback from Massive Stars. I. Implementation of Adaptive Ray-tracing Method and Tests [Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Algorithm]

Here, we present an implementation of an adaptive ray-tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses nonblocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems, including the propagation of radiation in vacuum and the expansions of various types of H II regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim {10}^{3}$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is 12% when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with those from the M 1 closure relation. Although the ART and M 1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Galactic Neutron Star Population for Use in Continuous Gravitational-wave Searches

Searches for continuous gravitational waves (GW) from unknown Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce GW if asymmetric deformations exist in its structure that are characterized by the star's ellipticity. Here, in this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to estimate the total number of neutron stars probed, using GW, to a given upper limit on the ellipticity. This may help optimize future searches with improved sensitivity. The improved sensitivity of third-generation GW detectors may increase the number of neutron stars probed, to a given ellipticity, by factors of 100 to 1000.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE STAR FORMATION LAW FOR SPIRAL AND IRREGULAR GALAXIES

A dynamical model for star formation on a galactic scale is proposed in which the interstellar medium is constantly condensing to star-forming clouds on the dynamical time of the average midplane density, and the clouds are constantly being disrupted on the dynamical timescale appropriate for their higher density. In this model, the areal star formation rate scales with the 1.5 power of the total gas column density throughout the main regions of spiral galaxies, and with a steeper power, 2, in the far outer regions and in dwarf irregular galaxies because of the flaring disks. At the same time, there is a molecular star formation law that is linear in the main and outer parts of disks and in dIrrs because the duration of individual structures in the molecular phase is also the dynamical timescale, canceling the additional 0.5 power of surface density. The total gas consumption time scales directly with the midplane dynamical time, quenching star formation in the inner regions if there is no accretion, and sustaining star formation for ∼100 Gyr or more in the outer regions with no qualitative change in gas stability or molecular cloud properties. The ULIRG track follows from high densities in galaxy collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

Optical and Near-infrared Radial Velocity Content of M Dwarfs: Testing Models with Barnard’s Star

High-precision radial velocity (RV) measurements have been central in the study of exoplanets during the last two decades, from the early discovery of hot Jupiters, to the recent mass measurements of Earth-sized planets uncovered by transit surveys. While optical RV is now a mature field, there is currently a strong effort to push the technique into the near-infrared domain (chiefly Y, J, H, and K bandpasses) to probe planetary systems around late-type stars. The combined lower mass and luminosity of M dwarfs leads to an increased reflex RV signal for planets in the habitable zone compared to Sun-like stars. The estimates on the detectability of planets rely on various instrumental characteristics but also on a prior knowledge of the stellar spectrum. While the overall properties of M dwarf spectra have been extensively tested against observations, the same is not true for their detailed line profiles, which leads to significant uncertainties when converting a given signal-to-noise ratio to a corresponding RV precision as attainable on a given spectrograph. By combining archival CRIRES and HARPS data with ESPaDOnS data of Barnard’s star, we show that state-of-the-art atmosphere models over-predict the Y- and J-band RV content by more than a factor of ∼2, while under-predicting the H- and K-band content by half.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Models of R Coronae Borealis Stars

We present an improved numerical method to model subsolar He+CO-WD merger progenitors of R Corona Borealis stars that builds on our previous work. These improvements include a smooth entropy transition from the core to the envelope of the post-merger, inclusion of single-zone nucleosynthesis to mimic the effects of burning during the merger event, and post-processing the models with a larger nuclear network for analysis of s-process nucleosynthesis. We perform a parameter study to understand the effects of the entropy transition, peak temperature, and overshooting on our models. The models that best agree with observations of R Corona Borealis stars are processed with a much larger nuclear network to investigate s-process nucleosynthesis and the dredge-up of s-process products into the outer envelope in detail. We present a model with a significant enhancement in s-process elements, which also agrees with observed surface abundances and isotopic ratios of 16 O/ 16 O and C/O between 1 and 10. Finally, we find that the neutron exposure and initial neutron densities this model requires to obtain such an enhancement are much more consistent with i-process nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetized Winds of M-type Stars and Star–Planet Magnetic Interactions: Uncertainties and Modeling Strategy

M-type stars are the most common stars in the Universe. They are ideal hosts for the search of exoplanets in the habitable zone (HZ), as their small size and low temperature make the HZ much closer-in than their solar twins. Harboring very deep convective layers, they also usually exhibit very intense magnetic fields. Understanding their environment, in particular their coronal and wind properties, is thus very important, as they might be very different from what is observed in the solar system. The mass-loss rate of M-type stars is poorly known observationally, and recent attempts to estimate it for some of them (e.g., TRAPPIST-1 and Proxima Centauri) can vary by an order of magnitude. In this work, we revisit the stellar wind properties of M dwarfs in the light of the latest estimates of $\dot{M}$ through Lyα absorption at the astropause and slingshot prominences. We outline a modeling strategy to estimate the mass-loss rate, radiative loss, and wind speed, with uncertainties, based on an Alfvén-wave-driven stellar wind model. We find that it is very likely that several TRAPPIST-1 planets lie within the Alfvén surface, which implies that these planets experience star–planet magnetic interactions (SPMIs). We also find that SPMIs between Proxima Cen b and its host star could be the reason for recently observed radio emissions.

M stars↗

GW170817 implications on the frequency and damping time of $f$-mode oscillations of neutron stars

Within a minimum model for neutron stars consisting of nucleons, electrons, and muons at $β$ equilibrium using about a dozen equations of state (EOS) from microscopic nuclear many-body theories and 40000 EOSs randomly generated using an explicitly isospin-dependent parametric EOS model for high-density neutron-rich nucleonic matter within its currently known uncertainty range, we study correlations among the $f$-mode frequency, its damping time and the tidal deformability as well as the compactness of neutron stars. Except for quark stars, both the $f$-mode frequency and damping time of canonical neutron stars are found to scale with the tidal deformability independent of the EOSs used. Applying the constraint on the tidal deformability of canonical neutron stars $Λ_{1.4} = 190^{+390}_{-120}$ extracted by the LIGO+VIRGO Collaborations from their improved analyses of the GW170817 event, the $f$-mode frequency and its damping time of canonical neutron stars are limited to 1.67–2.18 kHz and 0.155–0.255 s, respectively, providing a useful guidance for the ongoing search for gravitational waves from the $f$-mode oscillations of isolated neutron stars. Finally, assuming either or both the $f$-mode frequency and its damping time will be measured precisely in future observations with advanced gravitational wave detectors, we discuss how information about the mass and/or radius as well as the still rather elusive nuclear symmetry energies at suprasaturation densities may be extracted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Gaia DR2 parallax zero-point: hierarchical modelling of red clump stars

ABSTRACT The systematic offset of Gaia parallaxes has been widely reported with Gaia’s second data release, and it is expected to persist in future Gaia data. In order to use Gaia parallaxes to infer distances to high precision, we develop a hierarchical probabilistic model to determine the Gaia parallax zero-point offset along with the calibration of an empirical model for luminosity of red clump stars by combining astrometric and photometric measurements. Using a cross-matched sample of red clump stars from the Apache Point Observatory Galactic Evolution Experiment and Gaia Data Release 2 (DR2), we report the parallax zero-point offset in DR2 to be $\varpi _0 = -48 \pm 1\,\mu \text{as}$. We infer the red clump absolute magnitude to be −1.622 ± 0.004 in Ks, 0.435 ± 0.004 in G, −1.019 ± 0.004 in J, and −1.516 ± 0.004 in H. The intrinsic scatter of the red clump is ∼0.09 mag in J, H, and Ks, or ∼ 0.12 mag in G. We tailor our models to accommodate more complex analyses such as investigating the variations of the parallax zero-point with each source’s observed magnitude, observed colour, and sky position. In particular, we find fluctuations of the zero-point across the sky to be of order or less than a few 10s of $\mu \text{as}$.

Chan, Victor C.↗

STAR-CCM+ (CFD) Calculations and Validation L3:VVI.H2L.P15.02

This milestone presents a demonstration of the High-to-Low (Hi2Lo) process in the VVI focus area. Validation and additional calculations with the commercial computational fluid dynamics code, STAR-CCM+, were performed using a 5x5 fuel assembly with non-mixing geometry and spacer grids. This geometry was based on the benchmark experiment provided by Westinghouse. Results from the simulations were compared to existing experimental data and to the subchannel thermal-hydraulics code COBRA-TF (CTF). An uncertainty quantification (UQ) process was developed for the STAR-CCM+ model and results of the STAR UQ were communicated to CTF. Results from STAR-CCM+ simulations were used as experimental design points in CTF to calibrate the mixing parameter β and compared to results obtained using experimental data points. This demonstrated that CTF’s β parameter can be calibrated to match existing experimental data more closely. The Hi2Lo process for the STAR-CCM+/CTF code coupling was documented in this milestone and closely linked L3:VVI.H2LP15.01 milestone report.

97 MATHEMATICS AND COMPUTING↗

EG ANDROMEDAE: A NEW ORBIT AND ADDITIONAL EVIDENCE FOR A PHOTOIONIZED WIND

We analyze a roughly 20 yr set of spectroscopic observations for the symbiotic binary EG And. Radial velocities derived from echelle spectra are best fit with a circular orbit having an orbital period of P = 483.3 ± 1.6 days and semi-amplitude K = 7.34 ± 0.07 km s{sup −1}. Combined with previous data, these observations rule out an elliptical orbit at the 10 σ level. Equivalent widths of H i Balmer emission lines and various absorption features vary in phase with the orbital period. Relative to the radius of the red giant primary, the apparent size of the H ii region is consistent with a model where a hot secondary star with effective temperature T{sub h} ≈ 75,000 K ionizes the wind from the red giant.

79 ASTRONOMY AND ASTROPHYSICS↗

AN APPARENT PRECESSING HELICAL OUTFLOW FROM A MASSIVE EVOLVED STAR: EVIDENCE FOR BINARY INTERACTION

Massive, evolved stars play a crucial role in the metal enrichment, dust budget, and energetics of the interstellar medium; however, the details of their evolution are uncertain because of their rarity and short lifetimes before exploding as supernovae. Discrepancies between theoretical predictions from single-star evolutionary models and observations of massive stars have evoked a shifting paradigm that implicates the importance of binary interaction. We present mid- to far-infrared observations from the Stratospheric Observatory for Infrared Astronomy of a conical “helix” of warm dust (∼180 K) that appears to extend from the Wolf–Rayet star WR102c. Our interpretation of the helix is a precessing, collimated outflow that emerged from WR102c during a previous evolutionary phase as a rapidly rotating luminous blue variable. We attribute the precession of WR102c to gravitational interactions with an unseen compact binary companion whose orbital period can be constrained to 800 days < P < 1400 days from the inferred precession period, τ{sub p} ∼ 1.4 × 10{sup 4} yr, and limits imposed on the stellar and orbital parameters of the system. Our results concur with the range of orbital periods (P ≲ 1500 days) where spin-up via mass exchange is expected to occur for massive binary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Convection and Overshoot in Models of Doradus and Scuti Stars

We investigate the pulsation properties of stellar models that are representative of δ Scuti and γ Doradus variables. Here we have calculated a grid of stellar models from 1.2 to 2.2 M ⊙, including the effects of both rotation and convective overshoot using MESA, and we investigate the pulsation properties of these models using GYRE. We discuss the observable patterns in the frequency spacing for p modes and the period spacings for g modes. Using the observable patterns in the g mode period spacings, it may be possible to observationally constrain the convective overshoot and rotation of a model. We also calculate the pulsation constant (Q) for all models in our grid and investigate the variation with convective overshoot and rotation. The variation in the Q values of the radial modes can be used to place constraints on the convective overshoot and rotation of stars in this region. Finally, as a test case, we apply this method to a sample of 22 High-Amplitude δ Scuti stars (HADS) and provide estimates for the convective overshoot of the sample.

79 ASTRONOMY AND ASTROPHYSICS↗

What Sets the Massive Star Formation Rates and Efficiencies of Giant Molecular Clouds?

Galactic star formation scaling relations show increased scatter from kpc to sub-kpc scales. Investigating this scatter may hold important clues to how the star formation process evolves in time and space. Here, we combine different molecular gas tracers, different star formation indicators probing distinct populations of massive stars, and knowledge of the evolutionary state of each star-forming region to derive the star formation properties of ∼150 star-forming complexes over the face of the Large Magellanic Cloud (LMC). We find that the rate of massive star formation ramps up when stellar clusters emerge and boost the formation of subsequent generations of massive stars. In addition, we reveal that the star formation efficiency of individual giant molecular clouds (GMCs) declines with increasing cloud gas mass (M{sub cloud}). This trend persists in Galactic star-forming regions and implies higher molecular gas depletion times for larger GMCs. We compare the star formation efficiency per freefall time (ϵ{sub ff}) with predictions from various widely used analytical star formation models. While these models can produce large dispersions in ϵ{sub ff} similar to those in observations, the origin of the model-predicted scatter is inconsistent with observations. Moreover, all models fail to reproduce the observed decline of ϵ{sub ff} with increasing M{sub cloud} in the LMC and the Milky Way. We conclude that analytical star formation models idealizing global turbulence levels and cloud densities and assuming a stationary star formation rate (SFR) are inconsistent with observations from modern data sets tracing massive star formation on individual cloud scales. Instead, we reiterate the importance of local stellar feedback in shaping the properties of GMCs and setting their massive SFR.

79 ASTRONOMY AND ASTROPHYSICS↗

DETECTION OF SOLAR-LIKE OSCILLATIONS, OBSERVATIONAL CONSTRAINTS, AND STELLAR MODELS FOR θ CYG, THE BRIGHTEST STAR OBSERVED BY THE KEPLER MISSION

θ Cygni is an F3 spectral type magnitude V = 4.48 main-sequence star that was the brightest star observed by the original Kepler spacecraft mission. Short-cadence (58.8 s) photometric data using a custom aperture were first obtained during Quarter 6 (2010 June–September) and subsequently in Quarters 8 and 12–17. We present analyses of solar-like oscillations based on Q6 and Q8 data, identifying angular degree l = 0, 1, and 2 modes with frequencies of 1000–2700 μ Hz, a large frequency separation of 83.9 ± 0.4 μ Hz, and maximum oscillation amplitude at frequency ν {sub max} = 1829 ± 54 μ Hz. We also present analyses of new ground-based spectroscopic observations, which, combined with interferometric angular diameter measurements, give T {sub eff} = 6697 ± 78 K, radius 1.49 ± 0.03 R {sub ⊙}, [Fe/H] = -0.02 ± 0.06 dex, and log g = 4.23 ± 0.03. We calculate stellar models matching these constraints using the Yale Rotating Evolution Code and the Asteroseismic Modeling Portal. The best-fit models have masses of 1.35–1.39 M {sub ⊙} and ages of 1.0–1.6 Gyr. θ Cyg’s T {sub eff} and log g place it cooler than the red edge of the γ Doradus instability region established from pre- Kepler ground-based observations, but just at the red edge derived from pulsation modeling. The pulsation models show γ Dor gravity modes driven by the convective blocking mechanism, with frequencies of 1–3 cycles per day (11 to 33 μ Hz). However, gravity modes were not seen in Kepler data; one signal at 1.776 cycles per day (20.56 μ Hz) may be attributable to a faint, possibly background, binary.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dense-gas Mass versus Star Formation Rate Relation: A Misleading Linearity?

We model the star formation relation of molecular clumps in dependence of their dense-gas mass when their volume density profile is that of an isothermal sphere (i.e., ρ{sub clump}(r)∝r{sup −2}). Dense gas is defined as gas whose volume density is higher than a threshold ρ{sub th}=700 M{sub ⊙} pc{sup −3} (i.e., HCN(1-0)-mapped gas). We divide the clump into two regions: a dense inner region (where ρ{sub clump}(r)⩾ρ{sub th}), and low-density outskirts (where ρ{sub clump}(r)<ρ{sub th}). We find that the total star formation rate of clumps scales linearly with the mass of their dense inner region, even when more than half of the clump star formation activity takes place in the low-density outskirts. We therefore emphasize that a linear star formation relation does not necessarily imply that star formation takes place exclusively in the gas whose mass is given by the star formation relation. The linearity of the star formation relation is strengthened when we account for the mass of dense fragments (e.g., cores, fibers) seeding star formation in the low-density outskirts, and which our adopted clump density profile ρ{sub clump}(r) does not resolve. We also find that the star formation relation is significantly tighter when considering the dense gas than when considering all the clump gas, as observed for molecular clouds of the Galactic plane. When the clumps have no low-density outskirts (i.e., they consist of dense gas only), the star formation relation becomes superlinear and progressively wider.

79 ASTRONOMY AND ASTROPHYSICS↗

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗

Chiral Effective Field Theory and the High-Density Nuclear Equation of State

Born in the aftermath of core-collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yield novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons among nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗