Search NASA⌕ Search

SEARCH · Search NASA

Results for “ORBITAL ELEMENT”

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

Predicting Cislunar Orbit Lifetimes from Initial Orbital Elements

The volume of space between Earth’s geosynchronous orbit out to the Moon’s sphere of influence, including the lunar Lagrange points, is crucial for the successful planning and execution of space missions, but not fully understood dynamically. This region is a part of cislunar space. Trajectories through cislunar space are influenced by the gravitational forces of the Sun, Earth, Moon, and other Solar System planets leading to typically unpredictable and chaotic trajectory behavior. It is therefore difficult to predict the stability of an trajectory through cislunar space from a set of initial conditions or orbital elements. We simulate one million cislunar orbits to train a self-organizing map (SOM) to cluster the trajectories and orbits into families based on how long they remain stable within the cislunar space. Using the trained SOM, we are able to predict the stable lifetime of a trajectory through cislunar space from a set of initial orbital elements to within an accuracy of 10% for 8% of simulated trajectories and within 50% for 43% of the simulated trajectories. Clustering in the SOM suggests that a variety of trajectory morphologies have similar lifetimes. Once trained, the SOM can predict the stable lifetime of a given cislunar trajectory within milliseconds. The methods developed in this work enable the rapid identification of stable cislunar orbits and trajectories that could be used for future space exploration. Moreover, the developed SOM method can generate orbital and trajectory lifetime estimates from minimal observational data, such as a single two line element, making it useful for early warning systems and large-scale sensor network operations.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Orbital Elements of the Two-body Problem with Slowly Decreasing Mass: The Gyldén–Mestchersky Cases

We study the dynamics of two-body problems with gravitational parameters μ = μ(t) decreasing with the time. In particular we focus our attention on the behavior of the orbital elements in the first two Gyldén–Mestchersky cases: in the first case (GM-I) μ(t)=μ{sub 0}(1+αt){sup −1}, where μ{sub 0} is the constant initial value of the parameter and α a small positive constant. In the second case (GM-II) μ(t)=μ{sub 0}(1+2αt){sup −1/2}. In the GM-I problem, using a suitable transformation of the radius vector and the physical time, which reduces the original problem to a pure Kepler problem, we are able to obtain explicit analytical expressions of the orbital elements that hold true for all t⩾0 and therefore derive their asymptotic behavior. Moreover using the alternative scaled time these expressions show clearly their periodic and secular terms. For the GM-II problem a similar transformation to the GM-I problem leads to a perturbed Kepler problem with a radial direction perturbation. Although exact solutions in terms of elliptic functions may be obtained, the use of some integrals leads to analytical expressions of orbital elements that allow us to study their behavior. Finally, some comparisons between the time evolution of elements in both problems are also included.

79 ASTRONOMY AND ASTROPHYSICS↗

Orbital Elements and Stellar Parameters of the Active Binary UX Arietis

Stellar activity observed as large surface spots, radio flares, or emission lines is often found in binary systems. UX Arietis exhibits these signs of activity, originating on the K0 subgiant primary component. Our aim is to resolve the binary, measure the orbital motion, and provide accurate stellar parameters such as masses and luminosities to aid in the interpretation of the observed phenomena. Using the CHARA six-telescope optical long-baseline array on Mount Wilson, California, we obtained amplitudes and phases of the interferometric visibility on baselines up to 330 m in length, resolving the two components of the binary. We reanalyzed archival Center for Astrophysics spectra to disentangle the binary component spectra and the spectrum of the third component, which was resolved by speckle interferometry. We also obtained new spectra with the Nordic Optical Telescope, and we present new photometric data that we use to model stellar surface spot locations. Both interferometric visibilities and spectroscopic radial velocities are modeled with a spotted primary stellar surface using the Wilson–Devinney code. We fit the orbital elements to the apparent orbit and radial velocity data to derive the distance (52.1 ± 0.8 pc) and stellar masses (M{sub P}=1.30±0.06 M{sub ⊙}, M{sub S}=1.14±0.06 M{sub ⊙}). The radius of the primary can be determined to be R{sub P}=5.6±0.1 R{sub ⊙} and that of the secondary to be R{sub S}=1.6±0.2 R{sub ⊙}. The equivalent spot coverage of the primary component was found to be 62% with an effective temperature 20% below that of the unspotted surface.

