Magnetic Doppler imaging is currently the most powerful method of interpreting high-resolution spectropolarimetric observations of stars. This technique has revealed the presence of unexpected small-scale magnetic fields on the surfaces of Ap stars. These studies were recently criticisied by Stift et al. (2012), who claimed that magnetic inversions are not robust and are undermined by neglecting a feedback on the Stokes line profiles from the local atmospheric structure in the regions of enhanced metal abundance. We show that Stift et al. misinterpreted published magnetic Doppler imaging results and neglected some of the most fundamental principles behind magnetic mapping. We demonstrate that the variation of atmospheric structure across the surface of a star with chemical spots affects the local continuum intensity but is negligible for the normalised local Stokes profiles. For the disk-integrated spectra of an Ap star with extreme abundance variations, we find that the assumption of a mean model atmosphere leads to moderate errors in Stokes I but is negligible for polarisation spectra. Employing a new magnetic inversion code, which incorporates the horizontal variation of atmospheric structure, we reconstructed new maps of magnetic field and Fe abundance for the Ap star alpha^2 CVn. The resulting distribution of chemical spots changes insignificantly compared to the previous modelling based on a single model atmosphere, while the magnetic field geometry does not change at all. This shows that the assertions by Stift et al. are exaggerated as a consequence of unreasonable assumptions and extrapolations, as well as methodological flaws and inconsistencies of their analysis. Our discussion proves that published magnetic inversions based on a mean stellar atmosphere are robust and reliable, and that the presence of small-scale magnetic field structures on the surfaces of Ap stars is real.

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called “{\alpha}2″ LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs.

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called “{\alpha}2″ LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs.

We study the impact of the late time dynamical evolution of ejecta from core-collapse supernovae on \nu p-process nucleosynthesis. Our results are based on hydrodynamical simulations of neutrino wind ejecta. Motivated by recent two-dimensional wind simulations, we vary the dynamical evolution during the \nu p-process and show that final abundances strongly depend on the temperature evolution. When the expansion is very fast, there is not enough time for antineutrino absorption on protons to produce enough neutrons to overcome the \beta-decay waiting points and no heavy elements beyond A=64 are produced. The wind termination shock or reverse shock dramatically reduces the expansion speed of the ejecta. This extends the period during which matter remains at relatively high temperatures and is exposed to high neutrino fluxes, thus allowing for further (p,\gamma) and (n,p) reactions to occur and to synthesize elements beyond iron. We find that the \nu p-process starts to efficiently produce heavy elements only when the temperature drops below ~3 GK. At higher temperatures, due to the low alpha separation energy of 60Zn (S_{\alpha} = 2.7 MeV) the reaction 59Cu(p,\alpha)56Ni is faster than the reaction 59Cu(p,\gamma)60Zn. This results in the closed NiCu cycle that we identify and discuss here for the first time. We also investigate the late phase of the \nu p-process when the temperatures become too low to maintain proton captures. Depending on the late neutron density, the evolution to stability is dominated by \beta decays or by (n,\gamma) reactions. In the latter case, the matter flow can even reach the neutron-rich side of stability and the isotopic composition of a given element is then dominated by neutron-rich isotopes.

We discuss the peculiar nature of the nucleus of M94 (NGC 4736) in the context of new measurements of the broad H_alpha emission from HST-STIS observations. We show that this component is unambiguously associated with the high-resolution X-ray, radio, and variable UV sources detected at the optical nucleus of this galaxy. These multi-wavelength observations suggest that NGC 4736 is one of the least luminous broad-line (type 1) LINERs, with Lbol = 2.5 \times 10^40 erg/s. This LINER galaxy has also possibly the least luminous broad line region known (LH_alpha =2.2\times10^37 erg/s). We compare black hole mass estimates of this system to the recently measured ~7 \times 10^6 M_sun dynamical black hole mass measurement. The fundamental plane and M-sigma relationship roughly agree with the measured black hole mass, while other accretion based estimates (the M-FWHM(H_alpha) relation, empirical correlation of BH mass with high-ionization mid IR emission lines, and the X-ray excess variance) provide much lower estimates (~10^5 M_sun). An energy budget test shows that the AGN in this system may be deficient in ionizing radiation relative to the observed emission-line activity. This deficiency may result from source variability or the superposition of multiple sources including supernovae.

We discuss the peculiar nature of the nucleus of M94 (NGC 4736) in the context of new measurements of the broad H_alpha emission from HST-STIS observations. We show that this component is unambiguously associated with the high-resolution X-ray, radio, and variable UV sources detected at the optical nucleus of this galaxy. These multi-wavelength observations suggest that NGC 4736 is one of the least luminous broad-line (type 1) LINERs, with Lbol = 2.5 \times 10^40 erg/s. This LINER galaxy has also possibly the least luminous broad line region known (LH_alpha =2.2\times10^37 erg/s). We compare black hole mass estimates of this system to the recently measured ~7 \times 10^6 M_sun dynamical black hole mass measurement. The fundamental plane and M-sigma relationship roughly agree with the measured black hole mass, while other accretion based estimates (the M-FWHM(H_alpha) relation, empirical correlation of BH mass with high-ionization mid IR emission lines, and the X-ray excess variance) provide much lower estimates (~10^5 M_sun). An energy budget test shows that the AGN in this system may be deficient in ionizing radiation relative to the observed emission-line activity. This deficiency may result from source variability or the superposition of multiple sources including supernovae.

We study baryon asymmetry generation originated by the leptogenesis in the presence of hypermagnetic fields in the early Universe plasma before the electroweak phase transition (EWPT). For the simplest Chern-Simons (CS) wave configuration of hypermagnetic field we find baryon asymmetry growth when the hypermagnetic field value changes due to alpha^2-dynamo and the lepton asymmetry rises due to the Abelian anomaly. We solve the corresponding integro-differential equations for lepton asymmetries describing such self-consistent dynamics for lepto- and baryogenesis in the two scenarios : (i) when a primordial lepton asymmetry sits in right electrons e_R; and (ii) when, in addition to e_R, a left lepton asymmetry for e_L and nu_eL arises due to chirality flip reactions provided by inverse Higgs decays at the temperatures, T < T_RL ~ 10 TeV. We find that the baryon asymmetry of the Universe (BAU) rises very fast through such leptogenesis, especially, in strong hypermagnetic fields. Varying (decreasing) the CS wave number parameter k_0 < 10^-7 T_EW one can recover the observable value of BAU, eta_B ~ 10^-9, where k_0 = 10^-7 T_EW corresponds to the maximum value for CS wave number surviving ohmic dissipation of hypermagnetic field. In the scenario (ii) one predicts the essential difference of the lepton numbers of right- and left electrons at EWPT time, L_eR – L_eL ~ (mu_eR – mu_eL) / T_EW = Delta mu / T_EW ~ 10^-5 that can be used as an initial condition for chiral asymmetry after EWPT.

We determine the low-redshift field galaxy stellar mass function (GSMF) using an area of 143 deg^2 from the first three years of the Galaxy And Mass Assembly (GAMA) survey. The magnitude limits of this redshift survey are r < 19.4 mag over two thirds and 19.8 mag over one third of the area. The GSMF is determined from a sample of 5210 galaxies using a density-corrected maximum volume method. This efficiently overcomes the issue of fluctuations in the number density versus redshift. With H_0 = 70, the GSMF is well described between 10^8 and 10^11.5 Msun using a double Schechter function with mass^* = 10^10.66 Msun, phi_1^* = 3.96 x 10^-3 Mpc^-3, alpha_1 = -0.35, phi_2^* = 0.79 x 10^-3 Mpc^-3 and alpha_2 = -1.47. This result is more robust to uncertainties in the flow-model corrected redshifts than from the shallower Sloan Digital Sky Survey main sample (r < 17.8 mag). The upturn in the GSMF is also seen directly in the i-band and K-band galaxy luminosity functions. Accurately measuring the GSMF below 10^8 Msun is possible within the GAMA survey volume but as expected requires deeper imaging data to address the contribution from low surface-brightness galaxies.

Magnetic helicity fluxes in turbulently driven alpha^2 dynamos are studied to demonstrate their ability to alleviate catastrophic quenching. A one-dimensional mean-field formalism is used to achieve magnetic Reynolds numbers of the order of 10^5. We study both diffusive magnetic helicity fluxes through the mid-plane as well as those resulting from the recently proposed alternate dynamic quenching formalism. By adding shear we make a parameter scan for the critical values of the shear and forcing parameters for which dynamo action occurs. For this $\alpha\Omega$ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.

