Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

Simulating the formation of a protocluster at z∼ 2

We present results from two high-resolution hydrodynamical simulations of protocluster regions at   z ≃ 2.1  . The simulations have been compared to observational results for the so-called Spiderweb galaxy system, the core of a putative protocluster region at   z = 2.16  , found around a radio galaxy. The simulated regions have been chosen so as to form a poor cluster with   M ₂₀₀≃ 10¹⁴  h ⁻¹ M_⊙  (C1) and a rich cluster with   M ₂₀₀≃ 2 × 10¹⁵  h ⁻¹ M_⊙  (C2) at   z = 0  . The simulated protoclusters show evidence of ongoing assembly of a dominating central galaxy. The stellar mass of the brightest cluster galaxy of the C2 system is in excess with respect to observational estimates for the Spiderweb galaxy, with a total star formation rate which is also larger than indicated by observations. We find that the projected velocities of galaxies in the C2 cluster are consistent with observations, while those measured for the poorer cluster C1 are too low compared with the observed velocities. We argue that the Spiderweb complex resembles the high-redshift progenitor of a rich galaxy cluster. Our results indicate that the included supernovae feedback is not enough to suppress star formation in these systems, supporting the need of introducing active galactic nuclei feedback. According to our simulations, a diffuse atmosphere of hot gas in hydrostatic equilibrium should already be present at this redshift, and enriched at a level comparable to that of nearby galaxy clusters. The presence of this gas should be detectable with future deep X-ray observations.     

A circular statistical method for extracting rotation measures

We propose a new method for the extraction of rotation measures from spectral polarization data. The method is based on maximum likelihood analysis and takes into account the circular nature of the polarization data. The method is unbiased and statistically more efficient than the standard χ ² procedure. We also find that the method is computationally much more convenient than the standard χ ² procedure if the number of data points is very large. We find that for most sources the method gives results very close to the standard χ ² minimization procedure. We give results for all the cases for which the method gives significantly different results from χ ² minimization. We also make a λ ³ fit to the data in order to extract non-Faraday rotation behaviour for those sources for which a large number of data points are available.     

Hubble Space Telescope observations of SV Cam – II. First derivative light-curve modelling using phoenix and atlas model atmospheres

The variation of the specific intensity across the stellar disc is an essential input parameter in surface brightness reconstruction techniques such as Doppler imaging, where the relative intensity contributions of different surface elements are important in detecting star-spots. We use phoenix and atlas model atmospheres to model light curves derived from high precision (signal-to-noise ratio ≃ 5000) Hubble Space Telescope ( HST ) data of the eclipsing binary SV Cam (F9V+K4V), where the variation of specific intensity across the stellar disc will determine the contact points of the binary system light curve. For the first time, we use χ² comparison fits to the first derivative profiles to determine the best-fitting model atmosphere. We show the wavelength dependence of the limb darkening and that the first derivative profile is sensitive to the limb-darkening profile very close to the limb of the primary star. It is concluded that there is only a marginal difference (<1σ) between the χ² comparison fits of the two model atmospheres to the HST light curve at all wavelengths. The usefulness of the second derivative of the light curve for measuring the sharpness of the primary's limb is investigated, but we find that the data are too noisy to permit a quantitative analysis.     

A modified χ2-test for cosmic microwave background analyses

We present a new general procedure for determining a given set of quantities. To this end, we define a certain statistic, which we call 'modified  χ²' (χ²_M)  , because of its similarity to the standard  χ²  . The terms of this  χ²_M  are made up of the fluctuations of an unbiased estimator of some statistical quantities and certain weights. Only the diagonal terms of the covariance matrix appear explicitly in our statistic, while the full covariance matrix (and not its inverse) is included implicitly in the calculation of the weights. Choosing these weights, we may obtain, through minimizing  χ²_M  , the estimator that provides the minimum rms, either for those quantities or for the parameters on which these quantities depend. In this paper, we describe our method in the context of cosmic microwave background experiments, in order to obtain either the statistical properties of the maps or the cosmological parameters. The test here is constructed out of some estimator of the two-point correlation function at different angles. For the problem of one-parameter estimation, we show that our method has the same power as the maximum-likelihood method. We have also applied this method to Monte Carlo simulations of the COBE -DMR data, as well as to the actual 4-yr data, obtaining consistent results with previous analyses. We also provide a very good analytical approximation to the distribution function of our statistic, which could also be useful in other contexts.     

Gemini/GMOS spectra of globular clusters in the Virgo giant elliptical NGC 4649