79 ASTRONOMY AND ASTROPHYSICS↗

Probabilistic Data Association for Orbital-Element Estimation Using Multistage Expectation–Maximization

Tracking space objects is important for managing space traffic and predicting collisions, but is difficult in part due to data association and orbit model uncertainty. Expectation–maximization (EM) is a commonly used tracking method that has not been widely considered for tracking space objects. The technique consists of iteratively computing data association probabilities with a set of current element estimates, and updating estimates of the elements by solving a nonlinear weighted least-squares regression problem where the weights are the data association probabilities. In this paper, we demonstrate the use of EM for probabilistic data association and orbital-element estimation by applying the technique to simulated data from two angles-only tracking scenarios. In both scenarios, EM provides correct data associations and accurate maximum likelihood estimates of orbital elements. One scenario considers tracking a single object in clutter and quantifies the improvement of the orbital-element estimates and data associations as the detection probability increases. However, standard application of EM requires knowing the number of objects or may fail when a large number of objects are present. To address these issues, this paper employs a multistage version of EM that is applicable when there are a large and possibly unknown number of objects.

97 MATHEMATICS AND COMPUTING↗

Study of the Influence of Some NEA Orbit Elements on Sizes of Chaotic Zones of Resonances with the Earth and Jupiter

Results of investigation of regular and chaotic motions of near-Earth asteroids (NEAs) in the vicinity of the 1:2 mean motion resonance with the Earth and 2:1 mean motion resonance with Jupiter are discussed in the present work. The resonant zones, constructed from the set of initial values of the semimajor axis and the mean anomaly at some fixed values of other orbital elements, are numerically simulated. By means of the MEGNO parameter, the vicinities of the boundaries of these zones are investigated. Chaotic zones are represented graphically, and the widths of these zones are numerically estimated. The behavior of the resonant and chaotic zones is demonstrated for different values of the eccentricity and inclination.

79 ASTRONOMY AND ASTROPHYSICS↗

He IIλ4686 emission from the massive binary system in η Car: Constraints to the orbital elements and the nature of the periodic minima

Eta Carinae (η Car) is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He ii λ4686 line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He ii λ4686 emission line across the 2014.6 periastron passage using ground- and space-based telescopes. Comparison with previous data confirmed the overall repeatability of the line equivalent width (EW), radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW variations, combined with other measurements, yield an orbital period of 2022.7 ± 0.3 days. The observed variability of the EW was reproduced by a model in which the line flux primarily arises at the apex of the wind–wind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary’s atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to 135°–153°, and the longitude of periastron to 234°–252°. It also suggests that periastron passage occurred on T{sub 0}=2456874.4 (±1.3 days). Our model also reproduced EW variations from a polar view of the primary star as determined from the observed He ii λ4686 emission scattered off the Homunculus nebula.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolution of Binary Supermassive Black Holes in Rotating Nuclei

The interaction of a binary supermassive black hole with stars in a galactic nucleus can result in changes to all the elements of the binary’s orbit, including the angles that define its orientation. If the nucleus is rotating, the orientation changes can be large, causing large changes in the binary’s orbital eccentricity as well. We present a general treatment of this problem based on the Fokker–Planck equation for f , defined as the probability distribution for the binary’s orbital elements. First- and second-order diffusion coefficients are derived for the orbital elements of the binary using numerical scattering experiments, and analytic approximations are presented for some of these coefficients. Solutions of the Fokker–Planck equation are then derived under various assumptions about the initial rotational state of the nucleus and the binary hardening rate. We find that the evolution of the orbital elements can become qualitatively different when we introduce nuclear rotation: (1) the orientation of the binary’s orbit evolves toward alignment with the plane of rotation of the nucleus and (2) binary orbital eccentricity decreases for aligned binaries and increases for counteraligned ones. We find that the diffusive (random-walk) component of a binary’s evolution is small in nuclei with non-negligible rotation, and we derive the time-evolution equations for the semimajor axis, eccentricity, and inclination in that approximation. The aforementioned effects could influence gravitational wave production as well as the relative orientation of host galaxies and radio jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Spin–Orbit Matrix Elements for a Combined Spin-Flip and IP/EA approach

