Rapid Keck Spectroscopy of Cataclysmic Variables

We present an analysis of Rapid Keck Spectroscopy of the CVs AM Her (polar) and SS Cyg (dwarf nova). We decompose the spectra into constant and variable components and identify different types of variability in AM Her with different characteristic timescales. The variable flickering component of the accretion disc flux and the observational characteristics of a small flare in SS Cyg are isolated.

A regime-dependent parametrization of subgrid-scale cloud water content variability

The subgrid-scale spatial variability in cloud water content can be described by a parameter f called the fractional standard deviation. This is equal to the standard deviation of the cloud water content divided by the mean. This parameter is an input to schemes that calculate the impact of subgrid-scale cloud inhomogeneity on gridbox-mean radiative fluxes and microphysical process rates. A new regime-dependent parametrization of the spatial variability of cloud water content is derived from CloudSat observations of ice clouds. In addition to the dependencies on horizontal and vertical resolution and cloud fraction included in previous parametrizations, the new parametrization includes an explicit dependence on cloud type. The new parametrization is then implemented in the Global Atmosphere 6 (GA6) configuration of the Met Office Unified Model and used to model the effects of subgrid variability of both ice and liquid water content on radiative fluxes and autoconversion and accretion rates in three 20-year atmosphere-only climate simulations. These simulations show the impact of the new regime-dependent parametrization on diagnostic radiation calculations, interactive radiation calculations and both interactive radiation calculations and in a new warm microphysics scheme. The control simulation uses a globally constant f value of 0.75 to model the effect of cloud water content variability on radiative fluxes. The use of the new regime-dependent parametrization in the model results in a global mean which is higher than the control's fixed value and a global distribution of f which is closer to CloudSat observations. When the new regime-dependent parametrization is used in radiative transfer calculations only, the magnitudes of short-wave and long-wave top of atmosphere cloud radiative forcing are reduced, increasing the existing global mean biases in the control. When also applied in a new warm microphysics scheme, the short-wave global mean bias is reduced.

Spatial variability of liquid cloud and rain: observations and microphysical effects

Spatial variability of liquid cloud water content and rainwater content is analysed from three different observational platforms: in situ measurements from research aircraft, land-based remote sensing techniques using radar and lidar, and spaceborne remote sensing from CloudSat. The variance is found to increase with spatial scale, but also depends strongly on the cloud or rain fraction regime, with overcast regions containing less variability than broken cloud fields. This variability is shown to lead to large biases, up to a factor of 4, in both the autoconversion and accretion rates estimated at a model grid scale of ≈40 km by a typical microphysical parametrization using in-cloud mean values. A parametrization for the subgrid variability of liquid cloud and rainwater content is developed, based on the observations, which varies with both the grid scale and cloud or rain fraction, and is applicable for all model grid scales. It is then shown that if this parametrization of the variability is analytically incorporated into the autoconversion and accretion rate calculations, the bias is significantly reduced.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

Evolution of primordial black holes in a radiation and phantom energy environment

In this work we extend previous work on the evolution of a primordial black hole (PBH) to address the presence of a dark energy component with a super-negative equation of state as a background, investigating the competition between the radiation accretion, the Hawking evaporation and the phantom accretion, the latter two causing a decrease on black hole mass. It is found that there is an instant during the matter-dominated era after which the radiation accretion becomes negligible compared to the phantom accretion. The Hawking evaporation may become important again depending on a mass threshold. The evaporation of PBHs is quite modified at late times by these effects, but only if the generalized second law of thermodynamics is violated.

The old metal-poor open cluster ESO 92-SC05: accreted from a dwarf galaxy ?

The study of old open clusters outside the solar circle can bring constraints on formation scenarios of the outer disc. In particular, accretion of dwarf galaxies has been proposed as a likely mechanism in the area. We use BVI photometry for determining fundamental parameters of the faint open cluster ESO 92-SC05. Colour-magnitude diagrams are compared with Padova isochrones, in order to derive age, reddening and distance. We derive a reddening E(B - V) = 0.17, and an old age of similar to 6.0 Gyr. It is one of the rare open clusters known to be older than 5 Gyr. A metallicity of Z similar to 0.004 or [M/H] similar to -0.7 is found. The rather low metallicity suggests that this cluster might be the result of an accretion episode of a dwarf galaxy.

A multispectral view of the periodic events in eta Carinae

A full description of the 5.5-yr low excitation events in. Carinae is presented. We show that they are not as simple and brief as previously thought, but a combination of two components. The first, the slow variation component, is revealed by slow changes in the ionization level of circumstellar matter across the whole cycle and is caused by gradual changes in the wind wind collision shock-cone orientation, angular opening and gaseous content. The second, the collapse component, is restricted to around the minimum, and is due to a temporary global collapse of the wind-wind collision shock. High-energy photons (E > 16 eV) from the companion star are strongly shielded, leaving the Weigelt objects at low-ionization state for more than six months. High-energy phenomena are sensitive only to the collapse, low energy only to the slow variation and intermediate energies to both components. Simple eclipses and mechanisms effective only near periastron (e. g. shell ejection or accretion on to the secondary star) cannot account for the whole 5.5-yr cycle. We find anti-correlated changes in the intensity and the radial velocity of P Cygni absorption profiles in Fe II lambda 6455 and He I lambda 7065 lines, indicating that the former is associated to the primary and the latter to the secondary star. We present a set of light curves representative of the whole spectrum, useful for monitoring the next event (2009 January 11).

