Archive for category Solar and Stellar

Today's Postings

Magnetic Fields in Astrophysical Jets: From Launch to Termination

Long-lived, stable jets are observed in a wide variety of systems, from protostars, through Galactic compact objects to active galactic nuclei (AGN). Magnetic fields play a central role in launching, accelerating, and collimating the jets through various media. The termination of jets in molecular clouds or the interstellar medium deposits enormous amounts of mechanical energy and momentum, and their interactions with the external medium, as well, in many cases, as the radiation processes by which they are observed, are intimately connected with the magnetic fields they carry. This review focuses on the properties and structures of magnetic fields in long-lived jets, from their launch from rotating magnetized young stars, black holes, and their accretion discs, to termination and beyond. We compare the results of theory, numerical simulations, and observations of these diverse systems and address similarities and differences between relativistic and non-relativistic jets in protostellar versus AGN systems. On the observational side, we focus primarily on jets driven by AGN because of the strong observational constraints on their magnetic field properties, and we discuss the links between the physics of these jets on all scales.

Short Gas Dissipation Timescales: Diskless Stars in Taurus and Chamaeleon I

We present an Advanced Camera for Surveys/ Solar Blind Channel far-ultraviolet (FUV) study of \h2 gas in 12 weak T Tauri stars in nearby star-forming regions. The sample consists of sources which have no evidence of inner disk dust. Our new FUV spectra show that in addition to the dust, the gas is depleted from the inner disk. This sample is combined with a larger FUV sample of accretors and non-accretors with ages between 1 and 100 Myr, showing that as early as 1–3 Myr, systems both with and without gas are found. Possible mechanisms for depleting gas quickly include viscous evolution, planet formation and photoevaporation by stellar radiation fields. Since these mechanisms alone cannot account for the lack of gas at 1–3 Myr, it is likely that the initial conditions (e.g. initial disk mass or core angular momentum) contribute to the variety of disks observed at any age. We estimate the angular momentum of a cloud needed for most of the mass to fall very close to the central object and compare this to models of the expected distribution of angular momenta. Up to 20% of cloud cores have low enough angular momenta to form disks with the mass close to the star, which would then accrete quickly; this percentage is similar to the fraction of diskless stars in the youngest star forming regions. With our sample, we characterize the chromospheric contribution to the FUV luminosity and find that $L_{FUV}/L_{bol}$ saturates at $\sim10^{-4.1}$.

Really focused stellar winds in X-ray binaries

We investigate the anisotropy of stellar winds in binaries to improve the models of accretion in high-mass X-ray binaries. We model numerically the stellar wind from a supergiant component of a binary in radial and three-dimensional radiation hydrodynamic approximation taking into account the Roche potential, Coriolis force, and radiative pressure in the continuum and spectral lines. The Coriolis force influences substantially the mass loss and thus also the accretion rate. The focusing of the stellar wind by the gravitational field of the compact companion leads to the formation of a gaseous tail behind the companion.

The Great Eruption of Eta Carinae

During the years 1838-1858, the very massive star {\eta} Carinae became the prototype supernova impostor: it released nearly as much light as a supernova explosion and shed an impressive amount of mass, but survived as a star.1 Based on a light-echo spectrum of that event, Rest et al.2 conclude that “a new physical mechanism” is required to explain it, because the gas outflow appears cooler than theoretical expectations. Here we note that (1) theory predicted a substantially lower temperature than they quoted, and (2) their inferred observational value is quite uncertain. Therefore, analyses so far do not reveal any significant contradiction between the observed spectrum and most previous discussions of the Great Eruption and its physics.

Improved Variable Star Search in Large Photometric Data Sets -- New Variables in CoRoT Field LRa02 Detected by BEST II [Replacement]

The CoRoT field LRa02 has been observed with the Berlin Exoplanet Search Telescope II (BEST II) during the southern summer 2007/2008. A first analysis of stellar variability led to the publication of 345 newly discovered variable stars. Now, a deeper analysis of this data set was used to optimize the variability search procedure. Several methods and parameters have been tested in order to improve the selection process compared to the widely used J index for variability ranking. This paper describes an empirical approach to treat systematic trends in photometric data based upon the analysis of variance statistics that can significantly decrease the rate of false detections. Finally, the process of reanalysis and method improvement has virtually doubled the number of variable stars compared to the first analysis by Kabath et al. A supplementary catalog of 272 previously unknown periodic variables plus 52 stars with suspected variability is presented. Improved ephemerides are given for 19 known variables in the field. In addition, the BEST II results are compared with CoRoT data and its automatic variability classification.

Improved Variable Star Search in Large Photometric Data Sets -- New Variables in CoRoT Field LRa02 Detected ba BEST II

The CoRoT field LRa02 has been observed with the Berlin Exoplanet Search Telescope II (BEST II) during the southern summer 2007/2008. A first analysis of stellar variability led to the publication of 345 newly discovered variable stars. Now, a deeper analysis of this data set was used to optimize the variability search procedure. Several methods and parameters have been tested in order to improve the selection process compared to the widely used J index for variability ranking. This paper describes an empirical approach to treat systematic trends in photometric data based upon the analysis of variance statistics that can significantly decrease the rate of false detections. Finally, the process of reanalysis and method improvement has virtually doubled the number of variable stars compared to the first analysis by Kabath et al. A supplementary catalog of 272 previously unknown periodic variables plus 52 stars with suspected variability is presented. Improved ephemerides are given for 19 known variables in the field. In addition, the BEST II results are compared with CoRoT data and its automatic variability classification.

Stream Overflow in Z Cha and OY Car during Quiescence

Eclipses of the hot spot in Z Cha and OY Car observed by many authors during their quiescence are re-analyzed. Distances of the spot from the center of the disk r_s are determined from phases of ingress and egress. In the case of several eclipses of Z Cha and nearly all eclipses of OY Car it is found that r_s(egress) < r_s(ingress). This implies that they are not representative for the radius of the disk r_d and is interpreted as being due to the stream overflow. The r_s(ingress) – dt relations (where dt is the time since last outburst) are improved when points affected by the stream overflow are omitted.

Scaling of the growth rate of magnetic islands in the heliosheath

Current sheets thinner than the ion inertial length are unstable to the tearing instability and will develop magnetic islands that grow due to magnetic reconnection. We investigate whether the growth of magnetic islands in a current sheet can continue indefinitely, or in the case of the heliosheath until reaching a neighboring current sheet, and at what rate the islands grow. We investigate the development and growth of magnetic islands using a particle-in-cell code, starting from particle noise. Performing a scaling of the growth of magnetic islands versus the system size, we find that the growth rate is independent of the system size up to the largest simulation we were able to complete. The islands are able to continue growing as long as they merge with each other and maintain a high aspect ratio. Otherwise there is not enough magnetic tension to sustain reconnection. When applied to the sectored magnetic fields in the heliosheath, we show that the islands can continue growing until they reach the sector width and do so in much less time than it takes for the islands to convect through the heliosheath.

Constraining Tidal Dissipation in Stars from The Destruction Rates of Exoplanets

We use the distribution of extrasolar planets in circular orbits around stars with surface convective zones detected by ground based transit searches to constrain how efficiently tides raised by the planet are dissipated on the parent star. We parameterize this efficiency as a tidal quality factor (Q*). We conclude that the population of currently known planets is inconsistent with Q*<10^7 at the 99% level. Previous studies show that values of Q* between 10^5 and 10^7 are required in order to explain the orbital circularization of main sequence low mass binary stars in clusters, suggesting that different dissipation mechanisms might be acting in the two cases, most likely due to the very different tidal forcing frequencies relative to the stellar rotation frequency occurring for star–star versus planet–star systems.

Time series of high resolution photospheric spectra in a quiet region of the Sun. II. Analysis of the variation of physical quantities of granular structures

From the inversion of a time series of high resolution slit spectrograms obtained from the quiet sun, the spatial and temporal distribution of the thermodynamical quantities and the vertical flow velocity is derived as a function of logarithmic optical depth and geometrical height. Spatial coherence and phase shift analyzes between temperature and vertical velocity depict the height variation of these physical quantities for structures of different size. An average granular cell model is presented, showing the granule-intergranular lane stratification of temperature, vertical velocity, gas pressure and density as a function of logarithmic optical depth and geometrical height. Studies of a specific small and a specific large granular cell complement these results. A strong decay of the temperature fluctuations with increasing height together with a less efficient penetration of smaller cells is revealed. The T -T coherence at all granular scales is broken already at log tau =-1 or z~170 km. At the layers beyond, an inversion of the temperature contrast is revealed. Vertical velocities are in phase throughout the photosphere and penetrate into the highest layers under study.

Time series of high resolution photospheric spectra in a quiet region of the Sun. I. Analysis of global and spatial variations of line parameters

A 50 min time series of one-dimensional slit-spectrograms, taken in quiet sun at disk center, observed at the German Vacuum Tower Telescope (Observatorio del Teide), was used to study the global and spatial variations of different line parameters. In order to determine the vertical structure of the photosphere two lines with well separated formation heights have been considered. The data have been filtered of p-modes to isolate the pure convective phenomenon. From our studies of global correlation coefficients and coherence and phase shift analyzes between the several line parameters, the following results can be reported. The convective velocity pattern preserves structures larger than 1.0″ up to the highest layers of the photosphere (~ 435 km). However, at these layers, in the intensity pattern only structures larger than 2.0″ are still connected with those at the continuum level although showing inverted brightness contrast. This confirms an inversion of temperature that we have found at a height of ~140 km. A possible evidence of gravity waves superimposed to the convective motions is derived from the phase shift analysis. We interpret the behavior of the full width at half maximum and the equivalent width as a function of the distance to the granular borders, as a consequence of enhanced turbulence and/or strong velocity gradients in the intergranular lanes.

A new soft X-ray spectral model for polars with an application to AM Herculis

We present a simple heuristic model for the time-averaged soft X-ray temperature distribution in the accretion spot on the white dwarf in polars. The model is based on the analysis of the Chandra LETG spectrum of the prototype polar AM Her and involves an exponential distribution of the emitting area vs. blackbody temperature a(T) = a0 exp(-T/T0). With one free parameter besides the normalization, it is mathematically as simple as the single blackbody, but is physically more plausible and fits the soft X-ray and far-ultraviolet spectral fluxes much better. The model yields more reliable values of the wavelength-integrated flux of the soft X-ray component and the implied accretion rate than reported previously.

Period decrease in three SuperWASP eclipsing binary candidates near the short-period limit

SuperWASP light curves for 53 W UMa-type eclipsing binary (EB) candidates, identified in previous work as being close to the contact binary short-period limit, were studied for evidence of period change. The orbital periods of most of the stars were confirmed, and period decrease, significant at more than 5 sigma, was observed in three objects: 1SWASP J174310.98+432709.6 (-0.055 \pm0.003 s/yr), 1SWASP J133105.91+121538.0 (-0.075 \pm0.013 s/yr) and 1SWASP J234401.81-212229.1 (-0.313 \pm0.019 s/yr). The magnitudes of the observed period changes cannot be explained by magnetic braking or gravitational radiation effects, and are most likely primarily due to unstable mass transfer from primary to secondary components, possibly accompanied by unstable mass and angular momentum loss from the systems. If these period decreases persist, the systems could merge on a relatively short timescale.

Infant mortality without gas expulsion? The rapid dispersal of virialised low mass clusters