NGC 4649 (M60) is one of a handful of giant Virgo ellipticals. We have obtained Gemini/GMOS (Gemini North Multi-Object Spectrograph) spectra for 38 globular clusters (GCs) associated with this galaxy. Applying the multi-index  χ²  minimization technique of Proctor and Sansom with the single stellar population models of Thomas, Maraston and Korn, we derive ages, metallicities and α-element abundance ratios. We find several young (2–3 Gyr old) supersolar metallicity GCs, while the majority are old (>10 Gyr), spanning a range of metallicities from solar to  [Z/H]=−2  . At least two of these young GCs are at large projected radii of 17–20 kpc. The galaxy itself shows no obvious signs of a recent starburst, interaction or merger. A trend of decreasing α-element ratio with increasing metallicity is found.     

Bayesian photometric redshifts with empirical training sets

We combine in a single framework the two complementary benefits of  χ²  template fits and empirical training sets used e.g. in neural nets:  χ²  is more reliable when its probability density functions (PDFs) are inspected for multiple peaks, while empirical training is more accurate when calibration and priors of query data and training set match. We present a  χ²  empirical method that derives PDFs from empirical models as a subclass of kernel regression methods, and apply it to the Sloan Digital Sky Survey Data Release 5 sample of >75 000 quasi-stellar objects, which is full of ambiguities. Objects with single-peak PDFs show <1 per cent outliers, rms redshift errors <0.05 and vanishing redshift bias. At   z > 2.5  , these figures are two times better. Outliers result purely from the discrete nature and limited size of the model, and rms errors are dominated by the intrinsic variety of object colours. PDFs classed as ambiguous provide accurate probabilities for alternative solutions and thus weights for using both solutions and avoiding needless outliers. E.g. the PDFs predict 78.0 per cent of the stronger peaks to be correct, which is true for 77.9 per cent of them. Redshift incompleteness is common in faint spectroscopic surveys and turns into a massive undetectable outlier risk above other performance limitations, but we can quantify residual outlier risks stemming from size and completeness of the model. We propose a matched  χ²  error scale for noisy data and show that it produces correct error estimates and redshift distributions accurate within Poisson errors. Our method can easily be applied to future large galaxy surveys, which will benefit from the reliability in ambiguity detection and residual risk quantification.     

Lyman break galaxies at z= 4–6 in cosmological smoothed particle hydrodynamics simulations

We perform a spectrophotometric analysis of galaxies at redshifts z = 4–6 in cosmological smoothed particle hydrodynamics simulations of a Λ cold dark matter universe. Our models include radiative cooling and heating by a uniform ultraviolet (UV) background, star formation, supernova feedback, and a phenomenological model for galactic winds. Analysing a series of simulations of varying box size and particle number allows us to isolate the impact of numerical resolution on our results. Specifically, we determine the luminosity functions in B , V , R , i ' and z ' filters, and compare the results with observational surveys of Lyman break galaxies (LBGs) performed with the Subaru telescope and the Hubble Space Telescope . We find that the simulated galaxies have UV colours consistent with observations and fall in the expected region of the colour–colour diagrams used by the Subaru group. The stellar masses of the most massive galaxies in our largest simulation increase their stellar mass from   M _★∼ 10¹¹ M_⊙  at z = 6 to   M _★∼ 10^11.7 M_⊙  at z = 3. Assuming a uniform extinction of E ( B − V ) = 0.15, we also find reasonable agreement between simulations and observations in the space density of UV bright galaxies at z = 3–6, down to the magnitude limit of each survey. For the same moderate extinction level of E ( B − V ) ∼ 0.15, the simulated luminosity functions match observational data, but have a steep faint-end slope with α∼−2.0. We discuss the implications of the steep faint-end slope found in the simulations. Our results confirm the generic conclusion from earlier numerical studies that UV bright LBGs at z ≥ 3 are the most massive galaxies with E ( B − V ) ∼ 0.15 at each epoch.     

Gas-phase mean opacities for varying [M/H], N/O and C/O

We present a set of gas-phase Planck mean and Rosseland mean opacity tables applicable for simulations of star and planet formation, stellar evolution and disc modelling at various metallicities in hydrogen-rich environments. The tables are calculated for gas temperatures between 1000 and 10 000 K and total hydrogen number densities between 10² and  10¹⁷ cm⁻³  . The carbon-to-oxygen ratio is varied from 0.43 to well above 2.0, the nitrogen-to-oxygen ratio between 0.14 and 100.0. The tables are calculated for a range of metallicities down to  [M/H]'= log  N _M/ N _H=−7.0  . We demonstrate how the mean opacities and the abundances of the opacity species vary with C/O, N/O and  [M/H]'  . We use the element abundances from Grevesse et al., and we provide additional tables for the oxygen abundance value from Caffau et al. All tables will be available online at http://star-www.st-and.ac.uk/~ch80/datasources.html .     

Evidence for a massive dark object in NGC 4350

