The slow rotation of some young stars and the extreme rotation periods of some Ap stars have so far defied explanation. The absence of sufficiently efficient braking mechanisms for newly formed stars points to the star formation process itself as the origin. I find that a mode of star formation exists by which a protostar can form without accreting angular momentum. It depends on the survival of a magnetic connection between the accreting matter and the birth cloud. The conditions for this process to operate are analyzed, and illustrated with a generic (scale-free) model. Depending on the initial rotation rate of the accreting matter, either a Keplerian disk forms, or the gas ends up rotating with the rotation period of the cloud, even if it is at a large distance. The boundary in parameter space between these two outcomes is sharp.

The processes leading deformation and destruction of planets spiraling into the convective envelope of their host stars are described. The planet is compressed by the ram pressure, and deformed into a flattened shape for which a quantitative model is developed. Compression increases the planet's density contrast with the envelope and its gravitational binding energy. This increases the survivability especially of gas planets. An estimate is given for the depth of disruption by ram pressure, and for the subsequent fragmentation of the remnants. We show how the debris of rocky or iron planets, instead of mixing through the convection zone, sinks below the base of the convection zone. The time scale of the entire sequence of events is of the order of a few orbital times of the planet. If spiral-in of (partly) icy, rocky or iron planets has happened to the pre-main sequence Sun, it could account for the higher opacity below the base of the convection zone as inferred from helioseismology.

We review the current state of knowledge of magnetic fields inside stars, concentrating on recent developments concerning magnetic fields in stably stratified (zones of) stars, leaving out convective dynamo theories and observations of convective envelopes. We include the observational properties of A, B and O-type main-sequence stars, which have radiative envelopes, and the fossil field model which is normally invoked to explain the strong fields sometimes seen in these stars. Observations seem to show that Ap-type stable fields are excluded in stars with convective envelopes. Most stars contain both radiative and convective zones, and there are potentially important effects arising from the interaction of magnetic fields at the boundaries between them. Related to this, we discuss whether the Sun could harbour a magnetic field in its core. Recent developments regarding the various convective and radiative layers near the surfaces of early-type stars and their observational effects are examined. We look at possible dynamo mechanisms that run on differential rotation rather than convection. Finally we turn to neutron stars with a discussion of the possible origins for their magnetic fields.

We show that the angular momentum exchange mechanism governing the evolution of mass transferring binary stars does not apply to Roche-lobe filling planets, because most of the angular momentum of the mass transferring stream is absorbed by the host star. Apart from a correction for the difference in specific angular momentum of the stream and the centre of mass of the planet, the orbit does not expand much on Roche-lobe overflow. We explore the conditions for dynamically unstable Roche-lobe overflow as a function of planet mass and mass and radius (age) of host star and equation of state of planet. For a Sun-like host, gas giant planets in a range of mass and entropy can undergo dynamical mass transfer. Examples of the evolution of the mass transfer process are given. Dynamic mass transfer of rocky planets depend somewhat sensitively on equation of state used. Silicate planets in the range $1 M_{\oplus} <M_{\mathrm{p}} < 10 M_{\oplus} $ typically go through a phase of dynamical mass transfer before settling to slow overflow when their mass drops to less than $1 M_{\oplus}$.

This text is intended as an introduction to magnetohydrodynamics in astrophysics, emphasizing a fast path to the elements essential for physical understanding. It assumes experience with concepts from fluid mechanics: the fluid equation of motion and the Lagrangian and Eulerian descriptions of fluid flow. In addition, the basics of vector calculus and elementary special relativity are needed. Not much knowledge of electromagnetic theory is required. In fact, since MHD is much closer in spirit to fluid mechanics than to electromagnetism, an important part of the learning curve is to overcome intuitions based on the vacuum electrodynamics of one's high school days. The first chapter (only 39 pp) is meant as a practical introduction including exercises. This is the `essential' part. The exercises are important as illustrations of the points made in the text (especially the less intuitive ones). Almost all are mathematically unchallenging. The supplement in chapter 2 contains further explanations, more specialized topics and connections to the occasional topic somewhat outside MHD. The emphasis is on physical understanding by the visualization of MHD processes, as opposed to more formal approaches.

Helium burning in the convective cores of horizontal branch and `red clump' stars appears to involve a process of `ingestion' of unburnt helium into the core, the physics of which has not been identified yet. I show here that a limiting factor controlling the growth is the buoyancy of helium entering the denser C+O core. It yields a growth rate which scales directly with the convective luminosity of the core, and agrees with constraints on core size from current asteroseismology.

The process of convective settling is based on the assumption that a small fraction of the low-entropy downflows sink from the photosphere down to the bottom of the star's envelope convection zone retaining a substantial entropy contrast. We have previously shown that this process could explain the slow Li depletion observed in the Sun. We construct a parametric model of convective settling to investigate the dependence of Li and Be depletion on stellar mass and age. Our model is generally in good agreement with the Li abundances measured in open clusters and solar twins, although it seems to underestimate the Li depletion in the first ~1 Gyr. The model is also compatible with the Be abundances measured in a sample of field stars.

The question of possible small-scale dynamo action in the surface layers of the Sun is revisited with realistic 3D MHD simulations. As in other MHD problems, dynamo action is found to be a sensitive function of the magnetic Prandtl number ${\rm P_{\rm m} }=\nu/\eta$; it disappears below a critical value ${\rm P_{\rm c}}$ which is a function of the numerical resolution. At a grid spacing of 3.5 km, ${\rm P_{\rm c}}$ based on the hyperdiffusivities implemented in the code (STAGGER) is $\approx 1$, increasing with increasing grid spacing. As in other settings, it remains uncertain whether small scale dynamo action is present in the astrophysical limit where ${\rm P_{\rm m} }<<1$ and magnetic Reynolds number ${\rm R_m}\gg 1$. The question is discussed in the context of the strong effect that external stray fields are observed to have in generating and maintaining dynamo action in other numerical and laboratory systems, and in connection with the type-II hypertransient behavior of dynamo action observed in the absence of such external fields.

On the long nuclear time scale of stellar main-sequence evolution, even weak mixing processes can become relevant for redistributing chemical species in a star. We investigate a process of "differential heating," which occurs when a temperature fluctuation propagates by radiative diffusion from the boundary of a convection zone into the adjacent radiative zone. The resulting perturbation of the hydrostatic equilibrium causes a flow that extends some distance from the convection zone. We study a simplified differential-heating problem with a static temperature fluctuation imposed on a solid boundary. The astrophysically relevant limit of a high Reynolds number and a low P\'eclet number (high thermal diffusivity) turns out to be interestingly non-intuitive. We derive a set of scaling relations for the stationary differential heating flow. A numerical method adapted to a high dynamic range in flow amplitude needed to detect weak flows is presented. Our two-dimensional simulations show that the flow reaches a stationary state and confirm the analytic scaling relations. These imply that the flow speed drops abruptly to a negligible value at a finite height above the source of heating. We approximate the mixing rate due to the differential heating flow in a star by a height-dependent diffusion coefficient and show that this mixing extends about $4\%$ of the pressure scale height above the convective core of a $10\,M_\odot$ zero-age main sequence star.

In stably stratified stars, numerical magneto-hydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfv\'en times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

The net effect of the small scale magnetic field on the Sun's (bolometric) brightness is studied with realistic 3D MHD simulations. The direct effect of brightening within the magnetic field itself is consistent with measurements in high-resolution observations. The high 'photometric accuracy' of the simulations, however, reveal compensating brightness effects that are hard to detect observationally. The influence of magnetic concentrations on the surrounding nonmagnetic convective flows (a 'proximity effect') reduces the brightness by an amount exceeding the brightening by the magnetic concentrations themselves. The net photospheric effect of the small scale field (~ -0.34% at a mean flux density of 50 G) is thus negative. We conclude that the main contribution to the observed positive correlation between the magnetic field and total solar irradiance must be magnetic dissipation in layers around the temperature minimum and above (not included in the simulations). This agrees with existing inferences from observations.

We study a process of slow mixing in stars with convective envelopes, which is driven by the settling of cool downward plumes below the base of the convection zone. If a small fraction (of order $10^{-7}$) of the material cooled at the surface retains a significant entropy deficit while descending in plumes, it can reach the depth where lithium burning takes place. The model calculates the thermal response and mixing below the convection zone due to the settling process, assuming that the plumes arrive at the base of the convection zone with a broad range of entropy contrasts. We obtain a good fit to the observed lithium depletion in the Sun by assuming that the settling mass flux is distributed with respect to the entropy contrast as a power law with a slope around -2. We find convective settling to have a negligible influence on the stratification below the convection zone, although mixing induced by it could modify the gradient of helium concentration.

