Showing 1–50 of 158 results for author: Quillen, A C

The Limited Role of the Streaming Instability During Moon and Exomoon Formation

Authors: Miki Nakajima, Jeremy Atkins, Jacob B. Simon, Alice C. Quillen

Abstract: It is generally accepted that the Moon accreted from the disk formed by an impact between the proto-Earth and impactor, but its details are highly debated. Some models suggest that a Mars-sized impactor formed a silicate melt-rich (vapor-poor) disk around Earth, whereas other models suggest that a highly energetic impact produced a silicate vapor-rich disk. Such a vapor-rich disk, however, may not… ▽ More It is generally accepted that the Moon accreted from the disk formed by an impact between the proto-Earth and impactor, but its details are highly debated. Some models suggest that a Mars-sized impactor formed a silicate melt-rich (vapor-poor) disk around Earth, whereas other models suggest that a highly energetic impact produced a silicate vapor-rich disk. Such a vapor-rich disk, however, may not be suitable for the Moon formation, because moonlets, building blocks of the Moon, of 100 m-100 km may experience strong gas drag and fall onto Earth on a short timescale, failing to grow further. This problem may be avoided if large moonlets ($\gg 100$ km) form very quickly by streaming instability, which is a process to concentrate particles enough to cause gravitational collapse and rapid formation of planetesimals or moonlets. Here, we investigate the effect of the streaming instability in the Moon-forming disk for the first time and find that this instability can quickly form $\sim 100$ km-sized moonlets. However, these moonlets are not large enough to avoid strong drag and they still fall onto Earth quickly. This suggests that the vapor-rich disks may not form the large Moon, and therefore the models that produce vapor-poor disks are supported. This result is applicable to general impact-induced moon-forming disks, supporting the previous suggestion that small planets ($<1.6 R_\oplus$) are good candidates to host large moons because their impact-induced disks would be likely vapor-poor. We find a limited role of streaming instability in a satellite formation in an impact-induced disk, whereas it plays a key role during planet formation. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: Accepted for Publication in PSJ

Journal ref: 2024, The Planetary Science Journal, 5, 145

The azimuthal distribution of ejecta mass from oblique impacts into sand

Authors: Alice C. Quillen, Sean Doran

Abstract: We measure ejecta mass as a function of azimuthal and impact angle for 104 m/s oblique impacts into sand. We find that the ejecta mass distribution is strongly sensitive to azimuthal angle with up to 8 times more mass in ejecta on the downrange side compared to the uprange side. Crater radii, measured from the site of impact, are measured at different impact and azimuthal angles. Crater ejecta sca… ▽ More We measure ejecta mass as a function of azimuthal and impact angle for 104 m/s oblique impacts into sand. We find that the ejecta mass distribution is strongly sensitive to azimuthal angle with up to 8 times more mass in ejecta on the downrange side compared to the uprange side. Crater radii, measured from the site of impact, are measured at different impact and azimuthal angles. Crater ejecta scaling laws are modified to depend on azimuthal and impact angle. We find that crater radii are sensitive to both impact and azimuthal angle but the ejecta mass as a function of both angles can be estimated from the cube of the crater radius without an additional angular dependent function. The ejecta distributions are relevant for processes that depend upon the integrated properties of intermediate velocity impacts occurring in the outer solar system and possibly during planetesimal formation. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Galactic Structure From Binary Pulsar Accelerations: Beyond Smooth Models

Authors: Thomas Donlon II, Sukanya Chakrabarti, Lawrence M. Widrow, Michael T. Lam, Philip Chang, Alice C. Quillen

Abstract: We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρろー_0=0.062\pm0.017$ \m… ▽ More We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρろー_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($ρろー_{0,\textrm{DM}}=-0.010\pm0.018$ \msun/pc$^3$); these values are similar to, but have smaller uncertainties than previous pulsar timing measurements of these quantities. Here, we provide for the first time, values for the Oort constants and the slope of the rotation curve from direct acceleration measurements. We find that $A=15.4\pm2.6$ km/s/kpc and $B=-13.1\pm2.6$ km/s/kpc (consistent with results from \textit{Gaia}), and the slope of the rotation curve near the Sun is $-2\pm5$ km/s/kpc. We show that the Galactic acceleration field is clearly asymmetric, but due to data limitations it is not yet clear which physical processes drive this asymmetry. We provide updated models of the Galactic potential that account for various sources of disequilibrium; these models are incompatible with commonly used kinematic potentials. This indicates that use of kinematically derived Galactic potentials in precision tests (e.g., in tests of general relativity with pulsar timing) may be subject to larger uncertainties than reported. The acceleration data indicates that the mass of the Galaxy within the Solar circle is $2.3 \times 10^{11}$ M$_\odot$, roughly twice as large as currently accepted models. Additionally, the residuals of the acceleration data compared to existing Galactic models have a dependence on radial position; this trend can be explained if the Sun has an additional acceleration away from the Galactic center. △ Less

Submitted 20 May, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

Gaia DR3 data consistent with a short bar connected to a spiral arm

Authors: E. Vislosky, I. Minchev, S. Khoperskov, M. Martig, T. Buck, T. Hilmi, B. Ratcliffe, J. Bland-Hawthorn, A. C. Quillen, M. Steinmetz, R. de Jong

Abstract: We use numerical simulations to model Gaia DR3 data with the aim of constraining the Milky Way bar and spiral structure parameters. We show that both the morphology and the velocity field in Milky Way-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, v_R(x,y), we ca… ▽ More We use numerical simulations to model Gaia DR3 data with the aim of constraining the Milky Way bar and spiral structure parameters. We show that both the morphology and the velocity field in Milky Way-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, v_R(x,y), we can constrain the bar-spiral orientation. Incorporating uncertainties into our models is necessary to match the data; most importantly, a heliocentric distance uncertainty above 10-15% distorts the bar's shape and v_R quadrupole pattern morphology, and decreases its apparent angle with respect to the Sun-Galactocentric line. An excellent match to the Gaia DR3 v_R(x,y) field is found for a simulation with a bar length R_b\approx3.6 kpc. We argue that the data are consistent with a MW bar as short as ~3 kpc, for moderate strength inner disc spiral structure (A_2/A_0\approx0.25) or, alternatively, with a bar length up to ~5.2 kpc, provided that spiral arms are quite weak (A_2/A_0\approx0.1), and is most likely in the process of disconnecting from a spiral arm. We demonstrate that the bar angle and distance uncertainty can similarly affect the match between our models and the data - a smaller bar angle (20 deg instead of 30 deg) requires smaller distance uncertainty (20% instead of 30%) to explain the observations. Fourier components of the face-on density distribution of our models suggest that the MW does not have strong m=1 and/or m=3 spirals near the solar radius. △ Less

