A comparison between terrestrial, Cytherean and lunar impact cratering records

A comparison between the terrestrial, Cytherean and lunar cratering records indicates that the large craters (diameters   D > D ₀)  on these surfaces all have cumulative numbers that are proportional to   D ^-2.59±0.05  . Atmospheres have a negligible effect on the formation of   D > D ₀  craters. It is shown that this limiting diameter is  45±3 km  in the case of Venus, and  21.0±1.5 km  in the case of Earth. In this large-diameter range, there are about  1.51±0.34  times more craters, per unit area, on Venus than on the Earth, and about  1350±310  times more craters, per unit area, on the Moon than on the Earth.     

The dynamics of Neptune Trojan – I. The inclined orbits

The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the fast Fourier transform technique, we construct high-resolution dynamical maps on the plane of initial semimajor axis a ₀ versus inclination i ₀. These maps show three most stable regions, with i ₀ in the range of  (0°, 12°), (22°, 36°)  and  (51°, 59°),  respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points L ₄ and L ₅ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the  ν₈  secular resonance around   i ₀∼ 44°  pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination  (>60°)  and an unstable gap around  44°  on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of  ( a ₀, i ₀)  . The fine structures in the dynamical maps can be explained by these secular resonances.     

Simultaneous impact and lunar craters

The existence of large terrestrial impact crater doublets and Martian crater doublets that have been inferred to be impact craters demonstrates that simultaneous impact of two or more bodies occurs at nearly the same point on planetary surfaces. An experimental study of simultaneous impact of two projectiles near one another shows that doublet craters with ridges perpendicular to the bilateral axis of symmetry result when separation between impact points relative to individual crater diameter is large. When separation is progressively less, elliptical craters with central ridges and central peaks, circular craters with flat floors containing ridges and peaks, and circular craters with deep round bottoms are produced. These craters are similar in structure to many of the large lunar craters. Results suggest that the simultaneous impact of meteoroids near one another may be an important mechanism for the production of central peaks in large lunar craters.     

Principal components in active galactic nuclei variability data and the estimation of the flux contributions from different components

It has been found that the near-infrared flux variations of Seyfert galaxies satisfy relations of the form   F_i ≈α _{i j} +β _{i j} F_j   , where F_i , F_j are the fluxes in filters i and j ; and  α _{i , j} , β _{i , j}   are constants. These relations have been used to estimate the constant contributions of the non-variable underlying galaxies. The paper attempts a formal treatment of the estimation procedure, allowing for the possible presence of a third component, namely non-variable hot dust. In an analysis of a sample of 38 Seyfert galaxies, inclusion of the hot dust component improves the model fit in approximately half the cases. All derived dust temperatures are below 300 K, in the range 540–860 K or above 1300 K. A noteworthy feature is the estimation of confidence intervals for the component contributions: this is achieved by bootstrapping. It is also pointed out that the model implies that such data could be fruitfully analysed in terms of principal components.     

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.     

Contribution to the search for binaries among Am stars – VIII. New spectroscopic orbits of eight systems and statistical study of a sample of 91 Am stars

Red clump giant (RCG) stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration data base to constrain analytic triaxial models for the Galactic bar. We find the bar major-axis is oriented at an angle of 24°–27° to the Sun–Galactic Centre line-of-sight. The ratio of semimajor and semiminor bar axis scalelengths in the Galactic plane   x ₀, y ₀  , and vertical bar scalelength z ₀, is   x ₀ :  y ₀ :  z ₀= 10 : 3.5 : 2.6  , suggesting a slightly more prolate bar structure than the working model of Gerhard which gives the scalelength ratios as   x ₀ :  y ₀ :  z ₀= 10 : 4 : 3  .     

Degeneracy of semi-analytic solution for axially symmetric lens equations

