SEARCH · Search NASA

Results for “probability amplitudes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Theoretical calculations of ion acceleration in the vicinity of comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked up by the solar wind. The cometary ion pickup process for comet Giacobini-Zinner is studied in two ways: (1) with a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.

Mckenzie, M. L.↗

DOE OSTI2021

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Obliquity histories of Earth and Mars: Influence of inertial and dissipative core-mantle coupling

For both the Earth and Mars, secular variations in the angular separation of the spin axis from the orbit normal are suspected of driving major climatic changes. There is considerable interest in determining the amplitude and timing of these obliquity variations. If the orientation of the orbital plane were inertially fixed, and the planet were to act as a rigid body in it response to precessional torques, the spin axis would simply precess around the orbit at a fixed obliquity and at a uniform angular rate. The precession rate parameter depends on the principal moments of inertia and rotation rate of the perturbed body, and on the gravitational masses and semiminor axes of the perturbing bodies. For Mars, the precession rate is not well known, but probably lies in the interval 8 to 10 arcsec/year. Gravitational interactions between the planets lead to secular motions of the orbit planes. In the rigid body case, the spin axis still attempts to precess about the instantaneous orbit normal, but now the obliquity varies. The hydrostatic figure of a planet represents a compromise between gravitation, which attempts to attain spherical symmetry, and rotation, which prefers cylindrical symmetry. Due to their higher mean densities the cores of the Earth and Mars will be more nearly spherical than the outer layers of these planets. On short time scales it is appropriate to consider the core to be an inviscid fluid constrained to move with the ellipsoidal region bounded by the rigid mantle. The inertial coupling provided by this mechanism is effective whenever the ellipticicy of the container exceeds the ratio of precessional to rotational rates. If the mantle were actually rigid, this would be an extremely effective type of coupling. However, on sufficiently long time scales, the mantle will deform viscously and can accommodate the motions of the core fluid. A fundamentally different type of coupling is provided by electromagnetic or viscous torques. This type of coupling is likely to be most important on longer time scales. In each case, the mantle exerts an equal and opposite torque on the core.

Bills, Bruce G.↗

NASA NTRS1994

Interpretation of microlensing events in Q2237 + 0305

We point out the uncertainties in the potential models of the foreground galaxy in the gravitational lens Q2237 + 0305. The surface densities and external shears are uncertain by a factor of appproximately less than 2. Very long light curves are calculated for different models using the method developed by Witt (1993) and Lewis et al. (1993). We find that the probability of images being in a quiescent state depends strongly on the potential models and the direction of the source motion relative to the alignment of the (gravitational lens) shear term of the macroimages. For example, the probabilities of finding variations smaller than 0.1 mag and for a period of longer than 0.5 Einstein radii (approximately equals 3-5 yr) in image C and D are approximately 30%. This offers an explanation of why image C and D showed relatively small variations. The quiescent images are more likely to be demagnified, consistent with the observations. Further, the higher the magnification of the macroimages, the weaker are the microlensing events. In addition, we conclude that a simple Gaussian profile has some difficulties in explaining the asymmetric shape and the small amplitude of the observed high-magnification events if they are typical. Source models in which the emission is produced by a brightt compact core and a more extended halo can occur better explain the observed double peak events. More data are needed to clearly distinguish different galaxy potential models.

Witt, Hans J.↗

NASA NTRS2006

Relationship Between Vortex Meander and Ambient Turbulence

Efforts are currently underway to increase the capacity of airports by use of closely-spaced parallel runways. If such an objective is to be achieved safely and efficiently during both visual and instrument flight conditions, it will be necessary to develop more precise methods for the prediction of the motion and spread of the hazard posed by the lift-generated vortex-wakes of aircraft, and their uncertainties. The purpose of the present study is to relate the motion induced in vortex filaments by turbulence in the ambient flow field to the measured turbulence in the flow field. The problem came about when observations made in the two largest NASA wind tunnels indicated that extended exposure of vortex wakes to the turbulence in the wind tunnel air stream causes the centers of the vortices to meander about with time at a given downstream station where wake measurements are being made. Although such a behavior was expected, the turbulence level based on the maximum amplitude of meander was much less than the root-mean-squared value measured in the free-stream of the wind tunnel by use of hot-film anemometers. An analysis of the time-dependent motion of segments of vortex filaments as they interact with an eddy, indicates that the inertia of the filaments retards their motion enough in the early part of their travel to account for a large part of the difference in the two determinations of turbulence level. Migration of vortex filaments from one turbulent eddy to another (probably with a different orientation), is believed to account for the remainder of the difference. Methods that may possibly be developed for use in the measurement of the magnitude of the more intense eddies in turbulent flow fields and how they should be adjusted to predict vortex meander are then discussed.

