Old stellar population synthesis： new age and mass estimates for Mayall Ⅱ = G1

Mayall Ⅱ = G1 is one of the most luminous globular clusters （GCs） in M31. Here, we determine its age and mass by comparing multicolor photometry with theoretical stellar population synthesis models. Based on far- and near-ultraviolet GALEX photometry, broad-band UBVRI, and infrared JHKs 2MASS data, we construct the most extensive spectral energy distribution of G 1 to date, spanning the wavelength range from 1538 to 20 000A. A quantitative comparison with a variety of simple stellar population （SSP） models yields a mean age which is consistent with G1 being among the oldest building blocks of M31 and having formed within ～1.7 Gyr after the Big Bang. Irrespective of the SSP model or stellar initial mass function adopted, the resulting mass estimates （of order 10^7M⊙） indicate that GI is one of the most massive GCs in the Local Group. However, we speculate that the cluster＇s exceptionally high mass suggests that it may not be a genuine GC. Our results also suggest that G1 may contain, on average, （1.65±0.63） × 10^2L⊙ far-ultraviolet-bright, hot, extreme horizontal-branch stars, depending on the adopted SSP model. In addition, we demonstrate that extensive multi-passband photometry coupled with SSP analysis enables one to obtain age estimates for old SSPs that have similar accuracies as those from integrated spectroscopy or resolved stellar photometry, provided that some of the free parameters can be constrained independently.     

The environmental dependence of different colors in the CMASS sample of the SDSS DR9

I investigate the environmental dependence of galaxy colors in the CMASS sample of the Sloan Digital Sky Survey Data Release 9(SDSS DR9).To decrease the radial selection effect,I divide the CMASS sample into subsamples with a redshift binning size of?z=0.01 and analyze the environmental dependence of the u-r,u-g,g-r,r-i and i-z colors for these subsamples in each redshift bin.Statistical analysis shows that all five colors weakly correlate with the local environment,which may mean that the environmental processes responsible for a galaxy’s properties proceed slowly over cosmic time.     

Automated estimation of stellar fundamental parameters from low resolution spectra： the PLS method

PLS （Partial Least Squares regression） is introduced into an automatic estimation of fundamental stellar spectral parameters. It extracts the most correlative spectral component to the parameters （Teff, log g and [Fe/H]）, and sets up a linear regression function from spectra to the corresponding parameters. Considering the properties of stellar spectra and the PLS algorithm, we present a piecewise PLS regression method for estimation of stellar parameters, which is composed of one PLS model for Teff, and seven PLS models for log g and [Fe/H] estimation. Its performance is investigated by large experiments on flux calibrated spectra and continuum normalized spectra at different signal-to-noise ratios （SNRs） and resolutions. The results show that the piecewise PLS method is robust for spectra at the medium resolution of 0.23 nm. For low resolution 0.5 nm and 1 nm spectra, it achieves competitive results at higher SNR. Experiments using ELODIE spectra of 0.23 nm resolution illustrate that our piecewise PLS models trained with MILES spectra are efficient for O ～ G stars： for flux calibrated spectra, the systematic offsets are 3.8%, 0.14 dex, and -0.09 dex for Teff, log g and [Fe/H], with error scatters of 5.2%, 0.44 dex and 0.38 dex, respectively; for continuum normalized spectra, the systematic offsets are 3.8%, 0.12dex, and -0.13 dex for Teff, log g and [Fe/H], with error scatters of 5.2%, 0.49 dex and 0.41 dex, respectively. The PLS method is rapid, easy to use and does not rely as strongly on the tightness of a parameter grid of templates to reach high precision as Artificial Neural Networks or minimum distance methods do.     

Forecast daily indices of solar activity, F10.7, using support vector regression method

