On the detectability of lightcurves of Kuiper belt objects

We present a statistical study of the detectability of lightcurves of Kuiper belt objects (KBOs). Some Kuiper belt objects display lightcurves that appear “flat”; i.e., there are no significant brightness variations within the photometric uncertainties. Under the assumption that KBO lightcurves are mainly due to shape, the lack of brightness variations may be due to (1) the objects having very nearly spherical shapes or (2) their rotation axes coinciding with the line of sight. We investigate the relative importance of these two effects and relate it to the observed fraction of “flat” lightcurves. This study suggests that the fraction of KBOs with detectable brightness variations may provide clues about the shape distribution of these objects. Although the current database of rotational properties of KBOs is still insufficient to draw any statistically meaningful conclusions, we expect that, with a larger dataset, this method will provide a useful test for candidate KBO shape distributions.     

Oscillations in radiation-driven outflows

In this paper we investigate the dynamical behaviour of radiation-driven winds, specifically winds that arise when Compton scattering transfers momentum from the radiation field to the gas flow. Such winds occur during strong X-ray bursts from slowly accreting neutron stars, and also may be driven from the inner regions of a black hole or neutron star accretion disc when the mass transfer rate is very high. By linearizing the radiation hydrodynamic equations around steady spherical outflow, we evaluate the time-dependent response of these winds to perturbations introduced at their inner boundaries. We find that although radiation-driven winds are generally stable, they act as mechanical filters that should produce quasi-periodic oscillations or peaked noise in their radiation output when perturbations force them stochastically. This behaviour may underlie the photospheric oscillations observed in some strong Type I X-ray bursts.     

Intermittent Behaviour in Axisymmetric Mean-Field Dynamo Models in Spherical Shells

Axisymmetric mean-field dynamo models in spherical shells are shown to be capable of producing temporally intermittent behaviour. This is of potential importance since (i) it is, as far as we are aware, the first time such behaviour has been produced internally by a mean-field dynamo model in a spherical shell, without requiring any additional assumptions or truncations, and (ii) it may be characteristic of the type of behaviour observed in the long-term record of solar activity, such as Maunder minima. We also show that these types of behaviour persist when the functional form of the alpha quenching is altered and also occur over intervals of the shell thickness and the dynamo number.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Projected bispectrum in spherical harmonics and its application to angular galaxy catalogues

We present a theoretical and exact analysis of the bispectrum of projected galaxy catalogues. The result can be generalized to evaluate the projection in spherical harmonics of any 3D bispectrum and therefore has applications to cosmic microwave background and gravitational lensing studies.
By expanding the 2D distribution of galaxies on the sky in spherical harmonics, we show how the three-point function of the coefficients can be used in principle to determine the bias parameter of the galaxy sample. If this can be achieved, it would allow a lifting of the degeneracy between the bias and the matter density parameter of the Universe, which occurs in linear analysis of 3D galaxy catalogues. In previous papers, we have shown how a similar analysis can be done in three dimensions, and we show here through an error analysis and by implementing the method on a simulated projected catalogue that ongoing three-dimensional galaxy redshift surveys (even with all the additional uncertainties introduced by partial sky coverage, redshift-space distortions and smaller numbers) will do far better than all-sky projected catalogues with similar selection function.     

A possible long-lived asteroid population at the equilateral Lagrangian points of Saturn

The Lagrangian equilateral points of a planetary orbit are points of equilibrium that trail at 60°, ahead (L4) or behind (L5), the trajectory of a planet. Jupiter is the only major planet in our Solar system harbouring a known population of asteroids at those locations. Here we report the existence of orbits close to the Lagrangian points of Saturn, stable at time-scales comparable to the age of the Solar system. By scaling with respect to the Trojan population we have estimated the number of objects that would populate the regions, which gives a significant figure. Moreover, mutual physical collisions over the age of the Solar system would be very rare, so the evaporation rate of this swarm arising from mutual interactions would be very low. A population of asteroids not self-collisionally evolved after their formation stage would be the first to be observed in our planetary system. Our present estimations are based on the assumption that the capture efficiency at Saturn's equilateral points is comparable with the one corresponding to Jupiter, thus our figures may be taken as upper limits. In any case, observational constraints on their number would provide fundamental clues to our understanding of the history of the outer Solar system. If they existed, the surface properties and size distribution of those objects would represent unusually valuable fossil records of our early planetary system.     

