Showing papers in "arXiv: Solar and Stellar Astrophysics in 2013"

Modules for Experiments in Stellar Astrophysics (MESA): Giant Planets, Oscillations, Rotation, and Massive Stars

Abstract: We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]

Stellar Multiplicity

Abstract: Stellar multiplicity is an ubiquitous outcome of the star formation process Characterizing the frequency and main characteristics of multiple systems and their dependencies on primary mass and environment is therefore a powerful tool to probe this process While early attempts were fraught with selection biases and limited completeness, instrumentation breakthroughs in the last two decades now enable robust analyses In this review, we summarize our current empirical knowledge of stellar multiplicity for Main Sequence stars and brown dwarfs, as well as among populations of Pre-Main Sequence stars and embedded protostars Clear trends as a function of both primary mass and stellar evolutionary stage are identified that will serve as a comparison basis for numerical and analytical models of star formation

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

Abstract: We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. The synthetic spectra cover the wavelength range from 500AA to 50.000AA with resolutions of R=500.000 in the optical and near IR, R=100.000 in the IR and a step size of 0.1AA in the UV. The parameter space covers 2.300K<=Teff<=12.000K, 0.0<=log(g)<=+6.0, -4.0<=[Fe/H]<=+1.0, and -0.2<=[alpha/Fe]<=+1.2. The library is a work in progress and we expect to extend it up to Teff=25.000 K.

Intrinsic Colors, Temperatures, and Bolometric Corrections of Pre-Main Sequence Stars

Abstract: We present an analysis of the intrinsic colors and temperatures of 5-30 Myr old pre-main sequence (pre-MS) stars using the F0 through M9 type members of nearby, negligibly reddened groups: Eta Cha cluster, TW Hydra Association, Beta Pic Moving Group, and Tucana-Horologium Association. To check the consistency of spectral types from the literature, we estimate new spectral types for 52 nearby pre-MS stars with spectral types F3 through M4 using optical spectra taken with the SMARTS 1.5-m telescope. Combining these new types with published spectral types, and photometry from the literature (Johnson-Cousins BVIc, 2MASS JHKs and WISE W1, W2, W3, and W4), we derive a new empirical spectral type-color sequence for 5-30 Myr old pre-MS stars. Colors for pre-MS stars match dwarf colors for some spectral types and colors, but for other spectral types and colors, deviations can exceed 0.3 mag. We estimate effective temperatures (Teff) and bolometric corrections (BCs) for our pre-MS star sample through comparing their photometry to synthetic photometry generated using the BT-Settl grid of model atmosphere spectra. We derive a new Teff and BC scale for pre-MS stars, which should be a more appropriate match for T Tauri stars than often-adopted dwarf star scales. While our new Teff scale for pre-MS stars is within ~100 K of dwarfs at a given spectral type for stars

Alfven Wave Solar Model (AWSoM): Coronal Heating

Abstract: We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently.

X-Shooter spectroscopy of young stellar objects: IV -- Accretion in low-mass stars and sub-stellar objects in Lupus

Abstract: We present X-Shooter/VLT observations of a sample of 36 accreting low-mass stellar and sub-stellar objects (YSOs) in the Lupus star forming region, spanning a range in mass from ~0.03 to ~1.2Msun, but mostly with 0.1Msun < Mstar < 0.5Msun. Our aim is twofold: firstly, analyse the relationship between excess-continuum and line emission accretion diagnostics, and, secondly, to investigate the accretion properties in terms of the physical properties of the central object. The accretion luminosity (Lacc), and from it the accretion rate (Macc), is derived by modelling the excess emission, from the UV to the near-IR, as the continuum emission of a slab of hydrogen. The flux and luminosity (Ll) of a large number of emission lines of H, He, CaII, etc., observed simultaneously in the range from ~330nm to 2500nm, were computed. The luminosity of all the lines is well correlated with Lacc. We provide empirical relationships between Lacc and the luminosity of 39 emission lines, which have a lower dispersion as compared to previous relationships in the literature. Our measurements extend the Pab and Brg relationships to Lacc values about two orders of magnitude lower than those reported in previous studies. We confirm that different methodologies to measure Lacc and Macc yield significantly different results: Ha line profile modelling may underestimate Macc by 0.6 to 0.8dex with respect to Macc derived from continuum-excess measures. Such differences may explain the likely spurious bi-modal relationships between Macc and other YSOs properties reported in the literature. We derive Macc in the range 2e-12 -- 4e-8 Msun/yr and conclude that Macc is proportional to Mstar^1.8(+/-0.2), with a dispersion lower by a factor of about 2 than in previous studies. A number of properties indicate that the physical conditions of the accreting gas are similar over more than 5 orders of magnitude in Macc.

