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.

We present first results of a long-term project which aims at multi-scale, multi-physics simulations of wind accretion in microquasars and high-mass X-ray binaries. The 3D hydrodynamical simulations cover all scales, from the circum-binary environment down to the immediate vicinity of the black hole. We first introduce the numerical method and parallelization strategy of the AMR A-MAZE code. We then discuss some preliminary results of how, and on what scales, an accretion disk is formed around the black hole. We finally present some characteristics of this disk, which is far from Keplerian. We emphasize that on all scales shocks play a decisive role for the accretion process and the process of structure formation -- for the formation of the large scale, nearly coherent structure of the disk, but also for the formation of turbulent fluctuations.

Magnetic reconnection is a leading mechanism for magnetic energy conversion and high-energy non-thermal particle production in a variety of high-energy astrophysical objects, including ones with relativistic ion-electron plasmas (e.g., microquasars or AGNs) - a regime where first principle studies are scarce. We present 2D particle-in-cell (PIC) simulations of low $\beta$ ion-electron plasmas under relativistic conditions, i.e., with inflow magnetic energy exceeding the plasma rest-mass energy. We identify outstanding properties: (i) For relativistic inflow magnetizations (here $10 < \sigma_e < 360$), the reconnection outflows are dominated by thermal agitation instead of bulk kinetic energy. (ii) At large inflow electron magnetization ($\sigma_e > 80$), the reconnection electric field is sustained more by bulk inertia than by thermal inertia. It challenges the thermal-inertia-paradigm and its implications. (iii) The inflows feature sharp transitions at the entrance of the diffusion zones. These are not shocks but results from particle ballistic motions, all bouncing at the same location, provided that the thermal velocity in the inflow is far smaller than the inflow E cross B bulk velocity. (iv) Island centers are magnetically isolated from the rest of the flow, and can present a density depletion at their center. (v) The reconnection rates are slightly larger than in non-relativistic studies. They are best normalized by the inflow relativistic Alfv\'en speed projected in the outflow direction, which then leads to rates in a close range (0.14-0.25) thus allowing for an easy estimation of the reconnection electric field.

The PIC model relies on two building blocks. The first stems from the capability of computers to handle only up to $\sim10^{10}$ particles, while real plasmas contain from $10^4$ to $10^{20}$ particles per Debye sphere: a coarse-graining step must be used, whereby of the order of $p\sim10^{10}$ real particles are represented by a single computer superparticle. The second is field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e. physical quantities depending on the number $p$. They all derive from the plasma parameter $\Lambda$, which we show to be proportional to $1/p$. We explore three examples: the rapid collision- and fluctuation-induced thermalization of plasmas with different temperatures, that scale with the number of superparticles per grid cell and are a factor $p\sim10^{10}$ faster than in real plasmas; the high level of electrostatic fluctuations in a thermal plasma, with corrections due to the finite superparticle sizes; and the blurring of the linear spectrum of the filamentation instability, where the fastest growing modes do not dominate the total energy because of a high level of fluctuations. We stress that the enhanced collisions and correlations of PIC plasmas must be kept negligible toward kinetic physics.

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.

We report on the development of an implicit multi-D hydrodynamic code for stellar evolution. We present two test-cases relevant for the first scientific goal of the code: the simulation of convection in pulsating stars. First results on a realistic stellar model are also presented.

We revisit the tidal stability of extrasolar systems harboring a transiting planet and demonstrate that, independently of any tidal model, none but one (HAT-P-2b) of these planets has a tidal equilibrium state, which implies ultimately a collision of these objects with their host star. Consequently, conventional circularization and synchronization timescales cannot be defined because the corresponding states do not represent the endpoint of the tidal evolution. Using numerical simulations of the coupled tidal equations for the spin and orbital parameters of each transiting planetary system, we confirm these predictions and show that the orbital eccentricity and the stellar obliquity do not follow the usually assumed exponential relaxation but instead decrease significantly, reaching eventually a zero value, only during the final runaway merging of the planet with the star. The only characteristic evolution timescale of {\it all} rotational and orbital parameters is the lifetime of the system, which crucially depends on the magnitude of tidal dissipation within the star. These results imply that the nearly circular orbits of transiting planets and the alignment between the stellar spin axis and the planetary orbit are unlikely to be due to tidal dissipation. Other dissipative mechanisms, for instance interactions with the protoplanetary disk, must be invoked to explain these properties.