When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean-field equations are solved for oscillatory alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.

Observations in the near-infrared (NIR) allow a detailed study of young stellar clusters in grand-design spiral galaxies which in visual bands often are highly obscured by dust lanes along the arms. Deep JHK-maps of 10 spirals were obtained with HAWK-I/VLT. Data for NGC 2997 are presented here to illustrate the general results for the sample. The (H-K)-(J-H) diagrams suggest that most stellar clusters younger than 7 Myr are significantly attenuated by dust with visual extinctions reaching 7 mag. A gap between younger and older cluster complexes in the (J-K)-Mk diagram indicates a rapid reduction of extinction around 7 Myr possibly due to expulsion of dust and gas after supernovae explosions. The cluster luminosity function is consistent with a power law with an exponent alpha ~ 2. Cluster luminosities of Mk = -15 mag are reached, corresponding to masses close to 10^6 Mo, with no indication of a cut-off. Their azimuthal angles relative to the main spiral arms show that the most massive clusters are formed in the arm regions while fainter ones also are seen between the arms. Older clusters are more uniformly distribution with a weaker modulation relative to the arms.

We present global solutions of optically thin, two-temperature black hole accretion disks incorporating magnetic fields. We assume that the {\pi}{\phi}-component of Maxwell stress is proportional to the total pressure, and prescribe the radial dependence of the magnetic flux advection rate in order to complete the set of basic equations. We have obtained the magnetically supported (low-{\beta}) disk solutions whose luminosity exceeds the maximum luminosity for an advection-dominated accretion flow (ADAF), L > 0.4 {\alpha}^2 L_Edd. The accretion flow is composed of the outer ADAF, a luminous hot accretion flow (LHAF) inside the transition layer from the outer ADAF to the low-{\beta} disk, the low-{\beta} disk, and the inner ADAF. The low-{\beta} disk region becomes wider as the mass accretion rate increases further. In the low-{\beta} disk, the magnetic heating balances the radiative cooling, and the electron temperature decreases from ~ 10^9.5 K to ~ 10^8 K as the luminosity increases. These results are consistent with the anti-correlation between the energy cutoff in X-ray spectra (hence the electron temperature) and the luminosity when L > 0.1 L_Edd, observed in the bright/hard state during the bright hard-to-soft transitions of transient outbursts in galactic black hole candidates.

We constrain the densities of Earth- to Neptune-size planets around very cool (Te =3660-4660K) Kepler stars by comparing 1202 Keck/HIRES radial velocity measurements of 150 nearby stars to a model based on Kepler candidate planet radii and a power-law mass-radius relation. Our analysis is based on the presumption that the planet populations around the two sets of stars are the same. The model can reproduce the observed distribution of radial velocity variation over a range of parameter values, but, for the expected level of Doppler systematic error, the highest Kolmogorov-Smirnov probabilities occur for a power-law index alpha ~ 4, indicating that rocky-metal planets dominate the planet population in this size range. A single population of gas-rich, low-density planets with alpha = 2 is ruled out unless our Doppler errors are >= 5m/s, i.e., much larger than expected based on observations and stellar chromospheric emission. If small planets are a mix of gamma rocky planets (alpha = 3.85) and 1-gamma gas-rich planets (alpha = 2), then gamma > 0.5 unless Doppler errors are >=4 m/s. Our comparison also suggests that Kepler’s detection efficiency relative to ideal calculations is less than unity. One possible source of incompleteness is target stars that are misclassified subgiants or giants, for which the transits of small planets would be impossible to detect. Our results are robust to systematic effects, and plausible errors in the estimated radii of Kepler stars have only moderate impact.

We constrain the densities of Earth- to Neptune-size planets around very cool (3660-4660 K) Kepler stars by comparing 1332 Keck/HIRES radial velocity measurements of 167 nearby stars to a model based on Kepler candidate planet radii and a power-law mass-radius relation. Our analysis is based on the presumption that the planet populations around the two sets of stars are the same. The model can reproduce the observed distribution of radial velocity variation over a range of parameter values, but, for the expected magnitude of stellar jitter, the highest Kolmogorov-Smirnov probabilities occur for a power-law index alpha ~ 4, and our findings indicate that rocky-metal planets dominate the planet population in this size range. A dominant population of gas-rich, low-density planets with alpha 3 m/s) than expected and Kepler observations are severely (>80%) incomplete. The analysis limits the fraction of alpha=2 planets to <50% at 99% confidence. The comparison also suggests that that Kepler observations of these stars are incomplete compared to ideal calculations, and that many stars may have underestimated radii and/or are misclassified subgiants or giants, around which small planets would be almost impossible to detect by transit.

We compare the redshifts of neutral carbon and carbon monoxide in the redshifted sources in which the fine structure transition of neutral carbon, [CI], has been detected, in order to measure space-time variation of the fundamental constants. Comparison with the CO rotational lines measures gives the same combination of constants obtained from the comparison fine structure line of singly ionised carbon, [CII]. However, neutral carbon has the distinct advantage that it may be spatially coincident with the carbon monoxide, whereas [CII] could be located in the diffuse medium between molecular clouds, and so any comparison with CO could be dominated by intrinsic velocity differences. Using [CI], we obtain a mean variation of dF/F = (-3.6 +/- 8.5) x 10^-5, over z = 2.3 – 4.1, for the eight [CI] systems, which degrades to (-1.5+/- 11) x 10^-5, over z = 2.3 – 6.4 when the two [CII] systems are included. That is, zero variation over look-back times of 10.8-12.8 Gyr. However, the latest optical results indicate a spatial variation in alpha, which describes a dipole and we see the same direction in dF/F. This trend is, however, due to a single source for which the [CI] spectrum is of poor quality. This also applies to one of the two [CII] spectra previously used to find a zero variation in alpha^2/mu. Quantifying this, we find an anti-correlation between |dF/F| and the quality of the carbon detection, as measured by the spectral resolution, indicating that the typical values of >50 km/s, used to obtain a detection, are too coarse to reliably measure changes in the constants. From the fluxes of the known z > 1 CO systems, we predict that current instruments are incapable of the sensitivities required to measure changes in the constants through the comparison of CO and carbon lines. We therefore discuss in detail the use of ALMA for such an undertaking … ABRIDGED

We compare the redshifts of neutral carbon and carbon monoxide in the redshifted sources in which the fine structure transition of neutral carbon, [CI], has been detected, in order to measure space-time variation of the fundamental constants. Comparison with the CO rotational lines measures gives the same combination of constants obtained from the comparison fine structure line of singly ionised carbon, [CII]. However, neutral carbon has the distinct advantage that it may be spatially coincident with the carbon monoxide, whereas [CII] could be located in the diffuse medium between molecular clouds, and so any comparison with CO could be dominated by intrinsic velocity differences. Using [CI], we obtain a mean variation of dF/F = (-3.6 +/- 8.5) x 10^-5, over z = 2.3 – 4.1, for the eight [CI] systems, which degrades to (-1.5+/- 11) x 10^-5, over z = 2.3 – 6.4 when the two [CII] systems are included. That is, zero variation over look-back times of 10.8-12.8 Gyr. However, the latest optical results indicate a spatial variation in alpha, which describes a dipole and we see the same direction in dF/F. This trend is, however, due to a single source for which the [CI] spectrum is of poor quality. This also applies to one of the two [CII] spectra previously used to find a zero variation in alpha^2/mu. Quantifying this, we find an anti-correlation between |dF/F| and the quality of the carbon detection, as measured by the spectral resolution, indicating that the typical values of >50 km/s, used to obtain a detection, are too coarse to reliably measure changes in the constants. From the fluxes of the known z > 1 CO systems, we predict that current instruments are incapable of the sensitivities required to measure changes in the constants through the comparison of CO and carbon lines. We therefore discuss in detail the use of ALMA for such an undertaking … ABRIDGED

Theoretical predictions for the ensemble quasar structure function are tested using multi-epoch observations of Stripe 82 collected by the Sloan Digital Sky Survey. We reanalyze the entire available volume of the g-band imaging data using difference image photometry and build high quality light curves for 7562 spectroscopically confirmed quasars. Our structure function includes ~4.8×10^6 pairs of measurements and covers a wide range of time lags between 3 days and 6.9 years in the quasar rest frame. A broken power-law fit to this data shows the presence of two slopes alpha_1=0.33 and alpha_2=0.79 with the break at ~42 days. The structure function compiled using only flux increases is slightly lower than that for variations of the opposite sign, revealing a slight asymmetry between the leading and trailing edge of a typical flare. The reality of these features is confirmed with monte-carlo simulations. We give simple interpretation of the results in the frames of existing theoretical models.