Asteroseismic fundamental properties of solar-type stars observed by the NASA Kepler Mission

Abstract: We use asteroseismic data obtained by the NASA Kepler Mission to estimate the fundamental properties of more than 500 main-sequence and sub-giant stars. Data obtained during the first 10 months of Kepler science operations were used for this work, when these solar-type targets were observed for one month each in a survey mode. Stellar properties have been estimated using two global asteroseismic parameters and complementary photometric and spectroscopic data. Homogeneous sets of effective temperatures were available for the entire ensemble from complementary photometry; spectroscopic estimates of T_eff and [Fe/H] were available from a homogeneous analysis of ground-based data on a subset of 87 stars. [Abbreviated version... see paper for full abstract.]

The Mass Dependence Between Protoplanetary Disks and their Stellar Hosts

Abstract: We present a substantial extension of the mm-wave continuum photometry catalog for Taurus circumstellar dust disks. Combining new Submillimeter Array data with measurements in the literature, we construct a mm-wave luminosity distribution for Class II disks that is statistically complete for stellar hosts with spectral types earlier than M8.5 and has a (3-sigma) depth of ~3 mJy. The resulting census eliminates a longstanding bias against disks with late-type hosts, and thereby reveals a strong correlation between L_mm and the host spectral type. We confirm that this corresponds to a statistically robust relationship between the masses of dust disks and the stars that host them. A Bayesian regression technique is used to characterize these relationships: the results indicate a typical 1.3 mm flux density of 25 mJy for solar mass hosts and a power-law scaling L_mm \propto M_star^1.5-2.0. We suggest that a reasonable treatment of dust temperature in the conversion from L_mm to M_disk favors an inherently linear M_disk \propto M_star scaling, with a typical disk-to-star mass ratio of $\sim$0.2--0.6%. The RMS dispersion around this regression is 0.7 dex, suggesting that the combined effects of diverse evolutionary states, dust opacities, and temperatures in these disks imprint a FWHM range of a factor of 40 on the inferred M_disk (or L_mm) at any given host mass. We argue that this relationship between M_disk and M_star likely represents the origin of the inferred correlation between giant planet frequency and host star mass in the exoplanet population, and provides some basic support for the core accretion model for planet formation. Moreover, we caution that selection bias must be considered in comparative studies of disk evolution, and illustrate that fact with statistical comparisons of L_mm between Taurus and other clusters (abridged).

The Neutron Star Mass Distribution

Abstract: In recent years, the number of pulsars with secure mass measurements has increased to a level that allows us to probe the underlying neutron star mass distribution in detail. We critically review radio pulsar mass measurements and present a detailed examination through which we are able to put stringent constraints on the underlying neutron star mass distribution. For the first time, we are able to analyze a sizable population of neutron star-white dwarf systems in addition to double neutron star systems with a technique that accounts for systematically different measurement errors. We find that neutron stars that have evolved through different evolutionary paths reflect distinctive signatures through dissimilar distribution peak and mass cutoff values. Neutron stars in double neutron star and neutron star-white dwarf systems show consistent respective peaks at 1.35 Msun and 1.50 Msun which suggest significant mass accretion (Delta m~0.15 Msun) has occurred during the spin up phase. The width of the mass distribution implied by double neutron star systems is indicative of a tight initial mass function while the inferred mass range is significantly wider for neutron stars that have gone through recycling. We find a mass cutoff at 2 Msun for neutron stars with white dwarf companions which establishes a firm lower bound for the maximum neutron star mass. This rules out the majority of strange quark and soft equation of state models as viable configurations for neutron star matter. The lack of truncation close to the maximum mass cutoff suggests that the 2 Msun limit is set by evolutionary constraints rather than nuclear physics or general relativity, and the existence of rare super-massive neutron stars is possible.

