Determining bias with cumulant correlators

correlator of the galaxy density field Q ₂₁ is examined from the point of view of biasing. It is shown that, to leading order, it depends on two biasing parameters b b ₂, and on q ₂₁, the underlying cumulant correlator of the mass. As the skewness Q ₃ has analogous properties, the slope of the correlation function −γ, Q ₃ and Q ₂₁ uniquely determine the bias parameter on a particular scale to be b  = γ/6( Q ₂₁ −  Q ₃), when working in the context of gravitational instability with Gaussian initial conditions. Thus on large scales, easily accessible with the future Sloan Digital Sky Survey and the 2 Degree Field Survey, it will be possible to extract b b ₂ from simple counts-in-cells measurements. Moreover, the higher order cumulants, Q _N , successively determine the higher order biasing parameters. From these it is possible to predict higher order cumulant correlators as well. Comparison of the predictions with the measurements will provide internal consistency checks on the validity of the assumptions in the theory, most notably perturbation theory of the growth of fluctuations by gravity and Gaussian initial conditions. Since the method is insensitive Ω, it can be successfully combined with results from velocity fields, which determine Ω^0.6/b, to measure the total density parameter in the Universe.     

The shape of the dark matter halo in the early-type galaxy NGC 2974

We present H  i observations of the elliptical galaxy NGC 2974, obtained with the Very Large Array. These observations reveal that the previously detected H  i disc in this galaxy is in fact a ring. By studying the harmonic expansion of the velocity field along the ring, we constrain the elongation of the halo and find that the underlying gravitational potential is consistent with an axisymmetric shape.
We construct mass models of NGC 2974 by combining the H  i rotation curve with the central kinematics of the ionized gas, obtained with the integral-field spectrograph SAURON. We introduce a new way of correcting the observed velocities of the ionized gas for asymmetric drift, and hereby disentangle the random motions of the gas caused by gravitational interaction from those caused by turbulence. To reproduce the observed flat rotation curve of the H  i gas, we need to include a dark halo in our mass models. A pseudo-isothermal sphere provides the best model to fit our data, but we also tested an NFW halo and modified Newtonian dynamics, which fit the data marginally worse.
The mass-to-light ratio M / L _I increases in NGC 2974 from 4.3 M_⊙/L_⊙, _I at one effective radius to 8.5 M_⊙/L_{⊙, I} at 5  R _e. This increase of M / L already suggests the presence of dark matter: we find that within 5  R _e at least 55 per cent of the total mass is dark.     

Gaseous drag on a gravitational perturber in Modified Newtonian Dynamics and the structure of the wake

We calculate the structure of a wake generated by, and the dynamical friction force on, a gravitational perturber travelling through a gaseous medium of uniform density and constant background acceleration   g _ext  , in the context of Modified Newtonian Dynamics (MOND). The wake is described as a linear superposition of two terms. The dominant part displays the same structure as the wake generated in the Newtonian gravity scaled up by a factor  μ⁻¹( g _ext/ a ₀)  , where a ₀ is the constant MOND acceleration and μ the interpolating function. The structure of the second term depends greatly on the angle between   g _ext  and the velocity of the perturber. We evaluate the dynamical drag force numerically and compare our MOND results with the Newtonian case. We mention the relevance of our calculations to orbit evolution of globular clusters and satellites in a gaseous protogalaxy. Potential differences in the X-ray emission of gravitational galactic wakes in MOND and in Newtonian gravity with a dark halo are highlighted.     

On the relation between electron temperatures in the O+ and O++ zones in high-metallicity H ii regions

We suggest a new way to establish the relation between the electron temperature t ₃ within the [O  iii ] zone and the electron temperature t ₂ within the [O  ii ] zone in high-metallicity  (12 + log(O/H) > 8.25)  H  ii regions. The   t ₂– t ₃  diagram is constructed by applying our method to a sample of 372 H  ii regions. We find that the correlation between t ₂ and t ₃ is tight and can be approximated by a linear expression. The new   t ₂– t ₃  relation can be used to determine t ₂ and accurate abundances in high-metallicity H  ii regions with a measured t ₃. It can also be used in conjunction with the ff relation for the determination of t ₃ and t ₂ and oxygen abundances in high-metallicity H  ii regions, where the [O  iii ]λ4363 auroral line is not detected. The derived   t ₂– t ₃  relation is independent of photoionization models of H  ii regions.     