It is now believed that magnetohydrodynamic equilibria can exist in stably stratified stars due to the seminal works of Braithwaite & Spruit (2004) and Braithwaite & Nordlund (2006). What is still not known is whether magnetohydrodynamic equilibria can exist in a barotropic star, in which stable stratification is not present. It has been conjectured by Reisenegger (2009) that there will likely not exist any magnetohydrodynamical equilibria in barotropic stars. We aim to test this claim by presenting preliminary MHD simulations of barotropic stars using the three dimensional stagger code of Nordlund & Galsgaard (1995).

A grid of numerical simulations of double-diffusive convection is presented for the astrophysical case where viscosity (Prandtl number Pr) and solute diffusivity (Lewis number Le) are much smaller than the thermal diffusivity. As in laboratory and geophysical cases convection takes place in a layered form. The proper translation between subsonic flows in a stellar interior and an incompressible (Boussinesq) fluid is given, and the validity of the Boussinesq approximation for the semiconvection problem is checked by comparison with fully compressible simulations. The predictions of a simplified theory of mixing in semiconvection given in a companion paper are tested against the numerical results, and used to extrapolate these to astrophysical conditions. The predicted effective He-diffusion coefficient is nearly independent of the double-diffusive layering thickness $d$. For a fiducial main sequence model (15 $M_\odot$) the inferred mixing time scale is of the order $10^{10}$ yr. An estimate for the secular increase of $d$ during the semiconvective phase is given. It can potentially reach a significant fraction of a pressure scale height.

A model is developed for the transport of heat and solute in a system of double-diffusive layers under astrophysical conditions (viscosity and solute diffusivity low compared with the thermal diffusivity). The process of formation of the layers is not part of the model but, as observed in geophysical and laboratory settings, is assumed to be fast compared to the life time of the semiconvective zone. The thickness of the layers is \tbf{a} free parameter of the model. When the energy flux of the star is specified, the effective semiconvective diffusivities are only weakly dependent on this parameter. An estimate is given of the evolution of layer thickness with time in a semiconvective zone. The model predicts that the density ratio has a maximum for which a stationary layered state can exist, $R_\rho\la \mr{Le}^{-1/2}$. Comparison of the model predictions with a grid of numerical simulations is presented in a companion paper.

Large scale magnetic field threading an accretion disk is a key ingredient in the jet formation model. The most attractive scenario for the origin of such a large scale field is the advection of the field by the gas in the accretion disk from the interstellar medium or a companion star. However, it is realized that outward diffusion of the accreted field is fast compared to the inward accretion velocity in a geometrically thin accretion disk if the value of the Prandtl number Pm is around unity. In this work, we revisit this problem considering the angular momentum of the disk is removed predominantly by the magnetically driven outflows. The radial velocity of the disk is significantly increased due to the presence of the outflows. Using a simplified model for the vertical disk structure, we find that even moderately weak fields can cause sufficient angular momentum loss via a magnetic wind to balance outward diffusion. There are two equilibrium points, one at low field strengths corresponding to a plasma-beta at the midplane of order several hundred, and one for strong accreted fields, beta~1. We surmise that the first is relevant for the accretion of weak, possibly external, fields through the outer parts of the disk, while the latter one could explain the tendency, observed in full 3D numerical simulations, of strong flux bundles at the centers of disk to stay confined in spite of strong MRI turbulence surrounding them.

Use of the Hirsch-index ($h$) as measure of an author's visibility in the scientific literature has become popular as an alternative to a gross measure like total citations (c). I show that, at least in astrophysics, $h$ correlates tightly with overall citations. The mean relation is $h=0.5(\sqrt c+1)$. Outliers are few and not too far from the mean, especially if `normalized' ADS citations are used for $c$ and $h$. Whatever the theoretical reasoning behind it, the Hirsch index in practice does not appear to measure something significantly new.

We show that discs accreting onto the magnetosphere of a rotating star can end up in a 'trapped' state, in which the inner edge of the disc stays near the corotation radius, even at low and varying accretion rates. The accretion in these trapped states can be steady or cyclic; we explore these states over wide range of parameter space. We find two distinct regions of instability, one related to the buildup and release of mass in the disk outside corotation, the other to mass storage within the transition region near corotation. With a set of calculations over long time scales we show how trapped states evolve from both nonaccreting and fully accreting initial conditions, and also calculate the effects of cyclic accretion on the spin evolution of the star. Observations of cycles such as found here would provide important clues on the physics of magnetospheric accretion. Recent observations of cyclic and other unusual variability in T Tauri stars (EXors) and X-ray binaries are discussed in this context.

We hypothesize that hot Jupiters with inflated sizes represent a separate planet formation channel,the merging of two low-mass stars. We show that the abundance and properties of W UMa stars and low mass detached binaries are consistent with their being possible progenitors. The degree of inflation of the transiting hot Jupiters correlates with their expected spiral-in life time by tidal dissipation, and this could indicate youth if the stellar dissipation parameter Q'* is sufficiently low. Several Jupiter-mass planets can form in the massive compact disk formed in a merger event. Gravitational scattering between them can explain the high incidence of excentric, inclined, and retrograde orbits. If the population of inflated planets is indeed formed by a merger process, their frequency should be much higher around blue stragglers than around T Tauri stars.

In addition to the `facular' brightening of active regions, the quiet Sun also contains a small scale magnetic field with associated brightenings in continuum radiation. We measure this contribution of quiet regions to the Sun's brightness from high spatial resolution (0"16-0"32) observations of the Swedish 1-m Solar Telescope (SST) and Hinode satellite. The line-of-sight magnetic field and continuum intensity near \ion{Fe}{i} 6302.5 \AA\ are used to quantify the correlation between field strength and brightness. The data show that magnetic flux density contains a significant amount of intrinsically weak fields that contribute little to brightness. We show that with data of high spatial resolution a calibration of magnetic flux density as a proxy for brightness excess is possible. In the SST data, the magnetic brightening of a quiet region with an average (unsigned) flux density of 10 G is about 0.15%. In the Hinode data, and in SST data reduced to Hinode resolution, the measured brightening is some 40% lower. With appropriate correction for resolution, magnetic flux density can be used as a reliable proxy in regions of small scale mixed polarity. The measured brightness effect is larger than the variation of irradiance over a solar cycle. It is not clear, however, if this quiet Sun contribution actually varies significantly. \keywords{Sun: surface magnetism -- photosphere -- solar-terrestrial relations}

We investigate the evolution of a thin viscous disc surrounding magnetic star, including the spindown of the star by the magnetic torques it exerts on the disc. The transition from an accreting to a non-accreting state, and the change of the magnetic torque across the corotation radius are included in a generic way, the widths of the transition taken in the range suggested by numerical simulations. In addition to the standard accreting state, two more are found. An accreting state can develop into a 'dead' disc state, with inner edge well outside corotation. More often, a 'trapped' state develops, in which the inner disc edge stays close to corotation even at very low accretion rates. The long-term evolution of these two states is different. In the dead state the star spins down incompletely, retaining much of its initial spin. In the trapped state the star asymptotically can spin down to arbitarily low rates, its angular momentum transferred to the disc. We identify these outcomes with respectively the rapidly rotating and the very slowly rotating classes of Ap stars and magnetic white dwarfs.

A grid of numerical simulations of double-diffusive convection is presented for astrophysical conditions. As in laboratory and geophysical cases convection takes place in a layered form. A translation between the astrophysical fluid mechanics and incompressible (Boussinesq) approximation is given, valid for thin layers. Its validity is checked by comparison of the results of fully compressible and Boussinesq simulations of semiconvection. A fitting formula is given for the superadiabatic gradient as a function of this parameter. The superadiabaticity depends on the thickness $d$ of the double diffusive layers, for which no good theory is available, but the effective He-diffusion coefficient is nearly independent of $d$. For a fiducial main sequence model (15 $M_\odot$) the inferred mixing time scale is of the order $10^{10}$ yr.

Using simultaneous observations in X-rays and optical, we have performed a homogeneous analysis of the cross-correlation behaviours of four X-ray binaries: SWIFT J1753.5-0127, GX 339-4, Sco X-1, and Cyg X-2. With high time-resolution observations using ULTRACAM and RXTE, we concentrate on the short time-scale, dt<20 s, variability in these sources. Here we present our database of observations, with three simultaneous energy bands in both the optical and the X-ray, and multiple epochs of observation for each source, all with ~second or better time resolution. For the first time, we include a dynamical cross-correlation analysis, i.e., an investigation of how the cross-correlation function changes within an observation. We describe a number of trends which emerge. We include the full dataset of results, and pick a few striking relationships from among them for further discussion. We find, that the surprising form of X-ray/optical cross-correlation functions, a positive correlation signal preceded by an anti-correlation signal, is seen in all the sources at least some of the time. Such behaviour suggests a mechanism other than reprocessing as being the dominant driver of the short-term variability in the optical emission. This behaviour appears more pronounced when the X-ray spectrum is hard. Furthermore, we find that the cross-correlation relationships themselves are not stable in time, but vary significantly in strength and form. This all hints at dynamic interactions between the emitting components which could be modelled through non-linear or differential relationships.

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

Observations with the 1-m Swedish Solar Telescope of the flows seen in penumbral filaments are presented. Time sequences of bright filaments show overturning motions strikingly similar to those seen along the walls of small isolated structures in the active regions. The filaments show outward propagating striations with inclination angles suggesting that they are aligned with the local magnetic field. We interpret it as the equivalent of the striations seen in the walls of small isolated magnetic structures. Their origin is then a corrugation of the boundary between an overturning convective flow inside the filament and the magnetic field wrapping around it. The outward propagation is a combination of a pattern motion due to the downflow observed along the sides of bright filaments, and the Evershed flow. The observed short wavelength of the striation argues against the existence of a dynamically significant horizontal field inside the bright filaments. Its intensity contrast is explained by the same physical effect that causes the dark cores of filaments, light bridges and `canals'. In this way striation represents an important clue to the physics of penumbral structure and its relation with other magnetic structures on the solar surface. We put this in perspective with results from the recent 3-D radiative hydrodynamic simulations.

