Search

The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system., 13 pages, 5 figures. Accepted by MNRAS

The black holes (BH) in merging BH-BH binaries are likely progeny of binary O stars. Their properties, including their spins, will be strongly influenced by the evolution of their progenitor O stars. The remarkable observation that many single O stars spin very rapidly can be explained if they accreted angular momentum from a mass-transferring, O-type or Wolf-Rayet companion before that star blew up as a supernova. To test this prediction, we have measured the spin rates of eight O stars in Wolf-Rayet (WR) + O binaries, increasing the total sample size of such O stars' measured spins from two to ten. Polarimetric and other determinations of these systems' sin i allow us to determine an average equatorial rotation velocity from HeI (HeII) lines of $v_e$ = 348 (173) km/s for these O stars, with individual star's $v_e$ from HeI (HeII) lines ranging from 482 (237) to 290 (91) km/s. We argue that the $\sim$ 100\% difference between HeI and HeII speeds is due to gravity darkening. Super-synchronous spins, now observed in all 10 O stars in WR + O binaries where it has been measured, are strong observational evidence that Roche lobe overflow mass transfer from a WR progenitor companion has played a critical role in the evolution of WR+OB binaries. While theory predicts that this mass transfer rapidly spins-up the O-type mass gainer to a nearly break-up rotational velocity $v_e \sim 530 $ km/s, the observed average $v_e$ of the O-type stars in our sample is 65\% that speed. This demonstrates that, even over the relatively short WR-phase timescale, tidal and/or other effects causing rotational spin-down must be efficient. A challenge to tidal synchronization theory is that the two longest-period binaries in our sample (with periods of 29.7 and 78.5 days) unexpectedly display super-synchronous rotation., 10 pages, 10 figures, submitted to MNRAS

We examined variations in the 160/250 and 250/350 micron surface brightness ratios within 24 nearby (=250 microns (and sometimes at >=160 microns) originates from dust heated by the older stellar population. In the rest of the sample, either the results are indeterminate or both of these stellar populations may contribute equally to the global dust heating. The observed variations in dust heating mechanisms does not necessarily match what has been predicted by dust emission and radiative transfer models, which could lead to overestimated dust temperatures, underestimated dust masses, false detections of variability in dust emissivity, and inaccurate star formation rate measurements., Accepted for publication in MNRAS

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy., 12 pages, 4 figures, accepted for publication in MNRAS

In this paper, we consider the implications of the distributions of dust and metals in the disc of M31. We derive mean radial dust distributions using a dust map created from Herschel images of M31 sampling the entire far-infrared peak. Modified blackbodies are fit to approximately 4000 pixels with a varying, as well as a fixed, dust emissivity index (beta). An overall metal distribution is also derived using data collected from the literature. We use a simple analytical model of the evolution of the dust in a galaxy with dust contributed by stellar sources and interstellar grain growth, and fit this model to the radial dust-to-metals distribution across the galaxy. Our analysis shows that the dust-to-gas gradient in M31 is steeper than the metallicity gradient, suggesting interstellar dust growth is (or has been) important in M31. We argue that M31 helps build a case for cosmic dust in galaxies being the result of substantial interstellar grain growth, while the net dust production from stars may be limited.We note, however, that the efficiency of dust production in stars, e.g. in supernovae ejecta and/or stellar atmospheres, and grain destruction in the interstellar medium may be degenerate in our simple model. We can conclude that interstellar grain growth by accretion is likely at least as important as stellar dust production channels in building the cosmic dust component in M31.

Evidence is mounting that a significant fraction of the early-type galaxy population contains substantial reservoirs of cold interstellar gas and dust. We investigate the gas and dust in NGC 5485, an early-type galaxy with a prominent minor-axis dust lane. Using new Herschel PACS and SPIRE imaging data, we detect 3.8 × 106?M? of cool interstellar dust in NGC 5485, which is in stark contrast with the non-detection of the galaxy in sensitive H?i and CO observations from the ATLAS3D consortium. The resulting gas-to-dust ratio upper limit is Mgas/Md

We investigate the amount and spatial distribution of interstellar dust in edge-on spiral galaxies, using detailed radiative transfer modelling of a homogeneous sample of 12 galaxies selected from the Calar Alto Legacy Integral Field Area survey. Our automated fitting routine, fitskirt, was first validated against artificial data. This is done by simultaneously reproducing the Sloan Digital Sky Survey g-, r-, i- and z-band observations of a toy model in order to combine the information present in the different bands. We show that this combined, oligochromatic fitting has clear advantages over standard monochromatic fitting especially regarding constraints on the dust properties. We model all galaxies in our sample using a three-component model, consisting of a double-exponential disc to describe the stellar and dust discs and using a Sérsic profile to describe the central bulge. The full model contains 19 free parameters, and we are able to constrain all these parameters to a satisfactory level of accuracy without human intervention or strong boundary conditions. Apart from two galaxies, the entire sample can be accurately reproduced by our model. We find that the dust disc is about 75 per?cent more extended but only half as high as the stellar disc. The average face-on optical depth in the V band is 0.76 and the spread of 0.60 within our sample is quite substantial, which indicates that some spiral galaxies are relatively opaque even when seen face-on.

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina and Galactocentric orbits (consistent with proper motions), and compare the simulation line-of-sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed losVDs of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11-300 eV: always assuming that ??s=0.225 for H0 = 72?km?s?1?Mpc?1. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalizations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc?h?1, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low-mass (M200 1015.1?M?) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ? 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark haloes that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (also known as Modified Newtonian Dynamics or MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host and the size of the hole at the centre of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

In recent years it has become clear that the vertical structure of disk galaxies is a key ingredient for understanding galaxy evolution. In particular, the presence and struc- ture of extra-planar gas has been a focus of research. The Hydrogen Accretion in LOcal GAlaxieS (HALOGAS) survey aims to provide a census on the rate of cold neutral gas accretion in nearby galaxies as well as a statistically significant set of galaxies that can be investigated for their extra-planar gas properties. In order to better understand the the vertical structure of the neutral hydrogen in the two edge-on HALOGAS galaxies NGC 5023 and UGC 2082 we construct de- tailed tilted ring models. The addition of distortions resembling arcs or spiral arms significantly improves the fit of the models to these galaxies. In the case of UGC 2082 no vertical gradient in rotational velocity is required in either symmetric models nor non-symmetric models to match the observations. The best fitting model features two arcs of large vertical extent that may be due to accretion. In the case of NGC 5023 a vertical gradient is required in symmetric models (dV/dz =?14.9±3.8 km s?1 kpc?1) and its magnitude is significantly lowered when non-symmetric models are considered (dV/dz =?9.4 ± 3.8 km s?1 kpc?1). Additionally it is shown that the underlying disk of NGC 5023 can be made symmetric, in all parameters except the warp, in non- symmetric models. In comparison to the "classical" modeling these models fit the data significantly better with a limited addition of free parameters.