The virial ratio between kinetic and gravitational terms provides key insight into the balance of forces that confine a molecular cloud, but the clumpy and filamentary structures of resolved clouds make it difficult to evaluate this ratio in a consistent way. For clouds with resolved maps of column density as well as a line tracer, we demonstrate that the gravitational energy can be estimated directly from observations in a manner similar to the kinetic energy. This offers improved diagnostic power and consistency. Disentangling a cloud from foreground and background materials is a persistent challenge, for which we introduce a strategy based on Abel's transform. We provide proofs of principle using simulated clouds.

Dust barriers effectively capture the photon momentum of a central light source, but low-density channels, along with re-emission at longer wavelengths, enhance its escape. We use Monte Carlo simulations to study the effects of inhomogeneity on radiation forces imparted to a dust envelope around a central star. We survey the strength and scale of an inhomogeneous perturbation field, as well as the optical depth of its spherical reference state. We run at moderate numerical resolution, relying on our previous resolution study for calibration of the associated error. We find that inhomogeneities matter most when their scale exceeds the characteristic mean free path. As expected, they tend to reduce the net radiation force and extend its range; however, there is significant variance among realizations. Within our models, force integrals correlate with the emergent spectral energy distribution, given a specified set of dust properties. A critical issue is the choice of integral measures of the radiation force: for strong deviations from spherical symmetry the relevant measures assess radial forces relative to the cloud centre of mass. Of these, we find the virial term due to radiation to be the least stochastic of several integral measures in the presence of inhomogeneities.

We address in detail the radiation forces on spherical dust envelopes around luminous stars, and numerical solutions for these forces, as a first step toward more general dust geometries. Two physical quantities, a normalized force and a force-averaged radius, suffice to capture the overall effects of radiation forces; these combine to yield the radiation term in the virial theorem. In addition to the optically thin and thick regimes, the wavelength dependence of dust opacity allows for an intermediate case in which starlight is easily trapped but infrared radiation readily escapes. Scattering adds a non-negligible force in this intermediate regime. We address all three regimes analytically and provide approximate formulae for the force parameters, for arbitrary optical depth and inner dust temperature. Turning to numerical codes, we examine the convergence properties of the Monte Carlo code Hyperion run in Cartesian geometry. We calibrate both Hyperion and our analytical estimate using the DUSTY code, run in spherical geometry. We find that Monte Carlo codes tend to underestimate the radiation force when the mean free path of starlight is not well resolved, as this causes the inner dust temperature, and therefore the inner Rosseland opacity, to be too low. We briefly discuss implications for more complicated radiation transfer problems.

During star cluster formation, ongoing mass accretion is resisted by stellar feedback in the form of protostellar outflows from the low-mass stars and photo-ionization and radiation pressure feedback from the massive stars. We model the evolution of cluster-forming regions during a phase in which both accretion and feedback are present, and use these models to investigate how star cluster formation might terminate. Protostellar outflows are the strongest form of feedback in low-mass regions, but these cannot stop cluster formation if matter continues to flow in. In more massive clusters, radiation pressure and photo-ionization rapidly clear the cluster-forming gas when its column density is too small. We assess the rates of dynamical mass ejection and of evaporation, while accounting for the important effect of dust opacity on photo-ionization. Our models are consistent with the census of protostellar outflows in NGC 1333 and Serpens South, and with the dust temperatures observed in regions of massive star formation. Comparing observations of massive cluster-forming regions against our model parameter space, and against our expectations for accretion-driven evolution, we infer that massive-star feedback is a likely cause of gas disruption in regions with velocity dispersions less than a few kilometers per second, but that more massive and more turbulent regions are too strongly bound for stellar feedback to be disruptive.

The formation of brown dwarfs (BDs) poses a key challenge to star formation theory. The observed dearth of nearby ($\leq 5$ AU) brown dwarf companions to solar-mass stars, known as the brown dwarf desert, as well as the tendency for low-mass binary systems to be more tightly-bound than stellar binaries, have been cited as evidence for distinct formation mechanisms for brown dwarfs and stars. In this paper, we explore the implications of the minimal hypothesis that brown dwarfs in binary systems originate via the same fundamental fragmentation mechanism as stars, within isolated, turbulent giant molecular cloud cores. We demonstrate analytically that the scaling of specific angular momentum with turbulent core mass naturally gives rise to the brown dwarf desert, as well as wide brown-dwarf binary systems. Further, we show that the turbulent core fragmentation model also naturally predicts that very low-mass (VLM) binary and BD/BD systems are more tightly-bound than stellar systems. In addition, in order to capture the stochastic variation intrinsic to turbulence, we generate $10^4$ model turbulent cores with synthetic turbulent velocity fields to show that the turbulent fragmentation model accommodates a small fraction of binary brown dwarfs with wide separations, similar to observations. Indeed, the picture which emerges from the turbulent fragmentation model is that a single fragmentation mechanism may largely shape both stellar and brown dwarf binary distributions during formation.