Submitted 4 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 17p, 12 figs. Accepted to MNRAS on Dec 22, 2023

Ncorpi$\mathcal{O}$N : A $\mathcal{O}(N)$ software for N-body integration in collisional and fragmenting systems

Authors: Jérémy Couturier, Alice C. Quillen, Miki Nakajima

Abstract: Ncorpi$\mathcal{O}$N is a $N$-body software developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of Ncorpi$\mathcal{O}$N can choose between four different built-in modules to comp… ▽ More Ncorpi$\mathcal{O}$N is a $N$-body software developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of Ncorpi$\mathcal{O}$N can choose between four different built-in modules to compute self-gravity and detect collisions. One of these makes use of a mesh-based algorithm to treat mutual interactions in $\mathcal{O}(N)$ time. Another module, much more efficient than the standard Barnes-Hut tree code, is a $\mathcal{O}(N)$ tree-based algorithm called FalcON. It relies on fast multipole expansion for gravity computation and we adapted it to collision detection as well. Computation time is reduced by building the tree structure using a three-dimensional Hilbert curve. For the same precision in mutual gravity computation, Ncorpi$\mathcal{O}$N is found to be up to 25 times faster than the famous software REBOUND. Ncorpi$\mathcal{O}$N is written entirely in the C language and only needs a C compiler to run. A python add-on, that requires only basic python libraries, produces animations of the simulations from the output files. The name Ncorpi$\mathcal{O}$N, reminding of a scorpion, comes from the French $N$-corps, meaning $N$-body, and from the mathematical notation $\mathcal{O}(N)$, due to the running time of the software being almost linear in the total number $N$ of moonlets. Ncorpi$\mathcal{O}$N is designed for the study of accreting or fragmenting disks of planetesimal or moonlets. It detects collisions and computes mutual gravity faster than REBOUND, and unlike other $N$-body integrators, it can resolve a collision by fragmentation. The fast multipole expansions are implemented up to order six to allow for a high precision in mutual gravity computation. △ Less

Submitted 26 March, 2024; v1 submitted 31 October, 2023; originally announced October 2023.

Comments: 21 pages, 6 figures, 2 schemas

MSC Class: 70F10

Wind erosion and transport on planetesimals

Authors: Alice C. Quillen, Stephen Luniewski, Adam E. Rubinstein, Jeremy Couturier, Rachel Glade, Miki Nakajima

Abstract: We consider the possibility that aeolian (wind blown) processes occur on small, 1 to 100~km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AUえーゆー), and micron sized particles off the surf… ▽ More We consider the possibility that aeolian (wind blown) processes occur on small, 1 to 100~km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AUえーゆー), and micron sized particles off the surface of a 10 km diameter object in the Transneptunian region. The headwind is sufficiently strong to overcome surface cohesion in the inner solar system, but not in the outer solar system. However, in the outer solar system, surface particles can be redistributed or escape due to impacts from particles that are in the protosolar disk's wind. Based on scaling crater ejecta, we estimate that impacts from particles in the headwind will lead to erosion of mass rather than accretion for planetesimals below about 6 km in diameter. The erosion limit is independent of material strength but proportional to the wind velocity. We explore the sensitivity of splash particle trajectories to particle size, headwind velocity and Reynolds number. Winds from a protostellar disk could account for Kuiper Belt Object (486958) Arrokoth's smooth undulating terrain but only during an epoch of high particle flux and low wind velocity. These conditions could have been present during and just after coalescence of Arrokoth's building blocks. △ Less

Submitted 10 January, 2024; v1 submitted 14 October, 2023; originally announced October 2023.

Rearrangement of Granular Surfaces on Asteroids due to Thermal Cycling

Authors: Danielle Bovie, A. C. Quillen, Rachel Glade

Abstract: In granular systems, thermal cycling causes compaction, creep, penetration of dense objects, and ratcheting of grains against each other. On asteroid surfaces, thermal cycling is high amplitude and can happen billions of times in a few million years. We use a 1-dimensional thermophysical conductivity model to estimate the relative displacement of grains in proximity to one another, caused by varia… ▽ More In granular systems, thermal cycling causes compaction, creep, penetration of dense objects, and ratcheting of grains against each other. On asteroid surfaces, thermal cycling is high amplitude and can happen billions of times in a few million years. We use a 1-dimensional thermophysical conductivity model to estimate the relative displacement of grains in proximity to one another, caused by variations in thermal conductivity or shadows. We find that grains would experience relative displacements of order a few microns during each thermal cycle. If thermal cycling causes diffusive transport, then the asteroid's few centimeters deep thermal skin depth could flow a few centimeters in a million years. Thermal cycling could cause long-distance flows on sloped surfaces, allowing fine materials to collect in depressions. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Subsurface pulse, crater and ejecta asymmetry from oblique impacts into granular media

Authors: Bingcheng Suo, A. C. Quillen, Max Neiderbach, Luke O'Brient, Abobakar Sediq Miakhel, Nathan Skerrett, Jérémy Couturier, Victor Lherm, Jiaxin Wang, Hesam Askari, Esteban Wright, Paul Sánchez

Abstract: We carry out experiments of 104 m/s velocity oblique impacts into a granular medium (sand). Impact craters have nearly round rims even at a grazing angle of about $10^\circ$, however, the strength of seismic pulses excited by the impact is dependent upon impact angle, and the ratio between uprange and downrange velocity peaks can be as large as 5, particularly at shallow depths. Crater slope, an o… ▽ More We carry out experiments of 104 m/s velocity oblique impacts into a granular medium (sand). Impact craters have nearly round rims even at a grazing angle of about $10^\circ$, however, the strength of seismic pulses excited by the impact is dependent upon impact angle, and the ratio between uprange and downrange velocity peaks can be as large as 5, particularly at shallow depths. Crater slope, an offset between crater center and impact site, crater volume, azimuthal variation in ejection angle, seismic pulse shapes and subsurface flow direction are also sensitive to impact angle, but to a much lower degree than subsurface pulse strength. Uprange and downrange pulse peak amplitudes can be estimated from the horizontal and vertical components of the momentum imparted to the medium from the projectile △ Less

Submitted 22 September, 2023; v1 submitted 3 August, 2023; originally announced August 2023.

Comments: 28 pages, 22 figures

HOPS 361-C's Jet Decelerating and Precessing Through NGC 2071 IR

Authors: Adam E. Rubinstein, Nicole Karnath, Alice C. Quillen, Samuel Federman, Joel D. Green, Edward T. Chambers, Dan M. Watson, S. Thomas Megeath