We derive a simple semi-analytical approximation for lens equations with an arbitrary radially symmetric mass density ρ( r ), when   r /ξ₀≪ 1  and ξ₀ is the scalelength of the density profile. At the strong lensing regime, which is mostly constrained by the inner part of the mass density profile, we assume ρ∝ r ^α.
A dark matter (DM) haloes (GNFW model) are parametrized through a shape parameter α, a concentration parameter c ₁ and the total mass M . We apply our semi-analytical model to show how the solutions of the axially symmetric lens equations are degenerated in respect to the parameters α and c ₁.
In the case of an asymmetric dual image lens system, similar effective degeneracy is produced when the geometry of the lens is relaxed. Because it is impossible to determine the exact location of the source image, a family of solutions is acquired when the mass of the lens object and location of the observed images are fixed.
Our results indicate that the amount of degeneration is only weakly affected by the asymmetry in the lensing geometry set-up, e.g. the observational effective degeneracy is very close to the true physical degeneracy of the Einstein ring solutions. Basically with high-enough values for the concentration parameter, the degeneracy spawns the whole range for the shape parameter  α=[−2.0, −1.0]  .     

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.     

Electron density diagnostic potential of the Si x ion and its application in Procyon

Theoretical electron density sensitive line ratios   R ₁– R ₆  of Si  x soft X-ray emission lines are presented. We found that these line ratios are sensitive to electron density n _e, and the ratio R ₁ is insensitive to electron temperature T _e. For reliable determination of the electron density of laboratory and astrophysical plasmas, atomic data, such as electron impact excitation rates, are very important. Our results reveal that the discrepancy of the line ratios from different atomic data calculated with the distorted wave (DW) approximation and the R-matrix method is up to 19 per cent at   n _e= 2 × 10⁸ cm⁻³  . We applied the theoretical intensity ratio R ₁ to the Low Energy Transmission Grating Spectrometer (LETGS) spectrum of the solar-like star Procyon. By comparing the observed value (1.29) with the theoretical calculation, the derived electron density n _e is  2.6 × 10⁸ cm⁻³  , which is consistent with that derived from  (C  v < 8.3 × 10⁸ cm⁻³)  . When the temperature structure of the Procyon corona is taken into account, the derived electron density increases from   n _e= 2.6 × 10⁸  to  2.8 × 10⁸ cm⁻³  .     

Survival and disruption of subsystems during a cold collapse

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses under its own gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made a systematic study in two-parameter space: the number of the clumps N _c and the size of the clump r _v . We obtained the condition N _c ≫ 1 and n _k  ≥ 1, where n _k is related to r _v and the initial radius of the cluster R _ini through the relation R _ini/ r _v  = 2 N ^{( n k +5)/6}_c. A simple analytical argument supports the numerical result. This n _k corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than 2/Ω h ²)Mpc, the first violent collapse is strong enough to sweep away all the substructures that exist before the collapse.     

An asymmetric relativistic model for classical double radio sources

There is substantial observational evidence against the symmetric relativistic model of FR II radio sources. An asymmetric relativistic model is proposed which takes account of both relativistic effects and intrinsic/environmental asymmetries to explain the structural asymmetries of their radio lobes. A key parameter of the model is the jet-side of the double sources, which is estimated for 80 per cent of the FR II sources in the 3CRR complete sample. Statistical analyses of the properties of these sources show that the asymmetric model is in agreement with a wide range of observational data, and that the relativistic and intrinsic asymmetry effects are of comparable importance. Intrinsic/environmental asymmetry effects are more important at high radio luminosities and small physical scales. The mean translational speed of the lobes is found to be     consistent with the speeds found from spectral ageing arguments. According to a Gaussian model, the standard deviation of the distribution of v _lobe is σ _{v l}=0.04 c . The results are in agreement with an orientation-based unification scheme in which the critical angle separating the radio galaxies from quasars is about 45°.     

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.     

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.     

Maximum likelihood estimates of the two- and three-dimensional power spectra of the APM Galaxy Survey

