An asteroseismic study of the β Cephei star θ Ophiuchi: photometric results

We have carried out a three-site photometric campaign for the β Cephei star θ Oph from 2003 April to August. 245 h of differential photoelectric u v y photometry were obtained during 77 clear nights. The frequency analysis of our measurements has resulted in the detection of seven pulsation modes within a narrow frequency interval between 7.116 and 7.973 c d⁻¹. No combination or harmonic frequencies have been found. We have performed a mode identification of the individual pulsations from our colour photometry that shows the presence of one radial mode, one rotationally split  ℓ= 1  triplet and possibly three components of a rotationally split  ℓ= 2  quintuplet. We discuss the implications of our findings and point out the similarity of the pulsation spectrum of θ Oph to that of another β Cephei star, V836 Cen.     

Particle accretion and dissipation simulator: Collisional aggregation of icy particles

As a part of a project to study the cometary particle environment, a self-consistent ballistic rigid body particle colliding simulator, capable of producing fractal aggregates of icy particles modelled by hard or soft (metaball) spheres has been developed, to be used for creating test particles in future statistical studies of the cometary coma-particle interaction. Instead of a predefined sticking probability, actual surface interactions are modelled and a method for calculating internal stresses and fragmentation of an aggregate during the accretion phase is described. Simulations suggest that taking fragmentation into account has two major consequences, increase in the fractal dimension of resulting particles and a noticeable change in the mass spectrum of an ensemble.     

On the binarity of Herbig Ae/Be stars

We present high-resolution spectro-astrometry of a sample of 28 Herbig Ae/Be and three F-type pre-main-sequence stars. The spectro-astrometry, which is essentially the study of unresolved features in long-slit spectra, is shown from both empirical and simulated data to be capable of detecting binary companions that are fainter by up to 6 mag at separations larger than ∼0.1 arcsec. The nine targets that were previously known to be binary are all detected. In addition, we report the discovery of six new binaries and present five further possible binaries. The resulting binary fraction is 68 ± 11 per cent. This overall binary fraction is the largest reported for any observed sample of Herbig Ae/Be stars, presumably because of the exquisite sensitivity of spectro-astrometry for detecting binary systems. The data hint that the binary frequency of the Herbig Be stars is larger than that of the Herbig Ae stars. The Appendix presents model simulations to assess the capabilities of spectro-astrometry and reinforces the empirical findings. Most spectro-astrometric signatures in this sample of Herbig Ae/Be stars can be explained by the presence of a binary system. Two objects, HD 87643 and Z CMa, display evidence for asymmetric outflows. Finally, the position angles of the binary systems have been compared with available orientations of the circumprimary disc and these appear to be coplanar. The alignment between the circumprimary discs and the binary systems strongly suggests that the formation of binaries with intermediate-mass primaries is due to fragmentation as the alternative, stellar capture, does not naturally predict aligned discs. The alignment extends to the most massive B-type stars in our sample. This leads us to conclude that formation mechanisms that do result in massive stars, but predict random angles between the binaries and the circumprimary discs, such as stellar collisions, are also ruled out for the same reason.     

Simulations of multimode bolometric interferometers

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.     

Research on Parallel Algorithms for uv-faceting Imaging

The uv-faceting imaging is one of the widely used large field of view imaging technologies, and will be adopted for the data processing of the low-frequency array in the first stage of the Square Kilometre Array (SKA1). Due to the scale of the raw data of SKA1 is unprecedentedly large, the efficiency of data processing directly using the original uv-faceting imaging will be very low. Therefore, a uv-faceting imaging algorithm based on the MPI (Message Passing Interface)+OpenMP (Open Multi-Processing) and a uv-faceting imaging algorithm based on the MPI+CUDA (Compute Unified Device Architecture) are proposed. The most time-consuming data reading and gridding in the algorithm are optimized in parallel. The verification results show that the results of the proposed two algorithms are basically consistent with that obtained by the current mainstream data processing software CASA (Common Astronomy Software Applications), which indicates that the proposed two algorithms are basically correct. Further analysis of the accuracy and total running time shows that the MPI+CUDA method is better than the MPI+OpenMP method in both the correctness rate and running speed. The performance test results show that the proposed algorithms are effective and have certain extensibility.     

Study of Star/Galaxy Classification Based on the XGBoost Algorithm

Machine learning has achieved great success in many areas today. The lifting algorithm has a strong ability to adapt to various scenarios with a high accuracy, and has played a great role in many fields. But in astronomy, the application of lifting algorithms is still rare. In response to the low classification accuracy of the dark star/galaxy source set in the Sloan Digital Sky Survey (SDSS), a new research result of machine learning, eXtreme Gradient Boosting (XGBoost), has been introduced. The complete photometric data set is obtained from the SDSS-DR7, and divided into a bright source set and a dark source set according to the star magnitude. Firstly, the ten-fold cross-validation method is used for the bright source set and the dark source set respectively, and the XGBoost algorithm is used to establish the star/galaxy classification model. Then, the grid search and other methods are used to adjust the XGBoost parameters. Finally, based on the galaxy classification accuracy and other indicators, the classification results are analyzed, by comparing with the models of function tree (FT), Adaptive boosting (Adaboost), Random Forest (RF), Gradient Boosting Decision Tree (GBDT), Stacked Denoising AutoEncoders (SDAE), and Deep Belief Nets (DBN). The experimental results show that, the XGBoost improves the classification accuracy of galaxies in the dark source classification by nearly 10% as compared to the function tree algorithm, and improves the classification accuracy of sources with the darkest magnitudes in the dark source set by nearly 5% as compared to the function tree algorithm. Compared with other traditional machine learning algorithms and deep neural networks, the XGBoost also has different degrees of improvement.