Sub-structure formation in starless cores

Motivated by recent observational searches of sub-structure in starless molecular cloud cores, we investigate the evolution of density perturbations on scales smaller than the Jeans length embedded in contracting isothermal clouds, adopting the same formalism developed for the expanding Universe and the solar wind. We find that initially small amplitude, Jeans-stable perturbations (propagating as sound waves in the absence of a magnetic field), are amplified adiabatically during the contraction, approximately conserving the wave action density, until they either become nonlinear and steepen into shocks at a time $t_{\rm nl}$, or become gravitationally unstable when the Jeans length decreases below the scale of the perturbations at a time $t_{\rm gr}$. We evaluate analytically the time $t_{\rm nl}$ at which the perturbations enter the non-linear stage using a Burgers' equation approach, and we verify numerically that this time marks the beginning of the phase of rapid dissipation of the kinetic energy of the perturbations. We then show that for typical values of the rms Mach number in molecular cloud cores, $t_{\rm nl}$ is smaller than $t_{\rm gr}$, and therefore density perturbations likely dissipate before becoming gravitational unstable. Solenoidal modes grow at a faster rate than compressible modes, and may eventually promote fragmentation through the formation of vortical structures.Comment: 8 pages, 4 figure

A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows

In this paper we present a genuinely two-dimensional HLLC Riemann solver. On logically rectangular meshes, it accepts four input states that come together at an edge and outputs the multi-dimensionally upwinded fluxes in both directions. This work builds on, and improves, our prior work on two-dimensional HLL Riemann solvers. The HLL Riemann solver presented here achieves its stabilization by introducing a constant state in the region of strong interaction, where four one-dimensional Riemann problems interact vigorously with one another. A robust version of the HLL Riemann solver is presented here along with a strategy for introducing sub-structure in the strongly-interacting state. Introducing sub-structure turns the two-dimensional HLL Riemann solver into a two-dimensional HLLC Riemann solver. The sub-structure that we introduce represents a contact discontinuity which can be oriented in any direction relative to the mesh. The Riemann solver presented here is general and can work with any system of conservation laws. We also present a second order accurate Godunov scheme that works in three dimensions and is entirely based on the present multidimensional HLLC Riemann solver technology. The methods presented are cost-competitive with traditional higher order Godunov schemes

Sub-structural Niching in Estimation of Distribution Algorithms

We propose a sub-structural niching method that fully exploits the problem decomposition capability of linkage-learning methods such as the estimation of distribution algorithms and concentrate on maintaining diversity at the sub-structural level. The proposed method consists of three key components: (1) Problem decomposition and sub-structure identification, (2) sub-structure fitness estimation, and (3) sub-structural niche preservation. The sub-structural niching method is compared to restricted tournament selection (RTS)--a niching method used in hierarchical Bayesian optimization algorithm--with special emphasis on sustained preservation of multiple global solutions of a class of boundedly-difficult, additively-separable multimodal problems. The results show that sub-structural niching successfully maintains multiple global optima over large number of generations and does so with significantly less population than RTS. Additionally, the market share of each of the niche is much closer to the expected level in sub-structural niching when compared to RTS

J/psi dissociation by light mesons in an extended Nambu Jona-Lasinio model

An alternative model for the dissociation of the J/psi is proposed. Chiral symmetry is properly implemented. Abnormal parity interactions and mesonic form factors naturally arise from the underlying quark sub-structure. Analytic confinement for the light quarks is generated by appropriately chosen the quark interaction kernels. Dissociation cross sections of the J/psi by either a pion or a rho meson are then evaluated and discussed.Comment: 24 pages, 13 figures, final versio

3D Simulations of Betelgeuse's Bow Shock

Betelgeuse, the bright, cool red supergiant in Orion, is moving supersonically relative to the local interstellar medium. The star emits a powerful stellar wind which collides with this medium, forming a cometary structure, a bow shock, pointing in the direction of motion. We present the first 3D hydrodynamic simulations of the formation and evolution of Betelgeuse's bow shock. The models include realistic low temperature cooling and cover a range of plausible interstellar medium densities and stellar velocities between 0.3 - 1.9 cm-3 and 28 - 73 km/s. We show that the flow dynamics and morphology of the bow shock differ substantially due to the preferential growth of Rayleigh-Taylor or Kelvin-Helmholtz instabilities in the models. The former dominate the models with slow stellar velocities resulting in a clumpy bow shock sub-structure, whereas the latter produce a smoother, more layered sub-structure in the fast models. If the mass in the bow shock shell is low, as seems to be implied by the AKARI luminosities (~0.003 Msun), then Betelgeuse's bow shock is very young and is unlikely to have reached a steady state. The circular nature of the bow shock shell is consistent with this conclusion. Thus, our results suggest that Betelgeuse entered the red supergiant phase only recently.Comment: Minor revisions, replaced Fig. 1, 15, and 16, added movies. For a pdf version with higher resolution, see A&A: Forthcomin

The spiral structure of the Galaxy revealed by CS sources and evidence for the 4:1 resonance

We present a map of the spiral structure of the Galaxy, as traced by molecular CS emission associated with IRAS sources which are believed to be compact HII regions. The CS line velocities are used to determine the kinematic distances of the sources, in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical HII regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The well defined spiral arms are seen to be confined inside the co-rotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (4 epicycle oscillations in one turn around the galactic center). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inward curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky