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 1,513 records · Page 84

STAR CLUSTER PROPERTIES IN TWO LEGUS GALAXIES COMPUTED WITH STOCHASTIC STELLAR POPULATION SYNTHESIS MODELS

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete initial mass function (IMF) sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

79 ASTRONOMY AND ASTROPHYSICS↗

Simultaneous calibration of spectro-photometric distances and the Gaia DR2 parallax zero-point offset with deep learning

ABSTRACT Gaia measures the five astrometric parameters for stars in the Milky Way, but only four of them (positions and proper motion, but not distance) are well measured beyond a few kpc from the Sun. Modern spectroscopic surveys such as APOGEE cover a large area of the Milky Way disc and we can use the relation between spectra and luminosity to determine distances to stars beyond Gaia’s parallax reach. Here, we design a deep neural network trained on stars in common between Gaia and APOGEE that determines spectro-photometric distances to APOGEE stars, while including a flexible model to calibrate parallax zero-point biases in Gaia DR2. We determine the zero-point offset to be $-52.3 \pm 2.0\, \mu \mathrm{as}$ when modelling it as a global constant, but also train a multivariate zero-point offset model that depends on G, GBP − GRP colour, and Teff and that can be applied to all ≈58 million stars in Gaia DR2 within APOGEE’s colour–magnitude range and within APOGEE’s sky footprint. Our spectro-photometric distances are more precise than Gaia at distances ${\gtrsim} 2\, \mathrm{kpc}$ from the Sun. We release a catalogue of spectro-photometric distances for the entire APOGEE DR14 data set which covers Galactocentric radii $2\, \mathrm{kpc} \lesssim R \lesssim 19\, \mathrm{kpc}$; ${\approx} 150\, 000$ stars have ${\lt} 10{{\ \rm per\ cent}}$ uncertainty, making this a powerful sample to study the chemo-dynamical structure of the disc. We use this sample to map the mean [Fe/H] and 15 abundance ratios [X/Fe] from the Galactic Centre to the edge of the disc. Among many interesting trends, we find that the bulge and bar region at $R \lesssim 5\, \mathrm{kpc}$ clearly stands out in [Fe/H] and most abundance ratios.

Leung, Henry W.↗

Computer experiments on the structure and dynamics of spiral galaxies

The evolution of an initially balanced rotating disk of stars with an initial velocity dispersion given by Toomre's local criterion was investigated by means of a computer model for isolated disks of stars. It was found that the disk is unstable against very large scale modes. A stable axisymmetric disk with a velocity dispersion much larger than that given by Toomre's criterion was generated. The final mass distribution for the disk gives a high density central core and a disk population of stars that is closely approximated by an exponential variation. Various methods and rates of cooling the hot axisymmetric disks were investigated. It was found that the cooling resulted in the development of two-arm spiral structures which persisted as long as cooling continued. An experiment was performed to induce spiral structure in a galaxy by means of the close passage of a companion galaxy. Parameters similar to those expected for M51 and its companion were used. It was found that because of the high velocity dispersion of the disturbed disk galaxy, only a weak two-arm spiral structure appeared. The evolution of a uniformly rotating disk galaxy which is a stationary solution of the collisionless Boltzmann equation was investigated for various values of the initial rms velocity dispersion. It was found that the disk becomes stable at a value of the velocity dispersion predicted by theory.

Hohl, F.↗

Radiation mechanisms and magnetospheric structure of pulsars

A model of pulsars is outlined. Key early considerations were those which led to the identification of pulsars with neutron stars, the Goldreich-Julian model of pulsar magnetospheres, and the recognition that, in a pulsar magnetosphere, a high energy gamma ray may annihilate to produce an electron-positron pair. Arguments that suggest that pulsar magnetospheres may contain large masses of plasma, a suggestion which has important implications concerning the structure of the magnetosphere, are considered along with observational data which support a magnetosphere model based on that idea rather than the Goldreich-Julian model.

Sturrock, P. A.↗

Towards a Consolidated Approach for the Assessment of Evaluation Models of Nuclear Power Reactors

