This paper investigates the dynamical evolution of embedded stellar clusters from the protocluster stage, through the embedded star-forming phase, and out to ages of 10 Myr -- after the gas has been removed from the cluster. The relevant dynamical properties of young stellar clusters are explored over a wide range of possible star formation environments using N-body simulations. Many realizations of equivalent initial conditions are used to produce robust statistical descriptions of cluster evolution including the cluster bound fraction, radial probability distributions, as well as the distributions of close encounter distances and velocities. These cluster properties are presented as a function of parameters describing the initial configuration of the cluster, including the initial cluster membership N, initial stellar velocities, cluster radii, star formation efficiency, embedding gas dispersal time, and the degree of primordial mass segregation. The results of this parameter space survey, which includes about 25,000 simulations, provide a statistical description of cluster evolution as a function of the initial conditions. We also present a compilation of the FUV radiation fields provided by these same cluster environments. The output distributions from this study can be combined with other calculations, such as disk photoevaporation models and planetary scattering cross sections, to ascertain the effects of the cluster environment on the processes involved in planet formation.

Motivated by kinematic observations of young embedded clusters, this paper explores possible kinematic signatures produced by asphericity and departures from initial virial equilibrium in these systems. Specifically, the kinematic quantity that is measured and calculated in this study is the distribution of the line-of-sight velocities as a function of position along the cluster. Although clusters are found within a wide range of sizes, we focus on the regime with stellar membership $N \sim 10^3$. The gravitational potential of these young clusters is dominated by the gas, and the geometry of the gas distribution is generalized to include axisymmetric (and triaxial) forms. With this loss of symmetry, the kinematic results thus depend on viewing angle. This work also considers the stars to begin their trajectories with subvirial speeds, as indicated by observations of core motions in such clusters. Our results determine the conditions necessary for the kinematic signature to display interesting structure, i.e., a non-spherical potential, a viewing angle that is not along one of the principal axes, and subvirial starting conditions. We characterize the effects on this signature due to projection angle, initial stellar velocities, cluster elongation, and star formation efficiency. Finally, we compare our theoretical results to recent kinematic observations of the Orion Nebula Cluster; we find that the observations can be explained provided that the cluster is non-spherical, starts with subvirial initial velocities, and is not viewed along a principal axis.

We study the collapse of protoclusters within a giant molecular cloud (GMC) to determine the conditions under which collapse is significantly disrupted. Motivated by observations of star forming regions which exhibit flattened cloud structures, this study considers collapsing protoclusters with disk geometries. The collapse of a 10^3 Msun protocluster initially a distance of 2-10 pc from a 10^3 - 10^6 Msun point mass is numerically calculated. Simulations with zero initial relative velocity between the two are completed as well as simulations with relative velocities consistent with those observed in GMCs. The results allow us to define the conditions under which it is safe to assume protocluster collapse proceeds as if in isolation. For instance, we find the collapse of a 10^3 Msun protocluster will be significantly disrupted if it is within 2-4 pc of a 10^4 Msun point mass. Thus, the collapse of a 10^3 Msun protocluster can be considered to proceed as if in isolation if it is more than ~ 4 pc away from a 10^4 Msun compact object. In addition, in no portion of the sampled parameter space does the gravitational interaction between the protocluster disk and the massive particle significantly disperse the disk into the background GMC. We discuss the distribution of clusters of young stellar objects within the Perseus and Mon R2 star forming regions, which are consistent with the results of our simulations and the limitations of our results in gas dominated regions such as the Orion cloud.

This paper explores the stability of an Earth-like planet orbiting a solar-mass star in the presence of a stellar companion using ~ 400,000 numerical integrations. Given the chaotic nature of the systems being considered, we perform a statistical analysis of the ensuing dynamics for ~500 orbital configurations defined by the following set of orbital parameters: the companion mass; the companion eccentricity; the companion periastron; and the planet's inclination angle relative to the stellar binary plane. Specifically, we generate a large sample of survival times for each orbital configuration through the numerical integration of N >> 1 equivalent experiments (e.g., with the same orbital parameters but randomly selected initial orbital phases). We then construct distributions of survival time using the variable mu_s = log tau_s (where tau_s is in years) for each orbital configuration. The primary objective of this work is twofold. First, we use the mean of the distributions to gain a better understanding of what orbital configurations, while unstable, have sufficiently long survival times to make them interesting to the study of planet habitability. Second, we calculate the width, skew, and kurtosis of each mu_s distribution and look for general features that may aid further understanding and numerical exploration of these chaotic systems.

This paper studies the dynamical evolution of young stellar clusters with $N$ = 100 - 1000 members. We use N-body simulations to explore how evolution depends on system size $N$ and the initial conditions. Motivated by recent observations of extremely young systems, this study compares subvirial and virial starting states. Multiple realizations of equivalent cases (100 simulations per case) are used to build up a robust statistical description of these systems, e.g., distributions of closest approaches, mass profiles, and distributions of radial locations. These results provide a framework from which to assess the effects of these clusters on star and planet formation. The distributions of radial positions are used in conjunction with distributions of FUV luminosities (also calculated here) to determine the radiation exposure of circumstellar disks. The distributions of closest approaches are used in conjunction with scattering cross sections (calculated here from $10^5$ scattering experiments) to determine the probability of solar system disruption. We also use the nearby cluster NGC 1333 as a test case. Our main conclusion is that clusters in this size range have only a modest effect on forming planetary systems: Interaction rates are low so that the typical solar system experiences a single encounter within 1000 AU. Radiation exposure is low, with median FUV flux $G_0$ = 900, so that photoevaporation of disks is only important beyond 30 AU. Given the low interaction rates and modest radiation levels, we suggest that solar system disruption is a rare event in these clusters.