In this lecture the basic theory of accretion disks is reviewed, with emphasis on aspects relevant for X-ray binaries and Cataclysmic Variables. The text gives a general introduction as well as a selective discussion of a number of more recent topics.

A rapid timing analysis of VLT/ULTRACAM and RXTE observations of the black hole binary GX 339-4 in its 2007 low/hard state is presented. The optical light curves in the r', g' and u' filters show slow (~20 s) quasi-periodic variability. Upon this is superposed fast flaring activity on times approaching the best time resolution probed (~50 ms) and with maximum strengths of more than twice the local mean. Power spectral analysis over ~0.004-10 Hz is presented, and shows that although the average optical variability amplitude is lower than that in X-rays, the peak variability power emerges at a higher Fourier frequency in the optical. Energetically, we measure a large optical vs. X-ray flux ratio, higher than that seen when the source was fully jet-dominated. Such a large ratio cannot be easily explained with a disc alone. The optical:X-ray cross-spectrum shows a markedly different behaviour above and below ~0.2 Hz. The peak of the coherence function above this threshold is associated with a short optical time lag, also seen as the dominant feature in the time-domain cross-correlation at ~150 ms. The rms energy spectrum of these fast variations is best described by distinct physical components over the optical and X-ray regimes, and also suggests a maximal disc fraction of 20% at ~5000 A. If the constant time delay is due to propagation of fluctuations to (or within) the jet, this is the clearest optical evidence to date of the location of this component. The low-frequency QPO is seen in the optical but not in X-rays. Evidence of reprocessing emerges at the lowest Fourier frequencies, with optical lags at ~10 s and strong coherence in the blue u' filter. Simultaneous optical spectroscopy also shows the Bowen fluorescence blend, though its emission location is unclear. But canonical disc reprocessing cannot dominate the optical power easily, nor explain the fast variability. (abridged)

Some accreting neutron stars and young stars show unexplained episodic flares in the form of quasi-periodic oscillations or recurrent outbursts. In a series of two papers we present new work on an instability that can lead to episodic outbursts when the accretion disc is truncated by the star's strong magnetic field close to the corotation radius (where the Keplerian frequency matches the star's rotational frequency). In this paper we outline the physics of the instability and use a simple parameterization of the disc-field interaction to explore the instability numerically, which we show can lead to repeated bursts of accretion as well as steady-state solutions, as first suggested by Sunyaev and Shakura. The cycle time of these bursts increases with decreasing accretion rate. These solutions show that the usually assumed `propeller' state, in which mass is ejected from the system, does not need to occur even at very low accretion rates.

Some established views of the solar magnetic cycle are discussed critically, with focus on two aspects at the core of most models: the role of convective turbulence, and the role of the `tachocline' at the base of the convection zone. The standard view which treats the solar cycle as a manifestation of the interaction between convection and magnetic fields is shown to be misplaced. The main ingredient of the solar cycle, apart from differential rotation, is instead buoyant instability of the magnetic field itself. This view of the physics of the solar cycle was already established in the 1950s, but has been eclipsed mathematically by mean field turbulence formalisms which make poor contact with observations and have serious theoretical problems. The history of this development in the literature is discussed critically. The source of the magnetic field of the solar cycle is currently assumed to be located in the `tachocline': the shear zone at the base of the convection zone. While the azimuthal field of the cycle is indeed most likely located at the base of the convection zone, it cannot be powered by the radial shear of the tachocline as assumed in these models, since the radiative interior does not support significant shear stresses. Instead, it must be the powered by the latitudinal gradient in rotation rate in the convection zone, as in early models of the solar cycle. Possible future directions for research are briefly discussed.

The Sun supports a rich spectrum of internal waves that are continuously excited by turbulent convection. The GONG network and the MDI/SOHO space instrument provide an exceptional data base of spatially-resolved observations of solar oscillations, covering an entire sunspot cycle (11 years). Local helioseismology is a set of tools for probing the solar interior in three dimensions using measurements of wave travel times and local mode frequencies. Local helioseismology has discovered (i) near-surface vector flows associated with convection (ii) 250 m/s subsurface horizontal outflows around sunspots (iii) ~50 m/s extended horizontal flows around active regions (converging near the surface and diverging below), (iv) the effect of the Coriolis force on convective flows and active region flows (v) the subsurface signature of the 15 m/s poleward meridional flow, (vi) a +/-5 m/s time-varying depth-dependent component of the meridional circulation around the mean latitude of activity, and (vii) magnetic activity on the far side of the Sun.

The magnetic theory for the production of jets by accreting objects is reviewed with emphasis on outstanding problem areas. An effort is made to show the connections behind the occasionally diverging nomenclature in the literature, to contrast the different points of view about basic mechanisms, and to highlight concepts for interpreting the results of numerical simulations. The role of dissipation of magnetic energy in accelerating the flow is discussed, and its importance for explaining high Lorentz factors. The collimation of jets to the observed narrow angles is discussed, including a critical discussion of the role of `hoop stress'. The transition between disk and outflow is one of the least understood parts of the magnetic theory; its role in setting the mass flux in the wind, in possible modulations of the mass flux, and the uncertainties in treating it realistically are discussed. Current views on most of these problems are still strongly influenced by the restriction to 2 dimensions (axisymmetry) in previous analytical and numerical work; 3-D effects likely to be important are suggested. An interesting problem area is the nature and origin of the strong, preferably highly ordered magnetic fields known to work best for jet production. The observational evidence for such fields and their behavior in numerical simulations is discussed. I argue that the presence or absence of such fields may well be the `second parameter' governing not only the presence of jets but also the X-ray spectra and timing behavior of X-ray binaries.

We have performed 2.5D and 3D simulations of conical jets driven by the rotation of an ordered, large-scale magnetic field in a stratified atmosphere. The simulations cover about three orders of magnitude in distance to capture the centrifugal acceleration as well as the evolution past the Alfven surface. We find that the jets develop kink instabilities, the characteristics of which depend on the velocity profile imposed at the base of the flow. The instabilities are especially pronounced with a rigid rotation profile, which induces a shearless magnetic field. The jet's expansion appears to be limiting the growth of Alfven mode instabilities.

Recent observations of two black hole candidates (GX 339-4 and J1753.5-0127) in the low-hard state (L_X/L_Edd ~ 0.003-0.05) suggest the presence of a cool accretion disk very close to the innermost stable orbit of the black hole. This runs counter to models of the low-hard state in which the cool disk is truncated at a much larger radius. We study the interaction between a moderately truncated disk and a hot inner flow. Ion-bombardment heats the surface of the disk in the overlap region between a two-temperature advection-dominated accretion flow and standard accretion disk, producing a hot (kT_e ~70 keV) layer on the surface of the cool disk. The hard X-ray flux from this layer heats the inner parts of the underlying cool disk, producing a soft X-ray excess. Together with interstellar absorption these effects mimic the thermal spectrum from a disk extending to the last stable orbit. The results show that soft excesses in the low-hard state are a natural feature of truncated disk models.

We present the discovery of optical/X-ray flux correlations on rapid timescales in the low/hard state of the Galactic black hole GX 339-4. The source had recently emerged from outburst and was associated with a relatively-faint counterpart with mag V~17. The optical (VLT/ULTRACAM) and X-ray (RXTE/PCA) data show a clear positive cross-correlation function (CCF) signal, with the optical peak lagging X-rays by ~ 150 ms, preceded by a shallow rise and followed by a steep decline along with broad anti-correlation dips. Examination of the light curves shows that the main CCF features are reproduced in superpositions of flares and dips. The CCF peak is narrow and the X-ray auto-correlation function (ACF) is broader than the optical ACF, arguing against reprocessing as the origin for the rapid optical emission. X-ray flaring is associated with spectral hardening, but no corresponding changes are detected around optical peaks and dips. The variability may be explained in the context of synchrotron emission with interaction between a jet and a corona. The complex CCF structure in GX 339-4 has similarities to that of another remarkable X-ray binary XTE J1118+480, in spite of showing a weaker maximum strength. Such simultaneous multi-wavelength, rapid timing studies provide key constraints for modeling the inner regions of accreting stellar sources.