The STARS project at the Paul Scherrer Institut (PSI) has adopted the TRACE thermal-hydraulic (T-H) code for best-estimate system transient simulations of the Swiss Light Water Reactors (LWRs). For analyses involving interactions between system and core, a coupling of TRACE with the SIMULATE-3K (S3K) LWR core simulator has also been developed. In this configuration, the TRACE code and associated nuclear power reactor simulation models play a central role to achieve a comprehensive safety analysis capability. Thus, efforts have now been undertaken to consolidate the validation strategy by implementing a more rigorous and structured assessment approach for TRACE applications involving either only system T-H evaluations or requiring interfaces to e.g. detailed core or fuel behavior models. The first part of this paper presents the preliminary concepts of this validation strategy. The principle is to systematically track the evolution of a given set of predicted physical Quantities of Interest (QoIs) over a multidimensional parametric space where each of the dimensions represent the evolution of specific analysis aspects, including e.g. code version, transient specific simulation methodology and model "nodalisation". If properly set up, such environment should provide code developers and code users with persistent (less affected by user effect) and quantified information (sensitivity of QoIs) on the applicability of a simulation scheme (codes, input models, methodology) for steady state and transient analysis of full LWR systems. Through this, for each given transient/accident, critical paths of the validation process can be identified that could then translate into defining reference schemes to be applied for downstream predictive simulations. In order to illustrate this approach, the second part of this paper presents a first application of this validation strategy to an inadvertent blowdown event that occurred in a Swiss BWR/6. The transient was initiated by the spurious actuation of the Automatic Depressurization System (ADS). The validation approach progresses through a number of dimensions here: First, the same BWR system simulation model is assessed for different versions of the TRACE code, up to the most recent one. The second dimension is the "nodalisation" dimension, where changes to the input model are assessed. The third dimension is the "methodology" dimension. In this case imposed power and an updated TRACE core model are investigated. For each step in each validation dimension, a common set of QoIs are investigated. For the steady-state results, these include fuel temperatures distributions. For the transient part of the present study, the evaluated QoIs include the system pressure evolution and water carry-over into the steam line.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Standing on the shoulders of dwarfs: the kepler asteroseismic legacy sample. I. Oscillation mode parameters

The advent of space-based missions like Kepler has revolutionized the study of solar-type stars, particularly through the measurement and modeling of their resonant modes of oscillation. Here we analyze a sample of 66 Kepler main-sequence stars showing solar-like oscillations as part of the Kepler seismic LEGACY project. We use Kepler short-cadence data, of which each star has at least 12 months, to create frequency-power spectra optimized for asteroseismology. For each star, we identify its modes of oscillation and extract parameters such as frequency, amplitude, and line width using a Bayesian Markov chain Monte Carlo “peak-bagging” approach. We report the extracted mode parameters for all 66 stars, as well as derived quantities such as frequency difference ratios, the large and small separations Δν and δν{sub 02}; the behavior of line widths with frequency and line widths at ν{sub max} with T{sub eff}, for which we derive parametrizations; and behavior of mode visibilities. These average properties can be applied in future peak-bagging exercises to better constrain the parameters of the stellar oscillation spectra. The frequencies and frequency ratios can tightly constrain the fundamental parameters of these solar-type stars, and mode line widths and amplitudes can test models of mode damping and excitation.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Orbit for the Eclipsing Binary V577 Oph

Pulsating stars in eclipsing binary systems are unique objects for providing constraints on stellar models. To fully leverage the information available from the binary system, full orbital radial velocity curves must be obtained. We report 23 radial velocities for components of the eclipsing binary V577 Oph, whose primary star is a δ Sct variable. The velocities cover a nearly complete orbit and a time base of 20 years. We computed orbital elements for the binary and compared them to the ephemeris computed by Creevey et al. The comparison shows marginally different results. In particular, a change in the systemic velocity by −2 km s{sup −1} is suggested by our results. We compare this systemic velocity difference to that expected due to reflex motion of the binary in response to the third body in the system. The systemic velocity difference is consistent with reflex motion, given our mass determination for the eclipsing binary and the orbital parameters determined by Volkov and Volkova for the three-body orbit. We see no evidence for the third body in our spectra, but we do see strong interstellar Na D lines that are consistent in strength with the direction and expected distance of V577 Oph.

79 ASTRONOMY AND ASTROPHYSICS↗

Kepler Data Analysis: Non-Gaussian Noise and Fourier Gaussian Process Analysis of Stellar Variability

We develop a statistical analysis model of Kepler stellar flux data in the presence of planet transits, non-Gaussian noise, and stellar variability. We first develop a model for the Kepler noise probability distribution in the presence of outliers, which make the noise probability distribution non-Gaussian. We develop a signal likelihood analysis based on this probability distribution, in which we model the signal as a sum of the star variability and planetary transits. We argue that these components need to be modeled together if optimal signal is to be extracted from the data. For the stellar variability model we develop an optimal Gaussian process analysis using a Fourier-based Wiener filter approach, where the power spectrum is non-parametric and learned from the data. We develop high dimensional optimization of the objective function, where we jointly optimize all the model parameters, including thousands of star variability modes, and planet transit parameters. We apply the method to Kepler-90 data and show that it gives a better match to the stellar variability than the existing methods, and robustly handles noise outliers. As a consequence, the planet radii have a higher value than what the existing methods give, including splines and celerite.