Infant mortality brought about by the expulsion of a star cluster’s natal gas is widely invoked to explain cluster statistics at different ages. While a well studied problem, most recent studies of gas expulsion’s effect on a cluster have focused on massive clusters, with stellar counts of order $10^4$. Here we argue that the evolutionary timescales associated with the compact low-mass clusters typical of the median cluster in the Solar neighborhood are short enough that significant dynamical evolution can take place over the ages usually associated with gas expulsion. To test this we perform {\it N}-body simulations of the dynamics of a very young star forming region, with initial conditions drawn from a large-scale hydrodynamic simulation of gravitational collapse and fragmentation. The subclusters we analyse have high local star formation efficiencies and are roughly virialised, and have populations of a few hundred stars. Over 10 Myr they expand to a similar degree as would be expected from gas expulsion if they were initially gas-rich, but the expansion is purely due to the internal stellar dynamics of the young clusters. The expansion is such that the stellar densities at 2 Myr match those of YSOs in the Solar neighborhood. We argue that at the low-mass end of the cluster mass spectrum, a deficit of clusters at 10s of Myr does not necessarily imply gas expulsion as a disruption mechanism.

Long-term Variations of Solar Differential Rotation and Sunspot Activity: Revisited

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential rotation law obtained by Makarov, Tlatov, and Callebaut (1997), Javaraiah, Bertello, and Ulrich (2005a, b), and Javaraiah et al. (2009). Our results indicate that the solar surface rotation rate at the Equator (indicated by the A parameter of the standard solar rotation law) shows a secular decrease since cycle 12 onwards, given by about $1\,-\,1.5\times10^{-3}$($deg\ day^{-1} year^{-1}$). The B parameter of the standard differential rotation law seems to also show a secular decrease since cycle 12 onwards, but of weak statistical significance. The rotation rate averaged on latitudes ($0^{o}\,–\,40^{o}$) does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation in the long run.

Direct modeling of neutral helium in the heliosphere

Several years of neutral particle measurements by the NASA/IBEX mission have yielded direct observations of interstellar neutral helium and oxygen. The data indicate the presence of secondary neutral helium and oxygen, which are created within the heliosphere by charge exchange involving helium or oxygen ions. This contribution describes a detailed conserving calculation method based on Keplerian orbits that has been developed to characterize helium distribution functions throughout the heliosphere, in particular in the innermost heliosphere, while accounting for loss and production of neutral particles along their path. Coupled with global heliosphere models of plasma distributions, this code is useful for predicting the fluxes of heavy neutral atoms at spacecraft detectors, so enabling inferences on the characteristics of the interstellar medium.

On the Amplitude of Convective Velocities in the Deep Solar Interior

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5–30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

Circumbinary Gas Accretion onto a Central Binary: Infrared Molecular Hydrogen Emission from GG Tau A

We present high spatial resolution maps of ro-vibrational molecular hydrogen emission from the environment of the GG Tau A binary component in the GG Tau quadruple system. The H2 v= 1-0 S(1) emission is spatially resolved and encompasses the inner binary, with emission detected at locations that should be dynamically cleared on several hundred-year timescales. Extensions of H2 gas emission are seen to ~100 AU distances from the central stars. The v = 2-1 S(1) emission at 2.24 microns is also detected at ~30 AU from the central stars, with a line ratio of 0.05 +/- 0.01 with respect to the v = 1-0 S(1) emission. Assuming gas in LTE, this ratio corresponds to an emission environment at ~1700 K. We estimate that this temperature is too high for quiescent gas heated by X-ray or UV emission from the central stars. Surprisingly, we find that the brightest region of H2 emission arises from a spatial location that is exactly coincident with a recently revealed dust “streamer” which seems to be transferring material from the outer circumbinary ring around GG Tau A into the inner region. As a result, we identify a new excitation mechanism for ro-vibrational H2 stimulation in the environment of young stars. The H2 in the GG Tau A system appears to be stimulated by mass accretion infall as material in the circumbinary ring accretes onto the system to replenish the inner circumstellar disks. We postulate that H2 stimulated by accretion infall could be present in other systems, particularly binaries and “transition disk” systems which have dust cleared gaps in their circumstellar environments.

The star formation history of RCW 36

Recent studies of massive-star forming regions indicate that they can contain multiple generations of young stars. These observations suggest that star formation in these regions is sequential and/or triggered by a previous generation of (massive) stars. Here we present new observations of the star forming region RCW 36 in the Vela Molecular Ridge, hosting a young cluster of massive stars embedded in a molecular cloud complex. In the periphery of the cluster several young stellar objects (YSOs) are detected which produce bipolar jets (HH 1042 and HH 1043) demonstrating that these objects are still actively accreting. The VLT/X-shooter spectrum of the jet structure of HH 1042 provides detailed information on the physical conditions and kinematical properties of the jet plasma. From this information the YSO’s accretion history can be derived. Combining the photometric and spectroscopic observations of RCW 36 gives insight into the formation process of individual stars and the star formation history of this young massive-star forming region.

RX J0942.7-7726AB: an isolated pre-main sequence wide binary

We report the discovery of two young M-dwarfs, RX J0942.7-7726 (M1) and 2MASS J09424157-7727130 (M4.5), that were found only 42 arcsec apart in a survey for pre-main sequence stars surrounding the open cluster eta Chamaeleontis. Both stars have congruent proper motions and near-infrared photometry. Medium-resolution spectroscopy reveals that they are coeval (age 8-12 Myr), codistant (100-150 pc) and thus almost certainly form a true wide binary with a projected separation of 4000-6000 AU. The system appears too old and dynamically fragile to have originated in eta Cha and a traceback analysis argues for its birth in or near the Scorpius-Centaurus OB Association. RX J0942.7-7726AB joins a growing group of wide binaries kinematically linked to Sco-Cen, suggesting that such fragile systems can survive the turbulent environment of their natal molecular clouds while still being dispersed with large velocities. Conversely, the small radial velocity difference between the stars (2.7 \pm 1.0 km/s) could mean the system is unbound, a result of the coincidental ejection of two single stars with similar velocity vectors from the OB association early in its evolution.

DRAFTS: A Deep, Rapid Archival Flare Transient Search in the Galactic Bulge

We utilize the Sagittarius Window Eclipsing Extrasolar Planet Search (SWEEPS) HST/ACS dataset for a Deep Rapid Archival Flare Transient Search (DRAFTS) to constrain the flare rate toward the older stellar population in the Galactic bulge. During 7 days of monitoring 229,293 stars brighter than V=29.5, we find evidence for flaring activity in 105 stars between V=20 and V=28. We divided the sample into non-variable stars and variable stars whose light curves contain large-scale variability. The flare rate on variable stars is \sim 700 times that of non-variable stars, with a significant correlation between the amount of underlying stellar variability and peak flare amplitude. The flare energy loss rates are generally higher than those of nearby well-studied single dMe flare stars. The distribution of proper motions is consistent with the flaring stars being at the distance and age of the Galactic bulge. If they are single dwarfs, they span a range of \approx 1.0 – 0.25M\odot. A majority of the flaring stars exhibit periodic photometric modulations with P <3d. If these are tidally locked magnetically active binary systems, their fraction in the bulge is enhanced by a factor of \sim20 compared to the local value. These stars may be useful for placing constraints on the angular momentum evolution of cool close binary stars. Our results expand the type of stars studied for flares in the optical band, and suggest that future sensitive optical time-domain studies will have to contend with a larger sample of flaring stars than the M dwarf flare stars usually considered.

The Fundamental Parameters of Four Massive Eclipsing Binaries in Westerlund 1

Only a small number of high mass stars (> 30 Mo) have fundamental parameters measured with high enough accuracy from eclipsing binaries to constrain formation and evolutionary models of massive stars. This work aims to increase this limited sample, by studying the 4 massive eclipsing binary candidates discovered by Bonanos in the young massive cluster Westerlund 1. We present new follow-up echelle spectroscopy of these binaries and models of their light and radial velocity curves. We obtain fundamental parameters (i.e. masses, radii) for the 8 component stars, finding masses that span a range of 10-40 Mo, and contributing accurate fundamental parameters for 1 additional very massive star, the 33 Mo component of W13. WR77o is found to have a ~40 Mo companion, which provides a second dynamical constraint on the mass of the progenitor of the magnetar known in the cluster. We also use W13 to estimate the first, direct, eclipsing binary distance to Westerlund 1 and therefore the magnetar, and find it to be at 4.0 +/- 0.6 kpc. Our results confirm previous evidence for a high mass for the progenitor of the magnetar. In addition, the availability of eclipsing binaries with accurate parameters opens the way for direct, independent, high precision eclipsing binary distance measurements to Westerlund 1.

The complex case of V445 Lyr observed with Kepler: Two Blazhko modulations, a non-radial mode, possible triple mode RR Lyrae pulsation, and more

Rapid and strong changes in the Blazhko modulation of RR Lyrae stars, as they have recently been detected in high precision satellite data, have become a crucial topic in finding an explanation of the long-standing mystery of the Blazhko effect. We present here an analysis of the most extreme case detected so far, the RRab star V445 Lyr (KIC 6186029) which was observed with the Kepler space mission. V445 Lyr shows very strong cycle-to-cycle changes in its Blazhko modulation, which are caused both by a secondary long-term modulation period as well as irregular variations. In addition to the complex Blazhko modulation, V445 Lyr also shows a rich spectrum of additional peaks in the frequency range between the fundamental pulsation and the first harmonic. Among those peaks, the second radial overtone could be identified, which, combined with a metallicity estimate of [Fe/H]=-2.0 dex from spectroscopy, allowed to constrain the mass (0.55-0.65 M_sun) and luminosity (40-50 L_sun) of V445 Lyr through theoretical Petersen diagrams. A non-radial mode as well as possibly the first overtone are also excited. Furthermore, V445 Lyr shows signs of the period doubling phenomenon and a long term period change. A detailed Fourier analysis along with a study of the O-C variation of V445 Lyr is presented, and the origin of the additional peaks and possible causes of the changes in the Blazhko modulation are discussed. The results are then put into context with those of the only other star with a variable Blazhko effect for which a long enough set of high precision continuous satellite data has been published so far, the CoRoT star 105288363.

Thermohaline instabilities inside stars: a synthetic study including external turbulence and radiative levitation

We have derived a new expression for the thermohaline mixing coefficient in stars, including the effects of radiative levitation and external turbulence, by solving Boussinesq equations in a quasi-incompressible fluid with a linear approximation. It is well known that radiative levitation of individual elements can lead to their accumulation in specific stellar layers. In some cases, it can induce important effects on the stellar structure. Here we confirm that this accumulation is moderated by thermohaline convection due to the resulting inverse $\mu$-gradient. The new coefficient that we have derived shows that the effect of radiative accelerations on the thermohaline instability itself is small. This effect must however be checked in all computations. We also confirm that the presence of large horizontal turbulence can reduce or even suppress the thermohaline convection. These results are important as they concern all the cases of heavy element accumulation in stars. The computations of radiative diffusion have to be revisited including thermohaline convection and its consequences. It may be one of the basic reasons for the fact that the observed abundances are always smaller than those predicted by pure atomic diffusion. In any case, these processes have to compete with rotation-induced mixing, but this competition is more complex than previously thought due to their mutual interaction.

Prominence-cavity regions observed using SWAP 174A filtergrams and simultaneous eclipse flash spectra

Images from the SWAP (Proba 2 mission) taken at 174A in the Fe IX/X lines are compared to simultaneous slitless flash spectra taken during the last solar total eclipse of July, 11th 2010. Many faint low excitation emission lines together with the HeI and HeII Paschen Alpha chromospheric lines are recorded on eclipse spectra where regions of limb prominences are obtained with space-borne imagers. We consider a deep flash spectrum obtained by summing 80 individual spectra to show the intensity modulations of the continuum. Intensity depressions are observed around the prominences in both eclipse and SWAP images. The prominence cavities are interpreted as a relative depression of plasma density, produced inside the corona surrounding the prominences. Photometric measurements are shown at different scales and different, spectrally narrow, intervals for both the prominences and the coronal background.

The Period-Luminosity Relation of Red Supergiant Stars in the Small Magellanic Cloud

