Resonance scattering by N2

A study is made of the intensity distribution among the bands of the Meinel and first negative system of N₂⁺ due to resonance scattering of sunlight. Absolute transition probabilities are used to calculate the relative populations among the ion states under resonance scattering conditions; the mean lifetime for deactivation is the parameter which determines the amount of resonance scattering. Photon scattering rates are calculated for most of the ion bands and it is suggested that an appropriate value for the 3914 Å band would be 0·050 photons/ sec per ion. Observations of the Δυ = −1 sequence of the first negative system in the twilight spectrum are reported. Extended vibrational development is detected which indicates that only about 80 per cent of the emission is resonance scattered. Sunlit auroral spectra of N₂⁺, however, which have been generally considered to be due predominantly to resonance scattering, indicates only about 40 per cent of the emission is due to resonance scattering. Measurable effects resulting from a charge-transfer ion source (O⁺(²D)) are predicted.     

Druse clinopyroxene in D'Orbigny angritic meteorite studied by single‐crystal X‐ray diffraction,electron microprobe analysis,and Mössbauer spectroscopy

Abstract— The crystal structure of druse clinopyroxene from the D'Orbigny angrite, (Ca_0.944 Fe²⁺_0.042 Mg_0.010Mn_0.004) (Mg_0.469Fe²⁺_0.317Fe³⁺_0.035Al_0.125Cr_0.010Ti_0.044) (Si_1.742Al_0.258) O₆, a = 9.7684(2), b = 8.9124(2), c = 5.2859(1) Å, β = 105.903(1)°, V = 442.58 ų, space group C2/c, Z = 2, has been refined to an R₁ index of 1.92% using single‐crystal X‐ray diffraction data. The unit formula, calculated from electron microprobe analysis, and the refined site scattering values were used to assign site populations. The distribution of Fe²⁺and Mg over the M1 and M2 sites suggests a closure temperature of 1000 °C. Mössbauer spectroscopy measurements were done at room temperature on a single crystal and a powdered sample. The spectra are adequately fit by a Voigt‐based quadrupole‐splitting distribution model having two generalized sites, one for Fe²⁺with two Gaussian components and one for Fe³⁺with one Gaussian component. The two ferrous components are assigned to Fe²⁺at the M1 site, and arise from two different next‐nearest‐neighbor configurations of Ca and Fe cations at the M2 site: (3Ca,0Fe) and (2Ca,1Fe). The Fe³⁺/Fe_tot ratio determined by Mössbauer spectroscopy is in agreement with that calculated from the electron microprobe analysis. The results are discussed in connection with the redox and thermal history of D'Orbigny.     

The dust scattering model cannot explain the shallow X-ray decay in GRB afterglows

A dust scattering model was recently proposed to explain the shallow X-ray decay (plateau) observed prevalently in Gamma-Ray Burst (GRB) early afterglows. In this model, the plateau is the scattered prompt X-ray emission by the dust located close (about 10 to a few hundred pc) to the GRB site. In this paper, we carefully investigate the model and find that the scattered emission undergoes strong spectral softening with time, due to the model's essential ingredient that harder X-ray photons have smaller scattering angle thus arrive earlier, while softer photons suffer larger angle scattering and arrive later. The model predicts a significant change, that is  Δβ∼ 2–3  , in the X-ray spectral index from the beginning of the plateau towards the end of the plateau, while the observed data show close to zero softening during the plateau and the plateau-to-normal transition phase. The scattering model predicts a big difference between the harder X-ray light curve and the softer X-ray light curve, i.e. the plateau in harder X-rays ends much earlier than in softer X-rays. This feature is not seen in the data. The large scattering optical depths of the dust required by the model imply strong extinction in optical,   A _V ≳ 10  , which contradicts current findings of   A _V = 0.1–0.7  from optical and X-ray afterglow observations. We conclude that the dust scattering model cannot explain the X-ray plateaus.     

Rms-value and power spectrum of the photospheric intensity-fluctuations

The power spectrum and the rms-value of the granular intensity fluctuations were studied using granulation photographs of excellent quality obtained during the JOSO site testing campaign 1979 at Izaña. The observed power spectrum was corrected using various effective modulation transfer functions of the system: telescope+aberrations+atmospheric seeing, assuming different contributions of the atmospheric seeing. With this procedure a lower and upper limit for the ‘true’ power spectrum of the granular intensity fluctuations and thus for the rms-value could be derived: 7.2% <I _rms <12% at λ = 550 nm, with a most probable value of I _rms = 10.5%. We checked the validity of the upper limit by applying to our data a MTF (Deubner and Mattig, 1975), which certainly must lead to an overcorrection. This procedure lead to I _rms = 13.4%. Thus we can state that the true rms-value of the granular intensity fluctuations does certainly not exceed 13% at λ = 550 nm.     

A Monte Carlo study of polarization structures in the Thomson-scattered line radiation