We introduce a practical approach for computing the Breit–Pauli spin–orbit matrix elements of multiconfigurational systems with both spin and spatial degeneracies based on our recently developed RAS-nSF-IP/EA method (Houck, S. E.; et al. J. Chem. Theory Comput. 2019, 15, 2278). The spin–orbit matrix elements over all the multiplet components are computed using a single one-particle reduced density matrix as a result of the Wigner–Eckart theorem. A mean field spin–orbit approximation was used to account for the two-electron contributions. Basis set dependence as well as the effect of including additional excitations is presented. The effect of correlating the core and semicore orbitals is also examined. Surprisingly accurate results are obtained for spin–orbit coupling constants, despite the fact that the efficient wave function approximations we research neglect the bulk of dynamical correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TWO STARS TWO WAYS: CONFIRMING A MICROLENSING BINARY LENS SOLUTION WITH A SPECTROSCOPIC MEASUREMENT OF THE ORBIT

Light curves of microlensing events involving stellar binaries and planetary systems can provide information about the orbital elements of the system due to orbital modulations of the caustic structure. Accurately measuring the orbit in either the stellar or planetary case requires detailed modeling of subtle deviations in the light curve. At the same time, the natural, Cartesian parameterization of a microlensing binary is partially degenerate with the microlens parallax. Hence, it is desirable to perform independent tests of the predictions of microlens orbit models using radial velocity (RV) time series of the lens binary system. To this end, we present 3.5 years of RV monitoring of the binary lens system OGLE-2009-BLG-020 L, for which Skowron et al. constrained all internal parameters of the 200–700 day orbit. Our RV measurements reveal an orbit that is consistent with the predictions of the microlens light curve analysis, thereby providing the first confirmation of orbital elements inferred from microlensing events.

79 ASTRONOMY AND ASTROPHYSICS↗

Orbits for 18 Visual Binaries and Two Double-line Spectroscopic Binaries Observed with HRCAM on the CTIO SOAR 4 m Telescope, Using a New Bayesian Orbit Code Based on Markov Chain Monte Carlo

We present orbital elements and mass sums for 18 visual binary stars of spectral types B to K (five of which are new orbits) with periods ranging from 20 to more than 500 yr. For two double-line spectroscopic binaries with no previous orbits, the individual component masses, using combined astrometric and radial velocity data, have a formal uncertainty of ∼0.1 M{sub ⊙}. Adopting published photometry and trigonometric parallaxes, plus our own measurements, we place these objects on an H-R diagram and discuss their evolutionary status. These objects are part of a survey to characterize the binary population of stars in the Southern Hemisphere using the SOAR 4 m telescope+HRCAM at CTIO. Orbital elements are computed using a newly developed Markov chain Monte Carlo (MCMC) algorithm that delivers maximum-likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate the uncertainty of our derived parameters in a robust way. For spectroscopic binaries, using our approach, it is possible to derive a self-consistent parallax for the system from the combined astrometric and radial velocity data (“orbital parallax”), which compares well with the trigonometric parallaxes. We also present a mathematical formalism that allows a dimensionality reduction of the feature space from seven to three search parameters (or from 10 to seven dimensions—including parallax—in the case of spectroscopic binaries with astrometric data), which makes it possible to explore a smaller number of parameters in each case, improving the computational efficiency of our MCMC code.

79 ASTRONOMY AND ASTROPHYSICS↗

Colors of Centaurs observed by the Subaru/Hyper Suprime-Cam and implications for their origin

Abstract Centaurs have orbits between Jupiter and Neptune and are thought to originate from the trans-Neptunian region. Observations of surface properties of Centaurs and comparison with those of trans-Neptunian objects (TNOs) would provide constraints on their origin and evolution. We analyzed imaging data of nine known Centaurs observed by the Hyper Suprime-Cam (HSC) installed on the Subaru Telescope with the g- and i-band filters. Using the data available in the public HSC data archive, as well as those obtained by the HSC Subaru Strategic Program (HSC-SSP) by the end of 2017 June, we obtained the g − i colors of the nine Centaurs. We compared them with those of known TNOs in the HSC-SSP data obtained by T. Terai et al. (2018, PASJ, 70, S40). We found that the color distribution of the nine Centaurs is similar to that of those TNOs with high orbital inclinations, but distinct from those TNOs with low orbital inclinations. We also examined correlations between the colors of these Centaurs and their orbital elements and absolute magnitude. The Centaurs’ colors show a moderate positive correlation with semi-major axis, while no significant correlations between the color and other orbital elements or absolute magnitude were found for these Centaurs. On the other hand, recent studies on Centaurs with larger samples show interesting correlations between their color and absolute magnitude or orbital inclination. We discuss how our data fit in these previous studies, and also discuss implications of these results for their origin and evolution.

