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

Alice C. Quillen

Alice C. Quillen 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-ph.GAastro-phastro-ph.COastro-ph.SRBiological Physicscond-mat.softphysics.flu-dynphysics.geo-ph

Catalog footprint

What is connected

40works
9topics
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

40 published item(s)

preprint2022arXiv

Ricochets on Asteroids II: Sensitivity of laboratory experiments of low velocity grazing impacts on substrate grain size

We compare low velocity impacts that ricochet with the same impact velocity and impact angle into granular media with similar bulk density, porosity and friction coefficient but different mean grain size. The ratio of projectile diameter to mean grain length ranges from 4 in our coarsest medium to 50 in our finest sand. Using high speed video and fluorescent markers, we measure the ratio of pre- to post-impact horizontal and vertical velocity components, which we refer to as coefficients of restitution, and the angle of deflection caused by the impact in the horizontal plane. Coefficients of restitution are sensitive to mean grain size with the ratio associated with the horizontal velocity component about twice as large for our coarsest gravel as that for our finest sand. This implies that coefficients for hydro-static-like, drag-like and lift-like forces, used in empirical force laws, are sensitive to mean grain size. The coefficient that is most strongly sensitive to grain size is the lift coefficient which decreases by a factor of 3 between our coarsest and finest media. The deflection angles are largest in the coarser media and their size approximately depends on grain size to the 3/2 power. This scaling is matched with a model where momentum transfer takes place via collisions with individual grains. The dependence of impact mechanics on substrate size distribution should be considered in future models for populations of objects that impact granular asteroid surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Sub-surface granular dynamics in the context of oblique, low-velocity impacts into angular granular media

Oblique, low-velocity impacts onto extraterrestrial terrain are an inevitable occurrence during space exploration. We conduct two-dimensional discrete simulations to model such impacts into a bed of triangular grains. Finite element method provides the basis for simulation, enabling the angular grain geometry. Our findings re-create the three classes of impact behavior previously noted from experiments: full-stop, rollout, and ricochet \citep*{Wright2020}. An application of Set Voronoi tessellation assesses packing fraction at a high resolution, revealing how grains shift relative to each other during an impact event. Calculation of Von Mises strain distributions then reveal how grains shift relative to the overall system, leading to the notion of the 'skin zone'. Intuition would suggest that the region of perturbed grains would grow deeper with higher velocity impacts, results instead show that increasing velocity may actually evoke a change in the grains' dissipative response that boosts lateral perturbation. Finally, we consider as a whole how sub-surface response could link with impactor dynamics to deepen our understanding of oblique, low-velocity impact events and help to improve mission outcomes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Accretion of Ornamental Equatorial Ridges on Pan, Atlas and Daphnis

We explore scenarios for the accretion of ornamental ridges on Saturn's moons Pan, Atlas, and Daphnis from material in Saturn's rings. Accretion of complex shaped ridges from ring material should be possible when the torque from accreted material does not exceed the tidal torque from Saturn that ordinarily maintains tidal lock. This gives a limit on the maximum accretion rate and the minimum duration for equatorial ridge growth. We explore the longitude distribution of ridges accreted from ring material, initially in circular orbits, onto a moon that is on a circular, inclined or eccentric orbit. Sloped and lobed ridges can be accreted if the moon tidally realigns during accretion due to its change in shape or because the disk edge surface density profile allows ring material originating at different initial semi-major axes to impact the moon at different locations on its equatorial ridge. We find that accretion from an asymmetric gap might account for a depression on Atlas's equatorial ridge. Accretion from an asymmetric gap at orbital eccentricity similar to the Hill eccentricity, might allow accretion of multiple lobes, as seen on Pan. Two possibly connected scenarios are promising for growth of ornamental equatorial ridges. The moon migrates through the ring, narrowing its gap and facilitating accretion. The moon's orbital eccentricity could increase due to orbital resonance with another moon, pushing it into its gap edges and facilitating accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Synchronized oscillations in swarms of nematode Turbatrix aceti

There is a recent surge of interest in the behavior of active particles that can at the same time align their direction of movement and synchronize their oscillations, known as swarmalators. While theoretical and numerical models of such systems are now abundant, no real-life examples have been shown to date. We present an experimental investigation of the collective motion of the nematode Turbatrix aceti that self-propel by body undulation. We discover that these nematodes can synchronize their body oscillations, forming striking traveling metachronal waves, which produces strong fluid flows. We uncover that the location and strength of this collective state can be controlled through the shape of the confining structure; in our case the contact angle of a droplet. This opens a way for producing controlled work such as on-demand flows or displacement of objects. We illustrate this by showing that the force generated by this state is sufficient to change the physics of evaporation of fluid droplets, by counteracting the surface-tension force, which allow us to estimate its strength. The relatively large size and ease of culture make Turbatrix aceti a promising model organism for experimental investigation of swarming and oscillating active matter capable of producing controllable work.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Excitation of Tumbling in Phobos and Deimos

