Source author record

Anthony M. Bloch

Anthony M. Bloch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

On two notions of total positivity for partial flag varieties

Given integers $1 \le k_1 < \cdots < k_l \le n-1$, let $\text{Fl}_{k_1,\dots,k_l;n}$ denote the type $A$ partial flag variety consisting of all chains of subspaces $(V_{k_1}\subset\cdots\subset V_{k_l})$ inside $\mathbb{R}^n$, where each $V_k$ has dimension $k$. Lusztig (1994, 1998) introduced the totally positive part $\text{Fl}_{k_1,\dots,k_l;n}^{>0}$ as the subset of partial flags which can be represented by a totally positive $n\times n$ matrix, and defined the totally nonnegative part $\text{Fl}_{k_1,\dots,k_l;n}^{\ge 0}$ as the closure of $\text{Fl}_{k_1,\dots,k_l;n}^{>0}$. On the other hand, following Postnikov (2007), we define $\text{Fl}_{k_1,\dots,k_l;n}^{Δ>0}$ and $\text{Fl}_{k_1,\dots,k_l;n}^{Δ\ge 0}$ as the subsets of $\text{Fl}_{k_1,\dots,k_l;n}$ where all Plücker coordinates are positive and nonnegative, respectively. It follows from the definitions that Lusztig's total positivity implies Plücker positivity, and it is natural to ask when these two notions of positivity agree. Rietsch (2009) proved that they agree in the case of the Grassmannian $\text{Fl}_{k;n}$, and Chevalier (2011) showed that the two notions are distinct for $\text{Fl}_{1,3;4}$. We show that in general, the two notions agree if and only if $k_1, \dots, k_l$ are consecutive integers. We give an elementary proof of this result (including for the case of Grassmannians) based on classical results in linear algebra and the theory of total positivity. We also show that the cell decomposition of $\text{Fl}_{k_1,\dots,k_l;n}^{\ge 0}$ coincides with its matroid decomposition if and only if $k_1,\dots,k_l$ are consecutive integers, which was previously only known for complete flag varieties, Grassmannians, and $\text{Fl}_{1,3;4}$. Finally, we determine which notions of positivity are compatible with a natural action of the cyclic group of order $n$ that rotates the index set.

preprint2020arXiv

Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close Orbits

Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$.

preprint2020arXiv

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system.

preprint2016arXiv

Flag-based Control of Orbit Dynamics in Quantum Lindblad Systems

In this paper, we demonstrate that the dynamics of an $n$-dimensional Lindblad control system can be separated into its inter- and intra-orbit dynamics when there is fast controllability. This can be viewed as a control system on the simplex of density operator spectra, where the flag representing the eigenspaces is viewed as a control variable. The local controllability properties of this control system can be analyzed when the control-set of flags is limited to a finite subset. In particular, there is a natural finite subset of $n!$ flags that are effective for low-purity orbits.

preprint2016arXiv

Flag-Based Control of Quantum Purity for $n=2$ Systems

This paper investigates the fast Hamiltonian control of $n=2$ density operators by continuously varying the flag as one moves away from the completely mixed state. In general, the critical points and zeros of the purity derivative can only be solved analytically in the limit of minimal purity. We derive differential equations that maintain these features as the purity increases. In particular, there is a thread of points in the Bloch ball that locally maximizes the purity derivative, and a corresponding thread that minimizes it. Additionally, we show there is a closed surface of points inside of which the purity derivative is positive, and inside of which is negative. We argue that this approach may be useful in studying higher-dimensional systems.

preprint2016arXiv

The Stability of Tidal Equilibrium for Hierarchical Star-Planet-Moon Systems

Motivated by the current search for exomoons, this paper considers the stability of tidal equilibrium for hierarchical three-body systems containing a star, a planet, and a moon. In this treatment, the energy and angular momentum budgets include contributions from the planetary orbit, lunar orbit, stellar spin, planetary spin, and lunar spin. The goal is to determine the optimized energy state of the system subject to the constraint of constant angular momentum. Due to the lack of a closed form solution for the full three-body problem, however, we must use use an approximate description of the orbits. We first consider the Keplerian limit and find that the critical energy states are saddle points, rather than minima, so that these hierarchical systems have no stable tidal equilibrium states. We then generalize the calculation so that the lunar orbit is described by a time-averaged version of the circular restricted three-body problem. In this latter case, the critical energy state is a shallow minimum, so that a tidal equilibrium state exists. In both cases, however, the lunar orbit for the critical point lies outside the boundary (roughly half the Hill radius) where (previous) numerical simulations indicate dynamical instability. These results suggest that star-planet-moon systems have no viable long-term stable states analogous to those found for two-body systems.