Equation of state and neutron star properties constrained by nuclear physics and observation

Abstract: Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 +- 0.04 M_sun neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the nuclear symmetry energy. In particular, we show that the predicted range for neutron star radii is robust. For use in astrophysical simulations, we provide detailed numerical tables for a representative set of equations of state consistent with these constraints.

Populations of rotating stars. - I. Models from 1.7 to 15 Msun at Z = 0.014, 0.006, and 0.002 with {\Omega}/{\Omega}crit between 0 and 1

Abstract: B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars. In this paper, we provide stellar models covering the mass range between 1.7 to 15 Msun, which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence. We used the Geneva stellar-evolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 Msun present at the beginning of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. The values presently obtained are in agreement with the very few values obtained for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 Msun with a mean velocity on the MS of the order of 150 km/s show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.

Superflares on Solar Type Stars Observed with Kepler I. Statistical Properties of Superflares

Abstract: By extending our previous study by Maehara et al. (2012), we searched for superflares on G-type dwarfs (solar type stars) using Kepler data for a longer period (500 days) than that (120 days) in our previous study. As a result, we found 1547 superflares on 279 G-type dwarfs, which are much more than previous 365 superflares on 148 stars. Using these new data, we studied the statistical properties of occurrence frequency of superflares, and basically confirmed the previous results, i.e., the occurrence frequency (dN/dE) of superflares vs flare energy (E) shows power-law distribution with dN/dE \propto E^{-\alpha}, where \alpha ~ 2. It is interesting that this distribution is roughly on the same line as that for solar flares. In the case of the Sun-like stars (with surface temperature 5600-6000K and slowly rotating with a period longer than 10 days), the occurrence frequency of superflares with energy of 10^34 -10^35 erg is once in 800-5000 years. We also studied long term (500 days) stellar brightness variation of these superflare stars, and found that in some G-type dwarfs the occurrence frequency of superflares was extremely high, ~ 57 superflares in 500 days (i.e., once in 10 days). In the case of Sun-like stars, the most active stars show the frequency of one superflares (with 10^34 erg) in 100 days. There is evidence that these superflares have extremely large starspots with a size about 10 times larger than that of the largest sunspot. We argue that the physical origin of extremely high occurrence frequency of superflares in these stars may be attributed to the existence of extremely large starspots.

SASI Activity in Three-Dimensional Neutrino-Hydrodynamics Simulations of Supernova Cores

Abstract: The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.

Stellar Diameters and Temperatures III. Main Sequence A, F, G, & K Stars: Additional high-precision measurements and empirical relations

Abstract: Based on CHARA Array measurements, we present the angular diameters of 23 nearby, main- sequence stars, ranging from spectral type A7 to K0, five of which are exoplanet host stars. We derive linear radii, effective temperatures, and absolute luminosities of the stars using HIPPARCOS parallaxes and measured bolometric fluxes. The new data are combined with previously published values to create an Angular Diameter Anthology of measured angular diameters to main-sequence stars (luminosity class V and IV). This compilation consists of 125 stars with diameter uncertainties of less than 5%, ranging in spectral types from A to M. The large quantity of empirical data are used to derive color-temperature relations to an assortment of color indices in the Johnson (BVRIJHK), Cousins (RI), Kron (RI), Sloan (griz), and WISE (W3W4) photometric systems. These relations have an average standard deviation of ~3% and are valid for stars with spectral types A0 to M4. To derive even more accurate relations for Sun-like stars, we also determined these temperature relations omitting early-type stars (Teff > 6750 K) that may have biased luminosity estimates because of rapid rotation; for this subset the dispersion is only ~2.5%. We find effective temperatures in agreement within a couple percent for the interferometrically characterized sample of main sequence stars compared to those derived via the infrared-flux method and spectroscopic analysis.