79 ASTRONOMY AND ASTROPHYSICS↗

X-RAY BURST OSCILLATIONS: FROM FLAME SPREADING TO THE COOLING WAKE

Type I X-ray bursts are thermonuclear flashes observed from the surfaces of accreting neutron stars (NSs) in low mass X-ray binaries. Oscillations have been observed during the rise and/or decay of some of these X-ray bursts. Those seen during the rise can be well explained by a spreading hot spot model, but large amplitude oscillations in the decay phase remain mysterious because of the absence of a clear-cut source of asymmetry. To date there have not been any quantitative studies that consistently track the oscillation amplitude both during the rise and decay (cooling tail) of bursts. Here we compute the light curves and amplitudes of oscillations in X-ray burst models that realistically account for both flame spreading and subsequent cooling. We present results for several such “cooling wake” models, a “canonical” cooling model where each patch on the NS surface heats and cools identically, or with a latitude-dependent cooling timescale set by the local effective gravity, and an “asymmetric” model where parts of the star cool at significantly different rates. We show that while the canonical cooling models can generate oscillations in the tails of bursts, they cannot easily produce the highest observed modulation amplitudes. Alternatively, a simple phenomenological model with asymmetric cooling can achieve higher amplitudes consistent with the observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Using machine learning to predict the correlation of spectra using SDSS magnitudes as an improvement to the Locus Algorithm

The Locus Algorithm is a new technique to improve the quality of differential photometry by optimising the choices of reference stars. At the heart of this algorithm is a routine to assess how good each potential reference star is by comparing its SDSS magnitude values to those of the target star. In this way, the difference in wavelength-dependent effects of the Earth’s atmospheric scattering between target and reference can be minimised. This paper sets out a new way to estimate the quality of each reference star using machine learning. A random subset of stars from SDSS with spectra was chosen. For each one, a suitable reference star, also with a spectrum, was chosen. The correlation between the two spectra in the SDSS r band (between 550 nm and 700 nm) was taken to be the gold-standard measure of how well they match up for differential photometry. The five SDSS magnitude values for each of these stars were used as predictors. A number of supervised machine learning models were constructed on a training set of the stars and were each evaluated on a testing set. The model using Support Vector Regression had the best performance of these models. It was then tested on a final, hold-out, validation set of stars to get an unbiased measure of its performance. With an R 2 of 0.62, the SVR model presents enhanced performance for the Locus Algorithm technique.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on Quenching of Z ≲ 2 Massive Galaxies from the Evolution of the Average Sizes of Star-forming and Quenched Populations in COSMOS

We use >9400 log(m/M{sub ⊙})>10 quiescent and star-forming galaxies at z ≲ 2 in COSMOS/UltraVISTA to study the average size evolution of these systems, with focus on the rare ultra-massive population at log(m/M{sub ⊙})>11.4. The large 2 square degree survey area delivers a sample of ∼400 such ultra-massive systems. Accurate sizes are derived using a calibration based on high-resolution images from the Hubble Space Telescope. We find that at these very high masses, the size evolution of star-forming and quiescent galaxies is almost indistinguishable in terms of normalization and power-law slope. We use this result to investigate possible pathways of quenching massive m > M* galaxies at z < 2. We consistently model the size evolution of quiescent galaxies from the star-forming population by assuming different simple models for the suppression of star formation. These models include an instantaneous and delayed quenching without altering the structure of galaxies and a central starburst followed by compaction. We find that instantaneous quenching reproduces the observed mass-size relation of massive galaxies at z > 1 well. Our starburst+compaction model followed by individual growth of the galaxies by minor mergers is preferred over other models without structural change for log(m/M{sub ⊙})>11.0 galaxies at z > 0.5. None of our models is able to meet the observations at m > M* and z < 1 without significant contribution of post-quenching growth of individual galaxies via mergers. We conclude that quenching is a fast process in galaxies with m ≥ 10{sup 11} M {sub ⊙}, and that major mergers likely play a major role in the final steps of their evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Is There a Metallicity Ceiling to Form Carbon Stars? - A Novel Technique Reveals a Scarcity of C-Stars in the Inner M31 Disk

