A procedure to analyze nonlinear density waves in Saturn's rings using several occultation profiles

Cassini radio science experiments have provided multiple occultation optical depth profiles of Saturn's rings that can be used in combination to analyze density waves. This paper establishes an accurate procedure of inversion of the wave profiles to reconstruct the wave kinematic parameters as a function of semi-major axis, in the nonlinear regime. This procedure is achieved from simulated data in the presence of realistic noise perturbations, to control the reconstruction error. By way of illustration we have applied our procedure to the Mimas 5:3 density wave. We were able to recover precisely the kinematic parameters from the radio experiment occultation data in most of the propagation region; a preliminary analysis of the pressure-corrected dispersion allowed us to determine new but still uncertain values for the opacity ($K\simeq 0.02$ cm$^2$/g) and velocity dispersion of ($c_o\simeq 0.6$ cm/s) in the wave region. Our procedure constitutes the first step in our planned analysis of the density waves of Saturn's rings. It is very accurate and efficient in the far-wave region. However, improvements are required within the first wavelength. The ways in which this method can be used to establish diagnostics of ring physics are outlined.Comment: 50 pages,13 figures, 2 tables. Published in Icarus

Models of Ultraluminous X-Ray Sources with Intermediate-Mass Black Holes

We have computed models for ultraluminous X-ray sources ("ULXs") consisting of a black-hole accretor of intermediate mass ("IMBH"; e.g., ~1000 Msun) and a captured donor star. For each of four different sets of initial donor masses and orbital separations, we computed 30,000 binary evolution models using a full Henyey stellar evolution code. To our knowledge this is the first time that a population of X-ray binaries this large has been carried out with other than approximation methods, and it serves to demonstrate the feasibility of this approach to large-scale population studies of mass-transfer binaries. In the present study, we find that in order to have a plausible efficiency for producing active ULX systems with IMBHs having luminosities > 10^{40} ergs/sec, there are two basic requirements for the capture of companion/donor stars. First, the donor stars should be massive, i.e., > 8 Msun. Second, the initial orbital separations, after circularization, should be close, i.e., < 6-30 times the radius of the donor star when on the main sequence. Even under these optimistic conditions, we show that the production rate of IMBH-ULX systems may fall short of the observed values by factors of 10-100.Comment: 5 pages, 2 figures, submitted to Ap

Exploring autoionization and photo-induced proton-coupled electron transfer pathways of phenol in aqueous solution

The excited state dynamics of phenol in water have been investigated using transient absorption spectroscopy. Solvated electrons and vibrationally cold phenoxyl radicals are observed upon 200 and 267 nm excitation, but with formation time scales that differ by more than 4 orders of magnitude. The impact of these findings is assessed in terms of the relative importance of autoionization versus proton-coupled electron transfer mechanisms in this computationally tractable model system

Can Vertical Profiles of Tropospheric Methane on Titan Be Derived from Radio-Occultation Soundings?

The intensity of the received signal at Earth in the radio occultations of Titan is attenuated both by refractive defocusing and pressure-induced absorption from N2-N2 and CH4-N2 pairs. Because the absorption strength is different for the two sets of pairs, matching the retrieved absorptivity profile can in principle yield the vertical variation in gaseous methane in the troposphere. There are two factors that make this difficult. The first is the propagation of noise in the phase and amplitude of the received signal in the absorption retrieval. The phase data is first inverted to retrieve vertical profiles of refractivity, from which the refractive defocusing is calculated. This is then subtracted from the observed. intensity attenuation of the received signal to generate a profile of atmospheric absorption. The second problem is the uncertainty in the pressure-induced absorption coefficients. Laboratory data at radio wavelengths is only available near room temperature (see, e.g., [1] for N2-N2), and the extrapolation to the low temperatures in Titan's troposphere is not well established. Ab initio calculations by Borysow et al. [2, 3] provide absorption coefficients at low temperatures and long wavelengths, but their accuracy has come into question. We present examples from Cassini radio occultations of Titan to illustrate the difficulties. For methane mole fractions in the lower troposphere comparable to that inferred from the Huygens probe (approximately 0.05), it will be difficult to separate the contributions of N2-N2 collisions from those of N2-CH4, collisions to the retrieved absorption. However, higher concentrations of CH4 and/or a higher signal-to-noise ratio from a future uplink experiment could result in a successful separation of the two components. However, key to this are highly accurate estimates of the absorption from a combination of laboratory measurements at love temperatures and long wavelengths, and possibly improved theoretical calculations

Spatio-Temporal Pattern of Saturn's Equatorial Oscillation

Recent ground-based and Cassini CIRS thermal-infrared data have characterized the spatial and temporal characteristics of an equatorial oscillation in the middle atmosphere of Saturn above the 100-mbar level. The CIRS data [I] indicated a pattern of warm and cold anomalies near the equator, stacked vertically in alternating fashion. The ground-based observations s2, although not having the altitude range or vertical resolution of the CIRS observations, covered several years and indicated an oscillation cycle of approx.15 years, roughly half of Saturn's year. In Earth's middle atmosphere, both the quasi-biennial (approx.26 months) and semi-annual equatorial oscillations have been extensively observed and studied (see e.g., [3]), These exhibit a pattern of alternating warmer and cooler zonal-mean temperatures with altitude, relative to those at subtropical latitudes. Consistent with the thermal wind equation, this is also associated with an alternating pattern of westerly and easterly zonal winds. Moreover, the pattern of winds and temperatures descends with time. Momentum deposition by damped vertically propagating waves is thought to play a key role m forcing both types of oscillation, and it can plausibly account for the descent. Here we report the direct observation of this descent in Saturn's equatorial atmosphere from Cassini radio occultation soundings in 2005 and 2009. The retrieved temperatures are consistent with a descent of 0.7 x the pressure scale height. The descent rate is related to the magnitude of the wave forcing, radiative damping, and induced meridional circulations. We discuss possible implications

