BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Margaret Pan

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

ResearcherUnclaimed source record
astro-ph.EPastro-phastro-ph.SR

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

11 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade

Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λ$ = 450 $μ$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μ$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Apocenter glow in eccentric debris disks: implications for Fomalhaut and epsilon Eridani

Debris disks often take the form of eccentric rings with azimuthal asymmetries in surface brightness. Such disks are often described as showing "pericenter glow", an enhancement of the disk brightness in regions nearest the central star. At long wavelengths, however, the disk apocenters should appear brighter than their pericenters: in the long wavelength limit, we find the apocenter/pericenter flux ratio scales as 1+e for disk eccentricity e. We produce new models of this "apocenter glow" to explore its causes and wavelength dependence and study its potential as a probe of dust grain properties. Based on our models, we argue that several far-infrared and (sub)millimeter images of the Fomalhaut and epsilon Eridani debris rings obtained with Herschel, JCMT, SHARC II, ALMA, and ATCA should be reinterpreted as suggestions or examples of apocenter glow. This reinterpretation yields new constraints on the disks' dust grain properties and size distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the mass and origin of Chariklo's rings

Observations in 2013 and 2014 of the Centaur 10199 Chariklo and its ring system consistently indicated that the radial width of the inner, more massive ring varies with longitude. That strongly suggests that this ring has a finite eccentricity despite the fast differential precession that Chariklo's large quadrupole moment should induce. If the inferred apse alignment is maintained by the ring's self-gravity, as it is for the Uranian rings, we estimate a ring mass of a few times 10^16 g and a typical particle size of a few meters. These imply a short collisional spreading time of ~10^5 years, somewhat shorter than the typical Centaur dynamical lifetime of a few Myrs and much shorter than the age of the solar system. In light of this time constraint, we evaluate previously suggested ring formation pathways including collisional ejection and satellite disruption. We also investigate in detail a contrasting formation mechanism, the lofting of dust particles off Chariklo's surface into orbit via outflows of sublimating CO and/or N_2 triggered after Chariklo was scattered inward by giant planets. This latter scenario predicts that rings should be common among 100-km class Centaurs but rare among Kuiper belt objects and smaller Centaurs. It also predicts that Centaurs should show seasonal variations in cometary activity with activity maxima occurring shortly after equinox.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SMACK: A New Algorithm for Modeling Collisions and Dynamics of Planetesimals in Debris Disks

We present the Superparticle Model/Algorithm for Collisions in Kuiper belts and debris disks (SMACK), a new method for simultaneously modeling, in 3-D, the collisional and dynamical evolution of planetesimals in a debris disk with planets. SMACK can simulate azimuthal asymmetries and how these asymmetries evolve over time. We show that SMACK is stable to numerical viscosity and numerical heating over 10^7 yr, and that it can reproduce analytic models of disk evolution. We use SMACK to model the evolution of a debris ring containing a planet on an eccentric orbit. Differential precession creates a spiral structure as the ring evolves, but collisions subsequently break up the spiral, leaving a narrower eccentric ring.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Care and feeding of frogs

"Propellers" are features in Saturn's A ring associated with moonlets that open partial gaps. They exhibit non-Keplerian motion (Tiscareno 2010); the longitude residuals of the best-observed propeller, "Blériot," appear consistent with a sinusoid of period ~4 years. Pan and Chiang (2010) proposed that propeller moonlets librate in "frog resonances" with co-orbiting ring material. By analogy with the restricted three-body problem, they treated the co-orbital material as stationary in the rotating frame and neglected non-co-orbital material. Here we use simple numerical experiments to extend the frog model, including feedback due to the gap's motion, and drag associated with the Lindblad disk torques that cause Type I migration. Because the moonlet creates the gap, we expect the gap centroid to track the moonlet, but only after a time delay t_diff, the time for a ring particle to travel from conjunction with the moonlet to the end of the gap. We find that frog librations can persist only if t_diff exceeds the frog libration period P_lib, and if damping from Lindblad torques balances driving from co-orbital torques. If t_diff << P_lib, then the libration amplitude damps to zero. In the case of Blériot, the frog resonance model can reproduce the observed libration period P_lib ~ 4 yr. However, our simple feedback prescription suggests that Blériot's t_diff ~ 0.01P_lib, which is inconsistent with the observed libration amplitude of 260 km. We urge more accurate treatments of feedback to test the assumptions of our toy models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic flights of propellers

Kilometer-sized moonlets in Saturn's A ring create S-shaped wakes called "propellers" in surrounding material. The Cassini spacecraft has tracked the motions of propellers for several years and finds that they deviate from Keplerian orbits having constant semimajor axes. The inferred orbital migration is known to switch sign. We show using a statistical test that the time series of orbital longitudes of the propeller Blériot is consistent with that of a time-integrated Gaussian random walk. That is, Blériot's observed migration pattern is consistent with being stochastic. We further show, using a combination of analytic estimates and collisional N-body simulations, that stochastic migration of the right magnitude to explain the Cassini observations can be driven by encounters with ring particles 10-20 m in radius. That the local ring mass is concentrated in decameter-sized particles is supported on independent grounds by occultation analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-consistent size and velocity distributions of collisional cascades

