X-Shooter spectroscopy of FU Tau A

We have analyzed a broad-band optical and near-infrared spectrum of FU Tau A, a presumed young brown dwarf in the Taurus star forming region that has intrigued both theorists and observers by its over-luminosity in the HR diagram with respect to standard pre-main sequence evolutionary models. The new data, obtained with the X-Shooter spectrograph at the Very Large Telescope, include an unprecedented wealth of information on stellar parameters and simultaneously observed accretion and outflow indicators for FU Tau A. We present the first measurements of gravity (log g = 3.5 +- 0.5), radial velocity (RV = 22.5 +- 2.9 km/s), rotational velocity (v sin(i) = 20 +- 5 km/s) and lithium equivalent width (W_Li = 430 +- 20 mAA) for FUTau A. From the rotational velocity and the published period we infer a disk inclination of i ~ 50^deg. The lithium content is much lower than theoretically expected for such a young very low mass object, adding another puzzling feature to this object's properties. We determine the mass accretion rate of FU Tau A from comparison of the luminosities of 24 emission lines to empirical calibrations from the literature and find a mean of log (dM/dt)_acc [M_sun/yr] = -9.9 +- 0.2. The accretion rate determined independently from modeling of the excess emission in the Balmer and Paschen continua is consistent with this value. The corresponding accretion luminosity is too small to make a significant contribution to the bolometric luminosity. The existence of an outflow in FU Tau A is demonstrated through the first detection of forbidden emission lines from which we obtain an estimate for the mass loss rate, log (dM/dt)_out [M_sun/yr] < -10.4. The mass outflow and inflow rates can be combined to yield (dM/dt)_out / (dM/dt)_acc ~ 0.3, a value that is in agreement with jet launching models.Comment: 12 pages, accepted for publication in A&

Thermal States as Convex Combinations of Matrix Product States

We study thermal states of strongly interacting quantum spin chains and prove that those can be represented in terms of convex combinations of matrix product states. Apart from revealing new features of the entanglement structure of Gibbs states our results provide a theoretical justification for the use of White's algorithm of minimally entangled typical thermal states. Furthermore, we shed new light on time dependent matrix product state algorithms which yield hydrodynamical descriptions of the underlying dynamics.Comment: v3: 10 pages, 2 figures, final published versio

Measurement of Orbital Decay in the Double Neutron Star Binary PSR B2127+11C

We report the direct measurement of orbital period decay in the double neutron star pulsar system PSR B2127+11C in the globular cluster M15 at the rate of $(-3.95 \pm 0.13) \times 10^{-12}$, consistent with the prediction of general relativity at the $\sim 3$ level. We find the pulsar mass to be $m_p = (1.358 \pm 0.010) M_\odot$ and the companion mass $m_c = (1.354 \pm 0.010) M_\odot$. We also report long-term pulse timing results for the pulsars PSR B2127+11A and PSR B2127+11B, including confirmation of the cluster proper motion.Comment: 12 pages, 4 figures, accepted for publication in ApJ

Observations of Mira stars with the IOTA/FLUOR interferometer and comparison with Mira star models

We present K'-band observations of five Mira stars with the IOTA interferometer. The interferograms were obtained with the FLUOR fiber optics beam combiner, which provides high-accuracy visibility measurements in spite of time-variable atmospheric conditions. For the M-type Miras X Oph, R Aql, RU Her, R Ser, and the C-type Mira V CrB we derived the uniform-disk diameters 11.7mas, 10.9mas, 8.4mas, 8.1mas, and 7.9mas (+/- 0.3mas), respectively. Simultaneous photometric observations yielded the bolometric fluxes. The derived angular Rosseland radii and the bolometric fluxes allowed the determination of effective temperatures. For instance, the effective temperature of R Aql was determined to be 2970 +/- 110 K. A linear Rosseland radius for R Aql of (250 +100/-60) Rsun was derived from the angular Rosseland radius of 5.5mas +/- 0.2mas and the HIPPARCOS parallax of 4.73mas +/- 1.19mas. The observations were compared with theoretical Mira star models of Bessel et al. (1996) and Hofmann et al. (1998). The effective temperatures of the M-type Miras and the linear radius of R Aql indicate fundamental mode pulsation.Comment: 12 pages, 4 postscript figure

VLTI/AMBER spectro-interferometric imaging of VX Sgr's inhomogenous outer atmosphere

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.Comment: 9 Pages, Accepted for publication on Astronomy & Astrophysics, two references update

Validation of nonlinear PCA

Linear principal component analysis (PCA) can be extended to a nonlinear PCA by using artificial neural networks. But the benefit of curved components requires a careful control of the model complexity. Moreover, standard techniques for model selection, including cross-validation and more generally the use of an independent test set, fail when applied to nonlinear PCA because of its inherent unsupervised characteristics. This paper presents a new approach for validating the complexity of nonlinear PCA models by using the error in missing data estimation as a criterion for model selection. It is motivated by the idea that only the model of optimal complexity is able to predict missing values with the highest accuracy. While standard test set validation usually favours over-fitted nonlinear PCA models, the proposed model validation approach correctly selects the optimal model complexity.Comment: 12 pages, 5 figure

A path towards understanding the rotation-activity relation of M dwarfs with K2 mission, X-ray and UV data

We study the relation between stellar rotation and magnetic activity fora sample of 134 bright, nearby M dwarfs observed in the Kepler Two-Wheel (K2) mission during campaigns C0-C4. The K2 light curves yield photometrically derived rotation periods for 97 stars (79 of which without previous period measurement), as well as various measures for activity related to cool spots and flares. We find a clear difference between fast and slow rotators with a dividing line at a period of ~10 d at which the activity level changes abruptly. All photometric diagnostics of activity (spot cycle amplitude, flare peak amplitude and residual variability after subtraction of spot and flare variations) display the same dichotomy, pointing to a quick transition between a high-activity mode for fast rotators and a low-activity mode for slow rotators. This unexplained behavior is reminiscent of a dynamo mode-change seen in numerical simulations that separates a dipolar from a multipolar regime. A substantial number of the fast rotators a revisual binaries. A tentative explanation is accelerated disc evolution in binaries leading to higher initial rotation rates on the main sequence and associated longer spin-down and activity lifetimes. We combine the K2 rotation periods with archival X-ray and UV data. X-ray, FUV and NUV detections are found for 26, 41, and 11 stars from our sample, respectively. Separating the fast from the slow rotators, we determine for the first time the X-ray saturation level separately for early- and for mid-M stars.Publisher PDFPeer reviewe

Scaling and Correlation Functions in a Model of a Two-dimensional Earthquake Fault

We study numerically a two-dimensional version of the Burrige-Knopoff model. We calculate spatial and temporal correlation functions and compare their behavior with the results found for the one-dimensional model. The Gutenberg-Richter law is only obtained for special choices of parameters in the relaxation function. We find that the distribution of the fractal dimension of the slip zone exhibits two well-defined peaks coeersponding to intermediate size and large events.Comment: 14 pages, 23 Postscript figure