Abstract: We present a two-epoch Hubble Space Telescope (HST) study of NGC 2071 IR highlighting HOPS 361-C, a protostar producing an arced 0.2 parsec-scale jet. Proper motions for the brightest knots decrease from 350 to 100 km/s with increasing distance from the source. The [Fe II] and Pa$βべーた$ emission line intensity ratio gives a velocity jump through each knot of 40--50 km/s. A new [O I] 63 \mic\ spectrum,… ▽ More We present a two-epoch Hubble Space Telescope (HST) study of NGC 2071 IR highlighting HOPS 361-C, a protostar producing an arced 0.2 parsec-scale jet. Proper motions for the brightest knots decrease from 350 to 100 km/s with increasing distance from the source. The [Fe II] and Pa$βべーた$ emission line intensity ratio gives a velocity jump through each knot of 40--50 km/s. A new [O I] 63 \mic\ spectrum, taken with the German REciever for Astronomy at Terahertz frequencies (GREAT) instrument aboard Stratospheric Observatory for Infrared Astronomy (SOFIA), shows a low line-of-sight velocity indicative of high jet inclination. Proper motions and jump velocities then estimate 3D flow speed for knots. Subsequently, we model knot positions and speeds with a precessing jet that decelerates. Measurements are matched with a precession period of 1,000--3,000 years and half opening angle of $15^\circ$. The [Fe II] 1.26-to-1.64 \mic\ line intensity ratio determines visual extinction to each knot from 5--30 mag. Relative to $\sim$14 mag of extinction through the cloud from $\rm{C}^{18}$O emission maps, the jet is embedded at a 1/5 to 4/5 fractional cloud depth. Our model suggests the jet is dissipated over a 0.2 pc arc. This short distance may result from the jet sweeping through a wide angle, allowing the cloud time to fill cavities opened by the jet. Precessing jets contrast with nearly unidirectional protostellar jets that puncture host clouds and can propagate significantly further. △ Less

Submitted 24 March, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 25 pages, 10 figures, submitted to ApJ on 22-11-2022, revised until 8-3-2023, accepted 9-3-2023

Transmission of a Seismic Wave generated by impacts on Granular Asteroids

Authors: Paul Sánchez, Daniel J. Scheeres, Alice C. Quillen

Abstract: In this paper we use a Soft-Sphere Discrete Element method code to simulate the transmission and study the attenuation of a seismic wave. Then, we apply our findings to the different space missions that have had to touch the surface of different small bodies. Additionally, we do the same in regards to the seismic wave generated by the hypervelocity impacts produced by the DART and Hayabusa2 missio… ▽ More In this paper we use a Soft-Sphere Discrete Element method code to simulate the transmission and study the attenuation of a seismic wave. Then, we apply our findings to the different space missions that have had to touch the surface of different small bodies. Additionally, we do the same in regards to the seismic wave generated by the hypervelocity impacts produced by the DART and Hayabusa2 missions once the shock wave transforms into a seismic wave. We find that even at very low pressures, such as those present in the interior of asteroids, the seismic wave speed can still be on the order of hundreds of m/s depending on the velocity of the impact that produces the wave. As expected from experimental measurements, our results show that wave velocity is directly dependent on $P^{1/6}$, where $P$ is the total pressure (confining pressure plus wave induced pressure). Regardless of the pressure of the system and the velocity of the impact (in the investigated range), energy dissipation is extremely high. These results provide us with a way to anticipate the extent to which a seismic wave could have been capable of moving some small particles on the surface of a small body upon contact with a spacecraft. Additionally, this rapid energy dissipation would imply that even hypervelocity impacts should perturb only the external layer of a self-gravitating aggregate on which segregation and other phenomena could take place. This would in turn produce a layered structure of which some evidence has been observed △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: Accepted for publication in The Planetary Sciences Journal

Surface particle motions excited by a low velocity normal impact into a granular medium

Authors: Max Neiderbach, Bingcheng Suo, Esteban Wright, A. C. Quillen, Mokin Lee, Peter Miklavcic, Hesam Askari, Paul Sánchez

Abstract: In laboratory experiments, high speed videos are used to detect and track mm-size surface particle motions caused by a low velocity normal impact into sand. Outside the final crater radius and prior to the landing of the ejecta curtain, particle displacements are measured via particle tracking velocimetry and with a cross-correlation method. Surface particles rebound and are also permanently displ… ▽ More In laboratory experiments, high speed videos are used to detect and track mm-size surface particle motions caused by a low velocity normal impact into sand. Outside the final crater radius and prior to the landing of the ejecta curtain, particle displacements are measured via particle tracking velocimetry and with a cross-correlation method. Surface particles rebound and are also permanently displaced with both peak and permanent displacements rapidly decaying as a function of distance from the crater center. The surface begins to move before most of the ejecta curtain has landed, but continues to move after the subsurface seismic pulse has decayed. Ray angles for surface and subsurface velocities are similar to those described by a Maxwell's Z-model. This implies that the flow field outside the crater excavation region is a continuation of the crater excavation flow. The ratio of final particle displacement to crater radius resembles that measured for other impact craters. △ Less

Submitted 5 October, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

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

Authors: Peter M. Miklavčič, Hesam Askari, Paul Sánchez, Alice C. Quillen, Esteban Wright

Abstract: 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 experi… ▽ More 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. △ Less

Submitted 25 May, 2022; v1 submitted 7 January, 2022; originally announced January 2022.

Propagation and attenuation of pulses driven by low velocity normal impacts in granular media

Authors: A. C. Quillen, Max Neiderbach, Bingcheng Suo, Juliana South, Esteban Wright, Nathan Skerrett, Paul Sánchez, Fernando David Cúñez, Peter Miklavcic, Heesam Askari

Abstract: We carry out experiments of low velocity normal impacts into granular materials that fill an approximately cylindrical 42 litre tub. Motions in the granular medium are tracked with an array of 7 embedded accelerometers. Longitudinal pulses excited by the impact attenuate and their shapes broaden and become smoother as a function of travel distance from the site of impact. Pulse propagation is not… ▽ More We carry out experiments of low velocity normal impacts into granular materials that fill an approximately cylindrical 42 litre tub. Motions in the granular medium are tracked with an array of 7 embedded accelerometers. Longitudinal pulses excited by the impact attenuate and their shapes broaden and become smoother as a function of travel distance from the site of impact. Pulse propagation is not spherically symmetric about the site of impact. Peak amplitudes are about twice as large for the pulse propagating downward than at 45 degrees from vertical. An advection-diffusion model is used to estimate the dependence of pulse properties as a function of travel distance from the site of impact. The power law forms for pulse peak pressure, velocity and seismic energy depend on distance from impact to a power of -2.5 and this rapid decay is approximately consistent with our experimental measurements. Our experiments support a seismic jolt model, giving rapid attenuation of impact generated seismic energy into rubble asteroids, rather than a reverberation model, where seismic energy slowly decays. We apply our diffusive model to estimate physical properties of the seismic pulse that will be excited by the forthcoming DART mission impact onto the secondary, Dimorphos, of the asteroid binary (65803) Didymos system. We estimate that the pulse peak acceleration will exceed the surface gravity as it travels through the asteroid. △ Less

