A prominence eruption on 7 June 2011 produced spectacular curtains of plasma falling through the lower corona. At the solar surface they created an incredible display of extreme ultraviolet brightenings. The aim is to identify and analyze some of the local instabilities which produce structure in the falling plasma. The structures were investigated using SDO/AIA 171A and 193A images in which the falling plasma appeared dark against the bright coronal emission. Several instances of the Rayleigh-Taylor instability were investigated. In two cases the Alfven velocity associated with the dense plasma could be estimated from the separation of the Rayleigh-Taylor fingers. A second type of feature, which has the appearance of self-similar branching horns, is also discussed.

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

We present period measurements of a large sample of field stars in the solar neighbourhood, observed by CoRoT in two different directions of the Galaxy. The presence of a period was detected using the Scargle Lomb Normalized Periodogram technique and the autocorrelation analysis. The assessment of the results has been performed through a consistency verification supported by the folded light curve analysis. The data analysis procedure has discarded a non-negligible fraction of light curves due to instrumental artifacts, however it has allowed us to identify pulsators and binaries among a large number of field stars. We measure a wide range of periods, from 0.25 to 100 days, most of which are rotation periods. The final catalogue includes 1978 periods, with 1727 of them identified as rotational periods, 169 are classified as pulsations and 82 as orbital periods of binary systems. Our sample suffers from selection biases not easily corrected for, thus we do not use the distribution of rotation periods to derive the age distribution of the main-sequence population. Nevertheless, using rotation as a proxy for age, we can identify a sample of young stars (< 600 Myr), that will constitute a valuable sample, supported by further spectroscopic observations, to study the recent star formation history in the solar neighborhood.

Glycolaldehyde is a simple monosaccharide sugar linked to prebiotic chemistry. Recently it was detected in a molecular core in the star-forming region G31.41+0.31 at a reasonably high abundance. We investigate the formation of glycolaldehyde at 10K to determine whether it can form efficiently under typical dense core conditions. Using an astrochemical model, we test five different reaction mechanisms that have been proposed in the astrophysical literature, finding that a gas-phase formation route is unlikely. Of the grain-surface formation routes, only two are efficient enough at very low temperatures to produce sufficient glycolaldehyde to match the observational estimates, with the mechanism culminating in CH3OH + HCO being favoured. However, when we consider the feasibility of these mechanisms from a reaction chemistry perspective, the second grain-surface route looks more promising, H3CO + HCO.

We present some aspects and first results of the emission of sub-eV mass hidden photons from the Sun. The contribution from a resonant region below the photosphere can be quite significant, raising previous estimates. This is relevant for the Telescope for Hidden Photon Search, TSHIPS I, currently targeting at meV-mass hidden photons with O(10^-6) kinetic mixing with the photon. These particles could account for the large effective number of neutrinos pointed at by the cosmic microwave background and other large-scale structure probes, and are motivated in some scenarios of string theory.

CCD observations of V538 Cas have been made on nine nights during three weeks using the AAVSO Bright Star Monitor. No significant variations were found.

Our long-term photometric monitoring of southern nova-like cataclysmic variables with the 1.3-m SMARTS telescope found BB Doradus fading from V ~ 14.3 towards a deep low state at V ~ 19.3 in April 2008. Here we present time-resolved optical spectroscopy of BB Dor in this faint state in 2009. The optical spectrum in quiescence is a composite of a hot white dwarf with Teff = 30000 +- 5000 K and a M3-4 secondary star with narrow emission lines (mainly of the Balmer series and HeI) superposed. We associate these narrow profiles with an origin on the donor star. Analysis of the radial velocity curve of the H-alpha emission from the donor star allowed the measurement of an orbital period of 0.154095 +- 0.000003 d (3.69828 +- 0.00007 h), different from all previous estimates. We detected episodic accretion events which veiled the spectra of both stars and radically changed the line profiles within a timescale of tens of minutes. This shows that accretion is not completely quenched in the low state. During these accretion episodes the line wings are stronger and their radial velocity curve is delayed by ~ 0.2 cycle, similar to that observed in SW Sex and AM Her stars in the high state, with respect to the motion of the white dwarf. Two scenarios are proposed to explain the extra emission: impact of the material on the outer edge of a cold, remnant accretion disc, or the combined action of a moderately magnetic white dwarf (B1 <~ 5 MG) and the magnetic activity of the donor star.

We report two sympathetic solar eruptions, including a partial and a full flux rope eruption in a quadrupolar magnetic region, where a large and a small filament resided above the middle and the east neutral lines respectively. The large filament first rose slowly at a speed of 8 km/s for 23 minutes and then it was accelerated to 102 km/s. Finally, this filament erupted successfully and caused a coronal mass ejection. During the slow rising phase, various evidence for breakout-like external reconnection has been identified at high and low temperature lines. The eruption of the small filament started around the end of the large filament’s slow rising. This filament erupted partially and no associating coronal mass ejection could be detected. Based on a potential field extrapolation, we find that the topology of the three-dimensional coronal field above the source region is composed of three low-lying lobes and a large overlying flux system, and a null point located between the middle lobe and the overlying antiparallel flux system. We propose a possible mechanism within the framework of the magnetic breakout model to interpret the sympathetic filament eruptions, in which the magnetic implosion mechanism is thought to be a possible linkage between the sympathetic eruptions, and the external reconnection at the null point transfers field lines from the middle lobe to the lateral lobes and thereby leads to the full (partial) eruption of the observed large (small) filament. Other possible mechanisms are also discussed briefly. We conclude that the structural properties of coronal fields are important for producing sympathetic eruptions.

The last decade showed an impressive observational effort from the photometric and spectroscopic point of view for ancient stellar clusters in our Galaxy and beyond. The theoretical interpretation of these new observational results requires updated evolutionary models and isochrones spanning a wide range of chemical composition. With this aim we built the new “Pisa Stellar Evolution Database” of stellar models and isochrones by adopting a well-tested evolutionary code (FRANEC) implemented with updated physical and chemical inputs. In particular, our code adopts realistic atmosphere models and an updated equation of state, nuclear reaction rates and opacities calculated with recent solar elements mixture. A total of 32646 models have been computed in the range of initial masses 0.30 – 1.10 Msun for a grid of 216 chemical compositions with the fractional metal abundance in mass, Z, ranging from 0.0001 to 0.01, and the original helium content, Y, from 0.25 to 0.42. Models were computed for both solar-scaled and alpha-enhanced abundances with different external convection efficiencies. Correspondingly, 9720 isochrones were computed in the age range 8 – 15 Gyr, in time steps of 0.5 Gyr. The whole database is available to the scientific community on the web. Models and isochrones were compared with recent calculations available in the literature and with the color-magnitude diagram of selected Galactic globular clusters. The dependence of relevant evolutionary quantities on the chemical composition and convection efficiency were analyzed in a quantitative statistical way and analytical formulations were made available for reader’s convenience.

Recent studies have shown that stellar chromospheric activity, and its effect on convective energy transport in the envelope, is most likely the cause of significant radius and temperature discrepancies between theoretical evolution models and observations. We aim to determine absolute dimensions and abundances for the solar-type detached eclipsing binary EF Aqr, and to perform a detailed comparison with results from recent stellar evolutionary models. uvby-beta standard photometry was obtained with the Stromgren Automatic Telescope. The broadening function formalism was applied on spectra observed with HERMES at the Mercator telescope in La Palma, to obtain radial velocity curves. Masses and radii with a precision of 0.6% and 1.0% respectively have been established for both components of EF Aqr. The active 0.956 M_sol secondary shows star spots and strong Ca II H and K emission lines. The 1.224 M_sol primary shows signs of activity as well, but at a lower level. An [Fe/H] abundance of 0.00+-0.10 is derived with similar abundances for Si, Ca, Sc, Ti, V, Cr, Co, and Ni. Solar calibrated evolutionary models such as Yonsei-Yale, Victoria-Regina and BaSTI isochrones and evolutionary tracks are unable to reproduce EF Aqr, especially for the secondary, which is 9% larger and 400 K cooler than predicted. Models adopting significantly lower mixing length parameters l/H_p remove these discrepancies, as seen in other solar type binaries. For the observed metallicity, Granada models with a mixing length of l/H_p=1.30 (primary) and 1.05 (secondary) reproduce both components at a common age of 1.5+-0.6 Gyr. Observations of EF Aqr suggests that magnetic activity, and its effect on envelope convection, is likely to be the cause of discrepancies in both radius and temperature, which can be removed by adjusting the mixing length parameter of the models downwards.

We present very early UV to optical photometric and spectroscopic observations of the peculiar Type IIn supernova (SN) 2011ht in UGC 5460. The UV observations of the rise to peak are only the second ever recorded for a Type IIn SN and are by far the most complete. The SN, first classified as a SN impostor, slowly rose to a peak of M_V \sim -17 in \sim55 days. In contrast to the \sim2 magnitude increase in the v-band light curve from the first observation until peak, the UV flux increased by >7 magnitudes. The optical spectra are dominated by strong, Balmer emission with narrow peaks (FWHM\sim600 km/s), very broad asymmetric wings (FWHM\sim4200 km/s), and blue shifted absorption (\sim300 km/s) superposed on a strong blue continuum. The UV spectra are dominated by FeII, MgII, SiII, and SiIII absorption lines broadened by \sim1500 km/s. Merged X-ray observations reveal a L_(0.2-10)=(1.0+/-0.2)x10^(39) erg/s. Some properties of SN 2011ht are similar to SN impostors, while others are comparable to Type IIn SNe. Early spectra showed features typical of luminous blue variables at maximum and during giant eruptions. However, the broad emission profiles coupled with the strong UV flux have not been observed in previous SN impostors. The absolute magnitude and energetics (~2.5×10^(49) ergs in the first 112 days) are reminiscent of normal Type IIn SN, but the spectra are of a dense wind. We suggest that the mechanism for creating this unusual profile could be a shock interacting with a shell of material that was ejected a year before the discovery of the SN.

Since the launch of NASA’s Solar Dynamics Observatory on 2010 February 11, the Extreme ultraviolet Variability Experiment (EVE) has observed numerous flares. One interesting feature observed by EVE is that a subset of flares exhibit an additional enhancement of the 2-3 million K emission several hours after the flare’s soft X-ray emission. From the Atmospheric Imaging Assembly (AIA) images, we observe that this secondary emission, dubbed the EUV late phase, occurs in the same active region as the flare but not in the same coronal loops. Here, we examine the C8.8 flare that occurred on 2010 May 5 as a case study of EUV late phase flares. In addition to presenting detailed observations from both AIA and EVE, we develop a physical model of this flare and test it using the Enthalpy Based Thermal Evolution of Loops (EBTEL) model.