preprint2015arXiv

Planets in Other Universes: Habitability constraints on density fluctuations and galactic structure

Motivated by the possibility that different versions of the laws of physics could be realized within other universes, this paper delineates the galactic parameters that allow for habitable planets and revisits constraints on the amplitude $Q$ of the primordial density fluctuations. Previous work indicates that large values of $Q$ lead to galaxies so dense that planetary orbits cannot survive long enough for life to develop. Small values of $Q$ lead to delayed star formation, loosely bound galaxies, and compromised heavy element retention. This work generalizes previous treatments: [A] We consider models for the internal structure of galaxies and find the fraction of galactic real estate that allows stable, long-lived planetary orbits. [B] We perform a large ensemble of numerical simulations to estimate cross sections for the disruption of planetary orbits due to interactions with passing stars. [C] We consider disruption due to the background radiation fields produced by the galaxies. [D] One consequence of intense galactic background radiation fields is that some portion of the galaxy, denoted as the Galactic Habitable Zone, will provide the right flux levels to support habitable planets for essentially any planetary orbit. As $Q$ increases, the fraction of stars in a galaxy that allow for habitable planets decreases due to both orbital disruption and the intense background radiation. However, the outer parts of the galaxy always allow for habitable planets, so that the value of $Q$ does not have a well-defined upper limit. Moreover, some Galactic Habitable Zones are large enough to support more potentially habitable planets than the galaxies found in our universe. These results suggest that the possibilities for habitability in other universes are somewhat more favorable and far more diverse than previously imagined.

preprint2014arXiv

On the Stability of Extrasolar Planetary Systems and other Closely Orbiting Pairs

This paper considers the stability of tidal equilibria for planetary systems in which stellar rotation provides a significant contribution to the angular momentum budget. We begin by applying classic stability considerations for two bodies to planetary systems --- where one mass is much smaller than the other. The application of these stability criteria to a subset of the Kepler sample indicates that the majority of the systems are not in a stable equilibrium state. Motivated by this finding, we generalize the stability calculation to include the quadrupole moment for the host star. In general, a stable equilibrium requires that the total system angular momentum exceeds a minimum value (denoted here as $L_X$) and that the orbital angular momentum of the planet exceeds a minimum fraction of the total. Most, but not all, of the observed planetary systems in the sample have enough total angular momentum to allow an equilibrium state. Even with the generalizations of this paper, however, most systems have too little orbital angular momentum (relative to the total) and are not in an equilibrium configuration. Finally, we consider the time evolution of these planetary systems; the results constrain the tidal quality factor of the stars and suggest that $10^6\le{Q_\ast}\le10^7$.

preprint2013arXiv

Evolution of Planetary Orbits with Stellar Mass Loss and Tidal Dissipation

Intermediate mass stars and stellar remnants often host planets, and these dynamical systems evolve because of mass loss and tides. This paper considers the combined action of stellar mass loss and tidal dissipation on planetary orbits in order to determine the conditions required for planetary survival. Stellar mass loss is included using a so-called Jeans model, described by a dimensionless mass loss rate γand an index β. We use an analogous prescription to model tidal effects, described here by a dimensionless dissipation rate Γand two indices (q,p). The initial conditions are determined by the starting value of angular momentum parameter η(equivalently, the initial eccentricity) and the phase θof the orbit. Within the context of this model, we derive an analytic formula for the critical dissipation rate Γ, which marks the boundary between orbits that spiral outward due to stellar mass loss and those that spiral inward due to tidal dissipation. This analytic result Γ=Γ(γ,β,q,p,η,θ) is essentially exact for initially circular orbits and holds to within an accuracy of 50% over the entire multi-dimensional parameter space, where the individual parameters vary by several orders of magnitude. For stars that experience mass loss, the stellar radius often displays quasi-periodic variations, which produce corresponding variations in tidal forcing; we generalize the calculation to include such pulsations using a semi-analytic treatment that holds to the same accuracy as the non-pulsating case. These results can be used in many applications, e.g., to predict/constrain properties of planetary systems orbiting white dwarfs.

