Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.

Lithium in M 67: From the main sequence to the red giant branch

Context. Lithium abundances in open clusters are a very effective probe of mixing processes, and their study can help us to understand the large depletion of lithium that occurs in the Sun. Owing to its age and metallicity, the open cluster M 67 is especially interesting on this respect. Many studies of lithium abundances in M 67 have been performed, but a homogeneous global analysis of lithium in stars from subsolar masses and extending to the most massive members, has yet to be accomplished for a large sample based on high-quality spectra. Aims. We test our non-standard models, which were calibrated using the Sun with observational data. Methods. We collect literature data to analyze, for the first time in a homogeneous way, the non-local thermal equilibrium lithium abundances of all observed single stars in M 67 more massive than similar to 0.9 M-circle dot. Our grid of evolutionary models is computed assuming a non-standard mixing at metallicity [Fe/H] = 0.01, using the Toulouse-Geneva evolution code. Our analysis starts from the entrance into the zero-age main-sequence. Results. Lithium in M 67 is a tight function of mass for stars more massive than the Sun, apart from a few outliers. A plateau in lithium abundances is observed for turn-off stars. Both less massive (M >= 1.10 M-circle dot) and more massive (M >= 1.28 M-circle dot) stars are more depleted than those in the plateau. There is a significant scatter in lithium abundances for any given mass M <= 1.1 M-circle dot. Conclusions. Our models qualitatively reproduce most of the features described above, although the predicted depletion of lithium is 0.45 dex smaller than observed for masses in the plateau region, i.e. between 1.1 and 1.28 solar masses. More work is clearly needed to accurately reproduce the observations. Despite hints that chromospheric activity and rotation play a role in lithium depletion, no firm conclusion can be drawn with the presently available data.

Enhancement of optically thick to thin line intensities in solar and stellar coronal plasmas through radiative transfer effects: an angularly resolved study.

An analysis of radiative transfer effects present in the Fe XV ion stage of solar and stellar coronal plasmas provides a general explanation of line radiation intensity enhancement above the optically thin limit. Full linearization radiation transfer is compared with the escape factor method and found to be in good agreement at the lower column densities. An angular study of the enhancement shows that symmetry factors are of great importance. This gives a possible reason for the indeterminate status of opacity in relation to coronal lines of distant stellar sources, where only emission integrated across the whole surface is detected.

Propagation of an isothermal shock in stellar envelopes

We have studied in this paper the propagation of an isothermal shock in the radiative envelopes of the Bosman-Crespin model for a hot star and Boury’s model for a giant star. A spherically symmetric disturbance is supposed to be originated at or outside the surface of the convective core. We have used Whitham’s rule to study the variation in the shock strength and the shock velocity after modifying it for inclusion of pressure, energy and flux of radiation. We find the shock increases in strength as it propagates through the envelopes of decreasing density, pressure and temperature. The velocity of the shock decreases for very weak initial shock strengths, for intermediate initial shock strength it first decreases and then increases, while for large initial shock strength, it always increases. This aspect of the problem throws some light on the stability of the models under consideration.

Monsoon sensitivity to aerosol direct radiative forcing in the community atmosphere model

Aerosol forcing remains a dominant uncertainty in climate studies. The impact of aerosol direct radiative forcing on Indian monsoon is extremely complex and is strongly dependent on the model, aerosol distribution and characteristics specified in the model, modelling strategy employed as well as on spatial and temporal scales. The present study investigates (i) the aerosol direct radiative forcing impact on mean Indian summer monsoon when a combination of quasi-realistic mean annual cycles of scattering and absorbing aerosols derived from an aerosol transport model constrained with satellite observed Aerosol Optical Depth (AOD) is prescribed, (ii) the dominant feedback mechanism behind the simulated impact of all-aerosol direct radiative forcing on monsoon and (iii) the relative impacts of absorbing and scattering aerosols on mean Indian summer monsoon. We have used CAM3, an atmospheric GCM (AGCM) that has a comprehensive treatment of the aerosol-radiation interaction. This AGCM has been used to perform climate simulations with three different representations of aerosol direct radiative forcing due to the total, scattering aerosols and black carbon aerosols. We have also conducted experiments without any aerosol forcing. Aerosol direct impact due to scattering aerosols causes significant reduction in summer monsoon precipitation over India with a tendency for southward shift of Tropical Convergence Zones (TCZs) over the Indian region. Aerosol forcing reduces surface solar absorption over the primary rainbelt region of India and reduces the surface and lower tropospheric temperatures. Concurrent warming of the lower atmosphere over the warm oceanic region in the south reduces the land-ocean temperature contrast and weakens the monsoon overturning circulation and the advection of moisture into the landmass. This increases atmospheric convective stability, and decreases convection, clouds, precipitation and associated latent heat release. Our analysis reveals a defining negative moisture-advection feedback that acts as an internal damping mechanism spinning down the regional hydrological cycle and leading to significant circulation changes in response to external radiative forcing perturbations. When total aerosol loading (both absorbing and scattering aerosols) is prescribed, dust and black carbon aerosols are found to cause significant atmospheric heating over the monsoon region but the aerosol-induced weakening of meridional lower tropospheric temperature gradient (leading to weaker summer monsoon rainfall) more than offsets the increase in summer-time rainfall resulting from the atmospheric heating effect of absorbing aerosols, leading to a net decrease of summer monsoon rainfall. Further, we have carried out climate simulations with globally constant AODs and also with the constant AODs over the extended Indian region replaced by realistic AODs. Regional aerosol radiative forcing perturbations over the Indian region is found to have impact not only over the region of loading but over remote tropical regions as well. This warrants the need to prescribe realistic aerosol properties in strategic regions such as India in order to accurately assess the aerosol impact.

