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.     

Retrograde resonance in the planar three-body problem

We continue the investigation of the dynamics of retrograde resonances initiated in Morais and Giuppone (Mon Notices R Astron Soc 424:52–64, doi:10.1111/j.1365-2966.2012.21151.x, 2012). After deriving a procedure to deduce the retrograde resonance terms from the standard expansion of the three-dimensional disturbing function, we concentrate on the planar problem and construct surfaces of section that explore phase-space in the vicinity of the main retrograde resonances (2/ $-$ 1, 1/ $-$ 1 and 1/ $-$ 2). In the case of the 1/ $-$ 1 resonance for which the standard expansion is not adequate to describe the dynamics, we develop a semi-analytic model based on numerical averaging of the unexpanded disturbing function, and show that the predicted libration modes are in agreement with the behavior seen in the surfaces of section.     

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.     

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.     

Accreting He-rich material onto carbon-oxygen white dwarfs until explosive carbon ignition

Type Ia supernovae(SNe Ia) play an important role in studies of cosmology and galactic chemical evolution.They are believed to be thermonuclear explosions of carbon-oxygen white dwarfs(CO WDs)when their masses approach the Chandrasekar(Ch) mass limit.However,it is still not completely understood how a CO WD increases its mass to the Ch-mass limit in the classical single-degenerate(SD) model.In this paper,we studied the mass accretion process in the SD model to examine whether the WD can explode as an SN Ia.Employing the stellar evolution code called modules for experiments in stellar astrophysics(MESA),we simulated the He accretion process onto CO WDs.We found that the WD can increase its mass to the Ch-mass limit through the SD model and explosive carbon ignition finally occurs in its center,which will lead to an SN Ia explosion.Our results imply that SNe Ia can be produced from the SD model through steady helium accretion.Moreover,this work can provide initial input parameters for explosion models of SNe Ia.     

Stability of ion-acoustic waves in a pair-ion plasma with a third species of ions: application to cometary plasmas

A popular model of a cometary plasma is hydrogen (H⁺) with positively charged oxygen (O⁺) as a heavier ion component. However, the discovery of negatively charged oxygen (O^?) ions enables one to model a cometary plasma as a pair-ion plasma (of O⁺ and O^?) with hydrogen as a third ion constituent. We have, therefore, studied the stability of the ion-acoustic wave in such a pair-ion plasma with hydrogen and electrons streaming with velocities $V_{d\mathrm{H}^{+}}$ and V _de , respectively, relative to the oxygen ions. We find the calculated frequency of the ion-acoustic wave with this model to be in good agreement with the observed frequencies. The ion-acoustic wave can also be driven unstable by the streaming velocity of the hydrogen ions. The growth rate increases with increasing hydrogen density $n_{\mathrm{H}^{+}}$ , and streaming velocities $V_{d\mathrm{H}^{+}}$ and V _de . It, however, decreases with increasing oxygen ion densities $n_{\mathrm{O}^{+}}$ and $n_{\mathrm{O}^{-}}$ .     

Secular period decreasing of detached chromospherically active binaries

The long-term orbital period changes of a large sample of detached chromospherically active binaries (CABs) were studied. Eleven such systems were found to be undergoing secular period decreases with the rates of ?6.3×10^?9 to ?1.1×10^?6 days per year. The period decreasing rates are found to vary depending on the orbital period. The longer the orbital period is, the more rapidly the period decreases. Following Stepien (Mon. Not. R. Astron. Soc. 274:1019, 1995), the period decreasing rate predicted by angular momentum loss (AML) caused by magnetic wind is computed for each system. A comparison between the observed and calculated period decreasing rates shows that the former values are obviously larger than the latter by 1–3 orders of magnitude. It suggests that the magnetic wind is not likely the determinant mechanism driving the AML in these systems. Finally, the orbital angular momentum (AM) and the rate of AML, $\dot{J}$ , are computed for each system. It shows that the AM have a similar change with the orbital period like dP/dt does, but $\log\dot{J}/J$ presents no strict changing with the kinematical ages.     