In a recent study, Lee et al. presented new photometric follow-up timing observations of the semi-detached binary system SZ Herculis and proposed the existence of two hierarchical cirumbinary companions. Based on the light-travel time effect, the two low-mass M-dwarf companions are found to orbit the binary pair on moderate to high eccentric orbits. The derived periods of these two companions are close to a 2:1 mean-motion orbital resonance. We have studied the stability of the system using the osculating orbital elements as presented by Lee et al. Results indicate an orbit-crossing architecture exhibiting short-term dynamical instabilities leading to the escape of one of the proposed companions. We have examined the system’s underlying model parameter-space by following a Monte Carlo approach and found an improved fit to the timing data. A study of the stability of our best-fitting orbits also indicates that the proposed system is generally unstable. If the observed anomalous timing variations of the binary period is due to additional circumbinary companions, then the resulting system should exhibit a long-term stable orbital configuration much different from the orbits suggested by Lee et al. We, therefore, suggest that based on Newtonian-dynamical considerations, the proposed quadruple system cannot exist. To uncover the true nature of the observed period variations of this system, we recommend future photometric follow-up observations that could further constrain eclipse-timing variations and/or refine light-travel time models.

The BN/KL region of the Orion Nebula is the nearest region of high mass star formation in our galaxy. As such, it has been the subject of intense investigation at a variety of wavelengths, which have revealed it to be brightest in the infrared to sub-mm wavelength regime. Using the newly commissioned SOFIA airborne telescope and its 5-40 micron camera FORCAST, images of the entire BN/KL complex have been acquired. The 31.5 and 37.1 micron images represent the highest resolution observations (31.5 microns, and that this distinction goes instead to the source IRc4. It was determined from these images and derived dust color temperature maps that IRc4 is also likely to be self-luminous. A new source of emission has also been identified at wavelengths >31.5 microns that coincides with the northeastern outflow lobe from the protostellar disk associated with radio source I.

The luminous Young Stellar Object (YSO) IRAS 07422-2001 is studied in the infrared. We discover star forming activity in embedded clusters located in a cloud detected at mid-IR wavelengths in emission. Multiple outflows are discovered from these clusters in the H_2 ro-vibrational line at 2.122 micron. We detect at least six outflows from the cluster associated with the IRAS source and another outflow from a source located in a cluster detected ~2.7 arcmin NE of the IRAS source. Additional star formation is taking place in two other cluster candidates within the cloud. Three of the YSOs in the cluster associated with the IRAS source are detected at 11.2 micron at an angular resolution of ~0.8 arcsec. We have a tentative detection of a circumstellar disk in this cluster, seen as an extinction lane in the J and H-band images. The spectral energy distributions (SEDs) of the dominant YSOs in the cluster associated with the IRAS source and in the NE cluster are studied using radiative transfer models and the properties of the YSOs are estimated. The YSO associated with the IRAS source is probably in a very early Class I stage of formation. The source identified as the dominant YSO in the NE cluster appears to be older than the dominant YSO in the cluster associated with the IRAS source, but its observed flux seems to be contaminated by extra emission, which suggests the presence of a young source contributing to the SED at far-IR wavelengths. The star formation observed in the field of IRAS 07422-2001 supports the idea of hierarchical formation of massive star clusters and the growth of massive young stellar objects near the centres of multiple sub-clusters in a star forming clump through competitive accretion.

The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity fields are accounted for. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3-D nature must be taken into account.

We analyze the role of the symmetry energy slope parameter $L$ on the {\it r}-mode instability of neutron stars. Our study is performed using the microscopic Brueckner–Hartree–Fock approach of the nuclear equation of state, and several phenomenological Skyrme forces and relativistic mean field models. Our results show that the {\it r}-mode instability region is smaller for those models which give larger values of $L$. The reason is that both bulk ($\xi$) and shear ($\eta$) viscosities increase with $L$ and, therefore, the damping of the mode is more efficient for the models with larger $L$. We show also that the dependence of both viscosities on $L$ can be described at each density by simple power-laws of the type $\xi=A_{\xi}L^{B_\xi}$ and $\eta=A_{\eta}L^{B_\eta}$. Using the measured spin frequency and the estimated core temperature of the pulsar in the low-mass X-ray binary 4U 1608-52, we conclude that observational data seems to favor larger values of $L$ if the radius of this object is in the range $11.5-12$(10-12) km and its mass $1.4M_\odot$($2M_\odot$). Outside this range it is not possible to draw any conclusion on $L$ from this pulsar.

We report on results of fully consistent N-body simulations of globular cluster models with N = 100 000 members containing neutron stars and black holes. Using the improved `algorithmic regularization’ method of Hellstrom and Mikkola for compact subsystems, the new code NBODY7 enables for the first time general relativistic coalescence to be achieved for post-Newtonian terms and realistic parameters. Following an early stage of mass segregation, a few black holes form a small dense core which usually leads to the formation of one dominant binary. The subsequent evolution by dynamical shrinkage involves the competing processes of ejection and mergers by radiation energy loss. Unless the binary is ejected, long-lived triple systems often exhibit Kozai cycles with extremely high inner eccentricity (e > 0.999) which may terminate in coalescence at a few Schwarzschild radii. A characteristic feature is that ordinary stars as well as black holes and even BH binaries are ejected with high velocities. On the basis of the models studied so far, the results suggest a limited growth of a few remaining stellar mass black holes in globular clusters.

We report on results of fully consistent N-body simulations of globular cluster models with N = 100 000 members containing neutron stars and black holes. Using the improved `algorithmic regularization’ method of Hellstrom and Mikkola for compact subsystems, the new code NBODY7 enables for the first time general relativistic coalescence to be achieved for post-Newtonian terms and realistic parameters. Following an early stage of mass segregation, a few black holes form a small dense core which usually leads to the formation of one dominant binary. The subsequent evolution by dynamical shrinkage involves the competing processes of ejection and mergers by radiation energy loss. Unless the binary is ejected, long-lived triple systems often exhibit Kozai cycles with extremely high inner eccentricity (e > 0.999) which may terminate in coalescence at a few Schwarzschild radii. A characteristic feature is that ordinary stars as well as black holes and even BH binaries are ejected with high velocities. On the basis of the models studied so far, the results suggest a limited growth of a few remaining stellar mass black holes in globular clusters.

Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipole to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models obtained by direct numerical simulations. The control parameters, the aspect ratio and the mechanical boundary conditions have been varied to build up this sample of models. Both, strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes can in general be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case, the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos co-exist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated ’single mode’. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.

The symmetry energy effects on the location of the inner edge of neutron star crusts are studied. Three phenomenological models are employed in order to check the accuracy of the well known parabolic approximation of the equation of state for asymmetric nuclear matter in the determination of the transition density $n_t$ and transition pressure $P_t$. The results corroborate the statement that the error due to the assumption that a priori the equation of state is parabolic may introduce a large error in the determination of related properties of a neutron star as the crustal fraction of the moment of inertia.

In this study, we present the preliminary light curve analysis of the contact binary SW Lac, using B, V light curves of the system spanning 2 years (2009 – 2010). During the spot modeling process, we used the information coming from the Doppler maps of the system, which was performed using the high resolution and phase dependent spectra obtained at the 2.1 m Otto Struve Telescope of the McDonald Observatory, in 2009. The results showed that the spot modeling from the light curve analysis are in accordance with the Doppler maps, while the non-circular spot modeling technique is needed in order to obtain much better and reliable spot models.

We present a new set of zero metallicity models in the range 13-80 $\rm M_\odot$ together to the associated explosive nucleosynthesis. These models are fully homogeneous with the solar metallicity set we published in Limongi & Chieffi (2006) and will be freely available at the web site http://www.iasf-roma.inaf.it./orfeo/public{\_}html. A comparison between these yields and an average star that represents the average behavior of most of the very metal poor stars in the range $\rm -5.0<[Fe/H]<-2.5$ confirms previous findings that only a fraction of the elemental [X/Fe] may be fitted by the ejecta of $\it standard$ core collapse supernovae.

Vector magnetic fields of an active region filament in the photosphere and upper chromosphere are obtained from spectro-polarimetric observations recorded with the Tenerife Infrared Polarimeter (TIP II) at the German Vacuum Tower Telescope (VTT). We apply Milne-Eddington inversions on full Stokes vectors of the photospheric Si I 1082.7 nm and the upper chromospheric He I triplet at 1083.0 nm to obtain magnetic field vector and velocity maps in two atmosphere layers. We find that: (1)A complete filament was already present in H$\alpha$ at the beginning of the TIP II data acquisition. Only a partially formed one, composed of multiple small threads, was present in He I. (2) The AR filament comprises two sections. One shows strong magnetic field intensities, about 600 – 800 G in the upper chromosphere and 800 – 1000 G in the photosphere. The other exhibits only comparatively weak magnetic field strengths in both layers. (3) The Stokes V signal is indicative of a dip in the magnetic field strength close to the chromospheric PIL. (3) In the chromosphere consistent upflows are found along the PIL flanked by downflows. (4) The transversal magnetic field is nearly parallel to the PIL in the photosphere and inclined by 20 – 30 degree in the chromosphere. (5) The chromospheric magnetic field around the filament is found to be in normal configuration, while the photospheric field presents a concave magnetic topology. The observations are consistent with the emergence of a flux rope with a subsequent formation of a filament.

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R~80,000, sub milli-arcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with three-dimensional, NLTE radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws; Keperian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed SED, polarisation and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

Clumping in stellar winds of hot stars is a possible consequence of radiative-acoustic instability appearing in solutions of radiative-hydrodynamical equations. However, clumping is usually included to stellar atmosphere modeling and radiative transfer calculations in a highly approximate way via a global free parameter called the clumping factor. Using different values of clumping factors many researchers succeeded to fit the observed spectra better and to correct empirical mass loss rates. This usually leads to a conclusion that the stellar wind is clumped. To understand how clumping may influence theoretical predictions of mass-loss rates, different clumping properties have to be taken into account. If clumping appears already below the critical point, the mass-loss rates is changed.

We aim at investigating the morphology, kinematic and helicity evolution of a loop-like prominence during its eruption. We use multi-instrument observations from AIA/SDO, EUVI/STEREO and LASCO/SoHO. The kinematic, morphological, geometrical, and helicity evolution of a loop-like eruptive prominence are studied in the context of the magnetic flux rope model of solar prominences. The prominence eruption evolved as a height expanding twisted loop with both legs anchored in the chromosphere of a plage area. The eruption process consists of a prominence activation, acceleration, and a phase of constant velocity. The prominence body was composed of left-hand (counter-clockwise) twisted threads around the main prominence axis. The twist during the eruption was estimated at 6pi (3 turns). The prominence reached a maximum height of 526 Mm before contracting to its primary location and partially reformed in the same place two days after the eruption. This ejection, however, triggered a CME seen in LASCO C2. The prominence was located in the northern periphery of the CME magnetic field configuration and, therefore, the background magnetic field was asymmetric with respect to the filament position. The physical conditions of the falling plasma blobs were analysed with respect to the prominence kinematics. The same sign of the prominence body twist and writhe, as well as the amount of twisting above the critical value of 2pi after the activation phase indicate that possibly conditions for kink instability were present. No signature of magnetic reconnection was observed anywhere in the prominence body and its surroundings. The filament/prominence descent following the eruption and its partial reformation at the same place two days later suggest a confined type of eruption. The asymmetric background magnetic field possibly played an important role in the failed eruption.

We report the detection of carbon monoxide (CO) emission from the young supernova remnant Cassiopeia A (Cas A) at wavelengths corresponding to the fundamental vibrational mode at 4.65 micron. We obtained AKARI Infrared Camera spectra towards 4 positions which unambiguously reveal the broad characteristic CO ro-vibrational band profile. The observed positions include unshocked ejecta at the center, indicating that CO molecules form in the ejecta at an early phase. We extracted a dozen spectra across Cas A along the long 1 arcmin slits, and compared these to simple CO emission models in Local Thermodynamic Equilibrium to obtain first-order estimates of the excitation temperatures and CO masses involved. Our observations suggest that significant amounts of carbon may have been locked up in CO since the explosion 330 years ago. Surprisingly, CO has not been efficiently destroyed by reactions with ionized He or the energetic electrons created by the decay of the radiative nuclei. Our CO detection thus implies that less carbon is available to form carbonaceous dust in supernovae than is currently thought and that molecular gas could lock up a significant amount of heavy elements in supernova ejecta.

Here we summarize our recent results of high-resolution computer simulations on the turbulent amplification of weak magnetic seed fields showing that such fields will be exponentially amplified also during the gravitational collapse reminiscent to the situation during primordial star formation. The exponential magnetic field amplification is driven by the turbulent small-scale dynamo that can be only observed in computer simulations if the turbulent motions in the central core are sufficiently resolved. We find that the Jeans length, which determines the central core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We conclude from our studies that strong magnetic fields will be unavoidably created already during the formation of the first stars in the Universe, potentially influencing their evolution and mass distribution.

Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, multi-band BVRI polarimetry, broad-band and high-resolution interferometric visibilities and phases. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years. The broad H-band interferometry gives a disk diameter of (1.72+-0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50+-9) deg and a relatively low disk flattening of 1.66+-0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii).