Rossow, Vernon J.↗

NASA NTRS2006

System Would Detect Foreign-Object Damage in Turbofan Engine

A proposed data-fusion system, to be implemented mostly in software, would further process the digitized and preprocessed outputs of sensors in a turbofan engine to detect foreign-object damage (FOD) [more precisely, damage caused by impingement of such foreign objects as birds, pieces of ice, and runway debris]. The proposed system could help a flight crew to decide what, if any, response is necessary to complete a flight safely, and could aid mechanics in deciding what post-flight maintenance action might be needed. The sensory information to be utilized by the proposed system would consist of (1) the output of an accelerometer in an engine-vibration-monitoring subsystem and (2) features extracted from a gas path analysis. ["Gas path analysis" (GPA) is a term of art that denotes comprehensive analysis of engine performance derived from readings of fuel-flow meters, shaft-speed sensors, temperature sensors, and the like.] The acceleration signal would first be processed by a wavelet-transform-based algorithm, using a wavelet created for the specific purpose of finding abrupt FOD-induced changes in noisy accelerometer signals. Two additional features extracted would be the amplitude of vibration (determined via a single- frequency Fourier transform calculated at the rotational speed of the engine), and the rate of change in amplitude due to an FOD-induced rotor imbalance. This system would utilize two GPA features: the fan efficiency and the rate of change of fan efficiency with time. The selected GPA and vibrational features would be assessed by two fuzzy-logic inference engines, denoted the "Gas Path Expert" and the "Vibration Expert," respectively (see Figure 1). Each of these inference engines would generate a "possibility" distribution for occurrence of an FOD event: Each inference engine would assign, to its input information, degrees of membership, which would subsequently be transformed into basic probability assignments for the gas path and vibration components. The outputs of the inference engines would be fused by use of Dempster s combination algorithm (more precisely, an algorithm, based on the Dempster-Shafer-Yager theory of evidence, for fusing uncertain or imprecise information) to provide a reduced body of information to a human or computer decision maker. Figure 2 depicts some outputs generated in response to simulated accelerometer and GPA signals

Torso, James A.↗

NASA NTRS1984

A comparison of the properties of highly polarized QSOs versus low-polarization QSOs

An optical linear-polarization survey testing the relationships between polarization and other properties of QSOs is presented. Polarimetric observations of 239 QSOs have been made. Data from the literature are used to analyze the relationships between polarization and radio, optical, and X-ray properties. The highly polarized QSOs are generally compact radio sources, are associated with radio properties such as low-frequency variability and superluminal motion, exhibit large-amplitude rapid photometric variability, have steep nonthermal optical continua, and may show excess X-ray emission. No relationship is seen between polarization and redshift, optical luminosity, or the equivalent width of emission lines. The results are discussed in the context of both isotropic and anisotropic (beaming) models. While the associations between the optical and radio properties of highly polarized QSOs provide strong motivation for the anisotropic model, the lack of associations with redshift, optical luminosity, and emission-line strength is inconsistent with this model. It is concluded that the highly polarized QSOs are probably fundamentally different in their physical properties from low-polarization QSOs.

Moore, R. L.↗

DOE OSTI2022

Effective field theory analysis of 3 He–α scattering data

We treat low-energy 3 He–α elastic scattering in an effective field theory (EFT) that exploits the separation of scales in this reaction. We compute the amplitude up to next-to-next-to-leading order, developing a hierarchy of the effective-range parameters (ERPs) that contribute at various orders. We use the resulting formalism to analyse data for recent measurements at center-of-mass energies of 0.38–3.12 MeV using the scattering of nuclei in inverse kinematics (SONIK) gas target at TRIUMF as well as older data in this energy regime. We employ a likelihood function that incorporates the theoretical uncertainty due to truncation of the EFT and use Markov chain Monte Carlo sampling to obtain the resulting posterior probability distribution. We find that the inclusion of a small amount of data on the analysing power A y is crucial to determine the sign of the p-wave splitting in such an analysis. Here, the combination of A y and SONIK data constrains all ERPs up to O(p 4 ) in both s- and p-waves quite well. The asymptotic normalisation coefficients and s-wave scattering length are consistent with a recent EFT analysis of the capture reaction 3 He(α, γ) 7 Be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1987