Asteroseismic classification of stellar populations among 13000 red giants observed by Kepler

Abstract: Of the more than 150000 targets followed by the Kepler Mission, about 10% were selected as red giants. Due to their high scientific value, in particular for Galaxy population studies and stellar structure and evolution, their Kepler light curves were made public in late 2011. More than 13000 (over 85%) of these stars show intrinsic flux variability caused by solar-like oscillations making them ideal for large scale asteroseismic investigations. We automatically extracted individual frequencies and measured the period spacings of the dipole modes in nearly every red giant. These measurements naturally classify the stars into various populations, such as the red giant branch, the low-mass (M/Msol 1.8) secondary clump. The period spacings also reveal that a large fraction of the stars show rotationally induced frequency splittings. This sample of stars will undoubtedly provide an extremely valuable source for studying the stellar population in the direction of the Kepler field, in particular when combined with complementary spectroscopic surveys.

Reconstructing the Subsurface Three-Dimensional Magnetic Structure of A Solar Active Region Using SDO/HMI Observations

Abstract: A solar active region (AR) is a three-dimensional magnetic structure formed in the convection zone, whose property is fundamentally important for determining the coronal structure and solar activity when emerged. However, our knowledge on the detailed 3-D structure prior to its emergence is rather poor, largely limited by the low cadence and sensitivity of previous instruments. Here, using the 45-second high-cadence observations from the Helioseismic and Magnetic Imager (\emph{HMI}) onboard the Solar Dynamics Observatory (\emph{SDO}), we are able for the first time to reconstruct a 3-D datacube and infer the detailed subsurface magnetic structure of NOAA AR 11158 and to characterize its magnetic connectivity and topology. This task is accomplished with the aid of the image-stacking method and advanced 3-D visualization. We find that the AR consists of two major bipoles, or four major polarities. Each polarity in 3-D shows interesting tree-like structure, i.e. while the root of the polarity appears as a single tree-trunk-like tube, the top of the polarity has multiple branches consisting of smaller and thinner flux-tubes which connect to the branches of the opposite polarity that is similarly fragmented. The roots of the four polarities align well along a straight line, while the top branches are slightly non-coplanar. Our observations suggest that an active region, even appearing highly complicated on the surface, may originate from a simple straight flux-tube that undergoes both horizontal and vertical bifurcation processes during its rise through the convection zone.

Unveiling the gas and dust disk structure in HD 163296 using ALMA observations

Abstract: Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unveils for the first time at sub-millimeter frequencies the vertical structure details of a gaseous disk in Keplerian rotation, showing the back- and the front-side of a flared disk. Continuum emission at 850 microns reveals a compact dust disk with a 240 AU outer radius and a surface brightness profile that shows a very steep decline at radius larger than 125 AU. The gaseous disk is more than two times larger than the dust disk, with a similar critical radius but with a shallower radial profile. Radiative transfer models of the continuum data confirms the need for a sharp outer edge to the dust disk. The models for the CO(3-2) channel map require the disk to be slightly more geometrically thick than previous models suggested, and that the temperature at which CO gas becomes depleted (frozen-out) from the outer regions of the disk midplane is T < 20 K, in agreement with previous studies.

Local Enhancement of Surface Density in the Protoplanetary Ring Surrounding HD 142527

Abstract: We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

A WFC3/HST view of the three stellar populations in the Globular Cluster NGC6752

Abstract: Multi-band Hubble Space Telescope photometry reveals that the main sequence, sub-giant, and the red giant branch of the globular cluster NGC6752 splits into three main components in close analogy with the three distinct segments along its horizontal branch stars. These triple sequences are consistent with three stellar groups: a stellar population with a chemical composition similar to field halo stars (population a), a population (c) with enhanced sodium and nitrogen, depleted carbon and oxygen and enhanced helium abundance (Delta Y ~0.03), and a population (b) with an intermediate (between population a and c) chemical composition and slightly helium enhanced (Delta Y ~0.01). These components contain ~25% (population a), ~45% (population b), and ~30% (population c) of the stars. No radial gradient for the relative numbers of the three populations has been identified out to about 2.5 half mass radii.