Simultaneous multifrequency single-pulse properties of AXP XTE J1810–197

We have used the 76-m Lovell, 94-m equivalent Westerbork Synthesis Radio Telescope (WSRT) and 100-m Effelsberg radio telescopes to investigate the simultaneous single-pulse properties of the radio emitting magnetar Anomalous X-ray Pulsar (AXP) XTE J1810−197 at frequencies of 1.4, 4.8 and 8.35 GHz during 2006 May and July. We study the magnetar's pulse-energy distributions which are found to be very peculiar as they are changing on time-scales of days and cannot be fit by a single statistical model. The magnetar exhibits strong spiky single giant-pulse-like subpulses, but they do not fit the definition of the giant pulse or giant micropulse phenomena. Measurements of the longitude-resolved modulation index reveal a high degree of intensity fluctuations on day-to-day time-scales and dramatic changes across pulse phase. We find the frequency evolution of the modulation index values differs significantly from what is observed in normal radio pulsars. We find that no regular drifting subpulse phenomenon is present at any of the observed frequencies at any observing epoch. However, we find a quasi-periodicity of the subpulses present in the majority of the observing sessions. A correlation analysis indicates a relationship between components from different frequencies. We discuss the results of our analysis in light of the emission properties of normal radio pulsars and a recently proposed model which takes radio emission from magnetars into consideration.     

MHD waves at a spherical interface modelling EIT waves

Eruptive events such as flares and coronal mass ejections (CMEs) are known to generate global waves propagating over distances comparable to the solar radius in different layers of the solar atmosphere. Here we investigate the propagation of coronal EIT waves, modelled as fast magnetoacoustic modes propagating at a spherical interface in the presence of a purely radial magnetic field. Based on a simplified equilibrium we derive the dispersion relation of the waves. The generation and propagation of EIT waves at the spherical interface is studied numerically for different values of spherical degree and preliminary conclusions are reached regarding the properties of EIT waves. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Properties of spherical galaxies and clusters with an NFW density profile

Using the standard dynamical theory of spherical systems, we calculate the properties of spherical galaxies and clusters whose density profiles obey the universal form first obtained in high-resolution cosmological N -body simulations by Navarro, Frenk & White (NFW). We adopt three models for the internal kinematics: isotropic velocities, constant anisotropy and increasingly radial OsipkovMerritt anisotropy. Analytical solutions are found for the radial dependence of the mass, gravitational potential, velocity dispersion, energy and virial ratio and we test their variability with the concentration parameter describing the density profile and amount of velocity anisotropy. We also compute structural parameters, such as half-mass radius, effective radius and various measures of concentration. Finally, we derive projected quantities, the surface mass density and line-of-sight as well as aperture-velocity dispersion, all of which can be directly applied in observational tests of current scenarios of structure formation. On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW surface density profile is found to fit HubbleReynolds laws well. It is also well fitted by Sérsic R ^{1/ m} laws, for     but in a much narrower range of m and with much larger effective radii than are observed. Assuming in turn reasonable values of the effective radius, the mass density profiles imply a mass-to-light ratio that increases outwards at all radii.     

Toroidal versus poloidal magnetic fields in Sun-like stars: a rotation threshold

From a set of stellar spectropolarimetric observations, we report the detection of surface magnetic fields in a sample of four solar-type stars, namely HD 73350, HD 76151, HD 146233 (18 Sco) and HD 190771. Assuming that the observed variability of polarimetric signal is controlled by stellar rotation, we establish the rotation periods of our targets, with values ranging from 8.8 d (for HD 190771) to 22.7 d (for HD 146233). Apart from rotation, fundamental parameters of the selected objects are very close to the Sun's, making this sample a practical basis to investigate the specific impact of rotation on magnetic properties of Sun-like stars.
We reconstruct the large-scale magnetic geometry of the targets as a low-order  (ℓ < 10)  spherical harmonic expansion of the surface magnetic field. From the set of magnetic maps, we draw two main conclusions. (i) The magnetic energy of the large-scale field increases with rotation rate. The increase in chromospheric emission with the mean magnetic field is flatter than observed in the Sun. Since the chromospheric flux is also sensitive to magnetic elements smaller than those contributing to the polarimetric signal, this observation suggests that a larger fraction of the surface magnetic energy is stored in large scales as rotation increases. (ii) Whereas the magnetic field is mostly poloidal for low rotation rates, more rapid rotators host a large-scale toroidal component in their surface field. From our observations, we infer that a rotation period lower than ≈12 d is necessary for the toroidal magnetic energy to dominate over the poloidal component.     