Submitted 1 August, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

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

Authors: Esteban Wright, Alice C. Quillen, Paul Sanchez, Stephen R. Schwartz, Miki Nakajima, Hesam Askari, Peter Miklavcic

Abstract: 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- t… ▽ More 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. △ Less

Submitted 22 January, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Link to the high speed videos for the 4 experiments tracked and discussed: https://youtube.com/playlist?list=PLYnpHw29go3pCWwm554fwR6jNFvB5qVQn Submitted to the journal Icarus on September 12, 2021

Non-principal axis rotation in binary asteroid systems and how it weakens the BYORP effect

Authors: Alice C. Quillen, Anthony LaBarca, YuanYuan Chen

Abstract: Using viscoelastic mass/spring model simulations, we explore tidal evolution and migration of compact binary asteroid systems. We find that after the secondary is captured into a spin-synchronous state, non-principal axis rotation in the secondary can be long-lived. The secondary's long axis can remain approximately aligned along the vector connecting secondary to primary while the secondary rocks… ▽ More Using viscoelastic mass/spring model simulations, we explore tidal evolution and migration of compact binary asteroid systems. We find that after the secondary is captured into a spin-synchronous state, non-principal axis rotation in the secondary can be long-lived. The secondary's long axis can remain approximately aligned along the vector connecting secondary to primary while the secondary rocks back and forth about its long axis. Inward orbital semi-major axis migration can also resonantly excite non-principal axis rotation. By estimating solar radiation forces on triangular surface meshes, we show that the magnitude of the BYORP effect induced torque is sensitive to the secondary's spin state. Non-principal axis rotation within the 1:1 spin-orbit resonance can reduce the BYORP torque or cause frequent reversals in its direction. △ Less

Submitted 19 November, 2021; v1 submitted 30 July, 2021; originally announced July 2021.

Comments: to be published in Icarus

A measurement of the Galactic plane mass density from binary pulsar accelerations

Authors: Sukanya Chakrabarti, Philip Chang, Michael T. Lam, Sarah J. Vigeland, Alice C. Quillen

Abstract: We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the Solar System barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of… ▽ More We use compiled high-precision pulsar timing measurements to directly measure the Galactic acceleration of binary pulsars relative to the Solar System barycenter. Given the vertical accelerations, we use the Poisson equation to derive the Oort limit, i.e., the total volume mass density in the Galactic mid-plane. Our best-fitting model gives an Oort limit of $0.08^{0.05}_{-0.02} M_{\odot}/\rm pc^{3}$, which is close to estimates from recent Jeans analyses. Given the accounting of the baryon budget from McKee et al. (2015), we obtain a local dark matter density of $-0.004^{0.05}_{-0.02}~M_{\odot}/\rm pc^{3}$, which is slightly below other modern estimates but consistent within the current uncertainties of our method. While this first measurement of the Oort limit (and other Galactic parameters) has error bars that are currently several times larger than kinematical estimates, they should improve in the future. We also constrain the oblateness of the potential, finding it consistent with that expected from the disk and inconsistent with a potential dominated by a spherical halo, as is appropriate for our sample which is within a $\sim$ kpc of the Sun. We find that the slope of the rotation curve is not constrained by current measurements of binary pulsar accelerations. We give a fitting function for the vertical acceleration $a_{z}$: $a_{z} = -αあるふぁ_{1}z$; $\log_{10} (αあるふぁ_{1}/{\rm Gyr}^{-2})=3.69^{0.19}_{-0.12}$. By analyzing interacting simulations of the Milky Way, we find that large asymmetric variations in $da_{z}/dz$ as a function of vertical height may be a signature of sub-structure. We end by discussing the power of combining constraints from pulsar timing and high-precision radial velocity (RV) measurements towards lines-of-sight near pulsars, to test theories of gravity and constrain dark matter sub-structure. △ Less

Submitted 1 January, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: accepted to ApJ Letters, changes following referee comments

Accretion of Ornamental Equatorial Ridges on Pan, Atlas and Daphnis

Authors: Alice C. Quillen, Fatima Zaidouni, Miki Nakajima, Esteban Wright

Abstract: 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 durat… ▽ More 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. △ Less

Submitted 17 January, 2021; v1 submitted 26 July, 2020; originally announced July 2020.

Birth Sites of Young Stellar Associations and Recent Star Formation in a Flocculent Corrugated Disk

Authors: Alice C. Quillen, Alex R. Pettitt, Sukanya Chakrabarti, Yifan Zhang, Jonathan Gagné, Ivan Minchev

Abstract: With backwards orbit integration we estimate birth locations of young stellar associations and moving groups identified in the solar neighborhood that are younger than 70 Myr. The birth locations of most of these stellar associations are at smaller galactocentric radius than the Sun, implying that their stars moved radially outwards after birth. Exceptions to this rule include are the Argus and Oc… ▽ More With backwards orbit integration we estimate birth locations of young stellar associations and moving groups identified in the solar neighborhood that are younger than 70 Myr. The birth locations of most of these stellar associations are at smaller galactocentric radius than the Sun, implying that their stars moved radially outwards after birth. Exceptions to this rule include are the Argus and Octans associations which formed outside the Sun's Galactocentric radius. Variations in birth heights of the stellar associations suggest that they were born in a corrugated disk of molecular clouds, similar to that inferred from the current filamentary molecular cloud distribution and dust extinction maps. Multiple spiral arm features with different but near corotation pattern speeds and at different heights could account for the stellar association birth sites. We find that the young stellar associations are located in between peaks in the UV stellar velocity distribution for stars in the solar neighborhood. This would be expected if they were born in a spiral arm which perturbs stellar orbits that cross it. In contrast, stellar associations seem to be located near peaks in the vertical phase space distribution, suggesting that the gas in which stellar associations are born moves vertically together with the low velocity dispersion disk stars. △ Less

Submitted 3 August, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: The animated movie is available here: http://astro.pas.rochester.edu/~aquillen/figures/movie3D.gif

Fluctuations in galactic bar parameters due to bar-spiral interaction

Authors: T. Hilmi, I. Minchev, T. Buck, M. Martig, A. C. Quillen, G. Monari, B. Famaey, R. S. de Jong, C. F. P. Laporte, J. Read, J. L. Sanders, M. Steinmetz, C. Wegg

Abstract: We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlati… ▽ More We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value. △ Less

