Planet-planet perturbations can cause planets' orbital elements to change on secular timescales. Previous work has evaluated the nodal precession rate for planets in the limit of low $\alpha$ (semi-major axis ratio, 0$<$$\alpha$$\leq$1). Our simulations show that systems at high $\alpha$ (or low period ratio), similar to multiplanet systems found in the Kepler survey, have a nodal precession rate that is more strongly dependent on eccentricity and inclination. We present a complete expansion of the nodal precession rate to fourth order in the disturbing function and show that this analytical solution much better describes the simulated N-body behavior of high-$\alpha$ planet pairs; at $\alpha\approx$ 0.5, the fourth-order solution on average reduces the median analytical error by a factor of 7.5 from linear theory and 6.2 from a second-order expansion. We set limits on eccentricity and inclination where the theory is precisely validated by N-body integrations, which can be useful in future secular treatments of planetary systems.

Wide-orbit exoplanets are starting to be detected, and planetary formation models are under development to understand their properties. We propose a population of "Oort" planets around other stars, forming by a mechanism analogous to how the Solar System's Oort cloud of comets was populated. Gravitational scattering among planets is inferred from the eccentricity distribution of gas-giant exoplanets measured by the Doppler technique. This scattering is thought to commence while the protoplanetary disk is dissipating, $10^6-10^7$ yr after formation of the star, or perhaps soon thereafter, when the majority of stars are expected to be part of a natal cluster. Previous calculations of planet-planet scattering around isolated stars have one or more planets spending $10^4-10^7$ yr at distances >100 AU before ultimately being ejected. During that time, a close flyby of another star in the cluster may dynamically lift the periastron of the planet, ending further scattering with the inner planets. We present numerical simulations demonstrating this mechanism as well as an analysis of the efficiency. We estimate an occurrence of planets between 100 and 5000 AU by this mechanism to be <1% for gas giants and up to a few percent for Neptunes and super-Earths.

Previous studies show that the physical structures and kinematics of a region depend significantly on the ionisation fraction. In this paper, we extend our previous studies of the effect of ionisation fractions on star formation to clouds that include both non-ideal magnetohydrodynamics and turbulence. We aim to quantify the importance of a treatment of the ionisation fraction in turbulent magnetised media and investigate the effect of turbulence on shaping the clouds and filaments before star formation sets in. In particular, we investigate how the structure, mass and width of filamentary structures depend on the amount of turbulence in ionised media and the initial mass-to-flux ratio. We compare the resulting density and mass-to-flux ratio structures both qualitatively and quantitatively via filament and core masses and filament fitting techniques (Gaussian and Plummer profiles.) We find that even with almost no turbulence, filamentary structure still exists. Comparison of simulations show that for turbulent Mach numbers above 2, there is little structural difference between the Step-Like (SL) and Cosmic Ray only (CR-only) ionisation models, while below this threshold the ionisation structure significantly affects the formation of filaments. Analysis of the mass within cores and filaments show decrease in mass as the degree of turbulence increases. Finally, observed filaments within the Taurus L1495/B213 complex are best reproduced by models with supercritical mass-to-flux ratios and/or at least mildly supersonic turbulence, however, our models show that sterile fibres observed within Taurus may occur in highly ionised, subcritical environments. Based on this, we suggest that regions with fertile fibres likely indicate a trans- or supercritical mass-to-flux ratio within the region while sterile fibres are likely subcritical and transient.

Bailey & Basu (2014) show analysis of density and mass-to-flux ratio maps for simulations with either an ionisation profile which takes into account photoionisation (step-like profile) or a cosmic ray only ionisation profile. We extend this study to analyse the effect of these ionisation profiles on velocity structures, kinematics, and synthetic spectra. Clump regions are found to occur at the convergence of two flows with a low velocity region and velocity direction transition occurring at the junction. Models with evident substructure show that core formation occurs on the periphery of these velocity valleys. Analysis of synthetic spectra reveals the presence of large non-thermal components within low-density gas, especially for models with the step-like ionisation profile. All cores show small, sub-thermal relative motions compared to background gas. Large deviations within this analysis are due to the line of sight intersecting low- and high-density regions across the velocity switch transition. Positive deviations correspond to a foreground core moving away from the observer while negative deviations correspond to a background core moving toward the observer. Comparison of velocities resulting from different ionisation profiles suggest that high ionisation fractions yield supersonic velocities, up to two times the sound speed, while regions with low ionisation fractions tend to be subsonic or mildly supersonic. This suggests that the transition to coherence within cores could be a transition between high and low ionisation fractions within the gas.

We model molecular cloud fragmentation with thin disk non-ideal magnetohydrodynamic simulations that include ambipolar diffusion and partial ionization that transitions from primarily ultraviolet dominated to cosmic ray dominated regimes. These simulations are used to determine the conditions required for star clusters to form through a two-stage fragmentation scenario. Recent linear analyses have shown that the fragmentation length and time scales can undergo a dramatic drop across the column density boundary that separates the ultraviolet and cosmic ray dominated ionization regimes. As found in earlier studies, the absence of an ionization drop and regular perturbations leads to a single-stage fragmentation on parsec scales in transcritical clouds, so that the nonlinear evolution yields the same fragment sizes as predicted by linear theory. However, we find that a combination of initial transcritical mass-to-flux ratio, evolution through a column density regime in which the ionization drop takes place, and regular small perturbations to the mass-to-flux ratio are sufficient to cause a second stage of fragmentation during the nonlinear evolution. Cores of size ~0.1 pc are formed within an initial fragment of ~ pc size. Regular perturbations to the mass-to-flux ratio also accelerate the onset of runaway collapse.

We review recent results of non-ideal magnetohydrodynamic models for the fragmentation of molecular clouds and the collapse of cloud cores to form protostar-disk systems. Thin-disk models can elucidate many aspects of the physical problem and allow the calculation of large dynamic range of time and length scales.

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds yields information about length scales involved in star formation. Combining these length scales with various distributions of other environmental variables, (i.e., column density and mass-to-flux ratio) and applying Monte Carlo methods allow us to produce synthetic core mass functions (CMFs) for different environmental conditions. Our analysis shows that the shape of the CMF is directly dependent on the physical conditions of the cloud. Specifically, magnetic fields act to broaden the mass function and develop a high-mass tail while ambipolar diffusion will truncate this high-mass tail. In addition, we analyze the effect of small number statistics on the shape and high-mass slope of the synthetic CMFs. We find that observed core mass functions are severely statistically limited, which has a profound effect on the derived slope for the high-mass tail.

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds shows that molecular clouds exhibit a preferred length scale for collapse that depends on the mass-to-flux ratio and neutral-ion collision time within the cloud. We extend this linear analysis to the context of clustered star formation. By combining the results of the linear analysis with a realistic ionization profile for the cloud, we find that a molecular cloud may evolve through two fragmentation events in the evolution toward the formation of stars. Our model suggests that the initial fragmentation into clumps occurs for a transcritical cloud on parsec scales while the second fragmentation can occur for transcritical and supercritical cores on subparsec scales. Comparison of our results with several star forming regions (Perseus, Taurus, Pipe Nebula) shows support for a two-stage fragmentation model.