The 10.7 cm solar radio flux （F10.7）, the value of the solar radio emission flux density at a wavelength of 10.7 cm, is a useful index of solar activity as a proxy for solar extreme ultraviolet radiation. It is meaningful and important to predict F10.7 values accurately for both long-term （months-years） and short-term （days） forecasting, which are often used as inputs in space weather models. This study applies a novel neural network technique, support vector regression （SVR）, to forecasting daily values of F10.7. The aim of this study is to examine the feasibility of SVR in short-term F10.7 forecasting. The approach, based on SVR, reduces the dimension of feature space in the training process by using a kernel-based learning algorithm. Thus, the complexity of the calculation becomes lower and a small amount of training data will be sufficient. The time series of F10.7 from 2002 to 2006 are employed as the data sets. The performance of the approach is estimated by calculating the norm mean square error and mean absolute percentage error. It is shown that our approach can perform well by using fewer training data points than the traditional neural network.     

A logistic model for magnetic energy storage in solar active regions

Previous statistical analyses of a large number of SOHO/MDI full disk longitudinal magnetograms provided a result that demonstrated how responses of solar flares to photospheric magnetic properties can be fitted with sigmoid functions. A logistic model reveals that these fitted sigmoid functions might be related to the free energy storage process in solar active regions. Although this suggested model is rather simple, the free energy level of active regions can be estimated and the probability of a solar flare with importance over a threshold can be forecast within a given time window.     

A modified TreePM code

We discuss the performance characteristics of using the modification of the tree code suggested by Barnes in the context of the TreePM code. The optimization involves identifying groups of particles and using only one tree walk to compute the force for all the particles in the group. This modification has been in use in our implementation of the TreePM code for some time, and has also been used by others in codes that make use of tree structures. We present the first detailed study of the performance characteristics of this optimization. We show that the modification, if tuned properly, can speed up the TreePM code by a significant amount. We also combine this modification with the use of individual time steps and indicate how to combine these two schemes in an optimal fashion. We find that the combination is at least a factor of two faster than the modified TreePM without individual time steps. Overall performance is often faster by a larger factor because the scheme for the groups optimizes the use of cache for large simulations.     

The dynamical architecture and habitable zones of the quintuplet planetary system 55 Cancri

We perform numerical simulations to study the secular orbital evolution and dynamical structure of the quintuplet planetary system 55 Cancri with the self-consistent orbital solutions by Fischer and coworkers. In the simulations, we show that this system can be stable for at least 108 yr. In addition, we extensively investigate the planetary configuration of four outer companions with one terrestrial planet in the wide region of 0.790AU 〈 a 〈 5.900AU to examine the existence of potential asteroid structure and Habitable Zones （HZs）. We show that there are unstable regions for orbits about 4：1, 3：1 and 5：2 mean motion resonances （MMRs） of the outermost planet in the system, and several stable orbits can remain at 3：2 and 1：1 MMRs, which resembles the asteroid belt in the solar system. From a dynamical viewpoint, proper HZ candidates for the existence of more potential terrestrial planets reside in the wide area between 1.0 AU and 2.3 AU with relatively low eccentricities.     

Temporal variations and spectral properties of the Be/X-ray pulsar GRO J1008–57 studied by INTEGRAL

