Electric potential distributions at the interface between plasmasheet clouds

At the interface between two plasma clouds with different densities, temperatures, and/or bulk velocities, there are large charge separation electric fields which can be modeled in the framework of a collisionless theory for tangential discontinuities. Two different classes of layers were identified: the first one corresponds to (stable) ion layers which are thicker than one ion Lamor radius; the second one corresponds to (unstable) electron layers which are only a few electron Larmor radii thick. It is suggested that these thin electron layers with large electric potential gradients (up to 400 mV/m) are the regions where large-amplitude electrostatic waves are spontaneously generated. These waves scatter the pitch angles of the ambient plasmasheet electron into the atmospheric loss cone. The unstable electron layers can therefore be considered as the seat of strong pitch angle scattering for the primary auroral electrons

The nature of turbulence in OMC1 at the star forming scale: observations and simulations

Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions. Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence. Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.Comment: 12 pages, 8 figures, accepted A&A. Revised in response to referee. For higher resolution, see http://www.astro.phys.au.dk/~maikeng/sim_paper

Observations of spatial and velocity structure in the Orion Molecular Cloud

Observations are reported of H2 IR emission in the S(1) v=1-0 line at 2.121 microns in the Orion Molecular Cloud, OMC1, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF is a combination of adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a velocity discrimination as high as 1 km/s. Thanks to the high spatial and velocity resolution of the GriF data, 193 bright H2 emission regions can be identified in OMC1. The general characteristics of these features are described in terms of radial velocities, brightness and spatial displacement of maxima of velocity and brightness, the latter to yield the orientation of flows in the plane of the sky. Strong spatial correlation between velocity and bright H2 emission is found and serves to identify many features as shocks. Important results are: (i) velocities of the excited gas illustrate the presence of a zone to the south of BN-IRc2 and Peak 1, and the west of Peak 2, where there is a powerful blue-shifted outflow with an average velocity of -18 km/s. This is shown to be the NIR counterpart of an outflow identified in the radio from source I, a very young O-star. (ii) There is a band of weak velocity features (<5 km/s) in Peak 1 which may share a common origin through an explosive event, in the BN-IRc2 region, with the fast-moving fingers (or bullets) to the NW of OMC1. (iii) A proportion of the flows are likely to represent sites of low mass star formation and several regions show multiple outflows, probably indicative of multiple star formation within OMC1. The high spatial and velocity resolution of the GriF data show these and other features in more detail than has previously been possible.Comment: 27 pages, 19 figures, submitted to A&A Version 2: Several additions, including a section on protostellar candidates in OMC1, have been made based on the referee's suggestions v3: corrected typograph

Near-Infrared Spectroscopy of Molecular Hydrogen Emission in Four Reflection Nebulae: NGC 1333, NGC 2023, NGC 2068, and NGC 7023

We present near-infrared spectroscopy of fluorescent molecular hydrogen (H_2) emission from NGC 1333, NGC 2023, NGC 2068, and NGC 7023 and derive the physical properties of the molecular material in these reflection nebulae. Our observations of NGC 2023 and NGC 7023 and the physical parameters we derive for these nebulae are in good agreement with previous studies. Both NGC 1333 and NGC 2068 have no previously-published analysis of near-infrared spectra. Our study reveals that the rotational-vibrational states of molecular hydrogen in NGC 1333 are populated quite differently from NGC 2023 and NGC 7023. We determine that the relatively weak UV field illuminating NGC 1333 is the primary cause of the difference. Further, we find that the density of the emitting material in NGC 1333 is of much lower density, with n ~ 10^2 - 10^4 cm^-3. NGC 2068 has molecular hydrogen line ratios more similar to those of NGC 7023 and NGC 2023. Our model fits to this nebula show that the bright, H_2-emitting material may have a density as high as n ~ 10^5 cm^-3, similar to what we find for NGC 2023 and NGC 7023. Our spectra of NGC 2023 and NGC 7023 show significant changes in both the near-infrared continuum and H_2 intensity along the slit and offsets between the peaks of the H_2 and continuum emission. We find that these brightness changes may correspond to real changes in the density and temperatures of the emitting region, although uncertainties in the total column of emitting material along a given line of sight complicates the interpretation. The spatial difference in the peak of the H_2 and near-infrared continuum peaks in NGC 2023 and NGC 7023 shows that the near-infrared continuum is due to a material which can survive closer to the star than H_2 can.Comment: Submitted for publication in ApJ. 34 pages including 12 embedded postscript figures. Also available at http://www.astronomy.ohio-state.edu/~martini/pub