Theoretical predictions for the ensemble quasar structure function are tested using multi-epoch observations of Stripe 82 collected by the Sloan Digital Sky Survey. We reanalyze the entire available volume of the g-band imaging data using difference image photometry and build high quality light curves for 7562 spectroscopically confirmed quasars. Our structure function includes ~4.8×10^6 pairs of measurements and covers a wide range of time lags between 3 days and 6.9 years in the quasar rest frame. A broken power-law fit to this data shows the presence of two slopes alpha_1=0.33 and alpha_2=0.79 with the break at ~42 days. The structure function compiled using only flux increases is slightly lower than that for variations of the opposite sign, revealing a slight asymmetry between the leading and trailing edge of a typical flare. The reality of these features is confirmed with monte-carlo simulations. We give simple interpretation of the results in the frames of existing theoretical models.

The goals of this work are to provide NIR polarimetry of the stellar sources in the central pc at the resolution of an 8m telescope for the first time, along with new insights into the nature of the known bright bow-shock sources. We use AO assisted observations with the ESO VLT in the H- and Ks-band, applying high precision photometric methods developed for crowded fields and a new polarimetric calibration for NACO to produce polarization maps of the central 3″x19″, in addition to spatially resolved polarimetry and a variability analysis on the extended sources. We find foreground polarization parallel to the Galactic plane, with averages of (4.6+/-0.6)% at 26{\deg}+/-6{\deg} (Ks) and (9.3+/-1.3)% at 20{\deg}+/-6{\deg} (H) in the center. At larger distances from the center, we find different polarization parameters: (7.5+/-1.0)% at 11{\deg}+/-6{\deg} (Ks) and (12.1+/-2.1)% at 13{\deg}+/-6{\deg} (H). p_H/p_Ks peaks at 1.9+/-0.4, with a power law index for the wavelength dependency of alpha = 2.4+/-0.7. This also varies over the FOV, with higher values in the center, indicating local effects on the total polarization, possibly dichroic extinction by Northern Arm dust. The two extended sources IRS 21 and 1W show similar intrinsic polarization degrees of 6.1% resp. 7.8% (Ks) and 6.9% (H, only 1W) at polarization angles coincident with previous NIR and MIR findings, both in total and resolved. The spatial polarization pattern of both sources points to scattering on aligned elongated dust grains as the major source of intrinsic polarization, and matches the known orientation of the magnetic field. Our data also allows us to separate the bow shock of IRS 21 from the central source for the first time in the Ks-band, with the apex north of the center and a standoff distance of ~400 AU, matching previous estimates. This source also shows a ~50% increase in flux in the NIR over several years.

Helicity and \alpha effect driven by the nonaxisymmetric Tayler instability of toroidal magnetic fields in stellar radiation zones are computed. In the linear approximation a purely toroidal field always excites pairs of modes with identical growth rates but with opposite helicity so that the net helicity vanishes. If the magnetic background field has a helical structure by an extra (weak) poloidal component then one of the modes dominates producing a net kinetic helicity anticorrelated to the current helicity of the background field. The mean electromotive force is computed with the result that the \alpha effect by the most rapidly growing mode has the same sign as the current helicity of the background field. The \alpha effect is found as too small to drive an \alpha^{2} dynamo but the excitation conditions for an \alpha\Omega dynamo can be fulfilled for weak poloidal fields. Moreover, if the dynamo produces its own \alpha effect by the magnetic instability then problems with its sign do not arise. For all cases, however, the \alpha effect shows an extremely strong concentration to the poles so that a possible \alpha\Omega dynamo might only work at the polar regions. Hence, the results of our linear theory lead to a new topological problem for the existence of large-scale dynamos in stellar radiation zones on the basis of the current-driven instability of toroidal fields.

We consider mean-field dynamo models with fluctuating \alpha effect, both with and without shear. The \alpha effect is chosen to be Gaussian white noise with zero mean and given covariance. We show analytically that the mean magnetic field does not grow, but, in an infinitely large domain, the mean-squared magnetic field shows exponential growth of the fastest growing mode at a rate proportional to the shear rate, which agrees with earlier numerical results of Yousef et al (2008) and recent analytical treatment by Heinemann et al (2011) who use a method different from ours. In the absence of shear, an incoherent \alpha^2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third and fourth order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian.

We describe the serendipitous discovery of a radio point source in a 618 MHz image of the supernova remnant(SNR) G76.9+1.0. The SNR has a bipolar structure and the point source is located near a faint bridge of emission joining the two lobes of emission. The point source was also detected in follow-up higher frequency(1170 MHz) observations. The spectral index for the point source obtained from the GMRT observations is alpha = -2.1. The steep spectrum, as well as the location of the point source near the centre of the SNR establish the fact that it is indeed the pulsar J2022+3842 associated with this SNR. Consistent with this, subsequent analysis of archival Chandra X-ray data shows a point source coincident with the radio point source, as well as diffuse extended X-ray emission surrounding the unresolved source. However, no pulsed emission was detected despite deep searches at both 610 MHz and 1160 MHz although pulsed emission has been seen at 2 GHz with the GBT. It appears that the most likely reason for not detecting the pulsed signal at the GMRT is temporal broadening: for the estimated DM towards this SNR, the pulse broadening time could be as large as tens of milliseconds. The diffuse X-ray emission is elongated along the same direction as the bipolar structure seen in the radio. We interpret the radio lobes as having been formed from an equatorial wind. Although direct detection of pulsed signal has not been possible, we show convincingly that sensitive, high-resolution, radio imaging at multiple frequencies is a useful method to search for pulsar candidates.

We present new X-ray observations of the A0p stars alpha^2 CVn (log Lx < 26.0 erg/s) and IQ Aur (log Lx = 29.6 erg/s) and find that their X-ray luminosities differ by at least three orders of magnitude. IQ Aur possesses a strong cool plasma component with X-ray emitting regions located well above the stellar surface, but also significant amounts of hot plasma. Further, a large X-ray flare is detected from IQ Aur, implying the presence of magnetic reconnection. Our comparison study of similar stars indicates that the occurrence of X-ray emission generated by magnetically channelled wind shocks (MCWS) strongly depends on stellar properties. X-ray emission is preferably generated by more luminous and massive objects such as IQ Aur. The MCWS scenario can consistently describe the X-ray emission of these A0p stars, assuming that the very strong magnetic confinement of the stellar wind has led to the build up of a rigidly rotating disk around the star, where magnetic reconnection and centrifugal breakout events occur.