WD + He star systems as the progenitors of Type Ia supernovae and their surviving companion stars

Employing Eggleton’s stellar evolution code with the optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon–oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is ～0.3×10^?3 yr^?1, and that this channel can produce SNe Ia with short delay times (～45–140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity (>400 km/s) after the SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.     

Observational characteristics and possible asphericity of overluminous Type Ia supernovae

A few Type Ia supernovae (SNe Ia) have been suggested to be an explosion of a super-Chandrasekhar-mass white dwarf (WD) in order to account for their large luminosities, requiring a large amount of ⁵⁶Ni. However, the candidate overluminous SNe Ia 2003fg, 2006gz and (moderately overluminous) SN 1991T have very different observational features: the characteristic time-scale and velocity are very different. We examine if and how the diversity can be explained, by one-dimensional spherical radiation transport calculations covering a wide range of model parameters (e.g. WD mass). The observations of SN 2006gz are naturally explained by the super-Chandrasekhar-mass model. SN 1991T represents a marginal case, which may either be a Chandrasekhar or a super-Chandrasekhar-mass WD explosion. In contrast, the low velocity and short time-scale seen in SN 2003fg indicate that the ejecta mass is smaller than the Chandrasekhar mass, which is in apparent contradiction to the large luminosity. We suggest that the problem is solved if the progenitor WD, and thus the SN explosion, is aspherical. This may reflect a rapid rotation of the progenitor star, likely a consequence of the super-Chandrasekhar-mass WD progenitor. The observed differences between SNe 2003fg and 2006gz may be attributed to different viewing orientations.     

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.     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

Observational constraints on G-corrected holographic dark energy using a Markov chain Monte Carlo method

We constrain holographic dark energy (HDE) with time varying gravitational coupling constant in the framework of the modified Friedmann equations using cosmological data from type Ia supernovae, baryon acoustic oscillations, cosmic microwave background radiation and X-ray gas mass fraction. Applying a Markov Chain Monte Carlo (MCMC) simulation, we obtain the best fit values of the model and cosmological parameters within 1σ confidence level (CL) in a flat universe as: $\varOmega_{b}h^{2}=0.0222^{+0.0018}_{-0.0013}$ , $\varOmega_{c}h^{2}=0.1121^{+0.0110}_{-0.0079}$ , $\alpha_{G}\equiv \dot{G}/(HG) =0.1647^{+0.3547}_{-0.2971}$ and the HDE constant $c=0.9322^{+0.4569}_{-0.5447}$ . Using the best fit values, the equation of state of the dark component at the present time w _d0 at 1σ CL can cross the phantom boundary w=?1.     

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.     

Constraining the Angular Momentum of the Sun with Planetary Orbital Motions and General Relativity

The angular momentum of a star is an important astrophysical quantity related to its internal structure, formation, and evolution. Helioseismology yields $S_{\odot}= 1.92\times10^{41}\ \mathrm{kg\ m^{2}\ s^{-1}}$ for the angular momentum of the Sun. We show how it should be possible to constrain it in a near future by using the gravitomagnetic Lense?CThirring effect predicted by General Relativity for the orbit of a test particle moving around a central rotating body. We also discuss the present-day situation in view of the latest determinations of the supplementary perihelion precession of Mercury. A fit by Fienga et al. (Celestial Mech. Dynamical Astron. 111, 363, 2011) of the dynamical models of several standard forces acting on the planets of the solar system to a long data record yielded milliarcseconds per century. The modeled forces did not include the Lense?CThirring effect itself, which is expected to be as large as from helioseismology-based values of S _??. By assuming the validity of General Relativity, from its theoretical prediction for the gravitomagnetic perihelion precession of Mercury, one can straightforwardly infer $S_{\odot}\leq0.95\times10^{41}\ \mathrm{kg\, m^{2}\, s^{-1}}$ . It disagrees with the currently available values from helioseismology. Possible sources for the present discrepancy are examined. Given the current level of accuracy in the Mercury ephemerides, the gravitomagnetic force of the Sun should be included in their force models. MESSENGER, in orbit around Mercury since March 2011, will collect science data until 2013, while BepiColombo, to be launched in 2015, should reach Mercury in 2022 for a year-long science phase: the analysis of their data will be important in effectively constraining S _?? in about a decade or, perhaps, even less.     

Structure coefficients for use in stellar analysis

We present new values of the structural coefficients η _j , and related quantities, for realistic models of distorted stars in close binary systems. Our procedure involves numerical integration of Radau’s equation for detailed structural data and we verified our technique by referring to the 8-digit results of Brooker & Olle (Mon. Not. R. Astron. Soc. 115:101, 1955) for purely mathematical models. We provide tables of representative values of η _j , and related quantities, for j=2,3,…,7 for a selection of Zero Age Stellar Main Sequence (ZAMS) stellar models taken from the EZWeb compilation of the Dept. of Astronomy, University of Wisconsin-Madison. We include also some preliminary comparisons of our findings with the results of Claret and Gimenez (Astron. Astrophys. 519:A57 2010) for some observed stars.     

Perspectives on effectively constraining the location of a massive trans-Plutonian object with the New Horizons spacecraft: a sensitivity analysis

The radio tracking apparatus of the New Horizons spacecraft, currently traveling to the Pluto system where its arrival is scheduled for July 2015, should be able to reach an accuracy of 10 m (range) and 0.1  $\text{ mm } \text{ s }^{-1}$ mm s ? 1 (range-rate) over distances up to 50 au. This should allow to effectively constrain the location of a putative trans-Plutonian massive object, dubbed Planet X (PX) hereafter, whose existence has recently been postulated for a variety of reasons connected with, e.g., the architecture of the Kuiper belt and the cometary flux from the Oort cloud. Traditional scenarios involve a rock-ice planetoid with $m_\mathrm{X}\approx 0.7\,m_{\oplus }$ m X ≈ 0.7 m ⊕ at some 100–200 au, or a Jovian body with $m_\mathrm{X}\lesssim 5\,m_\mathrm{J}$ m X ? 5 m J at about 10,000–20,000 au; as a result of our preliminary sensitivity analysis, they should be detectable by New Horizons since they would impact its range at a km level or so over a time span 6 years long. Conversely, range residuals statistically compatible with zero having an amplitude of 10 m would imply that PX, if it exists, could not be located at less than about 4,500 au ( $m_\mathrm{X}=0.7\,m_{\oplus }$ m X = 0.7 m ⊕ ) or 60,000 au ( $m_\mathrm{X}=5\,m_\mathrm{J}$ m X = 5 m J ), thus making a direct detection quite demanding with the present-day technologies. As a consequence, it would be appropriate to rename such a remote body as Thelisto. Also fundamental physics would benefit from this analysis since certain subtle effects predicted by MOND for the deep Newtonian regions of our Solar System are just equivalent to those of a distant pointlike mass.     

Progenitors of type Ia supernovae

Type Ia supernovae (SNe Ia) play an important role in astrophysics and are crucial for the studies of stellar evolution, galaxy evolution and cosmology. They are generally thought to be thermonuclear explosions of accreting carbon–oxygen white dwarfs (CO WDs) in close binaries, however, the nature of the mass donor star is still unclear. In this article, we review various progenitor models proposed in the past years and summarize many observational results that can be used to put constraints on the nature of their progenitors. We also discuss the origin of SN Ia diversity and the impacts of SN Ia progenitors on some fields. The currently favourable progenitor model is the single-degenerate (SD) model, in which the WD accretes material from a non-degenerate companion star. This model may explain the similarities of most SNe Ia. It has long been argued that the double-degenerate (DD) model, which involves the merger of two CO WDs, may lead to an accretion-induced collapse rather than a thermonuclear explosion. However, recent observations of a few SNe Ia seem to support the DD model, and this model can produce normal SN Ia explosion under certain conditions. Additionally, the sub-luminous SNe Ia may be explained by the sub-Chandrasekhar mass model. At present, it seems likely that more than one progenitor model, including some variants of the SD and DD models, may be required to explain the observed diversity of SNe Ia.     

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