The mass ratio and the orbital parameters of the sdOB binary AA Doradus

The time sequence of 105 spectra covering one full orbital period of AA Dor has been analysed. Direct determination of   V  sin  i   for the sdOB component from 97 spectra outside of the eclipse for the lines Mg  ii 4481 Å and Si  iv 4089 Å clearly indicated a substantially smaller value than estimated before. Detailed modelling of line-profile variations for eight spectra during the eclipse for the Mg  ii 4481 Å line, combined with the out-of-eclipse fits, gave   V  sin  i = 31.8 ± 1.8 km s⁻¹  . The previous determinations of   V  sin  i   , based on the He  ii 4686 Å line, appear to be invalid because of the large natural broadening of the line. With the assumption of the solid-body, synchronous rotation of the sdOB primary, the measured values of the semi-amplitude K ₁ and   V  sin  i   lead to the mass ratio   q = 0.213 ± 0.013  which in turn gives K ₂ and thus the masses and radii of both components. The sdOB component appears to be less massive than assumed before,   M ₁= 0.25 ± 0.05 M_⊙  , but the secondary has its mass–radius parameters close to theoretically predicted for a brown dwarf,   M ₂= 0.054 ± 0.010 M_⊙  and   R ₂= 0.089 ± 0.005 R_⊙  . Our results do not agree with the recent determination of Vŭcković et al. based on a K ₂ estimate from line-profile asymmetries.     

Variance and skewness in the FIRST survey

We investigate the large-scale clustering of radio sources in the FIRST 1.4-GHz survey by analysing the distribution function ( counts in cells ). We select a reliable sample from the the FIRST catalogue, paying particular attention to the problem of how to define single radio sources from the multiple components listed. We also consider the incompleteness of the catalogue. We estimate the angular two-point correlation function w (θ), the variance Ψ₂ and skewness Ψ₃ of the distribution for the various subsamples chosen on different criteria. Both w (θ) and Ψ₂ show power-law behaviour with an amplitude corresponding to a spatial correlation length of r ₀ ∼ 10  h ⁻¹Mpc. We detect significant skewness in the distribution, the first such detection in radio surveys. This skewness is found to be related to the variance through Ψ₃ =  S ₃(Ψ₂)^α, with α = 1.9 ± 0.1, consistent with the non-linear gravitational growth of perturbations from primordial Gaussian initial conditions. We show that the amplitude of variance and the skewness are consistent with realistic models of galaxy clustering.     

The radial Tully–Fisher relation for spiral galaxies – I

We find a new Tully–Fisher-like relation for spiral galaxies holding at different galactocentric radii. This radial Tully–Fisher relation allows us to investigate the distribution of matter in the optical regions of spiral galaxies. This relation, applied to three different samples of rotation curves of spiral galaxies, directly proves that: (i) the rotation velocity of spirals is a good measure of their gravitational potential and both the rotation curve's amplitudes and profiles are well predicted by galaxy luminosity, (ii) the existence of a dark component, less concentrated than the luminous one, and (iii) a scaling law, according to which, inside the disc optical size:   M _dark/ M _lum= 0.5( L _B /10¹¹ L _{B ⊙})^−0.7  .     

Investigating cosmological weak lensing with the ray-bundle method

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.     

The minimum mass ratio of W UMa-type binary systems

When the total angular momentum of a binary system   J _tot= J _orb+ J _spin  is at a certain critical (minimum) value, a tidal instability occurs which eventually forces the stars to merge into a single, rapidly rotating object. The instability occurs when   J _orb= 3 J _spin  , which in the case of contact binaries corresponds to a minimum mass ratio   q _min≈  0.071–0.078. The minimum mass ratio is obtained under the assumption that stellar radii are fixed and independent. This is not the case with contact binaries where, according to the Roche model, we have   R ₂= R ₂( R ₁, a , q )  . By finding a new criterion for contact binaries, which arises from  d J _tot= 0  , and assuming   k ²₁≠ k ²₂  for the component's dimensionless gyration radii, a theoretical lower limit   q _min= 0.094–0.109  for overcontact degree   f = 0–1  is obtained.     

On the transfer of radiation at asteroidal surfaces in relation to their orbit deflection – I