(Abridged) We calculate the full stellar-structural evolution of donors in AM CVn systems formed through the WD channel coupled to the binary's evolution. Contrary to assumptions made in prior modelling, these donors are not fully convective over much of the AM CVn phase and do not evolve adiabatically under mass loss indefinitely. Instead, we identify three distinct phases of evolution: a mass transfer turn-on phase (during which the orbital period continues to decrease even after contact, the donor contracts, and the mass transfer rate accelerates to its maximum), a phase in which the donor expands adiabatically in response to mass loss, and a cooling phase beginning at orbital periods of approximately 45--55 minutes during which the donor contracts. The physics that determines the behaviour in the first and third phases, both of which are new outcomes of this study, are discussed in some detail. We find the overall duration of the turn-on phase to be between $\sim 10^4$-$10^6$ yrs, significantly longer than prior estimates. We predict the donor's luminosity and effective temperature. During the adiabatic expansion phase (ignoring irradiation effects), the luminosity is approximately $10^{-6}$--$10^{-4} L_\odot$ and the effective temperature is approximately 1000--1800 K. However, the flux generated in the accretion flow dominates the donor's intrinsic light at all times. The impact of irradiation on the donor extends the phase of adiabatic expansion to longer orbital periods and alters the donor's observational characteristics. Irradiated donors during the adiabatic phase can attain a surface luminosity of up to $\approx10^{-2} L_\odot$. We argue that the turn-on and cooling phases both will leave significant imprints on the AM CVn population's orbital period distribution.

We analyse the presence of nonradial oscillations in Cepheids, a problem which has not been theoretically revised since the work of Dziembowsky (1977) and Osaki (1977). Our analysis is motivated by a work of Moskalik et al. (2004) which reports the detec tion of low amplitude periodicities in a few Cepheids of the large Magellanic cloud. These newly discovered periodicities were interpreted as nonradial modes.} {Based on linear nonadiabatic stability analysis, our goal is to reanalyse the presence and stability of nonradial modes, taking into account improvement in the main input phys ics required for the modelling of Cepheids.} {We compare the results obtained from two different numerical methods used to solve the set of differential equations: a matrix method and the Ricatti method.} {We show the limitation of the matrix method to find low order p-modes ($l<6$), because of their dual character in evolved stars such as Cepheids. For higher order p-modes, we find an excellent agreement between the two methods.} {No nonradial instability is found below $l=5$, whereas many unstable nonradial modes exist for higher orders. We also find that nonradial modes remain unstable, even at hotter effective temperatures than the blue edge of the Cepheid instability strip, where no radial pulsations are expected.

We discuss results of our study on AM CVn binaries formed with donors that never ignited He before contact. For the first time, we treat the donor's in these systems in the context of a full stellar structure evolution theory and find that the binary's evolution can described in terms of 3 phases: contact, adiabatic donor expansion, and late-time donor cooling. Details of the first and third phase are new results from this study and we focus on generally characterizing these two phases. Finally, we present our predictions for the donor's light in these systems.

We present a new model of giant planet formation that extends the core-accretion model of Pollack etal (1996) to include migration, disc evolution and gap formation. We show that taking into account these effects can lead to a much more rapid formation of giant planets, making it compatible with the typical disc lifetimes inferred from observations of young circumstellar discs. This speed up is due to the fact that migration prevents the severe depletion of the feeding zone as observed in in situ calculations. Hence, the growing planet is never isolated and it can reach cross-over mass on a much shorter timescale. To illustrate the range of planets that can form in our model, we describe a set of simulations in which we have varied some of the initial parameters and compare the final masses and semi-major axes with those inferred from observed extra-solar planets.

We present a semi-analytical free-energy model aimed at characterizing the thermodynamic properties of dense fluid helium, from the low-density atomic phase to the high-density fully ionized regime. The model is based on a free-energy minimization method and includes various different contributions representative of the correlations between atomic and ionic species and electrons. This model allows the computation of the thermodynamic properties of dense helium over an extended range of density and temperature and leads to the computation of the phase diagram of dense fluid helium, with its various temperature and pressure ionization contours. One of the predictions of the model is that pressure ionization occurs abruptly at $\rho \simgr 10$ g cm$^{-3}$, {\it i.e.} $P\simgr 20$ Mbar, from atomic helium He to fully ionized helium He$^{2+}$, or at least to a strongly ionized state, without He$^{+}$ stage, except at high enough temperature for temperature ionization to become dominant. These predictions and this phase diagram provide a guide for future dynamical experiments or numerical first-principle calculations aimed at studying the properties of helium at very high density, in particular its metallization. Indeed, the characterization of the helium phase diagram bears important consequences for the thermodynamic, magnetic and transport properties of cool and dense astrophysical objects, among which the solar and the numerous recently discovered extrasolar giant planets.

We derive a new formulation to calculate the excess chemical potential of a fraction of $N_1$ particles interacting with $N_2$ particles of a different species. The excess chemical potential is calculated numerically from first principles by coupling molecular dynamics and Thomas-Fermi density functional theory to take into account the contribution arising from the quantum electrons on the forces acting on the ions. The choice of this simple functional is motivated by the fact that the present paper is devoted to the derivation and the validation of the method but more complicated functionals can and will be implemented in the future. This new method is applied in the microcanonical ensemble, the most natural ensemble for molecular dynamics simulations. This avoids the introduction of a thermostat in the simulation, and thus uncontroled modifications of the trajectories calculated from the forces between particles. The calculations are conducted for three values of the input thermodynamic quantities, energy and density, and for different total numbers of particles in order to examine the uncertainties due to finite size effects. This method and these calculations lie the basic foundation to study the thermodynamic stability of dense mixtures, without any {\it a priori} assumption on the degree of ionization of the different species.