Submitted 30 June, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 25 p., 19 figures, Accepted to MNRAS

Excitation of Tumbling in Phobos and Deimos

Authors: Alice C. Quillen, Mckenzie Lane, Miki Nakajima, Esteban Wright

Abstract: 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 t… ▽ More 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. △ Less

Submitted 13 January, 2020; v1 submitted 15 August, 2019; originally announced August 2019.

Star Formation Efficiencies at Giant Molecular Cloud Scales in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

Authors: D. Espada, S. Verley, R. E. Miura, F. P. Israel, C. Henkel, S. Matsushita, B. Vila-Vilaro, J. Ott, K. Morokuma-Matsui, A. B. Peck, A. Hirota, S. Aalto, A. C. Quillen, M. R. Hogerheijde, N. Neumayer, C. Vlahakis, D. Iono, K. Kohno

Abstract: We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previ… ▽ More We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $τたう=$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas. △ Less

Submitted 16 December, 2019; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: Minor update. 20 pages, 11 figures, accepted for publication in ApJ

Simulations of wobble damping in viscoelastic rotators

Authors: Alice C. Quillen, Katelyn J. Wagner, Paul Sanchez

Abstract: Using a damped mass-spring model, we simulate wobble of spinning homogeneous viscoelastic ellipsoids undergoing non-principal axis rotation. Energy damping rates are measured for oblate and prolate bodies with different spin rates, spin states, viscoelastic relaxation timescales, axis ratios, and strengths. Analytical models using a quality factor by Breiter et al. (2012) and for the Maxwell rheol… ▽ More Using a damped mass-spring model, we simulate wobble of spinning homogeneous viscoelastic ellipsoids undergoing non-principal axis rotation. Energy damping rates are measured for oblate and prolate bodies with different spin rates, spin states, viscoelastic relaxation timescales, axis ratios, and strengths. Analytical models using a quality factor by Breiter et al. (2012) and for the Maxwell rheology by Frouard & Efroimsky (2018) match our numerical measurements of the energy dissipation rate after we modify their predictions for the numerically simulated Kelvin-Voigt rheology. Simulations of nearly spherical but wobbling bodies with hard and soft cores show that the energy dissipation rate is more sensitive to the material properties in the core than near the surface. The sensitivity to viscoelastic model implies that inferred statistics of tumbling lifetimes in asteroids might be interpreted in terms of differences in their material properties. △ Less

Submitted 7 February, 2019; v1 submitted 5 January, 2019; originally announced January 2019.

Comments: accepted MNRAS

Boulder Stranding in Ejecta Launched by an Impact Generated Seismic Pulse

Authors: Esteban Wright, Alice C. Quillen, Juliana South, Randal C. Nelson, Paul Sanchez, Larkin Martini, Stephen Schwartz, Miki Nakajima, Erik Asphaug

Abstract: We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions… ▽ More We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa. △ Less

Submitted 13 December, 2019; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: Link to experiment videos: https://www.youtube.com/playlist?list=PLYnpHw29go3rqDIHmaC8BjWKojxr2wSu8

Journal ref: Icarus, Volume 337, 2020, 113424. ISSN 0019-1035

Near/Far Side Asymmetry in the Tidally Heated Moon

Authors: Alice C. Quillen, Larkin Martini, Miki Nakajima

Abstract: Using viscoelastic mass spring model simulations to track heat distribution inside a tidally perturbed body, we measure the asymmetry of heating in the crust of a spin synchronous Moon in eccentric orbit about the Earth. With the Moon within a few Earth radii of the Earth, we find that tidal heating per unit area in a lunar crustal shell is asymmetric due to the octupole order moment in the Earth'… ▽ More Using viscoelastic mass spring model simulations to track heat distribution inside a tidally perturbed body, we measure the asymmetry of heating in the crust of a spin synchronous Moon in eccentric orbit about the Earth. With the Moon within a few Earth radii of the Earth, we find that tidal heating per unit area in a lunar crustal shell is asymmetric due to the octupole order moment in the Earth's tidal field and is 10 to 20\% higher on its near side than on its far side. Tidal heating reduces the crustal basal heat flux and the rate of magma ocean crystallization. Assuming that the local crustal growth rate depends on the local basal heat flux and the distribution of tidal heating in latitude and longitude, a heat conductivity model illustrates that a moderately asymmetric and growing lunar crust could maintain its asymmetry. Near/far side asymmetric tidal heating could serve as a trigger for uneven lunar crustal growth. △ Less

Submitted 10 April, 2019; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: Accepted Icarus

Spiral arm crossings inferred from ridges in Gaia stellar velocity distributions

Authors: Alice C. Quillen, Ismael Carrillo, Friedrich Anders, Paul McMillan, Tariq Hilmi, Giacomo Monari, Ivan Minchev, Cristina Chiappini, Arman Khalatyan, Matthias Steinmetz

Abstract: The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spir… ▽ More The solar neighborhood contains disc stars that have recently crossed spiral arms in the Galaxy. We propose that boundaries in local velocity distributions separate stars that have recently crossed or been perturbed by a particular arm from those that haven't. Ridges in the stellar velocity distributions constructed from the second Gaia data release trace orbits that could have touched nearby spiral arms at apocentre or pericentre. The multiple ridges and arcs seen in local velocity distributions are consistent with the presence of multiple spiral features and different pattern speeds and imply that the outer Galaxy is flocculent rather than grand design. △ Less

Submitted 21 July, 2018; v1 submitted 25 May, 2018; originally announced May 2018.

Comments: For a movie https://youtu.be/5nbtNYDobcE

Coma Berenices: first evidence for incomplete vertical phase-mixing in local velocity space with RAVE - confirmed with Gaia DR2

Authors: G. Monari, B. Famaey, I. Minchev, T. Antoja, O. Bienayme, B. K. Gibson, E. K. Grebel, G. Kordopatis, P. McMillan, J. Navarro, Q. A. Parker, A. C. Quillen, W. Reid, G. Seabroke, A. Siebert, M. Steinmetz, R. F. G. Wyse, T. Zwitter

Abstract: Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mix… ▽ More Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mixing. We estimated for the first time from dynamical arguments that the moving group must have formed at most ~ 1.5 Gyr ago, and related this to a pericentric passage of the Sagittarius dwarf satellite galaxy. The present note is a rewritten version of the original arXiv post on this result now also including a confirmation of our finding with Gaia DR2. △ Less

Submitted 8 May, 2018; v1 submitted 20 April, 2018; originally announced April 2018.

Comments: AAS Research Notes

Impact Excitation of a Seismic Pulse and Vibrational Normal Modes on Asteroid Bennu and Associated Slumping of Regolith