Solar spicules are the fundamental magnetic structures in the chromosphere and considered to play a key role in channelling the chromosphere and corona. Recently, it was suggested by De Pontieu et al. that there were two types of spicules with very different dynamic properties, which were detected by space- time plot technique in the Ca ii H line (3968 A) wavelength from Hinode/SOT observations. ‘Type I’ spicule, with a 3-7 minute lifetime, undergoes a cycle of upward and downward motion; in contrast, ‘Type II’ spicule fades away within dozens of seconds, without descending phase. We are motivated by the fact that for a spicule with complicated 3D motion, the space-time plot, which is made through a slit on a fixed position, could not match the spicule behavior all the time and might lose its real life story. By revisiting the same data sets, we identify and trace 105 and 102 spicules in quiet sun (QS) and coronal hole (CH), respectively, and obtain their statistical dynamic properties. First, we have not found a single convincing example of ‘Type II’ spicules. Secondly, more than 60% of the identified spicules in each region show a complete cycle, i.e., majority spicules are ‘Type I’. Thirdly, the lifetime of spicules in QS and CH are 148 s and 112 s, respectively, but there is no fundamental lifetime difference between the spicules in QS and CH reported earlier. Therefore, the suggestion of coronal heating by ‘Type II’ spicules should be taken with cautions. Subject headings: Sun: chromosphere Sun:transition region Sun:corona

I discuss the role of self-gravity and radiative heating and cooling in shaping the nature of the turbulence in the interstellar medium (ISM) of our galaxy. The heating and cooling cause it to be highly compressible, and, in some regimes of density and temperature, to become thermally unstable, tending to spontaneously segregate into warm/diffuse and cold/dense phases. On the other hand, turbulence is an inherently mixing process, tending to replenish the density and temperature ranges that would be forbidden under thermal processes alone. The turbulence in the ionized ISM appears to be transonic (i.e, with Mach numbers $\Ms \sim 1$), and thus to behave essentially incompressibly. However, in the neutral medium, thermal instability causes the sound speed of the gas to fluctuate by up to factors of $\sim 30$, and thus the flow can be highly supersonic with respect to the dense/cold gas, although numerical simulations suggest that this behavior corresponds more to the ensemble of cold clumps than to the clumps’ internal velocity dispersion. Finally, coherent large-scale compressions in the warm neutral medium (induced by, say, the passage of spiral arms or by supernova shock waves) can produce large, dense molecular clouds that are subject to their own self-gravity, and begin to contract gravitationally. Because they are populated by nonlinear density fluctuations, whose local free-fall times are significantly smaller than that of the whole cloud, the fluctuations terminate their collapse earlier, giving rise to a regime of hierarchical gravitational fragmentation, with small-scale collapses occurring within larger-scale ones. Thus, the “turbulence” in molecular clouds may be dominated by a gravitationally contracting component at all scales.

We study the effects of optically thick clumps, non-void inter-clump medium, variation of the onset of clumping, and velocity dispersion inside clumps on the formation of resonance lines. For this purpose we developed a full 3-D Monte Carlo Radiative Transfer (MCRT) code that is able to handle 3-D shapes of clumps and arbitrary 3-D velocity fields. The method we developed allows us to take into account contributions from density and velocity wind inhomogeneities to the total opacity very precisely. The first comparison with observation shows that 3-D density and velocity wind inhomogeneities have a very strong influence on the resonance line formation, and that they have to be accounted for in order to obtain reliable mass-loss rate determinations.

We report chemical abundances of Na, Al, Zr, La, Nd, and Eu for 39 red giant branch (RGB) stars and 23 potential inner disk red clump stars located in Plaut-s low extinction window. We also measure lithium for a super Li-rich RGB star. The abundances were determined by spectrum synthesis of high resolution (R~25,000), high signal-to-noise (S/N~50-100 pixel-1) spectra obtained with the Blanco 4m telescope and Hydra multifiber spectrograph. For the bulge RGB stars, we find a general increase in the [Na/Fe] and [Na/Al] ratios with increasing metallicity, and a similar decrease in [La/Fe] and [Nd/Fe]. Additionally, the [Al/Fe] and [Eu/Fe] abundance trends almost identically follow those of the {\alpha}-elements, and the [Zr/Fe] ratios exhibit relatively little change with [Fe/H]. The consistently low [La/Eu] ratios of the RGB stars indicate that at least a majority of bulge stars formed rapidly (<1 Gyr) and before the main s-process could become a significant pollution source. In contrast, we find that the potential inner disk clump stars exhibit abundance patterns more similar to those of the thin and thick disks. Comparisons between the abundance trends at different bulge locations suggest that the inner and outer bulge formed on similar timescales. However, we find evidence of some abundance differences between the most metal-poor and metal-rich stars in various bulge fields. The data also indicate that the halo may have had a more significant impact on the outer bulge initial composition than the inner bulge composition. The [Na/Fe] and to a lesser extent [La/Fe] abundances further indicate that the metal-poor bulge, at least at ~1 kpc from the Galactic center, and thick disk may not share an identical chemistry.

HD 179949 is an F8V star, orbited by a giant planet at ~8 R* every 3.092514 days. The system was reported to undergo episodes of stellar activity enhancement modulated by the orbital period, interpreted as caused by Star-Planet Interactions (SPIs). One possible cause of SPIs is the large-scale magnetic field of the host star in which the close-in giant planet orbits. In this paper we present spectropolarimetric observations of HD 179949 during two observing campaigns (2009 September and 2007 June). We detect a weak large-scale magnetic field of a few Gauss at the surface of the star. The field configuration is mainly poloidal at both observing epochs. The star is found to rotate differentially, with a surface rotation shear of dOmega=0.216\pm0.061 rad/d, corresponding to equatorial and polar rotation periods of 7.62\pm0.07 and 10.3\pm0.8 d respectively. The coronal field estimated by extrapolating the surface maps resembles a dipole tilted at ~70 degrees. We also find that the chromospheric activity of HD 179949 is mainly modulated by the rotation of the star, with two clear maxima per rotation period as expected from a highly tilted magnetosphere. In September 2009, we find that the activity of HD 179949 shows hints of low amplitude fluctuations with a period close to the beat period of the system.

We analyse medium-resolution spectra (R\sim 18000) of 19 late type dwarfs in order to determine vsini-s using synthetic rather than observational template spectra. For this purpose observational data around 2.2 $\mu$m of stars with spectral classes from G8V to M9.5V were modelled. We find that the Na I (2.2062 and 2.2090 $\mu$m) and $^{12}$CO 2-0 band features are modelled well enough to use for vsini determination without the need for a suitable observational template spectra. Within the limit of the resolution of our spectra, we use synthetic spectra templates to derive vsini values consistent with those derived in the optical regime using observed templates. We quantify the errors in our vsini determination due to incorrect choice of model parameters \Teff, log $g$, $v_{\rm micro}$, [Fe/H] or FWHM and show that they are typically less than 10 per cent. We note that the spectral resolution of our data(\sim 16 km/s) limited this study to relatively fast rotators and that resolutions of 60000 will required to access most late-type dwarfs.

We present models of the components of the systems KOI-126 and CM Draconis, the two eclipsing binary systems known to date to contain stars with masses low enough to have fully convective interiors. We are able to model satisfactorily the system KOI-126, finding consistent solutions for the radii and surface temperatures of all three components, using a solar-like value of the mixing-length parameter \alpha in the convection zone, and PHOENIX NextGen 1D model atmospheres for the surface boundary conditions. Depending on the chemical composition, we estimate the age of the system to be in the range 3-5 Gyr. For CM Draconis, on the other hand, we cannot reconcile our models with the observed radii and T_eff using the current metal-poor composition estimate based on kinematics. Higher metallicities lessen but do not remove the discrepancy. We then explore the effect of varying the mixing length parameter \alpha. As previously noted in the literature, a reduced \alpha can be used as a simple measure of the lower convective efficiency due to rotation and induced magnetic fields. Our models show a sensitivity to \alpha (for \alpha < 1.0) sufficient to partially account for the radius discrepancies. It is, however, impossible to reconcile the models with the observations on the basis of the effect of the reduced \alpha alone. We therefore suggest that the combined effects of high metallicity and \alpha reduction could explain the observations of CM Draconis. For example, increasing the metallicity of the system towards super-solar values (i.e. Z = 2 Z_sun) yields an agreement within 2 \sigma with \alpha = 1.0.

The planetary nebula (PN) A 30 is believed to have undergone a very late thermal pulse resulting in the ejection of knots of hydrogen-poor material. Using multi-epoch HST images we have detected the angular expansion of these knots and derived an age of 850+280-150 yr. To investigate the spectral and spatial properties of the soft X-ray emission detected by ROSAT, we have obtained Chandra and XMM-Newton deep observations of A 30. The X-ray emission from A 30 can be separated into two components: a point-source at the central star and diffuse X-ray emission associated with the hydrogen-poor knots and the cloverleaf structure inside the nebular shell. To help us assess the role of the current stellar wind in powering this X-ray emission, we have determined the stellar parameters and wind properties of the central star of A 30 using a non-LTE model fit to its optical and UV spectrum. The spatial distribution and spectral properties of the diffuse X-ray emission is highly suggestive that it is generated by the post-born-again and present fast stellar winds interacting with the hydrogen-poor ejecta of the born-again event. Charge-exchange reactions between the ions of the stellar winds and neutral material of the born-again ejecta seem the most likely mechanism for the production of diffuse X-ray emission. Shock-heated plasma may also contribute to this emission, as the hydrogen-poor knots are ablated by the stellar winds, in which case efficient mass-loading of the stellar winds is needed to raise the density and damp the velocity of the stellar winds to reproduce the spectral properties of the diffuse emission. The origin of the X-ray emission from the central star of A 30 is puzzling: shocks in the present fast stellar wind and photospheric emission can be ruled out, while the development of a new, compact hot bubble confining the fast stellar wind seems implausible.

