Radial Velocity Jitter in Stars from the California and Carnegie Planet Search at Keck Observatory

I present an empirical model for predicting a star's radial velocity jitter from its B-V color, activity level, and absolute magnitude. This model is based on observations of 450 well- observed stars from Keck Observatory for the California and Carnegie Planet Search Program. The model includes noise from both astrophysical sources and systematic errors, and describes jitter as generally increasing with a star's activity and height above the main sequence.Comment: 16 pages, 7 figures, PASP in pres

A spiral structure in the disk of EX Draconis on the rise to outburst maximum

We report on the R-band eclipse mapping analysis of high-speed photometry of the dwarf nova EX Dra on the rise to the maximum of the November 1995 outburst. The eclipse map shows a one-armed spiral structure of ~180 degrees in azimuth, extending in radius from R ~0.2 to 0.43 R_{L1} (where R_{L1} is the distance from the disk center to the inner Lagrangian point), that contributes about 22 per cent of the total flux of the eclipse map. The spiral structure is stationary in a reference frame co-rotating with the binary and is stable for a timescale of at least 5 binary orbits. The comparison of the eclipse maps on the rise and in quiescence suggests that the outbursts of EX Dra may be driven by episodes of enhanced mass-transfer from the secondary star. Possible explanations for the nature of the spiral structure are discussed.Comment: To appear in the Astrophysical Journal Letters; 8 pages, 2 figures; coded with AAS latex styl

Ca II H and K Chromospheric Emission Lines in Late K and M Dwarfs

We have measured the profiles of the Ca II H and K chromospheric emission lines in 147 main sequence stars of spectral type M5-K7 (0.30-0.55 solar masses) using multiple high resolution spectra obtained during six years with the HIRES spectrometer on the Keck 1 telescope. Remarkably, the average FWHM, equivalent widths, and line luminosities of Ca II H and K increase by a factor of 3 with increasing stellar mass over this small range of stellar masses. We fit the H and K lines with a double Gaussian model to represent both the chromospheric emission and the non-LTE central absorption. Most of the sample stars display a central absorption that is typically redshifted by ~0.1 km/s relative to the emission, but the nature of this velocity gradient remains unknown. The FWHM of the H and K lines increase with stellar luminosity, reminiscent of the Wilson-Bappu effect in FGK-type stars. Both the equivalent widths and FWHM exhibit modest temporal variability in individual stars. At a given value of M_v, stars exhibit a spread in both the equivalent width and FWHM of Ca II H and K, due both to a spread in fundamental stellar parameters including rotation rate, age, and possibly metallicity, and to the spread in stellar mass at a given M_v. The K line is consistently wider than the H line, as expected, and its central absorption is more redshifted, indicating that the H and K lines form at slightly different heights in the chromosphere where the velocities are slightly different. The equivalent width of H-alpha correlates with Ca II H and K only for stars having Ca II equivalent widths above ~2 angstroms, suggesting the existence of a magnetic threshold above which the lower and upper chromospheres become thermally coupled.Comment: 40 pages including 12 figures and 17 pages of tables, accepted for publication in PAS

Spectropolarimetry of the nova-like variable V1315 Aql

We present spectropolarimetric observations of the eclipsing nova-like variable V1315 Aql, obtained with the aim of determining whether the single-peaked uneclipsed lines observed in this and related nova-likes are the result of disc emission scattered into the line of sight by the wind. The data show linear polarization with a mean value of 0.11+-0.02%. There are no significant differences between the continuum and line polarizations and no significant variations with wavelength or binary phase. We argue that the measured polarization may be attributed to scattering in the interstellar medium and hence conclude that there is no evidence of polarization intrinsic to V1315 Aql. We discuss alternative models which promise to resolve the controversy surrounding these objects

HST/FOS Eclipse Observations of the Nova-like Cataclysmic Variable UX Ursae Majoris

[abridged abstract] We present and analyze Hubble Space Telescope observations of the eclipsing nova-like cataclysmic variable UX UMa obtained with the Faint Object Spectrograph. Two eclipses each were observed with the G160L grating (covering the ultraviolet waveband) in August of 1994 and with the PRISM (covering the near-ultraviolet to near-infrared) in November of the same year. The system was 50% brighter in November than in August, which, if due to a change in the accretion rate, indicates a fairly substantial increase in Mdot_acc by >~ 50%. Model disk spectra constructed as ensembles of stellar atmospheres provide poor descriptions of the observed post-eclipse spectra, despite the fact that UX UMa's light should be dominated by the disk at this time. Suitably scaled single temperature model stellar atmospheres with T_eff = 12,500-14,500 K actually provide a better match to both the ultraviolet and optical post-eclipse spectra. Evidently, great care must be taken in attempts to derive accretion rates from comparisons of disk models to observations. One way to reconcile disk models with the observed post-eclipse spectra is to postulate the presence of a significant amount of optically thin material in the system. Such an optically thin component might be associated with the transition region (``chromosphere'') between the disk photosphere and the fast wind from the system, whose presence has been suggested by Knigge & Drew (1997).Comment: 35 pages, including 12 figures; to appear in the ApJ (Vol. 499

Non-reversible Gaussian processes for identifying latent dynamical structure in neural data

A common goal in the analysis of neural data is to compress large population recordings into sets of interpretable, low-dimensional latent trajectories. This problem can be approached using Gaussian process (GP)-based methods which provide uncertainty quantification and principled model selection. However, standard GP priors do not distinguish between underlying dynamical processes and other forms of temporal autocorrelation. Here, we propose a new family of “dynamical” priors over trajectories, in the form of GP covariance functions that express a property shared by most dynamical systems: temporal non-reversibility. Non-reversibility is a universal signature of autonomous dynamical systems whose state trajectories follow consistent flow fields, such that any observed trajectory could not occur in reverse. Our new multi-output GP kernels can be used as drop-in replacements for standard kernels in multivariate regression, but also in latent variable models such as Gaussian process factor analysis (GPFA). We therefore introduce GPFADS (Gaussian Process Factor Analysis with Dynamical Structure), which models single-trial neural population activity using low-dimensional, non-reversible latent processes. Unlike previously proposed non-reversible multi-output kernels, ours admits a Kronecker factorization enabling fast and memory-efficient learning and inference. We apply GPFADS to synthetic data and show that it correctly recovers ground truth phase portraits. GPFADS also provides a probabilistic generalization of jPCA, a method originally developed for identifying latent rotational dynamics in neural data. When applied to monkey M1 neural recordings, GPFADS discovers latent trajectories with strong dynamical structure in the form of rotations