Authors: Alice C. Quillen, Yuhui Zhao, YuanYuan Chen, Paul Sanchez, Randal C. Nelson, Stephen R. Schwartz

Abstract: We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal mo… ▽ More We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal modes. In this regime, the distribution of seismic reverberation is not evenly distributed across the body surface. With mass-spring model elastic simulations, we model impact excitation of seismic waves with a force pulse exerted on the surface and using three different asteroid shape models. The simulations exhibit antipodal focusing and normal mode excitation. If the impulse excited vibrational energy is long lasting, vibrations are highest at impact point, its antipode and at high surface elevations such as an equatorial ridge. A near equatorial impact launches a seismic impulse on a non-spherical body that can be focused on two additional points on an the equatorial ridge. We explore simple flow models for the morphology of vibration induced surface slumping. We find that the initial seismic pulse is unlikely to cause large shape changes. Long lasting seismic reverberation on Bennu caused by a near equatorial impact could have raised the height of its equatorial ridge by a few meters and raised two peaks on it, one near impact site and the other near its antipode. △ Less

Submitted 23 September, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: submitted to Icarus

The GALAH Survey: Stellar streams and how stellar velocity distributions vary with Galactic longitude, hemisphere and metallicity

Authors: Alice C. Quillen, Gayandhi De Silva, Sanjib Sharma, Michael Hayden, Ken Freeman, Joss Bland-Hawthorn, Maruša Žerjal, Martin Asplund, Sven Buder, Valentina D'Orazi, Ly Duong, Janez Kos, Jane Lin, Karin Lind, Sarah Martell, Katharine Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Borja Anguiano, Daniela Carollo, Luca Casagrande, Klemen Cotar, Peter L. Cottrell, Michael Ireland , et al. (8 additional authors not shown)

Abstract: Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to… ▽ More Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy. △ Less

Submitted 4 April, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: accepted for publication in MNRAS

Migration in the shearing sheet and estimates for young open cluster migration

Authors: Alice C. Quillen, Eric Nolting, Ivan Minchev, Gayandhi De Silva, Cristina Chiappini

Abstract: Using tracer particles embedded in self-gravitating shearing sheet N-body simulations, we investigate the distance in guiding centre radius that stars or star clusters can migrate in a few orbital periods. The standard deviations of guiding centre distributions and maximum migration distances depend on the Toomre or critical wavelength and the contrast in mass surface density caused by spiral stru… ▽ More Using tracer particles embedded in self-gravitating shearing sheet N-body simulations, we investigate the distance in guiding centre radius that stars or star clusters can migrate in a few orbital periods. The standard deviations of guiding centre distributions and maximum migration distances depend on the Toomre or critical wavelength and the contrast in mass surface density caused by spiral structure. Comparison between our simulations and estimated guiding radii for a few young super-solar metallicity open clusters, including NGC 6583, suggests that the contrast in mass surface density in the solar neighbourhood has standard deviation (in the surface density distribution) divided by mean of about 1/4 and larger than measured using COBE data by Drimmel and Spergel. Our estimate is consistent with a standard deviation of $\sim$0.07 dex in the metallicities measured from high-quality spectroscopic data for 38 young open clusters (<1 Gyr) with mean galactocentric radius 7-9 kpc. △ Less

Submitted 20 January, 2018; v1 submitted 14 September, 2017; originally announced September 2017.

Comments: accepted to MNRAS

Tilting Styx and Nix but not Uranus with a Spin-Precession-Mean-motion resonance

Authors: Alice C. Quillen, Yuan-Yuan Chen, Benoit Noyelles, Santiago Loane

Abstract: A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about star. We expand the planet-planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration… ▽ More A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about star. We expand the planet-planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration allows the spinning planet to be captured into spin resonance. With initial obliquity near zero, the spin resonance can lift the planet's obliquity to near 90 or 180 degrees depending upon whether the spin resonance is first or zero-th order in inclination. Past capture of Uranus into such a spin resonance could give an alternative non-collisional scenario accounting for Uranus's high obliquity. However we find that the time spent in spin resonance must be so long that this scenario cannot be responsible for Uranus's high obliquity. Our model can be used to study spin resonance in satellite systems. Our Hamiltonian model explains how Styx and Nix can be tilted to high obliquity via outward migration of Charon, a phenomenon previously seen in numerical simulations. △ Less

Submitted 28 September, 2017; v1 submitted 11 July, 2017; originally announced July 2017.

Comments: Submitted to Celestial Mechanics and Dynamical Astronomy

Disentangling the Circumnuclear Environs of Centaurus A: III. An Inner Molecular Ring, Nuclear Shocks and the CO to warm H2 interface

Authors: D. Espada, S. Matsushita, R. E. Miura, F. P. Israel, N. Neumayer, S. Martin, C. Henkel, T. Izumi, D. Iono, S. Aalto, J. Ott, A. B. Peck, A. C. Quillen, K. Kohno

Abstract: We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist with… ▽ More We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist within the CND, which form a ring-like structure with an unprojected diameter of 9 x 6 arcsec (162pc x 108pc) and a position angle PA = 155deg. Inside the nuclear ring, there are two leading and straight filamentary structures with lengths of about 30-60pc at PA = 120deg on opposite sides of the AGN, with a rotational symmetry of 180deg and steeper position-velocity diagrams, which are interpreted as nuclear shocks due to non-circular motions. Along the filaments, and unlike other nearby AGNs, several dense molecular clumps present low HCN/HCO+(4-3) ratios (~0.5). The filaments abruptly end in the probed transitions at r = 20pc from the AGN, but previous near-IR H2 (J=1-0) S(1) maps show that they continue in an even ~1000 K), winding up in the form of nuclear spirals, and forming an inner ring structure with another set of symmetric filaments along the N-S direction and within r = 10pc. The molecular gas is governed primarily by non-circular motions, being the successive shock fronts at different scales where loss of angular momentum occurs, a mechanism which may feed efficiently powerful radio galaxies down to parsec scales. △ Less

Submitted 18 June, 2017; originally announced June 2017.

Comments: 46 pages. Accepted for publication in ApJ

Obliquity evolution of the minor satellites of Pluto and Charon

Authors: Alice C. Quillen, Fiona Nichols-Fleming, Yuan-Yuan Chen, Benoit Noyelles