The extreme physical properties of the CoRoT-7b super-Earth

The search for rocky exoplanets plays an important role in our quest for extra-terrestrial life. Here, we discuss the extreme physical properties possible for the first characterised rocky super-Earth, CoRoT-7b (R(pl) = 1.58 +/- 0.10 R(Earth), M(pl) = 6.9 +/- 1.2 M(Earth)). It is extremely close to its star (a = 0.0171 AU = 4.48 R(st)), with its spin and orbital rotation likely synchronised. The comparison of its location in the (M(pl), R(pl)) plane with the predictions of planetary models for different compositions points to an Earth-like composition, even if the error bars of the measured quantities and the partial degeneracy of the models prevent a definitive conclusion. The proximity to its star provides an additional constraint on the model. It implies a high extreme-UV flux and particle wind, and the corresponding efficient erosion of the planetary atmosphere especially for volatile species including water. Consequently, we make the working hypothesis that the planet is rocky with no volatiles in its atmosphere, and derive the physical properties that result. As a consequence, the atmosphere is made of rocky vapours with a very low pressure (P <= 1.5 Pa), no cloud can be sustained, and no thermalisation of the planet is expected. The dayside is very hot (2474 +/- 71 K at the sub-stellar point) while the nightside is very cold (50-75 K). The sub-stellar point is as hot as the tungsten filament of an incandescent bulb, resulting in the melting and distillation of silicate rocks and the formation of a lava ocean. These possible features of CoRoT-7b could be common to many small and hot planets, including the recently discovered Kepler-10b. They define a new class of objects that we propose to name ""Lava-ocean planets"". (C) 2011 Elsevier Inc. All rights reserved.

Gas streaming motions towards the nucleus of M81

We present two-dimensional stellar and gaseous kinematics of the inner 120 x 250 pc2 of the LINER/Seyfert 1 galaxy M81, from optical spectra obtained with the Gemini Multi-Object Spectrograph (GMOS) integral field spectrograph on the Gemini-North telescope at a spatial resolution of approximate to 10 pc. The stellar velocity field shows circular rotation and, overall, is very similar to the published large-scale velocity field, but deviations are observed close to the minor axis which can be attributed to stellar motions possibly associated with a nuclear bar. The stellar velocity dispersion of the bulge is 162 +/- 15 km s-1, in good agreement with previous measurements and leading to a black hole mass of M(BH) = 5.5+3.6(-2.0) x 107 M(circle dot) based on the M(BH)-Sigma relationship. The gas kinematics is dominated by non-circular motions and the subtraction of the stellar velocity field reveals blueshifts of approximate to-100 km s-1 on the far side of the galaxy and a few redshifts on the near side. These characteristics can be interpreted in terms of streaming towards the centre if the gas is in the plane. On the basis of the observed velocities and geometry of the flow, we estimate a mass inflow rate in ionized gas of approximate to 4.0 x 10-3 M(circle dot) yr-1, which is of the order of the accretion rate necessary to power the LINER nucleus of M81. We have also applied the technique of principal component analysis (PCA) to our data, which reveals the presence of a rotating nuclear gas disc within approximate to 50 pc from the nucleus and a compact outflow, approximately perpendicular to the disc. The PCA combined with the observed gas velocity field shows that the nuclear disc is being fed by gas circulating in the galaxy plane. The presence of the outflow is supported by a compact jet seen in radio observations at a similar orientation, as well as by an enhancement of the [O i]/H alpha line ratio, probably resulting from shock excitation of the circumnuclear gas by the radio jet. With these observations we are thus resolving both the feeding - via the nuclear disc and observed gas inflow, and the feedback - via the outflow, around the low-luminosity active nucleus of M81.

Random primordial magnetic fields and the gas content of dark matter haloes

We recently predicted the existence of random primordial magnetic fields (RPMFs) in the form of randomly oriented cells with dipole-like structure with a cell size L(0) and an average magnetic field B(0). Here, we investigate models for primordial magnetic field with a similar web-like structure, and other geometries, differing perhaps in L(0) and B(0). The effect of RPMF on the formation of the first galaxies is investigated. The filtering mass, M(F), is the halo mass below which baryon accretion is severely depressed. We show that these RPMF could influence the formation of galaxies by altering the filtering mass and the baryon gas fraction of a halo, f(g). The effect is particularly strong in small galaxies. We find, for example, for a comoving B(0) = 0.1 mu G, and a reionization epoch that starts at z(s) = 11 and ends at z(e) = 8, for L(0) = 100 pc at z = 12, the f(g) becomes severely depressed for M < 10(7) M(circle dot), whereas for B(0) = 0 the f(g) becomes severely depressed only for much smaller masses, M < 10(5) M(circle dot). We suggest that the observation of M(F) and f(g) at high redshifts can give information on the intensity and structure of primordial magnetic fields.