The evolutionary properties of the old metal-rich Galactic open cluster NGC6791 are assessed, based on deep UB photometry and 2Mass JK data. For 4739 stars in the cluster, bolometric luminosity and effective temperature have been derived from theoretical (U-B) and (J-K) color fitting. The derived H-R diagram has been matched with the UVBLUE grid of synthetic stellar spectra to obtain the integrated SED of the system, together with a full set UV (Fanelli) and optical (Lick) narrow-band indices. The cluster appears to be a fairly good proxy of standard elliptical galaxies, although with significantly bluer infrared colors, a shallower 4000A Balmer break, and a lower Mg2 index. The confirmed presence of a dozen hot stars, along their EHB evolution, leads the cluster SED to consistently match the properties of the most active UV-upturn galaxies, with 1.7+/-0.4% of the total bolometric luminosity emitted shortward of 2500A. The cluster Helium abundance results Y=0.30 +/-0.04, while the Post-MS implied stellar lifetime from star number counts fairly agrees with the theoretical expectations from both the Padova and BASTI stellar tracks. A Post-MS fuel consumption of 0.43 +/- 0.01 M_sun is found for NGC6791 stars, in close agreement with the estimated mass of cluster He-rich white dwarfs. Such a tight figure may lead to suspect that a fraction of the cluster stellar population does actually not reach the minimum mass required to effectively ignite He in the stellar core.

We examine eight young stellar objects in the OMC-2 star forming region based on observations from the SOFIA/FORCAST early science phase, the Spitzer Space Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the literature. We show the spectral energy distributions of these objects from near-infrared to millimeter wavelengths, and compare the SEDs with those of sheet collapse models of protostars and circumstellar disks. Four of the objects can be modelled as protostars with infalling envelopes, two as young stars surrounded by disks, and the remaining two objects have double-peaked SEDs. We model the double-peaked sources as binaries containing a young star with a disk and a protostar. The six most luminous sources are found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities ranging from 300 L_sun to 20 L_sun. The most embedded source (OMC-2 FIR 4) can be fit by a class 0 protostar model having a luminosity of ~50 L_sun and mass infall rate of ~10^-4 solar masses per year.

Luminous blue variables (LBVs) are expected to play an important role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. In this paper we provide a multiwavelength study of the population of (candidate) LBVs identified within M33. New spectra provide an observational baseline of >4yr with respect to published data, which is well suited to identifying LBV outbursts. Multi-epoch optical and mid-IR surveys of M33 further constrain the variability of the sample and permit a search for dusty circumstellar ejecta. Spectroscopic and photometric variability appears common amongst the sample, although in many cases further observations will be required to determine its physical origin. Nevertheless, we report a new outburst of M33 Var C, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993-2010 confirms an LBV classification. Proof-of-concept quantitative analysis is provided for Romano’s star; the results being consistent with the finding that its bolometric luminosity varies during its LBV excursions. The combination of the temperature and luminosity of two stars, the B hypergiant [HS80] 110A and the cool hypergiant B324, appears to be in violation of the empirical Humphreys-Davidson limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects as extreme as Eta Carinae are identified. The combined multiwavelength dataset suggests that the population of LBVs studied is contaminated by stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity stars associated with hot dust are of particular interest since they appear similar to the likely progenitors of SN 2008S and the 2008 NGC300 transient, albeit suffering less intrinsic extinction. [ABRIDGED]

The analysis of recent, extended multicolour CCD and archive photoelectric, photographic and visual observations has revealed several important properties of RZ Lyr, an RRab-type variable exhibiting large-amplitude Blazhko modulation. On the time-base of \sim110 yr, a strict anticorrelation between the pulsation and modulation period changes is established. The light curve of RZ Lyr shows a remarkable bump on the descending branch in the small-amplitude phase of the modulation, similarly to the light curves of bump Cepheids. We speculate that the stellar structure temporally suits a 4:1 resonance between the periods of the fundamental and one of the higher-order radial modes in this modulation phase. The light-curve variation of RZ Lyr can be correctly fitted with a two-modulation-component solution; the 121 d period of the main modulation is nearly but not exactly four times longer than the period of the secondary modulation component. Using the inverse photometric method, the variations in the pulsation-averaged values of the physical parameters in different phases of both modulation components are determined.

New instrumental capabilities and the wealth of astrophysical information extractable from the near-infrared wavelength region have led to a growing interest in the field of high resolution spectroscopy at 1-5 mu. We aim to provide a library of observed high-resolution and high signal-to-noise-ratio near-infrared spectra of stars of various types throughout the Hertzsprung-Russell diagram. This is needed for the exploration of spectral features in this wavelength range and for comparison of reference targets with observations and models. High quality spectra were obtained using the CRIRES near-infrared spectrograph at ESO’s VLT covering the range from 0.97 to 5.3 mu at high spectral resolution. Accurate wavelength calibration and correction for of telluric lines were performed by fitting synthetic transmission spectra for the Earth’s atmosphere to each spectrum individually. We describe the observational strategy and the current status and content of the library which includes 13 objects. The first examples of finally reduced spectra are presented. This publication will serve as a reference paper to introduce the library to the community and explore the extensive amount of material.

Globally-propagating coronal bright fronts (CBFs) in the solar corona are among the most dramatic manifestations of solar activity, but are not well understood despite strong links with both solar flares and coronal mass ejections. Extreme UltraViolet (EUV) observations from the STEREO and SDO spacecraft are used here to study their kinematics and morphology. The first STEREO observations of a CBF are presented, with the pulse observed in all available EUV passbands (171, 195, 284 and for the first time, 304A). The pulse displayed similar kinematics in all passbands, although the derived pulse velocity and acceleration were found to be strongly influenced by the observing cadence, implying that previous kinematics may have been underestimated. Different techniques for identifying CBFs and deriving their true kinematics were tested, with traditional techniques shown to be prone to undefined user-dependent errors. This was overcome through the development of a statistically rigorous, semi-automated identification algorithm, which was then used to determine the kinematics of four CBF events observed by STEREO. All of the events studied exhibited clear deceleration as well as increases in both spatial and temporal pulse width, indicating that the CBFs are dispersive. A CBF pulse observed by both STEREO and SDO was also studied using the algorithm, with lower initial velocity and weaker deceleration noted in STEREO observations compared to SDO, reaffirming the effects of image cadence on the derived kinematics. The kinematics obtained using SDO were highly passband dependent, suggesting a compressive nature. Significant pulse broadening was also noted in observations from both spacecraft, allowing the dispersion rate of the pulse to be determined. These results indicate that coronal bright fronts are best interpreted as fast-mode magnetoacoustic waves propagating in an inhomogeneous medium.

In this contribution we describe some recent observations of high-speed magnetized flows in the quiet Sun granulation. These observations were carried out with the Imaging Magnetograph eXperiment (IMaX) onboard the stratospheric balloon {\sc Sunrise}, and possess an unprecedented spatial resolution and temporal cadence. These flows were identified as highly shifted circular polarization (Stokes $V$) signals. We estimate the LOS velocity responsible for these shifts to be larger than 6 km s$^{-1}$, and therefore we refer to them as {\it supersonic magnetic flows}. The average lifetime of the detected events is 81.3 s and they occupy an average area of about 23\,000 km$^2$. Most of the events occur within granular cells and correspond therefore to upflows. However some others occur in intergranular lanes or bear no clear relation to the convective velocity pattern. We analyze a number of representative examples and discuss them in terms of magnetic loops, reconnection events, and convective collapse.

Based on the turbulent convection model (TCM), we investigate chemical mixing in the bottom overshooting region of the convective envelope of intermediate-mass stars, focusing on its influence on the formation and extension of blue loops in the Hertzsprung-Russell (HR) diagram. A diffusive mixing model is adopted during the Red Giant Branch (RGB) phase. The properties of the blue loop are changed by modification of the element profiles above the H-burning shell, which results from the incomplete mixing in the bottom overshooting region when stellar model evolves up along the RGB. Such modification of the element profiles will lead to an increase of the opacity in the region just above the H-burning shell and a decrease of the opacity in the outer homogeneous convection zone, which will result in a quick decrease of the H-shell nuclear luminosity $L_{H}$ when the stellar model evolves from the RGB tip to its bottom and, finally, a much weaker and smaller convection zone will be obtained in the stellar envelope. This helps to form a longer blue loop. The extension of the blue loop is very sensitive to the parameters $(C_{X}$ and $\alpha_{TCM})$ of the diffusive mixing model and of the TCM. The results mainly show that: 1) comparing to the results of the classical model with the mixing-length theory, the lengths of the obtained blue loops with different combinations of the values of $C_{X}$ and $\alpha_{TCM}$ are all increased and the length of the blue loop increases with the values of parameters $C_{X}$ and $\alpha_{TCM}$; 2) the diffusive mixing model can significantly extend the time of stellar models lingering on the blue side of the HR diagram, even though the length of the blue loop for the 7$M_{\odot}$ star has a less prominent difference between the classical and diffusive mixing model;

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6″, corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation.

2011 was a successful year for the GREGOR project. The telescope was finally completed in May with the installation of the 1.5-meter primary mirror. The installation of the first-light focal plane instruments was completed by the end of the year. At the same time, the preparations for the installation of the high-order adaptive optics were finished, its integration to the telescope is scheduled for early 2012. This paper describes the telescope and its instrumentation in their present first light configuration, and provides a brief overview of the science goals of GREGOR.

The Kelvin-Helmholtz Instability (KHI) has been observed in the solar atmosphere. Ion-neutral collisions may play a relevant role for the growth rate and evolution of the KHI in solar partially ionized plasmas as in, e.g., solar prominences. Here, we investigate the linear phase of the KHI at an interface between two partially ionized magnetized plasmas in the presence of a shear flow. The effects of ion-neutral collisions and compressibility are included in the analysis. We obtain the dispersion relation of the linear modes and perform parametric studies of the unstable solutions. We find that in the incompressible case the KHI is present for any velocity shear regardless the value of the collision frequency. In the compressible case, the domain of instability depends strongly on the plasma parameters, specially the collision frequency and the density contrast. For high collision frequencies and low density contrasts the KHI is present for super-Alfvenic velocity shear only. For high density contrasts the threshold velocity shear can be reduced to sub-Alfvenic values. For the particular case of turbulent plumes in prominences, we conclude that sub-Alfvenic flow velocities can trigger the KHI thanks to the ion-neutral coupling.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

