Degree growth of meromorphic surface maps

We study the degree growth of iterates of meromorphic selfmaps of compact Kahler surfaces. Using cohomology classes on the Riemann-Zariski space we show that the degrees grow similarly to those of mappings that are algebraically stable on some birational model.Comment: 17 pages, final version, to appear in Duke Math Journa

Supersonic turbulence in 3D isothermal flow collision

Colliding supersonic bulk flows shape observable properties and internal physics of various astrophysical objects, like O-star winds, molecular clouds, galactic sheets, binaries, or gamma-ray bursts. Using numerical simulations, we show that the bulk flows leave a clear imprint on the collision zone, its mean properties and the turbulence it naturally develops. Our model setup consists of 3D head-on colliding isothermal hydrodynamical flows with Mach numbers between 2 and 43. Simulation results are in line with expectations from self-similarity: root mean square Mach numbers (Mrms) scale linearly with upstream Mach numbers, mean densities remain limited to a few times the upstream density. The density PDF is not log-normal. The turbulence is inhomogeneous: weaker in the zone center than close to the confining shocks. It is anisotropic: while Mrms is generally supersonic, Mrms transverse to the upstream flow is always subsonic. We argue that uniform, isothermal, head-on colliding flows generally disfavor isotropic, supersonic turbulence. The anisotropy carries over to other quantities like the density variance - Mach number relation. Structure functions differ depending on whether they are computed along a line-of-sight perpendicular or parallel to the upstream flow. We suggest that such line-of-sight effects should be kept in mind when interpreting turbulence characteristics derived from observations.Comment: 20 pages, 14 figures, 4 tables, accepted by Astronomy and Astrophysic

Tentative detection of ethylene glycol toward W51/e2 and G34.3+0.2

How complex organic - and potentially prebiotic - molecules are formed in regions of low- and high-mass star-formation remains a central question in astrochemistry. In particular, with just a few sources studied in detail, it is unclear what role environment plays in complex molecule formation. In this light, a comparison of relative abundances of related species between sources might be useful to explain observed differences. We seek to measure the relative abundance between three important complex organic molecules, ethylene glycol ((CH$_2$OH)$_2$), glycolaldehyde (CH$_2$OHCHO) and methyl formate (HCOOCH$_3$), toward high-mass protostars and thereby provide additional constraints on their formation pathways. We use IRAM 30-m single dish observations of the three species toward two high-mass star-forming regions - W51/e2 and G34.3+0.2 - and report a tentative detection of (CH2OH)2 toward both sources. Assuming that (CH$_2$OH)$_2$, CH$_2$OHCHO and HCOOCH$_3$ spatially coexist, relative abundance ratios, HCOOCH$_3$/(CH$_2$OH)$_2$, of 31 and 35 are derived for G34.3+0.2 and W51/e2, respectively. CH$_2$OHCHO is not detected, but the data provide lower limits to the HCOOCH$_3$/CH$_2$OHCHO abundance ratios of $\ge$193 for G34.3+0.2 and $\ge$550 for W51/e2. A comparison of these results to measurements from various sources in the literature indicates that the source luminosities may be correlated with the HCOOCH$_3$/(CH$_2$OH)$_2$ and HCOOCH$_3$/CH$_2$OHCHO ratios. This apparent correlation may be a consequence of the relative timescales each source spend at different temperatures-ranges in their evolution. Furthermore, we obtain lower limits to the ratio of (CH$_2$OH)$_2$/CH2OHCHO for G34.3+0.2 ($\ge$6) and W51/e2 ($\ge$16). This result confirms that a high (CH$_2$OH)$_2$/CH$_2$OHCHO abundance ratio is not a specific property of comets, as previously speculated.Comment: Accepted for publication by A&

The energetics of relativistic magnetic reconnection: ion-electron repartition and particle distribution hardness

