ALMA CO J=6-5 observations of IRAS16293-2422: Shocks and entrainment

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, IRAS16293-2422 (I16293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 9, which includes CO J=6-5 (E_up/k_B ~ 116 K), at an unprecedented spatial resolution (~0.2", 25 AU). I16293 itself is composed of two sources, A and B, with a projected distance of 5". CO J=6-5 emission is detected throughout the region, particularly in small, arcsecond-sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment takes place over large spatial scales, >~ 100 AU, and in small discrete knots. This unique dataset shows that the combination of a high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the warm inner environment of low-mass protostars.Comment: Accepted for publication in A&A Letter

Modelling diverse root density dynamics and deep nitrogen uptake — a simple approach

We present a 2-D model for simulation of root density and plant nitrogen (N) uptake for crops grown in agricultural systems, based on a modification of the root density equation originally proposed by Gerwitz and Page in J Appl Ecol 11:773–781, (1974). A root system form parameter was introduced to describe the distribution of root length vertically and horizontally in the soil profile. The form parameter can vary from 0 where root density is evenly distributed through the soil profile, to 8 where practically all roots are found near the surface. The root model has other components describing root features, such as specific root length and plant N uptake kinetics. The same approach is used to distribute root length horizontally, allowing simulation of root growth and plant N uptake in row crops. The rooting depth penetration rate and depth distribution of root density were found to be the most important parameters controlling crop N uptake from deeper soil layers. The validity of the root distribution model was tested with field data for white cabbage, red beet, and leek. The model was able to simulate very different root distributions, but it was not able to simulate increasing root density with depth as seen in the experimental results for white cabbage. The model was able to simulate N depletion in different soil layers in two field studies. One included vegetable crops with very different rooting depths and the other compared effects of spring wheat and winter wheat. In both experiments variation in spring soil N availability and depth distribution was varied by the use of cover crops. This shows the model sensitivity to the form parameter value and the ability of the model to reproduce N depletion in soil layers. This work shows that the relatively simple root model developed, driven by degree days and simulated crop growth, can be used to simulate crop soil N uptake and depletion appropriately in low N input crop production systems, with a requirement of few measured parameters

Influence of Pure Dephasing on Emission Spectra from Single Photon Sources

We investigate the light-matter interaction of a quantum dot with the electromagnetic field in a lossy microcavity and calculate emission spectra for non-zero detuning and dephasing. It is found that dephasing shifts the intensity of the emission peaks for non-zero detuning. We investigate the characteristics of this intensity shifting effect and offer it as an explanation for the non-vanishing emission peaks at the cavity frequency found in recent experimental work.Comment: Published version, minor change

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: II. Distribution and origin of warm molecular gas

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.Comment: 18 pages, accepted by A&A, A version with the figures in higher quality can be found on my website: http://www.cfa.harvard.edu/~tvankemp

Imaging Oxygen Distribution in Marine Sediments. The Importance of Bioturbation and Sediment Heterogeneity

The influence of sediment oxygen heterogeneity, due to bioturbation, on diffusive oxygen flux was investigated. Laboratory experiments were carried out with 3 macrobenthic species presenting different bioturbation behaviour patterns:the polychaetes Nereis diversicolor and Nereis virens, both constructing ventilated galleries in the sediment column, and the gastropod Cyclope neritea, a burrowing species which does not build any structure. Oxygen two-dimensional distribution in sediments was quantified by means of the optical planar optode technique. Diffusive oxygen fluxes (mean and integrated) and a variability index were calculated on the captured oxygen images. All species increased sediment oxygen heterogeneity compared to the controls without animals. This was particularly noticeable with the polychaetes because of the construction of more or less complex burrows. Integrated diffusive oxygen flux increased with oxygen heterogeneity due to the production of interface available for solute exchanges between overlying water and sediments. This work shows that sediment heterogeneity is an important feature of the control of oxygen exchanges at the sediment–water interface

Fractional differentiability for solutions of nonlinear elliptic equations

We study nonlinear elliptic equations in divergence form $${\operatorname{div}}{\mathcal A}(x,Du)={\operatorname{div}}G.$$ When ${\mathcal A}$ has linear growth in $Du$, and assuming that $x\mapsto{\mathcal A}(x,\xi)$ enjoys $B^\alpha_{\frac{n}\alpha, q}$ smoothness, local well-posedness is found in $B^\alpha_{p,q}$ for certain values of $p\in[2,\frac{n}{\alpha})$ and $q\in[1,\infty]$. In the particular case ${\mathcal A}(x,\xi)=A(x)\xi$, $G=0$ and $A\in B^\alpha_{\frac{n}\alpha,q}$, $1\leq q\leq\infty$, we obtain $Du\in B^\alpha_{p,q}$ for each $p<\frac{n}\alpha$. Our main tool in the proof is a more general result, that holds also if ${\mathcal A}$ has growth $s-1$ in $Du$, $2\leq s\leq n$, and asserts local well-posedness in $L^q$ for each $q>s$, provided that $x\mapsto{\mathcal A}(x,\xi)$ satisfies a locally uniform $VMO$ condition

Shock excitation of H2_2 in the James Webb Space Telescope era

(Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 excitation using such models, and to test over which parameter space distinct signatures are produced in H2 emission. We present simulated H2 emission using the Paris-Durham shock code over an extensive grid of 14,000 plane-parallel stationary shock models, a large subset of which are exposed to an external UV radiation field. The grid samples 6 input parameters: preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, cosmic-ray-ionization rate, and PAH abundance. Physical quantities, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. The strength of the transverse magnetic field, set by the scaling factor, b, plays a key role in the excitation of H2. At low values of b (<~ 0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; at higher values (b >~ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b >= 1 can be spatially resolved with JWST for nearby objects. When the input kinetic energy flux increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission.Comment: Published in A&