An investigation of Fe XVI emission lines in solar and stellar extreme-ultraviolet and soft X-ray spectra

New fully relativistic calculations of radiative rates and electron impact excitation cross-sections for Fe XVI are used to determine theoretical emission-line ratios applicable to the 251-361 and 32-77 angstrom portions of the extreme-ultraviolet (EUV) and soft X-ray spectral regions, respectively. A comparison of the EUV results with observations from the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS) reveals excellent agreement between theory and experiment. However, for emission lines in the 32-49 angstrom portion of the soft X-ray spectral region, there are large discrepancies between theory and measurement for both a solar flare spectrum obtained with the X-Ray Spectrometer/Spectrograph Telescope (XSST) and for observations of Capella from the Low- Energy Transmission Grating Spectrometer (LETGS) on the Chandra X-ray Observatory. These are probably due to blending in the solar flare and Capella data from both first-order lines and from shorter wavelength transitions detected in second and third order. By contrast, there is very good agreement between our theoretical results and the XSST and LETGS observations in the 50-77 angstrom wavelength range, contrary to previous results. In particular, there is no evidence that the Fe XVI emission from the XSST flare arises from plasma at a much higher temperature than that expected for Fe XVI in ionization equilibrium, as suggested by earlier work.

Stellar Surface Magneto-Convection as a Source of Astrophysical Noise. I. Multi-component Parameterisation of Absorption Line Profiles

We outline our techniques to characterise photospheric granulation as an astrophysical noise source. A four component parameterisation of granulation is developed that can be used to reconstruct stellar line asymmetries and radial velocity shifts due to photospheric convective motions. The four components are made up of absorption line profiles calculated for granules, magnetic intergranular lanes, non-magnetic intergranular lanes, and magnetic bright points at disc centre. These components are constructed by averaging Fe I $6302 \mathrm{\AA}$ magnetically sensitive absorption line profiles output from detailed radiative transport calculations of the solar photosphere. Each of the four categories adopted are based on magnetic field and continuum intensity limits determined from examining three-dimensional magnetohydrodynamic simulations with an average magnetic flux of $200 \mathrm{G}$. Using these four component line profiles we accurately reconstruct granulation profiles, produced from modelling 12 x 12 Mm$^2$ areas on the solar surface, to within $\sim \pm$ 20 cm s$^{-1}$ on a $\sim$ 100 m s$^{-1}$ granulation signal. We have also successfully reconstructed granulation profiles from a $50 \mathrm{G}$ simulation using the parameterised line profiles from the $200 \mathrm{G}$ average magnetic field simulation. This test demonstrates applicability of the characterisation to a range of magnetic stellar activity levels.

Two-dimensional Monte Carlo simulations of H i line formation in massive young stellar object disc winds

Massive young stellar objects (YSOs) are powerful infrared Hi line emitters. It has been suggested that these lines form in an outflow from a disc surrounding the YSO. Here, new two-dimensional Monte Carlo radiative transfer calculations are described which test this hypothesis. Infrared spectra are synthesized for a YSO disc wind model based on earlier hydrodynamical calculations. The model spectra are in qualitative agreement with the observed spectra from massive YSOs, and therefore provide support for a disc wind explanation for the Hi lines. However, there are some significant differences: the models tend to overpredict the Bra/Br? ratio of equivalent widths and produce line profiles which are slightly too broad and, in contrast to typical observations, are double-peaked. The interpretation of these differences within the context of the disc wind picture and suggestions for their resolution via modifications to the assumed disc and outflow structure are discussed. © 2005 RAS.

Line-driven Disk Winds in Activ Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

