Helical Fields and Filamentary Molecular Clouds

We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first derive a new virial equation appropriate for magnetized filamentary clouds, which includes the effects of non-thermal motions and the turbulent pressure of the surrounding ISM. When compared with the data, we find that many filamentary clouds have a mass per unit length that is significantly reduced by the effects of external pressure, and that toroidal fields play a significant role in squeezing such clouds. We also develop exact numerical MHD models of filamentary molecular clouds with more general helical field configurations than have previously been considered. We also examine the effects of the equation of state by comparing ``isothermal'' filaments, with constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed using a logatropic equation of state. We perform a Monte Carlo exploration of our parameter space to determine which choices of parameters result in models that agree with the available observational constraints. We find that both equations of state result in equilibria that agree with the observational results. Moreover, we find that models with helical fields have more realistic density profiles than either unmagnetized models or those with purely poloidal fields; we find that most isothermal models have density distributions that fall off as r^{-1.8} to r^{-2}, while logatropes have density profiles that range from r^{-1} to r^{-1.8}. We find that purely poloidal fields produce filaments with steep density gradients that not allowed by the observations.Comment: 21 pages, 8 eps figures, submitted to MNRAS. Significant streamlining of tex

Helical Fields and Filamentary Molecular Clouds II - Axisymmetric Stability and Fragmentation

In Paper I (Fiege & Pudritz, 1999), we constructed models of filamentary molecular clouds that are truncated by a realistic external pressure and contain a rather general helical magnetic field. We address the stability of our models to gravitational fragmentation and axisymmetric MHD-driven instabilities. By calculating the dominant modes of axisymmetric instability, we determine the dominant length scales and growth rates for fragmentation. We find that the role of pressure truncation is to decrease the growth rate of gravitational instabilities by decreasing the self-gravitating mass per unit length. Purely poloidal and toroidal fields also help to stabilize filamentary clouds against fragmentation. The overall effect of helical fields is to stabilize gravity-driven modes, so that the growth rates are significantly reduced below what is expected for unmagnetized clouds. However, MHD ``sausage'' instabilities are triggered in models whose toroidal flux to mass ratio exceeds the poloidal flux to mass ratio by more than a factor of $\sim 2$. We find that observed filaments appear to lie in a physical regime where the growth rates of both gravitational fragmentation and axisymmetric MHD-driven modes are at a minimum.Comment: 16 pages with 18 eps figures. Submitted to MNRA

Decorated Feynman Categories

In [KW14], the new concept of Feynman categories was introduced to simplify the discussion of operad--like objects. In this present paper, we demonstrate the usefulness of this approach, by introducing the concept of decorated Feynman categories. The procedure takes a Feynman category $\mathfrak F$ and a functor $\mathcal O$ to a monoidal category to produce a new Feynman category ${\mathfrak F}_{dec {\mathcal O}}$. This in one swat explains the existence of non--sigma operads, non--sigma cyclic operads, and the non--sigma--modular operads of Markl as well as all the usual candidates simply from the category $\mathfrak G$, which is a full subcategory of the category of graphs of [BM08]. Moreover, we explain the appearance of terminal objects noted in [Mar15]. We can then easily extend this for instance to the dihedral case. Furthermore, we obtain graph complexes and all other known operadic type notions from decorating and restricting the basic Feynman category $\mathfrak G$ of aggregates of corollas. We additionally show that the construction is functorial. There are further geometric and number theoretic applications, which will follow in a separate preprint.Comment: Updated version from 6/8/16 to appear in J. of Noncommutative Geometr

Resistive Magnetohydrodynamic Equilibria in a Torus

It was recently demonstrated that static, resistive, magnetohydrodynamic equilibria, in the presence of spatially-uniform electrical conductivity, do not exist in a torus under a standard set of assumed symmetries and boundary conditions. The difficulty, which goes away in the ``periodic straight cylinder approximation,'' is associated with the necessarily non-vanishing character of the curl of the Lorentz force, j x B. Here, we ask if there exists a spatial profile of electrical conductivity that permits the existence of zero-flow, axisymmetric r esistive equilibria in a torus, and answer the question in the affirmative. However, the physical properties of the conductivity profile are unusual (the conductivity cannot be constant on a magnetic surface, for example) and whether such equilibria are to be considered physically possible remains an open question.Comment: 17 pages, 4 figure

Tracking shocked dust: state estimation for a complex plasma during a shock wave

We consider a two-dimensional complex (dusty) plasma crystal excited by an electrostatically-induced shock wave. Dust particle kinematics in such a system are usually determined using particle tracking velocimetry. In this work we present a particle tracking algorithm which determines the dust particle kinematics with significantly higher accuracy than particle tracking velocimetry. The algorithm uses multiple extended Kalman filters to estimate the particle states and an interacting multiple model to assign probabilities to the different filters. This enables the determination of relevant physical properties of the dust, such as kinetic energy and kinetic temperature, with high precision. We use a Hugoniot shock-jump relation to calculate a pressure-volume diagram from the shocked dust kinematics. Calculation of the full pressure-volume diagram was possible with our tracking algorithm, but not with particle tracking velocimetry.Comment: 10 pages, 8 figures, accepted for publication in Physics of Plasma