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.     

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.     

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.     

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.     

Stability of triangular equilibrium points in the elliptic restricted problem of three bodies with radiating and triaxial primaries

This paper studies the stability of infinitesimal motions about the triangular equilibrium points in the elliptic restricted three body problem assuming bigger primary as a source of radiation and the smaller one a triaxial rigid body. The perturbation technique developed by Bennet (Icarus 4:177, 1965b) has been used for determination of characteristic exponents. This technique is based on Floquet’s Theory for determination of characteristic exponents in the system with periodic coefficients. The results of the study are analytical and numerical expressions are simulated for the transition curves bounding the region of stability in the μ–e plane, accurate to O(e ²). The unstable region is found to be divided into three parts. The effect of radiation parameter is significant. For small values of e, the results are in favor with the numerical analysis of Danby (Astron. J. 69:166, 1964), Bennet (Icarus 4:177, 1965b), Alfriend and Rand (AIAA J. 6:1024, 1969). The effect of radiation pressure is significant than the oblateness and triaxiality of the primaries.     

Effects of Random Flows on the Solar f Mode: II. Horizontal and Vertical Flow

We study the influence of horizontal and vertical random flows on the solar f mode in a plane-parallel, incompressible model that includes a static atmosphere. The incompressible limit is an adequate approximation for f-mode type of surface waves that are highly incompressible. The paper revisits and extends the problem investigated earlier by Murawski and Roberts (Astron. Astrophys. 272, 601, 1993). We show that the consideration of the proposed velocity profile requires several restrictive assumptions to be made. These constraints were not recognised in previous studies. The impact of the inconsistencies in earlier modelling is analysed in detail. Corrections to the dispersion relation are derived and the relevance of these corrections is analysed. Finally, the importance of the obtained results is investigated in the context of recent helioseismological data. Detailed comparison with our complementary studies on random horizontal flows (Mole, Kerekes, and Erdélyi, Solar Phys., accepted, 2008) and the random magnetic model of Erdélyi, Kerekes, and Mole (Astron. Astrophys. 431, 1083, 2005) is also given. In particular, for realistic solar parameters we find significant frequency reduction and wave damping, both of which increase with the characteristic thickness of the random layer.     

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}$ .     

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).     

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.     

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.     

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.     

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.     

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory/Helioseismic and Magnetic Imager Time?–?Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a?time?–?distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time?–?distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a?minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference travel-time perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.     

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.     

A family of exact solutions of Einstein-Maxwell field equations in isotropic coordinates: an application to optimization of quark star mass

In the present paper, we have obtained a class of charged super dense star models, starting with a static spherically symmetric metric in isotropic coordinates for perfect fluid by considering Hajj-Boutros (in J. Math. Phys. 27:1363, 1986) type metric potential and a specific choice of electrical intensity which involves a parameter K. The resulting solutions represent charged fluid spheres joining smoothly with the Reissner-Nordstrom metric at the pressure free interface. The solutions so obtained are utilized to construct the models for super-dense star like neutron stars (ρ _b =2 and 2.7×10¹⁴ g/cm³) and Quark stars (ρ _b =4.6888×10¹⁴ g/cm³). Our solution is well behaved for all values of n satisfying the inequalities \(4 < n \le4(4 + \sqrt{2} )\) and K satisfying the inequalities 0≤K≤0.24988, depending upon the value of n. Corresponding to n=4.001 and K=0.24988, we observe that the maximum mass of quark star M=2.335M _⊙ and radius R=10.04 km. Further, this maximum mass limit of quark star is in the order of maximum mass of stable Strange Quark Star established by Dong et al. (in arXiv:1207.0429v3, 2013). The robustness of our results is that the models are alike with the recent discoveries.