The characteristics of light variation of RSGs in SMC are analyzed based on the nearly 8-10 year long data collected by the ASAS and MACHO projects. The identified 126 RSGs are classified into five categories accordingly: 20 with poor photometry, 55 with no reliable period, 6 with semi-regular variation, 15 with Long Secondary Period (LSP) and distinguishable short period and 30 with only LSP. For the semi-regular variables and the LSP variables with distinguishable short period, the Ks band period-luminosity (P-L) relation is analyzed and compared with that of the Galaxy, LMC and M33. It is found that the RSGs in these galaxies obey similar P-L relation except the Galaxy. In addition, the P-L relations in the infrared bands, namely the 2MASS JHKs, Spitzer/IRAC and Spitzer/MIPS 24 {\mu}m bands, are derived with high reliability. The best P-L relation occurs in the Spitzer/IRAC [3.6] and [4.5] bands. Based on the comparison with the theoretical calculation of the P-L relation, the mode of pulsation of RSGs in SMC is suggested to be the first overtone radial mode.

Surface flux evolution constraints for flux transport dynamos

The surface flux transport (SFT) model of solar magnetic fields involves empirically well-constrained velocity and magnetic fields. The basic evolution of the Sun’s large-scale surface magnetic field is well described by this model. The azimuthally averaged evolution of the SFT model can be compared to the surface evolution of the flux transport dynamo (FTD), and the evolution of the SFT model can be used to constrain several near-surface properties of the FTD model. We compared the results of the FTD model with different upper boundary conditions and diffusivity profiles against the results of the SFT model. Among the ingredients of the FTD model, downward pumping of magnetic flux, related to a positive diffusivity gradient, has a significant effect in slowing down the diffusive radial transport of magnetic flux through the solar surface. Provided the pumping was strong enough to give rise to a downflow of a magnetic Reynolds number of 5 in the near-surface boundary layer, the FTD using a vertical boundary condition matches the SFT model based on the average velocities above the boundary layer. The FTD model with a potential field were unable to match the SFT results.

Free-free Emission and Radio Recombination Lines from Photoevaporating Disks

Recent infrared observations have demonstrated that photoevaporation driven by high-energy photons from the central star contributes to the dispersal of protoplanetary disks. Here, we show that photoevaporative winds should produce a detectable free-free continuum emission given the range of stellar ionizing photons and X-ray luminosities inferred for young sun-like stars. We point out that VLA observations of the nearby disk around TWHya might have already detected this emission at centimeter wavelengths and calculate the wind electron density and mass flow rate. We also estimate the intensities of H radio recombination lines tracing the wind and discuss which ones could be detected with current instrumentation. The detection and profiles of these recombination lines would unambiguously prove our inference of free-free emission from photoevaporating disks like TWHya. In addition, radio/millimeter data can help constraining wind parameters such as temperature and electron density that are fundamental in measuring mass flow rates.

Discovery of Super-Li Rich Red Giants in Dwarf Spheroidal Galaxies [Replacement]

Stars destroy lithium (Li) in their normal evolution. The convective envelopes of evolved red giants reach temperatures of millions of K, hot enough for the 7Li(p,alpha)4He reaction to burn Li efficiently. Only about 1% of first-ascent red giants more luminous than the luminosity function bump in the red giant branch exhibit A(Li) > 1.5. Nonetheless, Li-rich red giants do exist. We present 15 Li-rich red giants–14 of which are new discoveries–among a sample of 2054 red giants in Milky Way dwarf satellite galaxies. Our sample more than doubles the number of low-mass, metal-poor ([Fe/H] <~ -0.7) Li-rich red giants, and it includes the most-metal poor Li-enhanced star known ([Fe/H] = -2.82, A(Li)_NLTE = 3.15). Because most of these stars have Li abundances larger than the universe's primordial value, the Li in these stars must have been created rather than saved from destruction. These Li-rich stars appear like other stars in the same galaxies in every measurable regard other than Li abundance. We consider the possibility that Li enrichment is a universal phase of evolution that affects all stars, and it seems rare only because it is brief.

Discovery of Super-Li Rich Red Giants in Dwarf Spheroidal Galaxies

Stars destroy lithium (Li) in their normal evolution. The convective envelopes of evolved red giants reach temperatures of millions of K, hot enough for the 7Li(p,alpha)4He reaction to burn Li efficiently. Only about 1% of first-ascent red giants more luminous than the luminosity function bump in the red giant branch exhibit A(Li) > 1.5. Nonetheless, Li-rich red giants do exist. We present 15 Li-rich red giants–14 of which are new discoveries–among a sample of 2054 red giants in Milky Way dwarf satellite galaxies. Our sample more than doubles the number of low-mass, metal-poor ([Fe/H] <~ -0.7) Li-rich red giants, and it includes the most-metal poor Li-enhanced star known ([Fe/H] = -2.82, A(Li)_NLTE = 3.15). Because most of these stars have Li abundances larger than the universe's primordial value, the Li in these stars must have been created rather than saved from destruction. These Li-rich stars appear like other stars in the same galaxies in every measurable regard other than Li abundance. We consider the possibility that Li enrichment is a universal phase of evolution that affects all stars, and it seems rare only because it is brief.

New Variable Stars in the Globular Cluster NGC 6584

Using the Southeastern Association for Research in Astronomy 0.6 meter telescope located at Cerro Tololo, we searched for variable stars in the southern globular cluster NGC 6584. We obtained images during 8 nights between 28 May and 6 July of 2011. After processing the images, we used the image subtraction package ISIS developed by Alard (2000)to search for the variable stars. We identified a total of 69 variable stars in our 10×10 arcmin^2 field, including 43 variables cataloged by Millis & Liller (1980) and 26 hereto unknown variables. In total, we classified 46 of the variables as type RRab, with a mean period of 0.56776 days, 15 as type RRc with a mean period of 0.30886 days, perhaps one lower amplitude type RRe, with a period of 0.26482 days, 4 eclipsing binaries, and 3 long period (P > 2 days) variable stars. As many as 15 of the RRab Lyrae stars exhibited the Blazhko Effect. Furthermore, the mean periods of the RR Lyrae types, the exhibited period/amplitude relationship, and the ratio of N_c/(N_ab+N_c) of 0.25 are consistent with an Oosterhoff Type I cluster. Here we present refined periods, V-band light curves, and classifications for each of the 69 variables, as well as a color-magnitude diagram of the cluster.

Thermal signatures of tether-cutting reconnections in pre-eruption coronal flux ropes: hot central voids in coronal cavities

Using a 3D MHD simulation, we model the quasi-static evolution and the onset of eruption of a coronal flux rope. The simulation begins with a twisted flux rope emerging at the lower boundary and pushing into a pre-existing coronal potential arcade field. At a chosen time the emergence is stopped with the lower boundary taken to be rigid. Then the coronal flux rope settles into a quasi-static rise phase with an underlying, central sigmoid-shaped current layer developing. Reconnections in the dissipating current layer during the quasi-static phase effectively add twisted flux to the flux rope and thus allow it to rise quasi-statically, even though the magnetic energy is decreasing as the system relaxes. We examine the thermal features produced by such “tether-cutting” reconnections as a result of heating and field aligned thermal conduction. It is found that a central hot, low-density channel containing reconnected, twisted flux threading under the flux rope axis forms on top of the central current layer. When viewed in the line of sight roughly aligned with the hot channel (which is roughly along the neutral line), the central current layer appears as a high-density vertical column with upward extensions as a “U” shaped dense shell enclosing a central hot, low-density void. Such thermal features have been observed within coronal prominence cavities. Our MHD simulations suggest that they are the signatures and consequences of the tether-cutting reconnections, and that the central void grows and rises with the reconnections, until the flux rope reaches the critical height for the onset of the torus instability and dynamic eruption ensues.

New SX Phoenicis Variables in the Globular Cluster NGC 4833

We report the discovery of 6 SX Phoenicis stars in the southern globular cluster NGC 4833. Images were obtained from January through June 2011 with the Southeastern Association for Research in Astronomy 0.6 meter telescope located at Cerro Tololo Interamerican Observatory. The image subtraction method of Alard & Lupton (1998)was used to search for variable stars in the cluster. We confirmed 17 previously cataloged variables by Demers & Wehlau(1977). In addition to the previously known variables we have identified 10 new variables. Of the total number of confirmed variables in our 10×10 arcmin^2 field, we classified 10 RRab variables, with a mean period of 0.69591 days, 7 RRc, with a mean period of 0.39555 days, 2 possible RRe variables with a mean period of 0.30950 days, a W Ursae Majoris contact binary, an Algol-type binary, and the 6 SX Phoenicis stars with a mean period of 0.05847 days. The periods, relative numbers of RRab and RRc variables, and Bailey diagram are indicative of the cluster being of the Oosterhoff type II. We present the phased-light curves, periods of previously known variables and the periods and classifications of the newly discovered variables, and their location on the color-magnitude diagram.

Spectral Typing of Late Type Stellar Companions to Young Stars from Low Dispersion Near-Infrared Integral Field Unit Data

We used the Project 1640 near-infrared coronagraph and integral field spectrograph to observe 19 young solar type stars. Five of these stars are known binary stars and we detected the late-type secondaries and were able to measure their JH spectra with a resolution of R\sim30. The reduced, extracted, and calibrated spectra were compared to template spectra from the IRTF spectral library. With this comparison we test the accuracy and consistency of spectral type determination with the low-resolution near-infrared spectra from P1640. Additionally, we determine effective temperature and surface gravity of the companions by fitting synthetic spectra calculated with the PHOENIX model atmosphere code. We also present several new epochs of astrometry of each of the systems. Together these data increase our knowledge and understanding of the stellar make up of these systems. In addition to the astronomical results, the analysis presented helps validate the Project 1640 data reduction and spectral extraction processes and the utility of low-resolution, near-infrared spectra for characterizing late-type companions in multiple systems.

Primary neutral helium in the heliosphere

Two years of neutral measurements by IBEX-Lo have yielded several direct observations of interstellar neutral helium and oxygen during preferred viewing seasons. Besides the interstellar signal, there are indications of the presence of secondary neutral helium and oxygen created in the heliosphere. Detailed modeling of these particle species is necessary to connect the measured fluxes to the pristine local interstellar medium while accounting for loss and production of neutral particles during their path through the heliosphere. In this contribution, global heliosphere models are coupled to analytic calculations of neutral trajectories to obtain detailed estimates of the neutral distribution function of primary interstellar helium atoms in the heliosphere, in particular in the inner heliosphere.

On Simulating Type Ia Supernovae

Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwarf star. We describe the multi-scale physical processes that realistic models must incorporate and the numerical models for these that we employ. In particular, we describe a flame-capturing scheme that addresses the problem of turbulent thermonuclear combustion on unresolved scales. We present the results of our study of the systematics of type Ia supernovae including trends in brightness following from properties of the host galaxy that agree with observations. We also present performance results from simulations on leadership-class architectures.

Surface Alfven waves in solar flux tubes

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. Alfven waves and magneto-sonic waves are particular classes of MHD waves. These wave modes are clearly different and have pure properties in uniform plasmas of infinite extent only. Due to plasma non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is a quantity unequivocally related to Alfven waves as compression is for magneto-sonic waves. Here, we investigate MHD waves superimposed on a one-dimensional non-uniform straight cylinder with constant magnetic field. For a piecewise constant density profile we find that the fundamental radial modes of the non-axisymmetric waves have the same properties as surface Alfven waves at a true discontinuity in density. Contrary to the classic Alfveen waves in a uniform plasma of infinite extent, vorticity is zero everywhere except at the cylinder boundary. If the discontinuity in density is replaced with a continuous variation of density, vorticity is spread out over the whole interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves do not need compression to exist unlike the radial overtones. In thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves with phase velocities between the internal and the external Alfven velocities can be considered as surface Alfven waves. On the contrary, the radial overtones can be related to fast-like magneto-sonic modes.

The first confirmed superoutburst of the dwarf nova GALEX J215818.5+241924