Mass-spring model simulations are used to investigate past spin states of a viscoelastic Phobos and Deimos. From an initially tidally locked state, we find crossing of a spin-orbit resonance with Mars or a mean motion resonance with each other does not excite tumbling in Phobos or Deimos. However, once tumbling our simulations show that these moons can remain so for an extended period and during this time their orbital eccentricity can be substantially reduced. We attribute the tendency for simulations of an initially tumbling viscoelastic body to drop into spin-synchronous state at very low eccentricity to the insensitivity of the tumbling chaotic zone volume to eccentricity. After a tumbling body enters the spin synchronous resonance, it can exhibit long lived non-principal axis rotation and this too can prolong the period of time with enhanced tidally generated energy dissipation. The low orbital eccentricities of Phobos and Deimos could in part be due to spin excitation by nearly catastrophic impacts rather than tidal evolution following orbital resonance excitation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crustal Failure on Icy Moons from a Strong Tidal Encounter

Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress on the surface to cause brittle failure of an icy crust and simulated fractures can extend a large fraction of the radius of body. Strong tidal encounters may be responsible for the formation of long graben complexes and chasmata in ancient terrain of icy moons such as Dione, Tethys, Ariel and Charon.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Simulation of Tidal Evolution of a Viscoelastic Body Modelled with a Mass-Spring Network

We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Powerful Activity in the Bright Ages. I. A Visible/IR Survey of High Redshift 3C Radio Galaxies and Quasars

We present new rest frame UV and visible observations of 22 high-redshift (1 < z < 2.5) 3C radio galaxies and QSOs obtained with the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) instrument. Using a custom data reduction strategy in order to assure the removal of cosmic rays, persistence signal, and other data artifacts, we have produced high-quality science-ready images of the targets and their local environments. We observe targets with regions of UV emission suggestive of active star formation. In addition, several targets exhibit highly distorted host galaxy morphologies in the rest frame visible images. Photometric analyses reveals that brighter QSOs tend to be generally redder than their dimmer counterparts. Using emission line fluxes from the literature, we estimate that emission line contamination is relatively small in the rest frame UV images for the QSOs. Using archival VLA data, we have also created radio map overlays for each of our targets, allowing for analysis of the optical and radio axes alignment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tidal spin down rates of homogeneous triaxial viscoelastic bodies

We use numerical simulations to measure the sensitivity of the tidal spin down rate of a homogeneous triaxial ellipsoid to its axis ratios by comparing the drift rate in orbital semi-major axis to that of a spherical body with the same mass, volume and simulated rheology. We use a mass-spring model approximating a viscoelastic body spinning around its shortest body axis, with spin aligned with orbital spin axis, and in circular orbit about a point mass. The torque or drift rate can be estimated from that predicted for a sphere with equivalent volume if multiplied by $0.5 (1 + b^4/a^4)(b/a)^{-4/3} (c/a)^{-α_c}$ where $b/a$ and $c/a$ are the body axis ratios and index $α_c \approx 1.05$ is consistent with the random lattice mass spring model simulations but $α_c = 4/3$ suggested by scaling estimates. A homogeneous body with axis ratios 0.5 and and 0.8, like Haumea, has orbital semi-major axis drift rate about twice as fast as a spherical body with the same mass, volume and material properties. A simulation approximating a mostly rocky body but with 20\% of its mass as ice concentrated at its ends has a drift rate 10 times faster than the equivalent homogeneous rocky sphere. However, this increase in drift rate is not enough to allow Haumea's satellite, Hi'iaka, to have tidally drifted away from Haumea to its current orbital semi-major axis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Far Ultraviolet Morphology of Star Forming Filaments in Cool Core Brightest Cluster Galaxies

We present a multiwavelength morphological analysis of star forming clouds and filaments in the central ($< 50$ kpc) regions of 16 low redshift ($z<0.3$) cool core brightest cluster galaxies (BCGs). New Hubble Space Telescope (HST) imaging of far ultraviolet continuum emission from young ($\sim 10$ Myr), massive ($> 5$ \Msol) stars reveals filamentary and clumpy morphologies, which we quantify by means of structural indices. The FUV data are compared with X-ray, Ly$α$, narrowband H$α$, broadband optical/IR, and radio maps, providing a high spatial resolution atlas of star formation locales relative to the ambient hot ($\sim10^{7-8}$ K) and warm ionised ($\sim 10^4$ K) gas phases, as well as the old stellar population and radio-bright AGN outflows. Nearly half of the sample possesses kpc-scale filaments that, in projection, extend toward and around radio lobes and/or X-ray cavities. These filaments may have been uplifted by the propagating jet or buoyant X-ray bubble, or may have formed {\it in situ} by cloud collapse at the interface of a radio lobe or rapid cooling in a cavity's compressed shell. The morphological diversity of nearly the entire FUV sample is reproduced by recent hydrodynamical simulations in which the AGN powers a self-regulating rain of thermally unstable star forming clouds that precipitate from the hot atmosphere. In this model, precipitation triggers where the cooling-to- freefall time ratio is $t_{\mathrm{cool}}/t_{\mathrm{ff}}\sim 10$. This condition is roughly met at the maxmial projected FUV radius for more than half of our sample, and clustering about this ratio is stronger for sources with higher star formation rates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase Wrapping of Epicyclic Perturbations in the Wobbly Galaxy