One of the methods discussed in deflecting the orbit of an Earth-colliding asteroid is the use of nuclear explosives. In assessing its feasibility, apart from political considerations, it is important to quantify how effective it is in orbit deflection. The transfer of radiation incident at the surface is governed by a non-linear diffusion equation. For low-yield explosions with a slab geometry ( S ₀≃10⁸ kJ μs⁻¹), the temperature at depth x and time t is well approximated by a similarity solution of the form T ( x , t )= T ₀ f (ξ), ξ= x/ ( T ⁿ ₀ t )^1/2, with T ₀ given by ( S ₀/σ)^1/4, where σ is the Stefan–Boltzmann constant, n is an index that specifies the radiation transfer and f (ξ) is the solution of a non-linear differential equation subject to the condition f (0)=1 and lim_ξ→∞ f (ξ)=0. For high-yield explosions ( S ₀≃10¹⁰ kJ μs⁻¹), numerical solutions to the non-linear diffusion equation can be obtained. These solutions have properties similar to the case of low-yield explosions. If the duration of the explosion is d ×10⁻⁸ s, where d is close to 3, the fraction of energy absorbed by the surface is found to be 7, 12 and 23 per cent for S ₀=10⁸, 10⁹ and 10¹⁰ kJ μs⁻¹ respectively.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

Lensing, reddening and extinction effects of Mg ii absorbers from z= 0.4 to 2

Using a sample of almost 7000 strong Mg  ii absorbers with   W ₀ > 1 Å  and  0.4 < z < 2.2  detected in the SDSS DR4 data set, we investigate the gravitational lensing and dust extinction effects they induce on background quasars. After carefully quantifying several selection biases, we isolate the reddening effects as a function of redshift and absorber rest equivalent width, W ₀. We find the amount of dust to increase with cosmic time as  τ( z ) ∝ (1 + z )^−1.1±0.4  , following the evolution of cosmic star density or integrated star formation rate. We measure the reddening effects over a factor of 30 in E ( B − V ) and we find that  τ∝ ( W ₀)^1.9±0.1  , providing us with an important scaling for theoretical modelling of metal absorbers. We also measure the dust-to-metal ratio and find it similar to that of the Milky Way. In contrast to previous studies, we do not detect any gravitational magnification by Mg  ii systems. We measure the upper limit  μ < 1.10  and discuss the origin of the discrepancy. Finally, we estimate the fraction of absorbers missed due to extinction effects and show that it rises from 1 to 50 per cent in the range  1 < W ₀ < 6 Å  . We parametrize this effect and provide a correction for recovering the intrinsic  ∂ N /∂ W ₀  distribution.     

Multiple frequency analysis of the large-amplitude SX Phoenicis star BL Camelopardalis

A multifrequency analysis of the SX Phoenicis star BL Camelopardalis is presented on the basis of new high-speed photometry, along with fitting a total of 136 maxima. BL Cam is a multiple periodic pulsator. We find f ₀=25.5768, f ₁=25.2982, f ₂=25.8622, f ₃=31.5912, f ₄=25.1065, f ₅=25.5147 and f ₆=25.6188 cycle d⁻¹ together with the harmonics 51.1513 and 76.7268 cycle d⁻¹ and combination frequencies f ₀+ f ₁, f ₀+ f ₂ and f ₀+ f ₃. The new frequency solution represents the light curves of BL Cam quite well. The observed minus calculated (O-C) analysis indicates that the fundamental frequency is in good agreement with the results of Fourier analysis.     

Nodal and periastron precession of inclined orbits in the field of a rapidly rotating neutron star

We derive a formula for the nodal precession frequency and the Keplerian period of a particle at an arbitrary orbital inclination (with a minimum latitudinal angle reached at the orbit) in the post-Newtonian approximation in the external field of an oblate rotating neutron star (NS). We also derive formulas for the nodal precession and periastron rotation frequencies of slightly inclined low-eccentricity orbits in the field of a rapidly rotating NS in the form of asymptotic expansions whose first terms are given by the Okazaki-Kato formulas. The NS gravitational field is described by the exact solution of the Einstein equation that includes the NS quadrupole moment induced by rapid rotation. Convenient asymptotic formulas are given for the metric coefficients of the corresponding space-time in the form of Kerr metric perturbations in Boyer-Lindquist coordinates.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