In this work we build a detailed dynamic model for an S0 galaxy possibly hosting a central massive dark object (MDO). We show that the photometric profiles and the kinematics along the major and minor axes, including the h ₃ and h ₄ profiles, imply the presence of a central MDO of mass     i.e. 0.3–2.8 per cent of the mass derived for the stellar spheroidal component. Models without MDO are unable to reproduce the kinematic properties of the inner stars and of the rapidly rotating nuclear gas.
The stellar population consists of an exponential disc (27 per cent of the light) and a diffuse spheroidal component (73 per cent of the light) that cannot be represented by a simple de Vaucouleurs profile at any radius. The M L ratios we found for the stellar components (3.3 and 6.6 respectively) are typical of those of disc and elliptical galaxies.     

Clues for the origin of the fundamental metallicity relations – I. The hierarchical building up of the structure

We analyse the evolutionary history of galaxies formed in a hierarchical scenario consistent with the concordance Lambda cold dark matter (ΛCDM) model focusing on the study of the relation between their chemical and dynamical properties. Our simulations consistently describe the formation of the structure and its chemical enrichment within a cosmological context. Our results indicate that the luminosity–metallicity and the stellar mass–metallicity (LZR and MZR) relations are naturally generated in a hierarchical scenario. Both relations are found to evolve with redshift. In the case of the MZR, the estimated evolution is weaker than that deduced from observational works by approximately 0.10 dex. We also determine a characteristic stellar mass, M _c≈ 3 × 10¹⁰ M_⊙, which segregates the simulated galaxy population into two distinctive groups and which remains unchanged since z ∼ 3, with a very weak evolution of its metallicity content. The value and role played by M _c is consistent with the characteristic mass estimated from the SDSS galaxy survey by Kauffmann et al. Our findings suggest that systems with stellar masses smaller than M _c are responsible for the evolution of this relation at least from z ≈ 3. Larger systems are stellar dominated and have formed more than 50 per cent of their stars at   z ≥ 2  , showing very weak evolution since this epoch. We also found bimodal metallicity and age distributions from z ∼ 3, which reflects the existence of two different galaxy populations. Although SN feedback may affect the properties of galaxies and help to shape the MZR, it is unlikely that it will significantly modify M _c since, from   z = 3  this stellar mass is found in systems with circular velocities larger than 100 km s⁻¹.     

Active galactic nuclei jet mass loading and truncation by stellar winds

Active galactic nuclei can produce extremely powerful jets. While tightly collimated, the scale of these jets and the stellar density at galactic centres implies that there will be many jet/star interactions, which can mass load the jet through stellar winds. Previous work employed modest wind mass outflow rates, but this does not apply when mass loading is provided by a small number of high mass-loss stars. We construct a framework for jet mass loading by stellar winds for a broader spectrum of wind mass-loss rates than has previously been considered. Given the observed stellar mass distributions in galactic centres, we find that even highly efficient (0.1 Eddington luminosity) jets from supermassive black holes of masses M _BH≲ 10⁴ M_⊙ are rapidly mass loaded and quenched by stellar winds. For  10⁴ M_⊙ < M _BH < 10⁸ M_⊙  , the quenching length of highly efficient jets is independent of the jet's mechanical luminosity. Stellar wind mass loading is unable to quench efficient jets from more massive engines, but can account for the observed truncation of the inefficient M87 jet, and implies a baryon-dominated composition on scales ≳2 kpc therein even if the jet is initially pair plasma dominated.     

A Bayesian Kepler periodogram detects a second planet in HD 208487

An automatic Bayesian Kepler periodogram has been developed for identifying and characterizing multiple planetary orbits in precision radial velocity data. The periodogram is powered by a parallel tempering Markov chain Monte Carlo (MCMC) algorithm which is capable of efficiently exploring a multiplanet model parameter space. The periodogram employs an alternative method for converting the time of an observation to true anomaly that enables it to handle much larger data sets without a significant increase in computation time. Improvements in the periodogram and further tests using data from HD 208487 have resulted in the detection of a second planet with a period of 909⁸²₋₉₂ d, an eccentricity of 0.37^0.26_−0.20, a semimajor axis of 1.87^0.13_−0.14 au and an M sin  i = 0.45^+0.11_−0.13 M _J. The revised parameters of the first planet are period = 129.8 ± 0.4 d, eccentricity = 0.20 ± 0.09, semimajor axis = 0.51 ± 0.02 au and M sin  i = 0.41 ± 0.05  M _J. Particular attention is paid to several methods for calculating the model marginal likelihood which is used to compare the probabilities of models with different numbers of planets.     

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in modified Newtonian dynamics (MOND), using the N -body code n-mody , which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo & Nipoti. We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of 10⁴ to  10⁹ M_⊙  and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime  ( a _i, a _e≪ a ₀)  , where the motion of stars is either dominated by internal accelerations  ( a _i≫ a _e)  or constant external accelerations  ( a _e≫ a _i)  . In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime  ( a _i∼ a _e∼ a ₀)  . This allows for a much-improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.     