We use test-particle integrations to show that epicyclic motions excited by a pericentre passage of a dwarf galaxy could account for bulk vertical velocity streaming motions recently observed in the Galactic stellar disc near the Sun. We use fixed potential test-particle integrations to isolate the role of phase wrapping of epicyclic perturbations from bending and breathing waves or modes, which require self-gravity to oscillate. Perturbations from a fairly massive Sagittarius dwarf galaxy, $M_d \sim 2.5 \times 10^{10} M_\odot$, are required to account for the size of the observed streaming motions from its orbital pericentre approximately a Gyr ago. A previous passage of the dwarf through the Galactic disc approximately 2.2 Gyr ago (with a then more massive dwarf galaxy) is less effective. If phase wrapping of epicyclic perturbations is responsible for stellar streaming motions in the Galactic disc, then there should be variations in velocity gradients on scales of a few kpc in the vicinity of the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Parent Populations of 6 groups identified from Chemical Tagging in the Solar neighborhood

We estimate the size and distribution of the parent populations for the 6 largest (at least 20 stars in the Solar neighborhood) chemical groups identified in the Chemical Tagging experiment by Mitschang et al.~2014. Stars in the abundance groups tend to lie near a boundary in angular momentum versus eccentricity space where the probability is highest for a star to be found in the Solar neighborhood and where orbits have apocenter approximately equal to the Sun's galactocentric radius. Assuming that the parent populations are uniformly distributed at all azimuthal angles in the Galaxy, we estimate that the parent populations of these abundance groups contain at least 200,000 members. The spread in angular momentum of the groups implies that the assumption of a uniform azimuthal distribution only fails for the two youngest groups and only for the highest angular momentum stars in them. The parent populations of three thin disk groups have narrow angular momentum distributions, but tails in the eccentricity and angular momentum distributions suggest that only a small fraction of stars have migrated and increased in eccentricity. In contrast, the parent populations of the thick disk groups exhibit both wide angular momentum and eccentricity distributions implying that both heating and radial migration has taken place.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Torque on an exoplanet from an anisotropic evaporative wind

Winds from short-period Earth and Neptune mass exoplanets, driven by high energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super Earth mass planet that loses a large fraction of its mass in a wind could drift a few percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A search for eclipsing binaries that host discs

We search for systems hosting eclipsing discs using a complete sample of eclipsing binaries (EBs); those previously identified in the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). Within a subsample of 2,823 high-cadence, high-photometric precision and large eclipsing depth detached EBs previously identified in the Large Magellanic Cloud (LMC), we find that the skewness and kurtosis of the light curves magnitude distribution within the primary eclipse can distinguish EBs hosting a disc from those without. Two systems with previously identified eclipsing discs (OGLE-LMC-ECL-11893 and OGLE-LMC-ECL-17782) are identified with near zero skewness ($|S|<0.5$) and positive kurtosis. No additional eclipsing disc systems were found in the OGLE-III LMC, Small Magellanic Cloud (SMC) or Galactic Disc (GD) EB light curves. We estimate that the fraction of detached near main-sequence LMC EBs (which have a primary with an $I$-band magnitude brighter than $\simeq 19\,\rm{mag}$) that host a disc is approximately 1/1000. As circumstellar disc lifetimes are short, we expected to primarily find eclipsing discs around young stars. In addition, as there is more room for a disc in a widely separated binary and because a disk close to a luminous star would be above the dust sublimation temperature, we expected to primarily find eclipsing discs in long period binaries. However, OGLE-LMC-ECL-17782 is a 13.3 day period B star system with a transient and hot ($\sim 6000\,\rm{K}$, $\sim 0.1\,\rm{AU}$ radius) disc and Scott et al. (in prep.) estimate an age of $150\,\rm{Myr}$ for OGLE-LMC-ECL-11893. Both discs are unexpected in the EB sample and impel explanation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling Transiting Circumstellar Disks: Characterizing the Newly Discovered Eclipsing Disk System OGLE LMC-ECL-11893

We investigate the nature of the unusual eclipsing star OGLE LMC-ECL-11893 (OGLE J05172127-6900558) in the Large Magellanic Cloud recently reported by Dong et al. 2014. The eclipse period for this star is 468 days, and the eclipses exhibit a minimum of ~1.4 mag, preceded by a plateau of ~0.8 mag. Spectra and optical/IR photometry are consistent with the eclipsed star being a lightly reddened B9III star of inferred age ~150 Myr and mass of ~4 solar masses. The disk appears to have an outer radius of ~0.2 AU with predicted temperatures of ~1100-1400 K. We model the eclipses as being due to either a transiting geometrically thin dust disk or gaseous accretion disk around a secondary object; the debris disk produces a better fit. We speculate on the origin of such a dense circumstellar dust disk structure orbiting a relatively old low-mass companion, and on the similarities of this system to the previously discovered EE Cep.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Resonant Chains and Three-body Resonances in the Closely-Packed Inner Uranian Satellite System

