EXPRES. III. Revealing the Stellar Activity Radial Velocity Signature of ϵ Eridani with Photometry and Interferometry

The distortions of absorption line profiles caused by photospheric brightness variations on the surfaces of cool, main-sequence stars can mimic or overwhelm radial velocity (RV) shifts due to the presence of exoplanets. The latest generation of precision RV spectrographs aims to detect velocity amplitudes ≲ 10 cm s−1, but requires mitigation of stellar signals. Statistical techniques are being developed to differentiate between Keplerian and activity-related velocity perturbations. Two important challenges, however, are the interpretability of the stellar activity component as RV models become more sophisticated, and ensuring the lowest-amplitude Keplerian signatures are not inadvertently accounted for in flexible models of stellar activity. For the K2V exoplanet host ϵ Eridani, we separately used ground-based photometry to constrain Gaussian processes for modeling RVs and TESS photometry with a light-curve inversion algorithm to reconstruct the stellar surface. From the reconstructions of TESS photometry, we produced an activity model that reduced the rms scatter in RVs obtained with EXPRES from 4.72 to 1.98 m s−1. We present a pilot study using the CHARA Array and MIRC-X beam combiner to directly image the starspots seen in the TESS photometry. With the limited phase coverage, our spot detections are marginal with current data but a future dedicated observing campaign should allow for imaging, as well as allow the stellar inclination and orientation with respect to the debris disk to be definitively determined. This work shows that stellar surface maps obtained with high-cadence, time-series photometric and interferometric data can provide the constraints needed to accurately reduce RV scatter.

Measuring the density structure of an accretion hot spot

Magnetospheric accretion models predict that matter from protoplanetary disks accretes onto stars via funnel flows, which follow stellar magnetic field lines and shock on the stellar surfaces1–3, leaving hot spots with density gradients4–6. Previous work has provided observational evidence of varying density in hot spots7, but these observations were not sensitive to the radial density distribution. Attempts have been made to measure this distribution using X-ray observations8–10; however, X-ray emission traces only a fraction of the hot spot11,12 and also coronal emission13,14. Here we report periodic ultraviolet and optical light curves of the accreting star GM Aurigae, which have a time lag of about one day between their peaks. The periodicity arises because the source of the ultraviolet and optical emission moves into and out of view as it rotates along with the star. The time lag indicates a difference in the spatial distribution of ultraviolet and optical brightness over the stellar surface. Within the framework of a magnetospheric accretion model, this finding indicates the presence of a radial density gradient in a hot spot on the stellar surface, because regions of the hot spot with different densities have different temperatures and therefore emit radiation at different wavelengths.

Tangled magnetic field model of QPOs

The highly tangled magnetic field of a magnetar supports shear waves similar to Alfvén waves in an ordered magnetic field. Here we explore if torsional modes excited in the stellar interior and magnetosphere can explain the quasi-periodic oscillations (QPOs) observed in the tail of the giant flare of SGR 1900+14. We solve the initial value problem for a tangled magnetic field that couples interior shear waves to relativistic Alfvén shear waves in the magnetosphere. Assuming stellar oscillations arise from the sudden release of magnetic energy, we obtain constraints on the energetics and geometry of the process. If the flare energy is deposited initially inside the star, the wave energy propagates relatively slowly to the magnetosphere which is at odds with the observed rise time of the radiative event of ≲ 10 ms. Nor can the flare energy be deposited entirely outside the star, as most of the energy reflects off the stellar surface, giving surface oscillations of insufficient magnitude to produce detectable modulations of magnetospheric currents. Energy deposition in a volume that straddles the stellar surface gives agreement with the observed rise time and excites a range of modes with substantial amplitude at observed QPO frequencies. In general, localized energy deposition excites a broad range of modes that encompasses the observed QPOs, though many more modes are excited than the number of observed QPOs. If the flare energy is deposited axisymmetrically, as is possible for a certain class of MHD instabilities, the number of modes that is excited is considerably reduced.

Inferring possible magnetic field strength of accreting inflows in EXor-type objects from scaled laboratory experiments

Aims. EXor-type objects are protostars that display powerful UV-optical outbursts caused by intermittent and powerful events of magnetospheric accretion. These objects are not yet well investigated and are quite difficult to characterize. Several parameters, such as plasma stream velocities, characteristic densities, and temperatures, can be retrieved from present observations. As of yet, however, there is no information about the magnetic field values and the exact underlying accretion scenario is also under discussion. Methods. We use laboratory plasmas, created by a high power laser impacting a solid target or by a plasma gun injector, and make these plasmas propagate perpendicularly to a strong external magnetic field. The propagating plasmas are found to be well scaled to the presently inferred parameters of EXor-type accretion event, thus allowing us to study the behaviour of such episodic accretion processes in scaled conditions. Results. We propose a scenario of additional matter accretion in the equatorial plane, which claims to explain the increased accretion rates of the EXor objects, supported by the experimental demonstration of effective plasma propagation across the magnetic field. In particular, our laboratory investigation allows us to determine that the field strength in the accretion stream of EXor objects, in a position intermediate between the truncation radius and the stellar surface, should be of the order of 100 G. This, in turn, suggests a field strength of a few kilogausses on the stellar surface, which is similar to values inferred from observations of classical T Tauri stars.