Plasma wave evidence for reflected ions in front of subcritical shocks - ISEE 1 and 2 observations

Plasma wave profiles of bow shock crossings at ion acoustic frequencies exhibit 'feet' or 'pedestals', just as magnetic field profiles do. These ion-acoustic feet, consisting of increases in ion wave amplitude over background, are observed in front of virtually all quasi-perpendicular shocks, regardless of Mach number, but are closely associated with magnetic feet and reflected ions circulating around the ambient upstream field at supercritical shocks. Observed gyroradii of the reflected ions have been found to be commensurate with the thicknesses of the magnetic feet, which consist of increases in average field magnitude just outside the main shock ramps in front of supercritical shocks. Estimated scale lengths (gyroradii) of reflected protons are proportional to the thicknesses of acoustic wave feet in all cases, indicating that such protons are probably universally present at subcritical as well as supercritical shocks, although the protons are observed only sporadically because of instrumental limitations. The mere presence or absence of reflected ions in front of a shock may therefore be an insufficient diagnostic of the shock's status with respect to criticality and ion dissipation.

Greenstadt, E. W.↗

NASA NTRS1996

Effect of Cyclic Thermo-Mechanical Loads on Fatigue Reliability in Polymer Matrix Composites

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress dependent multi-factor interaction relationship developed at NASA Lewis Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability- based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/- 45/90)(sub s) graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Shah, A. R.↗

NASA NTRS1998

Ionospheric Signature of Surface Mine Blasts from Global Positioning System Measurements

Sources such as atmospheric or buried explosions and shallow earthquakes are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic and gravity waves induce variations of the ionospheric electron density. The Global Positioning System (GPS) provides a way of directly measuring the total electron content in the ionosphere and, therefore, of detecting such perturbations in the upper atmosphere. In July and August 1996, three large surface mine blasts (1.5 Kt each) were detonated at the Black Thunder coal mine in eastern Wyoming. As part of a seismic and acoustic monitoring- experiment, we deployed five dual-frequency GPS receivers at distances ranging from 50 to 200 km from the mine and were able to detect the ionospheric perturbation caused by the blasts. The perturbation starts 10 to 15 min after the blast, lasts for about 30 min, and propagates with an apparent horizontal velocity of 1200 meters per second. Its amplitude reaches 3 x 10 (exp 14) el per square meters in the 7-3 min period band, a value close to the ionospheric perturbation caused by the M = 6.7 Northridge earthquake. The small signal-to-noise ratio of the perturbation can be improved by slant-stacking the electron content time-series recorded by the different GPS receivers taking into account the horizontal propagation of the perturbation. The energy of the perturbation is concentrated in the 200 to 300 second period band, a result consistent with previous observations and numerical model predictions. The 300 second band probably corresponds to gravity modes and shorter periods to acoustic modes, respectively. Using a 1-D stratified velocity model of the atmosphere we show that linear acoustic ray tracing fits arrival times at all GPS receivers. We interpret the perturbation as a direct acoustic wave caused by the explosion itself. This study shows that even relatively small subsurface events can produce ionospheric perturbations that are above the detection threshold of the GPS technique. By sensing derivative signals, which can be detected over a relatively broad region, it appears that GPS might be particularly useful for source characterization within the first acoustic quiet zone where infrasound would probably be ineffective. This suggests that dual-frequency GPS monitoring could contribute to Comprehensive Nuclear Test Ban Treaty verification.

Calais, Eric↗

NASA NTRS2010

Probabilistic Simulation for Combined Cycle Fatigue in Composites

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress dependent multifactor interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/- 45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

DOE OSTI2020

Forecasting cosmological constraints from the weak lensing magnification of type Ia supernovae measured by the Nancy Grace Roman Space Telescope