We use medium-band near-infrared (NIR) Hubble Space Telescope WFC3 photometry with model NIR spectra of Asymptotic Giant Branch (AGB) stars to develop a new tool for efficiently distinguish- ing carbon-rich (C-type) AGB stars from oxygen-rich (M-type) AGB stars in galaxies at the edge of and outside the Local Group. We present the results of a test of this method on a region of the inner disk of M31, where we nd a surprising lack of C stars, contrary to the ndings of previous C star searches in other regions of M31. We nd only 1 candidate C star (plus up to 6 additional, less certain C stars candidates), resulting in an extremely low ratio of C to M stars (C=M = (3.3(sup +20)(sub - 0.1) x 10(sup -4)) that is 1-2 orders of magnitude lower than other C/M estimates in M31. The low C/M ratio is likely due to the high metallicity in this region which impedes stars from achieving C/O > 1 in their atmospheres. These observations provide stringent constraints to evolutionary models of metal-rich AGB stars and suggest that there is a metallicity threshold above which M stars are unable to make the transition to C stars, dramatically affecting AGB mass loss and dust production and, consequently, the observed global properties of metal-rich galaxies.

Metallicity↗

The metal-poor stellar halo in RAVE-TGAS and its implications for the velocity distribution of dark matter

The local velocity distribution of dark matter plays an integral role in interpreting the results from direct detection experiments. We previously showed that metal-poor halo stars serve as excellent tracers of the virialized dark matter velocity distribution using a high-resolution hydrodynamic simulation of a Milky Way-like halo. In this paper, we take advantage of the first Gaia data release, coupled with spectroscopic measurements from the RAdial Velocity Experiment (RAVE), to study the kinematics of stars belonging to the metal-poor halo within an average distance of ~5 kpc of the Sun. We study stars with iron abundances [Fe/H] < -1.5 and -1.8 that are located more than 1.5 kpc from the Galactic plane. Using a Gaussian mixture model analysis, we identify the stars that belong to the halo population, as well as some kinematic outliers. Here we find that both metallicity samples have similar velocity distributions for the halo component, within uncertainties. Assuming that the stellar halo velocities adequately trace the virialized dark matter, we study the implications for direct detection experiments. The Standard Halo Model, which is typically assumed for dark matter, is discrepant with the empirical distribution by ~6σ, predicts fewer high-speed particles, and is anisotropic. As a result, the Standard Halo Model overpredicts the nuclear scattering rate for dark matter masses below ~10 GeV. The kinematic outliers that we identify may potentially be correlated with dark matter substructure, though further study is needed to establish this correspondence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Corequake and shock heating model of the 5 March 1979 gamma-ray burst

Ramaty et al. (1980) have proposed a model of a neutron-star corequake and subsequent shock heating of the atmosphere to account for the gamma-ray burst of March 5, 1979. This model is elaborated by examining the overall energetics and characteristics of the radiation-dominated gas shocks under the assumption of thermodynamic equilibrium, taking into account the electron-positron pair production behind the shock. Using values for the density typical of those expected for neutron-star crusts (100-10,000 g/cu cm) and shock velocities characteristic of the sound speed (40,000-100,000), shock luminosities are obtained comparable to those required if the burst originated in the LMC. In addition, the fraction of energy deposited in electron-positron pairs is in good agreement with observation. Uneven shock heating between the polar and equatorial regions of the neutron star due to the presence of a dipole magnetic field, coupled with a rotation rate of 8 s, can also naturally account for the pulsating phase of the burst.

Ellison, D. C.↗

Stability of charged strange quark stars

We investigate the hydrostatic equilibrium and the stability of charged stars made of a charged perfect fluid. The matter contained in the star follows the MIT bag model equation of state and the charge distribution to a power-law of the radial coordinate. The hydrostatic equilibrium and the stability of charged strange stars are analyzed using the Tolman-Oppenheimer-Volkoff equation and the Chandrasekhar’s equation pulsation, respectively. These two equation are modified from their original form to the inclusion of the electric charge. We found that the stability of the star decreases with the increment of the central energy density and with the increment of the amount of charge.

79 ASTRONOMY AND ASTROPHYSICS↗

The Rapid Evolution of the Exciting Star of the Stingray Nebula

Context: SAO244567, the exciting star of the Stingray nebula, is rapidly evolving. Previous analyses suggested that it has heated up from an effective temperature of about 21 kK in 1971 to over 50 kK in the 1990s. Canonical post-asymptotic giant branch evolution suggests a relatively high mass while previous analyses indicate a low-mass star. Aims: A comprehensive model-atmosphere analysis of UV and optical spectra taken during 1988-2006 should reveal the detailed temporal evolution of its atmospheric parameters and provide explanations for the unusually fast evolution. Methods: Fitting line profiles from static and expanding non-LTE model atmospheres to the observed spectra allowed us to study the temporal change of effective temperature, surface gravity, mass-loss rate, and terminal wind velocity. In addition, we determined the chemical composition of the atmosphere. Results: We find that the central star has steadily increased its effective temperature from 38 kK in 1988 to a peak value of 60 kK in 2002. During the same time, the star was contracting, as concluded from an increase in surface gravity from log g = 4.8 to 6.0 and a drop in luminosity. Simultaneously, the mass-loss rate declined from log(M/M (solar mass) yr (exp −1)) = −9.0 to −11.6 and the terminal wind velocity increased from v (infinity) = 1800 km s (exp −1) to 2800 km s (exp −1). Since around 2002, the star stopped heating and has cooled down again to 55 kK by 2006. It has a largely solar surface composition with the exception of slightly subsolar carbon, phosphorus, and sulfur. The results are discussed by considering different evolutionary scenarios. Conclusions: The position of SAO244567 in the log T (sub eff) -log g plane places the star in the region of sdO stars. By comparison with stellar-evolution calculations, we confirm that SAO244567 must be a low-mass star (M < 0.55 solar mass). However, the slow evolution of the respective stellar evolutionary models is in strong contrast to the observed fast evolution and the young planetary nebula with a kinematical age of only about 1000 years. We speculate that the star could be a late He-shell flash object. Alternatively, it could be the outcome of close-binary evolution. Then SAO244567 would be a low-mass (0.354 solar mass) helium pre-white dwarf after the common-envelope phase, during which the planetary nebula was ejected.

Stingray↗

High-to-Low Activities across VVI, AMA and PHI

This milestone report documents work performed in CASL that was broadly focused on the use of validated high-fidelity simulation codes to inform, improve and calibrate low-fidelity codes to facilitate multiphysics coupling, design, performance analysis, and risk assessment. This high-to-low (Hi2Lo) initiative has the following objectives: employ physics encapsulated in high-fidelity codes to improve the functionality of low-fidelity codes, quantify biases in low-fidelity codes, develop and calibrate physical models and closure relations so that low-fidelity codes more accurately incorporate requisite physics, and develop low-fidelity codes that can be effectively coupled to simulate multiphysics behavior. In combination, these efforts are designed to impact both CIPS and DNB. For example, the development and experimental design-based Bayesian calibration of dynamic gas conductance models in CTFFuel will provide the capability to better model crud-induced power shifts. Whereas the present emphasis has focused on the high-fidelity codes STAR-CCM+ and BISON and low-fidelity codes CTF and CTFFuel, future emphasis will be expanded to include MAMBA. During the reporting period, Hi2Lo efforts focused on two broad components: demonstration of the experimental design-based Bayesian calibration framework to calibrate CTF using both STAR-generated synthetic data and experimental data, and the use of STAR to inform and improve the turbulence modeling and grid heat transfer modeling capabilities in CTF. A complementary milestone focused on the implementation and verification of fuel models in CTFFuel, for which we demonstrated initial calibration results. Future Hi2Lo work will focus on the use of the experimental design-based Bayesian calibration framework, demonstrated this year for STAR and CTF, to calibrate parameters in the fuel swelling, fuel irradiation-induced densification, and fuel relocation models in CTFFuel using both experimental and BISON-generated data. A second component will focus on the continued use of STAR to improve grid modeling capabilities in CTF.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

X-Ray Burst Oscillations: From Flame Spreading to the Cooling Wake

Type I X-ray bursts are thermonuclear flashes observed from the surfaces of accreting neutron stars (NSs) in low mass X-ray binaries. Oscillations have been observed during the rise and/or decay of some of these X-ray bursts. Those seen during the rise can be well explained by a spreading hot spot model, but large amplitude oscillations in the decay phase remain mysterious because of the absence of a clear-cut source of asymmetry. To date there have not been any quantitative studies that consistently track the oscillation amplitude both during the rise and decay (cooling tail) of bursts. Here we compute the light curves and amplitudes of oscillations in X-ray burst models that realistically account for both flame spreading and subsequent cooling. We present results for several such "cooling wake" models, a "canonical" cooling model where each patch on the NS surface heats and cools identically, or with a latitude-dependent cooling timescale set by the local effective gravity, and an "asymmetric" model where parts of the star cool at significantly different rates. We show that while the canonical cooling models can generate oscillations in the tails of bursts, they cannot easily produce the highest observed modulation amplitudes. Alternatively, a simple phenomenological model with asymmetric cooling can achieve higher amplitudes consistent with the observations.

ISM: individual objects (G32.8-0.1↗