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.

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.

Dependence of the large-scale vortex instability on latitude, stratification and domain size

In an earlier study, we reported on the excitation of large-scale vortices in Cartesian hydrodynamical convection models subject to rapid enough rotation. In that study, the conditions for the onset of the instability were investigated in terms of the Reynolds (Re) and Coriolis (Co) numbers in models located at the stellar North pole. In this study, we extend our investigation to varying domain sizes, increasing stratification, and place the box at different latitudes. The effect of the increasing box size is to increase the sizes of the generated structures, so that the principal vortex always fills roughly half of the computational domain. The instability becomes stronger in the sense that the temperature anomaly and change in the radial velocity are observed to be enhanced. The model with the smallest box size is found to be stable against the instability, suggesting that a sufficient scale separation between the convective eddies and the scale of the domain is required for the instability to work. The instability can be seen upto the colatitude of 30 degrees, above which value the flow becomes dominated by other types of mean flows. The instability can also be seen in a model with larger stratification. Unlike the weakly stratified cases, the temperature anomaly caused by the vortex structures is seen to depend on depth.

Lucky Imaging of transiting planet host stars with LuckyCam

We obtained high-resolution, high-contrast optical imaging in the Sloan Digital Sky Survey i′ band with the LuckyCam camera mounted on the 2.56 m Nordic Optical Telescope, to search for faint stellar companions to 16 stars harbouring transiting exoplanets. The Lucky imaging technique uses very short exposures to obtain near diffraction-limited images yielding sub-arcsecond sensitivity, allowing us to search for faint stellar companions within the seeing disc of the primary planet host. Here, we report the detection of two candidate stellar companions to the planet host TrES-1 at separations <6.5 arcsec and we confirm stellar companions to CoRoT-2, CoRoT-3, TrES-2, TrES-4 and HAT-P-7 already known in the literature. We do not confirm the candidate companions to HAT-P-8 found via Lucky imaging by Bergfors et al., however, most probably because HAT-P-8 was observed in poor seeing conditions. Our detection sensitivity limits allow us to place constraints on the spectral types and masses of the putative bound companions to the planet host stars in our sample. If bound, the stellar companions identified in this work would provide stringent observational constraints to models of planet formation and evolution. In addition, these companions could affect the derived physical properties of the exoplanets in these systems.

CI Aql: a Type Ia supernova progenitor?

If recurrent novae are progenitors of Type Ia supernovae, their white dwarfs must have masses close to the Chandrasekhar limit. The most reliable means of determining white dwarf masses in recurrent novae is dynamically, via radial-velocity and rotational-broadening measurements of the companion star. Such measurements require the system to be both eclipsing and to show absorption features from the secondary star. Prior to the work reported here, the only dynamical mass estimate of a recurrent nova was for U Sco, which has a white dwarf mass of 1.55 +/- 0.24 Msolar (Thoroughgood et al. 2001). We present new time-resolved, intermediate-resolution spectroscopy of the eclipsing recurrent nova CI Aquilae (CI Aql) during quiescence. We find the mass of the white dwarf to be 1.00 +/- 0.14 Msolar and the mass of the secondary star to be 2.32 +/- 0.19 Msolar. We estimate the radius of the secondary to be 2.07 +/- 0.06 Rsolar, implying that it is a slightly-evolved early A-type star. The high mass ratio of q = 2.35 +/- 0.24 and the high secondary-star mass implies that the mass transfer occurs on a thermal timescale. We suggest that CI Aql is rapidly evolving into a supersoft X-ray source, and ultimately may explode as a Type Ia supernova within 10 Myr.

Discovery of WASP-65b and WASP-75b: Two hot Jupiters without highly inflated radii

We report the discovery of two transiting hot Jupiters, WASP-65b (M_pl = 1.55 ± 0.16 M_J; R_pl = 1.11 ± 0.06 R_J), and WASP-75b (M_pl = 1.07 ± 0.05 M_J; R_pl = 1.27 ± 0.05 R_J). They orbit their host star every ∼2.311, and ∼2.484 days, respectively. The planet host WASP-65 is a G6 star (T_eff = 5600 K, [Fe/H] = −0.07 ± 0.07, age 8 Gyr); WASP-75 is an F9 star (T_eff = 6100 K, [Fe/H] = 0.07 ± 0.09, age ∼ 3 Gyr). WASP-65b is one of the densest known exoplanets in the mass range 0.1 and 2.0 MJ (ρpl = 1.13 ± 0.08 ρJ), a mass range where a large fraction of planets are found to be inflated with respect to theoretical planet models. WASP-65b is one of only a handful of planets with masses of ∼1.5 MJ, a mass regime surprisingly underrepresented among the currently known hot Jupiters. The radius of WASP-75b is slightly inflated (10%) as compared to theoretical planet models with no core, and has a density similar to that of Saturn (ρpl = 0.52 ± 0.06 ρJ). Key words. planetary systems – stars: individual:

WASP-54b, WASP-56b and WASP-57b: Three new sub-Jupiter mass planets from SuperWASP

We present three newly discovered sub-Jupiter mass planets from the SuperWASP survey: WASP-54b is a heavily bloated planet of mass 0.636$^{+0.025}_{-0.024}$ \mj and radius 1.653$^{+0.090}_{-0.083}$ \rj. It orbits a F9 star, evolving off the main sequence, every 3.69 days. Our MCMC fit of the system yields a slightly eccentric orbit ($e=0.067^{+0.033}_{-0.025}$) for WASP-54b. We investigated further the veracity of our detection of the eccentric orbit for WASP-54b, and we find that it could be real. However, given the brightness of WASP-54 V=10.42 magnitudes, we encourage observations of a secondary eclipse to draw robust conclusions on both the orbital eccentricity and the thermal structure of the planet. WASP-56b and WASP-57b have masses of 0.571$^{+0.034}_{-0.035}$ \mj and $0.672^{+0.049}_{-0.046}$ \mj, respectively; and radii of $1.092^{+0.035}_{-0.033}$ \rj for WASP-56b and $0.916^{+0.017}_{-0.014}$ \rj for WASP-57b. They orbit main sequence stars of spectral type G6 every 4.67 and 2.84 days, respectively. WASP-56b and WASP-57b show no radius anomaly and a high density possibly implying a large core of heavy elements; possibly as high as $\sim$50 M$_{\oplus}$ in the case of WASP-57b. However, the composition of the deep interior of exoplanets remain still undetermined. Thus, more exoplanet discoveries such as the ones presented in this paper, are needed to understand and constrain giant planets' physical properties.

A spectroscopically normal type Ic supernova from a very massive progenitor

We present observations of the Type Ic supernova (SN Ic) 2011bm spanning a period of about one year. The data establish that SN 2011bm is a spectroscopically normal SN Ic with moderately low ejecta velocities and with a very slow spectroscopic and photometric evolution (more than twice as slow as SN 1998bw). The Pan-STARRS1 retrospective detection shows that the rise time from explosion to peak was 40 days in the R band. Through an analysis of the light curve and the spectral sequence, we estimate a kinetic energy of 7-17 foe and a total ejected mass of 7-17 Mo, 5-10 Mo of which is oxygen and 0.6-0.7 Mo is 56Ni. The physical parameters obtained for SN 2011bm suggest that its progenitor was a massive star of initial mass 30-50 Mo. The profile of the forbidden oxygen lines in the nebular spectra show no evidence of a bi-polar geometry in the ejected material.