We report that the weak lensing magnification of type Ia supernovae (SNe Ia) is sensitive to the clustering of matter and provides an independent cosmological probe complementary to SN Ia distance measurements. The Nancy Grace Roman Space Telescope is uniquely sensitive to this measurement as it can discover high redshift SNe Ia and measure them with high precision. We present a methodology for reconstructing the probability distribution of the weak lensing magnification μ of SNe Ia, p(μ), from observational data, and using it to constrain cosmological parameters. We find that the reconstructed p(μ) can be fitted accurately by a stretched Gaussian distribution and used to measure the variance of μ, ξ μ , which can be compared to theoretical predictions in a likelihood analysis. Applying our methodology to a set of realistically simulated SNe Ia expected from the Roman Space Telescope, we find that using the weak lensing magnification of the SNe Ia constrains a combination of matter density Ω m and matter clustering amplitude σ 8 . SN Ia distances alone lead to a better than 1% measurement of Ωm. The combination of SN Ia weak lensing magnification and distance measurements result in a ~10 % measurement on σ 8 . The SNe Ia from Roman will be powerful in constraining the cosmological model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

High galactic latitude cataclysmic variables

A sample of stars which are either known or are good candidates for being halo cataclysmic variables has been isolated. This sample consists of 84 stars of which 57 are dwarf novae (two of these are probable DQ Her systems), 12 are novalike stars, seven are AM Her objects, and eight are classical novae. Periodicities are known for 27 of the systems. If these stars are indeed halo members, then their distances range from 350 pc to 8 kpc and they represent objects that should be significantly older than their disk counterparts. If they do not represent true halo members, then they are intrinsically very faint systems (Mv greater than +10). Comparing the gross properties of this sample to those for the disk cataclysmic variables, it is found that below the period gap, the mean outburst amplitude for the halo dwarf novae is 3 mag greater than for the disk dwarf novae. For systems with periods less than 4 h, the halo cataclysmic variables have a higher percentage of systems below the period gap than the disk cataclysmic variables. The implications of these findings are discussed, as well as possible sample biases and selection effects.

Howell, Steve B.↗

NASA NTRS2010

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

A methodology to compute probabilistically-combined thermo-mechanical fatigue life of polymer matrix laminated composites has been developed and is demonstrated. Matrix degradation effects caused by long-term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress-dependent multifactor-interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability-integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/-45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical-cyclic loads and low thermal-cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical-cyclic loads and high thermal-cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

NASA NTRS2011

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

Chamis, Christos C.↗

DOE OSTI2021

Temporal coherent control of resonant two-photon double ionization of the hydrogen molecule via doubly excited states

Here, we use time-delayed, counter-rotating, circularly polarized few-cycle attosecond nonoverlapping pulses to study the temporal coherent control of the resonant process of two-photon double ionization (TPDI) of hydrogen molecule via doubly excited states for pulse propagation direction along $\hat k$ either parallel or perpendicular to the molecular axis $\hat R$. For $\hat k$∥ $\hat R$ and a pulse carrier frequency of 36 eV resonantly populating the Q 2 1 Π$^+_u$ (1) doubly excited state as well as other 1 Π$^+_u$ doubly excited states, we find that the indirect ionization pathway through these doubly excited states changes the character of the kinematical vortex-shaped momentum distribution produced by the two direct ionization pathways from fourfold to twofold rotational symmetry. This result is similar to what found in TPDI of the He atom involving 1 P$^ o_ {±1}$ doubly excited states; however, angular distributions exhibiting a quantum beat effect between the ground state and a doubly excited state seen for the He atom are observed here for its molecular counterpart with an anomaly in shape and magnitude, not in frequency. The sixfold differential probability integrated over the azimuthal angle of the photoelectron pair shows that this anomaly is due to autoionization decays and quantum beats between doubly excited states. For $\hat k$⟂$\hat R$ and a broadband pulse carrier frequency of 30 eV populating the Q 1 1 Π$^+_u$(1), Q 1 1 Σ$^+_u$(1), Q 2 1 Π$^+_u$(1) , and Q 1 1 Σ$^+_u$ (2) doubly excited states, the momentum distribution is shown to exhibit dynamical electron vortices with four spiral arms, which originates from the interplay between the 1 Δ$^+_ g$ , 1 Π$^+_ g$, and 1 Σ$^+_ g$ dynamical ionization amplitudes. Our treatment within either the adiabatic-nuclei approximation or fixed-nuclei approximation shows that the latter provides a very good account for this correlated process.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1993

Mapping the stability region of the 3:2 Neptune-Pluto resonance

Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.

Levison, H. F.↗