Scaling of the electron dissipation range of solar wind turbulence

Abstract: Electron scale solar wind turbulence has attracted great interest in recent years. Clear evidences have been given from the Cluster data that turbulence is not fully dissipated near the proton scale but continues cascading down to the electron scales. However, the scaling of the energy spectra as well as the nature of the plasma modes involved at those small scales are still not fully determined. Here we survey 10 years of the Cluster search-coil magnetometer (SCM) waveforms measured in the solar wind and perform a statistical study of the magnetic energy spectra in the frequency range [$1, 180$]Hz. We show that a large fraction of the spectra exhibit clear breakpoints near the electon gyroscale $\rho_e$, followed by steeper power-law like spectra. We show that the scaling below the electron breakpoint cannot be determined unambiguously due to instrumental limitations that will be discussed in detail. We compare our results to recent ones reported in other studies and discuss their implication on the physical mechanisms and the theoretical modeling of energy dissipation in the SW.

Plasma Jets and Eruptions in Solar Coronal Holes: a 3D flux emergence experiment

Abstract: A three-dimensional numerical experiment of the launching of a hot and fast coronal jet followed by several violent eruptions is analyzed in detail. These events are initiated through the emergence of a magnetic flux rope from the solar interior into a coronal hole. We explore the evolution of the emerging magnetically-dominated plasma dome surmounted by a current sheet and the ensuing pattern of reconnection. A hot and fast coronal jet with inverted-Y shape is produced that shows properties comparable to those frequently observed with EUV and X-Ray detectors. We analyze its 3D shape, its inhomogeneous internal structure, and its rise and decay phases, lasting for some 15-20 min each. Particular attention is devoted to the field-line connectivities and the reconnection pattern. We also study the cool and high-density volume that appears encircling the emerged dome. The decay of the jet is followed by a violent phase with a total of five eruptions. The first of them seems to follow the general pattern of tether-cutting reconnection in a sheared arcade, although modified by the field topology created by the preceding reconnection evolution. The two following eruptions take place near and above the strong field-concentrations at the surface. They show a twisted, \Omega-loop like rope expanding in height, with twist being turned into writhe, thus hinting at a kink instability (perhaps combined with a torus-instability) as the cause of the eruption. The succession of a main jet ejection and a number of violent eruptions that resemble mini-CME's and their physical properties suggest that this experiment may provide a model for the blowout jets recently proposed in the literature.

An azimuthal asymmetry in the LkHa 330 disk

Abstract: Theory predicts that giant planets and low mass stellar companions shape circumstellar disks by opening annular gaps in the gas and dust spatial distribution. For more than a decade it has been debated whether this is the dominant process that leads to the formation of transitional disks. In this paper, we present millimeter-wave interferometric observations of the transitional disk around the young intermediate mass star LkHa330. These observations reveal a lopsided ring in the 1.3 mm dust thermal emission characterized by a radius of about 100 AU and an azimuthal intensity variation of a factor of 2. By comparing the observations with a Gaussian parametric model, we find that the observed asymmetry is consistent with a circular arc, that extends azimuthally by about 90 deg and emits about 1/3 of the total continuum flux at 1.3 mm. Hydrodynamic simulations show that this structure is similar to the azimuthal asymmetries in the disk surface density that might be produced by the dynamical interaction with unseen low mass companions orbiting within 70 AU from the central star. We argue that such asymmetries might lead to azimuthal variations in the millimeter-wave dust opacity and in the dust temperature, which will also affect the millimeter-wave continuum emission. Alternative explanations for the observed asymmetry that do not require the presence of companions cannot be ruled out with the existing data. Further observations of both the dust and molecular gas emission are required to derive firm conclusions on the origin of the asymmetry observed in the LkHa330 disk.