preprint2013arXiv

Evolution of Planetary Systems with Time Dependent Stellar Mass Loss

Observations indicate that intermediate mass stars, binary stars, and stellar remnants often host planets; a complete explanation of these systems requires an understanding of how planetary orbits evolve as their central stars lose mass. Motivated by these dynamical systems, this paper generalizes in two directions previous studies of orbital evolution in planetary systems with stellar mass loss: [1] Many previous treatments focus on constant mass loss rates and much of this work is carried out numerically. Here we study a class of single planet systems where the stellar mass loss rate is time dependent. The mass loss rate can be increasing or decreasing, but the stellar mass always decreases monotonically. For this class of models, we develop analytic approximations to specify the final orbital elements for planets that remain bound after the epoch of mass loss, and find the conditions required for the planets to become unbound. We also show that for some mass loss functions, planets become unbound only in the asymptotic limit where the stellar mass vanishes. [2] We consider the chaotic evolution for two planet systems with stellar mass loss. Here we focus on a model consisting of analogs of Jupiter, Saturn, and the Sun. By monitoring the divergence of initially similar trajectories through time, we calculate the Lyapunov exponents of the system. This analog solar system is chaotic in the absence of mass loss with Lyapunov time in the range 5 - 10 Myr; we find that the Lyapunov time decreases with increasing stellar mass loss rate, with a nearly linear relationship between the two time scales. Taken together, the results of this paper help provide an explanation for a wide range of dynamical evolution that occurs in solar systems with stellar mass loss.

preprint2013arXiv

Hill's Equation with Small Fluctuations: Cycle to Cycle Variations and Stochastic Processes

Hill's equations arise in a wide variety of physical problems, and are specified by a natural frequency, a periodic forcing function, and a forcing strength parameter. This classic problem is generalized here in two ways: [A] to Random Hill's equations which allow the forcing strength q_k, the oscillation frequency λ_k, and the period τ_k of the forcing function to vary from cycle to cycle, and [B] to Stochastic Hill's equations which contain (at least) one additional term that is a stochastic process ξ. This paper considers both random and stochastic Hill's equations with small parameter variations, so that p_k=q_k-<q_k>, \ell_k=λ_k-<λ_k>, and ξare all O(ε), where ε<<1. We show that random Hill's equations and stochastic Hill's equations have the same growth rates when the parameter variations p_k and \ell_k obey certain constraints given in terms of the moments of ξ. For random Hill's equations, the growth rates for the solutions are given by the growth rates of a matrix transformation, under matrix multiplication, where the matrix elements vary from cycle to cycle. Unlike classic Hill's equations where the parameter space (the λ-q plane) displays bands of stable solutions interlaced with bands of unstable solutions, random Hill's equations are generically unstable. We find analytic approximations for the growth rates of the instability; for the regime where Hill's equation is classically stable, and the parameter variations are small, the growth rate γ= O(ε^2). Using the relationship between the (\ell_k,p_k) and the ξ, this result for γcan be used to find growth rates for stochastic Hill's equations.

preprint2012arXiv

Gradient flows in the normal and Kähler metrics and triple bracket generated metriplectic systems

The dynamics of gradient and Hamiltonian flows with particular application to flows on adjoint orbits of a Lie group and the extension of this setting to flows on a loop group are discussed. Different types of gradient flows that arise from different metrics including the so-called normal metric on adjoint orbits of a Lie group and the Kähler metric are compared. It is discussed how a Kähler metric can arise from a complex structure induced by the Hilbert transform. Hybrid and metriplectic flows that arise when one has both Hamiltonian and gradient components are examined. A class of metriplectic systems that is generated by completely antisymmetric triple brackets is described and for finite-dimensional systems given a Lie algebraic interpretation. A variety of explicit examples of the several types of flows are given.