Sakugawa, Haruka↗

Evaluating the Dynamical Stability of Outer Solar System Objects in the Presence of Planet Nine

Here, we evaluate the dynamical stability of a selection of outer solar system objects in the presence of the proposed new solar system member Planet Nine. We use a Monte Carlo suite of numerical N-body integrations to construct a variety of orbital elements of the new planet and evaluate the dynamical stability of eight trans-Neptunian objects (TNOs) in the presence of Planet Nine. These simulations show that some combinations of orbital elements ($a,e$) result in Planet Nine acting as a stabilizing influence on the TNOs, which can otherwise be destabilized by interactions with Neptune. These simulations also suggest that some TNOs transition between several different mean-motion resonances during their lifetimes while still retaining approximate apsidal antialignment with Planet Nine. This behavior suggests that remaining in one particular orbit is not a requirement for orbital stability. As one product of our simulations, we present an a posteriori probability distribution for the semimajor axis and eccentricity of the proposed Planet Nine based on TNO stability. This result thus provides additional evidence that supports the existence of this proposed planet. We also predict that TNOs can be grouped into multiple populations of objects that interact with Planet Nine in different ways: one population may contain objects like Sedna and 2012 VP 113 , which do not migrate significantly in semimajor axis in the presence of Planet Nine and tend to stay in the same resonance; another population may contain objects like 2007 TG 422 and 2013 RF 98 , which may both migrate and transition between different resonances.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying the Influence of Jupiter on the Earth’s Orbital Cycles

A wealth of Earth-sized exoplanets will be discovered in the coming years, providing a large pool of candidates from which the targets for the search for life beyond the solar system will be chosen. The target selection process will require the leveraging of all available information in order to maximize the robustness of the target list and make the most productive use of follow-up resources. Here, we present the results of a suite of n-body simulations that demonstrate the degree to which the orbital architecture of the solar system impacts the variability of Earth’s orbital elements. By varying the orbit of Jupiter and keeping the initial orbits of the other planets constant, we demonstrate how subtle changes in solar system architecture could alter the Earth’s orbital evolution—a key factor in the Milankovitch cycles that alter the amount and distribution of solar insolation, thereby driving periodic climate change on our planet. The amplitudes and frequencies of Earth’s modern orbital cycles fall in the middle of the range seen in our runs for all parameters considered—neither unusually fast nor slow, neither large nor small. This finding runs counter to the “Rare Earth” hypothesis, which suggests that conditions on Earth are so unusual that life elsewhere is essentially impossible. Our results highlight how dynamical simulations of newly discovered exoplanetary systems could be used as an additional means to assess the potential targets of biosignature searches, and thereby help focus the search for life to the most promising targets.

79 ASTRONOMY AND ASTROPHYSICS↗

SPECTROSCOPIC ORBITS FOR 15 LATE-TYPE STARS

Spectroscopic orbital elements are determined for 15 stars with periods from 8 to 6528 days with six orbits computed for the first time. Improved astrometric orbits are computed for two stars and one new orbit is derived. Visual orbits were previously determined for four stars, four stars are members of multiple systems, and five stars have Hipparcos “G” designations or have been resolved by speckle interferometry. For the nine binaries with previous spectroscopic orbits, we determine improved or comparable elements. For HD 28271 and HD 200790, our spectroscopic results support the conclusions of previous authors that the large values of their mass functions and lack of detectable secondary spectrum argue for the secondary in each case being a pair of low-mass dwarfs. The orbits given here may be useful in combination with future interferometric and Gaia satellite observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating the Relativistic Motion of the Stars Near the Supermassive Black Hole in the Galactic Center