Possible origins of the magnetic fields of neutron stars include inheritance from the main sequence progenitor and dynamo action at some stage of evolution of progenitor. Inheritance is not sufficient to explain the fields of magnetars. Energetic considerations point to differential rotation in the final stages of core collapse process as the most likely source of field generation, at least for magnetars. A runaway phase of exponential growth is needed to achieve sufficient field amplification during relevant phase of core collapse; it can probably be provided by a some form of magnetorotational instability. Once formed in core collapse, the field is in danger of decaying again by magnetic instabilities. The evolution of a magnetic field in a newly formed neutron star is discussed, with emphasis on the existence of stable equilibrium configurations as end products of this evolution, and the role of magnetic helicity in their existence.

We explore the observational appearance of the photosphere of an ultrarelativistic flow with internal dissipation of energy as predicted by the magnetic reconnection model. Previous study of the radiative transfer in the photospheric region has shown that gradual dissipation of energy results in a hot photosphere. There, inverse Compton scattering of the thermal radiation advected with the flow leads to powerful photospheric emission with spectral properties close to those of the observed prompt GRB emission. Here, we build on that study by calculating the spectra for a large range of the characteristics of the flow. An accurate fitting formula is given that provides the photospheric spectral energy distribution in the ~10 keV to ~10 MeV energy range (in the central engine frame) as a function of the basic physical parameters of the flow. It facilitates the direct comparison of the model predictions with observations, including the variability properties of the lightcurves. We verify that the model naturally accounts for the observed clustering in peak energies of the E*f(E) spectrum. In this model, the Amati relation indicates a tendency for the most luminous bursts to have more energy per baryon. If this tendency also holds for individual GRB pulses, the model predicts the observed narrowing of the width of pulses with increasing photon energy.

We present results of numerical 3D MHD simulations with radiative energy transfer of fine structure in a small sunspot of about 4 Mm width. The simulations show the development of filamentary structures and flow patterns that are, except for the lengths of the filaments, very similar to those observed. The filamentary structures consist of gaps with reduced field strength relative to their surroundings. Calculated synthetic images show dark cores like those seen in the observations; the dark cores are the result of a locally elevated $\tau=1$ surface. The magnetic field in these cores is weaker and more horizontal than for adjacent brighter structures, and the core support a systematic outflow. Movies show migration of the dark-cored structures towards the umbra, and fragments of magnetic flux that are carried away from the spot by a large scale `moat flow'. We conclude that the simulations are in qualitative agreement with observed penumbra filamentary structures, Evershed flows and moving magnetic features.

We present numerical 2D magnetostatic models for sunspot penumbrae consisting of radially aligned field-free gaps in a potential magnetic field, as proposed by Spruit and Scharmer (2006). The shape of the gaps and the field configurations around them are computed consistently from the condition of magnetostatic pressure balance between the gap and the magnetic field. The results show that field-free gaps in the {\it inner} penumbra are cusp-shaped and bounded by a magnetic field inclined by about $70^\circ$ from the vertical. Here, the magnetic component has a Wilson depression on the order 200--300 km relative to the top of the field-free gap; the gaps should thus appear as noticeably elevated features. This structure explains the large variations in field strength in the inner penumbra inferred from magnetograms and two-component inversions, and the varying appearance of the inner penumbra with viewing angle. In the {\it outer} penumbra, on the other hand, the gaps are flat-topped with a horizontal magnetic field above the middle of the gap. The magnetic field has large inclination variations horizontally, but only small fluctuations in field strength, in agreement with observations.

Instabilities in magnetic fields wound up by differential rotation as reviewed in Spruit (1999) are discussed with some detail and new developments added. In stellar models which include magnetic torques, the differential rotation tends to accumulate in the gradients in composition. In view of this, instability in a $\mu$-gradient is studied in more detail here, resulting in the detection of a second instability. Its relevance for angular momentum transport is uncertain, however, since it requires high horizontal field gradients and would not operate near the pole. Finally, the possibility is discussed that magnetic instability in a $\mu$-gradient will lead to {\em layer formation}: the gradients breaking up into small steps of uniform composition and rotation rate. This would enhance the angular momentum transport across inhomogeneous zones, and decrease the rotation rates of the end products of stellar evolution. A recent {\tt astro-ph} submission by Dennisenkov and Pinsonneault proposing a modification of the instability conditions is shown to contain a mathematical error.

The role of kink instability in magnetically driven jets is explored through numerical one-dimensional steady relativistic MHD calculations. The instability is shown to have enough time to grow and influence the dynamics of Poynting-flux dominated jets. In the case of AGN jets, the flow becomes kinetic flux dominated at distances larger than ~1000 r_g because of the rapid dissipation of Poynting flux. When applied to GRB outflows, the model predicts more gradual Poynting dissipation and moderately magnetized flow at distances of ~10^{16} cm where the deceleration of the ejecta due to its interaction with the external medium is expected. The energy released by the instability can power the compact ``blazar zone'' emission and the prompt emission of GRB outflows with high radiative efficiencies.

We use numerical MHD to look at the stability of a possible poloidal field in neutron stars (Flowers & Ruderman 1977), and follow its unstable evolution, which leads to the complete decay of the field. We then model a neutron star after the formation of a solid crust of high conductivity. As the initial magnetic field we use the stable `twisted torus' field which was the result of our earlier work (Braithwaite & Nordlund 2005), since this field is likely to exist in the interior of the star at the time of crust formation. We follow the evolution of the field under the influence of diffusion, and find that large stresses build up in the crust, which will presumably lead to cracking. We put this forward as a model for outbursts in soft gamma repeaters.

We present numerical simulations of a self-sustaining magnetic field in a differentially rotating non-convective stellar interior. A weak initial field is wound up by the differential rotation; the resulting azimuthal field becomes unstable and produces a new meridional field component, which is then wound up anew, thus completing the `dynamo loop'. This effect is observed both with and without a stable stratification. A self-sustained field is actually obtained more easily in the presence of a stable stratification. The results confirm the analytical expectations of the role of Tayler instability.

We interpret penumbral filaments as due to convection in field-free, radially aligned gaps just below the visible surface of the penumbra, intruding into a nearly potential field above. This solves the classical discrepancy between the large heat flux and the low vertical velocities observed in the penumbra. The presence of the gaps causes strong small-scale fluctuations in inclination, azimuth angle and field strength, but without strong forces acting on the gas. The field is nearly horizontal in a region around the cusp-shaped top of the gap, thereby providing an environment for Evershed flows. We identify this region with the recently discovered dark penumbral cores. Its darkness has the same cause as the dark lanes in umbral light-bridges, reproduced in numerical simulations by Nordlund and Stein (2005). We predict that the large vertical and horizontal gradients of the magnetic field inclination and azimuth in the potential field model will produce the net circular polarization seen in observations. The model also explains the significant elevation of bright filaments above their surroundings. It predicts that dark areas in the penumbra are of two different kinds: dark filament cores containing the most inclined (horizontal) fields, and regions between bright filaments, containing the least inclined field lines.

We suggest a possible mechanism of efficient transport of the large-scale external magnetic flux inward through a turbulent accretion disk. The outward drift by turbulent diffusion which limits the effectiveness of this process is reduced if the external flux passes through the disk in concentrated patches. Angular momentum loss by a magnetocentrifugal wind from these patches of strong field adds to the inward drift. We propose that magnetic flux accumulating in this way at the center of the disk provides the `second parameter' determining X-ray spectral states of black hole binaries, and the presence of relativistic outflows.

Some main-sequence stars of spectral type A are observed to have a strong (300-30000 G), static, large-scale magnetic field, of a chiefly dipolar shape -- the `Ap stars' (for example Alioth, the fifth star in the Big Dipper). Following their discovery fifty years ago, it was speculated that these fields a remnants of the star's formation, a `fossil' field. Alternatively, they could be generated by a dynamo process in the star's convective core. The dynamo hypothesis has difficulty explaining high field strengths and the lack of a correlation with rotation. The weakness of the fossil-field theory has been the absence of known field configurations stable enough to survive in a star over its lifetime. We demonstrate here the formation of stable magnetic field configurations, with properties agreeing with those observed, by evolution from arbitrary, unstable initial fields. The results are applicable equally to Ap stars, magnetic white dwarfs and some highly magnetised neutron stars known as magnetars. This establishes fossil fields as the natural, unifying explanation for the magnetism of all these stars.