preprint2012arXiv

Hamel's Formalism and Variational Integrators on a Sphere

This paper discusses Hamel's formalism and its applications to structure-preserving integration of mechanical systems. It utilizes redundant coordinates in order to eliminate multiple charts on the configuration space as well as nonphysical artificial singularities induced by local coordinates, while keeping the minimal possible degree of redundancy and avoiding integration of differential-algebraic equations.

preprint2012arXiv

Mean Motion Resonances in Exoplanet Systems: An Investigation Into Nodding Behavior

Motivated by the large number of extrasolar planetary systems that are near mean motion resonances, this paper explores a related type of dynamical behavior known as "nodding". Here, the resonance angle of a planetary system executes libration (oscillatory motion) for several cycles, circulates for one or more cycles, and then enters once again into libration. This type of complicated dynamics can affect our interpretation of observed planetary systems that are in or near mean motion resonance. This work shows that planetary systems in (near) mean motion resonance can exhibit nodding behavior, and outlines the portion of parameter space where it occurs. This problem is addressed using both full numerical integrations of the planetary systems and via model equations obtained through expansions of the disturbing function. In the latter approach, we identify the relevant terms that allow for nodding. The two approaches are in agreement, and show that nodding often occurs when a small body is in an external mean motion resonance with a larger planet. As a result, the nodding phenomenon can be important for interpreting observations of transit timing variations, where the existence of smaller bodies is inferred through their effects on larger, observed transiting planets. For example, in actively nodding planetary systems, both the amplitude and frequency of the transit timing variations depend on the observational time window.

preprint2011arXiv

Accretion of Rocky Planets by Hot Jupiters

The observed population of Hot Jupiters displays a stunning variety of physical properties, including a wide range of densities and core sizes for a given planetary mass. Motivated by the observational sample, this paper studies the accretion of rocky planets by Hot Jupiters, after the Jovian planets have finished their principal migration epoch and become parked in $\sim4$-day orbits. In this scenario, rocky planets form later and then migrate inward due to torques from the remaining circumstellar disk, which also damps the orbital eccentricity. This mechanism thus represents one possible channel for increasing the core masses and metallicities of Hot Jupiters. This paper determines probabilities for the possible end states for the rocky planet: collisions with the Jovian planets, accretion onto the star, ejection from the system, and long-term survival of both planets. These probabilities depend on the mass of the Jovian planet and its starting orbital eccentricity, as well as the eccentricity damping rate for the rocky planet. Since these systems are highly chaotic, a large ensemble ($N\sim10^3$) of simulations with effectively equivalent starting conditions is required. Planetary collisions are common when the eccentricity damping rate is sufficiently low, but are rare otherwise. For systems that experience planetary collisions, this work determines the distributions of impact velocities -- both speeds and impact parameters -- for the collisions. These velocity distributions help determine the consequences of the impacts, e.g., where energy and heavy elements are deposited within the giant planets.

preprint2011arXiv

Discrete Hamilton-Jacobi Theory

We develop a discrete analogue of Hamilton-Jacobi theory in the framework of discrete Hamiltonian mechanics. The resulting discrete Hamilton-Jacobi equation is discrete only in time. We describe a discrete analogue of Jacobi's solution and also prove a discrete version of the geometric Hamilton-Jacobi theorem. The theory applied to discrete linear Hamiltonian systems yields the discrete Riccati equation as a special case of the discrete Hamilton-Jacobi equation. We also apply the theory to discrete optimal control problems, and recover some well-known results, such as the Bellman equation (discrete-time HJB equation) of dynamic programming and its relation to the costate variable in the Pontryagin maximum principle. This relationship between the discrete Hamilton-Jacobi equation and Bellman equation is exploited to derive a generalized form of the Bellman equation that has controls at internal stages.

preprint2011arXiv

Nonholonomic Hamilton-Jacobi Theory via Chaplygin Hamiltonization