Numerical integrations of the closely-packed inner Uranian satellite system show that variations in semi-major axes can take place simultaneously between three or four consecutive satellites. We find that the three-body Laplace angle values are distributed unevenly and have histograms showing structure, if the angle is associated with a resonant chain, with both pairs of bodies near first-order two-body resonances. Estimated three-body resonance libration frequencies can be only an order of magnitude lower than those of first-order resonances. Their strength arises from a small divisor from the distance to the first-order resonances and insensitivity to eccentricity, which make up for their dependence on moon mass. Three-body resonances associated with low-integer Laplace angles can also be comparatively strong due to the many multiples of the angle contributed from Fourier components of the interaction terms. We attribute small coupled variations in semi-major axis, seen throughout the simulation, to ubiquitous and weak three-body resonant couplings. We show that a system with two pairs of bodies in first-order mean-motion resonance can be transformed to resemble the well-studied periodically-forced pendulum with the frequency of a Laplace angle serving as a perturbation frequency. We identify trios of bodies and overlapping pairs of two-body resonances in each trio that have particularly short estimated Lyapunov timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variability in the 2MASS Calibration Fields: A Search for Transient Obscuration Events

We searched the light curves of over 40000 stars in the 2MASS calibration database, spanning approximately 4 years, for objects that have significant day long dimming events. We also searched the multi-color light curves for red dimming events that could be due to transient extinction. In the color independent sigma-limited search, we found 46 previously unknown eclipsing binaries, 6 periodic variable stars likely to be intrinsic pulsators and 21 young stellar objects in the $ρ$ Ophiucus star formation region previously studied by Parks et al. 2013. An additional 11 objects exhibited dimming events and most of these are unclassified. The search for red dimming events primarily reveals a population of low luminosity active galaxies that become bluer when they are brighter, and variable young stellar objects exhibiting high cross-correlation coefficients between color and brightness. The young stellar objects exhibit brightness and color variations in the direction of interstellar extinction whereas the active galaxies can have a bowed distribution in color and magnitude with reduced variation in color when the object is brightest. Among the objects that are usually quiescent (not strongly variable), we failed to find any dimming events deeper than 0.2 magnitude and lasting longer than a day. Two of the least embedded young stellar objects, however, dimmed by 0.2 mag for longer than a day without strong color variation. Having failed to find new exotic objects, we explored ways to eliminate commonly found objects so that a larger number of objects may be searched. We find that all but one eclipsing binary are excluded by requiring moderate color variation during a dimming event and most of the active galaxies and young stellar objects are excluded by placing a limit on the standard deviation of the magnitude distribution in the light curve.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Vertical Resonance Heating Model for X- or Peanut-Shaped Galactic Bulges

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of a planetesimal debris disk on stability scenarios for the extrasolar planetary system HR 8799

HR 8799 is a four planet system that also hosts a debris disk. By numerically integrating both planets and a planetesimal disk, we find interactions between an exterior planetesimal disk and the planets can influence the lifetime of the system. We first consider resonant planetary configurations that remained stable for at least 7 Myrs sans debris disk. An exterior debris disk with only 1 per cent the mass of the outermost planet (approximately a Neptune mass) was sufficiently large enough to pull the system out of resonance after 2 to 6 Myrs. Secondly, we consider configurations which are unstable in less than a few hundred thousand years. We find that these can be stabilized by a debris disk with a mass of more than 10 per cent that of the outermost planet. Our two sets of simulations suggest that estimates of the long term stability of a planetary system should take into account the role of the debris disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Origin Scenarios for the Kepler 36 Planetary System

We explore scenarios for the origin of two different density planets in the Kepler 36 system in adjacent orbits near the 7:6 mean motion resonance. We find that fine tuning is required in the stochastic forcing amplitude, the migration rate and planet eccentricities to allow two convergently migrating planets to bypass mean motion resonances such as the 4:3, 5:4 and 6:5, and yet allow capture into the 7:6 resonance. Stochastic forcing can eject the system from resonance causing a collision between the planets, unless the disk inducing migration and stochastic forcing is depleted soon after resonance capture. We explore a scenario with approximately Mars mass embryos originating exterior to the two planets and migrating inwards toward two planets. We find that gravitational interactions with embryos can nudge the system out of resonances. Numerical integrations with about a half dozen embryos can leave the two planets in the 7:6 resonance. Collisions between planets and embryos have a wide distribution of impact angles and velocities ranging from accretionary to disruptive. We find that impacts can occur at sufficiently high impact angle and velocity that the envelope of a planet could have been stripped, leaving behind a dense core. Some of our integrations show the two planets exchanging locations, allowing the outer planet that had experienced multiple collisions with embryos to become the innermost planet. A scenario involving gravitational interactions and collisions with embryos may account for both the proximity of the Kepler 36 planets and their large density contrast.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stability Boundaries for Resonant Migrating Planet Pairs