Magnetically Torqued Thin Accretion Disks

We compute the properties of a geometrically thin, steady accretion disk surrounding a central rotating, magnetized star. The magnetosphere is assumed to entrain the disk over a wide range of radii. The model is simplified in that we adopt two (alternate) ad hoc, but plausible, expressions for the azimuthal component of the magnetic field as a function of radial distance. We find a solution for the angular velocity profile tending to corotation close to the central star, and smoothly matching a Keplerian curve at a radius where the viscous stress vanishes. The value of this ''transition'' radius is nearly the same for both of our adopted B-field models. We then solve analytically for the torques on the central star and for the disk luminosity due to gravity and magnetic torques. When expressed in a dimensionless form, the resulting quantities depend on one parameter alone, the ratio of the transition radius to the corotation radius. For rapid rotators, the accretion disk may be powered mostly by spin-down of the central star. These results are independent of the viscosity prescription in the disk. We also solve for the disk structure for the special case of an optically thick alpha disk. Our results are applicable to a range of astrophysical systems including accreting neutron stars, intermediate polar cataclysmic variables, and T Tauri systems.Comment: 9 sharper figs, updated reference

EPIC 220204960: A Quadruple Star System Containing Two Strongly Interacting Eclipsing Binaries

We present a strongly interacting quadruple system associated with the K2 target EPIC 220204960. The K2 target itself is a Kp = 12.7 magnitude star at Teff ~ 6100 K which we designate as "B-N" (blue northerly image). The host of the quadruple system, however, is a Kp = 17 magnitude star with a composite M-star spectrum, which we designate as "R-S" (red southerly image). With a 3.2" separation and similar radial velocities and photometric distances, 'B-N' is likely physically associated with 'R-S', making this a quintuple system, but that is incidental to our main claim of a strongly interacting quadruple system in 'R-S'. The two binaries in 'R-S' have orbital periods of 13.27 d and 14.41 d, respectively, and each has an inclination angle of >89 degrees. From our analysis of radial velocity measurements, and of the photometric lightcurve, we conclude that all four stars are very similar with masses close to 0.4 Msun. Both of the binaries exhibit significant ETVs where those of the primary and secondary eclipses 'diverge' by 0.05 days over the course of the 80-day observations. Via a systematic set of numerical simulations of quadruple systems consisting of two interacting binaries, we conclude that the outer orbital period is very likely to be between 300 and 500 days. If sufficient time is devoted to RV studies of this faint target, the outer orbit should be measurable within a year.Comment: 20 pages, 18 figures, 7 tables; accepted for publication in MNRA

The Effects of Binary Evolution on the Dynamics of Core Collapse and Neutron-Star Kicks

We systematically examine how the presence in a binary affects the final core structure of a massive star and its consequences for the subsequent supernova explosion. Interactions with a companion star may change the final rate of rotation, the size of the helium core, the strength of carbon burning and the final iron core mass. Stars with initial masses larger than \sim 11\Ms that experiece core collapse will generally have smaller iron cores at the time of the explosion if they lost their envelopes due to a previous binary interaction. Stars below \sim 11\Ms, on the other hand, can end up with larger helium and metal cores if they have a close companion, since the second dredge-up phase which reduces the helium core mass dramatically in single stars does not occur once the hydrogen envelope is lost. We find that the initially more massive stars in binary systems with masses in the range 8 - 11\Ms are likely to undergo an electron-capture supernova, while single stars in the same mass range would end as ONeMg white dwarfs. We suggest that the core collapse in an electron-capture supernova (and possibly in the case of relatively small iron cores) leads to a prompt explosion rather than a delayed neutrino-driven explosion and that this naturally produces neutron stars with low-velocity kicks. This leads to a dichotomous distribution of neutron star kicks, as inferred previously, where neutron stars in relatively close binaries attain low kick velocities. We illustrate the consequences of such a dichotomous kick scenario using binary population synthesis simulations and discuss its implications. This scenario has also important consequences for the minimum initial mass of a massive star that becomes a neutron star. (Abbreviated.)Comment: 8 pages, 3 figures, submitted to ApJ, updated versio

The K2-ESPRINT Project. I. Discovery of the Disintegrating Rocky Planet K2-22b with a Cometary Head and Leading Tail

We present the discovery of a transiting exoplanet candidate in the K2 Field-1 with an orbital period of 9.1457 hr: K2-22b. The highly variable transit depths, ranging from $\sim$0\% to 1.3\%, are suggestive of a planet that is disintegrating via the emission of dusty effluents. We characterize the host star as an M-dwarf with $T_{\rm eff} \simeq 3800$ K. We have obtained ground-based transit measurements with several 1-m class telescopes and with the GTC. These observations (1) improve the transit ephemeris; (2) confirm the variable nature of the transit depths; (3) indicate variations in the transit shapes; and (4) demonstrate clearly that at least on one occasion the transit depths were significantly wavelength dependent. The latter three effects tend to indicate extinction of starlight by dust rather than by any combination of solid bodies. The K2 observations yield a folded light curve with lower time resolution but with substantially better statistical precision compared with the ground-based observations. We detect a significant "bump" just after the transit egress, and a less significant bump just prior to transit ingress. We interpret these bumps in the context of a planet that is not only likely streaming a dust tail behind it, but also has a more prominent leading dust trail that precedes it. This effect is modeled in terms of dust grains that can escape to beyond the planet's Hill sphere and effectively undergo `Roche lobe overflow,' even though the planet's surface is likely underfilling its Roche lobe by a factor of 2.Comment: 22 pages, 16 figures. Final version accepted to Ap