GALEX J215818.5+241924 was initially identified as a possible nova in Pegasus. We report unfiltered photometry of the object which revealed the presence of superhumps, with peak-to-peak amplitude of up to 0.22 magnitudes, diagnostic of it being a member of the SU UMa family of dwarf novae. The outburst amplitude was 4.3 magnitudes and it lasted at least 10 days, with a maximum brightness of magnitude 14.3. We determined the mean superhump period from our first 5 nights of observations as Psh = 0.06728(21) d. However analysis of the O-C residuals showed a dramatic evolution in Psh during the outburst. During the first part of the plateau phase the period increased with dPsh/dt = +2.67(15) \times 10-4. There was then an abrupt change following which the period decreased with dPsh/dt = -2.08(9)\times10-4. We found a signal in the power spectrum of the photometry which we interpret as the orbital signal with Porb = 0.06606(35) d. Thus the superhump period excess was epsilon = 0.020(8), such value being consistent with other SU UMa systems of similar orbital period.

New Bright Carbon Stars Found In The DFBS. (Research note, submitted to Astrophysics)

Carbon rich stars of Population II, such as CH giants, can provide direct information on the role of low to intermediate mass stars of the Halo in early Galactic evolution. Moreover, accurate knowledge of the CH stellar population is a critical requirement for building up scenarios for early Galactic chemical evolution. The first list of the faint high latitude C stars (FHLCs), found in the Digitized First Byurakan Survey {http://byurakan.phys.uniroma1.it and http://www.aras.am/Dfbs/dfbs.html} (DFBS) is given in a paper recently submitted for publication in Astronomy ad Astrophysics. In the present work, we report the recent discovery of two additional CH type stars (not previously catalogued), namely DFBS J075331.98+190344.3 and DFBS J111422.94+091442.7, detected on the DFBS plates with help of the image analysis softwares (FITSView and SAOImage DS9). Medium resolution spectra confirm the C rich nature for both of them. Using infrared color magnitudes relationship, we estimated the distances and K band absolute magnitudes of the new objects.

Thomas Hinsley Astbury: from an English market town schoolroom to the internal constitution of the stars [Cross-Listing]

T. H. Astbury (1858-1922) was for many years the much-respected headmaster of a boys’ junior school in the English market town of Wallingford. By night he was a dedicated amateur astronomer who enjoyed observing meteors, variable stars and many other objects. He began to search few new variable stars, his first discovery being the bright Cepheid variable, RT Aurigae. This, along with his discovery of 4 other variable stars, brought him to attention of some of the most famous professional astronomers of the age, including Herbert Hall Turner, Frank Dyson and Arthur Eddington.

Propagating Disturbances in Coronal Loops: A Detailed Analysis of Propagation Speeds [Replacement]

Quasi-periodic disturbances have been observed in the outer solar atmosphere for many years now. Although first interpreted as upflows (Schrijver et al. (1999)), they have been widely regarded as slow magnetoacoustic waves, due to observed velocities and periods. However, recent observations have questioned this interpretation, as periodic disturbances in Doppler velocity, line width and profile asymmetry were found to be in phase with the intensity oscillations (De Pontieu et al. (2010),Tian1 et al. (2011))}, suggesting the disturbances could be quasi-periodic upflows. Here we conduct a detailed analysis of the velocities of these disturbances across several wavelengths using the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We analysed 41 examples, including both sunspot and non sunspot regions of the Sun. We found that the velocities of propagating disturbances (PDs) located at sunspots are more likely to be temperature dependent, whereas the velocities of PDs at non sunspot locations do not show a clear temperature dependence. We also considered on what scale the underlying driver is affecting the properties of the PDs. Finally, we found that removing the contribution due to the cooler ions in the 193 A wavelength suggests that a substantial part of the 193 emission of sunspot PDs can be contributed to the cool component of 193\AA.

Propagating Disturbances in Coronal Loops: A Detailed Analysis of Propagation Speeds

Quasi-periodic disturbances have been observed in the outer solar atmosphere for many years now. Although first interpreted as upflows (Schrijver et al. (1999)), they have been widely regarded as slow magnetoacoustic waves, due to observed velocities and periods. However, recent observations have questioned this interpretation, as periodic disturbances in Doppler velocity, line width and profile asymmetry were found to be in phase with the intensity oscillations (De Pontieu et al. (2010),Tian1 et al. (2011))}, suggesting the disturbances could be quasi-periodic upflows. Here we conduct a detailed analysis of the velocities of these disturbances across several wavelengths using the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We analysed 41 examples, including both sunspot and non sunspot regions of the Sun. We found that the velocities of propagating disturbances (PDs) located at sunspots are more likely to be temperature dependent, whereas the velocities of PDs at non sunspot locations do not show a clear temperature dependence. We also considered on what scale the underlying driver is affecting the properties of the PDs. Finally, we found that removing the contribution due to the cooler ions in the 193 A wavelength suggests that a substantial part of the 193 emission of sunspot PDs can be contributed to the cool component of 193\AA.

Pure Neutron Matter Constraints on the Relativistic Mean-Field and Skyrme-Hartree-Fock Models [Cross-Listing]

Using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent \emph{ab initio} calculations, we present a general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs) while maintaining the broad quality of the predictions for binding energies and charge radii of nuclei. Such optimization leads to broadly consistent and tight predictions of the symmetry energy $J$ and its slope parameter $L$ at saturation density with associated joint 1$\sigma$ confidence ellipses in the $J-L$ plane. We demonstrate that given this optimization, the resulting neutron skin thicknesses are consistent with the experimental data so far, but clear model dependence shows up in (a) the slope of the correlation between $J$ and $L$ from the confidence ellipse, (b) the curvature parameter of the symmetry energy $K_{\rm sym}$, (c) the symmetry energy at supra-saturation densities, and (d) the neutron star radius. Notably, the RMF and SHF models are shown to be only marginally consistent at best with constraints on the isospin-dependent part of the incompressibility of neutron-rich nuclear matter $K_{\tau}$, and that the model dependence can lead to about 1 km difference of the neutron star radius given the same values of $J$, $L$ and symmetric nuclear matter (SNM) saturation properties.

Data Analysis of Gravitational Waves Signals from Millisecond Pulsars [Replacement]

The present work is devoted to the detection of monochromatic gravitational wave signals emitted by pulsars using ALLEGRO’s data detector. We will present the region (in frequency) of millisecond pulsars of the globular cluster 47 Tucanae (NGC 104) in the band of detector. With this result it was possible to analyse the data in the frequency ranges of the pulsars J1748-2446L and J1342+2822c, searching for annual Doppler variations using power spectrum estimates for the year 1999. We tested this method injecting a simulated signal in real data and we were able to detect it.

Data Analysis of Gravitational Waves Signals from Millisecond Pulsars [Cross-Listing]

The present work is devoted to the detection of monochromatic gravitational waves signals from pulsars using the ALLEGRO’s data detector. In this work we will present the region (in frequency) of millisecond pulsars of the 47 Tucanae (NGC 104) in the band of detector. According with this result was possible to analyse the frequency of the pulsar J1748-2446L and J1342+2822c, searching annual Doppler variations using power spectrum estimates for the year 1999. We tested this method injecting a simulated signal in real data and we were able to detect it.

Discovery of Two Very Wide Binaries with Ultracool Companions and a New Brown Dwarf at the L/T Transition

We present the discovery and spectroscopic follow-up of a nearby late-type L dwarf (2M0614+3950), and two extremely wide very-low-mass binary systems (2M0525-7425AB and 2M1348-1344AB), resulting from our search for common proper motion pairs containing ultracool components in the Two Micron All Sky Survey (2MASS) and the Wide-field Infrared Survey Explorer (WISE) catalogs. The near-infrared spectrum of 2M0614+3950 indicates a spectral type L9 \pm 1 object residing at a distance of 26.1 \pm 1.3 pc. The optical spectrum of the 2M0525-7425 primary reveals an M3.0 \pm 0.5 dwarf, accompanied by a secondary previously classified as L2. The system has an angular separation of ~44″, equivalent to ~2000 AU at the 45.7 \pm 2.5 pc distance. Using optical and infrared spectra, respectively, we classify the components of 2M1348-1344AB as M4.5 \pm 0.5 and T6 \pm 1. The angular separation of ~68″ is equivalent to ~1300 AU at the distance of 19.2 \pm 0.9 pc. 2M1348-1344AB is one of only five very wide (separation > 1000 AU) systems containing late T dwarfs known to date.

A Physical Model for SN 2001ay, a normal, bright, extremely slowly declining Type Ia supernova

We present a study of the peculiar Type Ia supernova 2001ay (SN 2001ay). The defining features of its peculiarity are: high velocity, broad lines, and a fast rising light curve, combined with the slowest known rate of decline. It is one magnitude dimmer than would be predicted from its observed value of Delta-m15, and shows broad spectral features. We base our analysis on detailed calculations for the explosion, light curves, and spectra. We demonstrate that consistency is key for both validating the models and probing the underlying physics. We show that this SN can be understood within the physics underlying the Delta-m15 relation, and in the framework of pulsating delayed detonation models originating from a Chandrasekhar mass, white dwarf, but with a progenitor core composed of 80% carbon. We suggest a possible scenario for stellar evolution which leads to such a progenitor. We show that the unusual light curve decline can be understood with the same physics as has been used to understand the Delta-m15 relation for normal SNe Ia. The decline relation can be explained by a combination of the temperature dependence of the opacity and excess or deficit of the peak luminosity, alpha, measured relative to the instantaneous rate of radiative decay energy generation. What differentiates SN 2001ay from normal SNe Ia is a higher explosion energy which leads to a shift of the Ni56 distribution towards higher velocity and alpha < 1. This result is responsible for the fast rise and slow decline. We define a class of SN 2001ay-like SNe Ia, which will show an anti-Phillips relation.

Precision Astrometry of the Exoplanet Host Candidate GD 66

The potential existence of a giant planet orbiting within a few AU of a stellar remnant has profound implications for both the survival and possible regeneration of planets during post-main sequence stellar evolution. This paper reports Hubble Space Telescope Fine Guidance Sensor and U.S. Naval Observatory relative astrometry of GD 66, a white dwarf thought to harbor a giant planet between 2 and 3 AU based on stellar pulsation arrival times. Combined with existing infrared data, the precision measurements here rule out all stellar-mass and brown dwarf companions, implying that only a planet remains plausible, if orbital motion is indeed the cause of the variations in pulsation timing.

Imaging the Algol Triple System in H Band with the CHARA Interferometer [Replacement]

Algol (Beta Per) is an extensively studied hierarchical triple system whose inner pair is a prototype semi-detached binary with mass transfer occurring from the sub-giant secondary to the main-sequence primary. We present here the results of our Algol observations made between 2006 and 2010 at the CHARA interferometer with the Michigan Infrared Combiner in the H band. The use of four telescopes with long baselines allows us to achieve better than 0.5 mas resolution and to unambiguously resolve the three stars. The inner and outer orbital elements, as well as the angular sizes and mass ratios for the three components are determined independently from previous studies. We report a significantly improved orbit for the inner stellar pair with the consequence of a 15% change in the primary mass compared to previous studies. We also determine the mutual inclination of the orbits to be much closer to perpendicularity than previously established. State-of-the-art image reconstruction algorithms are used to image the full triple system. In particular an image sequence of 55 distinct phases of the inner pair orbit is reconstructed, clearly showing the Roche-lobe-filling secondary revolving around the primary, with several epochs corresponding to the primary and secondary eclipses.

Imaging the Algol Triple System in H Band with the CHARA Interferometer