Convergent migration allows pairs of planet to become trapped into mean motion resonances. Once in resonance, the planets' eccentricities grow to an equilibrium value that depends on the ratio of migration time scale to the eccentricity damping timescale, $K=τ_a/τ_e$, with higher values of equilibrium eccentricity for lower values of $K$. For low equilibrium eccentricities, $e_{eq}\propto K^{-1/2}$. The stability of a planet pair depends on eccentricity so the system can become unstable before it reaches its equilibrium eccentricity. Using a resonant overlap criterion that takes into account the role of first and second order resonances and depends on eccentricity, we find a function $K_{min}(μ_p, j)$ that defines the lowest value for $K$, as a function of the ratio of total planet mass to stellar mass ($μ_p$) and the period ratio of the resonance defined as $P_1/P_2=j/(j+k)$, that allows two convergently migrating planets to remain stable in resonance at their equilibrium eccentricities. We scaled the functions $K_{min}$ for each resonance of the same order into a single function $K_c$. The function $K_{c}$ for planet pairs in first order resonances is linear with increasing planet mass and quadratic for pairs in second order resonances with a coefficient depending on the relative migration rate and strongly on the planet to planet mass ratio. The linear relation continues until the mass approaches a critical mass defined by the 2/7 resonance overlap instability law and $K_c \to \infty$. We compared our analytic boundary with an observed sample of resonant two planet systems. All but one of the first order resonant planet pair systems found by radial velocity measurements are well inside the stability region estimated by this model. We calculated $K_c$ for Kepler systems without well-constrained eccentricities and found only weak constraints on $K$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetary Construction Zones in Occultation: Discovery of an Extrasolar Ring System Transiting a Young Sun-like Star and Future Prospects for Detecting Eclipses by Circumsecondary and Circumplanetary Disks

