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

• Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructing a thermodynamic MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption's source region and initiate the eruption by boundary flows. More than 10^33 ergs of magnetic energy are released in the eruption within a few minutes, driving a flare, an EUV wave, and a coronal mass ejection (CME) that travels in the outer corona at about 1500 km/s, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic-equilibrium state. We find that the flux-rope center is very similar in character to the observed magnetic cloud, but arrives about 8.5 hours later and about 15 degrees too far to the North, with field strengths that are too weak by a factor of about 1.6. The front of the flux rope is highly distorted, exhibiting localized magnetic-field concentrations as it passes 1 AU. We discuss these properties with regard to the development of space-weather predictions based on MHD simulations of solar eruptions.
• Globular clusters are tracers of the history of star formation and chemical enrichment in the early Galaxy. Their abundance pattern can help understanding their chemical enrichment processes. In particular, the iron-peak elements have been relatively little studied so far in the Galactic bulge. The abundances of iron-peak elements are derived for the sample clusters, and compared with bulge field, and thick disk stars. We derive abundances of the iron-peak elements Sc, V, Mn, Cu, and Zn in individual stars of five bulge globular clusters (NGC 6528, NGC 6553, NGC 6522, NGC 6558, HP1), and of the reference thick disk/inner halo cluster 47 Tucanae (NGC 104). High resolution spectra were obtained with the UVES spectrograph at the Very Large Telescope over the years. The sample globular clusters studied span metallicities in the range -1.2<Fe/H]<0.0. V and Sc appear to vary in lockstep with Fe, indicating that they are produced in the same supernovae as Fe. We find that Mn is deficient in metal-poor stars, confirming that it is underproduced in massive stars; Mn-over-Fe steadily increases at the higher metallicities due to a metallicity-dependent enrichment by supernovae of type Ia. Cu behaves as a secondary element, indicating its production in a weak-s process in massive stars. Zn has an alpha-like behaviour at low metallicities, which can be explained in terms of nucleosynthesis in hypernovae. At the metal-rich end, Zn decreases with increasing metallicity, similarly to the alpha-elements.
• All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citepNaze and spectroscopic \citepWalborn variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing for the synthesis of photometric observables based on the Analytic Dynamical Magnetosphere (ADM) model of \citetOwocki. We apply our model to OGLE photometry in order to constrain their magnetic geometries and surface dipole strengths. We predict that the field strengths for some of these candidate extra-Galactic magnetic stars may be within the detection limits of the FORS2 instrument.
• We present the identification of WISE~J135501.90$-$825838.9 as a spectral binary system with a slight possibility of planetary-mass components in the $130-200$\u2009Myr AB~Doradus moving group. Peculiarities in the near-infrared spectrum of this source suggest it to be a blended-light binary with L6.0$\pm$1.0 and T3.0$\pm$1.8 or L7.0$\pm$0.6 and T7.5$\pm$0.4 components. Its proper motion and radial velocity as a combined-light source yield a high membership probability for AB~Doradus. While the young L6+T3 case is underluminous in a color-magnitude diagram at the AB~Doradus kinematic distance, the young L7+T7.5 case could be viable. Gravity-sensitive indicators are more consistent with a field-age binary. If confirmed as a young object, member of AB~Doradus, we estimate masses of 11$\pm$1\,\mjup and 9$\pm$1\,\mjup with both component masses below the Deuterium burning mass limit. Otherwise, we find masses of $72^{+4}_{-5}$ and $61^{+6}_{-8}$ for the field L6+T3 case and $70^{+2}_{-4}$ and $42^{+5}_{-6}$ for the field L7+T7.5 case. Our identification of WISE~J135501.90$-$825838.9 as a candidate young spectral binary introduces a new technique for detecting and characterizing planetary-mass companions to young brown dwarfs.
• Reduced fluid models including electron inertia and ion finite Larmor radius corrections are derived asymptotically, both from fluid basic equations and from a gyrofluid model. They apply to collisionless plasmas with small ion-to-electron equilibrium temperature ratio and low $\beta_e$, where $\beta_e$ indicates the ratio between the equilibrium electron pressure and the magnetic pressure exerted by a strong, constant and uniform magnetic guide field. The consistency between the fluid and gyrofluid approaches is ensured when choosing ion closure relations prescribed by the underlying ordering. A two-field reduction of the gyrofluid model valid for arbitrary equilibrium temperature ratio is also introduced, and is shown to have a noncanonical Hamiltonian structure. This model provides a convenient framework for studying kinetic Alfvén wave turbulence, from MHD to sub-$d_e$ scales (where $d_e$ holds for the electron skin depth). Magnetic energy spectra are phenomenologically determined within energy and generalized helicity cascades in the perpendicular spectral plane. Arguments based on absolute statistical equilibria are used to predict the direction of the transfers, pointing out that, within the sub-ion range associated with a $k_\perp^{-7/3}$ transverse magnetic spectrum, the generalized helicity could display an inverse cascade if injected at small scales, for example by reconnection processes.
• Atmospheric properties of F-,G- and K-type stars can be measured by spectral model fitting. These methods require data with good signal-to-noise ratio and reliable continuum normalisation. This is particularly challenging for the spectra we have obtained with the CORALIE échelle spectrograph for FGK stars with transiting M-dwarf companions. The spectra tend to have low signal-to-noise ratios, which makes it difficult to analyse them using existing methods. Our aim is to create a reliable automated spectral analysis routine to determine $T_{\rm eff}$, [Fe/H], $V \sin i$ from the CORALIE spectra of FGK stars. We use wavelet decomposition to distinguish between noise, continuum trends, and stellar spectral features in the CORALIE spectra. A subset of wavelet coefficients from the target spectrum are compared to those from a grid of models in a Bayesian framework to determine the posterior probability distributions of the atmospheric parameters. By testing our method using synthetic spectra we found that our method converges on the best fitting atmospheric parameters. We test the wavelet method on 20 FGK exoplanet host stars for which higher quality data have been independently analysed using equivalent width measurements. We find that we can determine $T_{\rm eff}$ to a precision of $85$K, [Fe/H] to a precision of 0.06 dex and $V \sin i$ to a precision of 1.35 kms$^{-1}$ for stars with $V \sin i$ $\geq$ 5 kms$^{-1}$. We find an offset in metallicity $\approx -$0.18 dex relative to the equivelent width fitting method. We can determine $\log g$ to a precision of $0.13$ dex but measurements are only reliable to confirm dwarf-like surface gravity. The wavelet method can be used to determine $T_{\rm eff}$, [Fe/H] and $V \sin i$ for FGK stars from échelle spectra. We find that our method is self consistent, and robust for spectra with SNR$\gtrapprox 40$.
• One of the objectives of the NASA LWS Institute on Nowcasting of Atmospheric Drag for LEO Spacecraft was to investigate whether and how to increase the accuracy of atmospheric drag models by improving the quality of the solar forcing inputs, namely Extreme Ultraviolet (EUV) irradiance information. In this focused review, we examine the status of and issues with EUV measurements and proxies, discuss recent promising developments, and suggest a number of ways to improve the reliability, availability, and forecast accuracy of EUV measurements in the next solar cycle.
• We present analysis of medium resolution optical spectra and long term V band photometry of four cool stars, BD+13 5000, BD+11 3024, TYC 3557-919-1 and TYC 5163-1764-1. Our spectroscopic analysis reveals that the stars are giant or sub-giant from K0 or K1 spectral type, and all of them exhibit emission features in their Ca ii H& K lines. These features appear to be modulated with the rotation of the stars. Except BD+11 3024, we observe that the radial velocities of the target stars are not stable, which suggests that each of them might be a member of a binary system. Global analysis of photometric data indicates clear cyclic variation for BD+13 5000 and TYC 5163-1764-1 with a period of 8.0$\pm$0.3 and 5.04$\pm$0.04 year, respectively. Besides that, we observe a dramatic increase ($\sim$0\fm7) in the mean brightness of BD+11 3024, accompanied with a 2.87$\pm$0.12 cyclic variation, embedded into the global brightening trend, which indicates possible multiple cycles on this star.
• Context. Despite the fact that VW Cephei is one of the well-studied contact binaries in the literature, there is no fully consistent model available that can explain every observed property of this system. Aims. Our motivation is to obtain new spectra along with photometric measurements, to analyze what kind of changes may have happened in the system in the past two decades, and to propose new ideas for explaining them. Methods. For the period analysis we determined 10 new times of minima from our light curves, and constructed a new O$-$C diagram of the system. Radial velocities of the components were determined using the cross-correlation technique. The light curves and radial velocities were modelled simultaneously with the PHOEBE code. All observed spectra were compared to synthetic spectra and equivalent widths of the H$\alpha$ line were measured on their differences. Results. We have re-determined the physical parameters of the system according to our new light curve and spectral models. We confirm that the primary component is more active than the secondary, and there is a correlation between spottedness and the chromospheric activity. We propose that flip-flop phenomenon occurring on the primary component could be a possible explanation of the observed nature of the activity. To explain the period variation of VW Cep, we test two previously suggested scenarios: presence of a fourth body in the system, and the Applegate-mechanism caused by periodic magnetic activity. We conclude that although none of these mechanisms can be ruled out entirely, the available data suggest that mass transfer with a slowly decreasing rate gives the most likely explanation for the period variation of VW Cep.
• LoTr 5 is a planetary nebula with an unusual long-period binary central star. As far as we know, the pair consists of a rapidly rotating G-type star and a hot star, which is responsible for the ionization of the nebula. The rotation period of the G-type star is 5.95 days and the orbital period of the binary is now known to be $\sim$2700 days, one of the longest in central star of planetary nebulae. The spectrum of the G central star shows a complex H$\alpha$ double-peaked profile which varies with very short time scales, also reported in other central stars of planetary nebulae and whose origin is still unknown. We present new radial velocity observations of the central star which allow us to confirm the orbital period for the long-period binary and discuss the possibility of a third component in the system at $\sim$129 days to the G star. This is complemented with the analysis of archival light curves from SuperWASP, ASAS and OMC. From the spectral fitting of the G-type star, we obtain a effective temperature of $T_{\rm eff}$ = 5410$\pm$250 K and surface gravity of $\log g$ = 2.7$\pm$0.5, consistent with both giant and subgiant stars. We also present a detailed analysis of the H$\alpha$ double-peaked profile and conclude that it does not present correlation with the rotation period and that the presence of an accretion disk via Roche lobe overflow is unlikely.
• We report the discovery of two massive eccentric systems with BRITE data, $\tau$ Ori and $\tau$ Lib, showing heartbeat effects close to the periastron passage. $\tau$ Lib exhibits shallow eclipses that will soon vanish due to the apsidal motion in the system. In neither system, tidally excited oscillations were detected.
• We present high-resolution observations of a magnetic reconnection event in the solar atmosphere taken with the New Vacuum Solar Telescope, AIA and HMI. The reconnection event occurred between the threads of a twisted arch filament system (AFS) and coronal loops. Our observations reveal that the relaxation of the twisted AFS drives some of its threads to encounter the coronal loops, providing inflows of the reconnection. The reconnection is evidenced by flared X-shape features in the AIA images, a current-sheet-like feature apparently connecting post-reconnection loops in the \halpha$+$1 Å images, small-scale magnetic cancellation in the HMI magnetograms and flows with speeds of 40--80 km/s along the coronal loops. The post-reconnection coronal loops seen in AIA 94 Å passband appear to remain bright for a relatively long time, suggesting that they have been heated and/or filled up by dense plasmas previously stored in the AFS threads. Our observations suggest that the twisted magnetic system could release its free magnetic energy into the upper solar atmosphere through reconnection processes. While the plasma pressure in the reconnecting flux tubes are significantly different, the reconfiguration of field lines could result in transferring of mass among them and induce heating therein.
• We report on high-resolution imaging and spectral observations of eruptions of a spiral structure in the transition region, which were taken with the Interface Region Imaging Spectrometer (IRIS), the Atmospheric Imaging Assembly (AIA) and the Helioseismic and Magnetic Imager (HMI). The eruption coincided with the appearance of two series of jets, with velocities comparable to the Alfvén speeds in their footpoints. Several pieces of evidence of magnetic braiding in the eruption are revealed, including localized bright knots, multiple well-separated jet threads, transition region explosive events and the fact that all these three are falling into the same locations within the eruptive structures. Through analysis of the extrapolated three-dimensional magnetic field in the region, we found that the eruptive spiral structure corresponded well to locations of twisted magnetic flux tubes with varying curl values along their lengths. The eruption occurred where strong parallel currents, high squashing factors, and large twist numbers were obtained. The electron number density of the eruptive structure is found to be $\sim3\times10^{12}$ cm$^{-3}$, indicating that significant amount of mass could be pumped into the corona by the jets. Following the eruption, the extrapolations revealed a set of seemingly relaxed loops, which were visible in the AIA 94 Å channel indicating temperatures of around 6.3 MK. With these observations, we suggest that magnetic braiding could be part of the mechanisms explaining the formation of solar eruption and the mass and energy supplement to the corona.
• Element abundance ratios of magnesium to neon (Mg/Ne) and neon to oxygen (Ne/O) in the transition region of the quiet Sun have been derived by re-assessing previously published data from the Coronal Diagnostic Spectrometer on board the Solar and Heliospheric Observatory in the light of new atomic data. The quiet Sun Mg/Ne ratio is important for assessing the effect of magnetic activity on the mechanism of the first ionization potential (FIP) effect, while the Ne/O ratio can be used to infer the solar photospheric abundance of neon, which can not be measured directly. The average Mg/Ne ratio is found to be $0.52\pm 0.11$, which applies over the temperature region 0.2--0.7~MK, and is consistent with the earlier study. The Ne/O ratio is, however, about 40\%\ larger, taking the value $0.24\pm 0.05$ that applies to the temperature range 0.08--0.40~MK. The increase is mostly due to changes in ionization and recombination rates that affect the equilibrium ionization balance. If the Ne/O ratio is interpreted as reflecting the photospheric ratio, then the photospheric neon abundance is $8.08\pm 0.09$ or $8.15\pm 0.10$ (on a logarithmic scale for which hydrogen is 12), according to whether the oxygen abundances of M.~Asplund et al.\ or E.~Caffau et al.\ are used. The updated photospheric neon abundance implies a Mg/Ne FIP bias for the quiet Sun of $1.6\pm 0.6$.
• We report the first estimate of the He abundance of the population of RR Lyrae stars in the Galactic bulge. This is done by comparing the recent observational data with the latest models. We use the large samples of ab type RR Lyrae stars found by OGLE IV in the inner bulge and by the VVV survey in the outer bulge. We present the result from the new models computed by Marconi et al. (2017), showing that the minimum Period for fundamental RR Lyrae pulsators depends on the He content. By comparing these models with the observations in a Period versus effective temperature plane, we find that the bulk of the bulge ab type RR Lyrae are consistent with primordial He abundance Y=0.245, ruling out a significant He-enriched population. This work demonstrates that the He content of the bulge RR Lyrae is different from that of the bulk of the bulge population as traced by the red clump giants, that appear to be significantly more He-rich.
• This paper presents a detailed analysis of the light and radial velocity curves of the semi-detached eclipsing binary system OGLE-LMC-ECL-09937. The system is composed of a hot, massive and luminous primary star of a late-O spectral type, and a more evolved, but less massive and luminous secondary, implying an Algol-type system that underwent a mass transfer episode. We derive masses of 21.04 +/- 0.34 M_Sun and 7.61 +/- 0.09 M_Sun and radii of 9.93 +/- 0.06 R_Sun and 9.18 +/- 0.04 R_Sun, for the primary and the secondary component, respectively, which make it the most massive known Algol-type system with masses and radii of the components measured with <2% accuracy. Consequently, the parameters of OGLE-LMC-ECL-09937 provide an important contribution to the sparsely populated high-mass end of the stellar mass distribution, and an interesting object for stellar evolution studies, being a possible progenitor of a binary system composed of two neutron stars.
• We find that the element abundances in solar energetic particles (SEPs) and in the slow solar wind (SSW), relative to those in the photosphere, show different patterns as a function of the first ionization potential (FIP) of the elements. Generally, the SEP and SSW abundances reflect abundance samples of the solar corona, where low-FIP elements, ionized in the photosphere, are more efficiently conveyed upward to the corona than high-FIP elements that are initially neutral atoms. Abundances of the elements, especially C, P, and S show a crossover from low to high FIP at ~10 eV in the SEPs but ~14 eV for the solar wind, suggesting that cooler photospheric plasma, perhaps from cool sunspots beneath active regions, supply the coronal source material that eventually will be shock-accelerated as SEPs, while the source of the SSW may lie at the base of diverging open-field lines surrounding but outside active regions. Meanwhile, energetic particles accelerated from the solar wind itself by shock waves at corotating interaction regions (CIRs), generally beyond 1 AU, confirm the FIP pattern of the solar wind.
• We study the evolution of protoplanetary discs that would have been precursors of a Trappist-1 like system under the action of accretion and external photoevaporation in different radiation environments. Dust grains swiftly grow above the critical size below which they are entrained in the photoevaporative wind, so although gas is continually depleted, dust is resilient to photoevaporation after only a short time. This means that the ratio of the mass in solids (dust plus planetary) to the mass in gas rises steadily over time. Dust is still stripped early on, and the initial disc mass required to produce the observed $4\,M_{\oplus}$ of Trappist-1 planets is high. For example, assuming a Fatuzzo & Adams (2008) distribution of UV fields, typical initial disc masses have to be $>30\,$per cent the stellar (which are still Toomre $Q$ stable) for the majority of similar mass M dwarfs to be viable hosts of the Trappist-1 planets. Even in the case of the lowest UV environments observed, there is a strong loss of dust due to photoevaporation at early times from the weakly bound outer regions of the disc. This minimum level of dust loss is a factor two higher than that which would be lost by accretion onto the star during 10 Myr of evolution. Consequently even in these least irradiated environments, discs that are viable Trappist-1 precursors need to be initially massive ($>10\,$per cent of the stellar mass).
• We conducted radio-interferometric observations of the well known binary cataclysmic system AM Herculis. This particular system is formed by the magnetic white dwarf (primary) and the red dwarf (secondary), and is the prototype of so-called polars. Our observations were conducted with the European VLBI Network (EVN) in the e-EVN mode at 5 GHz. We obtained six astrometric measurements spanning one year, which make it possible to update the annual parallax for this system with the best precision to date ($\pi=11.29\pm0.08$ mas), equivalent to the distance of $88.6\pm0.6$ pc. The system was observed mostly in the quiescent phase (visual magnitude $m_v\sim15.3$), when the radio emission was at the level of about $300~\mu$Jy. Our analysis suggests that the radio flux of AM Herculis is modulated with the orbital motion. Such specific properties of the radiation can be explained using the emission mechanism similar to the scenario proposed for V471 Tauri and, in general, for RS CVn type stars. In this scenario the radio emission arises near the surface of the red dwarf, where the global magnetic field strength may reach a few kG. We argue that the quiescent radio emission distinguishes AM Herculis together with AR Ursae Majoris (the second known persistent radio polar) from other polars, as the systems with a magnetized secondary star.
• We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.
• The nature of the progenitors of Type Ia supernovae (SNe Ia) remains a mystery. Comparing theoretical rates and delay-time distributions of SNe Ia with those inferred observationally can constrain their progenitor models. In this work, taking thermohaline mixing into account in the helium-enriched main-sequence (HEMS) donor scenario, we address rates and delay times of SNe Ia in this channel by combining the results of self-consistent binary evolution calculations with population synthesis models. We find that the Galactic SN Ia rate from the HEMS donor scenario is around 0.6-1.2 per century, which is about 30% of the observed rate. Delay times of SNe Ia in this scenario cover a wide range of 0.1-1.0 Gyr. We also present the pre-explosion properties of companion stars in the HEMS donor scenario, which will be helpful for placing constraints on SN Ia progenitors through analyzing their pre-explosion images.
• We determined the orbits and stellar parameters of three eccentric eclipsing binaries by light curve solutions of their Kepler data. KIC 8111622 and KIC 10518735 undergo total eclipses while KIC 8196180 reveals partial eclipses. The target components are G and K stars, excluding the primary of KIC 8196180 which is early F star. KIC 8196180 reveals well-visible tidally-induced feature at periastron, i.e. it is an eclipsing heartbeat star. The characteristics of the observed periastron feature (shape, width and amplitude) confirm the theoretical predictions. There are additional out-of-eclipse variations of KIC 8196180 with the orbital period which may be explained by spot activity of synchronously rotating component. Besides worse visible periastron feature KIC 811162 exhibits small amplitude light variations whose period is around 2.3 times shorter than the orbital one. These oscillations were attributed to spot(s) on asynchronously rotating component.