The nebular spectra of SN 2012aw and constraints on stellar nucleosynthesis from oxygen emission lines

We present nebular-phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with non-local thermodynamic equilibrium radiative transfer calculations applied to ejecta from stellar evolution/explosion models. Our spectral synthesis models generally show good agreement with the ejecta from a M_ZAMS = 15 M_⊙progenitor star. The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14-18 M_⊙ range.We also demonstrate how the evolution of the oxygen cooling lines of [O I] λ5577, [O I] λ6300, and [O I] λ6364 can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to < 1 M_⊙, independent of the mixing in the ejecta. This constraint implies that any progenitor model of initial mass greater than 20 M_⊙ would be difficult to reconcile with the observed line strengths. A stellar progenitor of around M_ZAMS = 15 M_⊙ can consistently explain the directly measured luminosity of the progenitor star, the observed nebular spectra, and the inferred pre-supernova mass-loss rate.We conclude that there is still no convincing example of a Type IIP supernova showing the nucleosynthesis products expected from an M_ZAMS > 20 M_⊙ progenitor. © 2014 The Author. Published by Oxford University Press on behalf of the Royal Astronomical Society.

Modeling the Rossiter-McLaughlin Effect: Impact of the Convective Center-to-limb Variations in the Stellar Photosphere

bservations of the Rossiter–McLaughlin (RM) effect provide information on star–planet alignments, which can inform planetary migration and evolution theories. Here, we go beyond the classical RM modeling and explore the impact of a convective blueshift that varies across the stellar disk and non-Gaussian stellar photospheric profiles. We simulated an aligned hot Jupiter with a four-day orbit about a Sun-like star and injected center-to-limb velocity (and profile shape) variations based on radiative 3D magnetohydrodynamic simulations of solar surface convection. The residuals between our modeling and classical RM modeling were dependent on the intrinsic profile width and v sin i; the amplitude of the residuals increased with increasing v sin i and with decreasing intrinsic profile width. For slowly rotating stars the center-to-limb convective variation dominated the residuals (with amplitudes of 10 s of cm s−1 to ~1 m s−1); however, for faster rotating stars the dominant residual signature was due a non-Gaussian intrinsic profile (with amplitudes from 0.5 to 9 m s−1). When the impact factor was 0, neglecting to account for the convective center-to-limb variation led to an uncertainty in the obliquity of ~10°–20°, even though the true v sin i was known. Additionally, neglecting to properly model an asymmetric intrinsic profile had a greater impact for more rapidly rotating stars (e.g., v sin i = 6 km s−1) and caused systematic errors on the order of ~20° in the measured obliquities. Hence, neglecting the impact of stellar surface convection may bias star–planet alignment measurements and consequently theories on planetary migration and evolution.

Étude de l'influence de la perte de masse sur l'évolution d'étoiles de plusieurs types

Cette thèse a été réalisée en cotutelle. Pour la forme, Gérard Jasniewicz était mon codirecteur 'officiel' en France, bien que mon codirecteur était plutôt Olivier Richard qui m'a encadré lorsque j'étais en France.

Radiative capture of nucleons at astrophysical energies with single-particle states

Radiative capture of nucleons at energies of astrophysical interest is one of the most important processes for nucleosynthesis. The nucleon capture can occur either by a compound nucleus reaction or by a direct process. The compound reaction cross sections are usually very small, especially for light nuclei. The direct capture proceeds either via the formation of a single-particle resonance or a non-resonant capture process. In this work we calculate radiative capture cross sections and astrophysical S-factors for nuclei in the mass region A < 20 using single-particle states. We carefully discuss the parameter fitting procedure adopted in the simplified two-body treatment of the capture process. Then we produce a detailed list of cases for which the model works well. Useful quantities, such as spectroscopic factors and asymptotic normalization coefficients, are obtained and compared to published data. (C) 2010 Elsevier Inc. All rights reserved.

Observational effects of stellar rotation on Be stars' disks

Be stars are known to be fast rotators. At high rotation rates a profound modification of the radiation field reaching the circumstellar environment is expected. The origin of this modification is the decrease of the effective gravity on stellar surface leading to the stellar geometrical flattening and the gravity darkening effect predicted by Von Zeipel. Making use of the radiative transfer code HDUST we discuss the consequences of such stellar rotation on the structure of Be star disks based on the Viscous Decretion Disk model. Observational predictions are also made, as SED, IR-excess and Hydrogen line profiles. The modified illumination of the circumstellar disk generates significant changes in these quantities. Ascertaining these changes is useful to set some of the fundamental parameters of the Be system and to unveil the role of stellar rotation over the stellar evolution.

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-Stream Approximation

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.