<abridged> The large relative sizes of circumstellar and circumplanetary disks imply that they might be seen in eclipse in stellar light curves. We estimate that a survey of ~10^4 young (~10 Myr old) post-accretion pre-MS stars monitored for ~10 years should yield at least a few deep eclipses from circumplanetary disks and disks surrounding low mass companion stars. We present photometric and spectroscopic data for a pre-MS K5 star (1SWASP J140747.93-394542.6), a newly discovered ~0.9 Msun member of the ~16 Myr-old Upper Cen-Lup subgroup of Sco-Cen at a kinematic distance of 128 pc. SuperWASP and ASAS light curves for this star show a remarkably long, deep, and complex eclipse event centered on 29 April 2007. At least 5 multi-day dimming events of >0.5 mag are identified, with a >3.3 mag deep eclipse bracketed by two pairs of ~1 mag eclipses symmetrically occurring +-12 days and +-26 days before and after. Hence, significant dimming of the star was taking place on and off over at least a ~54 day period in 2007, and a strong >1 mag dimming event occurred over a ~12 day span. We place a firm lower limit on the period of 850 days (i.e. the orbital radius of the eclipser must be >1.7 AU and orbital velocity must be <22 km/s). The shape of the light curve is similar to the lop-sided eclipses of the Be star EE Cep. We suspect that this new star is being eclipsed by a low-mass object orbited by a dense inner disk, girded by at least 3 dusty rings of lower optical depth. Between these rings are at least two annuli of near-zero optical depth (i.e. gaps), possibly cleared out by planets or moons, depending on the nature of the secondary. For possible periods in the range 2.33-200 yr, the estimated total ring mass is ~8-0.4 Mmoon (if the rings have optical opacity similar to Saturn's rings), and the edge of the outermost detected ring has orbital radius ~0.4-0.09 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stellar Migration by Short Lived Density Peaks Arising from Interference of Spiral Density Waves in an N-body Simulation

We identify migrating stars in an N-body simulation of a Milky-Way-like disk. Outward migration can occur when a star in a low eccentricity orbit lags a short-lived local spiral arm density peak. We interpret short lived local density peaks, that appear and fade on approximately an orbital period, as arising from positive interference between spiral density wave patterns. Stars near such a peak can migrate over a significant distance in galactocentric radius during the peak lifetime, providing that the peak is sufficiently dense. Using a Gaussian bar model for the potential perturbation associated with a narrow transient spiral feature, estimates of the migration rate, angular offset between particle and spiral feature, and maximum eccentricity for migrators roughly agrees with the values measured in our simulation. When multiple spiral density waves are present, local density peaks can appear and disappear on timescales faster than the timescale estimated for growth and decay of individual waves and the peak surface density can be larger than for any individual wave. Consequently, migration induced by transient density peaks may be more pervasive than that mediated by the growth and decay of individual patterns and occurring at their corotation resonance. We discuss interpretation of transient-like behavior in terms of interfering patterns, including estimating a coherence time for features that appear due to constructive interference, their effective angular rotation rates and the speed and direction that a density maximum would move across a galaxy inducing a localized and traveling burst of star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Capture of Irregular Satellites via Binary Planetesimal Exchange Reactions in Migrating Planetary Systems

By logging encounters between planetesimals and planets we compute the distribution of encounters in a numerically integrated two planet system that is migrating due to interactions with an exterior planetesimal belt. Capture of an irregular satellite in orbit about a planet through an exchange reaction with a binary planetesimal is only likely when the binary planetesimal undergoes a slow and close encounter with the planet. In our simulations we find that close and slow encounters between planetesimals and a planet primarily occur with the outermost and not innermost planet. Taking care to consider where a planet orbit crossing binary planetesimal would first be tidally disrupted, we estimate the probability of both tidal disruption and irregular satellite capture. We estimate that the probability that the secondary of a binary planetesimal is captured and becomes an irregular satellite about a Neptune mass outer planet is about 1/100 for binaries with masses and separations similar to transneptunian planetesimal binaries. If young exoplanetary debris disks host a binary planetesimal population then outwards migrating outer planets should host captured irregular satellite populations. We discuss interpretation of emission associated with the exoplanet Fomalhaut b in terms of collisional evolution of a captured irregular satellite population that is replenished due to planetary migration.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling disc non-axisymmetries: multiple patterns, radial migration, and thick discs

Disc non-axisymmetrc components, such as spirals and central bars, are nowadays known to play an important role in shaping galactic discs. Here we use Tree-SPH N-body simulations to examine the effect of these perturbers on two aspects: the occurrence of multiple patterns in discs and the effects of radial migration on disc thickening. We find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop in all models. Interaction among these asymmetric features results in a large scale stellar migration. However, we show that, despite the strong radial mixing, discs cannot be thickened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as stars migrate radially outwards. We also find that the bulge contribution to a thick-disc component for an Sa-type galaxy at ~2.5 disc scale-lengths is less than 1% and zero in the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick disc formation via stellar radial migration.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Structure in phase space associated with spiral and bar density waves in an N-body galactic disk

An N-body hybrid simulation, integrating both massive and tracer particles, of a Galactic disk is used to study the stellar phase space distribution or velocity distributions in different local neighborhoods. Pattern speeds identified in Fourier spectrograms suggest that two-armed and three-armed spiral density waves, a bar and a lopsided motion are coupled in this simulation, with resonances of one pattern lying near resonances of other patterns. We construct radial and tangential (uv) velocity distributions from particles in different local neighborhoods. More than one clump is common in these local velocity distributions regardless of the position in the disk. Features in the velocity distribution observed at one galactic radius are also seen in nearby neighborhoods (at larger and smaller radii) but with shifted mean v values. This is expected if the v velocity component of a clump sets the mean orbital galactic radius of its stars. We find that gaps in the velocity distribution are associated with the radii of kinks or discontinuities in the spiral arms. These gaps also seem to be associated with Lindblad resonances with spiral density waves and so denote boundaries between different dominant patterns in the disk. We discuss implications for interpretations of the Milky Way disk based on local velocity distributions. Velocity distributions created from regions just outside the bar's Outer Lindblad resonance and with the bar oriented at 45 degrees from the Sun-Galactic center line more closely resemble that seen in the solar neighborhood (triangular in shape at lower uv and with a Hercules like stream) when there is a strong nearby spiral arm, consistent with the observed Centaurus Arm tangent, just interior to the solar neighborhood.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Three Body Resonance Overlap in Closely Spaced Multiple Planet Systems

We compute the strengths of zero-th order (in eccentricity) three-body resonances for a co-planar and low eccentricity multiple planet system. In a numerical integration we illustrate that slowly moving Laplace angles are matched by variations in semi-major axes among three bodies with the outer two bodies moving in the same direction and the inner one moving in the opposite direction, as would be expected from the two quantities that are conserved in the three-body resonance. A resonance overlap criterion is derived for the closely and equally spaced, equal mass system with three-body resonances overlapping when interplanetary separation is less than an order unity factor times the planet mass to the one quarter power. We find that three-body resonances are sufficiently dense to account for wander in semi-major axis seen in numerical integrations of closely spaced systems and they are likely the cause of instability of these systems. For interplanetary separations outside the overlap region, stability timescales significantly increase. Crudely estimated diffusion coefficients in eccentricity and semi-major axis depend on a high power of planet mass and interplanetary spacing. An exponential dependence previously fit to stability or crossing timescales is likely due to the limited range of parameters and times possible in integration and the strong power law dependence of the diffusion rates on these quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

HST FUV Observations of Brightest Cluster Galaxies: The Role of Star Formation in Cooling Flows and BCG Evolution

Quillen et al. and O'Dea et al. carried out a Spitzer study of a sample of 62 brightest cluster galaxies (BCGs) from the ROSAT Brightest Cluster Sample chosen based on their elevated H-alpha flux. We present Hubble Space Telescope Advanced Camera for Surveys (ACS) far ultraviolet (FUV) images of the Ly-alpha and continuum emission of the luminous emission-line nebulae in 7 BCGs found to have an infrared excess. We confirm that the BCGs are actively forming stars suggesting that the IR excess seen in these BCGs is indeed associated with star formation. The FUV continuum emission extends over a region of ~7-28 kpc (largest linear size) and even larger in Ly-alpha. The young stellar population required by the FUV observations would produce a significant fraction of the ionizing photons required to power the emission line nebulae. Star formation rates estimated from the FUV continuum range from ~3 to ~14 times lower than those estimated from the IR, however both the Balmer decrement in the central few arcseconds and detection of CO in most of these galaxies imply that there are regions of high extinction that could have absorbed much of the FUV continuum. Analysis of archival VLA observations reveals compact radio sources in all 7 BCGs and kpc scale jets in Abell 1835 and RXJ 2129+00. The four galaxies with archival deep Chandra observations exhibit asymmetric X-ray emission, the peaks of which are offset from the center of the BCGs by ~10 kpc on average. A low feedback state for the AGN could allow increased condensation of the hot gas into the center of the galaxy and the feeding of star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jeans Instability in a Tidally Disrupted Halo Satellite Galaxy

We use a hybrid test particle/N-body simulation to integrate 4 million massless test particle trajectories within a fully self-consistent 10^5 particle N-body simulation. The number of massless particles allows us to resolve fine structure in the spatial distribution and phase space of a dwarf galaxy as it is disrupted in the tidal field of a Milky Way type galaxy. The tidal tails exhibit nearly periodic clumping or a smoke-like appearance. By running simulations with different satellite particle mass, halo particle mass, number of massive and massless particles and with and without a galaxy disk, we have determined that the instabilities are not due to numerical noise, amplification of structure in the halo, or shocking as the satellite passes through the disk of the Galaxy. We measure Jeans wavelengths and growth timescales in the tidal tail and show that the Jeans instability is a viable explanation for the clumps. We find that the instability causes velocity perturbations of order 10 km/s. Clumps in tidal tails present in the Milky Way could be seen in stellar radial velocity surveys as well as number counts. We find that the unstable wavelength growth is sensitive to the simulated mass of dark matter halo particles. A simulation with a smoother halo exhibits colder and thinner tidal tails with more closely spaced clumps than a simulation with more massive dark matter halo particles. Heating by the halo particles increases the Jeans wavelength in the tidal tail affecting substructure development, suggesting an intricate connection between tidal tails and dark matter halo substructure.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jeans Instability of Palomar 5's Tidal Tail

Tidal tails composed of stars should be unstable to the Jeans instability and this can cause them to look like beads on a string. The Jeans wavelength and tail diameter determine the wavelength and growth rate of the fastest growing unstable mode. Consequently the distance along the tail to the first clump and spacing between clumps can be used to estimate the mass density in the tail and its longitudinal velocity dispersion. Clumps in the tidal tails of the globular cluster Palomar 5 could be due to Jeans instability. We find that their spacing is consistent with the fastest growing mode if the velocity dispersion in the tail is similar to that in the cluster itself. While all tidal tails should exhibit gravitational instability, we find that clusters or galaxies with low concentration parameters are most likely to exhibit short wavelength rapidly growing Jeans modes in their tidal tails.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radial mixing in the outer Milky Way disk caused by an orbiting satellite

Using test particle simulations we examine the structure of the outer Galactic disk as it is perturbed by a satellite in a tight eccentric orbit about the Galaxy. A satellite of mass a few times 10^9 Msol can heat the outer Galactic disk, excite spiral structure and a warp and induce streams in the velocity distribution. We examine particle eccentricity versus the change in mean radius between initial and current orbits. Correlations between these quantities are reduced after a few satellite pericenter passages. Stars born in the outer galaxy can be moved in radius from their birth positions and be placed in low eccentricity orbits inside their birth radii. We propose that mergers and perturbations from satellite galaxies and subhalos can induce radial mixing in the stellar metallicity distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Warped Disk of Centaurus A from a radius of 2 to 6500 pc

We compile position and inclination angles for tilted ring fits to the warped dusty and gaseous disk spanning radius 1.8 to 6500 pc from recent observations. For radii exterior to a kpc, tilted ring orientations lie on an arc on a polar inclination versus position angle plot, suggesting that precession following a merger can account for the ring morphology. Three kinks in the ring orientations are seen on the polar plot, the one at radius of about 1.3 kpc we suspect corresponds to the location where self-gravity in the disk affects the ring precession rate. Another at a radius of about 600 pc may be associated with a gap in the gas distribution. A third kink is seen at a radius of 100 pc. A constant inclination tilted disk precessing about the jet axis may describe the disk between 100 and 20 pc but not interior to this. A model with disk orientation matching the molecular circumnuclear disk at 100 pc that decays at smaller radii to an inner flat disk perpendicular to the jet may account for disk orientations within 100 pc. Neither model would account for the cusps or changes in disk orientation at 100 or 600 pc.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Outflow Driven Turbulence in Molecular Clouds

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

When is star formation episodic? A delay differential equation negative feedback model

We introduce a differential equation for star formation in galaxies that incorporates negative feedback with a delay. When the feedback is instantaneous, solutions approach a self-limiting equilibrium state. When there is a delay, even though the feedback is negative, the solutions can exhibit cyclic and episodic solutions. We find that periodic or episodic star formation only occurs when two conditions are satisfied. Firstly the delay timescale must exceed a cloud consumption timescale. Secondly the feedback must be strong. This statement is quantitatively equivalent to requiring that the timescale to approach equilibrium be greater than approximately twice the cloud consumption timescale. The period of oscillations predicted is approximately 4 times the delay timescale. The amplitude of the oscillations increases with both feedback strength and delay time. We discuss applications of the delay differential equation (DDE) model to star formation in galaxies using the cloud density as a variable. The DDE model is most applicable to systems that recycle gas and only slowly remove gas from the system. We propose likely delay mechanisms based on the requirement that the delay time is related to the observationally estimated time between episodic events. The proposed delay timescale accounting for episodic star formation in galaxy centers on periods similar to P 10 Myrs, irregular galaxies with P 100 Myrs, and the Milky Way disk with P~ 2Gyr, could be that for exciting turbulence following creation of massive stars, that for gas pushed into the halo to return and interact with the disk and that for spiral density wave evolution, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Planetary embryos and planetesimals residing in thin debris disks

We consider constraints on the planetesimal population residing in the disks of AU Microscopii, Beta Pictoris and Fomalhaut taking into account their observed thicknesses and normal disk opacities. We estimate that bodies of radius 5, 180 and 70 km are responsible for initiating the collisional cascade accounting for the dust production for AU-Mic, Beta-Pic and Fomalhaut's disks, respectively, at break radii from the star where their surface brightness profiles change slope. Larger bodies, of radius 1000km and with surface density of order 0.01 g/cm^2, are required to explain the thickness of these disks assuming that they are heated by gravitational stirring. A comparison between the densities of the two sizes suggests the size distribution in the largest bodies is flatter than that observed in the Kuiper belt. AU Mic's disk requires the shallowest size distribution for bodies with radius greater than 10km suggesting that the disk contains planetary embryos experiencing a stage of runaway growth.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Spitzer IRS mapping of the central kpc of Centaurus A

We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dusty disk and associated with the oval dust shell. The molecular hydrogen transitions indicate the presence of warm gas at temperatures 250--720K. The ratio of the dust emission features at 7.7 and 11.3 micron and the ratio of the [NeII](12.8) and 11.3 dust emission feature are lower in the 500 pc dust shell than in the star forming disk. The clearer shell morphology at 11.3 micron, warm molecular hydrogen emission in the shell, and variation in line ratios in the shell compared to those in the disk, confirm spectroscopically that this shell is a separate coherent entity and is unlikely to be a chance superposition of dust filaments. The physical conditions in the shell are most similar to Galactic supernova remnants where blast waves encounter molecular clouds. The lines requiring the highest level of ionization, [NeV] and [OIV], are detected 20--25 arcsec north-east and south-west of the nucleus and at position angles near the radio jet axis. Fine structure line ratios and limits from this region suggest that the medium is low density and illuminated by a hard radiation field at low ionization parameter. These higher S molecular hydrogen pure rotational transitions are also particularly bright in the same region as the [OIV] and [NeV] emission. This suggests that the gas associated with the dust shell has been excited near the jet axis and is part of an ionization cone.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Chaotic zone boundary for low free eccentricity particles near an eccentric planet

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Via collisionless particle integration we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semi-major axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7 law has been previously explained by estimating the semi-major at which the first order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first order corotation resonances that have strength proportional to the planet's eccentricity. However, we find the corotation resonance width at low free eccentricity is small. Also the first order resonance width at zero free eccentricity is the same as that for a zero eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Discovery of a 500 pc shell in the nucleus of Centaurus A

Spitzer Space Telescope mid-infrared images of the radio galaxy Centaurus A reveal a shell-like, bipolar, structure 500 pc to the north and south of the nucleus. This shell is seen in 5.8, 8.0 and 24 micron broad-band images. Such a remarkable shell has not been previously detected in a radio galaxy and is the first extragalactic nuclear shell detected at mid-infrared wavelengths. We estimate that the shell is a few million years old and has a mass of order million solar masses. A conservative estimate for the mechanical energy in the wind driven bubble is 10^53 erg. The shell could have created by a small few thousand solar mass nuclear burst of star formation. Alternatively, the bolometric luminosity of the active nucleus is sufficiently large that it could power the shell. Constraints on the shell's velocity are lacking. However, if the shell is moving at 1000 km/s then the required mechanical energy would be 100 times larger.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Predictions for a planet just inside Fomalhaut's eccentric ring

We propose that the eccentricity and sharpness of the edge of Fomalhaut's disk are due to a planet just interior to the ring edge. The collision timescale consistent with the disk opacity is long enough that spiral density waves cannot be driven near the planet. The ring edge is likely to be located at the boundary of a chaotic zone in the corotation region of the planet. We find that this zone can open a gap in a particle disk as long as the collision timescale exceeds the removal or ejection timescale in the zone. We use the slope measured from the ring edge surface brightness profile to place an upper limit on the planet mass. The removal timescale in the chaotic zone is used to estimate a lower limit. The ring edge has eccentricity caused by secular perturbations from the planet. These arguments imply that the planet has a mass between that of Neptune and that of Saturn, a semi-major axis of approximately 119 AU and longitude of periastron and eccentricity, 0.1, the same as that of the ring edge.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

The decay of interplanetary coronal mass ejections and forbush decrease recovery times

We investigate the relation between Forbush cosmic ray decrease recovery time and coronal mass ejection transit time between the Sun and Earth. We identify 17 Forbush decreases from ground based neutron count rates between 1978 and 2003 that occur during the same phase in the solar cycle and can be associated with single coronal mass ejections (CMEs) in the SOHO LASCO CME Catalog or previously published reports, and with specific interplanetary coronal mass ejections (ICMEs) crossing the vicinity of Earth. We find an anti-correlation between Forbush recovery times and CME transit time that contradicts the predictions of simple cosmic ray diffusive barrier models. The anti-correlation suggests that the decay rate of ICMEs is anti-correlated with their travel speed. Forbush recovery times range from seven times the transit time for the fastest disturbance to a fifth the Sun-Earth transit time for the slowest. To account for the large range of measured recovery times we infer that the slowest disturbances must decay rapidly with radius whereas the fastest ones must remain strong. The longest recovery times suggest that the fastest disturbances in our sample decayed less rapidly with radius than the ambient solar wind magnetic field strength.

0 citations0 reviewsNo code linkedOpen access