The haloes of merger remnants

We perform collisionless N -body simulations of 1:1 galaxy mergers, using models which include a galaxy halo, disc and bulge, focusing on the behaviour of the halo component. The galaxy models are constructed without recourse to a Maxwellian approximation. We investigate the effect of varying the galaxies' orientation, their mutual orbit and the initial velocity anisotropy or cusp strength of the haloes upon the remnant halo density profiles and shape, as well as on the kinematics. We observe that the halo density profile (determined as a spherical average, an approximation we find appropriate) is exceptionally robust in mergers, and that the velocity anisotropy of our remnant haloes is nearly independent of the orbits or initial anisotropy of the haloes. The remnants follow the halo anisotropy – local density slope (β–γ) relation suggested by Hansen & Moore in the inner parts of the halo, but β is systematically lower than this relation predicts in the outer parts. Remnant halo axis ratios are strongly dependent on the initial parameters of the haloes and on their orbits. We also find that the remnant haloes are significantly less spherical than those described in studies of simulations which include gas cooling.     

Multiperiodic semiregular variable stars in the ASAS data base: A pilot study

Based on a search for multi‐periodic variability among the semi‐regular red variable stars in the database of the All Sky Automated Survey (ASAS), a sample of 72 typical examples is presented. Their period analysis was performed using the Discrete Fourier Transform. In 41 stars we identified two significant periods each, simultaneously present, while the remaining 31 cases revealed even three such periods per star. They occur in a range roughly between 50 and 3000 days. Inter‐relationships between these periods were analyzed using the “double period diagram” which compares adjacent periods, and the so‐called “Petersen diagram”, the period ratio vs. the shorter period. In both diagrams we could identify six sequences of accumulation of the period values. For five of these sequences (containing 97 % of all data points) we found an almost perfect coincidence with those of previous studies which were based on very different samples of semiregular red variables. Therefore, existence and locations of these sequences in the diagrams seem to be universal features, which appear in any data set of semi‐regularly variable red giants of the AGB; we conclude that they are caused by different pulsation modes as the typical and consistent properties of similar stellar AGB configurations. Stellar pulsations can be considered as the principal cause of the observed periodic variability of these stars, and not binary, rotation of a spotted surface or other possible reasons suggested in the literature. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-Gaussianity detections in the Bianchi VIIh corrected WMAP one-year data made with directional spherical wavelets

Many of the current anomalies reported in the WMAP text ( WMAP ) one-year data disappear after 'correcting' for the best-fitting embedded Bianchi type VII _h component, albeit assuming no dark energy component. We investigate the effect of this Bianchi correction on the detections of non-Gaussianity in the WMAP data that we previously made using directional spherical wavelets. We confirm that the deviations from Gaussianity in the kurtosis of spherical Mexican hat wavelet coefficients are eliminated once the data are corrected for the Bianchi component, as previously discovered by Jaffe et al. This is due to the reduction of the cold spot at Galactic coordinates  ( l , b ) = (209°, −57°)  , which Cruz et al. claimed to be the sole source of non-Gaussianity introduced in the kurtosis. Our previous detections of non-Gaussianity observed in the skewness of spherical wavelet coefficients are not reduced by the Bianchi correction. Indeed, the most significant detection of non-Gaussianity made with the spherical real Morlet wavelet at a significant level of 98.4 per cent remains (using a very conservative method to estimate the significance). Furthermore, we perform preliminary tests to determine if foregrounds or systematics are the source of this non-Gaussian signal, concluding that it is unlikely that these factors are responsible. We make our code to simulate Bianchi-induced temperature fluctuations publicly available.     