Due to their transitionary nature, yellow supergiants provide a critical challenge for evolutionary modeling. Previous studies within M31 and the SMC show that the Geneva evolutionary models do a poor job at predicting the lifetimes of these short-lived stars. Here we extend this study to the LMC while also investigating the galaxy’s red supergiant content. This task is complicated by contamination by Galactic foreground stars that color and magnitude criteria alone cannot weed out. Therefore, we use proper motions and the LMC’s large systemic radial velocity (\sim278 km/s) to separate out these foreground dwarfs. After observing nearly 2,000 stars, we identified 317 probable yellow supergiants, 6 possible yellow supergiants and 505 probable red supergiants. Foreground contamination of our yellow supergiant sample was \sim80%, while that of the the red supergiant sample was only 3%. By placing the yellow supergiants on the H-R diagram and comparing them against the evolutionary tracks, we find that new Geneva evolutionary models do an exemplary job at predicting both the locations and the lifetimes of these transitory objects.

We present a new numerical scheme for solving the advection equation and its application to the Vlasov simulation. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions compared to other existing schemes with the same memory usage. The two-dimensional scheme can solve the solid body rotation problem of a gaussian profile with little numerical diffusion. This is a very important property for Vlasov simulations of magnetized plasma. The application of the scheme to the electromagnetic Vlasov simulation of collisionless shock waves is presented as a benchmark test.