The S-star cluster in the Galactic center allows us to study the physics close to a supermassive black hole, including distinctive dynamical tests of general relativity. Our best estimates for the mass of and the distance to Sgr A* using the three stars with the shortest period (S2, S38, and S55/S0-102) and Newtonian models are M {sub BH} = (4.15 ± 0.13 ± 0.57) × 10{sup 6} M {sub ⊙} and R {sub 0} = 8.19 ± 0.11 ± 0.34 kpc. Additionally, we aim at a new and practical method to investigate the relativistic orbits of stars in the gravitational field near Sgr A*. We use a first-order post-Newtonian approximation to calculate the stellar orbits with a broad range of periapse distance r {sub p} . We present a method that employs the changes in orbital elements derived from elliptical fits to different sections of the orbit. These changes are correlated with the relativistic parameter defined as ϒ ≡ r {sub s} / r {sub p} (with r {sub s} being the Schwarzschild radius) and can be used to derive ϒ from observational data. For S2 we find a value of ϒ = 0.00088 ± 0.00080, which is consistent, within the uncertainty, with the expected value of ϒ = 0.00065 derived from M {sub BH} and the orbit of S2. We argue that the derived quantity is unlikely to be dominated by perturbing influences such as noise on the derived stellar positions, field rotation, and drifts in black hole mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Orbits and Masses of Binaries from Speckle Interferometry at SOAR

We present results from Speckle inteferometric observations of 15 visual binaries and one double-line spectroscopic binary, carried out with the HRCam Speckle camera of the SOAR 4.1 m telescope. These systems were observed as a part of an on-going survey to characterize the binary population in the solar vicinity, out to a distance of 250 pc. We obtained orbital elements and mass sums for our sample of visual binaries. The orbits were computed using a Markov Chain Monte Carlo algorithm that delivers maximum likelihood estimates of the parameters, as well as posterior probability density functions that allow us to evaluate their uncertainty. Their periods cover a range from 5 yr to more than 500 yr; and their spectral types go from early A to mid M, implying total system masses from slightly more than 4M{sub ⊙} down to 0.2M {sub ⊙}. They are located at distances between approximately 12 and 200 pc, mostly at low Galactic latitude. For the double-line spectroscopic binary YSC8, we present the first combined astrometric/radial-velocity orbit resulting from a self-consistent fit, leading to individual component masses of 0.897 ± 0.027 M {sub ⊙} and 0.857 ± 0.026 M {sub ⊙}; and an orbital parallax of 26.61 ± 0.29 mas, which compares very well with the Gaia DR2 trigonometric parallax (26.55 ± 0.27 mas). In combination with published photometry and trigonometric parallaxes, we place our objects on an H-R diagram and discuss their evolutionary status. We also present a thorough analysis of the precision and consistency of the photometry available for them.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detection of Water Vapor in the Thermal Spectrum of the Non-transiting Hot Jupiter Upsilon Andromedae b

The Upsilon Andromedae system was the first multi-planet system discovered orbiting a main-sequence star. We describe the detection of water vapor in the atmosphere of the innermost non-transiting gas giant ups And b by treating the star–planet system as a spectroscopic binary with high-resolution, ground-based spectroscopy. We resolve the signal of the planet’s motion and break the mass-inclination degeneracy for this non-transiting planet via deep combined flux observations of the star and the planet. In total, seven epochs of Keck NIRSPEC L band observations, three epochs of Keck NIRSPEC short-wavelength K band observations, and three epochs of Keck NIRSPEC long wavelength K band observations of the ups And system were obtained. We perform a multi-epoch cross-correlation of the full data set with an atmospheric model. We measure the radial projection of the Keplerian velocity (K {sub P} = 55 ± 9 km s{sup −1}), true mass (M{sub b}=1.7{sub −0.24}{sup +0.33} M {sub J}), and orbital inclination (i {sub b} 24° ± 4°), and determine that the planet’s opacity structure is dominated by water vapor at the probed wavelengths. Dynamical simulations of the planets in the ups And system with these orbital elements for ups And b show that stable, long-term (100 Myr) orbital configurations exist. These measurements will inform future studies of the stability and evolution of the ups And system, as well as the atmospheric structure and composition of the hot Jupiter.

79 ASTRONOMY AND ASTROPHYSICS↗