The standard theoretical treatment of collisional cascades derives a steady-state size distribution assuming a single constant velocity dispersion for all bodies regardless of size. Here we relax this assumption and solve self-consistently for the bodies' steady-state size and size-dependent velocity distributions. Specifically, we account for viscous stirring, dynamical friction, and collisional damping of the bodies' random velocities in addition to the mass conservation requirement typically applied to find the size distribution in a steady-state cascade. The resulting size distributions are significantly steeper than those derived without velocity evolution. For example, accounting self-consistently for the velocities can change the standard q=3.5 power-law index of the Dohnanyi (1969) differential size spectrum to an index as large as q=4. Similarly, for bodies held together by their own gravity, the corresponding power-law index range 2.88<q<3.14 of Pan & Sari (2005) can steepen to values as large as q=3.26. Our velocity results allow quantitative predictions of the bodies' scale heights as a function of size. Together with our predictions, observations of the scale heights for different sized bodies for the Kuiper belt, the asteroid belt, and extrasolar debris disks may constrain the mass and number of large bodies stirring the cascade as well as the colliding bodies' internal strengths.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Observational biases in determining extrasolar planet eccentricities in single-planet systems

We investigate potential biases in the measurements of exoplanet orbital parameters obtained from radial velocity observations for single-planet systems. We create a mock catalog of radial velocity data, choosing input planet masses, periods, and observing patterns from actual radial velocity surveys and varying input eccentricities. We apply Markov Chain Monte Carlo (MCMC) simulations and compare the resulting orbital parameters to the input values. We find that a combination of the effective signal-to-noise ratio of the data, the maximal gap in phase coverage, and the total number of periods covered by observations is a good predictor of the quality of derived orbit parameters. As eccentricity is positive definite, we find that eccentricities of planets on nearly circular orbits are preferentially overestimated, with typical bias of 1-2 times the median eccentricity uncertainty in a survey (e.g., 0.04 in the Butler et al. 2006 catalog). When performing population analysis, we recommend using the mode of the marginalized posterior eccentricity distribution to minimize potential biases. While the Butler et al. (2006) catalog reports eccentricities below 0.05 for just 17% of single-planet systems, we estimate that the true fraction of e < 0.05 orbits is about f(0.05)=38\pm 9%. For planets with P > 10 days, we find f(0.05)=28\pm 8% versus 10% from Butler et al. (2006). These planets either never acquired a large eccentricity or were circularized following any significant eccentricity excitation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The propeller and the frog

"Propellers" in planetary rings are disturbances in ring material excited by moonlets that open only partial gaps. We describe a new type of co-orbital resonance that can explain the observed non-Keplerian motions of propellers. The resonance is between the moonlet underlying the propeller, and co-orbiting ring particles downstream of the moonlet where the gap closes. The moonlet librates within the gap about an equilibrium point established by co-orbiting material and stabilized by the Coriolis force. In the limit of small libration amplitude, the libration period scales linearly with the gap azimuthal width and inversely as the square root of the co-orbital mass. The new resonance recalls but is distinct from conventional horseshoe and tadpole orbits; we call it the "frog" resonance, after the relevant term in equine hoof anatomy. For a ring surface density and gap geometry appropriate for the propeller Blériot in Saturn's A ring, our theory predicts a libration period of ~4 years, similar to the ~3.7 year period over which Blériot's orbital longitude is observed to vary. These librations should be subtracted from the longitude data before any inferences about moonlet migration are made.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Composite self-similar solutions for relativistic shocks: the transition to cold fluid temperatures

The flow resulting from a strong ultrarelativistic shock moving through a stellar envelope with a polytrope-like density profile has been studied analytically and numerically at early times while the fluid temperature is relativistic--that is, just before and just after the shock breaks out of the star. Such a flow should expand and accelerate as its internal energy is converted to bulk kinetic energy; at late enough times, the assumption of relativistic temperatures becomes invalid. Here we present a new self-similar solution for the post-breakout flow when the accelerating fluid has bulk kinetic Lorentz factors much larger than unity but is cooling through $p/n$ of order unity to subrelativistic temperatures. This solution gives a relation between a fluid element's terminal Lorentz factor and that element's Lorentz factor just after it is shocked. Our numerical integrations agree well with the solution. While our solution assumes a planar flow, we show that corrections due to spherical geometry are important only for extremely fast ejecta originating in a region very close to the stellar surface. This region grows if the shock becomes relativistic deeper in the star.

0 citations0 reviewsNo code linkedOpen access