Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stellar oscillations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Systematics in asteroseismic modelling: application of a correlated noise model for oscillation frequencies

ABSTRACT The detailed modelling of stellar oscillations is a powerful approach to characterizing stars. However, poor treatment of systematics in theoretical models leads to misinterpretations of stars. Here, we propose a more principled statistical treatment for the systematics to be applied to fitting individual mode frequencies with a typical stellar model grid. We introduce a correlated noise model based on a Gaussian process (GP) kernel to describe the systematics given that mode frequency systematics are expected to be highly correlated. We show that tuning the GP kernel can reproduce general features of frequency variations for changing model input physics and fundamental parameters. Fits with the correlated noise model better recover stellar parameters than traditional methods that either ignore the systematics or treat them as uncorrelated noise.

Li, Tanda (ORCID:0000000163962563)↗
DOE OSTI2025

Investigation of oscillations above the acoustic cut-off frequency for a sample of Kepler stars. II. Identifying temporal frequency shifts

ABSTRACT Stellar oscillations with frequencies greater than the acoustic cut-off frequency are not trapped within the star’s interior. Geometric interference between these high-frequency waves produces a peak-like structure in the power spectrum of the star, known as the pseudo-modes. In this study, we aim to analyse 51 stars previously identified (from a sample of 91 Kepler stars) to contain statistically significant evidence for pseudo-modes, to determine if the pseudo-mode frequencies vary in time. We obtained these frequency variations and the uncertainty by the resampled periodogram approach, where periodogram realisations were created from successive, overlapping 90-d time segments, and a cross-correlation function was generated between them. The computed shifts were compared against temporal p-mode frequency shifts and a stellar magnetic activity proxy, $S_{\mathrm{ph}}$. We found that for 16 per cent of the stars pseudo-mode frequency shifts were significantly anticorrelated with p-mode shifts, as is the case for the Sun. However, we also found pseudo-mode and p-mode shifts to be significantly correlated in-phase in 8 per cent of our sample. The remainder showed no significant correlation. We also searched for trends between the maximum variation of pseudo-mode frequency shifts and stellar parameters. Whilst we found no correlation between the maximum variation and stellar photospheric and chromospheric proxies, we found frequency variations to have a greater amplitude for cooler and slower rotating stars, which is in opposition to p-mode frequency shift variation. Temporal changes to pseudo-mode frequencies may be attributable to variations in stellar magnetic activity, and therefore can contribute to our understanding of stellar activity cycles.

Millson, L. J. (ORCID:0009000342542676)↗
DOE OSTI2021

Radial-velocity Precision of ESPRESSO Through the Analysis of the Solar Twin HIP 11915

Different stellar phenomena affect radial velocities (RVs), causing variations large enough to make it difficult to identify planet signals from the stellar variability. RV variations caused by stellar oscillations and granulation can be reduced through some methods, but the impact of rotationally modulated magnetic activity on RV, due to stellar active regions, is harder to correct. New instrumentation promises an improvement in precision of one order of magnitude, from about 1 m s{sup -1} to about 10 cm s{sup -1}. In this context, we report our first results from 24 spectroscopic Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO)/Very-Large Telescope observations of the solar twin star HIP 11915, spread over 60 nights. We used a Gaussian Process approach and found for HIP 11915 a RV residual rms scatter of about 20 cm s{sup -1}, representing an upper limit for the performance of ESPRESSO.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Fluid Mixing during Phase Separation in Crystallizing White Dwarfs

Abstract Accurate models of cooling white dwarfs must treat the energy released as their cores crystallize. This phase transition slows the cooling by releasing latent heat and also gravitational energy, which results from phase separation: liquid C is released from the solid C/O core, driving an outward carbon flux. The Gaia color–magnitude diagram provides striking confirmation of this theory by revealing a mass-dependent overdensity of white dwarfs, indicating slowed cooling at the expected location. However, the observed overdensity is enhanced relative to the models. Additionally, it is associated with increased magnetism, suggesting a link between crystallization and magnetic field generation. Recent works aimed at explaining an enhanced cooling delay and magnetic field generation employ a uniform mixing prescription that assumes large-scale turbulent motions; we show here that these calculations are not self-consistent. We also show that thermohaline mixing is most likely efficient enough to provide the required chemical redistribution during C/O phase separation, and that the resulting velocities and mixing lengths are much smaller than previous estimates. These reduced fluid motions cannot generate measurable magnetic fields, suggesting any link with crystallization needs to invoke a separate mechanism. Finally, this mixing alters the chemical profiles, which in turn affects the frequencies of the pulsation modes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Properties of F Stars with Stable Radial Velocity Timeseries: A Useful Metric for Selecting Low-jitter F Stars

