Residual noise covariance for Planck low-resolution data analysis

Aims: Develop and validate tools to estimate residual noise covariance in Planck frequency maps. Quantify signal error effects and compare different techniques to produce low-resolution maps. Methods: We derive analytical estimates of covariance of the residual noise contained in low-resolution maps produced using a number of map-making approaches. We test these analytical predictions using Monte Carlo simulations and their impact on angular power spectrum estimation. We use simulations to quantify the level of signal errors incurred in different resolution downgrading schemes considered in this work. Results: We find an excellent agreement between the optimal residual noise covariance matrices and Monte Carlo noise maps. For destriping map-makers, the extent of agreement is dictated by the knee frequency of the correlated noise component and the chosen baseline offset length. The significance of signal striping is shown to be insignificant when properly dealt with. In map resolution downgrading, we find that a carefully selected window function is required to reduce aliasing to the sub-percent level at multipoles, ell > 2Nside, where Nside is the HEALPix resolution parameter. We show that sufficient characterization of the residual noise is unavoidable if one is to draw reliable contraints on large scale anisotropy. Conclusions: We have described how to compute the low-resolution maps, with a controlled sky signal level, and a reliable estimate of covariance of the residual noise. We have also presented a method to smooth the residual noise covariance matrices to describe the noise correlations in smoothed, bandwidth limited maps.

Understanding Astrophysical Noise from Stellar Surface Magneto-Convection

To obtain cm/s precision, stellar surface magneto-convection must be disentangled from observed radial velocities (RVs). In order to understand and remove the convective signature, we create Sun-as-a-star model observations based on a 3D magnetohydrodynamic solar simulation. From these Sun-as-a-star model observations, we find several line characteristics are correlated with the induced RV shifts. The aim of this campaign is to feed directly into future high precision RV studies, such as the search for habitable, rocky worlds, with forthcoming spectrographs such as ESPRESSO.

Stout: Cloudy's Atomic and Molecular Database

We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures and densities spanning the low-to high-density limits. The radiation field between photon energies 10−8 Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~102 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon.

LOTUS: A low cost, ultraviolet spectrograph

We describe the design, construction and commissioning of LOTUS; a simple, low-cost long-slit spectrograph for the Liverpool Telescope. The design is optimized for near-UV and visible wavelengths and uses all transmitting optics. It exploits the instrument focal plane field curvature to partially correct axial chromatic aberration. A stepped slit provides narrow (2.5x95 arcsec) and wide (5x25 arcsec) options that are optimized for spectral resolution and flux calibration respectively. On sky testing shows a wavelength range of 3200-6300 Angstroms with a peak system throughput (including detector quantum efficiency) of 15 per cent and wavelength dependant spectral resolution of R=225-430. By repeated observations of the symbiotic emission line star AG Peg we demonstrate the wavelength stability of the system is less than 2 Angstroms rms and is limited by the positioning of the object in the slit. The spectrograph is now in routine operation monitoring the activity of comet 67P/Churyumov-Gerasimenko during its current post-perihelion apparition.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

Impact of cosmic inhomogeneities on SNe observations

We study the impact of cosmic inhomogeneities on the interpretation of SNe observations. We build an inhomogeneous universe model that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on the data, allow the model to mimic the concordance model.

On global solar dynamo simulations

Global dynamo simulations solving the equations of magnetohydrodynamics (MHD) have been a tool of astrophysicists who try to understand the magnetism of the Sun for several decades now. During recent years many fundamental issues in dynamo theory have been studied in detail by means of local numerical simulations that simplify the problem and allow the study of physical effects in isolation. Global simulations, however, continue to suffer from the age-old problem of too low spatial resolution, leading to much lower Reynolds numbers and scale separation than in the Sun. Reproducing the internal rotation of the Sun, which plays a crucual role in the dynamo process, has also turned out to be a very difficult problem. In the present paper the current status of global dynamo simulations of the Sun is reviewed. Emphasis is put on efforts to understand how the large-scale magnetic fields, i.e. whose length scale is greater than the scale of turbulence, are generated in the Sun. Some lessons from mean-field theory and local simulations are reviewed and their possible implications to the global models are discussed. Possible remedies to some of the current issues of the solar simulations are put forward.

Magnetorotational instability driven dynamos at low magnetic Prandtl numbers

Numerical simulations of the magnetorotational instability (MRI) with zero initial net flux in a non-stratified isothermal cubic domain are used to demonstrate the importance of magnetic boundary conditions. In fully periodic systems the level of turbulence generated by the MRI strongly decreases as the magnetic Prandtl number (Pm), which is the ratio of kinematic viscosity and magnetic diffusion, is decreased. No MRI or dynamo action below Pm=1 is found, agreeing with earlier investigations. Using vertical field conditions, which allow magnetic helicity fluxes out of the system, the MRI is found to be excited in the range 0.1

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.