We calculate the evaporation of a cool accretion disk around a black hole due to the ion-bombardment by an ion supported accretion flow (here ISAF, or optically thin ADAF). As first suggested by Spruit & Deufel (2002), this evaporation takes place in two stages: ion bombardment of the cool disk (Shakura-Sunyaev disk: SSD) produces an intermediate-temperature layer on top of the disk (`warm layer') which constitutes an independent accretion flow on both sides of the SSD. As this warm material accretes inward of the inner radius of the SSD, it becomes thermally unstable by lack of cooling of photons, and evaporates into the ISAF, thereby feeding the latter. Angular momentum conservation forces a certain fraction of the ISAF material to move outward, where it can bombard the SSD with its hot ions. The flow geometry is derived by computing stationary solutions of the continuity- and angular momentum equations for the three components (ISAF, warm flow and SSD). The overall radiative output is dominated by hard X-rays. They are produced mostly from the warm component, rather than the ISAF. The expected time dependence and stability of the flow, not computed here, is discussed briefly.

We investigate the production of the gamma-ray spectrum of a Poynting-flux dominated GRB ouflow. The very high magnetic field strengths (super-equipartition) in such a flow lead to very efficient synchrotron emission. In contrast with internal shocks, dissipation of magnetic energy by reconnection is gradual and does not produce the spectrum of cooling electrons associated with shock acceleration.We show that a spectrum with a break in the BATSE energy range is produced, instead, if the magnetic dissipation heats a small (\sim 10^{-4}) population of electrons.

As a massive star evolves through multiple stages of nuclear burning on its way to becoming a supernova, a complex, differentially rotating structure is set up. Angular momentum is transported by a variety of classic instabilities, and also by magnetic torques from fields generated by the differential rotation. We present the first stellar evolution calculations to follow the evolution of rotating massive stars including, at least approximately, all these effects, magnetic and non-magnetic, from the zero-age main sequence until the onset of iron-core collapse. The evolution and action of the magnetic fields is as described by Spruit 2002 and a range of uncertain parameters is explored. In general, we find that magnetic torques decrease the final rotation rate of the collapsing iron core by about a factor of 30 to 50 when compared with the non-magnetic counterparts. Angular momentum in that part of the presupernova star destined to become a neutron star is an increasing function of main sequence mass. That is, pulsars derived from more massive stars will rotate faster and rotation will play a more dominant role in the star's explosion. The final angular momentum of the core is determined - to within a factor of two - by the time the star ignites carbon burning. For the lighter stars studied, around 15 solar masses, we predict pulsar periods at birth near 15 ms, though a factor of two range is easily tolerated by the uncertainties. Several mechanisms for additional braking in a young neutron star, especially by fall back, are also explored.

We study possible mechanisms of excitation of quasiperiodic oscillations in the accretion flow of black hole accreters in their hard spectral states, in the context of the `truncated disk' model. Quasi-spherical oscillations of the inner ion-supported accretion flow (ISAF) can be excited by the interaction of this hot flow with the cool disk extending outward from it. The fundamental mode of (p-mode) oscillation is most easily excited, and has a frequency near the Kepler frequency at the inner edge of the cool disk. The strongest excitation mechanism is a feedback loop involving cooling of the ISAF by soft photons from the cool disk and heating of the cool disk by the ISAF, while synchrotron emission can be a relatively strong damping effect. Growth times are computed by detailed Comptonization calculations of the interaction of an idealized ISAF with a cool disk. Typical growth times as short as a few dynamical times are found, while amplitudes can be as large as 10%.

Coordinated observations X-ray+optical observations of two southern X-ray binaries, the black hole candidate XTE J1746-321 and the neutron star accreter Cir X-1 (a `microquasar') are reported. With a photon counting optical photometer on the 1.9m telescope at Sutherland, South Africa and the PCA detector on RXTE, 4h each of simultaneous data were obtained on XTE J1746 and Cir X-1. Cir X-1 showed no X-ray variability at the 2% level, XTE J1746 was variable at 5-7% with a 5Hz QPO. Cross-correlation yielded no correlated signals on either source, to a level of 1%. A problem with a recently published orbital ephemeris of Cir X-1 is pointed out.

We study the local stability of stratified, differentially-rotating fluids to axisymmetric perturbations in the presence of a weak magnetic field and of finite resistivity, viscosity and heat conductivity. This is a generalization of the Goldreich-Schubert-Fricke (GSF) double-diffusive analysis to the magnetized and resistive, triple-diffusive case. Our fifth-order dispersion relation admits a novel branch which describes a magnetized version of multi-diffusive modes. We derive necessary conditions for axisymmetric stability in the inviscid and perfect-conductor (double-diffusive) limits. In each case, rotation must be constant on cylinders and angular velocity must not decrease with distance from the rotation axis for stability, irrespective of the relative strength of viscous, resistive and heat diffusion. Therefore, in both double-diffusive limits, solid body rotation marginally satisfies our stability criteria. The role of weak magnetic fields is essential to reach these conclusions. The triple-diffusive situation is more complex, and its stability criteria are not easily stated. Numerical analysis of our general dispersion relation confirms our analytic double-diffusive criteria, but also shows that an unstable double-diffusive situation can be significantly stabilized by the addition of a third, ostensibly weaker, diffusion process. We describe a numerical application to the Sun's upper radiative zone and establish that it would be subject to unstable multi-diffusive modes if moderate or strong radial gradients of angular velocity were present.

We present a search for excess mid-IR emission due to circumbinary material in the orbital plane of cataclysmic variables (CVs). Our motivation stems from the fact that the strong braking exerted by a circumbinary (CB) disc on the binary system could explain several puzzles in our current understanding of CV evolution. Since theoretical estimates predict that the emission from a CB disc can dominate the spectral energy distribution (SED) of the system at wavelengths > 5 microns, we obtained simultaneous visible to mid-IR SEDs for eight systems. We report detections of SS Cyg at 11.7 microns and AE Aqr at 17.6 microns, both in excess of the contribution from the secondary star. In AE Aqr, the IR likely originates from synchrotron-emitting clouds propelled by the white dwarf. In SS Cyg, we argue that the observed mid-IR variability is difficult to reconcile with simple models of CB discs and we consider free-free emission from a wind. In the other systems, our mid-IR upper limits place strong constraints on the maximum temperature of a putative CB disc. The results show that if any sizeable CB disc are present in these systems, they must be self-shadowed or perhaps dust-free, with the peak thermal emission shifted to far-IR wavelengths.

The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron stars and their spin behavior during outbursts. In particular, it is unclear from the standard paradigm how these systems manage to accrete over such a wide range in dM/dt (i.e., >~ a factor of 50), and why the neutron stars exhibit a high rate of spindown in at least a number of cases. Following up on prior suggestions, we propose that `fast' X-ray pulsars can continue to accrete, and that their accretion disks terminate at approximately the corotation radius. We demonstrate the existence of such disk solutions by modifying the Shakura-Sunyaev equations with a simple magnetic torque prescription. The solutions are completely analytic, and have the same dependence on dM/dt and alpha (the viscosity parameter) as the original Shakura-Sunyaev solutions; but, the radial profiles can be considerably modified, depending on the degree of fastness. We apply these results to compute the torques expected during the outbursts of the transient millisecond pulsars, and find that we can explain the large spindown rates that are observed for quite plausible surface magnetic fields of ~10^9 G.

We use both time-domain and Fourier techniques to study the correlated optical and X-rays variability in the black hole X-ray nova XTE J1118+480. Its X-ray timing properties, such as the shape of the X-ray power spectrum and cross-correlation functions (CCFs) between X-ray energy bands, the slight decrease of RMS variability from 30 % in the 2-5.9 keV to 19 % in the 15.5--44.4 keV band, as well as the X-ray hardness/flux anti-correlation, are very similar to what is found in other black hole binaries in the hard state. The optical/X-ray CCF is virtually independent of the X-ray energies. The optical flux appears to be correlated not only with the X-ray flux but also with the X-ray spectral variability. Both the coherence function and the lags between optical and the X-rays are Fourier frequency dependent. The optical/X-ray coherence function reaches its maximum (~0.3) in the 0.1-1 Hz range and the time-lags decrease with frequency approximatively like f^-0.8. The correlation between X-ray and optical light curves appears to have time-scale-invariant properties. The X-ray/optical CCF maintains a similar but rescaled shape on time-scales ranging at least from 0.1 s to few 10 s. Using the event superposition method we show that the correlation is not triggered by a single type of event (dip or flare) in the light curves. Instead, optical and X-ray fluctuations of very different shapes, amplitudes and time-scales are correlated in a similar mode where the optical light curve is seemingly related to the time derivative of the X-rays.

The physics of GRB powered by a magnetic energy flux is reviewed. Magnetic fields are natural for transmitting the energy from the central compact object to the small amount of baryons required for a GRB. When dissipation of the flux of magnetic energy by reconnection inside the flow is taken into account, the magnetic model assumes several more convincing properties. For baryon loading typical of observed GRB, most of the dissipation takes place just outside photosphere, so that prompt emission is produced efficiently, and the magnetic field strength in this region is high, resulting in efficient synchrotron emission. Remarkably, the dissipation also causes very efficient acceleration of the bulk flow. This effect is illustrated with a classical hydrodynamic equivalent. In this context, the distinction between the flux of magnetic energy $cB^2/8\pi$ and the Poynting flux $cB^2/4\pi$ is important, and an interpretation of the Poynting flux as a `magnetic enthalpy flux' illuminating. Numerical and analytical results for flow acceleration and the relative contribution of photospheric (thermal) and nonthermal emission as functions of the asymptotic bulk Lorentz factor are given. The transition between X-ray flashes and true GRB is predicted at $\Gamma\approx 100$.

Rotation in massive stars has been studied on the main sequence and during helium burning for decades, but only recently have realistic numerical simulations followed the transport of angular momentum that occurs during more advanced stages of evolution. The results affect such interesting issues as whether rotation is important to the explosion mechanism, whether supernovae are strong sources of gravitational radiation, the star's nucleosynthesis, and the initial rotation rate of neutron stars and black holes. We find that when only hydrodynamic instabilities (shear, Eddington-Sweet, etc.) are included in the calculation, one obtains neutron stars spinning at close to critical rotation at their surface -- or even formally in excess of critical. When recent estimates of magnetic torques (Spruit 2002) are added, however, the evolved cores spin about an order of magnitude slower. This is still more angular momentum than observed in young pulsars, but too slow for the collapsar model for gamma-ray bursts.

A model is presented that explains the `torsional oscillation' pattern of deviations in the solar rotation rate as a geostrophic flow. The flow is driven by temperature variations near the surface due to the enhanced emission of radiation by the small scale magnetic field. The model explains the sign of the flow, its amplitude and the fact that the maxima occur near the boundaries of the main activity belts. The amplitude of the flow decreases with depth from its maximum at the surface but penetrates over much of the depth of the convection zone, in agreement with the data from helioseismology. It predicts that the flow is axisymmetric only on average, and in reality consists of a superposition of circulations around areas of enhanced magnetic activity. It must be accompanied by a meridional flow component, which declines more rapidly with depth.

We investigate the effects of magnetic energy release by local magnetic dissipation processes in Poynting flux-powered GRBs. For typical GRB parameters (energy and baryon loading) the dissipation takes place mainly outside the photosphere, producing non-thermal radiation. This process converts the total burst energy into prompt radiation at an efficiency of 10-50%. At the same time the dissipation has the effect of accelerating the flow to a large Lorentz factor. For higher baryon loading, the dissipation takes place mostly inside the photosphere, the efficiency of conversion of magnetic energy into radiation is lower, and an X-ray flash results instead of a GRB. We demonstrate these effects with numerical one-dimensional steady relativistic MHD calculations.

We present coordinated X-ray and optical observations of the X-ray transient KV UMa (= XTE J1118+480) at high time resolution. The optical variation associated with spikes in the X-rays consist of a dip preceding the spike by 0-5 s, and a peak following it. The cross-correlation function between X-rays and optical has the same shape. It rises sharply, within 30ms of the peak of the X-rays. The shape of the cross correlation shows significant intrinsic variability on time scales as short as 30s. Analyzing this variability with principal component analysis (calibrated with simulated data), we find two statistically independent components. The first is similar to the average cross-correlation. The second also consists of a dip followed by a peak, but on a 3x shorter time scale. These properties of the optical light, together with the high optical brightness of the transient, are not easily explained by reprocessing of the X-rays.

The X-ray spectrum from a cool accretion disk heated by virialized protons is computed. The cool disk is either embedded in a magnetically heated accretion disk corona or partly extends into an ion supported torus (or ADAF). We calculate the stationary equilibrium between proton heating, electron thermal conduction and the radiative losses by bremsstrahlung and Compton scattering. A heated surface layer on top of the accretion disk is produced with temperatures between 60--90 keV above a cool layer with temperatures of 0.01 keV (AGN) and 1keV (galactic black hole candidates). The spectra produced by the surface layer are reminiscent of hard state spectra, but a bit too steep, especially for AGN's. Near the inner edge of the disk, where the optical depth of the disk $\tau lesssim 1$, we find that the cool component of the disk disappears. Instead, the hot protons from the corona/ADAF heat the disk, on a dynamical time-scale, to temperatures of several 100 keV, limited by pair production. This region, here called a `warm disk', could contribute significantly to the hard X-ray spectra and could be important for feeding material into an ADAF.

We show that the inner regions of a cool accretion disk in an X-ray binary can transform into an advective, ion supported accretion flow (an optically thin ADAF, here called ISAF), through events involving only the known properties of the Coulomb interaction in a two-temperature plasma, standard radiation processes, and viscous heating. The optically thin inner edge of the disk is heated to a few 100 keV by the strong flux of hot ions from the surrounding hot ISAF. We show that he resident ions in this `warm' disk are thermally unstable due to internal viscous heating, and heat up to their virial temperature. The innermost disk regions thus evaporate and feed the ISAF. These processes are demonstrated with time dependent calculations of a two-temperature plasma in vertical hydrostatic equilibrium, including heating by external ions, internal proton--electron energy exchange, and viscous heating. The process complements the `coronal' evaporation mechanism which operates at larger distances from the central object.

We investigate the structure and evolution of a geometrically thin viscous Keplerian circumbinary (CB) disk, using detailed models of their radiative/convective vertical structure. We use a simplified description for the evolution of the cataclysmic binary and focus on cases where the circumbinary disk causes accelerated mass transfer (> 1e-8 Msun/yr). The inner edge of the disk is assumed to be determined by the tidal truncation radius and the mass input rate into the disk is assumed to be a small fraction (1e-5-0.01) of the mass transfer rate. Under the action of the viscous stresses in the disk the matter drifts outward with the optically thick region extending to several AU. The inner part of the disk is cool with maximum effective temperatures < 3,000 K while the outermost parts of the disk are < 30 K and optically thin. We calculate the effects of thermal instability on a sufficiently massive CB disk. It leads to outbursts reminiscent of those in thermally unstable accretion disks, with the instability remaining confined to the inner regions of the CB disk. However, for most of the evolutionary sequences the surface densities required to trigger instability are not reached. The spectral energy distributions from circumbinary disks are calculated, and the prospects for the detection of such disks in the infrared and submm wavelength regions are discussed.

Observations of Cataclysmic Variables show a number of phenomena that do not fit easily into the standard magnetic braking scenario. These include the large spread in mass transfer rates, the low surface temperatures of many of the companion stars, and evidence for material at low velocities. We propose that these anomalies have a common cause: the presence of a circumbinary (CB) disk. This may be a significant component of mass transferring binaries in general. Direct detection of such CB disks may be possible but difficult, because of their low luminosity Disks surrounding Cataclysmic Binaries and spectral energy distribution peaking in the mid-IR.

Magnetic fields can be created in stably stratified (non-convective) layers in a differentially rotating star. A magnetic instability in the toroidal field (wound up by differential rotation) replaces the role of convection in closing the field amplification loop. Tayler instability is likely to be the most relevant magnetic instability. A dynamo model is developed from these ingredients, and applied to the problem of angular momentum transport in stellar interiors. It produces a prodominantly horizontal field. This dynamo process might account for the observed pattern of rotation in the solar core.

We report on fast (msec) simultaneous observations in X-rays and optical light of the black hole candidate XTE J1118+480, showing correlation on time scales from 20msec to several seconds. The shape of the correlation is not easily explained with existing models. We suggest that is is caused by optical cyclosynchrotron emission in a dense, relatively slow outflow from the inner regions of the accretion.

We present detailed simulations of the interaction of a stream emanating from a mass-losing secondary with the core of a massive supergiant in the slow merger of the two stars inside a common envelope. The dynamics of the stream can be divided into a ballistic phase, starting at the L_1 point, and a hydrodynamical phase where the stream interacts strongly with the core. Considering the merger of a 1 and 5Msun star with a 20Msun evolved supergiant, we present two-dimensional hydrodynamical simulations using the PROMETHEUS code to demonstrate how the penetration depth and post-impact conditions depend on the initial properties of stream material (e.g. entropy, angular momentum, stream width) and the properties of the core (e.g. density structure and rotation rate). Using these results, we present a fitting formula for the entropy generated in the stream--core interaction and a recipe for the determination of the penetration depth based on a modified Bernoulli integral.

We investigate an interesting new class of high-mass X-ray binaries (HMXBs) with long orbital periods (P_orb > 30 days) and low eccentricities (e <~ 0.2). The orbital parameters suggest that the neutron stars in these systems did not receive a large impulse, or ``kick,'' at the time of formation. We develop a self-consistent phenomenological picture wherein the neutron stars born in the observed wide HMXBs receive only a small kick (<~ 50 km/s), while neutron stars born in isolation, in the majority of low-mass X-ray binaries, or in many of the well-known HMXBs with P_orb <~ 30 days receive the conventional large kicks, with a mean speed of ~ 300 km/s. We propose that the magnitude of the natal kick to a neutron star born in a binary system depends on the rotation rate of the pre-collapse core. We further suggest that the rotation rate of the core is a strong, well-defined function of the evolutionary path of the progenitor star.

We present a linear analysis of the stability of accretion disks in which angular momentum is removed by the magnetic torque exerted by a centrifugally driven wind. The effects of the dependence of the wind torque on field strength and inclination, the sub-Keplerian rotation due to magnetic forces, and the compression of the disk by the field are included. A WKB dispersion relation is derived for the stability problem. We find that the disk is always unstable if the wind torque is strong. At lower wind torques instability also occurs provided the rotation is close to Keplerian. The growth time scale of the instability can be as short as the orbital time scale. The instability is mainly the result of the sensitivity of the mass flux to changes in the inclination of the field at the disk surface. Magnetic diffusion in the disk stabilizes if the wind torque is small.

We consider the interaction of a slowly rotating unmagnetized neutron star with a hot (ion supported, ADAF) accretion flow. The virialized protons of the ADAF penetrate into the neutron star atmosphere, heating a surface layer. Detailed calculations are presented of the equilibrium between heating by the protons, electron thermal conduction, bremsstrahlung and multiple Compton scattering in this layer. Its temperature is of the order 40-70 keV. Its optical depth increases with the incident proton energy flux, and is of the order unity for accretion at $10^{-2}$--$10^{-1}$ of the Eddington rate. At these rates, the X-ray spectrum produced by the layer has a hard tail extending to 100 keV, and is similar to the observed spectra of accreting neutron stars in their hard states. The steep gradient at the base of the heated layer gives rise to an excess of photons at the soft end of the spectrum (compared to a blackbody) through an `inverse photosphere effect'. The differences with respect to previous studies of similar problems are discussed, they are due mostly to a more accurate treatment of the proton penetration process and the vertical structure of the heated layer.

We have performed Monte Carlo simulations of Compton upscattering of low-energy photons in an accretion disk around a Schwarzschild black hole. The photons gain energy from the rotational motion of the electrons in the disk. The upscattering occurs near the black hole horizon, where the flow velocity of the electrons approaches the speed of light. We show that this type of bulk-flow Comptonization can produce power-law X-ray spectra similar to the ones observed in black-hole X-ray transients in the high/soft state, i.e., a soft bump dominating the spectrum below ~ 10 keV and a power-law tail with photon index in the range 2-3. In order to reproduce the observed hard to soft flux ratio the disk has to have vertical optical depth above ~ 3 at the last stable orbit. We conclude that the power-law component of the high/soft state of black-hole transients may be due to an intrinsically cool disk extending all the way to the hole, without a separate hot plasma component.

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The 3--dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou (1995). In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outward and allows mass to accrete inward. The instability takes the form of a $m=1$ rotating `crescent', reminiscent of the purely hydrodynamic nonlinear instability previously found in pressure-supported discs. A model where the initial surface mass density $\Sigma(r)$ and $B_{\mathrm{z}}(r)$ decrease with radius as power laws shows transient mass accretion during about 6 orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.

We study the merging of massive stars inside a common envelope for binary systems consisting of a red supergiant with a mass of 15-20 Msun and a main-sequence companion of 1-5 Msun. We are particularly interested in the stage when the secondary, having overfilled its Roche lobe inside the common envelope, starts to transfer mass to the core of the primary at a very high mass-transfer rate and the subsequent nucleo-synthesis in the core-impact region. Using a parametrized model for the structure of the envelope at this stage, we perform 2-dimensional hydrodynamical calculations with the Munich Prometheus code to calculate the dynamics of the stream emanating from the secondary and its impact on the core of the primary. We find that, for the lower end of the estimated mass-transfer rate, low-entropy, hydrogen-rich material can penetrate deep into the primary core where nucleosynthesis through the hot CNO cycle can take place and that the associated neutron exposure may be sufficiently high for significant s-processing. For mass-transfer rates at the high end of our estimated range and higher densities in the stream, the stream impact can lead to the dredge-up of helium, but the neutron production is too low for significant s-processing.

The influence of a circumbinary (CB) disk on the evolution of cataclysmic variable (CV) binary systems is investigated. We show that CB mass surface densities sufficient to influence the evolution rate are plausibly provided by the outflows observed in CVs, if the net effect of these winds is to deliver $10^{-4}$--$10^{-3}$ of the mass transfer rate to the CB disk. The torque exerted by the CB disk provides a positive feedback between mass transfer rate and CB disk mass which can lead to mass transfer rates of $\sim 10^{-8} -10^{-7} \mpy$. This mechanism may be responsible for causing the range of variation of mass transfer rates in CV's. In particular, it may explain rates inferred for the novalike variables and the supersoft X-ray binary systems observed near the upper edge of the period gap ($P \sim 3 - 4$ hr), as well as the spread in mass transfer rates above and below the period gap. Consquences and the possible observability of such disks are discussed.

We consider possible geometries of magnetic fields in GRB outflows, and their evolution with distance from the source. For magnetically driven outflows, with an assumed ratio of magnetic to kinetic energy density of order unity, the field strengths are sufficient for efficient production of gamma-rays by synchrotron emission in the standard internal shock scenario, without the need for local generation of small scale fields. In these conditions, the MHD approximation is valid to large distances (>10^19cm). In outflows driven by nonaxisymmetric magnetic fields, changes of direction of the field cause dissipation of magnetic energy by reconnection. This dissipation takes place outside the photosphere of the outflow, and can convert a significant fraction of the magnetic energy flux into radiation.

Inhomogeneities in wave propagation conditions near and below the solar surface have been detected by means of time-distance helioseismology. Here we calculate the effect of temperature inhomogeneities on the travel times of sound waves. A temperature increase, e.g. in active regions, not only increases the sound speed but also lengthens the path along which the wave travels because the expansion of the heated layers shifts the upper turning of the waves upward. Using a ray tracing approximation we find that in many cases the net effect of a temperature enhancement is an {\it increase} of the travel times. We argue that the reduced travel times that are observed are caused by a combination of magnetic fields in the active region and {\it reduced} subsurface temperatures. Such a reduction may be related to the increased radiative energy loss from small magnetic flux tubes.

Radiatively inefficient (or advection dominated) disks are discussed at an introductory level. Ion supported and radiation supported flows are discussed, the different consequences of advection dominated flows onto black holes vs. solid surfaces (neutron stars, white dwarfs), hydrodynamics, the role of the ratio of specific heats, and the possible connection between ADAFs and outflows.

In this lecture the basic theory of accretion disks is introduced, with emphasis on aspects relevant for X-ray binaries and Cataclysmic Variables.

Observations at high spectral and temporal resolution are presented of the dwarf nova EX Dra in outburst. The disk seen in the He I line reconstructed by Doppler tomography shows a clear two-armed spiral pattern pointing to spiral shocks in the disk. The Balmer and He II maps also give evidence for the presence of spirals. The eclipse as seen in the red continuum indicates a disk radius of 0.31 times the orbital separation, which might be large enough to explain the observed spiral shocks through excitation by the tidal field of the secondary. The eclipse in the Balmer line profiles, well resolved in our observations, indicates a somewhat smaller disk size (0.25). We discuss the possibility that this is related to an optical depth effect in the lines.

Some topics in the theory of jets are reviewed. These include jet precession, unconfined jets, the origin of knots, the internal shock model as a unifying theme from protostellar jets to Gamma-ray bursts, relations between the Blandford-Znajek and MHD disk-wind models, and jet collimation in magnetic acceleration models.

The theory of stellar structure can be used to identify the most plausible mechanisms for the irradiance variations associated with the solar cycle. Changes in surface emissivity, i.e. the reduced cooling in spots and enhanced emission by small scale magnetic fields, are the most effective mechanisms and account for most of the observed variation. Helioseismology will soon be able to test the consequences of changes in surface emissivity, and distinguish them from other scenarios for irradiance variability.

We generalize a formula for the Goldreich-Julian charge density (\rho), originally derived for rotating neutron star, for arbitrary oscillations of a neutron star with arbitrary magnetic field configuration under assumption of low current density in the inner parts of the magnetosphere. As an application we consider toroidal oscillation of a neutron star with dipole magnetic field and calculate energy losses. For some oscillation modes the longitudinal electric field can not be canceled by putting charged particles in the magnetosphere without a presence of strong electric current j = (c/(\omega r))\rho c. It is shown that the energy losses are strongly affected by plasma in the magnetosphere, and cannot be described by vacuum formulas.

We present a model for the radio outbursts of microquasars based on the assumption of quasi-continuous jet ejection. The jets are `lit up' by shock fronts traveling along the jets during outbursts. The observed comparatively flat decay light curves combined with gradually steepening spectral slopes are explained by a superposition of the radiation of the aging relativistic particle population left behind by the shocks. This scenario is the low energy, time-resolved equivalent to the internal shock model for GRBs. We show that this model predicts energy contents of the radiating plasma similar to the plasmon model. At the same time, the jet model relaxes the severe requirements on the central source in terms of the rate at which this energy must be supplied to the jet. Observations of `mini-bursts' with flat spectral slopes and of infrared emission far from the source centre suggest two different states of jet ejections: (i) A `mini-burst' mode with relatively stable jet production and weak radio emission with flat spectra and (ii) an outburst mode with strong variations in the jet bulk velocities coupled with strong radio emission with steeper spectra. We also show that the continuous jets in microquasars should terminate in strong shocks and possibly inflate radio lobes similar to extragalactic jet sources. Finally, we investigate the possibility of testing the predictions of this model with resolved radio observations.

The X-ray spectrum from a cool disk embedded in an ion supported torus is computed. The interaction of the hot ions with the disk increases the hard X-ray luminosity of the system}. A surface layer of the disk is heated by the protons from the torus. The Comptonized spectrum produced by this layer has a shape that depends only weakly on the incident energy flux and the distance from the accreting compact object. It consists of a `blue bump' of unComptonized soft photons and a flat high energy tail, reminiscent of the observed spectra. The hard tail becomes flatter as the thermalization depth in the cool disk is increased. Further evidence for ion illumination are the Li abundance in the secondaries of low mass X-ray binaries and the 450 keV lines sometimes seen in black-hole transient spectra.

The processes contributing to the evolution of an initially weak magnetic field in a differentially rotating star are reviewed. These include rotational smoothing (akin to convective expulsion) and a list of about 5 instabilities, among them magnetorotational instability, byoyancy instability, and pinch-type instabilities. The important effects of thermal and magnetic diffusion on these instabilities are analyzed in some detail. The first instability to set in is a pinch-type instability. It becomes important in modifying the field configuration before magnetic buoyancy-driven instabilities set in. The evolution of an initially strong field remains a more open question, including the old problem whether dynamically stable magnetic equilibria exist in stars.

A weakly magnetized ($\sim 10^7$G) neutron star, slowly spun up by accretion in an X-ray binary, crosses the instability boundary for r-mode instability at P=1--2 msec. The amplitude of the oscillation, which initially increases only at the spinup time scale, is secularly unstable due to the negative temperature dependence of the viscosity in neutron star matter, and diverges after a few hundred years. Angular momentum loss by the gravitational wave causes strong differential rotation, in which the magnetic field is wound up to $10^{17}$G on a time scale of a few months. When this field becomes unstable to buoyancy instability, a surface field strength of a few $10^{16}$G is produced on a time scale of seconds, which then powers a GRB with energies of $\sim 1-8\times 10^{51}$ and duration of 1--100 sec.

A numerical code is described for constructing Doppler maps from the orbital variation of line profiles of (mass transfering) binaries. It uses an algorithm related to Richardson-Lucy iteration, and is much faster than the standard algorithm used for ME problems. The method has been tested on data of cataclysmic variables (including WZ Sge and SS Cyg), producing maps of up to 300X300 points. It includes an IDL-based set of routines for manipulating and plotting the input and output data, and can be downloaded from http://www.mpa-garching.mpg.de/~henk

RADIO pulsars are thought to born with spin periods of 0.02-0.5 s and space velocities of 100-1000 km/s, and they are inferred to have initial dipole magnetic fields of 10^{11}-10^{13}. The average space velocity of a normal star in the Milky Way is only 30 km/s, which means that pulsars must receive a substantial 'kick' at birth. Here we propose that the birth characteristics of pulsars have a simple physical connection with each other. Magnetic fields maintained by differential rotation between the core and envelope of the progenitor would keep the whole star in a state of approximately uniform rotation until 10 years before the explosion. Such a slowly rotating core has 1000 times less angular momentum than required to explain the rotation of pulsars. Although the specific physical process that 'kicks' the neutron star at birth has not been identified, unless its force is exerted exactly head-on, it will also cause the neutron star to rotate. We identify this process as the origin of the spin of pulsars. Such kicks will cause a correlation between the velocity and spin vectors of pulsars. We predict that many neutron stars are born with periods longer than 2 s, and never become radio pulsars.

I argue that the rotation of white dwarfs is not a remnant of the angular momentum of their main sequence progenitors but a result of the mass loss process on the AGB. Weak magnetic fields, if present in stellar interiors, are likely to maintain approximately uniform rotation in stars, both on the main sequence and on the giant branches. The nearly uniform rotation of the core of the Sun is evidence for the existence of such fields. Exactly axisymmetric mass loss on the AGB from uniformly rotating stars would lead lead to white dwarfs with very long rotation periods ($>$ 10 yr). Small random non-axisymmetries ($\sim 10^{-3}$) in the mass loss process, on the other hand, add sufficient angular momentum to explain the observed rotation periods around one day. The process illustrated with a computation of the probability distribution of the rotation periods under the combined influence of random forcing by weak nonaxisymmetries and angular momentum loss in the AGB superwind. Such asymmetries can in principle be observed by proper motion studies of the clumps in interferometric images of SiO maser emission.

The various effects leading to diversity in the bolometric light curves of supernovae are examined: nucleosynthesis, kinematic differences, ejected mass, degree of mixing, and configuration and intensity of the magnetic field are discussed. In Type Ia supernovae, a departure in the bolometric light curve from the full-trapping decline of $^{56}$Co can occur within the two and a half years after the explosion, depending on the evolutionary path followed by the WD during the accretion phase. If convection has developed in the WD core during the presupernova evolution, starting several thousand years before the explosion, a tangled magnetic field close to the equipartition value should have grown in the WD. Such an intense magnetic field would confine positrons where they originate from the $^{56}$Co decays, and preclude a strong departure from the full-trapping decline, as the supernova expands. This situation is expected to occur in C+O Chandrasekhar WDs as opposed to edge-lit detonated sub-Chandrasekhar WDs. If the pre-explosion magnetic field of the WD is less intense than 10$^{5-8}$G, a lack of confinement of the positrons emitted in the $^{56}$Co decay and a departure from full-trapping decline would occur. The time at which it takes place can provide estimates of the original magnetic field of the WD, its configuration, and also of the mass of the supernova ejecta. In SN 1991bg, the bolometric light curve suggests absence of a significant tangled magnetic field (intensity lower than $10^{3}$ G). Chandrasekhar-mass models do not reproduce the bolometric light curve of this supernova. For SN 1972E, on the contrary, there is evidence for a tangled configuration of the magnetic field and its light curve is well reproduced by a Chandrasekhar WD explosion.

Progress in the theory of stellar convection over the past decade is reviewed. The similarities and differences between convection in stellar envelopes and laboratory convection at high Rayleigh numbers are discussed. Direct numerical simulation of the solar surface layers, with no other input than atomic physics, the equations of hydrodynamics and radiative transfer is now capable of reproducing the observed heat flux, convection velocities, granulation patterns and line profiles with remarkably accuracy. These results show that convection in stellar envelopes is an essentially non-local process, being driven by cooling at the surface. This differs distinctly from the traditional view of stellar convection in terms of local concepts such as cascades of eddies in a mean superadiabatic gradient. The consequences this has for our physical picture of processes in the convective envelope are illustrated with the problems of sunspot heat flux blocking, the eruption of magnetic flux from the base of the convection zone, and the Lithium depletion problem.

The theory of magnetically accelerated outflows and jets from accretion disks is reviewed at an introductory level, with special attention to problem areas like the launching conditions of the flow at the disk surface, stability of the magnetic field, and collimation mechanisms. This text will appear in R.A.M.J. Wijers, M.B. Davies and C.A. Tout, eds., Physical processes in Binary Stars, Kluwer Dordrecht, 1996 (NATO ASI series).

We investigate the stability to nonaxisymmetric perturbations of an accretion disc in which a poloidal magnetic field provides part of the radial support against gravity. Interchange instability due to radial gradients in the magnetic field are strongly stabilized by the shear flow in the disc. For smooth field distributions this instability is restricted to discs in which the magnetic energy is comparable to the gravitational energy. An incompressible model for the instability akin to the Boussinesq approximation for convection is given which predicts the behaviour of the instability accurately. Global axisymmetric disturbances are also considered and found to be stable for a certain class of models. The results indicate that accretion discs may be able to support poloidal fields which are strong enough to suppress other forms of magnetic instability. These strong and stable field distributions are likely to be well suited for the magnetic acceleration of jets and winds.

This is an introduction to accretion disk theory, with emphasis on aspects relevant for X-ray Binaries and Cataclysmic Variables. The text corrects some mistakes in an earlier version, which appeared in 'Lives of Neutron Stars', A. Alpar, \"U. Kizilo\u glu and J. van Paradijs (eds.), Kluwer, Dordrecht (NATO ASI series, 1994).