The earliest phases of star formation - A Herschel key project. The thermal structure of low-mass molecular cloud cores

Abstract: The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

Advanced burning stages and fate of 8-10 Mo stars

Abstract: The stellar mass range 8

Analytical representations of unified equations of state for neutron-star matter

Abstract: Context. A unified equation of state (EoS) should describe the crust and the core of a neutron star using the same physical model. The Brussels-Montreal group has recently derived a family of such EoSs based on the nuclear energy-density functional theory with generalized Skyrme effective forces, fitted to the available mass data. At the same time, these forces were constrained to reproduce microscopic calculations of homogeneous neutron matter based on realistic two- and three-nucleon forces. Aims. We represent basic physical characteristics of the latest Brussels-Montreal EoS models by analytical expressions to facilitate their inclusion in astrophysical simulations. Methods. We consider three EoS models, which significantly differ by stiffness: BSk19, BSk20, and BSk21. For each of them we constructed two versions of the EoS parametrization. In the first version, pressure P and gravitational mass density \rho are given as functions of the baryon number density n_b. In the second version, P, \rho, and n_b are given as functions of pseudo-enthalpy, which is useful for two-dimensional calculations of stationary rotating configurations of neutron stars. In addition to the EoS, we derived analytical expressions for several related quantities that are required in neutron-star simulations: number fractions of electrons and muons in the stellar core, nucleon numbers per nucleus in the inner crust, and equivalent radii and shape parameters of the nuclei in the inner crust. Results. We obtain analytical representations for the basic characteristics of the models of cold dense matter, which are most important for studies of neutron stars. We demonstrate the usability of our results by applying them to calculations of neutron-star mass-radius relations, maximum and minimum masses, thresholds of direct Urca processes, and the electron conductivity in the neutron-star crust.

The radius discrepancy in low mass stars: single vs. binaries

Abstract: A long-standing issue in the theory of low mass stars is the discrepancy between predicted and observed radii and effective temperatures. In spite of the increasing availability of very precise radius determinations from eclipsing binaries and interferometric measurements of radii of single stars, there is no unanimous consensus on the extent (or even the existence) of the discrepancy and on its connection with other stellar properties (e.g. metallicity, magnetic activity). We investigate the radius discrepancy phenomenon using the best data currently available (accuracy about 5%). We have constructed a grid of stellar models covering the entire range of low mass stars (0.1-1.25 M_sun) and various choices of the metallicity and of the mixing length parameter \alpha. We used an improved version of the Yale Rotational stellar Evolution Code (YREC), implementing surface boundary conditions based on the most up-to-date PHOENIX atmosphere models. Our models are in good agreement with others in the literature and improve and extend the low mass end of the Yale-Yonsei isochrones. Our calculations include rotation-related quantities, such as moments of inertia and convective turnover time scales, useful in studies of magnetic activity and rotational evolution of solar-like stars. Consistently with previous works, we find that both binaries and single stars have radii inflated by about 3% with respect to the theoretical models; among binaries, the components of short orbital period systems are found to be the most deviant. We conclude that both binaries and single stars are comparably affected by the radius discrepancy phenomenon.

Discovery, Progenitor & Early Evolution of a Stripped Envelope Supernova iPTF13bvn

Abstract: The intermediate Palomar Transient Factory reports our discovery of a young supernova, iPTF13bvn, in the nearby galaxy, NGC5806 (22.5Mpc). Our spectral sequence in the optical and infrared suggests a likely Type Ib classification. We identify a single, blue progenitor candidate in deep pre-explosion imaging within a 2{\sigma} error circle of 80 mas (8.7 pc). The candidate has a MB luminosity of -5.2 +/- 0.4 mag and a B-I color of 0.1+/-0.3 mag. If confirmed by future observations, this would be the first direct detection for a progenitor of a Type Ib. Fitting a power law to the early light curve, we find an extrapolated explosion date around 1.1 days before our first detection. We see no evidence of shock cooling. The pre-explosion detection limits constrain the radius of the progenitor to be smaller than a few solar radii. iPTF13bvn is also detected in cm and mm-wavelengths. Fitting a synchrotron self-absorption model to our radio data, we find a mass loading parameter of 1.3*10^12 g/cm. Assuming a wind velocity of 10^3km/s, we derive a progenitor mass loss rate of 3*10^-5Msun/yr. Our observations, taken as a whole, are consistent with a Wolf Rayet progenitor of the supernova iPTF13bvn.