Thomson scattering is often invoked to explain broad wing features that are seen in various objects including active galactic nuclei and symbiotic stars. Despite the wavelength-independent scattering cross-section of Thomson scattering, the line flux may exhibit wavelength-dependent linear degree of polarization, because various parts of emission wings are contributed by photons with different scattering numbers. Specifically, more scattered and hence more weakly polarized photons tend to fill the farther wing parts from the line centre, while the neighbourhood of the line centre is dominated by less-scattered photons with higher degree of polarization. Using a Monte Carlo technique, we investigate the polarization structure of Thomson-scattered line radiation. A detailed analysis of polarization structure formation is conducted by investigating the dependence of the polarization and profile width on the scattering number for various finite electron scattering slabs. Significantly varying degree of polarization is obtained when the scattering medium has Thomson optical depth  τ_Th≥ 1  . We present our high-resolution spectrum of the symbiotic star V1016 Cyg obtained with the Bohyunsan Optical Echelle Spectrograph (BOES) in order to fit the broad profile around Hα by electron scattering wings adopting an oblate spheroidal geometry with Thomson optical depth  τ_Th= 0.5  and electron temperature   T _e= 6.2 × 10⁴  K  . Local maxima in the linear degree of polarization of Thomson-scattered line radiation are expected to appear in the spectral regions characterized by the average scattering number ≃1.     

Exact solution of the equation of transfer in a finite exponential atmosphere by the method of laplace transform and linear singular operators

An approximate solution of the transfer equation for coherent scattering in stellar atmospheres with Planck's function as a nonlinear function of optical depth, viz.,

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.     

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.     

Color gradients of spiral disks in the SIoan Digital Sky Survey

We investigate the radial color gradients of galactic disks using a sample of - 20 000 face-on spiral galaxies selected from the fourth data release of the Sloan Digital Sky Survey （SDSS-DR4）. We combine galaxies with similar concentrations, sizes and luminosities to construct composite galaxies, and then measure their color profiles by stacking the azimuthally averaged radial color profiles of all the member galaxies. Except for the smallest galaxies （R50 〈 3 kpc）, almost all galaxies show negative disk color gradients with mean 9 - r gradient Ggr = -0.006 magkpc-1 and r - z gradient Grz = -0.018 mag kpc^-1. The disk color gradients are independent of the morphological types of galaxies and strongly dependent on the disk surface brightness μd, with lower surface brightness galactic disks having steeper color gradients. We quantify the intrinsic correlation between color gradients and surface brightness as Ggr = -0.011μd ＋ 0.233 and Grz - -0.015μd ＋ 0.324. These quantified correlations provide tight observational constraints on the formation and evolution models of spiral galaxies.     

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.     

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.     

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.     

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.     

Solution of the equation of transfer for coherent scattering in an exponential atmosphere by the method of discrete ordinates

A solution of the transfer equation for coherent scattering in stellar atmosphere with Planck's function as a nonlinear function of optical depth, viz., $$B_v (T) = b_0 + b_1 {\text{ }}e^{ - \beta \tau } $$ is obtained by the method of discrete ordinates originally due to Chandrasekhar.     

Voyager disk resolved photometry of the Saturnian satellites

Photometric analysis of Voyager images of the medium-sized icy satellites of Saturn shows that their surfaces exhibit a wide range of scattering properties. At low phase angles, Rhea and Dione closely follow lunar behavior with almost no limb darkening. Mimas, Tethys, and especially Enceladus shiw significant limb darkening at low phase angles, which suggests multiple scattering is important for their surfaces. A simple photometric function of the form I/F = f(α)Aμ₀/(μ + μ₀) + (1 ? A)μ₀ has been fit to the observations. For normal reflectances <0.6, we find lunar-like scattering properties (A = 1). No satellite's surface can be described by Lambert's Law (A = 0). Dione exhibits the widest albedo variations (about 50%). A longitudinal dark stripe which represents a 15% decrease in albedo is situated near the center of the trailing side of Tethys. A correlation is found between the albedo and color of the satellites: the darker objects are redder. Similarly, darker areas of each satellite are redder. Spectral reflectances of Mimas and Enceladus can be derived for the first time. After the proper calibrations to the Voyager color images are made, it is found that both satellites have remarkably flat spectra into the ultraviolet.     

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.     

LO Pegasi: an investigation of multiband optical polarization

We present BVR polarimetric study of the cool active star LO Pegasi (LO Peg) for the first time. LO Peg was found to be highly polarized among the cool active stars. Our observations yield average values of polarization in LO Peg:   P_B = 0.387 ± 0.004 per cent, θ_B= 88°± 1°; P_V = 0.351 ± 0.004 per cent, θ_V= 91°± 1°  and   P_R = 0.335 ± 0.003 per cent, θ_R= 91°± 1°  . Both the degree of polarization and the position angle are found to be variable. The semi-amplitude of the polarization variability in B, V and R bands is found to be  0.18 ± 0.02, 0.13 ± 0.01  and  0.10 ± 0.02  per cent, respectively. We suggest that the levels of polarization observed in LO Peg could be the result of scattering of an anisotropic stellar radiation field by an optically thin circumstellar envelope or scattering of the stellar radiation by prominence-like structures.     

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.     

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.