The spin period variations and hard X-ray spectral properties of the Be/Xray pulsar GRO J1008–57 are studied with INTEGRAL observations during two outbursts in 2004 June and 2009 March.The pulsation periods of～93.66 s in 2004and～93.73 s in 2009 are determined.Pulse profiles of GRO J1008–57 during outbursts are strongly energy dependent with a double-peaked profile from 3–7 keV and a single-peaked profile in hard X-rays above 7 keV.Combined with previous measurements,we find that GRO J1008–57 has undergone a spin-down trend from 1993–2009 with a rate of～4.1×10-5s d-1,and could have changed into a spin-up trend after 2009.We find a relatively soft spectrum in the early phase of the 2009 outburst with cutoff energy～13 keV.Above a hard X-ray flux of～10-9erg cm-2s-1,the spectra of GRO J1008–57 during outbursts need an enhanced hydrogen absorption with column density～6×1022cm-2.The observed dip-like pulse profile of GRO J1008–57 in soft X-ray bands could be caused by this intrinsic absorption.Around the outburst peaks,a possible cyclotron resonance scattering feature at～74 keV is detected in the spectra of GRO J1008–57 which is consistent with the feature that was reported in MAXI/GSC observations,making the source a neutron star with the highest known magnetic field(～6.6×1012G)among accreting X-ray pulsars.This marginal feature is supported by the present detections in GRO J1008–57 following the correlation between the fundamental line energies and cutoff energies in accreting X-ray pulsars.Finally we discovered two modulation periods at～124.38 d and～248.78 d using RXTE/ASM light curves of GRO J1008–57.Two flare peaks appearing in the folded light curve had different spectral properties.The normal outburst lasting 0.1 of an orbital phase had a hard spectrum and could not be significantly detected below 3 keV.The second flare lasting ten days showed a very soft spectrum without significant detections above 5 keV.GRO J1008–57 is a good candidate of an accreting system with an equatorial circumstellar disk around the companion star.The neutron star passing the disk of the Be star near periastron and apastron produces two X-ray flares.The soft spectral properties in the secondary flares still need further detailed studies with soft X-ray spectroscopy.     

The spatial distribution of dark matter annihilation originating from a gamma-ray line signal

The GeV-TeVγ-ray line signal is the smoking gun signature of dark matter annihilation or decay.The detection of such a signal is one of the main targets of some space-based telescopes,including Fermi-LAT and the upcoming missions CALET,DAMPE and Gamma-400.An important feature ofγ-ray line photons that originate from dark-matter-annihilation is that they are concentrated at the center of the Galaxy.So far,no reliableγ-ray line has been detected by Fermi-LAT,and the upper limits on the cross section of annihilation intoγ-rays have been reported.We use these upper limits to estimate the"maximal"number ofγ-ray line photons detectable for FermiLAT,DAMPE and Gamma-400,and then investigate the spatial distribution of these photons.We show that the center of the distribution will usually be offset from the Galactic center(Sgr A)due to the limited statistics.Such a result is almost independent of models of the dark matter distribution,and will render the reconstruction of the dark matter distribution with theγ-ray line signal very challenging for foreseeable space-based detectors.     

On properties of mass transfer in an Algol-type binary system： an example of VV Vir

Some properties of mass transfer within a semi-detached binary system are discussed based on an example of VV Vir. The observations and analysis of VV Vir show that it is a semi-detached binary system with the less massive cooler component filling its Roche lobe, which is also called an Algol-type binary system. Based on the parameters of this semi-detached binary, a theoretical study on the mass flow is carded out, including the trajectory of the mass flow, the position, radius and temperature of the impact spot caused by the mass flow, the inconsistent form of the mass flow and the possibly huge rate of energy transfer carried by the mass flow. Humps and distortions in light curves of VV Vir are deemed to be weak evidence for the theoretical results.     

New upper limit on the cosmological constant from solar system dynamics

The cosmological constantΛis the simplest model for explaining the dark energy which supposedly drives the observed accelerated expansion rate of the Universe.Together with the concept of cold dark matter,it satisfactorily accommodates a wealth of observations related to cosmology.Due to its assumed constancy throughout the Universe,Λmight also affect the dynamics of the planets in the solar system,although with extremely small effects.However,modern high-precision ephemerides provide a promising tool for constraining it.Using the supplementary advances in the perihelia provided by current INPOP10a and EPM2011 ephemerides,we obtain a new upper limit onΛin the solar system when the Lense-Thirring effect due to the Sun’s angular momentum and the uncertainty of the Sun’s quadrupole moment are properly taken into account.These two factors were mostly absent in previous works dealing withΛ.We find that INPOP10a yields an upper limit ofΛ=(0.26±1.45)×10-43m-2and EPM2011 givesΛ=(-0.44±8.93)×10-43m-2.Such bounds are about 10 times less than previously estimated results.     

Foreground removal of 21 cm fluctuation with multifrequency fitting

