Engineering topics

Malhotra, Renu

Publications and source records attributed to Malhotra, Renu.

At least 19 records

DOE OSTI2020

Dynamical Instabilities in Systems of Multiple Short-period Planets Are Likely Driven by Secular Chaos: A Case Study of Kepler-102

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in ∼20% of the cases on a wide range of timescales (up to 5 × 10{sup 9} orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multiplanet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean-motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet–planet collisions in ∼80% of Kepler-102 simulations. Kepler-102 b likely needs a mass ≳0.1 M {sub ⊕}, hence a bulk density exceeding about half Earth’s, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system’s AMD evolution with a “spectral fraction” calculated from the power spectrum of short integrations (∼5 × 10{sup 6} orbits). We find that small spectral fractions (≲0.01) are strongly associated with dynamical stability on long timescales (5 × 10{sup 9} orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many nonresonant multiplanet systems and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long-term (in)stability of multiplanet systems from short simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Extreme Debris Disk Variability: Exploring the Diverse Outcomes of Large Asteroid Impacts During the Era of Terrestrial Planet Formation

The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100–200 Myr, when violent impacts occur between large planetesimals of sizes up to protoplanets. Such events are marked by the production of large amounts of debris, as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-sized bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system suggest two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

THE MASS DISTRIBUTION FUNCTION OF PLANETS

The distribution of orbital period ratios of adjacent planets in extrasolar planetary systems discovered by the Kepler space telescope exhibits a peak near ∼1.5–2, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below ∼1.5. We find from these data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill’s criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed 10{sup −3.2}m{sub ∗}, or about two-thirds of Jupiter’s mass for solar-mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill’s criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of logm/M{sub ⊕} of ∼(0.6−1.0) and ∼(1.1−1.2), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

DEBRIS DISTRIBUTION IN HD 95086—A YOUNG ANALOG OF HR 8799

HD 95086 is a young early-type star that hosts (1) a 5 M{sub J} planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 μm crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 μm. Due to the low resolution of the MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing ∼5% of the total dust mass. We also present detailed analysis of the disk spectral energy distribution and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (∼175 K) belt, a cold (∼55 K) disk, and an extended disk halo (up to ∼800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet, and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass ∼ 5 M{sub J} , could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

PLANETARY CHAOTIC ZONE CLEARING: DESTINATIONS AND TIMESCALES

We investigate the orbital evolution of particles in a planet's chaotic zone to determine their final destinations and their timescales of clearing. There are four possible final states of chaotic particles: collision with the planet, collision with the star, escape, or bounded but non-collision orbits. In our investigations, within the framework of the planar circular restricted three body problem for planet-star mass ratio μ in the range 10{sup –9} to 10{sup –1.5}, we find no particles hitting the star. The relative frequencies of escape and collision with the planet are not scale-free, as they depend upon the size of the planet. For planet radius R{sub p} ≥ 0.001 R{sub H} where R{sub H} is the planet's Hill radius, we find that most chaotic zone particles collide with the planet for μ ≲ 10{sup –5}; particle scattering to large distances is significant only for higher mass planets. For fixed ratio R{sub p} /R{sub H} , the particle clearing timescale, T {sub cl}, has a broken power-law dependence on μ. A shallower power law, T {sub cl} ∼ μ{sup –1/3}, prevails at small μ where particles are cleared primarily by collisions with the planet; a steeper power law, T {sub cl} ∼ μ{sup –3/2}, prevails at larger μ where scattering dominates the particle loss. In the limit of vanishing planet radius, we find T {sub cl} ≈ 0.024 μ{sup –3/2}. The interior and exterior boundaries of the annular zone in which chaotic particles are cleared are increasingly asymmetric about the planet's orbit for larger planet masses; the inner boundary coincides well with the classical first order resonance overlap zone, Δa {sub cl,} {sub int} ≅ 1.2 μ{sup 0.28} a{sub p} ; the outer boundary is better described by Δa {sub cl,} {sub ext} ≅ 1.7 μ{sup 0.31} a{sub p} , where a{sub p} is the planet-star separation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Predictions for shepherding planets in scattered light images of debris disks

Planets can affect debris disk structure by creating gaps, sharp edges, warps, and other potentially observable signatures. However, there is currently no simple way for observers to deduce a disk-shepherding planet's properties from the observed features of the disk. Here we present a single equation that relates a shepherding planet's maximum mass to the debris ring's observed width in scattered light, along with a procedure to estimate the planet's eccentricity and minimum semimajor axis. We accomplish this by performing dynamical N-body simulations of model systems containing a star, a single planet, and an exterior disk of parent bodies and dust grains to determine the resulting debris disk properties over a wide range of input parameters. We find that the relationship between planet mass and debris disk width is linear, with increasing planet mass producing broader debris rings. We apply our methods to five imaged debris rings to constrain the putative planet masses and orbits in each system. Observers can use our empirically derived equation as a guide for future direct imaging searches for planets in debris disk systems. In the fortuitous case of an imaged planet orbiting interior to an imaged disk, the planet's maximum mass can be estimated independent of atmospheric models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

PLANETS NEAR MEAN-MOTION RESONANCES

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases, the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the three-body problem. We show that for the 3 : 2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100 M{sub Circled-Plus }. This predicted mass range is higher-by at least a factor of three-than the range expected from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

ASTEROID BELTS IN DEBRIS DISK TWINS: VEGA AND FOMALHAUT

Vega and Fomalhaut are similar in terms of mass, ages, and global debris disk properties; therefore, they are often referred to as 'debris disk twins'. We present Spitzer 10-35 {mu}m spectroscopic data centered at both stars and identify warm, unresolved excess emission in the close vicinity of Vega for the first time. The properties of the warm excess in Vega are further characterized with ancillary photometry in the mid-infrared and resolved images in the far-infrared and submillimeter wavelengths. The Vega warm excess shares many similar properties with the one found around Fomalhaut. The emission shortward of {approx}30 {mu}m from both warm components is well described as a blackbody emission of {approx}170 K. Interestingly, two other systems, {epsilon} Eri and HR 8799, also show such an unresolved warm dust using the same approach. These warm components may be analogous to the solar system's zodiacal dust cloud, but of far greater mass (fractional luminosity of {approx}10{sup -5} to 10{sup -6} compared to 10{sup -8} to 10{sup -7}). The dust temperature and tentative detections in the submillimeter suggest that the warm excess arises from dust associated with a planetesimal ring located near the water-frost line and presumably created by processes occurring at similar locations in other debris systems as well. We also review the properties of the 2 {mu}m hot excess around Vega and Fomalhaut, showing that the dust responsible for the hot excess is not spatially associated with the dust we detected in the warm belt. We suggest it may arise from hot nano grains trapped in the magnetic field of the star. Finally, the separation between the warm and cold belt is rather large with an orbital ratio {approx}>10 in all four systems. In light of the current upper limits on the masses of planetary objects and the large gap, we discuss the possible implications for their underlying planetary architecture and suggest that multiple, low-mass planets likely reside between the two belts in Vega and Fomalhaut.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2000

The nu Andromedae System: Models and Stability

Radial velocity observations of the F8 V star nu Andromedae taken at Lick and at Whipple Observatories have revealed evidence of three periodicities in the line-of-sight velocity of the star. These periodicities have been interpreted as evidence for at least three low-mass companions (LMCs) revolving around nu Andromedae. The mass and orbital parameters inferred for these companions raise questions about the dynamical stability of the system. We report here results from our independent analysis of the published radial velocity data, as well as new unpublished data taken at Lick Observatory. Our results confirm the finding of three periods in the data. Our best fits to the data, on the assumption that these periods arise from the gravitational perturbations of companions in Keplerian orbits, are also generally in agreement but with some differences from the earlier findings. We find that the available data do not constrain well the orbital eccentricity of the middle companion in a three-companion model of the data. We also find that in order for our best-fit model to the Lick data to be dynamically stable over the lifetime of the star (approximately 2 billion years), the system must have a mean inclination to the plane of the sky greater than 13 deg. The corresponding minimum inclination for the best fit to the Whipple data set is 19 deg. These values imply that the maximum mass for the outer companion can be no greater than about 20 Jupiter masses. Our analysis of the stability of the putative systems also places constraints on the relative inclinations of the orbital planes of the companions. We comment on global versus local (i.e., method of steepest descent) means of finding best-fit orbits from radial velocity data sets.

Stepinski, Tomasz F.↗

NASA NTRS1997

Depletion of the Outer Asteroid Belt

During the early history of the solar system, it is likely that the outer planets changed their distance from the sun, and hence, their influence on the asteroid belt evolved with time. The gravitational influence of Jupiter and Saturn on the orbital evolution of asteroids in the outer asteroid belt was calculated. The results show that the sweeping of mean motion resonances associated with planetary migration efficiently destabilizes orbits in the outer asteroid belt on a time scale of 10 million years. This mechanism provides an explanation for the observed depletion of asteroids in that region.

Liou, Jer-Chyi↗

NASA NTRS1996

The Phase Space Structure Near Neptune Resonances in the Kuiper Belt

The Solar system beyond Neptune is believed to house a population of small primordial bodies left over from the planet formation process. The region up to heliocentric distance -50 AU (a.k.a. the Kuiper Belt) may be the source of the observed short-period comets. In this region, the phase space structure near orbital resonances with Neptune is of special interest for the long-term stability of orbits. There is reason to believe that a significant fraction (perhaps most) of the Kuiper Belt objects reside preferentially in these resonance locations. This paper describes the dynamics of small objects near the major orbital resonances with Neptune. Estimates of the widths of stable resonance zones as well as the properties of resonant orbits are obtained from the circular, planar restricted three-body model. Although this model does not contain the full complexity of the long-term orbital dynamics of Kuiper Belt objects subject to the full N-body perturbations of all the planets, it does provide a baseline for the phase space structure and properties of resonant orbits in the trans-Neptunian Solar system.

Malhotra, Renu↗

NASA NTRS1995

The Origin of Pluto's Orbit: Implications for the Solar System Beyond Neptune

The origin of the highly eccentric, inclined, and resonance-locked orbit of Pluto has long been a puzzle. A possible explanation has been proposed recently which suggests that these extraordinary orbital properties may be a natural consequence of the formation and early dynamical evolution of the outer solar system. A resonance capture mechanism is possible during the clearing of the residual planetesimal debris and the formation of the Oort Cloud of comets by planetesimal mass loss from the vicinity of the giant planets. If this mechanism were in operation during the early history of the planetary system, the entire region between the orbit of Neptune and approximately 50 AU would have been swept by first-order mean motion resonances. Thus, resonance capture could occur not only for Pluto, but quite generally for other trans-Neptunian small bodies. Some consequences of this evolution for the present-day dynamical structure of the trans-Neptunian region are (1) most of the objects in the region beyond Neptune and up to approximately 50 AU exist in very narrow zones located at orbital resonances with Neptune (particularly the 3:2 and the 2:1 resonances); and (2) these resonant objects would have significantly large eccentricities. The distribution of objects in the Kuiper Belt as predicted by this theory is presented here.

Malhotra, Renu↗

NASA NTRS1994

A mapping method for the gravitational few-body problem with dissipation

Recently a new class of numerical integration methods - 'mixed variable symplectic integrators' - has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ordinary differential equations (ODE) integration methods. Here we present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.

Malhotra, Renu↗

NASA NTRS1994

Nonlinear Resonances in the Solar System

Orbital resonances are ubiquitous in the Solar system. They play a decisive role in the long term dynamics, and in some cases the physical evolution, of the planets and of their natural satellites, as well as the evolution of small bodies (including dust) in the planetary system. The few-body gravitational problem of hierarchical planetary-type systems allows for a complex range of dynamical timescales, from the fast orbital periods to the very slow orbit precession rates. The interaction of fast and slow degrees of freedom produces a rich diversity of resonance phenomena. Weak dissipative effects - such as tides or radiation drag forces - also produce unexpectedly rich dynamical behaviors. This paper provides a mostly qualitative discussion of simple dynamical models for the commonly encountered orbital resonance phenomena in the Solar system.

Malhotra, Renu↗

NASA NTRS1994

Pluto's Heliocentric Orbit

We review the current state of knowledge regarding Pluto's orbital motion.

Pluto planetary orbit orbit charon↗

NASA NTRS1993

Orbital Resonances in the Solar Nebula: Strengths and Weaknesses

A planetesimal moving in the Solar Nebula experiences an aerodynamic drag which causes its orbit to circularize and shrink. However, resonant perturbations from a protoplanet interior to the planetesimal's orbit can counteract both the orbital decay and the damping of the eccentricity: the planetesimal can be captured into an orbital resonance and its eccentricity pumped up to a modestly high equilibrium value. Thus, orbital resonances constitute (partial) barriers to the delivery of planetesimals into the feeding zone of the protoplanet. We have established the characteristics of the phenomenon of resonance capture by gas drag in the circular restricted three-body approximation. We have determined the strengths of the equilibrium resonant orbits with respect to impulsive velocity perturbations. We conclude that planetesimals captured in orbital resonances are quite vulnerable to being dislocated from these orbits by mutual planetesimal interactions, but that the resonances are effective in slowing down the rate of orbital decay of planetesimals. Only very small bodies, less than or equal to 100 m, are able to reach a approx. 1 Earth mass protoplanet without being slowed down by resonances.

Malhotra, Renu↗

NASA NTRS1993

Three-body effects in the PSR 1257+12 planetary system

A detailed theoretical analysis of the three-body effects in the putative planetary system of PSR 1257+12 is presented. Ways in which these effects are manifested in the pattern of pulse arrival times are discussed; the dominant perturbation is described as a modulation of the phases of the near-sinusoidal signals of the two planetary companions. Explicit formulas for the time dependence of the oscillating orbital elements that are needed for an improved timing model for this system are provided. If a timing model with fixed, independent Keplerian orbits continues to be used for the timing analysis, and if two planets are indeed orbiting this pulsar, then the three-body effects should become detectable by means of a growth in the postfit residuals as more observations are accumulated. If the typical error in the pulse arrival time measurements is about 10 microsec, the amplitude of the postfit residuals will increase beyond this level with three to five years of timing observations.

Malhotra, Renu↗

NASA NTRS1993

Orbital dynamics of PSR1257+12 and its two planetary companions

I present a review of the orbital dynamics of the three body problem of PSR1257+12 and its two putative planetary companions. The orbital dynamics provides two 'reality checks' on the planetary interpretation of the pulse arrival time residuals observed by Wolszczan and Frail (1992). First, by examining the dynamical stability of the system for the deduced orbital periods, one can place upper limits on the companion masses of about two times the mass of Jupiter. Since the lower limits of the companion masses deduced from the observations are well below this upper limit, the planetary interpretation passes this first 'reality check'. Second, the 3:2 near-commensurability of the orbital periods effects periodic variations of the Keplerian orbital parameters. The amplitude of these variations provides a direct measure of the masses of the companions relative to the mass of the pulsar. The signature of these variations in pulse time-of-arrival data is not straightforward, but, if the planetary interpretation is correct, it should become detectable as an additional quasi-periodic phase residual over a period of a few years.

Malhotra, Renu↗