Chromospheric activity as age indicator

Abstract: Chromospheric activity has been calibrated and widely used as an age indicator. However, it has been suggested that the viability of this age indicator is, in the best case, limited to stars younger than about 1.5 Gyr. I aim to define the age range for which chromospheric activity is a robust astrophysical clock. I collected literature measurements of the S-index in field stars, which is a measure of the strength of the H and K lines of the Ca II and a proxy for chromospheric activity, and exploited the homogeneous database of temperature and age determinations for field stars provided by the Geneva-Copenhagen survey of the solar neighbourhood. Field data, inclusive data previously used to calibrate chromospheric ages, confirm the result found using open cluster data, i.e. there is no decay of chromospheric activity after about 2 Gyr. The only existing indication supporting the viability of chromospheric ages older than 2 Gyr is the similarity of chromospheric activity levels in the components of 35 dwarf binaries. However, even in the most optimistic scenario, uncertainty in age determination for field stars and lack of sufficient data in open clusters make any attempt of calibrating an age activity relationship for old stars premature. The hypothesis that chromospheric activity follows the Skumanich law, i.e. that it is proportional to the inverse square root of the age, should be relaxed.

The Compactness of Presupernova Stellar Cores

Abstract: The success or failure of the neutrino-transport mechanism for producing a supernova in an evolved massive star is known to be sensitive not only to the mass of the iron core that collapses, but also to the density gradient in the silicon and oxygen shells surrounding that core. Here we study the systematics of a presupernova core's "compactness" (O'Connor & Ott 2011) as a function of the mass of the star and the physics used in its calculation. Fine-meshed surveys of presupernova evolution are calculated for stars from 15 to 65 Msun. The metallicity and the efficiency of semiconvection and overshoot mixing are both varied and bare carbon-oxygen cores are explored as well as full hydrogenic stars. Two different codes, KEPLER and MESA, are used for the study. A complex interplay of carbon and oxygen burning, especially in shells, can cause rapid variations in the compactness for stars of very nearly the same mass. On larger scales, the distribution of compactness with main sequence mass is found to be robustly non-monotonic, implying islands of "explodability", particularly around 8 to 20 Msun and 25 to 30 Msun. The carbon-oxygen (CO) core mass of a presupernova star is a better, though still ambiguous discriminant of its core structure than the main sequence mass.

Hydrogen Burning on Accreting White Dwarfs: Stability, Recurrent Novae, and the Post-Novae Supersoft Source

Abstract: We examine the properties of white dwarfs (WDs) accreting hydrogen-rich matter in and near the stable burning regime of accretion rates as modeled by time-dependent calculations done with Modules for Experiments in Stellar Astrophysics (MESA). We report the stability boundary for WDs of masses between 0.51 solar masses and 1.34 solar masses as found via time-dependent calculations. We also examine recurrent novae that are accreting at rates close to, but below, the stable burning limit and report their recurrence times and ignition masses. Our dense grid in accretion rates finds the expected minimum possible recurrence times as a function of the WD mass. This enables inferences to be made about the minimum WD mass possible to reach a specific recurrence time. We compare our computational models of post-outburst novae to the stably burning WDs and explicitly calculate the duration and effective temperature (Teff) of the post-novae WD in the supersoft phase. We agree with the measured turnoff time - Teff relation in M31 by Henze and collaborators, infer WD masses in the 1.0-1.3 solar masses range, and predict ejection masses consistent with those observed. We close by commenting on the importance of the hot helium layer generated by stable or unstable hydrogen burning for the short- and long-term evolution of accreting white dwarfs.