The 21 centimeter （21 cm） line emission from neutral hydrogen in the intergalactic medium （IGM） at high redshifts is strongly contaminated by foreground sources such as the diffuse Galactic synchrotron emission and free-free emission from the Galaxy, as well as emission from extragalactic radio sources, thus making its observation very complicated. However, the 21 cm signal can be recovered through its structure in frequency space, as the power spectrum of the foreground contamination is expected to be smooth over a wide band in frequency space while the 21 cm fluctuations vary significantly. We use a simple polynomial fitting to reconstruct the 21 cm signal around four frequencies 50, 100, 150 and 200 MHz with an especially small channel width of 20 kHz. Our calculations show that this multifrequency fitting approach can effectively recover the 21 cm signal in the frequency range 100 - 200 MHz. However, this method doesn＇t work well around 50 MHz because of the low intensity of the 21 cm signal at this frequency. We also show that the fluctuation of detector noise can be suppressed to a very low level by taking long integration times, which means that we can reach a sensitivity of ≈ 10 mK at 150 MHz with 40 antennas in 120 hours of observations.     

Ionization state of cosmic hydrogen by early stars and quasars

Cosmic hydrogen is reionized and maintained in its highly ionized state by the ultraviolet emission attributed to an early generation of stars and quasars. The Lyα opacity observed in absorption spectra of high-redshift quasars permits more stringent constraints on the ionization state of cosmic hydrogen. Based on density perturbation and structure formation theory, we develop an analytic model to trace the evolution of the ionization state in the post-overlap epoch of reionization, in which the bias factor is taken into account. For quasars, we represent an improved luminosity function by utilizing a hybrid approach for the halo formation rate that is in reasonable agreement with the published measurements at 2 ≤ z ≤ 6. Comparison with the classic Press-Schechter mass function of dark matter halos, we demonstrate that the biased mass distribution indeed enhances star formation efficiency in the overdense environment by more than 25 per cent following the overlap of ionized bubbles. In addition, an alternative way is introduced to derive robust estimates of the mean free path for ionizing photons. In our model, star-forming galaxies are likely to dominate the ionizing background radiation beyond z = 3, and quasars contribute equally above a redshift of z ～ 2.5. From 5 ≤ z ≤ 6, the lack of evolution in photoionization rate can thus be explained by the relatively flat evolution in star formation efficiency, although the mean free path of ionizing photons increases rapidly. Moreover, in the redshift interval z ～ 2 - 6, the expected mean free path and Gunn-Peterson optical depth obviously evolve by a factor of ～ 500 and ～50 respectively. We find that the relative values of critical overdensities for hydrogen ionization and collapse could be 430% at z ≈2 and 2% at z ≈6, suggesting a rapid overlap process in the overdense regions around instant quasars following reionization. We further illustrate that the absolute estimates of the fraction of neutral hydrogen computed from theoretical models may n     

Relations between stellar mass and electron temperature-based metallicity for star-forming galaxies in a wide mass range

We select 947 star-forming galaxies from SDSS-DR7 with [O III]λ4363emission lines detected at a signal-to-noise ratio larger than 5σ. Their electron temperatures and direct oxygen abundances are then determined. We compare the results from different methods. t2, the electron temperature in the low ionization region, estimated from t3, that in the high ionization region, is compared using three analysis relations between t2- t3. These show obvious differences, which result in some different ionic oxygen abundances. The results of t3, t2, O++/H+and O+/H+derived by using methods from IRAF and literature are also compared. The ionic abundances O++/H+are higher than O+/H+for most cases. The different oxygen abundances derived from Teand the strong-line ratios show a clear discrepancy, which is more obvious following increasing stellar mass and strong-line ratio R23. The sample of galaxies from SDSS with detected [O III]λ4363 have lower metallicites and higher star formation rates, so they may not be typical representatives of the whole population of galaxies. Adopting data objects from Andrews Martini, Liang et al. and Lee et al. data, we derive new relations of stellar mass and metallicity for star-forming galaxies in a much wider stellar mass range: from 106 M to 1011 M.     

What does the Sun tell and hint now？