Algol (Beta Per) is an extensively studied hierarchical triple system whose inner pair is a prototype semi-detached binary with mass transfer occurring from the main-sequence primary to the sub-giant secondary. We present here the results of our Algol observations made between 2006 and 2010 at the CHARA interferometer with the Michigan Infrared Combiner in the H band. The use of four telescopes with long baselines allows us to achieve better than 0.5 mas resolution and to unambiguously resolve the three stars. The inner and outer orbital elements, as well as the angular sizes and mass ratios for the three components are determined independently from previous studies. We report a significantly improved orbit for the inner stellar pair with the consequence of a 15% change in the primary mass compared to previous studies. We also determine the mutual inclination of the orbits to be much closer to perpendicularity than previously established. State-of-the-art image reconstruction algorithms are used to image the full triple system. In particular an image sequence of 55 distinct phases of the inner pair orbit is reconstructed, clearly showing the Roche-lobe-filling secondary revolving around the primary, with several epochs corresponding to the primary and secondary eclipses.

Interpreting Magnetic Variance Anisotropy Measurements in the Solar Wind [Cross-Listing]

The magnetic variance anisotropy ($\mathcal{A}_m$) of the solar wind has been used widely as a method to identify the nature of solar wind turbulent fluctuations; however, a thorough discussion of the meaning and interpretation of the $\mathcal{A}_m$ has not appeared in the literature. This paper explores the implications and limitations of using the $\mathcal{A}_m$ as a method for constraining the solar wind fluctuation mode composition and presents a more informative method for interpreting spacecraft data. The paper also compares predictions of the $\mathcal{A}_m$ from linear theory to nonlinear turbulence simulations and solar wind measurements. In both cases, linear theory compares well and suggests the solar wind for the interval studied is dominantly Alfv\’{e}nic in the inertial and dissipation ranges to scales $k \rho_i \simeq 5$.

The Relationship Between the Expansion Speed and Radial Speed of CMEs Confirmed Using Quadrature Observations of the 2011 February 15 CME

It is difficult to measure the true speed of Earth-directed CMEs from a coronagraph located along the Sun-Earth line because of the occulting disk. However, the expansion speed (the speed with which the CME appears to spread in the sky plane) can be measured by such a coronagraph. In order to convert the expansion speed to radial speed (which is important for space weather applications) one can use an empirical relationship between the two that assumes an average width for all CMEs. If we have the width information from quadrature observations, we can confirm the relationship between expansion and radial speeds derived by Gopalswamy et al. (2009a). The STEREO spacecraft were in qudrature with SOHO (STEREO-A ahead of Earth by 87o and STEREO-B 94o behind Earth) on 2011 February 15, when a fast Earth-directed CME occurred. The CME was observed as a halo by the Large-Angle and Spectrometric Coronagraph (LASCO) on board SOHO. The sky-plane speed was measured by SOHO/LASCO as the expansion speed, while the radial speed was measured by STEREO-A and STEREO-B. In addition, STEREO-A and STEREO-B images provided the width of the CME, which is unknown from Earth view. From the SOHO and STEREO measurements, we confirm the relationship between the expansion speed (Vexp) and radial speed (Vrad) derived previously from geometrical considerations (Gopalswamy et al. 2009a): Vrad = 1/2 (1 + cot w)Vexp, where w is the half width of the CME. STEREO-B images of the CME, we found that CME had a full width of 76o, so w = 38o. This gives the relation as Vrad = 1.14 Vexp. From LASCO observations, we measured Vexp = 897 km/s, so we get the radial speed as 1023 km/s. Direct measurement of radial speed yields 945 km/s (STEREO-A) and 1058 km/s (STEREO-B). These numbers are different only by 7.6% and 3.4% (for STEREO-A and STEREO-B, respectively) from the computed value.

Properties of Ground Level Enhancement Events and the Associated Solar Eruptions during Solar Cycle 23

We present an overview of the observed properties of the GLEs and those of the associated flares and CMEs. The solar eruptions are very intense involving X-class flares and extreme CME speeds (average ~2000 km/s). The active regions in which the GLE events originate are generally large: 1290 msh (median 1010 msh) compared to 934 msh (median: 790 msh) for SEP-producing active regions. The initial acceleration of GLE-associated CMEs is much larger (by a factor of 2) than that of ordinary CMEs (2.3 km/s2 vs.1 km/s2). The GLE particle release is delayed with respect to the onset of all electromagnetic signatures of the eruptions: type II bursts, low frequency type III bursts, soft X-ray flares and CMEs. The presence of metric type II radio bursts some 17 min (median: 16 min; range: 3 to 48 min) before the GLE onset indicates shock formation well before the particle release. The release of GLE particles occurs when the CMEs reach an average height of ~3.09 Rs for well-connected events. For poorly connected events, the average CME height at GLE particle release is ~66% larger (mean: 5.18 Rs). The longitudinal dependence is consistent with shock accelerations because the shocks from poorly connected events need to expand more to cross the field lines connecting to an Earth observer. The CME height at metric type II burst onset is in the narrow range 1.29 to 1.8 Rs, with A mean of 1.53 Rs. The CME heights at metric type II burst onset and GLE particle release correspond to the minimum and maximum in the Alfven speed profile. The CME heights at GLE particle release are in good agreement with those obtained from the velocity dispersion analysis (Reames, 2009a,b) including the source longitude dependence. We also discuss the implications of the delay of GLE particle release with respect to complex type III bursts and hard X-ray emission.

Sensitivity of astrophysical reaction rates to nuclear uncertainties

Sensitivities of nuclear reaction rates to a variation of nuclear properties are studied. Target nuclei range from proton- to neutron-dripline for 10<=Z<=83. Reactions considered are nucleon- and alpha-induced reactions mediated by the strong interaction. The contribution of reactions proceeding on the target ground state to the total stellar rate is also given. General dependences on various input quantities are discussed. Additionally, sensitivities of laboratory cross sections of nucleon-, alpha-, and gamma-induced reactions are shown, allowing to estimate the impact of cross section measurements. Finally, recommended procedures to explore and improve reaction rate uncertainties using the present sensitivity data are outlined.

A unifying picture of helical and azimuthal MRI, and the universal significance of the Liu limit

The magnetorotational instability (MRI) plays a key role for cosmic structure formation by triggering turbulence in the rotating flows of accretion disks that would be otherwise hydrodynamically stable. In the limit of small magnetic Prandtl number, the helical and the azimuthal version of MRI are known to be governed by a quite different scaling behaviour than the standard MRI with a vertical applied magnetic field. Using the short-wavelength approximation for an incompressible, resistive, and viscous rotating fluid we present a unified description of helical and azimuthal MRI, and we identify the universal character of the Liu limit Ro= -0.8284 for the critical Rossby number. From this universal behaviour we are also lead to the prediction of higher azimuthal wavenumber for rather small ratios of azimuthal to axial applied fields.

Modeling the irregularities of solar cycle using flux transport dynamo models

The sunspot number varies roughly periodically with time. However the individual cycle durations and the amplitudes are found to vary in an irregular manner. It is observed that the stronger cycles are having shorter rise times and vice versa. This leads to an important effect know as the Waldmeier effect. Another important feature of the solar cycle irregularity are the grand minima during which the activity level is strongly reduced. We explore whether these solar cycle irregularities can be studied with the help of the flux transport dynamo model of the solar cycle. We show that with a suitable stochastic fluctuations in a regular dynamo model, we are able to reproduce many irregular features of the solar cycle including the Waldmeier effect and the grand minimum. However, we get all these results only if the value of the turbulent diffusivity in the convection zone is reasonably high.

First Simultaneous Detection of Moving Magnetic Features in Photospheric Intensity and Magnetic Field Data

The formation and the temporal evolution of a bipolar moving magnetic feature (MMF) was studied with high spatial and temporal resolution. The photometric properties were observed with the New Solar Telescope at Big Bear Solar Observatory using a broadband TiO filter (705.7 nm), while the magnetic field was analyzed using the spectropolarimetric data obtained by Hinode. For the first time, we observed a bipolar MMF simultaneously in intensity images and magnetic field data, and studied the details of its structure. The vector magnetic field and the Doppler velocity of the MMF were also studied. A bipolar MMF having its positive polarity closer to the negative penumbra formed being accompanied by a bright, filamentary structure in the TiO data connecting the MMF and a dark penumbral filament. A fast downflow (<2km/s) was detected at the positive polarity. The vector magnetic field obtained from the full Stokes inversion revealed that a bipolar MMF has a U-shaped magnetic field configuration. Our observations provide a clear intensity counterpart of the observed MMF in the photosphere, and strong evidence of the connection between the MMF and the penumbral filament as a serpentine field.

The locations of SNe Ib/c and their comparison to those of WR stars and GRBs

The locations of long GRBs and stripped supernovae are compared to those of their favored progenitors, WR stars, and their sub-classes. Compared to Leloudas et al. (2010), we have doubled the number of galaxies with suitable WR data. In the combined sample, WC stars are found, on average, in brighter locations than WN stars. The WN distribution is fully consistent with the one of SNe Ib, while it is inconsistent with those of SNe II, Ic and GRBs. The WC distribution is both consistent with SNe Ib and Ic. It is inconsistent with SNe II, and marginally consistent with GRBs. Furthermore, we present a spectroscopic study of the locations of SNe Ib/c. The average metallicity in the environments of SNe Ic is found to be a little higher than for SNe Ib, but the difference is small and not significant within our sample. Under the assumption that the SN regions were formed in an instantaneous burst of star formation, we find that a fraction of them appear older than what is allowed in order to host SNe Ib/c from single massive stars. Within this framework, these SNe must come from lower mass binaries.

Detection of gravity modes in the massive binary V380 Cyg from Kepler spacebased photometry and high-resolution spectroscopy

We report the discovery of low-amplitude gravity-mode oscillations in the massive binary star V380 Cyg, from 180 d of Kepler custom-aperture space photometry and 5 months of high-resolution high signal-to-noise spectroscopy. The new data are of unprecedented quality and allowed to improve the orbital and fundamental parameters for this binary. The orbital solution was subtracted from the photometric data and led to the detection of periodic intrinsic variability with frequencies of which some are multiples of the orbital frequency and others are not. Spectral disentangling allowed the detection of line-profile variability in the primary. With our discovery of intrinsic variability interpreted as gravity mode oscillations, V380 Cyg becomes an important laboratory for future seismic tuning of the near-core physics in massive B-type stars.

Solving the mode identification problem in asteroseismology of F stars observed with Kepler

Asteroseismology of F-type stars has been hindered by an ambiguity in identification of their oscillation modes. The regular mode pattern that makes this task trivial in cooler stars is masked by increased linewidths. The absolute mode frequencies, encapsulated in the asteroseismic variable epsilon, can help solve this impasse because the values of epsilon implied by the two possible mode identifications are distinct. We find that the correct epsilon can be deduced from the effective temperature and the linewidths and we apply these methods to a sample of solar-like oscillators observed with Kepler.

The Coronal Abundance Anomalies of M Dwarfs

We analyze Chandra X-ray spectra of the M0 V+M0 V binary GJ 338. As quantified by X-ray surface flux, these are the most inactive M dwarfs ever observed with X-ray grating spectroscopy. We focus on measuring coronal abundances, in particular searching for evidence of abundance anomalies related to First Ionization Potential (FIP). In the solar corona and wind, low FIP elements are overabundant, which is the so-called “FIP effect.” For other stars, particularly very active ones, an “inverse FIP effect” is often observed, with low FIP elements being underabundant. For both members of the GJ 338 binary, we find evidence for a modest inverse FIP effect, consistent with expectations from a previously reported correlation between spectral type and FIP bias. This amounts to strong evidence that all M dwarfs should exhibit the inverse FIP effect phenomenon, not just the active ones. We take the first step towards modeling the inverse FIP phenomenon in M dwarfs, building on past work that has demonstrated that MHD waves coursing through coronal loops can lead to a ponderomotive force that fractionates elements in a manner consistent with the FIP effect. We demonstrate that in certain circumstances this model can also lead to an inverse FIP effect, pointing the way to more detailed modeling of M dwarf coronal abundances in the future.

