Molecular Line Profile Fitting with Analytic Radiative Transfer Models

We present a study of analytic models of starless cores whose line profiles have ``infall asymmetry,'' or blue-skewed shapes indicative of contracting motions. We compare the ability of two types of analytical radiative transfer models to reproduce the line profiles and infall speeds of centrally condensed starless cores whose infall speeds are spatially constant and range between 0 and 0.2 km s-1. The model line profiles of HCO+ (J=1-0) and HCO+ (J=3-2) are produced by a self-consistent Monte Carlo radiative transfer code. The analytic models assume that the excitation temperature in the front of the cloud is either constant (``two-layer'' model) or increases inward as a linear function of optical depth (``hill'' model). Each analytic model is matched to the line profile by rapid least-squares fitting. The blue-asymmetric line profiles with two peaks, or with a blue shifted peak and a red shifted shoulder, can be well fit by the ``HILL5'' model (a five parameter version of the hill model), with an RMS error of 0.02 km s-1. A peak signal to noise ratio of at least 30 in the molecular line observations is required for performing these analytic radiative transfer fits to the line profiles.Comment: 48 pages, 20 figures, accepted for publication in Ap

Pair Correlations, Short Range Order and Dispersive Excitations in the Quasi-Kagome Quantum Magnet Volborthite

We present spatial and dynamic information on the s=1/2 distorted kagome antiferromagnet volborthite, Cu3V2O7(OD)2.2D2O, obtained by polarized and inelastic neutron scattering. The instantaneous structure factor, S(Q), is dominated by nearest neighbor pair correlations, with short range order at wave vectors Q1=0.65(3) {\AA}^-1 and Q2=1.15(5) {\AA}^-1 emerging below 5 K. The excitation spectrum, S(Q,{\omega}), reveals two steep branches dispersing from Q1 and Q2, and a flat mode at {\omega}=5.0(2) meV. The results allow us to identify the cross-over at T*=1 K in 51V NMR and specific heat measurements as the build-up of correlations at Q_1. We compare our data to theoretical models proposed for volborthite, and demonstrate that the excitation spectrum can be explained by spin-wave-like excitations with anisotropic exchange parameters, as also suggested by recent local density calculations.Comment: Rewritten article resubmitted to Phys. Rev. Lett. 021

Cantor type functions in non-integer bases

Cantor's ternary function is generalized to arbitrary base-change functions in non-integer bases. Some of them share the curious properties of Cantor's function, while others behave quite differently

Dzyaloshinsky-Moriya Anisotropy in the Spin-1/2 Kagom\'e Compound ZnCu3_{3}(OH)6_{6}Cl2_{2}

We report the determination of the Dzyaloshinsky-Moriya interaction, the dominant magnetic anisotropy term in the \kagome spin-1/2 compound {\herbert}. Based on the analysis of the high-temperature electron spin resonance (ESR) spectra, we find its main component $|D_z|=15(1)$ K to be perpendicular to the \kagome planes. Through the temperature dependent ESR line-width we observe a building up of nearest-neighbor spin-spin correlations below $\sim$150 K.Comment: 4 pages, 3 figures, minor modification

The problematically short superwind of OH/IR stars - Probing the outflow with the 69 {\mu}m spectral band of forsterite

Spectra of OH/IR stars show prominent spectral bands of crystalline olivine (Mg$_{(2-2x)}$Fe$_{(2x)}$SiO$_{4}$). To learn more about the timescale of the outflows of OH/IR stars, we study the spectral band of crystalline olivine at 69 {\mu}m. The 69 {\mu}m band is of interest because its width and peak wavelength position are sensitive to the grain temperature and to the exact composition of the crystalline olivine. With Herschel/PACS, we observed the 69 {\mu}m band in the outflow of 14 OH/IR stars. By comparing the crystalline olivine features of our sample with those of model spectra, we determined the size of the outflow and its crystalline olivine abundance. The temperature indicated by the observed 69 {\mu}m bands can only be reproduced by models with a geometrically compact superwind ($R_{\rm{SW}}\lesssim$ 2500 AU = 1400 R$_{*}$).This means that the superwind started less than 1200 years ago (assuming an outflow velocity of 10 km/s). The small amount of mass lost in one superwind and the high progenitor mass of the OH/IR stars introduce a mass loss and thus evolutionary problem for these objects, which has not yet been understood.Comment: Accepted by A&

Measurement and physical interpretation of the mean motion of turbulent density patterns detected by the BES system on MAST

The mean motion of turbulent patterns detected by a two-dimensional (2D) beam emission spectroscopy (BES) diagnostic on the Mega Amp Spherical Tokamak (MAST) is determined using a cross-correlation time delay (CCTD) method. Statistical reliability of the method is studied by means of synthetic data analysis. The experimental measurements on MAST indicate that the apparent mean poloidal motion of the turbulent density patterns in the lab frame arises because the longest correlation direction of the patterns (parallel to the local background magnetic fields) is not parallel to the direction of the fastest mean plasma flows (usually toroidal when strong neutral beam injection is present). The experimental measurements are consistent with the mean motion of plasma being toroidal. The sum of all other contributions (mean poloidal plasma flow, phase velocity of the density patterns in the plasma frame, non-linear effects, etc.) to the apparent mean poloidal velocity of the density patterns is found to be negligible. These results hold in all investigated L-mode, H-mode and internal transport barrier (ITB) discharges. The one exception is a high-poloidal-beta (the ratio of the plasma pressure to the poloidal magnetic field energy density) discharge, where a large magnetic island exists. In this case BES detects very little motion. This effect is currently theoretically unexplained.Comment: 28 pages, 15 figures, submitted to PPC

