# Solar and Stellar Astrophysics (astro-ph.SR)

• We report on subarcsecond observations of complex organic molecules (COMs) in the high-mass protostar IRAS20126+4104 with the Plateau de Bure Interferometer in its most extended configurations. In addition to the simple molecules SO, HNCO and H2-13CO, we detect emission from CH3CN, CH3OH, HCOOH, HCOOCH3, CH3OCH3, CH3CH2CN, CH3COCH3, NH2CN, and (CH2OH)2. SO and HNCO present a X-shaped morphology consistent with tracing the outflow cavity walls. Most of the COMs have their peak emission at the putative position of the protostar, but also show an extension towards the south(east), coinciding with an H2 knot from the jet at about 800-1000 au from the protostar. This is especially clear in the case of H2-13CO and CH3OCH3. We fitted the spectra at representative positions for the disc and the outflow, and found that the abundances of most COMs are comparable at both positions, suggesting that COMs are enhanced in shocks as a result of the passage of the outflow. By coupling a parametric shock model to a large gas-grain chemical network including COMs, we find that the observed COMs should survive in the gas phase for about 2000 yr, comparable to the shock lifetime estimated from the water masers at the outflow position. Overall, our data indicate that COMs in IRAS20126+4104 may arise not only from the disc, but also from dense and hot regions associated with the outflow.
• Asteroseismic parameters allow us to measure the basic stellar properties of field giants observed far across the Galaxy. Most of such determinations are, up to now, based on simple scaling relations involving the large frequency separation, ∆\nu, and the frequency of maximum power, \nu$_{max}$. In this work, we implement ∆\nu and the period spacing, ∆P, computed along detailed grids of stellar evolutionary tracks, into stellar isochrones and hence in a Bayesian method of parameter estimation. Tests with synthetic data reveal that masses and ages can be determined with typical precision of 5 and 19 per cent, respectively, provided precise seismic parameters are available. Adding independent information on the stellar luminosity, these values can decrease down to 3 and 10 per cent respectively. The application of these methods to NGC 6819 giants produces a mean age in agreement with those derived from isochrone fitting, and no evidence of systematic differences between RGB and RC stars. The age dispersion of NGC 6819 stars, however, is larger than expected, with at least part of the spread ascribable to stars that underwent mass-transfer events.
• The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model \sc fastwind with the genetic fitting algorithm \sc pikaia to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and Köhler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]
• Grid-based hydrodynamics simulations of circumstellar disks are often performed in the curvilinear coordinate system, in which the center of the computational domain coincides with the motionless star. However, the center of mass may be shifted from the star due to the presence of any non-axisymmetric mass distribution. As a result, the system exerts a gravity force on the star, causing the star to move in response, which can affect the evolution of the circumstellar disk. We aim at studying the effects of stellar motion on the evolution of protostellar and protoplanetary disks. In protostellar disks, a non-axisymmetric distribution of matter in the form of spiral arms or massive clumps can form due to gravitational instability. Protoplanetary disks can also feature non-axisymmetric structures caused by a high-mass planet or a large-scale vortex. We use 2D grid-based hydrodynamic simulations to explore the effect of stellar motion. We adopt a non-inertial polar coordinate system centered on the star, in which the stellar motion is taken into account by calculating the indirect potential caused by the non-axisymmetric disk, a high-mass planet, or a large-scale vortex. We found that the stellar motion has a moderate effect on the evolution history in protostellar disks, reducing somewhat the disk size and mass, while having a profound effect on the collapsing envelope, changing its inner shape from an initially axisymmetric to a non-axisymmetric configuration. Protoplanetary disk simulations show that the stellar motion slightly reduces the width of the gap opened by a high-mass planet, decreases the planet migration rate, and strengthens the large-scale vortices formed at the viscosity transition. We conclude that the inclusion of the indirect potential is recommended in grid-based hydrodynamics simulations of circumstellar disks which use the curvilinear coordinate system.
• Photometric observations in Sloan g' and i' bands of four W UMa binaries, V796 Cep, V797 Cep, CSS J015341.9+381641 and NSVS 3853195, are presented. Our observations showed that CSS J015404.1+382805 and NSVS 3853195 are the same star. We determined the initial epochs T_0 of all targets and improved the period of NSVS 3853195. The light curve solutions of our data revealed that the components of each target are almost the same in mass, temperature, radius and luminosity. The stellar components are of G and K spectral types and undergo partial eclipses. All systems have barely-overcontact configurations and belong to H subtype W UMa binaries. We established that the relation between the luminosity ratio l_2/l_1 and mass ratio q of our targets is approximately l_2/l_1 = q^1.5.
• [Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textitHerschel HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\u2009pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,\sc i] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,\sc i] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,\sc i] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources.
• In order to provide a better basis for the study of mechanisms of nucleosynthesis of the light elements beyond hydrogen and helium in the oldest stars, the abundances of C, O, Mg, Si, P, S, K, and Ca have been derived from UV-HST and visible-ESO high resolution spectra in the old, very metal-poor star HD 84937, at a metallicity that is 1/200 that of the Sun's. For this halo main-sequence turnoff star, the abundance determination of P and S are the first published determinations. The LTE profiles of the lines were fitted to the observed spectra. Wherever possible, we corrected the derived abundances for non-LTE effects. Three-dimensional (3D) CO5BOLD model atmospheres have been used to determine the abundances of C and O from molecular CH and OH bands. The abundances of these light elements in HD 84937 are found to agree well with the abundances in classical metal-poor stars. Our HD 84937 carbon abundance determination points toward a solar (or mildly enhanced) value of [C/Fe]. The modest overabundance of the alpha elements of even atomic number Z, typical of halo turnoff stars, is confirmed in this example. The odd-Z element P is found to be somewhat deficient in HD 84937, at [P/Fe]=-0.32, which is again consistent with the handful of existing determinations for turnoff stars of such low metallicity. We show that the abundance of oxygen, deduced from the OH band from 3D computations, is not compatible with the abundance deduced from the red oxygen triplet. This incompatibility is explained by the existence of a chromosphere heating the shallow layers of the atmosphere where the OH band, in 3D computations, is mainly formed. The abundance ratios are compared to the predictions of models of galactic nucleosynthesis and evolution
• We investigate numerically Alfvén waves propagating along an axisymmetric and non-isothermal solar flux tube embedded in the solar atmosphere. The tube magnetic field is current-free and diverges with height, and the waves are excited by a periodic driver along the tube magnetic field lines. The main results are that the two wave variables, the velocity and magnetic field perturbations in the azimuthal direction, behave differently as a result of gradients of physical parameters along the tube. To explain these differences in the wave behavior, the time evolution of the wave variables and the resulting cutoff period for each wave variable are calculated, and used to determine regions in the solar chromosphere where strong wave reflection may occur.
• XMM-Newton has deeply changed our picture of X-ray emission of hot, massive stars. High-resolution X-ray spectroscopy as well as monitoring of these objects helped us gain a deeper insight into the physics of single massive stars with or without magnetic fields, as well as of massive binary systems, where the stellar winds of both stars interact. These observations also revealed a number of previously unexpected features that challenge our understanding of the dynamics of the stellar winds of massive stars. I briefly summarize the results obtained over the past 15 years and highlight the perspectives for the next decade. It is anticipated that coordinated (X-ray and optical or UV) monitoring and time-critical observations of either single or binary massive stars will become the most important topics in this field over the coming years. Synergies with existing or forthcoming X-ray observatories (NuSTAR, Swift, eROSITA) will also play a major role and will further enhance the importance of XMM-Newton in our quest for understanding the physics of hot, massive stars.
• A new collection of photodissociation and photoionisation cross sections for 102 atoms and molecules of astrochemical interest has been assembled, along with a brief review of the basic processes involved. These have been used to calculate dissociation and ionisation rates, with uncertainties, in a standard ultraviolet interstellar radiation field (ISRF) and wavelength-dependent radiation fields. The new ISRF rates generally agree within 30% with our previous compilations, with a few notable exceptions. The reduction of rates in shielded regions was calculated as a function of dust, molecular and atomic hydrogen, atomic C, and self-shielding column densities. The relative importance of shielding types depends on the species in question and the dust optical properties. The new data are publicly available from the Leiden photodissociation and ionisation database. Sensitivity of rates to variation of temperature and isotope, and cross section uncertainties, are tested. Tests were conducted with an interstellar-cloud chemical model, and find general agreement (within a factor of two) with the previous iteration of the Leiden database for the ISRF, and order-of-magnitude variations assuming various kinds of stellar radiation. The newly parameterised dust-shielding factors makes a factor-of-two difference to many atomic and molecular abundances relative to parameters currently in the UDfA and KIDA astrochemical reaction databases. The newly-calculated cosmic-ray induced photodissociation and ionisation rates differ from current standard values up to a factor of 5. Under high temperature and cosmic-ray-flux conditions the new rates alter the equilibrium abundances of abundant dark cloud abundances by up to a factor of two. The partial cross sections for H2O and NH3 photodissociation forming OH, O, NH2 and NH are also evaluated and lead to radiation-field-dependent branching ratios.
• We present here the application of a model for a mass circulation mechanism in between the H-burning shell and the base of the convective envelope of low mass AGB stars, aimed at studying the isotopic composition of those presolar grains showing the most extreme levels of 18O depletion and high concentration of 26Mg from the decay of 26Al. The mixing scheme we present is based on a previously suggested magnetic-buoyancy process, already shown to account adequately for the formation of the main neutron source for slow neutron captures in AGB stars. We find that this scenario is also capable of reproducing for the first time the extreme values of the 17O/16O, the 18O/16O, and the 26Al/27Al isotopic ratios found in the mentioned oxide grains, including the highest amounts of 26Al there measured.
• This is a manual for the MOSiC package. MOSiC is a collection of IDL programs for profile analysis and Gaussian fitting of the Mg II h/k lines along with Gaussian fitting of the C II 133.5 nm line pair, the O I 135.6, the Cl I 135.2, the Si IV 139.7 and 140.3 and the O IV 140.0 nm lines observed with the IRIS near UV and far UV spectrograph. It was tested by analyzing over a hundred different IRIS data sets (quiet Sun, sunspot, ...). It works for off limb data, although it is still experimental. MOSiC analyzes different spectral lines separately and returns line intensity, width, and velocity for each line. A few sample profiles and maps are included in this manual.
• We investigate how star formation efficiency can be significantly decreased by the removal of a molecular cloud's envelope by feedback from an external source. Feedback from star formation has difficulties halting the process in dense gas but can easily remove the less dense and warmer envelopes where star formation does not occur. However, the envelopes can play an important role keeping their host clouds bound by deepening the gravitational potential and providing a constraining pressure boundary. We use numerical simulations to show that removal of the cloud envelopes results in all cases in a fall in the star formation efficiency (SFE). At 1.38 free-fall times our 4 pc cloud simulation experienced a drop in the SFE from 16 to six percent, while our 5 pc cloud fell from 27 to 16 per cent. At the same time, our 3 pc cloud (the least bound) fell from an SFE of 5.67 per cent to zero when the envelope was lost. The star formation efficiency per free-fall time varied from zero to $\approx$ 0.25 according to $\alpha$, defined to be the ratio of the kinetic plus thermal to gravitational energy, and irrespective of the absolute star forming mass available. Furthermore the fall in SFE associated with the loss of the envelope is found to even occur at later times. We conclude that the SFE will always fall should a star forming cloud lose its envelope due to stellar feedback, with less bound clouds suffering the greatest decrease.
• We analyse archival CGRO-BATSE X-ray flux and spin frequency measurements of GX 1+4 over a time span of 3000 days. We systematically search for time dependent variations of torque luminosity correlation. Our preliminary results indicate that the correlation shifts from being positive to negative on time scales of few 100 days.
• The characteristics of the gravitational collapse of a supernova and the fluxes of active and sterile neutrinos produced during the formation of its protoneutron core have been calculated numerically. The relative yields of active and sterile neutrinos in core matter with different degrees of neutronization have been calculated for various input parameters and various initial conditions. A significant increase in the fraction of sterile neutrinos produced in superdense core matter at the resonant degree of neutronization has been confirmed. The contributions of sterile neutrinos to the collapse dynamics and the total flux of neutrinos produced during collapse have been shown to be relatively small. The total luminosity of sterile neutrinos is considerably lower than the luminosity of electron neutrinos, but their spectrum is considerably harder at high energies.
• Power spectra of segmentation-cell length (a dominant length scale of EUV emission in the transition region) from full-disk HeII extreme ultraviolet (EUV) images observed by the Extreme ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory (SOHO) and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) during periods of quiet Sun conditions for a time interval from 1996 to 2015 were analyzed. The spatial power as a function of the spatial frequency from about 0.04 to 0.27 (EIT) or up to 0.48 (AIA) 1/Mm depends on the distribution of the observed segmentation-cell dimensions, -- a structure of the solar EUV network. The temporal variations of the spatial power reported by Didkovsky and Gurman (Solar Phys., 289, 153) were suggested as decreases at the mid-spatial frequencies for the compared spectra when the power curves at the highest spatial frequencies of 0.5 1/pix were adjusted to match each other. That approach has been extended in this work to compare spectral ratios at high spatial frequencies expressed in the solar spatial frequency units of 1/Mm. Based on this approach we represent these ratio changes as a long-term network transformation which may be interpreted as a continuous dissipation of mid-size network structures to the smaller-size structures in the transition region. In contrast to expected cycling of the segmentation-cell dimension structures and associated spatial power in the spectra with the solar cycle, the spectra demonstrate a significant and steady change of the EUV network. The temporal trend across these structural spectra is not critically sensitive to any long-term instrumental changes, e.g. degradation of sensitivity, but to the change of the segmentation-cell dimensions of the EUV network structure.