An asteroseismic study of the β Cephei star 12 Lacertae: multisite spectroscopic observations, mode identification and seismic modelling

We present the results of a spectroscopic multisite campaign for the β Cephei star 12 (DD) Lacertae. Our study is based on more than thousand high-resolution high S/N spectra gathered with eight different telescopes in a time span of 11 months. In addition, we make use of numerous archival spectroscopic measurements. We confirm 10 independent frequencies recently discovered from photometry, as well as harmonics and combination frequencies. In particular, the slowly pulsating B-stars (SPB)-like g -mode with frequency 0.3428 d⁻¹ reported before is detected in our spectroscopy. We identify the four main modes as  (ℓ₁, m ₁) = (1, 1), (ℓ₂, m ₂) = (0, 0), (ℓ₃, m ₃) = (1, 0)  and  (ℓ₄, m ₄) = (2, 1)  for   f ₁= 5.178 964 d⁻¹, f ₂= 5.334 224 d⁻¹, f ₃= 5.066 316 d⁻¹  and   f ₄= 5.490 133 d⁻¹  , respectively. Our seismic modelling shows that f ₂ is likely the radial first overtone and that the core overshooting parameter  α_ov  is lower than 0.4 local pressure scale heights.     

VLT spectroscopy of XTE J2123−058 during quiescence: the masses of the two components

We present Very Large Telescope (VLT) low-resolution spectroscopy of the neutron star X-ray transient XTE J2123−058 during its quiescent state. Our data reveal the presence of a K7V companion which contributes 77 per cent to the total flux at λ 6300 and orbits the neutron star at     . Contrary to other soft X-ray transients (SXTs), the H _α emission is almost exactly in antiphase with the velocity curve of the optical companion. Using the light-centre technique we obtain     and hence     This, combined with a previous determination of the inclination angle     yields     and     . M ₂ agrees well with the observed spectral type. Doppler tomography of the H α emission shows a non-symmetric accretion disc distribution mimicking that seen in SW Sex stars. Although we find a large systemic velocity of −     this value is consistent with the galactic rotation velocity at the position of J2123−058, and hence a halo origin. The formation scenario of J2123−058 is still unresolved.     

Dynamics of the boxy elliptical galaxy NGC 1600

We use three-integral models to infer the distribution function (DF) of the boxy E3–E4 galaxy NGC 1600 from surface brightness and line-profile data on the minor and major axes. We assume axisymmetry and that the mass-to-light ratio is constant in the central ∼1 R _e. Stars in the resulting gravitational potential move mainly on regular orbits. We use an approximate third integral K from perturbation theory and write the DF as a sum of basis functions in the three integrals E , L _z and K . We then fit the projected moments of these basis functions to the kinematic observables and deprojected density, using a non-parametric algorithm. The deduced dynamical structure is radially anisotropic, with σ _θ σ _r ≈ σ _φ σ _r ≈0.7 on the major axis. Both on the minor axis and near the centre the velocity distribution is more isotropic; thus the model is flattened by equatorial radial orbits. The kinematic data are fitted without the need for a central black hole; the central mass determined previously from ground-based data therefore overestimates the actual black-hole mass. The mass-to-light ratio of the stars is M L _V =6  h ₅₀. The anisotropy structure of NGC 1600 with a radially anisotropic main body and more nearly isotropic centre is similar to that found recently in NGC 1399, 2434, 3379 and 6703, suggesting that this pattern may be common amongst massive elliptical galaxies. We discuss a possible merger origin of NGC 1600 in the light of these results.     

Bulk viscosity of mixed nucleon–hyperon–quark matter in neutron stars

We calculate the coefficient of bulk viscosity by considering the non-leptonic weak interactions in the cores of hybrid stars with both hyperons and quarks. We first determine the dependence of the production rate of neutrons on the reaction rate of quarks in the non-leptonic processes, that is,  Γ _n = K _s Γ _s +Γ_Λ+ 2Γ_Σ−  . The conversion rate,   K _s   , in our scenario is a complicated function of baryon number density. We also consider the medium effect of quark matter on bulk viscosity. Using these results, we estimate the limiting rotation of the hybrid stars, which may suppress the r-mode instability more effectively. Hybrid stars should be the candidates for the extremely rapid rotators.