A Comprehensive X-ray Absorption Model for Atomic Oxygen

An analytical formula is developed to represent accurately the photoabsorption cross section of O I for all energies of interest in X-ray spectral modeling. In the vicinity of the Kedge, a Rydberg series expression is used to fit R-matrix results, including important orbital relaxation effects, that accurately predict the absorption oscillator strengths below threshold and merge consistently and continuously to the above-threshold cross section. Further minor adjustments are made to the threshold energies in order to reliably align the atomic Rydberg resonances after consideration of both experimental and observed line positions. At energies far below or above the K-edge region, the formulation is based on both outer- and inner-shell direct photoionization, including significant shake-up and shake-off processes that result in photoionization-excitation and double photoionization contributions to the total cross section. The ultimate purpose for developing a definitive model for oxygen absorption is to resolve standing discrepancies between the astronomically observed and laboratory measured line positions, and between the inferred atomic and molecular oxygen abundances in the interstellar medium from XSTAR and SPEX spectral models

Variational calculations of the Λ\Lambda-seperation energy of the Λ17_{\Lambda}^{17}O hypernucleus

Variational Monte Carlo calculations have been made for the $_{ \Lambda}^{17}$O hypernucleus using realistic two- and three-baryon interactions. A two pion exchange potential with spin- and space-exchange components is used for the $\Lambda$N potential. Three-body two-pion exchange and strongly repulsive dispersive $\Lambda$NN interactions are also included. The trial wave function is constructed from pair- and triplet-correlation operators acting on a single particle determinant. These operators consist of central, spin, isospin, tensor and three- baryon potential components. A cluster Monte Carlo method is developed for noncentral correlations and is used with up to four-baryon clusters in our calculations. The three-baryon $\Lambda$NN force is discussed.Comment: 24 pages, 2 figs available by fax., for publication in Phys. Rev.

Quasar Selection Based on Photometric Variability

We develop a method for separating quasars from other variable point sources using SDSS Stripe 82 light curve data for ~10,000 variable objects. To statistically describe quasar variability, we use a damped random walk model parametrized by a damping time scale, tau, and an asymptotic amplitude (structure function), SF_inf. With the aid of an SDSS spectroscopically confirmed quasar sample, we demonstrate that variability selection in typical extragalactic fields with low stellar density can deliver complete samples with reasonable purity (or efficiency, E). Compared to a selection method based solely on the slope of the structure function, the inclusion of the tau information boosts E from 60% to 75% while maintaining a highly complete sample (98%) even in the absence of color information. For a completeness of C=90%, E is boosted from 80% to 85%. Conversely, C improves from 90% to 97% while maintaining E=80% when imposing a lower limit on tau. With the aid of color selection, the purity can be further boosted to 96%, with C= 93%. Hence, selection methods based on variability will play an important role in the selection of quasars with data provided by upcoming large sky surveys, such as Pan-STARRS and the Large Synoptic Survey Telescope (LSST). For a typical (simulated) LSST cadence over 10 years and a photometric accuracy of 0.03 mag (achieved at i~22), C is expected to be 88% for a simple sample selection criterion of tau>100 days. In summary, given an adequate survey cadence, photometric variability provides an even better method than color selection for separating quasars from stars.Comment: (v2) 50 pages, accepted to Ap

Relativistic Beaming and the Intrinsic Properties of Extragalactic Radio Jets

Relations between the observed quantities for a beamed radio jet, apparent transverse speed and apparent luminosity (beta_app,L), and the intrinsic quantities, Lorentz factor and intrinsic luminosity (gamma,L_o), are investigated. The inversion from measured to intrinsic values is not unique, but approximate limits to gamma and L_o can be found using probability arguments. Roughly half the sources in a flux density--limited, beamed sample have a value of gamma close to the measured beta_app. The methods are applied to observations of 119 AGN jets made with the VLBA at 15 GHz during 1994-2002. The results strongly support the common relativistic beam model for an extragalactic radio jet. The (beta_app,L) data are closely bounded by a theoretical envelope, an aspect curve for gamma=32, L_o= 10^25 W/Hz. This gives limits to the maximum values of gamma and L_o in the sample: gamma_max about 32, and L_o,max ~ 10^26 W/Hz. No sources with both high beta_app and low L are observed. This is not the result of selection effects due to the observing limits, which are flux density S>0.5 Jy, and angular velocity mu<4 mas/yr. Many of the fastest quasars have a pattern Lorentz factor gamma_p close to that of the beam, gamma_b, but some of the slow quasars must have gamma_p<<gamma_b. Three of the 10 galaxies in the sample have a superluminal feature, with speeds up to beta_app about 6. The others are at most mildly relativistic. The galaxies are not off-axis versions of the powerful quasars, but Cygnus A might be an exception.Comment: 12 pages, 9 figures, 1 table, accepted for publication in the Astrophysical Journa