Some historical records, which have held since the beginning of modern solar activity cycles, are being broken by the present Sun： cycle 23 records the longest cycle length and fall time; latitudes of high-latitude sunspots belonging to a new cycle around the minimum time of the cycle are statistically the lowest at present, compared with those of other cycles; there are only one or no sunspots in a month appearing at high latitudes for 58 months, which is the first time that such a long duration has been observed. The solar dynamo is believed to be slowing down due to： （1） the minimum smoothed monthly mean sunspot number is the smallest since cycle 16 onwards, and even probably among all modern solar cycles; and （2） once the time interval between the first observations of two neighboring sunspot groups is larger than 14 d, it should be approximately regarded as an observation of no sunspots on the visible solar disk, called a spotless event. Spotless events occur with the highest frequency around the minimum time of cycle 24, and the longest spotless event also appears around the minimum time for observations of the Sun since cycle 16. Cycle 24 is expected to have the lowest level of sunspot activity from cycle 16 onwards and even probably for all of the modern solar cycles.     

Equilibria of a charged artificial satellite subject to gravitational and Lorentz torques

The attitude dynamics of a rigid artificial satellite subject to a gravity gradient and Lorentz torques in a circular orbit are considered. Lorentz torque is developed on the basis of the electrodynamic effects of the Lorentz force acting on the charged satellite's surface. We assume that the satellite is moving in a Low Earth Orbit in the geomagnetic field, which is considered to be a dipole. Our model of torque due to the Lorentz force is developed for an artificial satellite with a general shape, and the nonlinear differential equations of Euler are used to describe its attitude orientation. All equilibrium positions are determined and conditions for their existence are obtained.The numerical results show that the charge q and radius ρ0of the center of charge for the satellite provide a certain type of semi-passive control for the attitude of the satellite. The technique for this kind of control would be to increase or decrease the electrostatic screening on the satellite. The results obtained confirm that the change in charge can affect the magnitude of the Lorentz torque, which can also affect control of the satellite. Moreover, the relationship between magnitude of the Lorentz torque and inclination of the orbit is investigated.     

A review of solar type III radio bursts

Solar type III radio bursts are an important diagnostic tool in the understanding of solar accelerated electron beams. They are a signature of propagating beams of nonthermal electrons in the solar atmosphere and the solar system. Consequently, they provide information on electron acceleration and transport, and the conditions of the background ambient plasma they travel through. We review the observational properties of type III bursts with an emphasis on recent results and how each property can help identify attributes of electron beams and the ambient background plasma. We also review some of the theoretical aspects of type III radio bursts and cover a number of numerical efforts that simulate electron beam transport through the solar corona and the heliosphere.     

Numerical simulations of three-dimensional magnetic swirls in a solar flux-tube

We aim to numerically study evolution of Alfv′en waves that accompany short-lasting swirl events in a solar magnetic flux-tube that can be a simple model of a magnetic pore or a sunspot. With the use of the FLASH code we numerically solve three-dimensional ideal magnetohydrodynamic equations to simulate twists which are implemented at the top of the photosphere in magnetic field lines of the flux-tube. Our numerical results exhibit swirl events and Alfv′en waves with associated clockwise and counterclockwise rotation of magnetic lines, with the largest values of vorticity at the bottom of the chromosphere, and a certain amount of energy flux.     

A 96-antenna radioheliograph

Here we briefly present some design approaches for a multifrequency 96-antenna radioheliograph. The configuration of the array antenna, transmission lines and digital receivers is the main focus of this work. The radioheliograph is a T-shaped centrally condensed radiointerferometer operating in the frequency range 4–8 GHz.The justification for the choice of such a configuration is discussed. The signals from antennas are transmitted to a workroom by analog optical links. The dynamic range and phase errors of the microwave-over-optical signal are considered. The signals after downconverting are processed by digital receivers for delay tracking and fringe stopping. The required step of delay tracking and data rates are considered. Two 3-bit data streams(I and Q) are transmitted to a correlator with the transceivers embedded in Field Programmed Gate Array chips and with PCI Express cables.