Generalized Forward-Backward Splitting

This paper introduces the generalized forward-backward splitting algorithm for minimizing convex functions of the form $F + \sum_{i=1}^n G_i$, where $F$ has a Lipschitz-continuous gradient and the $G_i$'s are simple in the sense that their Moreau proximity operators are easy to compute. While the forward-backward algorithm cannot deal with more than $n = 1$ non-smooth function, our method generalizes it to the case of arbitrary $n$. Our method makes an explicit use of the regularity of $F$ in the forward step, and the proximity operators of the $G_i$'s are applied in parallel in the backward step. This allows the generalized forward backward to efficiently address an important class of convex problems. We prove its convergence in infinite dimension, and its robustness to errors on the computation of the proximity operators and of the gradient of $F$. Examples on inverse problems in imaging demonstrate the advantage of the proposed methods in comparison to other splitting algorithms.Comment: 24 pages, 4 figure

Short-range repulsion and isospin dependence in the KN system

The short-range properties of the KN interaction are studied within the meson-exchange model of the Juelich group. Specifically, dynamical explanations for the phenomenological short-range repulsion, required in this model for achieving agreement with the empirical KN data, are explored. Evidence is found that contributions from the exchange of a heavy scalar-isovector meson (a_0(980)) as well as from genuine quark-gluon exchange processes are needed. Taking both mechanisms into account a satisfactory description of the KN phase shifts can be obtained without resorting to phenomenological pieces.Comment: 26 pages, 5 figure

Thomas Decomposition of Algebraic and Differential Systems

In this paper we consider disjoint decomposition of algebraic and non-linear partial differential systems of equations and inequations into so-called simple subsystems. We exploit Thomas decomposition ideas and develop them into a new algorithm. For algebraic systems simplicity means triangularity, squarefreeness and non-vanishing initials. For differential systems the algorithm provides not only algebraic simplicity but also involutivity. The algorithm has been implemented in Maple

Observations and modeling of H_2 fluorescence with partial frequency redistribution in giant planet atmospheres

Partial frequency redistribution (PRD), describing the formation of the line profile, has negligible observational effects for optical depths smaller than ~10^3, at the resolving power of most current instruments. However, when the spectral resolution is sufficiently high, PRD modeling becomes essential in interpreting the line shapes and determining the total line fluxes. We demonstrate the effects of PRD on the H_2 line profiles observed at high spectral resolution by the Far-Ultraviolet Spectroscopic Explorer (FUSE) in the atmospheres of Jupiter and Saturn. In these spectra, the asymmetric shapes of the lines in the Lyman (v"- 6) progression pumped by the solar Ly-beta are explained by coherent scattering of the photons in the line wings. We introduce a simple computational approximation to mitigate the numerical difficulties of radiative transfer with PRD, and show that it reproduces the exact radiative transfer solution to better than 10%. The lines predicted by our radiative transfer model with PRD, including the H_2 density and temperature distribution as a function of height in the atmosphere, are in agreement with the line profiles observed by FUSE. We discuss the observational consequences of PRD, and show that this computational method also allows us to include PRD in modeling the continuum pumped H_2 fluorescence, treating about 4000 lines simultaneously.Comment: 17 pages, accepted for publication in Ap