Turbulent convection models are thought to be good tools to deal with the convective overshooting in the stellar interior. However, they are too complex to be applied in calculations of stellar structure and evolution. In order to understand the physical processes of the convective overshooting and to simplify the application of turbulent convection models, a semi-analytic solution is necessary. We obtain the approximate solution and asymptotic solution of the turbulent convection model in the overshooting region, and find some important properties of the convective overshooting: I. The overshooting region can be partitioned into three parts: a thin region just outside the convective boundary with high efficiency of turbulent heat transfer, a power law dissipation region of turbulent kinetic energy in the middle, and a thermal dissipation area with rapidly decreasing turbulent kinetic energy. The decaying indices of the turbulent correlations $k$, $\bar{u_{r}’T'}$, and $\bar{T’T'}$ are only determined by the parameters of the TCM, and there is an equilibrium value of the anisotropic degree $\omega$. II. The overshooting length of the turbulent heat flux $\bar{u_{r}’T'}$ is about $1H_k$($H_k=|\frac{dr}{dlnk}|$). III. The value of the turbulent kinetic energy at the convective boundary $k_C$ can be estimated by a method called \textsl{the maximum of diffusion}. Turbulent correlations in the overshooting region can be estimated by using $k_C$ and exponentially decreasing functions with the decaying indices.

Turbulent convection models are thought to be good tools to deal with the convective overshooting in the stellar interior. However, they are too complex to be applied in calculations of stellar structure and evolution. In order to understand the physical processes of the convective overshooting and to simplify the application of turbulent convection models, a semi-analytic solution is necessary. We obtain the approximate solution and asymptotic solution of the turbulent convection model in the overshooting region, and find some important properties of the convective overshooting: I. The overshooting region can be partitioned into three parts: a thin region just outside the convective boundary with high efficiency of turbulent heat transfer, a power law dissipation region of turbulent kinetic energy in the middle, and a thermal dissipation area with rapidly decreasing turbulent kinetic energy. The decaying indices of the turbulent correlations $k$, $\bar{u_{r}’T'}$, and $\bar{T’T'}$ are only determined by the parameters of the TCM, and there is an equilibrium value of the anisotropic degree $\omega$. II. The overshooting length of the turbulent heat flux $\bar{u_{r}’T'}$ is about $1H_k$($H_k=|\frac{dr}{dlnk}|$). III. The value of the turbulent kinetic energy at the convective boundary $k_C$ can be estimated by a method called \textsl{the maximum of diffusion}. Turbulent correlations in the overshooting region can be estimated by using $k_C$ and exponentially decreasing functions with the decaying indices.

Most of the observational studies of supernova (SN) explosions are limited to early phases (<a> 100 yr) in our Galaxy or very nearby galaxies. SNe at the epoch between these two, which we call “transitional” phase, have not been explored in detail except for several extragalactic SNe including SN 1987A in the Large Magellanic Cloud. We present theoretical predictions for the infrared (IR) dust emissions by several mechanisms; emission from dust formed in the SN ejecta, light echo by circumstellar and interstellar dust, and emission from shocked circumstellar dust. We search for IR emission from 6 core-collapse SNe at the transitional phase in the nearby galaxies NGC 1313, NGC 6946, and M101 by using the data taken with the AKARI satellite and Spitzer. Among 6 targets, we detect the emission from SN 1978K in NGC 1313. SN 1978K is associated with 1.3 x 10^{-3} Msun of silicate dust. We show that, among several mechanisms, the shocked circumstellar dust is the most probable emission source to explain the IR emission observed for CSM-rich SN 1978K. IR emission from the other 5 objects is not detected. Our current observations are sensitive to IR luminosity of > 10^{38} erg s^{-1}, and the non-detection of SN 1962M excludes the existence of the shocked circumstellar dust for a high gas mass-loss rate of sim 10^{-4} Msun yr^{-1}. Observations of SNe at the transitional phase with future IR satellites will fill the gap of IR observations of SNe with the age of 10-100 years, and give a new opportunity to study the circumstellar and interstellar environments of the progenitor, and possibly dust formation in SNe.

We identify two new tidally distorted white dwarfs (WDs), SDSS J174140.49+652638.7 and J211921.96-001825.8 (hereafter J1741 and J2119). Both stars are extremely low mass (ELM, < 0.2 Msun) WDs in short-period, detached binary systems. High-speed photometric observations obtained at the McDonald Observatory reveal ellipsoidal variations and Doppler beaming in both systems; J1741, with a minimum companion mass of 1.1 Msun, has one of the strongest Doppler beaming signals ever observed in a binary system (0.59 \pm 0.06% amplitude). We use the observed ellipsoidal variations to constrain the radius of each WD. For J1741, the star's radius must exceed 0.074 Rsun. For J2119, the radius exceeds 0.10 Rsun. These indirect radius measurements are comparable to the radius measurements for the bloated WD companions to A-stars found by the Kepler spacecraft, and they constitute some of the largest radii inferred for any WD. Surprisingly, J1741 also appears to show a 0.23 \pm 0.06% reflection effect, and we discuss possible sources for this excess heating. Both J1741 and J2119 are strong gravitational wave sources, and the time-of-minimum of the ellipsoidal variations can be used to detect the orbital period decay. This may be possible on a timescale of a decade or less.

The active region NOAA 11158 produced the first X-class flare of Solar Cycle 24, an X2.2 flare at 01:44 UT on 2011 February 15. Here we analyze SDO/HMI magnetograms covering a 12-hour interval centered at the time of this flare. We describe the spatial distributions of the photospheric magnetic changes associated with this flare, including the abrupt changes in the field vector, vertical electric current and Lorentz force vector. We also trace these parameters’ temporal evolution. The abrupt magnetic changes were concentrated near the neutral line and in two neighboring sunspots. Near the neutral line, the field vectors became more horizontal during the flare and the shear increased. This was due to an increase in strength of the horizontal field components near the neutral line, most significant in the horizontal component parallel to the neutral line but the perpendicular component also increased in strength. The vertical component did not show a significant, permanent overall change at the neutral line. The increase in total flux at the neutral line was accompanied by a compensating flux decrease in the surrounding volume. In both of the sunspots near the neutral line the azimuthal flux abruptly decreased during the flare but this change was permanent in only one of the spots. There was a large, abrupt, downward vertical Lorentz force change during the flare, consistent with results of past analyses and recent theoretical work. The horizontal Lorentz force acted in opposite directions on each side of neutral line, with the two sunspots at each end subject to abrupt torsional forces. The shearing forces were consistent with a decrease of shear near the neutral line, whereas the field itself became more sheared as a result of the flux collapsing towards the neutral line from the surrounding volume.

We investigate the use of the gravity sensitive neutral sodium (NaI) doublet at 8183 Angstroms 8195 Angstroms (Na 8200 Angstrom doublet) as an age indicator for M dwarfs. We measured the Na doublet equivalent width (EW) in giants, old dwarfs, young dwarfs, and candidate members of the Beta Pic moving group using medium resolution spectra. Our Na 8200 Angstrom doublet EW analysis shows that the feature is useful as an approximate age indicator in M-type dwarfs with (V-K_s) >= 5.0, reliably distinguishing stars older and younger than 100 Myr. A simple derivation of the dependence of the Na EW on temperature and gravity supports the observational results. An analysis of the effects of metallicity show that this youth indicator is best used on samples with similar metallicity. The age estimation technique presented here becomes useful in a mass regime where traditional youth indicators are increasingly less reliable, is applicable to other alkali lines, and will help identify new-low mass members in other young clusters and associations.

Using dimensional analysis techniques we present an extension of Newton’s gravitational theory built under the assumption that Milgrom’s acceleration constant is a fundamental quantity of nature. The gravitational force converges to Newton’s gravity and to a MOND-like description in two different mass and length regimes. It is shown that a modification on the force sector (and not in the dynamical one as MOND does) is more convenient and can reproduce and predict different phenomena usually ascribed to dark matter at the non-relativistic level.

We present 37\micron\ imaging of the S140 complex of infrared sources centered on IRS1 made with the FORCAST camera on SOFIA. These observations are the longest wavelength imaging to resolve clearly the three main sources seen at shorter wavelengths, IRS 1, 2 and 3, and are nearly at the diffraction limit of the 2.5-m telescope. We also obtained a small number of images at 11 and 31\micron\ that are useful for flux measurement. Our images cover the area of several strong sub-mm sources seen in the area — SMM 1, 2, and 3 — that are not coincident with any mid-infrared sources and are not visible in our longer wavelength imaging either. Our new observations confirm previous estimates of the relative dust optical depth and source luminosity for the components in this likely cluster of early B stars. We also investigate the use of super-resolution to go beyond the basic diffraction limit in imaging on SOFIA and find that the van Cittert algorithm, together with the “multi-resolution” technique, provides excellent results.

We present the results of a comparison between classical and newly identified EXor based on literature data and aimed at recognizing possible differences or similarities of both categories. Optical and near-IR two-color diagrams, modalities of fluctuations, and derived values of the mass accretion rates are indicative of strong similarities between the two samples. We demonstrate how the difference between the outburst and the quiescence spectral energy distribution of all the EXor can be well fitted with a single blackbody, as if an additional thermal component appears during the outbursting phase. Temperatures of this additional component span between 1000 and 4500 K, while the radii of the emitting regions (assumed to be a uniform disk) span between 0.01 and 0.1 AU, sizes typical of the inner portions of the circumstellar disk. Spots persisting up to 50% of the outburst duration, not exceeding the 10% of the stellar surface, and with temperatures compatible with the EXor mass accretion rates, are able to account for both the appearance of the additional thermal component and the dust sublimation in the inner structures of the disk. We also compare the EXor events with the most significant color and magnitude fluctuations of active T Tauri stars finding that (i} burst accretion phenomena should also be important for this latter class; (ii} EXor events could be more frequent then those accidentally discovered. Remarkable is the case of the source V2493 Cyg, a T Tauri star recently identified as a strong outbursting object: new optical and near-IR photometric and spectroscopic data are presented trying to clarify its EXor or FUor nature.

We review (i) observations and numerical simulations of vortical flows in the solar atmosphere and (ii) measurements of the horizontal magnetic field in quiet Sun regions. First, we discuss various manifestations of vortical flows and emphasize the role of magnetic fields in mediating swirling motion created near the solar surface to the higher layers of the photosphere and to the chromosphere. We reexamine existing simulation runs of solar surface magnetoconvection with regard to vortical flows and compare to previously obtained results. Second, we review contradictory results and problems associated with measuring the angular distribution of the magnetic field in quiet Sun regions. Furthermore, we review the Stokes-V-amplitude ratio method for the lines Fe I 630.15 and 630.25 nm. We come to the conclusion that the recently discovered two distinct populations in scatter plots of this ratio must not bee interpreted in terms of ‘uncollapsed’ and ‘collapsed’ fields but they stem from weak granular magnetic fields and weak canopy fields located at the boundaries between granules and the intergranular space. Based on new simulation runs, we reaffirm earlier findings of a predominance of the horizontal field components over the vertical one, particularly in the upper photosphere and at the base of the chromosphere.

Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy stored in the magnetic field by permitting a change in the magnetic topology. In this article we present a review of the current state of understanding of magnetic reconnection. We discuss theoretical results regarding the formation of current sheets in complex 3D magnetic fields, and describe the fundamental differences between reconnection in two and three dimensions. We go on to outline recent developments in modelling of reconnection with kinetic theory, as well as in the MHD framework where a number of new 3D reconnection regimes have been identified. We discuss evidence from observations and simulations of solar system plasmas that support this theory, and summarise some prominent locations in which this new reconnection theory is relevant in astrophysical plasmas.

Solar coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic field from the Sun into the corona and interplanetary space. They are the most significant drivers of adverse space weather at Earth and other locations in the heliosphere, so it is important to understand the physics governing their eruption and propagation. However the diffuse morphology and transient nature of CMEs makes them difficult to identify and track using traditional image processing techniques. In this thesis the implementation of multiscale image processing techniques to identify and track the CME front through coronagraph images is detailed. An ellipse characterisation of the CME front is used to determine the CME kinematics and morphology with increased precision as compared to techniques used in current CME catalogues, and efforts are underway to automate this procedure for applying to a large number of CME observations for future analysis. It was found that CMEs do not simply undergo constant acceleration, but rather tend to show a higher acceleration early in their propagation. The angular width of CMEs was also found to change as they propagate, normally increasing with height from the Sun. However these results were derived from plane-of-sky measurements with no correction for how the true CME geometry and direction affect the kinematics and morphology observed. With the advent of the unique dual perspectives of the STEREO spacecraft, the multiscale methods were extended to an elliptical tie-pointing technique in order reconstruct the front of a CME in three-dimensions. Applying this technique to the Earth-directed CME of 12 December 2008 allowed an accurate determination of its true kinematics and morphology, and the CME was found to undergo early acceleration, non-radial motion, angular width expansion, and aerodynamic drag in the solar wind as it propagated towards Earth.

A well known behavior of EUV light curves of discrete coronal loops is that the peak intensities of cooler channels or spectral lines are reached at progressively later times than hotter channels. This time lag is understood to be the result of hot coronal loop plasma cooling through these lower respective temperatures. However, loops typically comprise only a minority of the total emission in active regions. Is this cooling pattern a common property of active region coronal plasma, or does it only occur in unique circumstances, locations, and times? The new SDO/AIA data provide a wonderful opportunity to answer this question systematically for an entire active region. We measure the time lag between pairs of SDO/AIA EUV channels using 24 hours of images of AR 11082 observed on 19 June 2010. We find that there is a time-lag signal consistent with cooling plasma, just as is usually found for loops, throughout the active region including the diffuse emission between loops for the entire 24 hour duration. The pattern persists consistently for all channel pairs and choice of window length within the 24 hour time period, giving us confidence that the plasma is cooling from temperatures of greater than 3 MK, and sometimes exceeding 7 MK, down to temperatures lower than ~ 0.8 MK. This suggests that the bulk of the emitting coronal plasma in this active region is not steady; rather, it is dynamic and constantly evolving. These measurements provide crucial constraints on any model which seeks to describe coronal heating.

We describe the detonation mechanism comprising the “Pulsationally Assisted” Gravitationally Confined Detonation (GCD) model of Type Ia supernovae (SNe Ia). This model is analogous to the previous GCD model reported in Jordan (2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final energy releases and nickel-56 yields conform better to observational values than is the case for the “classical” GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel) the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that produce a detonation. We demonstrate this mechanism with three 3-dimensional (3D), full star simulations of this model using the FLASH code, varying the initial offset of the ignition points for each model. The simulations are characterized by nuclear energy releases ranging from 38% to 78% of the binding energy of the white dwarf during the deflagration phase. We show that the conditions for detonation are achieved in all three of the models.

Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. The long-term datasets revealed that they are remarkably stable over long periods (>40yr), with the possible exception of zeta Sco prior to the 20th century, in contrast to the typical excursions that characterise luminous blue variables (LBVs). Zeta Sco, HD190603 and BP Cru possess physical properties intermediate between B supergiants and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars (for initial masses in the range ~30-60Msun). In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110Msun) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, the properties of Cyg OB2 #12 are difficult to understand under current evolutionary paradigms. [ABRIDGED]

A Bayesian approach to calibrating period-luminosity (PL) relations has substantial benefits over generic least-squares fits. In particular, the Bayesian approach takes into account the full prior distribution of the model parameters, such as the a priori distances, and refits these parameters as part of the process of settling on the most highly-constrained final fit. Additionally, the Bayesian approach can naturally ingest data from multiple wavebands and simultaneously fit the parameters of PL relations for each waveband in a procedure that constrains the parameter posterior distributions so as to minimize the scatter of the final fits appropriately in all wavebands. Here we describe the generalized approach to Bayesian model fitting and then specialize to a detailed description of applying Bayesian linear model fitting to the mid-infrared PL relations of RR Lyrae variable stars. For this example application we quantify the improvement afforded by using a Bayesian model fit. We also compare distances previously predicted in our example application to recently published parallax distances measured with the Hubble Space Telescope and find their agreement to be a vindication of our methodology. Our intent with this article is to spread awareness of the benefits and applicability of this Bayesian approach and encourage future PL relation investigations to consider employing this powerful analysis method.

We present the results of our photometric and spectroscopic studies of the new eclipsing cataclysmic variable star 1RXS J180834.7+101041. Its spectrum exhibits double-peaked hydrogen and helium emission lines. The Doppler maps constructed from hydrogen lines show a nonuniform distribution of emission in the disk similar to that observed in IP Peg. This suggests that the object can be a cataclysmic variable with tidal density waves in the disk. We have determined the component masses (M_WD =0.8 \pm 0.22 M_sun and M_RD =0.14 \pm 0.02 M_sun) and the binary inclination (i =78 \pm 1.5 deg) based on well-known relations between parameters for cataclysmic variable stars. We have modeled the binary light curves and showed that the model of a disk with two spots is capable of explaining the main observed features of the light curves.

This paper presents a detailed assessment of the ability of the 240 Skyrme interaction parameter sets in the literature to satisfy a series of criteria derived from macroscopic properties of nuclear matter in the vicinity of nuclear saturation density at zero temperature and their density dependence, derived by the liquid drop model, experiments with giant resonances and heavy-ion collisions. The objective is to identify those parameterizations which best satisfy the current understanding of the physics of nuclear matter over a wide range of applications. Out of the 240 models, only 16 are shown to satisfy all these constraints. Additional, more microscopic, constraints on density dependence of the neutron and proton effective mass beta-equilibrium matter, Landau parameters of symmetric and pure neutron nuclear matter, and observational data on high- and low-mass cold neutron stars further reduce this number to 5, a very small group of recommended Skyrme parameterizations to be used in future applications of the Skyrme interaction of nuclear matter related observables. Full information on partial fulfillment of individual constraints by all Skyrme models considered is given. The results are discussed in terms of the physical interpretation of the Skyrme interaction and the validity of its use in mean-field models. Future work on application of the Skyrme forces, selected on the basis of variables of nuclear matter, in Hartree-Fock calculation of properties of finite nuclei, is outlined.

Planetary transits provide a unique opportunity to investigate the surface distributions of star spots. Our aim is to determine if, with continuous observation (such as the data that will be provided by the Kepler mission), we can in addition measure the rate of drift of the spot belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, incorporating both flux emergence and surface transport. This provides the radial magnetic field distribution on the stellar surface as a function of time. We then model the transit of a planet whose orbital axis is misaligned with the stellar rotation axis. Such a planet could occult spots at a range of latitudes. This allows us to complete the forward modelling of the shape of the transit lightcurve. We then attempt the inverse problem of recovering spot locations from the transit alone. From this we determine if transit lightcurves can be used to measure spot belt locations as a function of time. We find that for low-activity stars such as the Sun, the 3.5 year Kepler window is insufficient to determine this drift rate. For more active stars, it may be difficult to distinguish subtle differences in the nature of flux emergence, such as the degree of overlap of the “butterfly wings”. The rate and direction of drift of the spot belts can however be determined for these stars. This would provide a critical test of dynamo theory.

Spitzer IRAC observations of 15 metal-polluted white dwarfs reveal infrared excesses in the spectral energy distributions of HE 0110-5630, GD 61, and HE 1349-2305. All three of these stars have helium-dominated atmospheres, and their infrared emissions are consistent with warm dust produced by the tidal destruction of (minor) planetary bodies. This study brings the number of metal-polluted, helium and hydrogen atmosphere white dwarfs surveyed with IRAC to 53 and 38 respectively. It also nearly doubles the number of metal-polluted helium-rich white dwarfs found to have closely orbiting dust by Spitzer. From the increased statistics for both atmospheric types with circumstellar dust, we derive a typical disk lifetime of log[t_{disk} (yr)] = 5.6+-1.1 (ranging from 3*10^4 – 5*10^6 yr). This assumes a relatively constant rate of accretion over the timescale where dust persists, which is uncertain. We find that the fraction of highly metal-polluted helium-rich white dwarfs that have an infrared excess detected by Spitzer is only 23 per cent, compared to 48 per cent for metal-polluted hydrogen-rich white dwarfs, and we conclude from this difference that the typical lifetime of dusty disks is somewhat shorter than the diffusion time scales of helium-rich white dwarf. We also find evidence for higher time-averaged accretion rates onto helium-rich stars compared to the instantaneous accretion rates onto hydrogen-rich stars; this is an indication that our picture of evolved star-planetary system interactions is incomplete. We discuss some speculative scenarios that can explain the observations.

Starting from June 2011, the activity of a 137Cs source has been measured by means of a HPGe detector installed deep underground in the Gran Sasso Laboratory. In total about 5100 energy spectra, one hour measuring time each, have been collected. These data allowed the search for time variations of the decay constant with periods from a few hours to 1 year. No signal with amplitude larger than 9.6{\cdot}10-5 at 95% C.L. has been detected. These limits are more than one order of magnitude lower than the values on the oscillation amplitude reported in literature. In particular, for 1 year period an oscillation amplitude larger than 8.5{\cdot}10-5 has been excluded at 95% C.L., independently of the phase. The same data give a value of 29.96{\pm}0.08 years for the 137Cs half life, in good agreement with the world mean value of 30.05{\pm}0.08 years.

We report on extended photometry of two pulsating sdB stars in close binaries. For both cases, we use rotational splitting of the pulsation frequencies to show that the sdB component rotates much too slowly to be in synchronous rotation. We use a theory of tidal interaction in binary stars to place limits on the mass ratios that are independent of estimates based on the radial velocity curves. The companions have masses below 0.26 M\odot. The pulsation spectra show the signature of high-overtone g-mode pulsation. One star, KIC 11179657, has a clear sequence of g-modes with equal period spacings as well as several periodicities that depart from that trend. KIC 02991403 shows a similar sequence, but has many more modes that do not fit the simple pattern.

Context. This paper investigates the role of the Hall term in the propagation and dissipation of waves which interact with 2D magnetic X-points and considers the effect of the Hall term on the nature of the resulting reconnection. Aims. The goal is to determine how the evolution of a nonlinear fast magnetoacoustic wave pulse, and the behaviour of the oscillatory reconnection which results from the interaction of the pulse with a line-tied 2D magnetic X-point, is affected by the Hall term in the generalised Ohm’s law. Methods. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth (di) in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Results. In general, the fast wave is coupled to a shear wave and, for finite di, to whistler and ion cyclotron waves. Dispersive whistler effects cause rapid oscillations of the X-point, which (in combination with the arrival of the main body of the pulse) leads to the creation of magnetic islands and multiple null points under the influence of the Hall term. At later times, competition of local Lorentz and gas pressure forces return the system to a near-equilibrium state. The rate of oscillatory reconnection recovered during this latter phase appears to be unaffected by the value of di.

Aims: Jets are excellent signposts for very young embedded protostars, so we want to identify jet-driving protostars as a tracer of the currently forming generation of stars in the Carina Nebula, which is one of the most massive galactic star-forming regions and which is characterised by particularly high levels of massive-star feedback on the surrounding clouds. Methods: We used archive data to construct large (> 2 deg x 2 deg) Spitzer IRAC mosaics of the Carina Nebula and performed a spatially complete search for objects with excesses in the 4.5 micron band, typical of shock-excited molecular hydrogen emission. We also identified the mid-infrared point sources that are the likely drivers of previously discovered Herbig-Haro jets and molecular hydrogen emission line objects. We combined the Spitzer photometry with our recent Herschel far-infrared data to construct the spectral energy distributions, and used the Robitaille radiative-transfer modelling tool to infer the properties of the objects. Results: The radiative-transfer modelling suggests that the jet sources are protostars with masses between ~1 M_sol and ~10 M_sol that are surrounded by circumstellar disks and embedded in circumstellar envelopes. Conclusions: The estimated protostar masses < 10 M_sol suggest that the current star-formation activity in the Carina Nebula is restricted to low- and intermediate-mass stars. More optical than infrared jets can be observed, indicating that star formation predominantly takes place close to the surfaces of clouds.

We monitored the Ca II H and K lines of HD 179949, a notable star in the southern hemisphere, to observe and confirm previously identified planet induced emission (PIE) as an effect of star-planet interaction. We obtained high resolution spectra (R ~ 53,000) with a signal-to-noise ratio S/N >~ 50 in the Ca II H and K cores during 10 nights of observation at the McDonald Observatory. Wide band echelle spectra were taken using the 2.7 m telescope. Detailed statistical analysis of Ca II K revealed fluctuations in the Ca II K core attributable to planet induced chromospheric emission. This result is consistent with previous studies by Shkolnik et al. (2003). Additionally, we were able to confirm the reality and temporal evolution of the phase shift of the maximum of star-planet interaction previously found. However, no identifiable fluctuations were detected in the Ca II H core. The Al I lambda 3944 A line was also monitored to gauge if the expected activity enhancements are confined to the chromospheric layer. Our observations revealed some variability, which is apparently unassociated with planet induced activity.

One of the most challenging difficulties that precedes the frequency analysis of {\it Kepler} data for a Blazhko star is stitching together the data from different seasons (quarters). We discuss the preliminary steps in the stitching, detrending and rescaling process using the data for long-term Blazhko stars. We present the process on {\it Kepler} data of a Blazhko star with a variable Blazhko cycle and some first results of our analysis.

We present BVR light curves of pulsating stars, ST Boo and RR Leo, obtained between March and September 2007 at the Ankara University Observatory (AUG) and the T\”{U}B{\.I}TAK National Observatory (TUG). Although these observational data are insufficient to obtain the reliable results for a frequency analysis of ST Boo and RR Leo stars, in this study, we tried to investigate the pulsation phenomena of these two stars, as an overview, using the Period04 software package. As preliminary results, we present the possible frequencies for ST Boo and RR Leo.

The SDO/HMI instruments provide photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semi-analytic and numeric equilibria and been applied before to vector magnetograms from Hinode and ground based observations. Recently we implemented a new version which takes measurement errors in photospheric vector magnetograms into account. Photospheric field measurements are often due to measurement errors and finite nonmagnetic forces inconsistent as a boundary for a force-free field in the corona. In order to deal with these uncertainties, we developed two improvements: 1.) Preprocessing of the surface measurements in order to make them compatible with a force-free field 2.) The new code keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which have to be optimized for use with data from SDO/HMI. Within this work we describe the corresponding analysis method and evaluate the force-free equilibria by means of how well force-freeness and solenoidal conditions are fulfilled, the angle between magnetic field and electric current and by comparing projections of magnetic field lines with coronal images from SDO/AIA. We also compute the available free magnetic energy and discuss the potential influence of control parameters.