Magnetic torques on T Tauri stars: Accreting versus non-accreting systems

Context. Classical T Tauri stars (CTTs) magnetically interact with their surrounding disks, a process that is thought to regulate their rotational evolution. Aims. We compute torques acting on the stellar surface of CTTs that arise from different accreting (accretion funnels) and ejecting (stellar winds and magnetospheric ejections) flow components. Furthermore, we compare the magnetic braking due to stellar winds in two different systems: isolated (i.e., weak-line T Tauri and main-sequence) and accreting (i.e., classical T Tauri) stars. Methods. We use 2.5D magnetohydrodynamic, time-dependent, axisymmetric simulations that were computed with the PLUTO code. For both systems, the stellar wind is thermally driven. In the star-disk-interaction (SDI) simulations, the accretion disk is Keplerian, viscous, and resistive, and is modeled with an alpha prescription. Two series of simulations are presented, one for each system (i.e., isolated and accreting stars). Results. In classical T Tauri systems, the presence of magnetospheric ejections confines the stellar-wind expansion, resulting in an hourglass-shaped geometry of the outflow, and the formation of the accretion columns modifies the amount of open magnetic flux exploited by the stellar wind. These effects have a strong impact on the stellar-wind properties, and we show that the stellar-wind braking is more efficient in the SDI systems than in the isolated ones. We further derive torque scalings over a wide range of magnetic field strengths for each flow component in an SDI system (i.e., magnetospheric accretion and ejections, and stellar winds), which directly applies a torque on the stellar surface. Conclusions. In all the performed SDI simulations, the stellar wind extracts less than 2% of the mass accretion rate and the disk is truncated by up to 66% of the corotation radius. All simulations show a net spin-up torque. We conclude that in order to achieve a stellar-spin equilibrium, we need either more massive stellar winds or disks that are truncated closer to the corotation radius, which increases the torque efficiency of the magnetospheric ejections.

Can complex T Tauri star light curves be modelled with star-spots?

ABSTRACT ‘Transiting Exoplanet Survey Satellite’ photometry of 10 young stars with very complex periodic light curves is considered. Previous findings that these cannot be due to dark surface spots are re-evaluated by allowing arbitrarily shaped areas on the stellar surface to have sub-photospheric fluxes. This is done by approximating flux integrals by sums over surface elements. The unknown ratios of spot to photospheric fluxes are determined by lasso or ridge regression procedures. It is found that almost all light curves can be modelled very accurately in this way. The usual, if rarely stated, caveat applies – star-spot models presented in the paper are not unique.

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch (AGB) stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics, which amplifies the uncertainties on stellar parameters and distances. Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics (RHD) simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius using the average 3D intensity profile and, finally, compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 μm, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths (3.299 ± 0.005 mas and 3.053 ± 0.006 mas at 1.52 and 1.70 μm, respectively) and recovered the radius (R = 307 ± 41 and R = 284 ± 38 R⊙) using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

What determines the maximum stellar surface density of galaxies?

ABSTRACT Observationally, it has been reported that the densest stellar system in the Universe does not exceed a maximum stellar surface density, $\Sigma ^{\max }_{*}$ = $3\times 10^5\, {\rm M}_{\odot }\,{\rm pc}^{-2}$, throughout a wide physical scale ranging from star cluster to galaxy. This suggests that there exists a fundamental physics that regulates the star formation and stellar density. However, factors that determine this maximum limit are not clear. In this study, we show that $\Sigma ^{\max }_{*}$ of galaxies is not a constant as previous work reported, but actually it depends on the stellar mass. We select galaxy sample from the Sloan Digital Sky Survey Data Release 12 at z = 0.01–0.5. In contrast to a constant maximum predicted by theoretical models, $\Sigma ^{\max }_{*}$ strongly depends on stellar mass, especially for less massive galaxies with $\text{$\sim$}10^{10}\, {\rm M}_{\odot }$. We also found that a majority of high-Σ* galaxies show red colours and low star formation rates. These galaxies probably reach the $\Sigma ^{\max }_{*}$ as a consequence of the galaxy evolution from blue star forming to red quiescent by quenching star formation. One possible explanation of the stellar-mass dependence of $\Sigma ^{\max }_{*}$ is a mass-dependent efficiency of stellar feedback. The stellar feedback could be relatively more efficient in a shallower gravitational potential, which terminates star formation quickly before the stellar system reaches a high stellar density.