We estimate the two- and three-dimensional power spectra, P ₂( K ) and P ₃( k ), of the galaxy distribution by applying a maximum likelihood estimator to pixel maps of the APM Galaxy Survey. The analysis provides optimal estimates of the power spectra and of their covariance matrices if the fluctuations are assumed to be Gaussian. Our estimates of P ₂( K ) and P ₃( k ) are in good agreement with previous work, but we find that the errors at low wavenumbers have been underestimated in some earlier studies. If the galaxy power spectrum is assumed to have the same shape as the mass power spectrum, then the APM maximum likelihood P ₃( k ) estimates at k ≤0.19  h  Mpc⁻¹ constrain the amplitude and shape parameter of a scale-invariant CDM model to lie within the 2 σ ranges 0.74≤( σ ₈)_g≤1.28 and 0.06≤Γ≤0.46 . Using the Galactic extinction estimates of Schlegel, Finkbeiner & Davis, we show that Galactic obscuration has a negligible effect on galaxy clustering over most of the area of the APM Galaxy Survey.     

Spectroscopic solution of the star QX Cas

High-resolution spectroscopic observations around the Hα line of the binary star QX Cas covering the whole orbital period are presented. Our radial velocity solution, the first ever determined, requires an eccentric orbit with the following orbital parameters: eccentricity,   e = 0.22 ± 0.01  ; longitude of periastron,  ω= 45°± 5°  ; semi-amplitudes of the radial velocity curves of the primary and secondary stars,   K ₁ sin  i = 125.8 ± 0.9 km s⁻¹  and   K ₂ sin  i = 144.8 ± 1.1 km s⁻¹  ; gamma velocity,   V ₀= 65.1 ± 0.5 km s⁻¹  ; and mass ratio,   q = 0.869 ± 0.013  . The corresponding lower limits of the masses of the components and their separation are         , and   a sin  i = 31.34 ± 0.48 R_⊙  .     

Testing tidal-torque theory – I. Spin amplitude and direction

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N -body simulation with thousands of well-resolved haloes. The protohaloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the protohaloes match, on average, the spin amplitudes of the virialized haloes of today, if linear growth is assumed until ∼ t ₀/3, or  55–70  per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ∼50° in the TTT prediction at t ₀, such that the linear spatial correlations of spins on scales ≥1  h ⁻¹ Mpc are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second-order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.     

Monte Carlo simulations of star clusters – II. Tidally limited, multimass systems with stellar evolution

A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W ₀=3). Models with larger W ₀ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.     

Roche tomography of cataclysmic variables – III. Star-spots on AE Aqr

We present a Roche tomography reconstruction of the secondary star in the cataclysmic variable AE Aqr. The tomogram reveals several surface inhomogeneities that are due to the presence of large, cool star-spots. In addition to a number of lower latitude spots, the maps also show the presence of a large, high-latitude spot similar to that seen in Doppler images of rapidly rotating isolated stars, and a relative paucity of spots at a latitude of 40°. In total, we estimate that some 18 per cent of the Northern hemisphere of AE Aqr is spotted.
We have also applied the entropy landscape technique to determine accurate parameters for the binary system. We obtain optimal masses   M ₁= 0.74 M_⊙, M ₂= 0.50 M_⊙  , a systemic velocity  γ=−63 km s⁻¹  and an orbital inclination   i = 66°  .
Given that this is the first study to successfully image star-spots on the secondary star in a cataclysmic variable, we discuss the role that further studies of this kind may play in our understanding of these binaries.     

On the relation between orbital structure and observed bar morphology

We investigate the morphological relation between the orbits of the central family of periodic orbits ( x ₁ family) and the bar itself using models of test particles moving in a barred potential. We show that different bar morphologies may have as a backbone the same set of x ₁ periodic orbits. We point out that by populating initially axisymmetric stellar discs exponentially with test particles in circular, or almost circular motion, we may end up with a response bar which reveals a shape different in crucial details from that of the individual stable x ₁ orbits. For example, a bar model in which the x ₁ orbits are pure ellipses may have a much more complicated response morphology. This depends on the particular invariant curves around x ₁, which are populated in each model.     

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.