Advances in understanding young high-mass stars using optical interferometry

The closest examples of high-mass star birth occurs in deeply embedded environments at kiloparsec distances. Although much progress has been made, an observationally validated picture of the dominant processes which allows the central hydrostatic object to grow in mass has yet to be established. The observational technique of optical interferometry has demonstrated its potential in the field of high-mass star formation by delivering a milli-arcsecond infrared view on the complex accretion environment. We provide an overview of the scientific results obtained with multi-aperture telescope arrays and briefly discuss future instruments and their anticipated impact on our understanding of massive young stellar objects.

The transition mass-loss rate: Calibrating the role of line-driven winds in massive star evolution

A debate has arisen regarding the importance of stationary versus eruptive mass loss for massive star evolution. The reason is that stellar winds have been found to be clumped, which results in the reduction of unclumped empirical mass-loss rates. Most stellar evolution models employ theoretical mass-loss rates which are already reduced by a moderate factor of ~2-3 compared to non-corrected empirical rates. A key question is whether these reduced rates are of the correct order of magnitude, or if they should be reduced even further, which would mean that the alternative of eruptive mass loss becomes necessary. Here we introduce the transition mass-loss rate (dM/dt)_trans between O and Wolf-Rayet (WR) stars. Its novelty is that it is model independent. All that is required is postulating the spectroscopic transition point in a given data-set, and determining the stellar luminosity, which is far less model dependent than the mass-loss rate. The transition mass-loss rate is subsequently used to calibrate stellar wind strength by its application to the Of/WNh stars in the Arches cluster. Good agreement is found with two alternative modelling/theoretical results, suggesting that the rates provided by current theoretical models are of the right order of magnitude in the ~50Msun mass range. Our results do not confirm the specific need for eruptive mass loss as Luminous Blue Variables, and current stellar evolution modelling for Galactic massive stars seems sound. Mass loss through alternative mechanisms might still become necessary at lower masses, and/or metallicities, and the quantification of alternative mass loss is desirable.

SN 2009kn - The Twin of the Type IIn SN 1994W

We present an optical and near-infrared photometric and spectroscopic study of supernova (SN) 2009kn spanning ~1.5 yr from the discovery. The optical spectra are dominated by the narrow (FWHM ~1000 km s^-1) Balmer lines distinctive of a Type IIn SN with P-Cygni profiles. Contrarily the photometric evolution resembles more that of a Type IIP SN with a large drop in luminosity at the end of the plateau phase. These characteristics are similar to those of SN 1994W, whose nature has been explained with two different models with different approaches. The well-sampled dataset on SN 2009kn offers the possibility to test these models, both in the case of SN 2009kn and SN 1994W. We associate the narrow P-Cygni lines with a swept-up shell composed of circumstellar matter and SN ejecta. The broad emission line wings, seen during the plateau phase arise from internal electron scattering in this shell. The slope of the light curve after the post-plateau drop is fairly consistent with that expected from the radioactive decay of 56Co, suggesting a supernova origin for SN 2009kn. Assuming radioactivity to be the main source powering the light curve of SN 2009kn in the tail phase, we infer an upper limit for 56Ni mass of 0.023 M_sun. This is significantly higher than that estimated for SN 1994W, which also showed a much steeper decline of the light curve after the post-plateau drop. We also observe late-time near-infrared emission which most likely arises from newly-formed dust produced by SN 2009kn. As with SN 1994W, no broad lines are observed in the spectra of SN 2009kn, not even in the late time tail phase.

Effects of Thomson-Scattering Geometry on White-Light Imaging of an Interplanetary Shock: Synthetic Observations from Forward Magnetohydrodynamic Modelling

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time-elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east-west asymmetry in its brightness and degree of polarization. Therefore, measurement of the interplanetary polarized brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

The Mass-Loss Return From Evolved Stars to The Large Magellanic Cloud VI: Luminosities and Mass-Loss Rates on Population Scales

