Search

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., Comment: Accepted for Publication in PSJ

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.

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 ($\rho_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($\rho_{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.

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., Comment: 17p, 12 figs. Accepted to MNRAS on Dec 22, 2023

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., Comment: 21 pages, 6 figures, 2 schemas

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.

Using the recent ability of quantum computers to initialize quantum states rapidly with high fidelity, we use a function operating on a discrete set to create a simple class of quantum channels. Fixed points and periodic orbits, that are present in the function, generate fixed points and periodic orbits in the associated quantum channel. Phenomenology such as periodic doubling is visible in a 6 qubit dephasing channel constructed from a truncated version of the logistic map. Using disjoint subsets, discrete function-generated channels can be constructed that preserve coherence within subspaces. Error correction procedures can be in this class as syndrome detection uses an initialized quantum register. A possible application for function-generated channels is in hybrid classical/quantum algorithms. We illustrate how these channels can aid in carrying out classical computations involving iteration of non-invertible functions on a quantum computer with the Euclidean algorithm for finding the greatest common divisor of two integers.

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.

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, Comment: 28 pages, 22 figures

The interactions of solid objects with granular media is countered by a resistance force that stems from frictional forces between the grains and the media's resistance to inertia imposed by the intruder. Earlier theories of granular intrusion have suggested an additive contribution of these two families of forces and had tremendous success in predicting resistive forces on arbitrary shaped objects. However, it remains unclear how these forces are influenced by gravitational conditions. We examine the role of gravity on surface impact behavior into cohesionless granular media using hundreds of soft-sphere discrete element simulations, we demonstrate that the outcome of impacts remain qualitatively similar under varying gravitational conditions if initial velocities are scaled with the Froude number, suggesting an underlying law. Using theoretical arguments, we provide reasoning for the observed universality and show that there is a hidden dependency in resistive forces into granular media on Froude number. Following the theoretical framework, we show that Froude number scaling precisely collapses impact trajectories across gravitational conditions, setting the foundation for explorations in granular behavior beyond Earth.

Biological systems can rely on collective formation of a metachronal wave in an ensemble of oscillators for locomotion and for fluid transport. We consider one-dimensional chains of phase oscillators with nearest neighbor interactions, connected in a loop and with rotational symmetry, so each oscillator resembles every other oscillator in the chain. Numerical integrations of the discrete phase oscillator systems and a continuum approximation show that directional models (those that do not obey reversal symmetry), can exhibit instability to short wavelength perturbations but only in regions where the slope in phase has a particular sign. This causes short wavelength perturbations to develop that can vary the winding number that describes the sum of phase differences across the loop and the resulting metachronal wave speed. Numerical integrations of stochastic directional phase oscillator models show that even a weak level of noise can seed instabilities that resolve into metachronal wave states.

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$\beta$ 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., Comment: 25 pages, 10 figures, submitted to ApJ on 22-11-2022, revised until 8-3-2023, accepted 9-3-2023

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, Comment: Accepted for publication in The Planetary Sciences Journal

Recent experiments have shown that the nematode {\it T. aceti} can assemble into collectively undulating groups at the edge of fluid drops. This coordinated state consists of metachronal waves and drives fluid circulation inside the drop. We find that the circulation velocity is about 2 mm/s and nearly half the speed of the metachronal wave. We develop a quasi two-dimensional hydrodynamics model using the Stokes flow approximation. The periodic motion of the nematodes constitute our moving boundary condition that drives the flow. Our model suggests that large amplitude excursions of the nematodes tails produce the fluid circulation. We discuss the constraints on containers that would enhance fluid motion, which could be used in the future design of on demand flow generating systems.

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.

"What Will They Learn?" has been documenting for the last 13 years the crisis of weak core curriculum requirements. In addition to helping students choose the right college for the right reasons, American Council of Trustees and Alumni (ACTA) works with college and university leaders on campuses around the country to strengthen curricula and raise academic standards. The printed report and its companion website WhatWillTheyLearn.com, represent the only assessment of U.S. colleges and universities that focuses on what schools are actually teaching. The report assigns grades to over 1,100 colleges and universities on an "A" through "F" scale so that families know, in advance, just how serious an institution is about providing a rigorous core curriculum. Trustees, faculty members, administrators, policymakers, donors, students, parents, and high school counselors all play crucial roles in shaping American higher education. "What Will They Learn?"® is designed to help each constituency better understand the importance of a robust core curriculum, how to identify schools dedicated to preparing students for career and citizenship, and what each group can do to promote academic excellence at the nation's colleges and universities. [For the 2020 report, see ED616591.]

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.

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.

We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago., Comment: We welcome comments and missing references

This special report looks at how arts-based charter schools in the District of Columbia, Louisiana, Maryland and North Carolina provide access to arts education and integrate the arts schoolwide across curricula, community engagement and teacher professional development.

On a given day, more than 43,000 youths who largely identify as students of color and are disproportionately male are held in residential placement facilities as a result of involvement with the juvenile justice system. More than two-thirds of youths living in placement facilities who participated in a national survey shared aspirations to continue their education, with nearly half wanting to pursue postsecondary education. Despite these aspirations, youth involved in the juvenile justice system encounter educational barriers before they enter the system: Nearly one third are diagnosed with a learning disability, nearly half demonstrate academic proficiency below their grade level and close to a quarter are not enrolled in school. A disruption in their academic trajectory because of involvement with the juvenile justice system may prevent a student from continuing their education upon reentry. Sixty-six percent of youths do not return to school after leaving placement. High-quality arts education participation has lasting, positive impacts for youth by supporting development from early childhood into adulthood. The arts provide opportunities to build self-efficacy and achieve personal goals by helping youth develop ownership of their learning, determine individual criteria for success and track personal progress. Research shows that these effects most strongly benefit youth who have limited access to opportunities for academic, personal and social success and are thus at higher risk for low academic achievement, dropping out of high school or entering the juvenile justice system. Arts participation can also reduce the likelihood that they will engage in delinquent, risky or violent behavior. Additionally, the arts support civic outcomes, including political engagement and volunteering, among youth. This Special Report explores research and programs that employ the arts along the juvenile justice continuum: in prevention, intervention, transition and healing.

In this 12th edition of its flagship annual survey of college general education curricula, "What Will They Learn?," American Council of Trustees and Alumni (ACTA) shines a light on the wide disparity in academic standards at institutions of higher learning across the country. "What Will They Learn?" reviews the core curricula of over 1,100 four-year public and private institutions to determine the rigor of each school's general education program. Ratings are based on how many of the following seven core subjects are included in the curriculum: Composition, Literature, Foreign Language (intermediate-level), U.S. Government or History, Economics, Mathematics, and Natural Science. Grades are awarded on an "A" through "F" scale, depending on how many of the seven core subjects institutions require. This year, 365 schools earned an "A" or "B" grade, while over 400 received a "D" or "F." Too many colleges and universities continue to offer unstructured and chaotic curricula, causing students to sort through a plethora of course options in the name of "choice" and "self-discovery." "What Will They Learn?" 2020-2021 and its interactive companion website WhatWillTheyLearn.com can help orient students, parents, and trustees as they consider what educational quality looks like in 21st century America. [For the 2019 report, see ED602432.]

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., Comment: 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

When a sufficiently large far-field training data is presented, jointly optimizing a multichannel frontend and an end-to-end (E2E) Automatic Speech Recognition (ASR) backend shows promising results. Recent literature has shown traditional beamformer designs, such as MVDR (Minimum Variance Distortionless Response) or fixed beamformers can be successfully integrated as the frontend into an E2E ASR system with learnable parameters. In this work, we propose the self-attention channel combinator (SACC) ASR frontend, which leverages the self-attention mechanism to combine multichannel audio signals in the magnitude spectral domain. Experiments conducted on a multichannel playback test data shows that the SACC achieved a 9.3% WERR compared to a state-of-the-art fixed beamformer-based frontend, both jointly optimized with a ContextNet-based ASR backend. We also demonstrate the connection between the SACC and the traditional beamformers, and analyze the intermediate outputs of the SACC., Comment: In Proceedings of Interspeech 2021

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., Comment: to be published in Icarus

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

At high concentration, free swimming nematodes known as vinegar eels ({\it Turbatrix aceti}), collectively exhibit metachronal waves near a boundary. We find that the frequency of the collective traveling wave is lower than that of the freely swimming organisms. We explore models based on a chain of oscillators with nearest neighbor interactions that inhibit oscillator phase velocity. The phase of each oscillator represents the phase of the motion of the eel's head back and forth about its mean position. A strongly interacting directed chain model mimicking steric repulsion between organisms robustly gives traveling wave states and can approximately match the observed wavelength and oscillation frequency of the observed traveling wave. We predict body shapes assuming that waves propagate down the eel body at a constant speed. The phase oscillator model that impedes eel head overlaps also reduces close interactions throughout the eel bodies., Comment: Movies https://youtu.be/jETmIRah-mo and https://youtu.be/_HzUbkXZRKw

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} = -\alpha_{1}z$; $\log_{10} (\alpha_{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., Comment: accepted to ApJ Letters, changes following referee comments

High precision spectrographs can enable not only the discovery of exoplanets, but can also provide a fundamental measurement in Galactic dynamics. Over about ten year baselines, the expected change in the line-of-sight velocity due to the Galaxy's gravitational field for stars at $\sim$ kpc scale distances above the Galactic mid-plane is $\sim$ few - 10 cm/s, and may be detectable by the current generation of high precision spectrographs. Here, we provide theoretical expectations for this measurement based on both static models of the Milky Way and isolated Milky Way simulations, as well from controlled dynamical simulations of the Milky Way interacting with dwarf galaxies. We simulate a population synthesis model to analyze the contribution of planets and binaries to the Galactic acceleration signal. We find that while low-mass, long-period planetary companions are a contaminant to the Galactic acceleration signal, their contribution is very small. Our analysis of $\sim$ ten years of data from the LCES HIRES/Keck precision radial velocity (RV) survey shows that slopes of the RV curves of standard RV stars agree with expectations of the local Galactic acceleration near the Sun within the errors, and that the error in the slope scales inversely as the square root of the number of observations. Thus, we demonstrate that a survey of stars with low intrinsic stellar jitter at kpc distances above the Galactic mid-plane for realistic sample sizes can enable a direct determination of the dark matter density., Comment: accepted to ApJ Letters

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.

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., Comment: The animated movie is available here: http://astro.pas.rochester.edu/~aquillen/figures/movie3D.gif

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., Comment: 25 p., 19 figures, Accepted to MNRAS

Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites., Comment: Link to the high speed videos for the 3 experiments tracked and discussed: https://www.youtube.com/playlist?list=PLYnpHw29go3rKYOvM4WEQkHWYYIGgRC-E Preliminary results presented at American Geophysical Union (AGU) 2019 Fall Meeting (P43D-3497). Submitted to the journal Icarus on February 3, 2020

We numerically explore the behavior of repelling and aligning self-propelled polar particles (boids) in 2D enclosed by a damped flexible and elastic loop-shaped boundary. We observe disordered, polar ordered (or jammed) and circulating states. The latter produce a rich variety of boundary shapes including; circles, ovals, irregulars, ruffles, or sprockets, depending upon the bending moment of the boundary and the boundary to particle mass ratio. With the exception of the circulating states with non-round boundaries, states resemble those exhibited by attracting self-propelled particles, but here the confining boundary acts in place of a cohesive force. We attribute the formation of ruffles to instability mediated by pressure on the boundary when the speed of waves on the boundary approximately matches the self-propelled particle's swim speed.

We show that Smoothed Particle Hydrodynamics (SPH) simulations of dwarf galaxies interacting with a Milky Way-like disk produce moving groups in the simulated stellar disk. We analyze three different simulations: one that includes dwarf galaxies that mimic the Large Magellanic Cloud, Small Magellanic Cloud and the Sagittarius dwarf spheroidal; another with a dwarf galaxy that orbits nearly in the plane of the Milky Way disk; and a null case that does not include a dwarf galaxy interaction. We present a new algorithm to find large moving groups in the $V_R, V_\phi$ plane in an automated fashion that allows us to compare velocity sub-structure in different simulations, at different locations, and at different times. We find that there are significantly more moving groups formed in the interacting simulations than in the isolated simulation. A number of dwarf galaxies are known to orbit the Milky Way, with at least one known to have had a close pericenter approach. Our analysis of simulations here indicates that dwarf galaxies like those orbiting our Galaxy produce large moving groups in the disk. Our analysis also suggests that some of the moving groups in the Milky Way may have formed due to dynamical interactions with perturbing dwarf satellites. The groups identified in the simulations by our algorithm have similar properties to those found in the Milky Way, including similar fractions of the total stellar population included in the groups, as well as similar average velocities and velocity dispersions., Comment: 15 pages, 12 figures, Accepted for publication in MNRAS. Major revisions in the accepted version of the paper