Infrared helium–hydrogen line ratios as a measure of stellar effective temperature

We have observed a large sample of compact planetary nebulae in the near-infrared to determine how the 2¹P–2¹S He  i line at 2.058 μm varies as a function of stellar effective temperature, T _eff. The ratio of this line with H  i Br γ at 2.166 μm has often been used as a measure of the highest T _eff present in a stellar cluster, and hence of whether there is a cut-off in the stellar initial mass function at high masses. However, recent photoionization modelling has revealed that the behaviour of this line is more complex than previously anticipated. Our work shows that in most aspects the photoionization models are correct. In particular, we confirm the weakening of the 2¹P–2¹S line as T _eff increases beyond 40 000 K. However, in many cases the model underpredicts the observed ratio when we consider the detailed physical conditions in the individual planetary nebulae. Furthermore, there is evidence that there is still significant 2¹P–2¹S He  i line emission even in the planetary nebulae with very hot     central stars. It is clear from our work that this ratio cannot be considered as a reliable measure of effective temperature on its own.     

Compton dragged gamma-ray bursts: the spectrum

Recently several studies have jointly analysed data from different cosmological probes with the motivation of estimating cosmological parameters. Here we generalize this procedure to allow freedom in the relative weights of various probes. This is done by including in the joint χ ² function a set of 'hyper-parameters', which are dealt with using Bayesian considerations. The resulting algorithm, which assumes uniform priors on the log of the hyper-parameters, is very simple: instead of minimizing     (where     is per data set j ) we propose to minimize     (where N _j is the number of data points per data set j ). We illustrate the method by estimating the Hubble constant H ₀ from different sets of recent cosmic microwave background (CMB) experiments (including Saskatoon, Python V, MSAM1, TOCO and Boomerang ). The approach can be generalized for combinations of cosmic probes, and for other priors on the hyper-parameters.     

Quantifying the coexistence of massive black holes and dense nuclear star clusters

In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  10⁶–10⁹ M_⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  10⁵–10⁷ M_⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  10⁸–10¹¹ M_⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M _BH+ M _NC)/ M _sph]=−(0.39 ± 0.07) log[ M _sph/10¹⁰ M_⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.     

Wilkinson Microwave Anisotropy Probe 5-yr constraints on fnl with wavelets

We present a Gaussianity analysis of the Wilkinson Microwave Anisotropy Probe ( WMAP ) 5-yr cosmic microwave background (CMB) temperature anisotropy data maps. We use several third-order estimators based on the spherical Mexican hat wavelet. We impose constraints on the local non-linear coupling parameter f _nl using well-motivated non-Gaussian simulations. We analyse the WMAP maps at resolution of 6.9 arcmin for the Q , V , and W frequency bands. We use the KQ 75 mask recommended by the WMAP team which masks out 28 per cent of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 150 arcmin. With these coefficients, we compute the third order estimators which are used to perform a  χ²  analysis. The  χ²  statistic is used to test the Gaussianity of the WMAP data as well as to constrain the f _nl parameter. Our results indicate that the WMAP data are compatible with the Gaussian simulations, and the f _nl parameter is constrained to  −8 < f _nl < +111  at 95 per cent confidence level (CL) for the combined   V + W   map. This value has been corrected for the presence of undetected point sources, which add a positive contribution of  Δ f _nl= 3 ± 5  in the   V + W   map. Our results are very similar to those obtained by the WMAP team using the bispectrum.     

Resonant excitation of white dwarf oscillations in compact object binaries – I. The no back reaction approximation

We consider the evolution of white dwarfs with compact object companions (specifically black holes with masses up to  ∼10⁶  M_⊙  , neutron stars, and other white dwarfs). We suppose that the orbits are initially quite elliptical and then shrink and circularize under the action of gravitational radiation. During this evolution, the white dwarfs will pass through resonances when harmonics of the orbital frequency match the stellar oscillation eigenfrequencies. As a star passes through these resonances, the associated modes will be excited and can be driven to amplitudes that are so large that there is a back reaction on the orbit which, in turn, limits the growth of the modes. A formalism is presented for describing this dynamical interaction for a non-rotating star in the linear approximation when the orbit can be treated as non-relativistic. A semi-analytical expression is found for computing the resonant energy transfer as a function of stellar and orbital parameters for the regime where back reaction may be neglected. This is used to calculate the results of passage through a sequence of resonances for several hypothetical systems. It is found that the amplitude of the  ℓ= m = 2  f -mode can be driven into the non-linear regime for appropriate initial conditions. We also discuss where the no back reaction approximation is expected to fail, and the qualitative effects of back reaction.