Inclusion of sdBs in evolutionary population synthesis for binary stellar populations and the application: the determination of photo-z and galaxy morphology

Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t～1 Gyr, and have been even included in the evolutionary population synthesis (EPS) models by Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007) and those of the Yunnan group (Zhang et al. in Astron. Astrophys. 415:117, 2004; Mon. Not. R. Astron. Soc. 357:1088, 2005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is available upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0–Sd and Irr types of galaxies by using the package of Bruzual and Charlot (Mon. Not. R. Astron. Soc. 344:1000, 2003). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (Astron. Astrophys. 363:476, 2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (Astron. J. 134:579, 2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (Astron. J. 134:579, 2007), we find that the photo-z fluctuate with the SDSS spectroscopic redshifts, while the Sa–Sc galaxies in the catalogue of Fukugita et al. (Astron. J. 134:579, 2007) are classified earlier as Burst-E galaxies.     

Absolute magnitude calibration for red clump stars

We combined the (K _s , J?K _s ) data in Laney et al. (Mon. Not. R. Astron. Soc. 419:1637, 2012) with the V apparent magnitudes and trigonometric parallaxes taken from the Hipparcos catalogue and used them to fit the $M_{K_{s}}$ absolute magnitude to a linear polynomial in terms of V?K _s colour. The mean and standard deviation of the absolute magnitude residuals, ?0.001 and 0.195 mag, respectively, estimated for 224 red clump stars in Laney et al. (2012) are (absolutely) smaller than the corresponding ones estimated by the procedure which adopts a mean $M_{K_{s}}=-1.613~\mbox{mag}$ absolute magnitude for all red clump stars, ?0.053 and 0.218 mag, respectively. The statistics estimated by applying the linear equation to the data of 282 red clump stars in Alves (Astrophys. J. 539:732, 2000) are larger, $\Delta M_{K_{s}}=0.209$ and σ=0.524 mag, which can be explained by a different absolute magnitude trend, i.e. condensation along a horizontal distribution.     

Very Fast and Accurate Azimuth Disambiguation of Vector Magnetograms

We present a method for fast and accurate azimuth disambiguation of vector magnetogram data regardless of the location of the analyzed region on the solar disk. The direction of the transverse field is determined with the principle of minimum deviation of the field from the reference (potential) field. The new disambiguation (NDA) code is examined on the well-known models of Metcalf et al. (Solar Phys. 237, 267, 2006) and Leka et al. (Solar Phys. 260, 83, 2009), and on an artificial model based on the observed magnetic field of AR 10930 (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We compare Hinode/SOT-SP vector magnetograms of AR 10930 disambiguated with three codes: the NDA code, the nonpotential magnetic-field calculation (NPFC: Georgoulis, Astrophys. J. Lett. 629, L69, 2005), and the spherical minimum-energy method (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We then illustrate the performance of NDA on SDO/HMI full-disk magnetic-field observations. We show that our new algorithm is more than four times faster than the fastest algorithm that provides the disambiguation with a satisfactory accuracy (NPFC). At the same time, its accuracy is similar to that of the minimum-energy method (a very slow algorithm). In contrast to other codes, the NDA code maintains high accuracy when the region to be analyzed is very close to the limb.     

The two-body interaction potential in the STF tensor formalism: an application to binary asteroids

The symmetric trace free (STF) tensor formalism, developed by Hartmann et al. (Celest Mech Dyn Astron 60:139–159. doi:10.1007/BF00693097, 1994), is a nice tool, not much used in Celestial Mechanics. It is fully equivalent to the usual spherical harmonics but permits more elegant and compact formulations. The coupling between the gravitational fields of extended bodies with this formalism has been used in Mathis and Le Poncin-Lafitte (Astron Astrophys 497:889–910. doi:10.1051/0004-6361/20079054, 2009) for binary stars or planetary systems, but not yet applied to binary asteroids. However, binary asteroids are common in the Solar System and usually their study requires a full two rigid body approach. The formulation of the two-body interaction potential in the STF formalism in the full two rigid body problem is detailed and completed in this article. An application to the binary asteroid (66391) 1999 KW4 is presented with a comparison of our results with other results of the literature for validation.     

Spectral behavior of AM Her and QQ Vul in high and intermediate states in the UV

We present low resolution UV spectra of two polar systems, AM Her and QQ Vul from the observations taken by the IUE (International Ultraviolet Explorer) of the period between 1978–1996 and 1983–1996 for both systems respectively, to accomplish a large scale study of what happens to the ultraviolet flux of C IV 1550 Å spectral line during different orbital phases. Two spectra for both systems showing the variations in line fluxes and line widths at different orbital phases in high and intermediate states are presented. We concentrated on calculating the line fluxes and line widths of C IV 1550 Å emission line originating in the accretion stream. Our results show that there is spectral variability for the aformentioned physical parameters at different times, similar to that known for the light curve (Heise and Verbunt, Astron. Astrophys. 189:112, 1988; Gansicke et al., Astron. Astrophys. 303:127, 1995; Kafka and Honeycutt, Astron. J. 125:2188K, 2003). We attribute it to the variations of both density and temperature as a result of changing the mass transfer rate (Hutchings et al., Astron. J. 123:2841H, 2002; King and Lasota, Astron. Astrophys. 140L:16K, 1984) which is responsible for this spectral variability. Also we found that the line fluxes of AM Her are greater than the line fluxes of QQ Vul, while the line widths of both systems are approximately the same.     

Nonsingular charged analogues of Schwarzschild’s interior solution

The problem of finding nonsingular charged analogue of Schwarzschild’s interior solutions has been reduced to that of finding a monotonically decreasing function f. The models are discussed in generality by imposing reality condition on f. It is shown that the physical solutions are possible only for surface density to central density ratio greater than or equal to 2/3 i.e. $\frac{\rho_{a}}{\rho_{0}}\ge2/3$ . The unphysical nature of solutions with linear equation state has been proved. A generalization procedure has been utilized to generalize solutions by Guilfoyle (1999). Recently found solutions by Gupta and Kumar (2005a, 2005b, 2005c) are generalized by taking particular form of f and seen to have higher mass and more stable. The maximum mass is found to be 1.59482 M _Θ . The models have been found to be stable once the physical requirements are established due to mass to radius less than 4/9, total charge to total mass ratio less than 1 and redshift quite low.     

The masses of sdB and sdO stars

Mass is a fundamental parameter, but the masses are not well known for most hot subdwarfs. We propose a method of determining the masses of hot subdwarfs. Using this method, we studied the masses of hot subdwarfs from the ESO supernova Ia progenitor survey and the Hamburg quasar survey. The study shows that most of the subdwarf B stars have masses between 0.42 and 0.54 M _⊙, whilst most sdO stars are in the range 0.40～0.55 M _⊙. Comparing our study to the theoretical mass distributions of Han et al. (Mon. Not. R. Astron. Soc. 341:669, 2003), we found that sdO stars with mass less than ～0.5 M _⊙ may evolve from sdB stars, whilst most high-mass (>0.5 M _⊙) sdO stars result from mergers directly.     

On the Use of Roche Equipotentials in Analysing the Problems of Binary and Rotating Stars

Kopal (Adv. Astron. Astrophys., 9, 1, 1972) introduced the concept of Roche equipotentials to analyse the effects of rotational and tidal distortions in case of stars in binary systems. In this approach a mathematical expression for the potential of a star in a binary system is obtained by approximating its inner structure with Roche model. This expression for the potential has been used in subsequent analysis by various authors to analyse the problems of structures and oscillations of synchronous and nonsynchronus binary stars as well as single rotating stars. Occasionally, doubts have been expressed regarding the validity of the use of this approach for analysing nonsynchronous binaries and rotationally and tidally distorted single stars. In this paper we have tried to clarify these doubts.     

Global spectral energy distributions of the Large Magellanic Cloud with interstellar dust

The effects of dust on infrared emission vary among galaxies of different morphological types. We investigated integrated spectral energy distributions (SEDs) in infrared and submillimeter/millimeter emissions from the Large Magellanic Cloud (LMC) based on observations from the Herschel Space Observatory (HSO) and near- to mid-infrared observations from the Spitzer Space Telescope (SST). We also used IRAS and WMAP observations to constrain the SEDs and present the results of radiative transfer calculations using the spectrophotometric galaxy model. We explain the observations by using dust models with different grain size distributions in the interstellar medium of the LMC, noting that the LMC has undergone processes that differ from those in the Milky Way. We determined a spectral index and a normalization factor in the range of ?3.5 to ?3.45 with grain radii in the range of 1 nm–300 nm for the silicate grain and 2 nm–1 μm for the graphite grain. The best fit to the observed SED was obtained with a spectral index of ?3.47, similar to the value derived by Piovan et al. (Mon. Not. R. Astron. Soc. 366(3):923, 2006a). The grain size distribution is described using a power law but with a break that is introduced below a _b , where a larger exponent is used. Changing the graphite grain size distribution significantly changed the SED pattern within the observational uncertainties. Based on the SED fits to the observations from submillimeter wavelengths to infrared radiation from the LMC using GRASIL (Silva et al., Astrophys. J. 509(1):103, 1998), we obtained a reasonable set of parameter values in chemical and geometric space together with the grain size distributions (Weingartner and Draine, Astrophys. J. 548(1):296, 2001) and a modified MRN model with the LMC extinction curve (Piovan et al., Mon. Not. R. Astron. Soc. 366(3):923, 2006a). For a given set of parameters including the disc scale height, synthesis of the starlight spectrum, optical depth, escape time scale, dust model, and star formation efficiency, the adopted dust-to-gas ratio for modeling the observed SEDs, ～1/300 (from the literature) yields a reasonable fit to the observed SEDs and similar results with the metallicity of the LMC as those reported in Russell and Bessell (Astrophys. J. Suppl. Ser. 70:865, 1989). The dust-to-gas ratios that are given as the metallicity caused the variation in the model fits. The difference mainly appears at the wavelengths near 100 μm.     

The local standard of rest and the well in the velocity distribution

It is now recognised that the traditional method of calculating the LSR fails. We find an improved estimate of the LSR by making use of the larger and more accurate database provided by XHIP and repeating our preferred analysis from Francis and Anderson (New Astron 14:615–629, 2009a). We confirm an unexpected high value of $U_0$ by calculating the mean for stars with orbits sufficiently inclined to the galactic plane that they do not participate in bulk streaming motions. Our best estimate of the solar motion with respect to the LSR $(U_0, V_0, W_0) = (14.1\, \pm \, 1.1, 14.6\, \pm \, 0.4, 6.9\, \pm \, 0.1)$ km s $^{-1}$ .     

Primordial magnetic fields constrained from CMB anisotropies on dynamo cosmology

Magneto-curvature stresses could deform magnetic field lines giving rise to back reaction and restoring magnetic stresses (Tsagas in Phys. Rev. Lett., 2001). Barrow and Tsagas (Phys. Rev. D, 2008) have shown that in Friedman universe the expansion slows down in its spatial section of negative Riemann curvature. Earlier, Chicone and Latushkin (Proc. Am. Math. Soc. 125(11):3391, 1995) proved that fast dynamos in compact 2D manifold implies negatively constant Riemannian curvature. Here one applies the Barrow-Tsagas ideas to cosmic dynamos of negative curvature. Fast dynamo, covariant stretching of Riemann slices of cosmic Lobachevsky plane is given. Inclusion of advection term on dynamo equations (Clarkson and Marklund in Mon. Not. R. Astron. Soc., 2005) is considered. In advection absence, slow dynamos are also obtained. It is shown the viscous and restoring forces on stretching particles decrease, as magnetic rates increase. From COBE data ( $\frac{{\delta}B}{B}\approx{10^{-5}}$ ), one is able to compute the stretching $\frac{{\delta}V^{y}}{V^{y}}=1.5\frac{{\delta}B}{B}\approx{1.5{\times}10^{-5}}$ . Zeldovich et al. have computed the maximum magnetic growth rate as γ _max ≈8.0×10^?1 t ^?1. From COBE data a lower growth rate as γ _COBE ≈6.0×10^?6 t ^?1, is well-within Zeldovich et al estimate. Instead of Harrison value $B\approx{t^{\frac{4}{3}}}$ one obtains a lower primordial field B≈10^?6 t which yields B≈10^?6 G at 1 s Big Bang time.     

On the Origin of Pulsations of Sub-THz Emission from Solar Flares

We propose a model to explain fast pulsations in sub-THz emission from solar flares. The model is based on the approach of a flaring loop as an equivalent electric circuit and explains the pulse-repetition rate, the high-quality factor, Q≥10³, low modulation depth, pulse synchronism at different frequencies, and the dependence of the pulse-repetition rate on the emission flux, observed by Kaufmann et al. (Astrophys. J. 697, 420, 2009). We solved the nonlinear equation for electric current oscillations using a Van der Pol method and found the steady-state value for the amplitude of the current oscillations. Using the pulse rate variation during the flare on 4 November 2003, we found a decrease of the electric current from 1.7×10¹² A in the flare maximum to 4×10¹⁰ A just after the burst. Our model is consistent with the plasma mechanism of sub-THz emission suggested recently by Zaitsev, Stepanov, and Melnikov (Astron. Lett. 39, 650, 2013).     

Capture probability in the 3:1 mean motion resonance with Jupiter: an application to the Vesta family

We study the capture and crossing probabilities in the 3:1 mean motion resonance with Jupiter for a small asteroid that migrates from the inner to the middle Main Belt under the action of the Yarkovsky effect. We use an algebraic mapping of the averaged planar restricted three-body problem based on the symplectic mapping of Hadjidemetriou (Celest Mech Dyn Astron 56:563–599, 1993), adding the secular variations of the orbit of Jupiter and non-symplectic terms to simulate the migration. We found that, for fast migration rates, the captures occur at discrete windows of initial eccentricities whose specific locations depend on the initial resonant angles, indicating that the capture phenomenon is not probabilistic. For slow migration rates, these windows become narrower and start to accumulate at low eccentricities, generating a region of mutual overlap where the capture probability tends to 100 %, in agreement with the theoretical predictions for the adiabatic regime. Our simulations allow us to predict the capture probabilities in both the adiabatic and non-adiabatic cases, in good agreement with results of Gomes (Celest Mech Dyn Astron 61:97–113, 1995) and Quillen (Mon Not RAS 365:1367–1382, 2006). We apply our model to the case of the Vesta asteroid family in the same context as Roig et al. (Icarus 194:125–136, 2008), and found results indicating that the high capture probability of Vesta family members into the 3:1 mean motion resonance is basically governed by the eccentricity of Jupiter and its secular variations.     

Validation of the 3D AMR SIP–CESE Solar Wind Model for Four Carrington Rotations

We carry out the adaptive mesh refinement (AMR) implementation of our solar–interplanetary space-time conservation element and solution element (CESE) magnetohydrodynamic model (SIP–CESE MHD model) using a six-component grid system (Feng, Zhou, and Wu, Astrophys. J. 655, 1110, 2007; Feng et al., Astrophys. J. 723, 300, 2010). By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation (Jiang et al., Solar Phys. 267, 463, 2010), the SIP–AMR–CESE MHD model is implemented in the reference coordinates with the aid of the parallel AMR package PARAMESH available at http://sourceforge.net/projects/paramesh/ . Meanwhile, the volumetric heating source terms derived from the topology of the magnetic-field expansion factor and the minimum angular separation (at the photosphere) between an open-field foot point and its nearest coronal-hole boundary are also included. We show the preliminary results of applying the SIP–AMR–CESE MHD model for simulating the solar-wind background of different solar-activity phases by comparison with SOHO observations and other spacecraft data from OMNI. Our numerical results show overall good agreements in the solar corona and in interplanetary space with these multiple-spacecraft observations.     

Enhanced luminosity of young stellar objects in cometary globules

In this study we investigated the effects of external trigger on the characteristics of young stellar objects (YSOs) associated with cometary globules (CGs). We made optical spectroscopy of stars associated with star-forming CGs. We find that the masses of the most massive stars associated with CGs are correlated with the masses of the parent cloud but they are systematically larger than expected for clouds of similar mass from the relation M _max-star=0.33M _cl ^0.43 given by Larson (Mon. Not. R. Astron. Soc. 200:159, 1982). We have also estimated the luminosities of the IRAS sources found associated with CGs as a function of cloud mass and then compared them with those of the IRAS sources found associated with isolated opacity class 6 clouds (isolated and relatively away from large star forming regions). We find that the luminosities of IRAS sources associated with CGs are larger than those of the opacity class 6 clouds. These findings support results from recent simulations in which it was shown that the Radiation Driven Implosion (RDI) process, believed to be responsible for the cometary morphology and star formation, can increase the luminosity 1–2 orders of magnitudes higher than those of protostars formed without external triggering due to an increase in accretion rates. Thus implying that the massive stars can have profound influence on the star formation in clouds located in their vicinity.     

On the Use of Meteor Camera Systems in the Detection of Kuiper Belt Objects Through Serendipitous Stellar Occultations

The direct detection of Kuiper Belt Objects (KBOs) by telescopic imaging is not currently practical for objects much less than 100 km in diameter. However, indirect methods such as serendipitous stellar occultations might still be employed to detect these bodies. The method of serendipitous stellar occultations has been previously used with some success in detecting KBOs—Roques et al. (Astron J 132(2):819–822, 2006) detected three Trans-Neptunian objects; Schlichting et al. (Nature 462(7275):895–897, 2009) and Schlichting et al. (Astrophys J 761:150, 2012) each detected a single object in archival Hubble Space Telescope data. However, previous assessments of KBO occultation detection rates have been calculated only for telescopes—we extend this method to video camera systems, and we apply this derivation to the automated meteor camera systems currently in use at the University of Western Ontario. We find that in a typical scenario we can expect one occultation per month. However recent studies such as those of Shankman et al. (Astrophys. J. Lett. 764. doi:10.1088/2041-8205/764/1/L2, 2013) and Gladman et al. (AAS/Division for Planetary Sciences Meeting Abstracts, 2012) which indicate that the population of small KBOs may be smaller than has been assumed in the past may result in a sharp reduction of these rates. Nonetheless, a survey for KBO occultations using existing meteor camera systems may provide valuable information about the number density of KBOs.     

Eigenfrequencies of rotationally and tidally distorted white dwarf models of stars

In the present paper we have studied the eigenfrequencies of small adiabatic barotropic pseudo-radial and nonradial modes of oscillations of the white dwarf models of rotating stars in binary systems. In this work the methodology of Mohan and Saxena (in Astrophys. Space Sci. 113:155, 1985) has been used that utilizes the averaging technique of Kippenhahn and Thomas (in Proc. IAU Colloq., vol. 4, p. 20, 1970) and certain results on Roche equipotential as that given by Kopal (in Advances in Astronomy and Astrophysics, Academic Press, 1972). The objective of this study is to investigate the effects of rotation and/or tidal distortion on the periods of oscillations of rotationally and/or tidally distorted white dwarf models of stars assuming it to be the primary component of the binary system and rotating uniformly. The results of present study show that the eigenfrequencies (both radial and nonradial modes) of the rotationally distorted and rotationally and tidally distorted white dwarf model of stars in binary systems tend to decrease under the influence of rotational distortions and rotational and tidal distortions, respectively. However, results are contrary for tidally distorted white dwarf model of stars.     

Electron capture strength on odd-A nucleus 59Co in explosive astrophysical environment

The Gamow-Teller (GT) transitions within massive stars play sumptuous role in the dynamics of core collapse supernovae. GT strength distributions and electron capture rates have been calculated for odd-A nucleus ⁵⁹Co within the proton-neutron quasiparticles random phase approximation (pn-QRPA) formalism. The pn-QRPA results are compared with other model calculations and (n,p) reaction experiment carried out at TRIUMF charge-exchange facility. The pn-QRPA calculated a total B(GT ₊) strength of 3.3 for ⁵⁹Co to be compared with the shell model value of 2.5 and the 1.9±0.1 in the (n,p) charge-exchange reaction. Aufderheide et al. (1993) extracted total strength equaling 2.4±0.3. The placement of GT centroid at 5.6 MeV by the pn-QRPA model is in reasonable agreement with the shell model centroid at 5.1 MeV whereas the measured GT centroid was placed at 4.4±0.3 MeV in the (n,p) experiment. Fuller, Fowler and Newman (FFN) (1980, 1982a, 1982b), placed the GT centroid at too low excitation energy of 2.0 MeV in the daughter nucleus ⁵⁹Fe, and this misplacement led to the enhancement of FFN rates. The suppressed pn-QRPA and shell model electron capture rates are in good agreement with each other. The rates are suggestive of higher value of Y _e (electron-to-baryon ratio) and may contribute to a more massive homologously collapsing core resulting in a more energetic shock. It might be interesting for the simulators to check the effect of these suppressed rates on the fine-tuning of the time rate of Y _e , the concomitant heavy element nucleosynthesis, and, on the energetics of the subsequent shock wave.     

Modified Coronal Index of the Solar Activity

The original coronal index of the solar activity (CI) has been constructed on the basis of ground-based measurements of the intensities of the coronal line of 530.3 nm (Rybanský in Bull. Astron. Inst. Czechoslov., 28, 367, 1975; Rybanský et al. in J. Geophys. Res., 110, A08106, 2005). In this paper, CI is compared with the EUV measurements on the CELIAS/SEM equipment based on the same idea as the original idea of the coronal index. The correlation is very good for the period 1996?–?2005 (r=0.94 for daily values). The principal result of this paper is the introduction of the modified coronal index (MCI) which in all uses and contexts can replace the existing CI index. Daily MCI values extend over a time period of six solar activity cycles. Future MCI measurements will be derived from more reliable measurements made by space-based observatories that are not influenced by the weather. MCI measurements are and will continue to be archived at the web site of the Slovak Central Observatory in Hurbanovo ( http://www.suh.sk/obs/vysl/MCI.htm ).     

Developing an Optical Lunar Occultation Measurement Reduction System for Observations at Kaau Observatory

We present here a detailed explanation of the reduction method that we use to determine the angular diameters of the stars occulted by the dark limb of the moon. This is a main part of the lunar occultation observation program running at King Abdul Aziz University observatory since late 1993. The process is based on the least square model fitting method of analyzing occultation data, first introduced by Nather et al. (Astron. J. 75:963, 1970).