Analytical galactic models with mild stellar cusps

In the past two decades, it has been established by high-resolution observations of early-type galaxies that their nuclear surface brightness and corresponding stellar mass densities are characterized by cusps. In this paper, we present a new spherical analytical model family describing mild cuspy centres. We study isotropic and anisotropic models of Osipkov–Merritt type. It is shown that the associated distribution functions and intrinsic velocity dispersions can be represented analytically in a unified way in terms of hypergeometric series, allowing thus a straightforward comparison of these important global quantities for galaxies having underlying mass densities which may differ significantly in their degree of central cuspiness or radial falloff.     

Optical spectroscopic classification and membership of young M dwarfs in star-forming regions

The spectral type is a key parameter in calibrating the temperature which is required to estimate the mass of young stars and brown dwarfs. We describe an approach developed to classify low-mass stars and brown dwarfs in the Trapezium Cluster using red optical spectra, which can be applied to other star-forming regions. The classification uses two methods for greater accuracy: the use of narrow-band spectral indices which rely on the variation of the strength of molecular lines with spectral type and a comparison with other previously classified young, low-mass objects in the Chamaeleon I star-forming region. We have investigated and compared many different molecular indices and have identified a small number of indices which work well for classifying M-type objects in nebular regions. The indices are calibrated for young, pre-main-sequence objects whose spectra are affected by their lower surface gravities compared with those on the main sequence. Spectral types obtained are essentially independent of both reddening and nebular emission lines.
Confirmation of candidate young stars and brown dwarfs as bona fide cluster members may be accomplished with moderate resolution spectra in the optical region by an analysis of the strength of the gravity-sensitive Na doublet. It has been established that this feature is much weaker in these very young objects than in field dwarfs. A sodium spectral index is used to estimate the surface gravity and to demonstrate quantitatively the difference between young (1–2 Myr) objects, and dwarf and giant field stars.     

Detecting Sunyaev–Zel'dovich clusters with Planck– II. Foreground components and optimized filtering schemes

The Planck mission is the most sensitive all-sky cosmic microwave background (CMB) experiment currently planned. The High-Frequency Instrument (HFI) will be especially suited for observing clusters of galaxies by their thermal Sunyaev–Zel'dovich (SZ) effect. In order to assess Planck 's SZ capabilities in the presence of spurious signals, a simulation is presented that combines maps of the thermal and kinetic SZ effects with a realization of the CMB, in addition to Galactic foregrounds (synchrotron emission, free–free emission, thermal emission from dust, CO-line radiation) as well as the submillimetric emission from celestial bodies of our Solar system. Additionally, observational issues such as the finite angular resolution and spatially non-uniform instrumental noise of Planck 's sky maps are taken into account, yielding a set of all-sky flux maps, the autocorrelation and cross-correlation properties of which are examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive and matched filtering are extended to spherical data sets, that enable the amplification of the weak SZ signal in the presence of all contaminations stated above. The theory of scale-adaptive and matched filtering in the framework of spherical maps is developed, the resulting filter kernel shapes are discussed and their functionality is verified.     

Magnetized, mass-loaded, rotating accretion flows

We present a semi-analytical investigation of a simple one-dimensional, steady-state model for a mass-loaded, rotating, magnetized, hydrodynamical flow. Our approach is analogous to one used in early studies of magnetized winds. The model represents the infall towards a central point mass of the gas generated in a cluster of stars surrounding it, as is likely to occur in some active nuclei and starburst galaxies. We describe the properties of the different classes of infall solutions. We find that the flow becomes faster than the fast-mode speed, and hence decoupled from the centre, only for a limited range of parameter values, and when magnetic stresses are ineffective. Such flow is slowed as it approaches a centrifugal barrier, implying the existence of an accretion disc. When the flow does not become super-fast and the magnetic torque is insufficient, no steady solution extending inward to the centre exists. Finally, with a larger magnetic torque, solutions representing steady sub-Alfvénic flows are found, which can resemble spherical hydrodynamical infall. Such solutions, if applicable, would imply that rotation is not important and that any accretion disc formed would be of very limited size.     

