Planetesimal disk evolution driven by embryo-planetesimal gravitational scattering

The process of gravitational scattering of planetesimals by a massive protoplanetary embryo is explored theoretically. We propose a method to describe the evolution of the disk surface density, eccentricity, and inclination caused by the embryo-planetesimal interaction. It relies on the analytical treatment of the scattering in two extreme regimes of the planetesimal epicyclic velocities: shear-dominated (dynamically ``cold'') and dispersion-dominated (dynamically ``hot''). In the former, planetesimal scattering can be treated as a deterministic process. In the latter, scattering is mostly weak because of the large relative velocities of interacting bodies. This allows one to use the Fokker-Planck approximation and the two-body approximation to explore the disk evolution. We compare the results obtained by this method with the outcomes of the direct numerical integrations of planetesimal orbits and they agree quite well. In the intermediate velocity regime an approximate treatment of the disk evolution is proposed based on interpolation between the two extreme regimes. We also calculate the rate of embryo's mass growth in an inhomogeneous planetesimal disk and demonstrate that it is in agreement with both the simulations and earlier calculations. Finally we discuss the question of the direction of the embryo-planetesimal interaction in the dispersion-dominated regime and demonstrate that it is repulsive. This means that the embryo always forms a gap in the disk around it, which is in contrast with the results of other authors. The machinery developed here will be applied to realistic protoplanetary systems in future papers.Comment: 40 pages, 9 figures, submitted to A

Predicting spinor condensate dynamics from simple principles

We study the spin dynamics of quasi-one-dimensional F=1 condensates both at zero and finite temperatures for arbitrary initial spin configurations. The rich dynamical evolution exhibited by these non-linear systems is explained by surprisingly simple principles: minimization of energy at zero temperature, and maximization of entropy at high temperature. Our analytical results for the homogeneous case are corroborated by numerical simulations for confined condensates in a wide variety of initial conditions. These predictions compare qualitatively well with recent experimental observations and can, therefore, serve as a guidance for on-going experiments.Comment: 4 pages, 2 figures. v3: matches version appeared in PR

On the Fairness of Centralised Decision-Making Strategies in multi-TSO Power Systems

In this paper, we consider an interconnected power system, where the different Transmission System Operators (TSOs) have agreed to transferring some of their competences to a Centralised Control Center (CCC). In such a context, a recurrent difficulty for the CCC is to define decision-making strategies which are fair enough to every TSO of the interconnected system. We address this multiobjective problem when the objective of every TSO can be represented by a real-valued function. We propose an algorithm to elect the solution that leads to the minimisation of the distance with the utopian minimum after having normalised the different objectives. We analyse the fairness of this solution in the sense of economics. We illustrate the approach with the IEEE 118 bus system partitioned in 3 areas having as local objective the minimisation of active power losses, the maximisation of reactive power reserves, or a combination of both criteria.multi-area power system, centralised control, multi-objective optimisation, fairness.

Good rotations

Numerical integrations in celestial mechanics often involve the repeated computation of a rotation with a constant angle. A direct evaluation of these rotations yields a linear drift of the distance to the origin. This is due to roundoff in the representation of the sine s and cosine c of the angle theta. In a computer, one generally gets c^2 + s^2 1, resulting in a mapping that is slightly contracting or expanding. In the present paper we present a method to find pairs of representable real numbers s and c such that c^2 + s^2 is as close to 1 as possible. We show that this results in a drastic decrease of the systematic error, making it negligible compared to the random error of other operations. We also verify that this approach gives good results in a realistic celestial mechanics integration.Comment: 24 pages, 3 figure

Spot distribution and fast surface evolution on Vega

Spectral signatures of surface spots were recently discovered from high cadence observations of the A star Vega. We aim at constraining the surface distribution of these photospheric inhomogeneities, and investigating a possible short term evolution of the spot pattern. Using data collected over five consecutive nights, we employ the Doppler Imaging method to reconstruct three different maps of the stellar surface, from three consecutive subsets of the whole time-series. The surface maps display a complex distribution of dark and bright spots, covering most of the visible fraction of the stellar surface. A number of surface features are consistently recovered in all three maps, but other features seem to evolve over the time span of observations, suggesting that fast changes can affect the surface of Vega within a few days at most. The short-term evolution is observed as emergence or disappearance of individual spots, and may also show up as zonal flows, with low-latitude and high latitude belts rotating faster than intermediate latitudes. It is tempting to relate the surface brightness activity to the complex magnetic field topology previously reconstructed for Vega, although strictly simultaneous brightness and magnetic maps will be necessary to assess this potential link.Comment: MNRAS Letters, in pres

A polarity reversal in the large-scale magnetic field of the rapidly rotating Sun HD 190771

Aims. We investigate the long-term evolution of the large-scale photospheric magnetic field geometry of the solar-type star HD 190771. With fundamental parameters very close to those of the Sun except for a shorter rotation period of 8.8 d, HD 190771 provides us with a first insight into the specific impact of the rotation rate in the dynamo generation of magnetic fields in 1 $M_\odot$ stars. Methods. We use circularly polarized, high-resolution spectra obtained with the NARVAL spectropolarimeter (Observatoire du Pic du Midi, France) and compute cross-correlation line profiles with high signal-to-noise ratio to detect polarized Zeeman signatures. From three phase-resolved data sets collected during the summers of 2007, 2008, and 2009, we model the large-scale photospheric magnetic field of the star by means of Zeeman-Doppler imaging and follow its temporal evolution. Results. The comparison of the magnetic maps shows that a polarity reversal of the axisymmetric component of the large-scale magnetic field occurred between 2007 and 2008, this evolution being observed in both the poloidal and toroidal magnetic components. Between 2008 and 2009, another type of global evolution occured, characterized by a sharp decrease of the fraction of magnetic energy stored in the toroidal component. These changes were not accompanied by significant evolution in the total photospheric magnetic energy. Using our spectra to perform radial velocity measurements, we also detect a very low-mass stellar companion to HD 190771.Comment: Accepted by Astronomy and Astrophysics (Letter to the Editor

Transient behavior of surface plasmon polaritons scattered at a subwavelength groove

We present a numerical study and analytical model of the optical near-field diffracted in the vicinity of subwavelength grooves milled in silver surfaces. The Green's tensor approach permits computation of the phase and amplitude dependence of the diffracted wave as a function of the groove geometry. It is shown that the field diffracted along the interface by the groove is equivalent to replacing the groove by an oscillating dipolar line source. An analytic expression is derived from the Green's function formalism, that reproduces well the asymptotic surface plasmon polariton (SPP) wave as well as the transient surface wave in the near-zone close to the groove. The agreement between this model and the full simulation is very good, showing that the transient "near-zone" regime does not depend on the precise shape of the groove. Finally, it is shown that a composite diffractive evanescent wave model that includes the asymptotic SPP can describe the wavelength evolution in this transient near-zone. Such a semi-analytical model may be useful for the design and optimization of more elaborate photonic circuits whose behavior in large part will be controlled by surface waves.Comment: 12 pages, 10 figure