In a companion paper, we have conducted an in-depth analysis of radial velocity jitter of over 600 stars, examining the astrophysical origins including stellar granulation, oscillation, and magnetic activity. In this paper, we highlight a subsample of those stars, specifically the main sequence and “retired” F stars—which we refer to as “MSRF” stars—that show low levels of RV jitter (<10 m s{sup −1}). We describe the observational signatures of these stars that allow them to be identified in radial velocity planet programs, for instance, those performing follow-up of transiting planets discovered by TESS. We introduce a “jitter metric” that combines the two competing effects of RV jitter with age: activity and convection. Using thresholds in the jitter metric, we can select both “complete” and “pure” samples of low jitter F stars. We also provide recipes for identifying these stars using only Gaia colors and magnitudes. Finally, we describe a region in the Gaia color–magnitude diagram where low jitter F stars are most highly concentrated. By fitting a ninth-order polynomial to the Gaia main sequence, we use the height above the main sequence as a proxy for evolution, allowing for a crude selection of low jitter MSRF stars when activity measurements are otherwise unavailable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Limits on Mode Coherence in Pulsating DA White Dwarfs Due to a Nonstatic Convection Zone

The standard theory of pulsations deals with the frequencies and growth rates of infinitesimal perturbations in a stellar model. Modes which are calculated to be linearly driven should increase their amplitudes exponentially with time; the fact that nearly constant amplitudes are usually observed is evidence that nonlinear mechanisms inhibit the growth of finite amplitude pulsations. Models predict that the mass of convection zones in pulsating hydrogen-atmosphere (DAV) white dwarfs is very sensitive to temperature (i.e., M CZ α T –90 eff ), leading to the possibility that even low-amplitude pulsators may experience significant nonlinear effects. In particular, the outer turning point of finite-amplitude g-mode pulsations can vary with the local surface temperature, producing a reflected wave that is out of phase with what is required for a standing wave. This can lead to a lack of coherence of the mode and a reduction in its global amplitude. In this paper we show that: (1) whether a mode is calculated to propagate to the base of the convection zone is an accurate predictor of its width in the Fourier spectrum, (2) the phase shifts produced by reflection from the outer turning point are large enough to produce significant damping, and (3) amplitudes and periods are predicted to increase from the blue edge to the middle of the instability strip, and subsequently decrease as the red edge is approached. Furthermore, this amplitude decrease is in agreement with the observational data while the period decrease has not yet been systematically studied.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Pulsating Helium-atmosphere White Dwarfs. I. New DBVs from the Sloan Digital Sky Survey

Abstract We present a dedicated search for new pulsating helium-atmosphere (DBV) white dwarfs from the Sloan Digital Sky Survey using the McDonald 2.1 m Otto Struve Telescope. In total we observed 55 DB and DBA white dwarfs with spectroscopic temperatures between 19,000 and 35,000 K. We find 19 new DBVs and place upper limits on variability for the remaining 36 objects. In combination with previously known DBVs, we use these objects to provide an update to the empirical extent of the DB instability strip. With our sample of new DBVs, the red edge is better constrained, as we nearly double the number of DBVs known between 20,000 and 24,000 K. We do not find any new DBVs hotter than PG 0112+104, the current hottest DBV is at T eff ≈ 31,000 K, but do find pulsations in four DBVs with temperatures between 27,000 and 30,000 K, improving empirical constraints on the poorly defined blue edge. We investigate the ensemble pulsation properties of all currently known DBVs, finding that the weighted mean period and total pulsation power exhibit trends with effective temperature that are qualitatively similar to the pulsating hydrogen-atmosphere white dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Confrontation of Observation and Theory: High-frequency QPOs in X-Ray Binaries, Tidal Disruption Events, and Active Galactic Nuclei

In this paper, we compile observations of high-frequency quasi-periodic oscillations (QPOs) around black holes, both stellar and supermassive, and compare their positions in the parameter space of black hole mass, spin, and oscillation frequency. We find that supermassive black holes occupy a separate region of parameter space than stellar, and further, that QPOs seen around tidal disruption events rather than Seyfert-type active galactic nuclei (AGN) occupy an entirely different space. We then compare these results to the orbital resonance, diskoseismic, warped disk, and disk–jet coupling theoretical models for the origin of high-frequency QPOs. We find that while oscillations around stellar mass black holes are generally consistent with the above models, supermassive black holes are decidedly not. Oscillations seen in tidal disruption events are consistent with oscillations near the frequency of the innermost stable circular orbit, while QPOs in AGN are not accounted for by any of the physical models in consideration. This indicates that despite the scale invariance of accretion processes implied by a decades-wide correlation between QPO frequency and black hole mass, any theory of high-frequency QPOs must relate the frequency to more than just the mass and spin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators

The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate kx, but the oscillation and damping time scales are longer in local representations, particularly for small kx approaching the system size.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Review of uncertainties in the cosmic supernova relic neutrino background

We review the computation of and associated uncertainties in the current understanding of the relic neutrino background due to core-collapse supernovae, black hole formation and neutron star merger events. We consider the current status of uncertainties due to the nuclear equation of state (EoS), the progenitor masses, the source supernova neutrino spectrum, the cosmological star formation rate, the stellar initial mass function, neutrino oscillations, and neutrino self-interactions. We summarize the current viability of future neutrino detectors to distinguish the nuclear EoS and the temperature of supernova neutrinos via the detected relic supernova neutrino spectrum.

Astronomy & Astrophysics↗
DOE OSTI2020

Probing non-standard neutrino interactions with supernova neutrinos at Hyper-K

Non-standard neutrino self interactions (NSSI) could be stronger than Fermi interactions. We investigate the ability to constrain these four-neutrino interactions by their effect on the flux of neutrinos originating from a galactic supernova. In the dense medium of a core collapse supernova, these new self interactions can have a significant impact on neutrino oscillations, leading to changes at the flavor evolution and spectra level. We use simulations of the neutrino flux from a 13 solar mass, core collapse supernova at 10 kpc away, and numerically propagate these neutrinos through the stellar medium taking into account vacuum/MSW oscillations, SM ν – ν scattering as well as ν – ν interactions that arise from NSSI. We pass the resulting neutrino flux to a simulation of the future Hyper-Kamiokande detector to see what constraints on NSSI parameters are possible when the next galactic supernova becomes visible. We find that these constraints depend strongly on the neutrino mass hierarchy and if the NSSI is flavor-violating or preserving. Sensitivity to NSSI in the normal hierarchy (NH) at Hyper-K is limited by the experiment’s ability to efficiently detect ν e , but deviations from no NSSI could be seen if the NSSI is particularly strong. In the inverted hierarchy (IH) scenario, Hyper-K can significantly improve constraints on flavor-violating NSSI down to O(10 –1 )G F .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

helioLuna: Prospect for cislunar spacecraft and near-earth asteroid detection using heliostat fields at night

We experimentally and computationally investigate a proposed frequency-domain method for detecting and tracking cislunar spacecraft and near-earth asteroids using heliostat fields at night. Unlike imaging, which detects spacecraft and asteroids by their streak in sidereally-fixed long-exposure photographs, our proposed detection method oscillates the orientation of heliostats concentrating light from the stellar field and measures the light’s photocurrent power spectrum at sub-milliHertz resolution. If heliostat oscillation traces out an ellipse fixed in the galactic coordinate system, spacecraft or asteroids produce a peak in photocurrent power spectrum at a frequency slightly shifted from the starlight peak. The frequency shift is on the scale of milliHertz and proportional to apparent angular rate relative to sidereal. Relative phase corresponds to relative angular position, enabling tracking. A potential advantage of this frequency-domain method over imaging is that detectivity improves with apparent angular rate and number of heliostats. Since heliostats are inexpensive compared to an astronomical observatory and otherwise unused at night, the proposed method may cost-effectively augment observatory systems such as NASA’s Asteroid Terrestrial-impact Last Alert System (ATLAS).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2020

Chromium Nucleosynthesis and Silicon–Carbon Shell Mergers in Massive Stars

In this work, we analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] ≈ 0 reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 ${M}_{\odot }$ model at metallicity Z = 0.01, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including 51 V and 52 Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 ${M}_{\odot }$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 M⊙ model at Z = 0.01 in the yields interpolation scheme of our GCE model for stars between 15 and 25 ${M}_{\odot }$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] ≈ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si–C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the presupernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si shell could launch a convective–reactive global oscillation if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Generalized Mercier stability criterion for stellarators

The Mercier criterion is a well-known stability criterion for tokamaks. It is derived from a 2 × 2 matrix eigenvalue problem arising from the expansion of resonant solutions about a singular surface where m−nq=0, with m and n being the poloidal and toroidal mode numbers, respectively, and q being the safety factor. The stability criterion is that the eigenvalues must be real, otherwise, the solution oscillates, violating the Newcomb crossing criterion. Because of the non-axisymmetry of stellarators, different toroidal as well as poloidal harmonics couple to each other. It follows that each singular surface can have multiple resonant harmonics, with multiplicity M≥1. The corresponding matrix eigenvalue problem involves a 2M×2M matrix, resulting in M pairs of positive and negative eigenvalues. The generalized stability criterion is that all eigenvalues must be real. While the original Mercier criterion can be expressed in terms of quadratures of equilibrium quantities over the singular surface, which can be evaluated anywhere, the generalized Mercier criterion can only be evaluated on rational q surfaces with a given set of resonant harmonics.

Physics↗