We present an analysis of the relationship between the orientation of magnetic fields and filaments that form in 3D magnetohydrodynamic simulations of cluster-forming, turbulent molecular cloud clumps. We examine simulated cloud clumps with size scales of L ~ 2-4 pc and densities of n ~ 400-1000 cm^-3 with Alfven Mach numbers near unity. We simulated two cloud clumps of different masses, one in virial equilibrium, the other strongly gravitationally bound, but with the same initial turbulent velocity field and similar mass-to-flux ratio. We apply various techniques to analyze the filamentary and magnetic structure of the resulting cloud, including the DisPerSE filament-finding algorithm in 3D. The largest structure that forms is a 1-2 parsec-long filament, with smaller connecting sub-filaments. We find that our simulated clouds, wherein magnetic forces and turbulence are comparable, coherent orientation of the magnetic field depends on the virial parameter. Subvirial clumps undergo strong gravitational collapse and magnetic field lines are dragged with the accretion flow. We see evidence of filament-aligned flow and accretion flow onto the filament in the subvirial cloud. Magnetic fields oriented more parallel in the subvirial cloud and more perpendicular in the denser, marginally bound cloud. Radiative feedback from a 16 Msun star forming in a cluster in one of our simulations ultimately results in the destruction of the main filament, the formation of an HII region, and the sweeping up of magnetic fields within an expanding shell at the edges of the HII region.

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.

We use numerical simulations of turbulent cluster-forming regions to study the nature of dense filamentary structures in star formation. Using four hydrodynamic and magnetohydrodynamic simulations chosen to match observations, we identify filaments in the resulting column density maps and analyze their properties. We calculate the radial column density profiles of the filaments every 0.05 Myr and fit the profiles with the modified isothermal and pressure confined isothermal cylinder models, finding reasonable fits for either model. The filaments formed in the simulations have similar radial column density profiles to those observed. Magnetic fields provide additional pressure support to the filaments, making `puffier' filaments less prone to fragmentation than in the pure hydrodynamic case, which continue to condense at a slower rate. In the higher density simulations, the filaments grow faster through the increased importance of gravity. Not all of the filaments identified in the simulations will evolve to form stars: some expand and disperse. Given these different filament evolutionary paths, the trends in bulk filament width as a function of time, magnetic field strength, or density, are weak, and all cases are reasonably consistent with the finding of a constant filament width in different star-forming regions. In the simulations, the mean FWHM lies between 0.06 and 0.26 pc for all times and initial conditions, with most lying between 0.1 to 0.15 pc; the range in FWHMs are, however, larger than seen in typical Herschel analyses. Finally, the filaments display a wealth of substructure similar to the recent discovery of filament bundles in Taurus.

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

Recent observations and simulations have suggested that HII regions around massive stars may vary in their size and emitted flux on timescales short enough to be observed. This variability can have a number of causes, ranging from environmental causes to variability of the ionizing source itself. We explore the latter possibility by considering the pre-main-sequence evolution of massive protostars and conducting numerical simulations with ionizing radiation feedback using the FLASH AMR hydrodynamics code. We investigate three different models: a simple ZAMS model, a self-consistent one-zone model by Offner et al. (2009), and a model fit to the tracks computed by Hosokawa & Omukai (2009). The protostellar models show that hypercompact HII regions around massive, isolated protostars collapse or shrink from diameters of 80 or 300 AU, depending on the model choice, down to near absence during the swelling of stellar radius that accompanies the protostar's transition from a convective to a radiative internal structure. This occurs on timescales as short as ~3000 years.

Radiative feedback is among the most important consequences of clustered star formation inside molecular clouds. At the onset of star formation, radiation from massive stars heats the surrounding gas, which suppresses the formation of many low-mass stars. When simulating pre-main-sequence stars, their stellar properties must be defined by a prestellar model. Different approaches to prestellar modeling may yield quantitatively different results. In this paper, we compare two existing prestellar models under identical initial conditions to gauge whether the choice of model has any significant effects on the final population of stars. The first model treats stellar radii and luminosities with a ZAMS model, while separately estimating the accretion luminosity by interpolating to published prestellar tracks. The second, more accurate prestellar model self-consistently evolves the radius and luminosity of each star under highly variable accretion conditions. Each is coupled to a raytracing-based radiative feedback code that also treats ionization. The impact of the self-consistent model is less ionizing radiation and less heating during the early stages of star formation. This may affect final mass distributions. We noted a peak stellar mass reduced by 8% from 47.3 Msun to 43.5 Msun in the evolutionary model, relative to the track-fit model. Also, the difference in mass between the two largest stars in each case is reduced from 14 Msun to 7.5 Msun. The HII regions produced by these massive stars were also seen to flicker on timescales down to the limit imposed by our timestep (< 560 years), rapidly changing in size and shape, confirming previous cluster simulations using ZAMS-based estimates for prestellar ionizing flux.