Analytical solutions of accreting black holes immersed in a Lambda CDM model

The evolution of the mass of a black hole embedded in a universe filled with dark energy and cold dark matter is calculated in a closed form within a test fluid model in a Schwarzschild metric, taking into account the cosmological evolution of both fluids. The result describes exactly how accretion asymptotically switches from the matter-dominated to the Lambda-dominated regime. For early epochs, the black hole mass increases due to dark matter accretion, and on later epochs the increase in mass stops as dark energy accretion takes over. Thus, the unphysical behaviour of previous analyses is improved in this simple exact model. (C) 2010 Elsevier B.V. All rights reserved.

Modelling the secular evolution of migrating planet pairs

The subject of this paper is the secular behaviour of a pair of planets evolving under dissipative forces. In particular, we investigate the case when dissipative forces affect the planetary semimajor axes and the planets move inwards/outwards the central star, in a process known as planet migration. To perform this investigation, we introduce fundamental concepts of conservative and dissipative dynamics of the three-body problem. Based on these concepts, we develop a qualitative model of the secular evolution of the migrating planetary pair. Our approach is based on the analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces is invoked. We show that, under the assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the Mode I and Mode II stationary solutions of the conservative secular problem. The ultimate convergence and the evolution of the system along one of these secular modes of motion are determined uniquely by the condition that the dissipation rate is sufficiently smaller than the proper secular frequency of the system. We show that it is possible to reassemble the starting configurations and the migration history of the systems on the basis of their final states and consequently to constrain the parameters of the physical processes involved.

DIFFUSION OF MAGNETIC FIELD AND REMOVAL OF MAGNETIC FLUX FROM CLOUDS VIA TURBULENT RECONNECTION

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence reassures that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti-correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic field and density, which corresponds well to the observations of the interstellar media. While earlier studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for de-correlating magnetic field and density, we get the effect of permanent de-correlation with one fluid code, i.e., without invoking ambipolar diffusion. In addition, in the presence of gravity and turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus, the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent advection of heat obtained in an earlier study.

Precession and nutation in the eta Carinae binary system: evidence from the X-ray light curve

It is believed that eta Carinae is actually a massive binary system, with the wind-wind interaction responsible for the strong X-ray emission. Although the overall shape of the X-ray light curve can be explained by the high eccentricity of the binary orbit, other features like the asymmetry near periastron passage and the short quasi-periodic oscillations seen at those epochs have not yet been accounted for. In this paper we explain these features assuming that the rotation axis of eta Carinae is not perpendicular to the orbital plane of the binary system. As a consequence, the companion star will face eta Carinae on the orbital plane at different latitudes for different orbital phases and, since both the mass-loss rate and the wind velocity are latitude dependent, they would produce the observed asymmetries in the X-ray flux. We were able to reproduce the main features of the X-ray light curve assuming that the rotation axis of eta Carinae forms an angle of 29 degrees +/- 4 degrees with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with an amplitude of about 5 degrees and a period of about 22 days. The nutation parameters, as well as the precession of the apsis, with a period of about 274 years, are consistent with what is expected from the torques induced by the companion star.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS: THE MAGNETIC FIELD GEOMETRY AND THE INFLUENCE OF THE WIND ON GIANT PLANET MIGRATION

By means of numerical simulations, we investigate magnetized stellar winds of pre-main-sequence stars. In particular, we analyze under which circumstances these stars will present elongated magnetic features (e.g., helmet streamers, slingshot prominences, etc). We focus on weak-lined T Tauri stars, as the presence of the tenuous accretion disk is not expected to have strong influence on the structure of the stellar wind. We show that the plasma-beta parameter (the ratio of thermal to magnetic energy densities) is a decisive factor in defining the magnetic configuration of the stellar wind. Using initial parameters within the observed range for these stars, we show that the coronal magnetic field configuration can vary between a dipole-like configuration and a configuration with strong collimated polar lines and closed streamers at the equator (multicomponent configuration for the magnetic field). We show that elongated magnetic features will only be present if the plasma-beta parameter at the coronal base is beta(0) << 1. Using our self-consistent three-dimensional magnetohydrodynamics model, we estimate for these stellar winds the timescale of planet migration due to drag forces exerted by the stellar wind on a hot-Jupiter. In contrast to the findings of Lovelace et al., who estimated such timescales using the Weber and Davis model, our model suggests that the stellar wind of these multicomponent coronae are not expected to have significant influence on hot-Jupiters migration. Further simulations are necessary to investigate this result under more intense surface magnetic field strengths (similar to 2-3 kG) and higher coronal base densities, as well as in a tilted stellar magnetosphere.