We use $N$-body simulations to investigate the excitation of bending waves in a Milky Way-like disc-bulge-halo system. The dark matter halo consists of a smooth component and a population of subhaloes while the disc is composed of thin and thick components. Also considered is a control simulation where all of the halo mass is smoothly distributed. We find that bending waves are more vigorously excited in the thin disc than the thick one and that they are strongest in the outer regions of the disc, especially at late times. By way of a Fourier decomposition, we find that the complicated pattern of bending across the disc can be described as a superposition of waves, which concentrate along two branches in the radius-rotational frequency plane. These branches correspond to vertical resonance curves as predicted by a WKB analysis. Bending waves in the simulation with substructure have a higher amplitude than those in the smooth-halo simulation, though the frequency-radius characteristics of the waves in the two simulations are very similar. A cross correlation analysis of vertical displacement and bulk vertical velocity suggests that the waves oscillate largely as simple plane waves. We suggest that the wave-like features in astrometric surveys such as the Second Data Release from \textit{Gaia} may be due to long-lived waves of a dynamically active disc rather than, or in addition to, perturbations from a recent satellite-disc encounter.

We study the spontaneous generation and evolution of bending waves in $N$-body simulations of two isolated Milky Way-like galaxy models. The models differ by their disc-to-halo mass ratios, and hence by their susceptibility to the formation of a bar and spiral structure. Seeded from shot noise in the particle distribution, bending waves rapidly form in both models and persist for many billions of years. Waves at intermediate radii manifest as corrugated structures in vertical position and velocity that are tightly wound, morphologically leading, and dominated by the $m=1$ azimuthal Fourier component. A spectral analysis of the waves suggests they are a superposition of modes from two continuous branches in the Galactocentric radius-rotational frequency plane. The lower-frequency branch is dominant and is responsible for the corrugated, leading, and warped structure. Over time, power in this branch migrates outward, lending credence to an inside-out formation scenario for the warp. Our power spectra qualitatively agree with results from linear perturbation theory and a WKB analysis, both of which include self-gravity. Thus, we conclude that the waves in our simulations are self-gravitating and not purely kinematic. These waves are reminiscent of the wave-like pattern recently found in Galactic star counts from the Sloan Digital Sky Survey and smoothly transition to a warp near the disc's edge. Velocity measurements from \textit{Gaia} data will be instrumental in testing the true wave nature of the corrugations. We also compile a list of "minimum requirements" needed to observe bending waves in external galaxies.

We present a disc-halo N-body model of the low surface brightness galaxy UGC 628, one of the few systems that harbours a "slow" bar with a ratio of corotation radius to bar length of $\mathcal{R} \equiv R_c/a_b \sim 2$. We select our initial conditions using SDSS DR10 photometry, a physically motivated radially variable mass-to-light ratio profile, and rotation curve data from the literature. A global bar instability grows in our submaximal disc model, and the disc morphology and dynamics agree broadly with the photometry and kinematics of UGC 628 at times between peak bar strength and the onset of buckling. Prior to bar formation, the disc and halo contribute roughly equally to the potential in the galaxy's inner region, giving the disc enough self gravity for bar modes to grow. After bar formation there is significant mass redistribution, creating a baryon dominated inner and dark matter dominated outer disc. This implies that, unlike most other low surface brightness galaxies, UGC 628 is not dark matter dominated everywhere. Our model nonetheless implies that UGC 628 falls on same the relationship between dark matter fraction and rotation velocity found for high surface brightness galaxies, and lends credence to the argument that the disc mass fraction measured at the location where its contribution to the potential peaks is not a reliable indicator of its dynamical importance at all radii.

We explore the hypothesis that a passing satellite or dark matter subhalo has excited coherent oscillations of the Milky Way's stellar disk in the direction perpendicular to the Galactic midplane. This work is motivated by recent observations of spatially dependent bulk vertical motions within ~ kpc of the Sun. A satellite can transfer a fraction of its orbital energy to the disk stars as it plunges through the Galactic midplane thereby heating and thickening the disk. Bulk motions arise during the early stages of such an event when the disk is still in an unrelaxed state. We present simple toy-model calculations and simulations of disk-satellite interactions, which show that the response of the disk depends on the relative velocity of the satellite. When the component of the satellite's velocity perpendicular to the disk is small compared with that of the stars, the perturbation is predominantly a bending mode. Conversely, breathing and higher order modes are excited when the vertical velocity of the satellite is larger than that of the stars. We argue that the compression and rarefaction motions seen in three different surveys are in fact breathing mode perturbations of the Galactic disk.