Solar and Stellar Astrophysics (astro-ph.SR)

  • PDF
    Four different approximate approaches used to model the stressing of coronal magnetic fields due to an imposed photospheric motion are compared with each other and the results from a full time-dependent magnetohydrodynamic (MHD) code. The assumptions used for each of the approximate methods are tested by considering large photospheric footpoint displacements. We consider a simple model problem, comparing the full nonlinear magnetohydrodynamic evolution, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time-dependent. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements but has significant errors when the plasma $\beta_0$ is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low $\beta_0$ plasmas if the viscosity is too large.
  • PDF
    Observations of protostars are often compared with synthetic observations of models in order to infer the underlying physical properties of the protostars. The majority of these models have a single protostar, attended by a disc and an envelope. However, observational and numerical evidence suggests that a large fraction of protostars form as multiple systems. This means that fitting models of single protostars to observations may be inappropriate. We produce synthetic observations of protostellar multiple systems undergoing realistic, non-continuous accretion. These systems consist of multiple protostars with episodic luminosities, embedded self-consistently in discs and envelopes. We model the gas dynamics of these systems using smoothed particle hydrodynamics and we generate synthetic observations by post-processing the snapshots using the \textscspamcart Monte Carlo radiative transfer code. We present simulation results of three model protostellar multiple systems. For each of these, we generate $4\times10^4$ synthetic spectra at different points in time and from different viewing angles. We propose a Bayesian method, using similar calculations to those presented here, but in greater numbers, to infer the physical properties of protostellar multiple systems from observations.
  • PDF
    We present near infrared high-precision photometry for eight transiting hot Jupiters observed during their predicted secondary eclipses. Our observations were carried out using the staring mode of the WIRCam instrument on the Canada-France-Hawaii Telescope (CFHT). We present the observing strategies and data reduction methods which delivered time series photometry with statistical photometric precisionas low as 0.11%. We performed a Bayesian analysis to model the eclipse parameters and systematics simultaneously. The measured planet-to-star flux ratios allowed us to constrain the thermal emission from the day side of these hot Jupiters, as we derived the planet brightness temperatures. Our results combined with previously observed eclipses reveal an excess in the brightness temperatures relative to the blackbody prediction for the equilibrium temperatures of the planets for a wide range of heat redistribution factors. We find a trend that this excess appears to be larger for planets with lower equilibrium temperatures. This may imply some additional sources of radiation, such as reflected light from the host star and/or thermal emission from residual internal heat from the formation of the planet.
  • PDF
    The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure.
  • PDF
    This review summarises results of recent magnetic and chemical abundance surface mapping studies of early-type stars. We discuss main trends uncovered by observational investigations and consider reliability of spectropolarimetric inversion techniques used to infer these results. A critical assessment of theoretical attempts to interpret empirical magnetic and chemical maps in the framework of, respectively, the fossil field and atomic diffusion theories is also presented. This confrontation of theory and observations demonstrates that 3D MHD models of fossil field relaxation are successful in matching the observed range of surface magnetic field geometries. At the same time, even the most recent time-dependent atomic diffusion calculations fail to reproduce diverse horizontal abundance distributions found in real magnetic hot stars.
  • PDF
    A series of recent theoretical atomic diffusion studies has address the challenging problem of predicting inhomogeneous vertical and horizontal chemical element distributions in the atmospheres of magnetic ApBp stars. Here we critically assess the most sophisticated of such diffusion models - based on a time-dependent treatment of the atomic diffusion in a magnetised stellar atmosphere - by direct comparison with observations as well by testing the widely used surface mapping tools with the spectral line profiles predicted by this theory. We show that the mean abundances of Fe and Cr are grossly underestimated by the time-dependent theoretical diffusion model, with discrepancies reaching a factor of 1000 for Cr. We also demonstrate that Doppler imaging inversion codes, based either on modelling of individual metal lines or line-averaged profiles simulated according to theoretical three-dimensional abundance distribution, are able to reconstruct correct horizontal chemical spot maps despite ignoring the vertical abundance variation. These numerical experiments justify a direct comparison of the empirical two-dimensional Doppler maps with theoretical diffusion calculations. This comparison is generally unfavourable for the current diffusion theory, as very few chemical elements are observed to form overabundance rings in the horizontal field regions as predicted by the theory and there are numerous examples of element accumulations in the vicinity of radial field zones, which cannot be explained by diffusion calculations.
  • PDF
    We present results from the investigation of 5-min umbral oscillations in a single-polarity sunspot of active region NOAA 12132. The spectra of TiO, H$\alpha$, and 304 Å are used for corresponding atmospheric heights from the photosphere to lower corona. Power spectrum analysis at the formation height of H$\alpha$ - 0.6 Å to H$\alpha$ center resulted in the detection of 5-min oscillation signals in intensity interpreted as running waves outside the umbral center, mostly with vertical magnetic field inclination $>15\deg$. A phase-speed filter is used to extract the running wave signals with speed $v_{ph}> 4$ km s$^{-1}$, from the time series of H$\alpha$ - 0.4 Å images, and found twenty-four 3-min umbral oscillatory events in a duration of one hour. Interestingly, the initial emergence of the 3-min umbral oscillatory events are noticed closer to or at umbral boundaries. These 3-min umbral oscillatory events are observed for the first time as propagating from a fraction of preceding Running Penumbral Waves (RPWs). These fractional wavefronts rapidly separates from RPWs and move towards umbral center, wherein they expand radially outwards suggesting the beginning of a new umbral oscillatory event. We found that most of these umbral oscillatory events develop further into RPWs. We speculate that the waveguides of running waves are twisted in spiral structures and hence the wavefronts are first seen at high latitudes of umbral boundaries and later at lower latitudes of the umbral center.
  • PDF
    This short review aims at clarifying upon the origins of the well-known Eddington-Barbier relationships. The later relate indeed both the emergent specific intensity and the flux to the stellar photospheric source function at certain specific optical depths. Hereafter we shall discuss about the very assumptions used for the original derivation of Barbier (1943). We also point to the fact that Milne (1921) had already formulated these two relations.
  • PDF
    Classical Cepheids are well-known and widely used distance indicators. Since distance and extinction are usually degenerate, it is important to develop suitable methods to robustly anchor the distance scale. Here, we introduce a near-infrared (near-IR) optimal distance method to determine both the extinction values of and distances to a large sample of 289 Galactic classical Cepheids. The overall uncertainty in the derived distances is less than 4.9%. We compare our newly determined distances to the Cepheids in our sample with previously published distances to the same Cepheids with Hubble Space Telescope parallax measurements and distances based on the IR surface brightness method, Wesenheit functions, and the main-sequence fitting method. The systematic deviations in the distances determined here with respect to those of previous publications is less than 1-2%. We hence constructed Galactic mid-IR period-luminosity (PL) relations for classical Cepheids in the four Wide-Field Infrared Survey Explorer (WISE) bands (W1, W2, W3, and W4) and the four Spitzer Space Telescope IRAC bands ([3.6], [4.5], [5.8] and [8.0]). Based on our sample of hundreds of Cepheids, the WISE PL relations have been determined for the first time; their dispersion is approximately 0.10 mag. Using the currently most complete sample, our Spitzer PL relations represent a significant improvement in accuracy, especially in the [3.6] band which has the smallest dispersion (0.066 mag). In addition, the average mid-IR extinction curve for Cepheids has been obtained: A_W1/A_Ks=0.560, A_W2/A_Ks=0.479, A_W3/A_Ks=0.507, A_W4/A_Ks=0.406, A_[3.6]/A_Ks=0.481, A_[4.5]/A_Ks=0.469, A_[5.8]/A_Ks=0.427, and A_[8.0]/A_Ks=0.427 mag.
  • PDF
    The Cygnus complex is one of the most powerful star forming regions at a close distance from the Sun (~1.4 kpc). Its richest OB association Cygnus OB2 is known to harbour many tens of O-type stars and hundreds of B-type stars, providing a large homogeneous population of OB stars that can be analysed. Many studies of its massive population have been developed in the last decades, although the total number of OB stars is still incomplete. We have obtained new blue intermediate-resolution spectra suitable for spectral classification of the 61 candidates presented by Comerón & Pasquali (2012) in Cygnus OB2 and surroundings. We thus performed a spectral classification of the sample using He I-II and metal lines rates, as well as the Marxist Ghost Buster (MGB) software for O-type stars and the IACOB standards catalogue for B-type stars. Out of the 61 candidates, we have classified 42 stars as new massive OB-type stars, earlier than B3, in Cygnus OB2 and surroundings, including 11 O-type stars. The other candidates are discarded as they display later spectral types inconsistent with membership in the association. However, the magnitude cutoff and dust extinction introduce an incompleteness. Many O and early B stars at B > 16 mag are still undiscovered in the region. Finally, we have studied the age and extinction distribution of our sample within the region, placing them in the Hertzsprung-Russell Diagram using both rotating and non-rotating models calculated by Ekström et al. (2012) and Brott et al. (2011) in order to assess age uncertainties. Massive star formation in Cygnus OB2 seems to have proceeded from lower to higher Galactic longitudes, regardless of the details of the models used. The correlation between age and Galactic longitude previously found in the region is now confirmed.
  • PDF
    Magnetic confinement of the winds of hot, massive stars has far-reaching consequences on timescales ranging from hours to Myr. Understanding the long-term effects of this interplay has already led to the identification of two new evolutionary pathways to form `heavy' stellar mass black holes and pair-instability supernova even at galactic metallicity. We are performing 1D stellar evolution model calculations that, for the first time, account for the surface effects and the time evolution of fossil magnetic fields. These models will be thoroughly confronted with observations and will potentially lead to a significant revision of the derived parameters of observed magnetic massive stars.
  • PDF
    We used Fourier transform techniques to determine macroturbulent velocity. The analysis is done with mictoturbulent velocity and rotation velocity as an unknown quantities. In order to distinguish the effects of rotation from macroturbulence effects in slowly rotating stars, primarily the main lobe of residual Fourier transforms of the observed lines, which were taken from the solar spectrum and the spectra of two other stars, was used. This case of Fourier analysis of spectral lines is the most complicated one. The end results were in a satisfactory agreement with the data obtained using different methods. We have obtained the average values of microturbulent, macroturbulent, and rotation velocities were 0.85, 2.22, and 1.75 km/s for the Sun as the star, 0.58, 1.73, and 0.78 km/s for HD 10700, and 1.16, 3.56, and 6.24 km/s for HD 1835. It was found that the macroturbulent velocity decreases with height in the atmosphere of the Sun and HD 1835. In the case of HD 10700, the macroturbulent velocity did not change with height, and the determined rotation velocity was two times lower than the one obtained using other methods. It was concluded that Fourier transform techniques are suitable for determining the velocities in atmospheres of solar-type stars with very slow rotation.
  • PDF
    Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.
  • PDF
    The age and evolutionary status of MWC349A, the unique emission line star with maser and laser radiation in hydrogen recombination lines, remain unknown because the spectrum of the star is veiled by bright emission from the ionized disk and wind. The major argument for this massive ($>$10 M_sun) star being evolved is its association with a close-by (2.4 arcsec) companion, MWC349B, whose B0III spectrum implies an age of a few Myrs. However, newly obtained high-resolution spectra of MWC349B reveal a difference $\approx$ 35 km/s in the radial velocities of the two stars, which makes their being gravitationally bound highly improbable. An estimate of the relative proper motion of the two stars seems to confirm this conclusion. This reopens the previously suggested possibility that MWC349A is a young massive star in a region of active star formation close to the Cyg OB2 association. MWC349B, which moves with a speed $\geq$ 35 km/s relative to Cyg OB2, may be a runaway star from this association.
  • PDF
    A prediction of the steady-state reconnection electric field in asymmetric reconnection is obtained by maximizing the reconnection rate as a function of the opening angle made by the upstream magnetic field on the weak magnetic field (magnetosheath) side. The prediction is within a factor of two of the widely examined asymmetric reconnection model [Cassak and Shay, Phys. Plasmas 14, 102114, 2007] in the collisionless limit, and they scale the same over a wide parameter regime. The previous model had the effective aspect ratio of the diffusion region as a free parameter, which simulations and observations suggest is on the order of 0.1, but the present model has no free parameters. In conjunction with the symmetric case [Liu et al., Phys. Rev. Lett. 118, 085101, 2017], this work further suggests that this nearly universal number 0.1, essentially the normalized fast reconnection rate, is a geometrical factor arising from maximizing the reconnection rate within magnetohydrodynamic (MHD)-scale constraints.
  • PDF
    We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through one- and two-dimensional test problems.