Abstract: New Horizons mission observations show that the small satellites Styx, Nix, Kerberos and Hydra, of the Pluto-Charon system, have not tidally spun-down to near synchronous spin states and have high obliquities with respect to their orbit about the Pluto-Charon binary (Weaver et al. 2016). We use a damped mass-spring model within an N-body simulation to study spin and obliquity evolution for single… ▽ More New Horizons mission observations show that the small satellites Styx, Nix, Kerberos and Hydra, of the Pluto-Charon system, have not tidally spun-down to near synchronous spin states and have high obliquities with respect to their orbit about the Pluto-Charon binary (Weaver et al. 2016). We use a damped mass-spring model within an N-body simulation to study spin and obliquity evolution for single spinning non-round bodies in circumbinary orbit. Simulations with tidal dissipation alone do not show strong obliquity variations from tidally induced spin-orbit resonance crossing and this we attribute to the high satellite spin rates and low orbital eccentricities. However, a tidally evolving Styx exhibits intermittent obliquity variations and episodes of tumbling. During a previous epoch where Charon migrated away from Pluto, the minor satellites could have been trapped in orbital mean motion inclination resonances. An outward migrating Charon induces large variations in Nix and Styx's obliquities. The cause is a commensurability between the mean motion resonance frequency and the spin precession rate of the spinning body. As the minor satellites are near mean motion resonances, this mechanism could have lifted the obliquities of all four minor satellites. The high obliquities need not be primordial if the minor satellites were at one time captured into mean motion resonances. △ Less

Submitted 20 April, 2017; v1 submitted 19 January, 2017; originally announced January 2017.

Comments: Icarus, in press

Rings beyond the giant planets

Authors: Bruno Sicardy, Maryame El Moutamid, Alice C. Quillen, Paul M. Schenk, Mark R. Showalter, Kevin Walsh

Abstract: Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected… ▽ More Until 2013, only the giant planets were known to host ring systems. In June 2013, a stellar occulation revealed the presence of narrow and dense rings around Chariklo, a small Centaur object that orbits between Saturn and Uranus. Meanwhile, the Cassini spacecraft revealed evidence for the possible past presence of rings around the Saturnian satellites Rhea and Iapetus. Mars and Pluto are expected to have tenuous dusty rings, though they have so far evaded detection. More remotely, transit events observed around a star in 2007 may have revealed for the first time exoplanetary rings around a giant planet orbiting that star. So, evidence is building to show that rings are more common features in the universe than previously thought. Several interesting issues arising from the discovery (or suspicion) of new ring systems are described in this chapter. △ Less

Submitted 11 April, 2017; v1 submitted 10 December, 2016; originally announced December 2016.

Comments: 23 pages, 8 figures, version 2. This material will be published as Chapter 7 of the book Planetary Ring Systems, edited by Matthew S. Tiscareno and Carl D. Murray. See: www.cambridge.org/9781107113824

Tidal spin down rates of homogeneous triaxial viscoelastic bodies

Authors: Alice C. Quillen, Andrea Kueter-Young, Julien Frouard, Darin Ragozzine

Abstract: 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 orb… ▽ More 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. △ Less

Submitted 23 August, 2016; v1 submitted 28 July, 2016; originally announced July 2016.

Comments: accepted for publication in MNRAS

Cold, clumpy accretion onto an active supermassive black hole

Authors: Grant R. Tremblay, J. B. Raymond Oonk, Françoise Combes, Philippe Salomé, Christopher P. O'Dea, Stefi A. Baum, G. Mark Voit, Megan Donahue, Brian R. McNamara, Timothy A. Davis, Michael A. McDonald, Alastair C. Edge, Tracy E. Clarke, Roberto Galván-Madrid, Malcolm N. Bremer, Louise O. V. Edwards, Andrew C. Fabian, Stephen L. Hamer, Yuan Li, Anaëlle Maury, Helen R. Russell, Alice C. Quillen, C. Megan Urry, Jeremy S. Sanders, Michael Wise

Abstract: Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead… ▽ More Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the "hot mode" accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing "shadows" cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it. △ Less

Submitted 10 June, 2016; v1 submitted 7 June, 2016; originally announced June 2016.

Comments: Published in the June 9th, 2016 issue of Nature

Journal ref: Nature, 534, 218-221 (2016)

Dippers and Dusty Disks Edges: A Unified Model

Authors: Eva H. L. Bodman, Alice C. Quillen, Megan Ansdell, Michael Hippke, Tabetha S. Boyajian, Eric E. Mamajek, Eric G. Blackman, Aaron Rizzuto, Joel H. Kastner

Abstract: We revisit the nature of large dips in flux from extinction by dusty circumstellar material that is observed by Kepler for many young stars in the Upper Sco and $ρろー$ Oph star formation regions. These young, low-mass "dipper" stars are known to have low accretion rates and primarily host moderately evolved dusty circumstellar disks. Young low mass stars often exhibit rotating star spots that cause q… ▽ More We revisit the nature of large dips in flux from extinction by dusty circumstellar material that is observed by Kepler for many young stars in the Upper Sco and $ρろー$ Oph star formation regions. These young, low-mass "dipper" stars are known to have low accretion rates and primarily host moderately evolved dusty circumstellar disks. Young low mass stars often exhibit rotating star spots that cause quasi-periodic photometric variations. We found no evidence for periods associated with the dips that are different from the star spot rotation period in spectrograms constructed from the light curves. The material causing the dips in most of these light curves must be approximately corotating with the star. We find that disk temperatures computed at the disk corotation radius are cool enough that dust should not sublime. Crude estimates for stellar magnetic field strengths and accretion rates are consistent with magnetospheric truncation near the corotation radius. Magnetospheric truncation models can explain why the dips are associated with material near corotation and how dusty material is lifted out of the midplane to obscure the star which would account for the large fraction of young low mass stars that are dippers. We propose that variations in disk orientation angle, stellar magnetic field dipole tilt axis, and disk accretion rate are underlying parameters accounting for differences in the dipper light curves. △ Less

Submitted 21 February, 2017; v1 submitted 12 May, 2016; originally announced May 2016.

Comments: submitted to MNRAS, updated for referee revisions

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

Authors: Bryan Hilbert, Marco Chiaberge, JohnPaul Kotyla, Grant R. Tremblay, Carlo Stanghellini, William B. Sparks, Stefi A. Baum, Alessandro Capetti, F. Duccio Macchetto, George K. Miley, Christopher P. O'Dea, Eric S. Perlman, Alice C. Quillen

Abstract: 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… ▽ More 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. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 25 pages, 54 figures; Accepted for publication in ApJS

Hot planetary winds near a star: dynamics, wind-wind interactions, and observational signatures

Authors: Jonathan Carroll-Nellenback, Adam Frank, Baowei Liu, Alice C. Quillen, Eric G. Blackman, Ian Dobbs-Dixon