We derive the extinction curve towards the Galactic Center from 1 to 19 micron. We use hydrogen emission lines of the minispiral observed by ISO-SWS and SINFONI. The extinction free flux reference is the 2 cm continuum emission observed by the VLA. Towards the inner 14″ * 20″ we find an extinction of A(2.166 micron)=2.62 +/- 0.11, with a power-law slope of alpha=-2.11 +/- 0.06 shortward of 2.8 micron, consistent with the average near infrared slope from the recent literature. At longer wavelengths, however, we find that the extinction is grayer than shortward of 2.8 micron. We find it is not possible to fit the observed extinction curve with a dust model consisting of pure carbonaceous and silicate grains only, and the addition of composite particles, including ices, is needed to explain the observations. Combining a distance dependent extinction with our distance independent extinction we derive the distance to the GC to be R_0=7.94 +/- 0.65 kpc. Towards Sgr A* (r<0.5") we obtain A_H=4.21 +/- 0.10, A_Ks=2.42 +/- 0.10 and A_L'=1.09 +/- 0.13.

We derive the extinction curve towards the Galactic Center from 1 to 19 micron. We use hydrogen emission lines of the minispiral observed by ISO-SWS and SINFONI. The extinction free flux reference is the 2 cm continuum emission observed by the VLA. Towards the inner 14″ * 20″ we find an extinction of A(2.166 micron)=2.62 +/- 0.11, with a power-law slope of alpha=-2.11 +/- 0.06 shortward of 2.8 micron, consistent with the average near infrared slope from the recent literature. At longer wavelengths, however, we find that the extinction is grayer than shortward of 2.8 micron. We find it is not possible to fit the observed extinction curve with a dust model consisting of pure carbonaceous and silicate grains only, and the addition of composite particles, including ices, is needed to explain the observations. Combining a distance dependent extinction with our distance independent extinction we derive the distance to the GC to be R_0=7.94 +/- 0.65 kpc. Towards Sgr A* (r<0.5") we obtain A_H=4.21 +/- 0.10, A_Ks=2.42 +/- 0.10 and A_L'=1.09 +/- 0.13.

We develop a general method to fit the planetary distribution function (PLDF) to exoplanet survey data. This maximum likelihood method accommodates more than one planet per star and any number of planet or target star properties. Application to \Kepler data relies on estimates of the efficiency of discovering transits around Solar type stars by Howard et al. (2011). These estimates are shown to agree with theoretical predictions for an ideal transit survey. Using announced \Kepler planet candidates, we fit the PLDF as a joint powerlaw in planet radius, down to 0.5 R_Eart, and orbital period, up to 50 days. The estimated number of planets per star in this sample is ~ 0.7 — 1.4, where the broad range covers systematic uncertainties in the detection efficiency. To analyze trends in the PLDF we consider four planet samples, divided between shorter and longer periods at 7 days and between large and small radii at 3 R_Earth. At longer periods, the size distribution of the small planets, with index \alpha = -1.2 \pm 0.2 steepens to \alpha = -2.0 \pm 0.2 for the larger planet sample. For shorter periods, the opposite is seen: smaller planets follow a steep powerlaw, \alpha = -1.9 \pm 0.2 that is much shallower, \alpha = -0.7 \pm 0.2 at large radii. The observed deficit of intermediate-sized planets at the shortest periods may arise from the evaporation and sublimation of Neptune and Saturn-like planets. If the trend and explanation hold, it would be spectacular observational confirmation of the core accretion and migration hypotheses, and allow refinement of these theories.

We present a study of HII regions in M51 using the Hubble Space Telescope ACS images taken as part of the Hubble Heritage Program. We have catalogued about 19,600 HII regions in M51 with Ha luminosity in the range of L = 10^{35.5} erg/s to 10^{39.0} erg/s. The Ha luminosity function of HII regions (HII LF) in M51 is well represented by a double power law with its index alpha=-2.25\pm0.02 for the bright part and alpha=-1.42\pm0.01 for the faint part, separated at a break point L= 10^{37.1} erg/s. This break was not found in previous studies of M51 HII regions. Comparison with simulated HII LFs suggests that this break is caused by the transition of HII region ionizing sources, from low mass clusters (with ~ 10^3 M_sun, including several OB stars) to more massive clusters (including several tens of OB stars). The HII LFs with L < 10^{37.1} erg/s are found to have different slopes for different parts in M51: the HII LF for the interarm region is steeper than those for the arm and the nuclear regions. This observed difference in HII LFs can be explained by evolutionary effects that HII regions in the interarm region are relatively older than those in the other parts of M51.

Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar dynamos.

We have discovered a new rapidly oscillating Ap star among the Kepler Mission target stars, KIC 10195926. This star shows two pulsation modes with periods that are amongst the longest known for roAp stars at 17.1 min and 18.1 min, indicating that the star is near the terminal age main sequence. The principal pulsation mode is an oblique dipole mode that shows a rotationally split frequency septuplet that provides information on the geometry of the mode. The secondary mode also appears to be a dipole mode with a rotationally split triplet, but we are able to show within the improved oblique pulsator model that these two modes cannot have the same axis of pulsation. This is the first time for any pulsating star that evidence has been found for separate pulsation axes for different modes. The two modes are separated in frequency by 55 microHz, which we model as the large separation. The star is an alpha^2 CVn spotted magnetic variable that shows a complex rotational light variation with a period of Prot = 5.68459 d. For the first time for any spotted magnetic star of the upper main sequence, we find clear evidence of light variation with a period of twice the rotation period; i.e. a subharmonic frequency of $\nu_{\rm rot}/2$. We propose that this and other subharmonics are the first observed manifestation of torsional modes in an roAp star. From high resolution spectra we determine Teff = 7400 K, log g = 3.6 and v sin i = 21 km/s. We have found a magnetic pulsation model with fundamental parameters close to these values that reproduces the rotational variations of the two obliquely pulsating modes with different pulsation axes. The star shows overabundances of the rare earth elements, but these are not as extreme as most other roAp stars. The spectrum is variable with rotation, indicating surface abundance patches.

Recent theoretical work has suggested that Lyman-alpha nebulae could be substantially polarized in the Lyman-alpha emission line, depending on the geometry, kinematics, and powering mechanism at work. Polarization observations can therefore provide a useful constraint on the source of ionization in these systems. In this Letter, we present the first Lyman-alpha polarization measurements for a giant Lyman-alpha nebula at z~2.656. We do not detect any significant linear polarization of the Lyman-alpha emission: P_{Lyman-alpha}=2.6+/-2.8% (corrected for statistical bias) within a single large aperture. The current data also do not show evidence for the radial polarization gradient predicted by some theoretical models. These results rule out singly scattered Lyman-alpha (e.g., from the nearby AGN) and may be inconsistent with some models of backscattering in a spherical outflow. However, the effects of seeing, diminished signal-to-noise ratio, and angle averaging within radial bins make it difficult to put strong constraints on the radial polarization profile. The current constraints may be consistent with higher density outflow models, spherically symmetric infall models, photoionization by star formation within the nebula or the nearby AGN, resonant scattering, or non-spherically symmetric cold accretion (i.e., along filaments). Higher signal-to-noise ratio data probing to higher spatial resolution will allow us to harness the full diagnostic power of polarization observations in distinguishing between theoretical models of giant Lyman-alpha nebulae.

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the {\alpha} – {\Omega} mode is generally expected to dominate because it is the fastest growing. In an {\alpha} – {\Omega} dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory {\alpha} – {\Omega} type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both {\alpha} – {\Omega} and {\alpha}2 dynamo action is possible. Results. After initial dominance and saturation, the {\alpha} – {\Omega} mode is found to be destroyed by an orthogonal {\alpha}2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of {\alpha} – {\Omega} dynamos by its {\alpha}2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of {\alpha} – {\Omega} dynamos from an {\alpha}2 dynamo has not yet been found.

Context: Large-scale magnetic fields resulting from hydromagnetic dynamo action may differ substantially in their time dependence. Cyclic field variations, characteristic for the solar magnetic field, are often explained by an important omega-effect, i.e. by the stretching of field lines due to strong differential rotation. Aims: The dynamo mechanism of a convective, oscillatory dynamo model is investigated. Methods: We solve the MHD-equations for a conducting Boussinesq fluid in a rotating spherical shell. For a resulting oscillatory model, dynamo coefficients have been computed with the help of the so-called test-field method. Subsequently, these coefficients have been used in a mean-field calculation in order to explore the underlying dynamo mechanism. Results: Although the rather strong differential rotation present in this model influences the magnetic field, the omega-effect alone is not responsible for its cyclic time variation. If the omega-effect is suppressed, the resulting alpha^2-dynamo remains oscillatory. Surprisingly, the corresponding alpha-omega dynamo leads to a non-oscillatory magnetic field. Conclusions: The assumption of an alpha-omega mechanism does not explain the occurrence of magnetic cycles satisfactorily.

Over the past years observations of young and populous star clusters have shown that the stellar initial mass function (IMF) can be conveniently described by a two-part power-law with an exponent alpha_2 = 2.3 for stars more massive than about 0.5 Msol and an exponent of alpha_1 = 1.3 for less massive stars. A consensus has also emerged that most, if not all, stars form in stellar groups and star clusters, and that the mass function of these can be described as a power-law (the embedded cluster mass function, ECMF) with an exponent beta ~2. These two results imply that the integrated galactic IMF (IGIMF) for early-type stars cannot be a Salpeter power-law, but that they must have a steeper exponent. An application to star-burst galaxies shows that the IGIMF can become top-heavy. This has important consequences for the distribution of stellar remnants and for the chemo-dynamical and photometric evolution of galaxies. In this contribution the IGIMF theory is described, and the accompanying contribution by Pflamm-Altenburg, Weidner & Kroupa (this volume) documents the applications of the IGIMF theory to galactic astrophysics.

Aims: GRB 050502B is well known for the very bright flare displayed in its X-ray light curve. Despite extensive studies, however, the optical light curve has never been discussed and its redshift is unconstrained. Possible correlations between optical and X-ray data are analysed. Methods: Photometric data from TNG in the R and I bands were used to compare the optical afterglow with the X-ray light curve. The HyperZ package and a late time VLT host observation were used to derive redshift estimates. Results: The I-band afterglow decay followed a power-law of index {\alpha} = 2.1 $\pm$ 0.6, after a late break at ~ $1.3 \times 10^5$ s. The R – I color is remarkably red and the broadband spectral index {\beta}_OX = $0.9 \pm 0.1$ is consistent with the X-ray spectral slope {\beta}_X. Although a photometric redshift of z > 4 is the most conservative result to consider, a photometric redshift of z = 5.2 $\pm$ 0.3 is suggested with no extinction in the host, based on which an isotropic energy E_{\gamma},iso = $(3.8 \pm 0.7) \times 10^52$ erg and a jet opening angle {\theta} ~ $3.7^{\circ}$ are subsequently derived. Conclusions: The combined X-ray and optical data suggest an achromatic break, which we interpret as a jet break. The post jet break slope obeys roughly the closure relation for the jet slow cooling model. Because of the afterglow’s very red color, in order for the redshift to be low (z < 1), extinction, if present in the host, must be significantly high. Since the optical-to-X-ray index is consistent with the X-ray spectrum, and there is no XRT evidence for excess N_H, GRB 050502B was likely at high redshift.

The detection of a spatial variation of the fine-structure constant, alpha, based on study of quasar absorption systems has recently been reported. The physics that causes this alpha-variation should have other observable manifestations, and this motivates us to look for complementary astrophysical effects. In this paper we propose a method to test whether spatial variation of fundamental constants existed during the epoch of big bang nucleosynthesis. Using existing measurements of primordial deuterium abundance we find very weak indications that such a signature might exist, but the paucity of measurements precludes any firm conclusion. We also examine existing quasar absorption spectra data that are sensitive to variation of the electron-to-proton mass ratio, mu, and x = (alpha^2 mu g_p) for spatial variation.

We introduce a new class of models of chaotic inflation inspired by the superconformal approach to supergravity. This class of models allows a functional freedom of choice of the inflaton potential V = |f(\phi)|^2. The simplest model of this type has a quadratic potential m^2\phi^2/2. Another model describes an inflaton field with the standard symmetry breaking potential \lambda^2 (\phi^2-v^2)^2. Depending on the value of v and on initial conditions for inflation, the spectral index n_s may take any value from 0.97 to 0.93, and the tensor-to-scalar ratio r may span the interval form 0.3 to 0.01. A generalized version of this model has a potential \lambda^2 (\phi^\alpha-v^\alpha)^2. At large \phi and \alpha > 0, this model describes chaotic inflation with the power law potential \phi^{2\alpha}. For \alpha < 0, this potential describes chaotic inflation with a potential which becomes flat in the large field limit. We further generalize these models by introducing a nonminimal coupling of the inflaton field to gravity. The mechanism of moduli stabilization used in these models allows to improve and generalize several previously considered models of chaotic inflation in supergravity.

The physical mechanism behind the TeV gamma-ray source observed at the centre of the Galaxy is still unknown. One intriguing possibility is that the accretion flow onto the central supermassive black hole is responsible for accelerating protons to TeV energies which then diffuse outward to interact with molecular gas at distances of ~1 pc. Here, we build on our earlier detailed calculations of the proton transport to consider the time and energy dependence of the TeV signal following a burst of particle acceleration at Sgr A*. We find that, due to the strong energy dependence of the proton diffusion, any variability in the particle acceleration rate will only be visible in the TeV signal after a delay of ~ 10 yrs, and only at energies >~ 10 TeV. If the accelerator is long-lived it must have been running for at least 10^4 yrs and have a hard proton injection spectrum of \alpha=0.7 (where dn/dE_inj \propto E^{-\alpha}) in order to produce the correct amount of high energy gamma-ray flux. This rapid diffusion of high energy protons also rules out the possibility that the observed TeV source is directly related to the period of increased activity of Sgr A* that ended ~ 100 yrs ago. However, a good fit to the observed H.E.S.S. data was found with \alpha=2.7 for the scenario of a brief (~few year long) burst of particle acceleration that occurred ~ 10 yrs ago. If such bursts are common then they will keep the TeV source energised and will likely produce spectral variability at >~ 10 TeV on <~ 5 yr timescales. This model also implies that particle acceleration may be an important mechanism in reducing the radiative efficiency of weakly accreting black holes.

The newly installed Wide Field Camera 3 (WFC3) on the Hubble Space Telescope has been used to obtain multi-band images of the nearby spiral galaxy M83. These new observations are the deepest and highest resolution images ever taken of a grand-design spiral, particularly in the near ultraviolet, and allow us to better differentiate compact star clusters from individual stars and to measure the luminosities of even faint clusters in the U band. We find that the luminosity function for clusters outside of the very crowded starburst nucleus can be approximated by a power law, dN/dL \propto L^{alpha}, with alpha = -2.04 +/- 0.08, down to M_V ~ -5.5. We test the sensitivity of the luminosity function to different selection techniques, filters, binning, and aperture correction determinations, and find that none of these contribute significantly to uncertainties in alpha. We estimate ages and masses for the clusters by comparing their measured UBVI,Halpha colors with predictions from single stellar population models. The age distribution of the clusters can be approximated by a power-law, dN/dt propto t^{gamma}, with gamma=-0.9 +/- 0.2, for M > few x 10^3 Msun and t < 4×10^8 yr. This indicates that clusters are disrupted quickly, with ~80-90% disrupted each decade in age over this time. The mass function of clusters over the same M-t range is a power law, dN/dM propto M^{beta}, with beta=-1.94 +/- 0.16, and does not have bends or show curvature at either high or low masses. Therefore, we do not find evidence for a physical upper mass limit, M_C, or for the earlier disruption of lower mass clusters when compared with higher mass clusters, i.e. mass-dependent disruption. We briefly discuss these implications for the formation and disruption of the clusters.

Radial velocity (RV) observations of an exoplanet system giving a value of M_T sin(i) condition (ie. give information about) not only the planet’s true mass M_T but also the value of sin(i) for that system (where i is the orbital inclination angle). Thus the value of sin(i) for a system with any particular observed value of M_T sin(i) cannot be assumed to be drawn randomly from a distribution corresponding to an isotropic i distribution, i.e. the presumptive prior distribution . Rather, the posterior distribution from which it is drawn depends on the intrinsic distribution of M_T for the exoplanet population being studied. We give a simple Bayesian derivation of this relationship and apply it to several "toy models" for the (currently unknown) intrinsic distribution of M_T. The results show that the effect can be an important one. For example, even for simple power-law distributions of M_T, the median value of sin(i) in an observed RV sample can vary between 0.860 and 0.023 (as compared to the 0.866 value for an isotropic i distribution) for indices (alpha) of the power-law in the range between -2 and +1, respectively. Over the same range of indicies, the 95% confidence interval on M_T varies from 1.002-4.566 (alpha = -2) to 1.13-94.34 (alpha = +1) times larger than M_T sin(i) due to sin(i) uncertainty alone. Our qualitative conclusion is that RV studies of exoplanets, both individual objects and statistical samples, should regard the sin(i) factor as more than a "numerical constant of order unity" with simple and well understood statistical properties. We argue that reports of M_T sin(i) determinations should be accompanied by a statement of the corresponding confidence bounds on M_T at, say, the 95% level based on an explicitly stated assumed form of the true M_T distribution in order to more accurately reflect the mass uncertainties associated with RV studies.

It is suggested an expanding locally anisotropic metric (ELA) ansatz describing matter in a flat expanding universe which interpolates between the Schwarzschild (SC) metric near point-like central bodies of mass ‘M’ and the Robertson-Walker (RW) metric for large radial coordinate: ‘ds^2=Z(cdt)2 – 1/Z (dr1-(Hr1/c) Z^(alpha/2+1/2)(cdt))^2-r1^2 dOmega’, where ‘Z=1-U’ with ‘U=2GM/(c^2r1)’, ‘G’ is the Newton constant, ‘c’ is the speed of light, ‘H=H(t)=\dot(a)/a’ is the time-dependent Hubble rate, ‘dOmega=dtheta^2+sin^2(theta) dvarphi^2′ is the solid angle element, ‘a’ is the universe scale factor and we are employing the coordinates ‘r1=ar’, being ‘r’ the radial coordinate for which the RW metric is diagonal. For constant exponent ‘alpha=alpha0=0′ it is retrieved the isotropic McVittie (McV) metric and for ‘alpha=alpha0=1′ it is retrieved the locally anisotropic Cosmological-Schwarzschild (SCS) metric, both already discussed in the literature. However it is shown that only for constant exponent ‘alpha=alpha0> 1′ exists an event horizon at the SC radius ‘r1=2GM/c^2′ and only for ‘alpha=alpha0>= 3′ space-time is singularity free for this value of the radius. These bounds exclude the previous existing metrics, for which the SC radius is a naked extended singularity. In addition it is shown that for ‘alpha=alpha0>5′ space-time is approximately Ricci flat in a neighborhood of the event horizon such that the SC metric is a good approximation in this neighborhood. It is further shown that to strictly maintain the SC mass pole at the origin ‘r1=0′ without the presence of more severe singularities it is required a radial coordinate dependent correction to the exponent ‘alpha(r1)=alpha0+alpha1 ‘2GM/(c^2 r1)’ with a negative coefficient ‘alpha1<0′. The energy-momentum density, pressures and equation of state are discussed.

The ACS and NICMOS have been used to obtain new HST images of NGC 4038/4039 ("The Antennae"). These new observations allow us to better differentiate compact star clusters from individual stars, based on both size and color. We use this ability to extend the cluster luminosity function by approximately two magnitudes over our previous WFPC2 results, and find that it continues as a single power law, dN/dL propto L^alpha with alpha=-2.13+/-0.07, down to the observational limit of Mv~-7. Similarly, the mass function is a single power law dN/dM propto M^beta with beta=-2.10+/-0.20 for clusters with ages t<3×10^8 yr, corresponding to lower mass limits that range from 10^4 to 10^5 Msun, depending on the age range of the subsample. Hence the power law indices for the luminosity and mass functions are essentially the same. The luminosity function for intermediate-age clusters (i.e., ~100-300 Myr old objects found in the loops, tails, and outer areas) shows no bend or turnover down to Mv~-6, consistent with relaxation-driven cluster disruption models which predict the turnover should not be observed until Mv~-4. An analysis of individual ~0.5-kpc sized areas over diverse environments shows good agreement between values of alpha and beta, similar to the results for the total population of clusters in the system. Several of the areas studied show evidence for age gradients, with somewhat older clusters appearing to have triggered the formation of younger clusters. The area around Knot B is a particularly interesting example, with an ~10-50 Myr old cluster of estimated mass ~10^6 Msun having apparently triggered the formation of several younger, more massive (up to 5×10^6 Msun) clusters along a dust lane.

(Abridged) A procedure is suggested to explore the value of F = alpha^2/mu, where mu = m_e/m_p is the electron-to-proton mass ratio, and alpha is the fine-structure constant. The fundamental physical constants, which are measured in different physical environments of high (terrestrial) and low (interstellar) densities of baryonic matter are supposed to vary in chameleon-like scalar field models, which predict that both masses and coupling constant may depend on the local matter density. The parameter Delta F/F = (F_obs – F_lab)/F_lab can be estimated from the radial velocity offset, Delta V = V_rot-V_fs, between the low-laying rotational transitions in carbon monoxide 13CO and the fine-structure transitions in atomic carbon [CI]. A model-dependent constraint on Delta alpha/alpha can be obtained from Delta F/F using Delta mu/mu independently measured from the ammonia method. Currently available radio astronomical datasets provide an upper limit on |Delta V| < 110 m/s (1sigma). When interpreted in terms of the spatial variation of F, this gives |Delta F/F| < 3.7*10^-{7}. An order of magnitude improvement of this limit will allow us to test independently a non-zero value of Delta mu/mu = (2.2 +/- 0.4_stat +/- 0.3_sys)*10^{-8} recently found with the ammonia method. Taking into account that the ammonia method restricts the spatial variation of mu at the level of |Delta mu/mu| <= 3*10^{-8} and assuming that Delta F/F is the same in the entire interstellar medium, one obtains that the spatial variation of alpha does not exceed the value |Delta alpha/alpha| < 2*10^{-7}. Since extragalactic gas clouds have densities similar to those in the interstellar medium, the bound on Delta alpha/alpha is also expected to be less than 2*10^{-7} at high redshift if no significant temporal dependence of alpha is present.

The existence of magnetohydrodynamic mean-field alpha^2-dynamos with spherically symmetric, isotropic helical turbulence function alpha is related to a non-self-adjoint spectral problem for a coupled system of two singular second order ordinary differential equations. We establish global estimates for the eigenvalues of this system in terms of the turbulence function alpha and its derivative alpha’. They allow us to formulate an anti-dynamo theorem and a non-oscillation theorem. The conditions of these theorems, which again involve alpha and alpha’, must be violated in order to reach supercritical or oscillatory regimes.

In 2008 the blazar Markarian 421 entered a very active phase and was one of the brightest sources in the sky at TeV energies, showing frequent flaring episodes. Using the data of ARGO-YBJ, a full coverage air shower detector located at Yangbajing (4300 m a.s.l., Tibet, China), we monitored the source at gamma ray energies E > 0.3 TeV during the whole year. The observed flux was variable, with the strongest flares in March and June, in correlation with X-ray enhanced activity. While during specific episodes the TeV flux could be several times larger than the Crab Nebula one, the average emission from day 41 to 180 was almost twice the Crab level, with an integral flux of (3.6 +-0.6) 10^-11 photons cm^-2 s^-1 for energies E > 1 TeV, and decreased afterwards. This paper concentrates on the flares occurred in the first half of June. This period has been deeply studied from optical to 100 MeV gamma rays, and partially up to TeV energies, since the moonlight hampered the Cherenkov telescope observations during the most intense part of the emission. Our data complete these observations, with the detection of a signal with a statistical significance of 3.8 standard deviations on June 11-13, corresponding to a gamma ray flux about 6 times larger than the Crab one above 1 TeV. The reconstructed differential spectrum, corrected for the intergalactic absorption, can be represented by a power law with an index alpha = -2.1 extending up to several TeV. The spectrum slope is fully consistent with previous observations reporting a correlation between the flux and the spectral index, suggesting that this property is maintained in different epochs and characterizes the source emission processes.

The optical spectra of objects classified as QSOs in the SDSS DR6 are analyzed with the aim of determining the value of the fine structure constant in the past and then check for possible changes in the constant over cosmological timescales. The analysis is done by measuring the position of the fine structure lines of the [OIII] doublet (4959 and 5008) in QSO nebular emission. From the sample of QSOs at redshifts z < 0.8 a subsample was selected on the basis of the amplitude and width of the [OIII] lines. Two different method were used to determine the position of the lines of the [OIII] doublet, both giving similar results. Using a clean sample containing 1568 of such spectra, a value of Delta alpha /alpha=(+2.4 +-2.5) x 10^{-5} (in the range of redshifts 0-0.8) was determined. The use of a larger number of spectra allows a factor ~5 improvement on previous constraints based on the same method. On the whole, we find no evidence of changes in alpha on such cosmological timescales. The mean variation compatible with our results is 1/ <t> Delta alpha/alpha=(+0.7 +- 0.7) x 10^{-14} yr^{-1}. The analysis was extended to the [NeIII] and [SII] doublets, although their usefulness is limited due to the fact that all these doublets in QSOs tend to be fainter than [OIII], and that some of them are affected by systematics.

We report the detection of 21-cm and H2 absorption lines in the same DLA system (log N(HI)=21.36+-0.10) at zabs=3.17447 towards SDSSJ133724+315254 (z=3.174). We estimate the spin temperature of the gas to be, Ts~600 K, intermediate between the expected values for cold and warm neutral media. This suggests that the HI absorption originates from a mixture of different phases. The total molecular fraction is low, f=10^-7, and H2 rotational level populations are not in equilibrium. The average abundance of the alpha-elements is, [S/H]=-1.45. N and Fe are found underabundant with respect to alpha-elements by ~1.0 dex and ~0.5 dex respectively. Using photoionization models we conclude that the gas is located more than 270 kpc away from the QSO. While the position of 21-cm absorption line coincides with the H2 velocity profile, their centroid are shifted by 2.7+-1.0 km/s from each other. However, the position of the strongest metal absorption component matches the position of the 21-cm absorption line within 0.5 km/s. From this, we constrain the variation of the combination of fundamental constants x=alpha^2 Gp/mu, Delta x/x=-(1.7+-1.7)x10^-6. This system is unique as we can at the same time have an independent constrain on mu using H2 lines. However only Werner band absorption lines are seen and the range of sensitivity coefficients is too narrow to provide a stringent constraint: Delta mu/mu <= 4.0×10^-4. The VLT/UVES spectrum reveals another DLA at zabs=3.16768 with log N(HI)=20.41+-0.15 and low metallicity, [Si/H]=-2.68+-0.11. We derive log N(DI)/N(HI)=-(4.93+-0.15) in this system. This is a factor of two smaller than the value expected from the best fitted value of Omega_b from the WMAP 5 yr data. This confirms the presence of astration of deuterium even at very low metallicity. [abridged]

The angular power spectrum of HI 21 cm opacity fluctuations is a useful statistic for quantifying the observed opacity fluctuations as well as for comparing these with theoretical models. We present here the HI 21 cm opacity fluctuation power spectrum towards the supernova remnant Cas A from interferometric data with spacial resolution of 5" and spectral resolution of 0.4 km/s. The power spectrum has been estimated using a simple but robust visibility based technique. We find that the power spectrum is well fit by a power law P_tau(U) = U^{alpha} with a power law index of alpha ~ -2.86 +/- 0.10 (3 sigma error) over the scales of 0.07 – 2.3 pc for the gas in the Perseus spiral arm and 0.002 – 0.07 pc (480 – 15730 au) for that in the Local arm. This estimated power law index is consistent with earlier observational results based on both HI emission over larger scales and absorption studies over a similar range of scales. We do not detect any statistically significant change in the power law index with the velocity width of the frequency channels. This constrains the power law index of the velocity structure function to be beta = 0.2 +/- 0.6 (3 sigma error).

The angular power spectrum of HI 21 cm opacity fluctuations is a useful statistic for quantifying the observed opacity fluctuations as well as for comparing these with theoretical models. We present here the HI 21 cm opacity fluctuation power spectrum towards the supernova remnant Cas A from interferometric data with spacial resolution of 5" and spectral resolution of 0.4 km/s. The power spectrum has been estimated using a simple but robust visibility based technique. We find that the power spectrum is well fit by a power law P_tau(U) = U^{alpha} with a power law index of alpha ~ -2.86 +/- 0.10 (3 sigma error) over the scales of 0.07 – 2.3 pc for the gas in the Perseus spiral arm and 0.002 – 0.07 pc (480 – 15730 au) for that in the Local arm. This estimated power law index is consistent with earlier observational results based on both HI emission over larger scales and absorption studies over a similar range of scales. We do not detect any statistically significant change in the power law index with the velocity width of the frequency channels. This constrains the power law index of the velocity structure function to be beta = 0.2 +/- 0.6 (3 sigma error).

Strong organized magnetic fields have been studied in the upper main sequence chemically peculiar stars for more than half a century. However, only recently have observational methods and numerical techniques become sufficiently mature to allow us to record and interpret high-resolution four Stokes parameter spectra, leading to the first assumption-free magnetic field models of these stars. Here we present a detailed magnetic Doppler imaging analysis of the spectropolarimetric observations of the prototypical magnetic Ap star alpha^2 CVn. The surface abundance distributions of Fe and Cr and a full vector map of the stellar magnetic field are reconstructed in a self-consistent inversion using our state-of-the-art magnetic Doppler imaging code Invers10. We succeeded in reproducing most of the details of the available spectropolarimetric observations of alpha^2 CVn with a magnetic map which combines a global dipolar-like field topology with localized spots of higher field intensity. We demonstrate that these small-scale magnetic structures are inevitably required to fit the linear polarization spectra; however, their presence cannot be inferred from the Stokes I and V observations alone. Our magnetic Doppler imaging analysis of alpha^2 CVn and previous results for 53 Cam support the view that the upper main sequence stars can harbour fairly complex surface magnetic fields which resemble oblique dipoles only at the largest spatial scales. Spectra in all four Stokes parameters are absolutely essential to unveil and meaningfully characterize this field complexity in Ap stars. We therefore suggest that understanding magnetism of stars in other parts of the H-R diagram is similarly incomplete without investigation of their linear polarization spectra.

We present a survey, using the Chandra X-ray observatory, of the central gravitating mass profiles in a sample of 10 galaxies, groups and clusters, spanning ~2 orders of magnitude in virial mass. We find the total mass distributions from ~0.2–10Re, where Re is the optical effective radius of the central galaxy, are remarkably similar to powerlaw density profiles. The negative logarithmic slope of the mass density profiles, alpha, systematically varies with Re, from alpha=2, for systems with Re~4kpc to alpha=1.2 for systems with Re>30kpc. Departures from hydrostatic equilibrium are likely to be small and cannot easily explain this trend. We show that the conspiracy between the baryonic (Sersic) and dark matter (NFW/ Einasto) components required to maintain a powerlaw total mass distribution naturally predicts an anti-correlation between alpha and Re that is very close to what is observed. The systematic variation of alpha with Re implies a dark matter fraction within Re that varies systematically with the properties of the galaxy in such a manner as to reproduce, without fine tuning, the observed tilt of the fundamental plane. We speculate that establishing a nearly powerlaw total mass distribution is therefore a fundamental feature of galaxy formation and the primary factor which determines the tilt of the fundamental plane.

We determine the slope of the near infrared extinction power law (A$_{\lambda} \propto \lambda^{-\alpha}$) for 8 regions of the Galaxy between l$\sim27^{\circ}$ and $\sim100^{\circ}$. UKIDSS Galactic Plane Survey data are compared, in colour-colour space, with Galactic population synthesis model data reddened using a series of power laws and convolved through the UKIDSS filter profiles. Monte Carlo simulations allow us to determine the best fit value of $\alpha$ and evaluate the uncertainty. All values are consistent with each other giving an average extinction power law of $\alpha$=2.14$^{+0.04}_{-0.05}$. This is much steeper than most laws previously derived in the literature from colour excess ratios, which are typically between 1.6 and 1.8. We show that this discrepancy is due to an inappropriate choice of filter wavelength in conversion from colour excess ratios to $\alpha$ and that effective rather than isophotal wavelengths are more appropriate. In addition, curved reddening tracks, which depend on spectral type and filter system, should be used instead of straight vectors.

We determine the slope of the near infrared extinction power law (A$_{\lambda} \propto \lambda^{-\alpha}$) for 8 regions of the Galaxy between l$\sim27^{\circ}$ and $\sim100^{\circ}$. UKIDSS Galactic Plane Survey data are compared, in colour-colour space, with Galactic population synthesis model data reddened using a series of power laws and convolved through the UKIDSS filter profiles. Monte Carlo simulations allow us to determine the best fit value of $\alpha$ and evaluate the uncertainty. All values are consistent with each other giving an average extinction power law of $\alpha$=2.14$^{+0.04}_{-0.05}$. This is much steeper than most laws previously derived in the literature from colour excess ratios, which are typically between 1.6 and 1.8. We show that this discrepancy is due to an inappropriate choice of filter wavelength in conversion from colour excess ratios to $\alpha$ and that effective rather than isophotal wavelengths are more appropriate. In addition, curved reddening tracks, which depend on spectral type and filter system, should be used instead of straight vectors.

(Note: this is a shortened version of the original "structured" A&A format abstract.) We performed a large optical spectroscopic and photometric survey of the Lynds~1630N and 1641 clouds. We provide a catalog of 132 confirmed young stars in L1630N and 267 such objects in L1641. We identify 28 transition disk systems, 20 of which were previously unknown, as well as 42 new transition disk candidates for which we have broad-band photometry but no optical spectroscopy. We estimate mass accretion rates M_acc from the equivalent widths of the H_alpha, H_beta, and HeI 5876\AA emission lines, and find a dependence on stellar mass of M_acc propto Mstar^alpha, with alpha~3.1 in the subsolar mass range that we probe. An investigation of a large literature sample of mass accretion rate estimates yields a similar slope of alpha~2.8 in the subsolar regime, but a shallower slope of alpha~2.0 if the whole mass range of 0.04 M_sun-5 Msun is included. Among the transition disk objects, the fraction of stars that show significant accretion activity is relatively low compared to stars with still optically thick disks (26\pm11% vs. 57\pm6%, respectively). However, those transition disks that do show significant accretion have the same median accretion rate as normal optically thick disks of 3-4*10^{-9} M_sun/yr. We find that the ages of the transition disks and the WTTSs without disks are statistically indistinguishable, and both groups are significantly older than the CTTSs. These results argue against disk-binary interaction or gravitational instability as mechanisms causing a transition disk appearance. Our observations indicate that disk lifetimes in the clustered population are shorter than in the distributed population. We propose refined Halpha equivalent width criteria to distinguish WTTSs from CTTSs.

The Tensor-Vector-Scalar theory of gravity, which was designed as a relativistic implementation to the modified dynamics paradigm, has fared quite well as an alternative to dark matter, on both galactic and cosmological scales. However, its performance in the solar system, as embodied in the post-Newtonian formalism, has not yet been fully investigated. Tamaki has recently attempted to calculate the preferred frame parameters for TeVeS, but ignored the cosmological value of the scalar field, thus concluding that the Newtonian potential must be static in order to be consistent with the vector equation. We show that when the cosmological value of the scalar field is taken into account, there is no constraint on the Newtonian potential; however, the cosmological value of the scalar field is tightly linked to the vector field coupling constant K, preventing the former from evolving as predicted by its equation of motion. We then proceed to investigate the post-Newtonian limit of a generalized version of TeVeS, with {\AE}ther type vector action, and show that its \beta,\gamma and \xi parameters are as in GR, while solar system constraints on the preferred frame parameters \alpha_1 and \alpha_2 can be satisfied within a modest range of small values of the scalar and vector fields coupling parameters, and for values of the cosmological scalar field consistent with evolution within the framework of existing models.

We summarize the attempts by our group and others to derive constraints on variations of fundamental constants over cosmic time using quasar absorption lines. Most upper limits reside in the range 0.5-1.5×10-5 at the 3sigma level over a redshift range of approximately 0.5-2.5 for the fine-structure constant, alpha, the proton-to-electron mass ratio, mu, and a combination of the proton gyromagnetic factor and the two previous constants, gp(alpha^2/mu)^nu, for only one claimed variation of alpha. It is therefore very important to perform new measurements to improve the sensitivity of the numerous methods to at least <0.1×10-5 which should be possible in the next few years. Future instrumentations on ELTs in the optical and/or ALMA, EVLA and SKA pathfinders in the radio will undoutedly boost this field by allowing to reach much better signal-to-noise ratios at higher spectral resolution and to perform measurements on molecules in the ISM of high redshift galaxies.

If a dynamical system is long-lived and non-resonant (that is, if there is a set of tracers that have evolved independently through many orbital times), and if the system is observed at any non-special time, it is possible to infer the dynamical properties of the system (such as the gravitational force or acceleration law) from a snapshot of the positions and velocities of the tracer population at a single moment in time. In this paper we describe a general inference technique that solves this problem while allowing (1) the unknown distribution function of the tracer population to be simultaneously inferred and marginalized over, and (2) prior information about the gravitational field and distribution function to be taken into account. As an example, we consider the simplest problem of this kind: We infer the force law in the Solar System using only an instantaneous kinematic snapshot (valid at 2009 April 1.0) for the eight major planets. We consider purely radial acceleration laws of the form a_r = -A [r/r_0]^{-\alpha}, where r is the distance from the Sun. Using a probabilistic inference technique, we infer 1.989 < \alpha < 2.052 (95 percent interval), largely independent of any assumptions about the distribution of energies and eccentricities in the system beyond the assumption that the system is phase-mixed. Generalizations of the methods used here will permit, among other things, inference of Milky Way dynamics from Gaia-like observations.

If a dynamical system is long-lived and non-resonant (that is, if there is a set of tracers that have evolved independently through many orbital times), and if the system is observed at any non-special time, it is possible to infer the dynamical properties of the system (such as the gravitational force or acceleration law) from a snapshot of the positions and velocities of the tracer population at a single moment in time. In this paper we describe a general inference technique that solves this problem while allowing (1) the unknown distribution function of the tracer population to be simultaneously inferred and marginalized over, and (2) prior information about the gravitational field and distribution function to be taken into account. As an example, we consider the simplest problem of this kind: We infer the force law in the Solar System using only an instantaneous kinematic snapshot (valid at 2009 April 1.0) for the eight major planets. We consider purely radial acceleration laws of the form a_r = -A [r/r_0]^{-\alpha}, where r is the distance from the Sun. Using a probabilistic inference technique, we infer 1.989 < \alpha < 2.052 (95 percent interval), largely independent of any assumptions about the distribution of energies and eccentricities in the system beyond the assumption that the system is phase-mixed. Generalizations of the methods used here will permit, among other things, inference of Milky Way dynamics from Gaia-like observations.

We calculate heavy element enrichment in a Jupiter-mass protoplanet formed by disk instability at various radial distances from the star, considering different disk masses and surface density distributions. Although the available mass for accretion increases with radial distance (a) for disk solid surface density (sigma) functions sigma=sigma_0*a^(-alpha) with alpha < 2, the accretion timescale is significantly longer at larger radial distances. Efficient accretion is limited to the first ~ 1E5 years of planetary evolution, when the planet is extended and before gap opening and type II migration take place. The accreted mass is calculated for disk masses of 0.01, 0.05 and 0.1 M_sun with alpha = 1/2, 1, and 3/2. We show that a Jupiter-mass protoplanet can accrete 1 to 110 M_earth of heavy elements, depending on the disk properties. Our results explain the large variation in heavy element enrichment found in extra-solar giant planets. Since higher disk surface density is found to lead to larger heavy element enrichment, our model results are consistent with the correlation between heavy element enrichment and stellar metallicity. Our calculations also suggest that Jupiter could have formed at a larger radial distance than its current location while still accreting the mass of heavy elements predicted by interior models. We conclude that in the disk instability model the final composition of a giant planet is strongly determined by its formation environment. The heavy element abundance of a giant planet does not discriminate between its origin by either disk instability or core accretion.

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.

The study of young stellar populations has revealed that most stars are in binary or higher order multiple systems. In this study the influence on the stellar initial mass function (IMF) of large quantities of unresolved multiple massive stars is investigated by taking into account stellar evolution and photometrically determined system masses. The models where initial masses are derived from the luminosity and colour of unresolved multiple systems show that even under extreme circumstances (100% binaries or higher order multiples) the difference between the power-law index of the mass function of all stars and the observed mass function is small (~0.1). Thus, if the observed IMF has the Salpeter index alpha = 2.35 then the true stellar IMF has an index not flatter than alpha = 2.25. Additionally, unresolved multiple systems may hide between 15 and 60% of the underlying true mass of a star cluster. While already a known result, it is important to point out that the presence of a large number of unresolved binaries amongst pre-main-sequence (PMS) stars induces a significant spread in the measured ages of these stars even if there is none. Also, lower-mass stars in a single-age binary-rich cluster appear older than the massive stars by about 0.6 Myr.

We present new results on the frequency distribution of projected HI column densities f(N,X), total comoving covering fraction, and integrated mass densities rho_HI of high redshift, HI `disks’ from a survey of damped Lya systems (DLAs) in the Sloan Digital Sky Survey, Data Release 5. For the full sample spanning z=2.2 to 5 [738 DLAs], f(N,X) is well fitted by a double power-law with a break column density N_d = 10^(21.55 +/- 0.04) and low/high-end exponents alpha = -2.00 +/- 0.05, -6.4^{+1.1}_{-1.6}. The shape of f(N,X) is invariant during this redshift interval and also follows the projected surface density distribution of present-day HI disks as inferred from 21cm observations. We conclude that HI gas has been distributed in a self-similar fashion for the past 12Gyr. The normalization of f(N,X), in contrast, decreases by a factor of two during the ~2Gyr interval from z=4 to 2.2 giving corresponding decreases in both the total covering fraction and rho_HI. At z~2, these quantities match the present-day values suggesting no evolution during the past ~10Gyr. We argue that the evolution at early times is driven by `violent’ processes that removes gas from nearly half the galaxies at z~3 establishing the antecedants of current early-type galaxies. The perceived constancy of rho_HI, meanwhile, implies that HI gas is a necessary but insufficient pre-condition for star formation and that the global star-formation rate is driven by the accretion and condensation of fresh gas from the intergalactic medium.