Collisionless magnetic reconnection is a prime candidate to account for flare-like or steady emission, outflow launching, or plasma heating, in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas. But the fate of the initial magnetic energy in a reconnection event remains poorly known: what is the amount given to kinetic energy, the ion/electron repartition, and the hardness of the particle distributions? We explore these questions with 2D particle-in-cell simulations of ion-electron plasmas. We find that 45 to 75% of the total initial magnetic energy ends up in kinetic energy, this fraction increasing with the inflow magnetization. Depending on the guide field strength, ions get from 30 to 60% of the total kinetic energy. Particles can be separated into two populations that only weakly mix: (i) particles initially in the current sheet, heated by its initial tearing and subsequent contraction of the islands; and (ii) particles from the background plasma that primarily gain energy via the reconnection electric field when passing near the X-point. Particles (ii) tend to form a power-law with an index $p=-d\log n(\gamma)/d\log\gamma$, that depends mostly on the inflow Alfv\'en speed $V_A$ and magnetization $\sigma_s$ of species $s$, with for electrons $p=5$ to $1.2$ for increasing $\sigma_e$. The highest particle Lorentz factor, for ions or electrons, increases roughly linearly with time for all the relativistic simulations. This is faster, and the spectra can be harder, than for collisionless shock acceleration. We discuss applications to microquasar and AGN coronae, to extragalactic jets, and to radio lobes. We point out situations where effects such as Compton drag or pair creation are important.Comment: 15 pages, submitted to A&

Gait analysis of patients with knee osteoarthritis highlights a pathological mechanical pathway and provides a basis for therapeutic interventions

Knee osteoarthritis (OA) is a painful and incapacitating disease affecting a large portion of the elderly population, for which no cure exists. There is a critical need to enhance our understanding of OA pathogenesis, as a means to improve therapeutic options. Knee OA is a complex disease influenced by many factors, including the loading environment. Analysing knee biomechanics during walking - the primary cyclic load-bearing activity - is therefore particularly relevant. There is evidence of meaningful differences in the knee adduction moment, flexion moment and flexion angle during walking between non-OA individuals and patients with medial knee OA. Furthermore, these kinetic and kinematic gait variables have been associated with OA progression. Gait analysis provides the critical information needed to understand the role of ambulatory biomechanics in OA development, and to design therapeutic interventions. Multidisciplinary research is necessary to relate the biomechanical alterations to the structural and biological components of OA. Cite this article: Favre J, Jolles BM. Analysis of gait, knee biomechanics and the physiopathology of knee osteoarthritis in the development of therapeutic interventions. EFORT Open Rev 2016;1:368-374. DOI: 10.1302/2058-5241.1.000051

Apar-T: code, validation, and physical interpretation of particle-in-cell results

We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and non-linear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. Second, we detail a new method for initial loading of Maxwell-J\"uttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. Third, we scrutinize the question of what description of physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p~10^10 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter Lambda^{PIC}, which we show to scale as 1/p. We explore two implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p~10^10 smaller than in real plasmas. The other is that the level of electric field fluctuations scales as 1/Lambda^{PIC} ~ p. We provide a corresponding exact expression. Fourth, we compare the Vlasov-Maxwell theory, which describes a phase-space fluid with infinite Lambda, to the PIC model and its relatively small Lambda.Comment: 24 pages, 14 figures, accepted in Astronomy & Astrophysic

HCOOCH3 as a probe of temperature and structure of Orion-KL

We studied the O-bearing molecule HCOOCH3 to characterize the physical conditions of the different molecular source components in Orion-KL. We identify 28 methyl formate emission peaks throughout the 50" field of observations. The two strongest peaks are in the Compact Ridge (MF1) and in the SouthWest of the Hot Core (MF2). Spectral confusion is still prevailing as half of the expected transitions are blended over the region. Assuming that the transitions are thermalized, we derive the temperature at the five main emission peaks. At the MF1 position we find a temperature of 80K in a 1.8"x0.8" beam size and 120K on a larger scale (3.6" x2.2"), suggesting an external source of heating, whereas the temperature is about 130K at the MF2 position on both scales. Transitions of HCOOCH3 in vt=1 are detected as well and the good agreement of the positions on the rotational diagrams between the vt=0 and the vt=1 transitions suggests a similar temperature. The velocity of the gas is between 7.5 and 8.0km/s depending on the positions and column density peaks vary from 1.6x10^16 to 1.6x10^17cm^-2. A second velocity component is observed around 9-10 km/s in a North-South structure stretching from the Compact Ridge up to the BN object; this component is warmer at the MF1 peak. The two other C2H4O2 isomers are not detected and the derived upper limit for the column density is <3x10^14cm^-2 for glycolaldehyde and <2x10^15cm^-2 for acetic acid. From the 223GHz continuum map, we identify several dust clumps with associated gas masses in the range 0.8 to 5.8Msun. Assuming that the HCOOCH3 is spatially distributed as the dust, we find relative abundances of HCOOCH3 in the range <0.1x10^-8 to 5.2x10^-8. We suggest a relation between the methyl formate distribution and shocks as traced by 2.12 mum H2 emission.Comment: Accepted for publication in A&