Context: The study of variations in total solar irradiance (TSI) is important for understanding how the Sun affects the Earth’s climate. Aims: Full-disk continuum images and magnetograms are now available for three full solar cycles. We investigate how modelled TSI compares with direct observations by building a consistent modelled TSI dataset. The model, based only on changes in the photospheric magnetic flux can then be tested on rotational, cyclical and secular timescales. Methods: We use Kitt Peak and SoHO/MDI continuum images and magnetograms in the SATIRE-S model to reconstruct TSI over cycles 21-23. To maximise independence from TSI composites, SORCE/TIM TSI data are used to fix the one free parameter of the model. We compare and combine the separate data sources for the model to estimate an uncertainty on the reconstruction and prevent any additional free parameters entering the model. Results: The reconstruction supports the PMOD composite as being the best historical record of TSI observations, although on timescales of the solar rotation the IRMB composite provides somewhat better agreement. Further to this, the model is able to account for 92% of TSI variations from 1978 to 2009 in the PMOD composite and over 96% during cycle 23. The reconstruction also displays an inter-cycle, secular decline of 0.20 (+0.12 / -0.09) Wm-2 between cycle 23 minima, in agreement with the PMOD composite. Conclusions: SATIRE-S is able to recreate TSI observations on all timescales of a day and longer over 31 years from 1978. This is strong evidence that changes in photospheric magnetic flux alone are responsible for almost all solar irradiance variations over the last three solar cycles.

Le Borgne et al. (2007) report on the determination of evolutionary changes in the periods of field RR Lyr stars. Thanks to the extension of the GEOS database, we could analyze a sample twice larger than the previous one. We obtained a different picture of the period changes, with a number of stars showing an increasing period greater than that of stars showing a decreasing period.

CoRoT photometric measurements of asteroseismic targets need complementary ground-based spectroscopic observations. We are using the planet-hunter HARPS spectrograph attached to the 3.6m-ESO telescope in the framework of two consecutive Large Programmes. We discuss its use to study line-profile variations and we report on a specific result obtained for the Delta Sct star HD 170699.

We report time-resolved photometry of the cataclysmic variable V849 Her, and measure a period of 0.1414 \pm 0.0030 days (3.394 \pm 0.072 hours). We also present photometry taken over several weeks in 2010 and 2011, as well as light curves from 1995 to 2011 by the American Association of Variable Star Observers. The spectra, absolute magnitude derived from infrared magnitudes, and variability all suggest that V849 Her is a nova-like variable. The shallow (0.5-magnitude) low states we observe resemble the erratic low states of the VY Sculptoris stars, although they may recur quasi-periodically over an average cycle of 12.462 \pm 0.074 days.

Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for the turbulent dynamics can be built. Numerical simulations of MHD turbulence are widely used to guide and test the theoretical predictions; however, simulating MHD turbulence and accurately measuring its scaling properties is far from straightforward. Computational power limits the calculations to moderate Reynolds numbers and often simplifying assumptions are made in order that a wider range of scales can be accessed. After describing the theoretical predictions and the numerical approaches that are often employed in studying strong incompressible MHD turbulence, we present the findings of a series of high-resolution direct numerical simulations. We discuss the effects that insufficiencies in the computational approach can have on the solution and its physical interpretation.

We examine the observational consequences of partial gaps being opened by planets in protoplanetary disks. We model the disk using a static alpha-disk model with detailed radiative transfer, parametrizing the shape and size of the partially cleared gaps based on the results of hydrodynamic simulations. Shadowing and illumination by stellar irradiation at the surface of the gap leads to increased contrast as the gap trough is deepened by shadowing and cooling and the far gap wall is puffed up by illumination and heating. In calculating observables, we find that multiple scattering is important and derive an approximation to include these effects. A gap produced by a 200 M_Earth (70 M_Earth) planet at 10 AU can lower/raise the midplane temperature of the disk by up to ~-25/+29% (~-11/+19%) by shadowing in the gap trough and illumination on the far shoulder of the gap. At the distance of Taurus, this gap would be resolvable with ~0.01″ angular resolution. The gap contrast is most significant in scattered light and at thermal continuum wavelengths characteristic of the surface temperature, reducing or raising the surface brightness by up to order of magnitude. Since gaps sizes are correlated to planet mass, this is a promising way of finding and determining the masses of planets embedded in protoplanetary disks.

Elemental and isotopic abundances are the fossils of galactic archaeology. The observed [X/Fe]-[Fe/H] relations in the Galactic bulge and disk and the mass-metallicity relation of galaxies are roughly reproduced with chemodynamical simulations of galaxies under the standard \Lambda-CDM picture and standard stellar physics. The isotopic ratios such as ^{17,18}O and ^{25,26}Mg may require a refinement of modelling of supernova and asymptotic giant branch stars. The recent observation of the Carbon-rich damped Lyman \alpha system can be reproduced only with faint core-collapse supernovae. This suggests that chemical enrichment by the first stars in the first galaxies is driven not by pair-instability supernovae but by core-collapse supernovae (\sim 20-50M_\odot). The observed F abundances can be reproduced with the neutrino processes of core-collapse supernovae. As in F, the observations of elemental abundances in small systems may requires further complications of chemical enrichment. In globular clusters the relative contribution from low-mass supernovae is likely to be smaller than in the field, while the contribution from massive supernovae seems smaller in dwarf spheroidal galaxies than in the solar neighbourhood.

Aims: We demonstrate the universal character of the quiet-Sun chromosphere among inactive stars (solar-type and giants). By assessing the main physical processes, we shed new light on some common observational phenomena. Methods: We discuss measurements of the solar Mt. Wilson S-index, obtained by the Hamburg Robotic Telescope around the extreme minimum year 2009, and compare the established chromospheric basal Ca II K line flux to the Mt. Wilson S-index data of inactive (“flat activity”) stars, including giants. Results: During the unusually deep and extended activity minimum of 2009, the Sun reached S-index values considerably deeper than in any of its previously observed minima. In several brief periods, the Sun coincided exactly with the S-indices of inactive (“flat”, presumed Maunder Minimum-type) solar analogues of the Mt. Wilson sample; at the same time, the solar visible surface was also free of any plages or remaining weak activity regions. The corresponding minimum Ca II K flux of the quiet Sun and of the presumed Maunder Minimum-type stars in the Mt. Wilson sample are found to be identical to the corresponding Ca II K chromospheric basal flux limit. Conclusions: We conclude that the quiet-Sun chromosphere is a universal phenomenon among inactive stars. Its mixed-polarity magnetic field, generated by a local, “fast” turbulent dynamo finally provides a natural explanation for the minimal soft X-ray emission observed for inactive stars. Given such a local dynamo also works for giant chromospheres, albeit on larger length scales, i.e., l ~ R/g, with R and g as stellar radius and surface gravity, respectively, the existence of giant spicular phenomena and the guidance of mechanical energy toward the acceleration zone of cool stellar winds along flux-tubes have now become traceable.

Much of our knowledge of galaxies comes from analysing the radiation emitted by their stars. It depends on the stellar initial mass function (IMF) describing the distribution of stellar masses when the population formed. Consequently knowledge of the IMF is critical to virtually every aspect of galaxy evolution. More than half a century after the first IMF determination, no consensus has emerged on whether it is universal in different galaxies. Previous studies indicated that the IMF and the dark matter fraction in galaxy centres cannot be both universal, but they could not break the degeneracy between the two effects. Only recently indications were found that massive elliptical galaxies may not have the same IMF as our Milky Way. Here we report unambiguous evidence for a strong systematic variation of the IMF in early-type galaxies as a function of their stellar mass-to-light ratio, producing differences up to a factor of three in mass. This was inferred from detailed dynamical models of the two-dimensional stellar kinematics for the large Atlas3D representative sample of nearby early-type galaxies spanning two orders of magnitude in stellar mass. Our finding indicates that the IMF depends intimately on a galaxy’s formation history.

Calculating the microscopic dissociation rate of a bound state, such as a classical diatomic molecule, has been difficult so far. The problem was that standard theories require an energy barrier over which the bound particle (or state) escapes into the preferred low-energy state. This is not the case when the long-range repulsion responsible for the barrier is either absent or screened (as in Cooper pairs, ionized plasma, or biomolecular complexes). We solve this classical problem by accounting for entropic memory at the microscopic level. The theory predicts dissociation rates for arbitrary potentials and is successfully tested on the example of plasma, where it yields an estimate of ionization in the core of Sun in excellent agreement with experiments. In biology, the new theory accounts for crowding in receptor-ligand kinetics and protein aggregation.

Mennickent et al.and Sabogal et al.identified a large number of Classical Be (CBe) candidates in the L&SMC based on their photometric variability using the OGLEII database. They classified these stars into four different groups based on the appearance of their variability. We studied the infrared properties of the sample as well as the spectroscopic properties of a subsample. We cross-correlated the optical sample with the IRSF catalog to obtain the J, H, Ks magnitudes of all the four types of stars in the L&SMC. Spectra of 120 stars belonging to the types 1, 2 and 3 were analysed to study their spectral properties. Among the four types, the type 4 stars is the dominant group. The NIR colour-colour diagrams suggest that the type 4 stars in the LMC have a subclass, which is not found in our Galaxy or in the SMC. The main type 4 sample which is \sim 49% of the total sample has NIR properties similar to the Galactic CBe stars and the SMC type 4 stars. Though the new subclass of type 4 stars have high E(B – V) \sim 0.75, they are not located close to regions with high reddening. The type 3 stars (\sim 6% & 7.3% in the L&SMC) are found to have large H{\alpha} EW in the SMC and some are found to have large NIR excess. This small fraction of stars are unlikely to be CBe stars. The type 2 stars are found in larger fraction in the SMC, when compared to the LMC. The spectroscopic and the NIR properties suggest that these could be CBe stars. The spectroscopic sample of type 1 stars which show H{\alpha} in emission and confirmed as CBe stars are more abundant in the SMC by a factor of 2.6. If the effect of metallicity is to cause more CBe stars in the SMC, when compared to the LMC, then type 1, type 2 and type 4 stars follow this rule, with an enhancement of 2.6, 2.4 and 1.3 respectively.

We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetised gas cloud. We focus on the dependence of the growth and saturation level of the dynamo generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose an initial field strength of $\simeq 1\,\mkG$ with the magnetic to kinetic energy ratio, $E_{\rm m}/E_{\rm k} \sim 10^{-4}$. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of $B/\rho^{2/3}$ and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from $v_{\rm rms}/c_{\rm s}\sim 0.2$ to 0.4 and then starts decreasing from $v_{\rm rms}/c_{\rm s}\sim 1.0$. We obtain saturation values of $E_{\rm m}/E_{\rm k} = 0.2 – 0.3$ for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at $E_{\rm m}/E_{\rm k}\sim 0.2 – 0.4$ of the equipartition value. (Abridged)

Line ratios in “fir” triplets of helium-like ions have proven to be a powerful diagnostic of conditions in X-ray emitting gas surrounding massive stars. Recent observations indicate that these ratios can be variable with time. The possible causes of variation in line ratios are limited: changes in the radiation field or changes in density, and changes in mass-loss or geometry. In this paper, we investigate the ability of changes in the radiation field to induce variability in the ratio R=f/i. To isolate the radiative effect, we use a heuristic model of temperature and radius changes in variable stars in the B and O range with low-density, steady-state winds. We model the changes in emissivity of X-ray emitting gas close to the star due to differences in level-pumping from available UV photons at the location of the gas. We find that under these conditions, variability in R is dominated by the stellar temperature. Although the relative amplitude of variability is roughly comparable for most lines at most temperatures, detectable variations are limited to a few lines for each spectral type. We predict that variable values in R due to stellar variability must follow predictable trends found in our simulations. Our model uses radial pulsations as a mode of stellar variability that maximizes the amplitude of variation in R. This model is robust enough to show which ions will provide the best opportunity for observing variability in the f/i ratio at different stellar temperatures, and the correlation of that variability with other observable parameters. In real systems, the effects would be more complex than in our model, with differences in phase and suppressed amplitude in the presence of non-radial pulsations. This suggests that changes in R across many lines concurrently are not likely to be produced by a variable radiation field.

Extrasolar planets and belts of debris orbiting post-main-sequence single stars may become unbound as the evolving star loses mass. In multiple star systems, the presence or co-evolution of the additional stars can significantly complicate the prospects for orbital excitation and escape. Here, we investigate the dynamical consequences of multi-phasic, nonlinear mass loss and establish a criterion for a system of any stellar multiplicity to retain a planet whose orbit surrounds all of the parent stars. For single stars which become white dwarfs, this criterion can be combined with the Chandrasekhar Limit to establish the maximum allowable mass loss rate for planet retention. We then apply the criterion to circumbinary planets in evolving binary systems over the entire stellar mass phase space. Through about 10^5 stellar evolutionary track realizations, we characterize planetary ejection prospects as a function of binary separation, stellar mass and metallicity. This investigation reveals that planets residing at just a few tens of AU from a central concentration of stars are susceptible to escape in a wide variety of multiple systems. Further, planets are significantly more susceptible to ejection from multiple star systems than from single star systems for a given system mass. For system masses greater than about 2 Solar masses, multiple star systems represent the greater source of free-floating planets.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula HaTr 4, one of few known to contain a post-common-envelope central star system. Common envelope evolution is believed to play an important role in the shaping of planetary nebulae, but the exact nature of this role is yet to be understood. High spatial- and spectral- resolution spectroscopy of the [OIII]5007 nebular line obtained with VLT-UVES are presented alongside deep narrowband Ha+[NII]6584 imagery obtained using EMMI-NTT, and together the two are used to derive the three-dimensional morphology of HaTr 4. The nebula is found to display an extended ovoid morphology with an enhanced equatorial region consistent with a toroidal waist – a feature believed to be typical amongst planetary nebulae with post-common-envelope central stars. The nebular symmetry axis is found to lie perpendicular to the orbital plane of the central binary, concordant with the idea that the formation and evolution of HaTr 4 has been strongly influenced by its central binary.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula HaTr 4, one of few known to contain a post-common-envelope central star system. Common envelope evolution is believed to play an important role in the shaping of planetary nebulae, but the exact nature of this role is yet to be understood. High spatial- and spectral- resolution spectroscopy of the [OIII]5007 nebular line obtained with VLT-UVES are presented alongside deep narrowband Ha+[NII]6584 imagery obtained using EMMI-NTT, and together the two are used to derive the three-dimensional morphology of HaTr 4. The nebula is found to display an extended ovoid morphology with an enhanced equatorial region consistent with a toroidal waist – a feature believed to be typical amongst planetary nebulae with post-common-envelope central stars. The nebular symmetry axis is found to lie perpendicular to the orbital plane of the central binary, concordant with the idea that the formation and evolution of HaTr 4 has been strongly influenced by its central binary.

We report on a BVRcIc survey of field W Ursae Majoris binary stars and present accurate colors for 606 systems that have been observed on at least three photometric nights from a robotic observatory in southern Arizona. Comparison with earlier photometry for a subset of the systems shows good agreement. We investigate two independent methods of determining the interstellar reddening, although both have limitations that can render them less effective than desired. A subset of 101 systems shows good agreement between the two reddening methods.

We present evidence that there is an M dwarf problem similar to the previously identified G dwarf and K dwarf problems: the number of low-metallicity M dwarfs is not sufficient to match simple closed-box models of local Galactic chemical evolution. We estimated the metallicity of 4141 M dwarf stars with spectra from the Sloan Digital Sky Survey (SDSS) using a molecular band strength versus metallicity calibration developed using high resolution spectra of nearby M dwarfs. Using a sample of M dwarfs with measured magnitudes, parallaxes, and metallicities, we derived a relation that describes the absolute magnitude variation as a function of metallicity. When we examined the metallicity distribution of SDSS stars, after correcting for the different volumes sampled by the magnitude-limited survey, we found that there is an M dwarf problem, with the number of M dwarfs at [Fe/H] ~ -0.5 less than 1% the number at [Fe/H] = 0, where a simple model of Galactic chemical evolution predicts a more gradual drop in star numbers with decreasing metallicity.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.