We present results from the first application of the Grid of Red Supergiant and Asymptotic Giant Branch ModelS (GRAMS) model grid to the entire evolved stellar population of the Large Magellanic Cloud (LMC). GRAMS is a pre-computed grid of 80,843 radiative transfer (RT) models of evolved stars and circumstellar dust shells composed of either silicate or carbonaceous dust. We fit GRAMS models to ~30,000 Asymptotic Giant Branch (AGB) and Red Supergiant (RSG) stars in the LMC, using 12 bands of photometry from the optical to the mid-infrared. Our published dataset consists of thousands of evolved stars with individually determined evolutionary parameters such as luminosity and mass-loss rate. The GRAMS grid has a greater than 80% accuracy rate discriminating between Oxygen- and Carbon-rich chemistry. The global dust injection rate to the interstellar medium (ISM) of the LMC from RSGs and AGB stars is on the order of 1.5×10^(-5) solar masses/yr, equivalent to a total mass injection rate (including the gas) into the ISM of ~5×10^(-3) solar masses/yr. Carbon stars inject two and a half times as much dust into the ISM as do O-rich AGB stars, but the same amount of mass. We determine a bolometric correction factor for C-rich AGB stars in the K band as a function of J – K color, BC(K) = -0.40(J-K)^2 + 1.83(J-K) + 1.29. We determine several IR color proxies for the dust mass-loss rate (MLR) from C-rich AGB stars, such as log (MLR) = (-18.90)/((K-[8.0])+3.37)-5.93. We find that a larger fraction of AGB stars exhibiting the `long-secondary period’ phenomenon are O-rich than stars dominated by radial pulsations, and AGB stars without detectable mass-loss do not appear on either the first-overtone or fundamental-mode pulsation sequences.

A geometrical model for polar plumes observed by STEREO SECCH

Solar polar plumes are bright radial rays rooted at the sun’s polar areas. They are widely believed to have the structure of expanding tube. A four degree polynomial function was assumed to represent the change of cross section diameter as a function of height and four unknown parameters were calculated with measured average widths of total 31 plumes at 4 heights (1.04$R_s$,1.10$R_s$,1.16$R_s$,1.20$R_s$).

The SDSS-HET Survey of Kepler Eclipsing Binaries: Spectroscopic Dynamical Masses of the Kepler-16 Circumbinary Planet Hosts

We have used high-resolution spectroscopy to observe the Kepler-16 eclipsing binary as a double-lined system, and measure precise radial velocities for both stellar components. These velocities yield a dynamical mass-ratio of q=0.2994+-0.0031. When combined with the inclination, i=90.3401+0.0016-0.0019 deg, measured from the Kepler photometric data by Doyle et al. 2011, we derive dynamical masses for the Kepler-16 components of M_A=0.654+-0.017 M_sun and M_B=0.1959+-0.0031 M_sun, a precision of 2.5% and 1.5% respectively. Our results confirm at the ~2% level the mass-ratio derived by Doyle et al. with their photometric-dynamical model, q=0.2937+-0.0006. These are among the most precise spectroscopic dynamical masses ever measured for low-mass stars, and provide an important direct test of the results from the photometric-dynamical modeling technique.

Observing Simulated Protostars with Outflows: How Accurate are Protostellar Properties Inferred from SEDs?

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. We perform synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte-Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against over-interpreting properties inferred from SEDs for unresolved protostellar sources. (Abridged)

Analysis of Gamma Radiation from a Radon Source: Indications of a Solar Influence

This article presents an analysis of about 29,000 measurements of gamma radiation associated with the decay of radon in a sealed container at the Geological Survey of Israel (GSI) Laboratory in Jerusalem between 28 January 2007 and 10 May 2010. These measurements exhibit strong variations in time of year and time of day, which may be due in part to environmental influences. However, time-series analysis reveals a number of periodicities, including two at approximately 11.2 year$^{-1}$ and 12.5 year$^{-1}$. We have previously found these oscillations in nuclear-decay data acquired at the Brookhaven National Laboratory (BNL) and at the Physikalisch-Technische Bundesanstalt (PTB), and we have suggested that these oscillations are attributable to some form of solar radiation that has its origin in the deep solar interior. A curious property of the GSI data is that the annual oscillation is much stronger in daytime data than in nighttime data, but the opposite is true for all other oscillations. This may be a systematic effect but, if it is not, this property should help narrow the theoretical options for the mechanism responsible for decay-rate variability.

International observational campaigns of the last two eclipses in EE Cephei: 2003 and 2008/9

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations in low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete observational data which were collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous to a “bump” that is very clear in all the UBVRI light curves. Variations of the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and H_alpha emission is superimposed with the same strong absorption. These observations show that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the “bump” observed during the last two eclipses in terms of different times of obscuration of the hot polar regions of the Be star by the disk; (iii) our suggested period of the disk precession ~11-12 P_orb and predicted depth of about 2 mag the forthcoming eclipse in 2014.

Anchoring the Distance Scale via X-ray/IR Data for Cepheid Clusters: SU Cas

New X-ray (XMM-Newton) and JHKs (OMM) observations for members of the star cluster Alessi 95, which Turner et al.(2012) discovered hosts the classical Cepheid SU Cas, were used in tandem with UCAC3 (proper motion) and 2MASS observations to determine precise cluster parameters: E(J-H)=0.08+-0.02 and d=405+-15 pc. The ensuing consensus among cluster, pulsation, IUE, and trigonometric distances (d=414+-5(se)+-10(sd) pc) places SU Cas in a select group of nearby fundamental Cepheid calibrators (Delta Cep, Zeta Gem). High-resolution X-ray observations may be employed to expand that sample as the data proved pertinent for identifying numerous stars associated with SU Cas. Acquiring X-ray observations of additional fields may foster efforts to refine Cepheid calibrations used to constrain H_0.

Scars of Intense Accretion Episodes at Metal-Rich White Dwarfs

A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points to historical, time-averaged rates significantly higher than the currently observed episodes at their DAZ counterparts. The difference between the ongoing, instantaneous accretion rates witnessed at DAZ white dwarfs, which often exceed 1e8 g/s, and those inferred over the past 1e5-1e6 yr for the DBZ stars can be a few orders of magnitude, and therefore must result from high-rate episodes of tens to hundreds of years so they remain undetected to date. This paper explores the likelihood that such brief, intense accretion episodes of gas-phase material can account for existing data. For reasonable assumptions about the circumstellar gas, accretion rates approaching or exceeding 1e15 g/s are possible, similar to rates observed in quiescent cataclysmic variables, and potentially detectable with future x-ray missions or wide-field monitoring facilities. Gaseous debris that is prone to such rapid accretion may be abundant immediately following a tidal disruption event via collisions and sublimation, or if additional bodies impinge upon an extant disk. Particulate disk matter accretes at or near the Poynting-Robertson drag rate for long periods between gas-producing events, consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ stars without infrared excesses have rates consistent with accretion from particulate disks that remain undetected. This overall picture has implications for quasi-steady state models of accretion and the derived chemical composition of asteroidal debris in DBZ white dwarfs.

Observational Study Of the Quasi-Periodic Fast Propagating Magnetosonic Waves and the Associated Flare on 2011 May 30 [Replacement]

On 2011 May 30, quasi-periodic fast propagating (QFP) magnetosonic waves accompanied by a C2.8 flare were directly imaged by the Atomospheric Imaging Assembly instrument on board the Solar Dynamics Observatory. The QFP waves successively emanated from the flare kernel, they propagated along a cluster of open coronal loops with a phase speed of 834 km/s during the flare’s rising phase, and the multiple arc-shaped wave trains can be fitted with a series of concentric circles. We generate the k-omega diagram of the Fourier power and find a straight ridge that represents the dispersion relation of the waves. Along the ridge, we find a lot of prominent nodes which represent the available frequencies of the QFP waves. On the other hand, the frequencies of the flare are also obtained by analyzing the flare light curves using the wavelet technique. The results indicate that almost all the main frequencies of the flare are consistent with those of the QFP waves. This suggests that the flare and the QFP waves were possibly excited by a common physical origin. On the other hand, a few low frequencies revealed by the k-omega diagram can not be found in the accompanying flare. We propose that these low frequencies were possibly due to the leakage of the pressure-driven p-mode oscillations from the photosphere into the low corona, which should be a noticeable mechanism for driving the QFP waves observed in the corona.

Observational Study Of the Quasi-Periodic Fast Propagating Magnetosonic Waves and the Associated Flare on 2011 May 30

On 2011 May 30, quasi-periodic fast propagating (QFP) magnetosonic waves accompanied by a C2.8 flare were directly imaged by the Atomospheric Imaging Assembly instrument on board the Solar Dynamics Observatory. The QFP waves successively emanated from the flare kernel, they propagated along a cluster of open coronal loops with a phase speed of 834 km/s during the flare’s rising phase, and the multiple arc-shaped wave trains can be fitted with a series of concentric circles. We generate the k-omega diagram of the Fourier power and find a straight ridge that represents the dispersion relation of the waves. Along the ridge, we find a lot of prominent nodes which represent the available frequencies of the QFP waves. On the other hand, the frequencies of the flare are also obtained by analyzing the flare light curves using the wavelet technique. The results indicate that almost all the main frequencies of the flare are consistent with those of the QFP waves. This suggests that the flare and the QFP waves were possibly excited by a common physical origin. On the other hand, a few low frequencies revealed by the k-omega diagram can not be found in the accompanying flare. We propose that these low frequencies were possibly due to the leakage of the pressure-driven p-mode oscillations from the photosphere into the low corona, which should be a noticeable mechanism for driving the QFP waves observed in the corona.

The Diagnostic Potential of Transition Region Lines under-going Transient Ionization in Dynamic Events

We discuss the diagnostic potential of high cadence ultraviolet spectral data when transient ionization is considered. For this we use high cadence UV spectra taken during the impulsive phase of a solar flares (observed with instruments on-board the Solar Maximum Mission) which showed excellent correspondence with hard X-ray pulses. The ionization fraction of the transition region ion O V and in particular the contribution function for the O V 1371A line are computed within the Atomic Data and Analysis Structure, which is a collection of fundamental and derived atomic data and codes which manipulate them. Due to transient ionization, the O V 1371A line is enhanced in the first fraction of a second with the peak in the line contribution function occurring initially at a higher electron temperature than in ionization equilibrium. The rise time and enhancement factor depend mostly on the electron density. The fractional increase in the O V 1371A emissivity due to transient ionization can reach a factor of 2–4 and can explain the fast response in the line flux of transition regions ions during the impulsive phase of flares solely as a result of transient ionization. This technique can be used to diagnostic the electron temperature and density of solar flares observed with the forth-coming Interface Region Imaging Spectrograph.

On the energy release in solar flares

High-resolution observations show the fine structure of the global equilibrium magnetic field configuration in solar atmosphere to be essentially different from that assumed in the traditional ‘potential + force-free’ field scenarios. The interacting large-scale structures of fine field elements are separated by numerous non-force-free elements (tangential discontinuities) which are neglected in the traditional field picture. An incorporation of these elements into the model implies a dynamical rather than statical character of equilibrium of the field configuration. A transition of the system into flaring can be triggered by the ballooning mode of flute instability of prominences or/and coronal condensations. Tearing-mode and MHD instabilities as well as the effects of overheating of the turbulent current sheet prevent the field from stationary reconnection as it is adopted in the traditional scenario. We speculate around the assumption that the energy release in active regions is governed by the same scenario as dynamical current percolation through a random network of resistors in which the saltatory resistance is controlled by local current.

Active Region Moss: Doppler Shifts from Hinode/EIS Observations

Studying the Doppler shifts and the temperature dependence of Doppler shifts in moss regions can help us understand the heating processes in the core of the active regions. In this paper we have used an active region observation recorded by the Extreme-ultraviolet Imaging Spectrometer (EIS) onboard Hinode on 12-Dec-2007 to measure the Doppler shifts in the moss regions. We have distinguished the moss regions from the rest of the active region by defining a low density cut-off as derived by Tripathi et al. (2010). We have carried out a very careful analysis of the EIS wavelength calibration based on the method described in Young et al. (2012). For spectral lines having maximum sensitivity between log T = 5.85 and log T = 6.25 K, we find that the velocity distribution peaks at around 0 km/s with an estimated error of 4-5 km/s. The width of the distribution decreases with temperature. The mean of the distribution shows a blue shift which increases with increasing temperature and the distribution also shows asymmetries towards blue-shift. Comparing these results with observables predicted from different coronal heating models, we find that these results are consistent with both steady and impulsive heating scenarios. However, the fact that there are a significant number of pixels showing velocity amplitudes that exceed the uncertainty of 5 km s$^{-1}$ is suggestive of impulsive heating. Clearly, further observational constraints are needed to distinguish between these two heating scenarios.

Exploring the Morphology of RAVE Stellar Spectra

The RAdial Velocity Experiment (RAVE) is a medium resolution R~7500 spectroscopic survey of the Milky Way which already obtained over half a million stellar spectra. They present a randomly selected magnitude-limited sample, so it is important to use a reliable and automated classification scheme which identifies normal single stars and discovers different types of peculiar stars. To this end we present a morphological classification of 350,000 RAVE survey stellar spectra using locally linear embedding, a dimensionality reduction method which enables representing the complex spectral morphology in a low dimensional projected space while still preserving the properties of the local neighborhoods of spectra. We find that the majority of all spectra in the database ~90-95% belong to normal single stars, but there is also a significant population of several types of peculiars. Among them the most populated groups are those of various types of spectroscopic binary and chromospherically active stars. Both of them include several thousands of spectra. Particularly the latter group offers significant further investigation opportunities since activity of stars is a known proxy of stellar ages. Applying the same classification procedure to the sample of normal single stars alone shows that the shape of the projected manifold in two dimensional space correlates with stellar temperature, surface gravity and metallicity.

STEREO-A and PROBA2 Quadrature Observations of Reflections of three EUV Waves from a Coronal Hole

We investigate the interaction of three consecutive large-scale coronal waves with a polar coronal hole, simultaneously observed on-disk by the Solar TErrestrial Relations Observatory (STEREO)-A spacecraft and on the limb by the PRoject for On-Board Autonomy 2 (PROBA2) spacecraft on January 27, 2011. All three extreme-ultraviolet(EUV) waves originate from the same active region NOAA 11149 positioned at N30E15 in the STEREO-A field-of-view and on the limb in PROBA2. We derive for the three primary EUV waves start velocities in the range of ~310 km/s for the weakest up to ~500 km/s for the strongest event. Each large-scale wave is reflected at the border of the extended coronal hole at the southern polar region. The average velocities of the reflected waves are found to be smaller than the mean velocities of their associated direct waves. However, the kinematical study also reveals that in each case the end velocity of the primary wave matches the initial velocity of the reflected wave. In all three events the primary and reflected waves obey the Huygens-Fresnel principle, as the incident angle with ~10{\deg} to the normal is of the same size as the angle of reflection. The correlation between the speed and the strength of the primary EUV waves, the homologous appearance of both the primary and the reflected waves, and in particular the EUV wave reflections themselves implicate that the observed EUV transients are indeed nonlinear large-amplitude MHD waves.

Matching the Spectral Energy Distribution and p Mode Oscillation Frequencies of the Rapidly Rotating Delta Scuti Star ? Ophiuchi with a 2D Rotating Stellar Model

Spectral energy distributions are computed using 2D rotating stellar models and NLTE plane parallel model atmospheres. A rotating, 2D stellar model has been found which matches the observed ultraviolet and visible spectrum of ? Oph. The SED match occurs for the interferometrically deduced surface shape and inclination, and is different from the SED produced by spherical models. The p mode oscillation frequencies in which the latitudinal variation is modelled by a linear combination of eight Legendre polynomials were computed for this model. The five highest and seven of the nine highest amplitude modes show agreement between computed axisymmetric, equatorially symmetric mode frequencies and the mode frequencies observed by MOST to within the observational error. Including nonaxisymmetric modes up through |m| = 2 and allowing the possibility that the eight lowest amplitude modes could be produced by modes which are not equatorially symmetric produces matches for 24 out of the 35 MOST modes to within the observational error and another eight modes to within twice the observational error. The remaining three observed modes can be fit within 4.2 times the observational error, but even these may be fit to within the observational error if the criteria for computed modes are expanded.

Chemical segregation toward massive hot cores: The AFGL2591 star forming region

We present high angular resolution observations (0.5″x0.3″) carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

Discovery of Three Distant, Cold Brown Dwarfs in the WFC3 Infrared Spectroscopic Parallels Survey

We present the discovery of three late type (>T4) brown dwarfs, including a probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. We use the G141 grism spectra to determine the spectral types of the dwarfs and derive distance estimates based on a comparison with nearby T dwarfs with known parallaxes. These are the most distant spectroscopically confirmed T/Y dwarfs, with the farthest at an estimated distance of ~400 pc. We compare the number of cold dwarfs found in the WISP survey with simulations of the brown dwarf mass function. The number found is generally consistent with an initial stellar mass function dN/dM \propto M^{-\alpha} with \alpha = 0.0–0.5, although the identification of a Y dwarf is somewhat surprising and may be indicative of either a flatter absolute magnitude/spectral type relation than previously reported or an upturn in the number of very late type brown dwarfs in the observed volume.

The M16 molecular complex under the influence of NGC6611. Herschel's perspective of the heating effect on the Eagle Nebula

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70{\mu}m, 160{\mu}m, 250{\mu}m, 350{\mu}m, 500{\mu}m. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 – 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

Fast Differential Emission Measure Inversion of Solar Coronal Data

We present a fast method for reconstructing Differential Emission Measures (DEMs) using solar coronal data. The method averages over 1000 DEMs per second for an example active region observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO), and achieves reduced chi-squared of order unity with no negative emission in all but a few test cases. The high performance of this method is especially relevant in the context of AIA, which images of order one million solar pixels per second. This paper describes the method, analyzes its fidelity, compares its performance and results with other DEM methods, and applies it to an active region and loop observed by AIA and by the Extreme-ultraviolet Imaging Spectrometer (EIS) on Hinode.

John Goodricke, Edward Pigott, and Their Study of Variable Stars [Cross-Listing]

John Goodricke and Edward Pigott, working in York, England, between 1781 and 1786, determined the periods of variation of eclipsing binaries such as Algol and Beta Lyrae and speculated that the eclipses of Algol might be caused by a “dark body,” perhaps even a planet. They also determined the periods of variation of the first two known Cepheid variables, the stars whose period-luminosity relation today enables astronomers to determine distances to distant galaxies. Goodricke holds special interest because he was completely deaf and because he died at the age of 21. The lives and work of these two astronomers are described.

Electron-ion Recombination of Fe XII forming Fe XI: Laboratory Measurements and Theoretical Calculations

We have measured electron-ion recombination for Fe XII forming Fe XI using a merged beams configuration at the heavy-ion storage ring TSR located at the Max Planck Institute for Nuclear Physics in Heidelberg, Germany. The measured merged beams recombination rate coefficient (MBRRC) for collision energies from 0 to 1500 eV is presented. This work uses a new method for determining the absolute MBRRC based on a comparison of the ion beam decay rate with and without the electron beam on. For energies below 75 eV, the spectrum is dominated by dielectronic recombination (DR) resonances associated with 3s-3p and 3p-3d core excitations. At higher energies we observe contributions from 3-N’ and 2-N’ core excitations DR. We compare our experimental results to state-of-the-art multi-configuration Breit-Pauli (MCBP) calculations and find significant differences, both in resonance energies and strengths. We have extracted the DR contributions from the measured MBRRC data and transformed them into a plasma recombination rate coefficient (PRRC) for temperatures in the range of 10^3 to 10^7 K. We show that the previously recommended DR data for Fe XII significantly underestimate the PRRC at temperatures relevant for both photoionized plasmas (PPs) and collisionaly ionized plasmas (CPs). This is to be contrasted with our MCBP PRRC results which agree with the experiment to within 30% at PP temperatures and even better at CP temperatures. We find this agreement despite the disagreement shown by the detailed comparison between our MCBP and experimental MBRRC results. Lastly, we present a simple parameterized form of the experimentally derived PRRC for easy use in astrophysical modelling codes.

Did the ancient egyptians record the period of the eclipsing binary Algol - the Raging one?

The eclipses in binary stars give precise information of orbital period changes. Goodricke discovered the 2.867 days period in the eclipses of Algol in the year 1783. The irregular orbital period changes of this longest known eclipsing binary continue to puzzle astronomers. The mass transfer between the two members of this binary should cause a long-term increase of the orbital period, but observations over two centuries have not confirmed this effect. Here, we present evidence indicating that the period of Algol was 2.850 days three millenia ago. For religious reasons, the ancient Egyptians have recorded this period into the Cairo Calendar, which describes the repetitive changes of the Raging one. Cairo Calendar may be the oldest preserved historical document of the discovery of a variable star.

The scaler mode in the Pierre Auger Observatory to study heliospheric modulation of cosmic rays

The impact of the solar activity on the heliosphere has a strong influence on the modulation of the flux of low energy galactic cosmic rays arriving at Earth. Different instruments, such as neutron monitors or muon detectors, have been recording the variability of the cosmic ray flux at ground level for several decades. Although the Pierre Auger Observatory was designed to observe cosmic rays at the highest energies, it also records the count rates of low energy secondary particles (the scaler mode) for the self-calibration of its surface detector array. From observations using the scaler mode at the Pierre Auger Observatory, modulation of galactic cosmic rays due to solar transient activity has been observed (e.g., Forbush decreases). Due to the high total count rate coming from the combined area of its detectors, the Pierre Auger Observatory (its detectors have a total area greater than $16\,000$\,m$^2$) detects a flux of secondary particles of the order of $\sim 10^8$\,counts per minute. Time variations of the cosmic ray flux related to the activity of the heliosphere can be determined with high accuracy. In this paper we briefly describe the scaler mode and analyze a Forbush decrease together with the interplanetary coronal mass ejection that originated it. The Auger scaler data are now publicly available.

Comments on "Experimental study of radiative shocks at PALS facility" by C.Stehl\'e et al. (arXiv:1003.2739v1 [astro-ph.SR]) [Replacement]

Discussion on paper “Experimental study of radiative shocks at PALS facility” by C.Stehle\’et al (arXiv:1003.2739v1 [astro-ph.SR]) is presented, findings are questioned.

Comments on "Experimental study of radiative shocks at PALS facility" by C.Stehl\'e et al. (arXiv:1003.2739v1 [astro-ph.SR])

Discussion on paper “Experimental study of radiative shocks at PALS facility” by C.Stehle\’et al (arXiv:1003.2739v1 [astro-ph.SR]) is presented, findings are questioned.

Optimal Interpolation and Prediction in Pulsar Timing

For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.

Optical flares and flaring oscillations on the M-type eclipsing binary CU Cnc

We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.

A public database for white dwarf asteroseismology with fully evolutionary models: I. Chemical profiles and pulsation periods of ZZ Ceti (DAV) stars

We present a large bank of chemical profiles and pulsation periods suited for asteroseismological studies of ZZ Ceti (or DAV) variable stars. Our background equilibrium DA white dwarf models are the result of fully evolutionary computations that take into account the complete history of the progenitor stars from the ZAMS. The models are characterized by self-consistent chemical structures from the centre to the surface, and cover a wide range of stellar masses, effective temperatures, and H envelope thicknesses. We present dipole and quadrupole pulsation $g$-mode periods comfortably covering the interval of periods observed in ZZ Ceti stars. Complete tabulations of chemical profiles and pulsation periods to be used in asteroseismological period fits, as well as other quantities of interest, can be freely downloaded from our website (\url{http://www.fcaglp.unlp.edu.ar/evolgroup}

Search for ionized jets towards high-mass young stellar objects

We are carrying out multi-frequency radio continuum observations, using the Australia Telescope Compact Array, to systematically search for collimated ionized jets towards high-mass young stellar objects (HMYSOs). Here we report observations at 1.4, 2.4, 4.8 and 8.6 GHz, made with angular resolutions of about 7, 4, 2, and 1 arcsec, respectively, towards six objects of a sample of 33 southern HMYSOs thought to be in very early stages of evolution. The objects in the sample were selected from radio and infrared catalogs by having positive radio spectral indices and being luminous (L_bol > 20,000 L_sun), but underluminous in radio emission compared to that expected from its bolometric luminosity. This criteria makes the radio sources good candidates for being ionized jets. As part of this systematic search, two ionized jets have been discovered: one previously published and the other reported here. The rest of the observed candidates correspond to three hypercompact hii regions and two ultracompact hii regions. The two jets discovered are associated with two of the most luminous (70,000 and 100,000 Lsun) HMYSOs known to harbor this type of objects, showing that the phenomena of collimated ionized winds appears in the formation process of stars at least up to masses of ~ 20 M_sun and provides strong evidence for a disk-mediated accretion scenario for the formation of high-mass stars. From the incidence of jets in our sample, we estimate that the jet phase in high-mass protostars lasts for 40,000 yr.

The search for habitable worlds: 1. The viability of a starshade mission

As part of NASA’s mission to explore habitable planets orbiting nearby stars, this paper explores the detection and characterization capabilities of a 4-m space telescope plus 50-m starshade located at the Earth-Sun L2 point, a.k.a. the New Worlds Observer (NWO). Our calculations include the true spectral types and distribution of stars on the sky, an iterative target selection protocol designed to maximize efficiency based on prior detections, and realistic mission constraints. We carry out both analytical calculations and simulated observing runs for a wide range in exozodiacal background levels ({\epsilon} = 1 – 100 times the local zodi brightness) and overall prevalence of Earth-like terrestrial planets ({\eta}\oplus = 0.1 – 1). We find that even without any return visits, the NWO baseline architecture (IWA = 65 mas, limiting FPB = 4\times10-11) can achieve a 95% probability of detecting and spectrally characterizing at least one habitable Earth-like planet, and an expectation value of ~3 planets found, within the mission lifetime and {\Delta}V budgets, even in the worst-case scenario ({\eta}\oplus = 0.1 and {\epsilon} = 100 zodis for every target). This achievement requires about one year of integration time spread over the 5 year mission, leaving the remainder of the telescope time for UV-NIR General Astrophysics. Cost and technical feasibility considerations point to a “sweet spot” in starshade design near a 50-m starshade effective diameter, with 12 or 16 petals, at a distance of 70,000-100,000 km from the telescope.

The Chandra X-ray Survey of Planetary Nebulae (ChanPlaNS): Probing Binarity, Magnetic Fields, and Wind Collisions

We present an overview of the initial results from the Chandra Planetary Nebula Survey (ChanPlaNS), the first systematic (volume-limited) Chandra X-ray Observatory survey of planetary nebulae (PNe) in the solar neighborhood. The first phase of ChanPlaNS targeted 21 mostly high-excitation PNe within ~1.5 kpc of Earth, yielding 3 detections of diffuse X-ray emission and 9 detections of X-ray-luminous point sources at the central stars (CSPNe) of these objects. Combining these results with those obtained from Chandra archival data for all (14) other PNe within ~1.5 kpc that have been observed to date, we find an overall X-ray detection rate of 68%. Roughly 50% of the PNe observed by Chandra harbor X-ray-luminous CSPNe, while soft, diffuse X-ray emission tracing shocks formed by energetic wind collisions is detected in ~30%; five objects display both diffuse and point-like emission components. The presence of X-ray sources appears correlated with PN density structure, in that molecule-poor, elliptical nebulae are more likely to display X-ray emission (either point-like or diffuse) than molecule-rich, bipolar or Ring-like nebulae. All but one of the X-ray point sources detected at CSPNe display X-ray spectra that are harder than expected from hot (~100 kK) central stars emitting as simple blackbodies, possibly indicating a high frequency of binary companions to CSPNe. Other potential explanations include self-shocking winds, non-LTE photospheric emission from the CSPN, or PN mass fallback. Most PNe detected as diffuse X-ray sources are elliptical nebulae that display a nested shell/halo structure and bright ansae; the diffuse X-ray emission regions are confined within inner, sharp-rimmed shells. All sample PNe that display diffuse X-ray emission have inner shell dynamical ages <~5×10^3 yr, placing firm constraints on the timescale for strong shocks due to wind interactions in PNe.

2D Radiative MHD Simulations of the Importance of Partial Ionization in the Chromosphere

The solar chromosphere is weakly ionized and interactions between ionized particles and neutral particles likely have significant consequences for the thermodynamics of the plasma. We investigate the importance of introducing neutral particles using numerical 2.5D radiative MHD simulations obtained with the Bifrost code. The models span from the upper layers of the convection zone to the low corona, and solve the full MHD equations with non-grey and NLTE radiative transfer, and thermal conduction. The effects of partial ionization are implemented using the generalized Ohm’s law. The approximations required in going from three fluids to the generalized Ohm’s law are tested in our simulations. The Ohmic diffusion, the Hall term, and ambipolar diffusion show strong variations in the chromosphere. These strong variations of the various magnetic diffusivities are absent or significantly underestimated when, as has been common for these types of studies, using the VAL-C model as a basis for estimates. In addition, we find that differences in estimating the magnitude of ambipolar diffusion arise depending on which method is used to calculate the ion-neutral collision frequency. These differences cause uncertainties in the different magnetic diffusivity terms. In the chromosphere, we find that the ambipolar diffusion is of the same order of magnitude or even larger than the numerical diffusion used to stabilize our code. As a consequence, ambipolar diffusion produces a strong impact on the modeled atmosphere. Perhaps more importantly, it suggests that at least in the chromospheric domain, self-consistent simulations of the solar atmosphere driven by magneto-convection can accurately describe the impact of the dominant form of resistivity (ambipolar diffusion). This suggests that such simulations may be more realistic in their approach to the lower solar atmosphere than previously assumed.

Enthalpy-based Thermal Evolution of Loops: II. Improvements to the Model

This paper develops the zero-dimensional (0D) hydrodynamic coronal loop model “Enthalpy-based Thermal Evolution of Loops” (EBTEL) proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating, especially impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. The new features in EBTEL are important for accurate tracking of, in particular, the density. The 0D results are compared to a 1D hydro code (Hydrad) with generally good agreement. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function, (b) extensive parameter surveys and (c) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Justification of the two-bulge method in the theory of bodily tides

On various occasions, several authors suggested to model bodily tides with superposition of two symmetrical bulges. One bulge is always aimed at the secondary, and thus implements the instantaneous reaction of the primary’s shape and potential to the tide-rising gravitational pull exerted on it by the secondary. This portion of the tide is called “adiabatic tide” (Zahn 1966a,b) or “elastic tide” (Ferraz Mello 2012; Krasinsky 2006). The second bulge is assumed to align orthogonally to the direction to the tide-raising secondary. So this bulge is set to implement the entire nonelastic portion of the primary’s deformation. This, second bulge is called “dissipative tide” (Zahn 1966a,b; Krasinsky 2006) or “creep tide” (Ferraz Mello 2012). We demonstrate that the two-bulge method is not a separate approximation, but ensues directly from the Fourier expansion of a linear tidal theory equipped with an arbitrary rheological model involving a departure from elasticity. While less economical mathematically, the two-bulge approach has a good illustrative power, and may be employed (like, e.g., in Remus et al. 2012a,b) on a par with a more concise method of complex amplitudes.

 

You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.

Powered by Vote It Up

^ Return to the top of page ^