We develop Hamilton-Jacobi theory for Chaplygin systems, a certain class of nonholonomic mechanical systems with symmetries, using a technique called Hamiltonization, which transforms nonholonomic systems into Hamiltonian systems. We give a geometric account of the Hamiltonization, identify necessary and sufficient conditions for Hamiltonization, and apply the conventional Hamilton-Jacobi theory to the Hamiltonized systems. We show, under a certain sufficient condition for Hamiltonization, that the solutions to the Hamilton-Jacobi equation associated with the Hamiltonized system also solve the nonholonomic Hamilton-Jacobi equation associated with the original Chaplygin system. The results are illustrated through several examples.

preprint2010arXiv

Effects of Turbulence, Eccentricity Damping, and Migration Rate on the Capture of Planets into Mean Motion Resonance

Pairs of migrating extrasolar planets often lock into mean motion resonance as they drift inward. This paper studies the convergent migration of giant planets (driven by a circumstellar disk) and determines the probability that they are captured into mean motion resonance. The probability that such planets enter resonance depends on the type of resonance, the migration rate, the eccentricity damping rate, and the amplitude of the turbulent fluctuations. This problem is studied both through direct integrations of the full 3-body problem, and via semi-analytic model equations. In general, the probability of resonance decreases with increasing migration rate, and with increasing levels of turbulence, but increases with eccentricity damping. Previous work has shown that the distributions of orbital elements (eccentricity and semimajor axis) for observed extrasolar planets can be reproduced by migration models with multiple planets. However, these results depend on resonance locking, and this study shows that entry into -- and maintenance of -- mean motion resonance depends sensitively on migration rate, eccentricity damping, and turbulence.

preprint2010arXiv

Hill's Equation with Random Forcing Parameters: Determination of Growth Rates through Random Matrices

This paper derives expressions for the growth rates for the random 2 x 2 matrices that result from solutions to the random Hill's equation. The parameters that appear in Hill's equation include the forcing strength and oscillation frequency. The development of the solutions to this periodic differential equation can be described by a discrete map, where the matrix elements are given by the principal solutions for each cycle. Variations in the forcing strength and oscillation frequency lead to matrix elements that vary from cycle to cycle. This paper presents an analysis of the growth rates including cases where all of the cycles are highly unstable, where some cycles are near the stability border, and where the map would be stable in the absence of fluctuations. For all of these regimes, we provide expressions for the growth rates of the matrices that describe the solutions.

preprint2009arXiv

General Analysis of Type I Planetary Migration with Stochastic Perturbations

This paper presents a generalized treatment of Type I planetary migration in the presence of stochastic perturbations. In many planet-forming disks, the Type I migration mechanism, driven by asymmetric torques, acts on a short time scale and compromises planet formation. If the disk also supports MHD instabilities, however, the corresponding turbulent fluctuations produce additional stochastic torques that modify the steady inward migration scenario. This work studies the migration of planetary cores in the presence of stochastic fluctuations using complementary methods, including a Fokker-Planck approach and iterative maps. Stochastic torques have two main effects: [1] Through outward diffusion, a small fraction of the planetary cores can survive in the face of Type I inward migration. [2] For a given starting condition, the result of any particular realization of migration is uncertain, so that results must be described in terms of the distributions of outcomes. In addition to exploring different regimes of parameter space, this paper considers the effects of the outer disk boundary condition, varying initial conditions, and time-dependence of the torque parameters. For disks with finite radii, the fraction of surviving planets decreases exponentially with time. We find the survival fractions and decay rates for a range of disk models, and find the expected distribution of locations for surviving planets. For expected disk properties, the survival fraction lies in the range $0.01 < p_S < 0.1$.

preprint2009arXiv

Hill's Equation with Random Forcing Parameters: The Limit of Delta Function Barriers

This paper considers random Hill's equations in the limit where the periodic forcing function becomes a Dirac delta function. For this class of equations, the forcing strength $q_k$, the oscillation frequency $\af_k$, and the period are allowed to vary from cycle to cycle. Such equations arise in astrophysical orbital problems in extended mass distributions, in the reheating problem for inflationary cosmologies, and in periodic Schr{ö}dinger equations. The growth rates for solutions to the periodic differential equation can be described by a matrix transformation, where the matrix elements vary from cycle to cycle. Working in the delta function limit, this paper addresses several coupled issues: We find the growth rates for the $2 \times 2$ matrices that describe the solutions. This analysis is carried out in the limiting regimes of both large $q_k \gg 1$ and small $q_k \ll 1$ forcing strength parameters. For the latter case, we present an alternate treatment of the dynamics in terms of a Fokker-Planck equation, which allows for a comparison of the two approaches. Finally, we elucidate the relationship between the fundamental parameters $(\af_k,q_k)$ appearing in the stochastic differential equation and the matrix elements that specify the corresponding discrete map. This work provides analytic -- and accurate -- expressions for the growth rates of these stochastic differential equations in both the $q_k \gg1 $ and the $q_k \ll 1$ limits.