Harmonic analysis of cosmic microwave background data – I. Ring reductions and point-source catalogue

We present a harmonic model for the data analysis of an all-sky cosmic microwave background survey, such as Planck , where the survey is obtained through ring-scans of the sky. In this model, resampling and pixelization of the data are avoided. The spherical transforms of the sky at each frequency, in total intensity and polarization, as well as the bright-point-source catalogue, are derived directly from the data reduced on to the rings. Formal errors and the most significant correlation coefficients for the spherical transforms of the frequency maps are preserved. A clean and transparent path from the original samplings in the time domain to the final scientific products is thus obtained. The data analysis is largely based on Fourier analysis of rings; the positional stability of the instrument's spin axis during these scans is a requirement for the data model and is investigated here for the Planck satellite. Brighter point sources are recognized and extracted as part of the ring reductions and, on the basis of accumulated data, used to build the bright-point-source catalogue. The analysis of the rings is performed in an iterative loop, involving a range of geometric and detector response calibrations. The geometric calibrations are used to reconstruct the paths of the detectors over the sky during a scan and the phase offsets between scans of different detectors; the response calibrations eliminate short- and long-term variations in detector response. Point-source information may allow the reconstruction of the beam profile. The reconstructed spherical transforms of the sky in each frequency channel form the input to the subsequent analysis stages. Although the methods in this paper were developed with the data processing for the Planck satellite in mind, there are many aspects which have wider implementation possibilities, including the construction of real-space pixelized maps.     

Bursting and quenching in massive galaxies without major mergers or AGNs

In this work, we study the formation and evolution of dark matter haloes by means of the spherical infalling model with shell-crossing. We present a framework to tackle this effect properly based on the numerical follow-up, with time, of that individual shell of matter that contains always the same fraction of mass with respect to the total mass. In this first step, we do not include angular momentum, velocity dispersion or triaxiality. Within this framework – named as the spherical shell tracker – we investigate the dependence of the evolution of the halo with virial mass, with the adopted mass fraction of the shell, and for different cosmologies. We find that our results are very sensitive to a variation of the halo virial mass or the mass fraction of the shell that we consider. However, we obtain a negligible dependence on cosmology. Furthermore, we show that the effect of shell-crossing plays a crucial role in the way that the halo reaches the stabilization in radius and the virial equilibrium. We find that the values currently adopted in the literature for the actual density contrast at the moment of virialization,  δ_vir  , may not be accurate enough. In this context, we stress the problems related to the definition of a virial mass and a virial radius for the halo. The question of whether the results found here may be obtained by tracking the shells with an analytic approximation remains to be explored.     

Dark matter halo creation in moving barrier models

In hierarchical models of structure formation, the time derivative of the halo mass function may be thought of as the difference of two terms – a creation term, which describes the increase in the number of haloes of mass m from mergers of less massive objects, and a destruction term, which describes the decrease in the number of m -haloes as these merge with other haloes, creating more massive haloes as a result. The first part of this paper focuses on estimating the distribution of times when these creation events take place. In models where haloes form from a spherical collapse, this distribution can be estimated from the same formalism which is used to estimate halo abundances: the constant-barrier excursion-set approach. In the excursion-set approach, moving rather than constant barriers are necessary for estimating halo abundances when the collapse is triaxial. First, we generalize the excursion-set estimate of the creation time distribution by incorporating ellipsoidal collapse. Then, we show that these moving barrier based predictions are in better agreement with measurements in numerical simulations than are the corresponding predictions of the spherical collapse model. In the second part of the paper, we link the creation time distribution to the creation term mentioned above. For this quantity, the improvement provided by the ellipsoidal collapse model is more evident. These results should be useful for studies of merger-driven star formation rates and active galactic nucleus activity. We also present a similar study of the creation of haloes conditioned on belonging to an object of a certain mass today, and reach similar conclusions – the moving barrier based estimates are in substantially better agreement with the simulations. This part of the study may be useful for understanding the tendency for the oldest stars to exist in the most massive objects, and for star formation to only occur in lower mass objects at late times.