Abstract: Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind l… ▽ More Signatures of "evaporative" winds from exo-planets on short (hot) orbits around their host star have been observed in a number of systems. In this paper we present global AMR simulations that track the launching of the winds, their expansion through the circumstellar environment, and their interaction with a stellar wind. We focus on purely hydrodynamic flows including the anisotropy of the wind launching and explore the orbital/fluid dynamics of the resulting flows in detail. In particular we find that a combination of the tidal and Coriolis forces strongly distorts the planetary "Parker" wind creating "up-orbit" and "down-orbit" streams. We characterize the flows in terms of their orbital elements which change depending on their launch position on the planet. We find that the anisotropy in the atmospheric temperature leads to significant backflow on to the planet. The planetary wind interacts strongly with the stellar wind creating instabilities that cause eventual deposition of planetary gas onto the star. We present synthetic observations of both transit and absorption line-structure for our simulations. For our initial conditions, we find that the orbiting wind material produces absorption signatures at significant distances from the planet and substantial orbit to orbit variability. Ly-αあるふぁ absorption shows red and blueshifted features out to 70 km/s. Finally, using semi-analytic models we constrain the effect of radiation pressure, given the approximation of uniform stellar absorption. △ Less

Submitted 28 April, 2016; originally announced April 2016.

Comments: 17 pages, 14 figures

MSC Class: 85A20

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

Authors: Julien Frouard, Alice C. Quillen, Michael Efroimsky, David Giannella

Abstract: 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 th… ▽ More 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. △ Less

Submitted 3 March, 2016; v1 submitted 29 January, 2016; originally announced January 2016.

Comments: 14 pages, 4 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 458, Issue 3, p.2890-2901, 2016

Identification of Globular Cluster Stars in RAVE data II: Extended tidal debris around NGC 3201

Authors: B. Anguiano, G. M. De Silva, K. Freeman, G. S. Da Costa, T. Zwitter, A. C. Quillen, D. B. Zucker, J. F. Navarro, A. Kunder, A. Siebert, R. F. G. Wyse, E. K. Grebel, G. Kordopatis, B. K. Gibson, G. Seabroke, S. Sharma, J. Wojno, J. Bland-Hawthorn, Q. A. Parker, M. Steinmetz, C. Boeche, G. Gilmore, O. Bienayme, W. Reid, F. Watson

Abstract: We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Thre… ▽ More We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution. △ Less

Submitted 10 January, 2016; originally announced January 2016.

Comments: Accepted 2016 January 8. Received 2015 December 10; in original form 2015 November 13 MNRAS

Crustal Failure on Icy Moons from a Strong Tidal Encounter

Authors: Alice C. Quillen, David Giannella, John G. Shaw, Cynthia Ebinger

Abstract: 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… ▽ More 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. △ Less

Submitted 5 April, 2016; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: accepted for publication in Icarus, for a movie see https://youtu.be/W5GyTEyXTss

Planetary Evaporation and the Dynamics of Planet Wind/Stellar Wind Bow Shocks

Authors: A. Frank, B. Liu, J. Carroll-Nellenback, A. C. Quillen, J. F. Kasting, I. Dobbs-Dixon, E. G. Blackman

Abstract: We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a stren… ▽ More We present initial results of a new campaign of simulations focusing on the interaction of planetary winds with stellar environments using Adaptive Mesh Refinement methods. We have confirmed the results of Stone & Proga 2009 that an azimuthal flow structure is created in the planetary wind due to day/night temperatures differences. We show that a backflow towards the planet will occur with a strength that depends on the escape parameter. When a stellar outflow is included, we see unstable bow waves forming through the outflow's interaction with the planetary wind. △ Less

Submitted 4 August, 2015; originally announced August 2015.

Comments: To appear in the proceedings of IAU Symposium 314 "Young Stars and Planets Near the Sun"

Phase Wrapping of Epicyclic Perturbations in the Wobbly Galaxy

Authors: Alexander de la Vega, Alice C. Quillen, Jeffrey L. Carlin, Sukanya Chakrabarti, Elena D'Onghia

Abstract: 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 r… ▽ More 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. △ Less

Submitted 2 September, 2015; v1 submitted 27 July, 2015; originally announced July 2015.

Comments: 15 pages, 8 figures, accepted for publication in MNRAS

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

Authors: Grant R. Tremblay, Christopher P. O'Dea, Stefi A. Baum, Rupal Mittal, Michael McDonald, Françoise Combes, Yuan Li, Brian McNamara, Malcolm N. Bremer, Tracy E. Clarke, Megan Donahue, Alastair C. Edge, Andrew C. Fabian, Stephen L. Hamer, Michael T. Hogan, Raymond Oonk, Alice C. Quillen, Jeremy S. Sanders, Philippe Salomé, G. Mark Voit

Abstract: 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… ▽ More 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. △ Less

Submitted 13 May, 2015; originally announced May 2015.

Comments: 36 pages, 30 figures, accepted for publication in MNRAS

Torque on an exoplanet from an anisotropic evaporative wind

Authors: Jean Teyssandier, James E. Owen, Fred C. Adams, Alice C. Quillen

Abstract: 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… ▽ More 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. △ Less

Submitted 7 April, 2015; originally announced April 2015.

Comments: 12 pages, 7 figures, submitted to MNRAS

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

Authors: Alice C. Quillen, Borja Anguiano, Gayandhi De Silva, Ken Freeman, Dan B. Zucker, Ivan Minchev, Joss Bland-Hawthorn

Abstract: 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 neighborhoo… ▽ More 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. △ Less

Submitted 9 April, 2015; v1 submitted 10 March, 2015; originally announced March 2015.

Comments: accepted for publication in MNRAS

Journal ref: 2015, MNRAS, 450, 2354

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

Authors: Erin L. Scott, Eric E. Mamajek, Mark J. Pecaut, Alice C. Quillen, Fred Moolekamp, Cameron P. M. Bell

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: 7 pages, 3 figures, 2 tables; accepted for publication in the Astrophysical Journal

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

Authors: Alice C. Quillen, Robert S. French

Abstract: 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 tw… ▽ More 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. △ Less

Submitted 5 August, 2014; originally announced August 2014.

Journal ref: 2014, MNRAS, 445, 3959

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

Authors: Alice C. Quillen, Marco Ciocca, Jeffrey L. Carlin, Zeyang Meng, Cameron P. M. Bell

Abstract: 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 vari… ▽ More 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. △ Less

Submitted 5 February, 2014; originally announced February 2014.

A search for eclipsing binaries that host discs

Authors: Zeyang Meng, Alice C. Quillen, Cameron P. M. Bell, Eric E. Mamajek, Erin L. Scott, Ji-Lin Zhou

Abstract: 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 a… ▽ More 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. △ Less

Submitted 30 April, 2014; v1 submitted 6 January, 2014; originally announced January 2014.

Comments: 10 pages, 11 figures, 2 tables, accepted for publication in MNRAS

Journal ref: 2014, MNRAS, 441, 3733