preprint2009arXiv

Nonholonomic Hamilton-Jacobi equation and Integrability

We discuss an extension of the Hamilton-Jacobi theory to nonholonomic mechanics with a particular interest in its application to exactly integrating the equations of motion. We give an intrinsic proof of a nonholonomic analogue of the Hamilton--Jacobi theorem. Our intrinsic proof clarifies the difference from the conventional Hamilton-Jacobi theory for unconstrained systems. The proof also helps us identify a geometric meaning of the conditions on the solutions of the Hamilton-Jacobi equation that arise from nonholonomic constraints. The major advantage of our result is that it provides us with a method of integrating the equations of motion just as the unconstrained Hamilton--Jacobi theory does. In particular, we build on the work by Iglesias-Ponte, de Leon, and Martin de Diego so that the conventional method of separation of variables applies to some nonholonomic mechanical systems. We also show a way to apply our result to systems to which separation of variables does not apply.

preprint2008arXiv

Mean Motion Resonances in Extrasolar Planetary Systems with Turbulence, Interactions, and Damping

This paper continues previous work on the effects of turbulence on mean motion resonances in extrasolar planetary systems. Turbulence is expected to arise in the disks that form planets, and these fluctuations act to compromise resonant configurations. This paper extends previous work by considering how interactions between the planets and possible damping effects imposed by the disk affect the outcomes. These physical processes are studied using three approaches: numerical integrations of the 3-body problem with additional forcing due to turbulence, model equations that reduce the problem to stochastically driven oscillators, and Fokker-Planck equations that describe the time evolution of an ensemble of systems. With this combined approach, we elucidate the physics of how turbulence can remove extrasolar planetary systems from mean motion resonance. As expected, systems with sufficiently large damping (dissipation) can maintain resonance, in spite of turbulent forcing. In the absence of strong damping, ensembles of these systems exhibit two regimes of behavior, where the fraction of the bound states decreases as a power-law or as an exponential. Both types of behavior can be understood through the model developed herein. For systems with weak interactions between planets, the model reduces to a stochastic pendulum, and the fraction of bound states decreases as a power-law. For highly interactive systems, the dynamics are more complicated and the fraction of bound states decreases exponentially. We show how planetary interactions lead to drift terms in the Fokker-Planck equation and account for this exponential behavior. In addition to clarifying the physical processes involved, this paper strengthens the finding that turbulence implies that mean motions resonances should be rare.

preprint2006arXiv

A Discrete Variational Integrator for Optimal Control Problems on SO(3)

In this paper we study a discrete variational optimal control problem for the rigid body. The cost to be minimized is the external torque applied to move the rigid body from an initial condition to a pre-specified terminal condition. Instead of discretizing the equations of motion, we use the discrete equations obtained from the discrete Lagrange--d'Alembert principle, a process that better approximates the equations of motion. Within the discrete-time setting, these two approaches are not equivalent in general. The kinematics are discretized using a natural Lie-algebraic formulation that guarantees that the flow remains on the Lie group SO(3) and its algebra so(3). We use Lagrange's method for constrained problems in the calculus of variations to derive the discrete-time necessary conditions. We give a numerical example for a three-dimensional rigid body maneuver.

preprint2006arXiv

Controlled Lagrangians and Potential Shaping for Stabilization of Discrete Mechanical Systems

The method of controlled Lagrangians for discrete mechanical systems is extended to include potential shaping in order to achieve complete state-space asymptotic stabilization. New terms in the controlled shape equation that are necessary for matching in the discrete context are introduced. The theory is illustrated with the problem of stabilization of the cart-pendulum system on an incline. We also discuss digital and model predictive control.