Orbital evolution of binary systems in dense stellar clusters is important in a variety of contexts: origin of blue stragglers, progenitors of compact object mergers, millisecond pulsars, and so on. Here we consider the general problem of secular evolution of the orbital elements of a binary system driven by the smooth tidal field of an axisymmetric stellar cluster (globular, nuclear, etc.) in which the binary orbits. We derive a secular Hamiltonian (averaged over both the inner Keplerian orbit of the binary and its outer orbit within the cluster) valid to quadrupole order for an arbitrary cluster potential and explore its characteristics. This doubly-averaged 'tidal' Hamiltonian depends on just two parameters, which fully absorb the information about the background cluster potential and the binary's orbit within it: a dimensional parameter $A$ setting the secular timescale, and a dimensionless parameter $\Gamma$ which determines the phase portrait of the binary's inner orbital evolution. We examine the dependence of $A$ and $\Gamma$ on cluster potential (both spherical and axisymmetric) and on the binary orbit within the cluster. Our theory reproduces known secular results - such as Lidov-Kozai evolution and the effect of the Galactic tide on Oort Cloud comets - in appropriate limits, but is more general. It provides a universal framework for understanding dynamical evolution of various types of binaries driven by the smooth tidal field of any axisymmetric potential. In a companion paper (Hamilton & Rafikov 2019b) we provide a detailed exploration of the resulting orbital dynamics.

The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive black hole) until gravitational wave emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH-BH, NS-BH and NS-NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several pc. For NS-NS binaries we find merger rates in the range $0.01-0.07$ Gpc$^{-3}$ yr$^{-1}$ from globular clusters and $0.1-0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters. For NS-BH and BH-BH binaries we find small merger rates from globular clusters, but a rate of $0.1 - 0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several per cent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future gravitational wave detectors.

Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $\Gamma$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $\Gamma$ is varied. We find that for $\Gamma > 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $\Gamma = 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $\gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters.

The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

We examine the light and color evolution of the T Tauri binary KH 15D through photometry obtained at wavelengths between 0.55 and 8.0 $\mu$m. The data were collected with ANDICAM on the 1.3 m SMARTS telescope at Cerro-Tololo Inter-American Observatory and with IRAC on the Spitzer Space Telescope. We show that the system's circumbinary ring, which acts as a screen that covers and uncovers different portions of the binary orbit as the ring precesses, has reached an orientation where the brighter component (star B) fully or nearly fully emerges during each orbital cycle. The fainter component (star A) remains fully occulted by the screen at all phases. The leading and trailing edges of the screen move across the sky at the same rate of $\sim$15 meters per second, consistent with expectation for a ring with a radius and width of $\sim$4 AU and a precession period of $\sim$6500 years. Light and color variations continue to indicate that the screen is sharp edged and opaque at \emph{VRIJH} wavelengths. However, we find an increasing transparency of the ring edge at 2.2, 3.6, and 4.5 $\mu$m. Reddening seen at the beginning of the eclipse that occurred during the CSI 2264 campaign particularly suggests selective extinction by a population of large dust grains. Meanwhile, the gradual bluing observed while star B is setting is indicative of forward scattering effects at the edge of the ring. The SED of the system at its bright phase shows no evidence of infrared excess emission that can be attributed to radiation from the ring or other dust component out to 8 microns.

Following two years of complete occultation of both stars by its opaque circumbinary ring, the binary T Tauri star KH 15D has abruptly brightened again during apastron phases, reaching I = 15 mag. Here, we show that the brightening is accompanied by a change in spectral class from K6/K7 (the spectral class of star A) to ~K1, and a bluing of the system in V-I by about 0.3 mag. A radial velocity measurement confirms that, at apastron, we are now seeing direct light from star B, which is more luminous and of earlier spectral class than star A. Evidently, the trailing edge of the occulting screen has just become tangent to one anse of star B's projected orbit. This confirms a prediction of the precession models, supports the view that the tilted ring is self-gravitating, and ushers in a new era of the system's evolution that should be accompanied by the same kind of dramatic phenomena observed from 1995-2009. It also promotes KH 15D from a single-lined to a double-lined eclipsing binary, greatly enhancing its value for testing pre-main sequence models. The results of our study strengthen the case for truncation of the outer ring at around 4 AU by a sub-stellar object such as an extremely young giant planet. The system is currently at an optimal configuration for detecting the putative planet and we urge expedient follow-up observations.

We test analytic predictions from different models of magnetospheric accretion, which invoke disk-locking, using stellar and accretion parameters derived from models of low resolution optical spectra of 36 T Tauri stars (TTSs) in NGC 2264 (age~3 Myrs). Little evidence is found for models that assume purely dipolar field geometries; however, strong support is found in the data for a modified version of the X-wind model (Shu et al. 1994) which allows for non-dipolar field geometries. The trapped flux concept in the X-wind model is key to making the analytic predictions which appear supported in the data. By extension, our analysis provides support for the outflows predicted by the X-wind as these also originate in the trapped flux region. In addition, we find no support in the data for accretion powered stellar winds from young stars. By comparing the analysis presented here of NGC 2264 with a similar analysis of stars in Taurus (age~1-2 Myr), we find evidence that the equilibrium interaction between the magnetic field and accretion disk in TTS systems evolves as the stars grow older, perhaps as the result of evolution of the stellar magnetic field geometry. We compare the accretion rates we derive with accretion rates based on U-band excess, finding good agreement. In addition, we use our accretion parameters to determine the relationship between accretion and H-beta luminosity, again finding good agreement with previously published results; however, we also find that care must be used when applying this relationship due to strong chromospheric emission in young stars which can lead to erroneous results in some cases.

We have obtained 48 high resolution echelle spectra of the pre-main sequence eclipsing binary system KH~15D (V582 Mon, P = 48.37 d, $e$ $\sim$ 0.6, M$_{A}$ = 0.6 M$_{\odot}$, M$_{B}$ = 0.7 M$_{\odot}$). The eclipses are caused by a circumbinary disk seen nearly edge on, which at the epoch of these observations completely obscured the orbit of star B and a large portion of the orbit of star A. The spectra were obtained over five contiguous observing seasons from 2001/2002 to 2005/2006 while star A was fully visible, fully occulted, and during several ingress and egress events. The H$\alpha$ line profile shows dramatic changes in these time series data over timescales ranging from days to years. A fraction of the variations are due to "edge effects" and depend only on the height of star A above or below the razor sharp edge of the occulting disk. Other observed variations depend on the orbital phase: the H$\alpha$ emission line profile changes from an inverse P Cygni type profile during ingress to an enhanced double-peaked profile, with both a blue and red emission component, during egress. Each of these interpreted variations are complicated by the fact that there is also a chaotic, irregular component present in these profiles. We find that the complex data set can be largely understood in the context of accretion onto the stars from a circumbinary disk with gas flows as predicted by the models of eccentric T Tauri binaries put forward by Artymowicz & Lubow, G\"{u}nther & Kley, and de Val-Borro et al. In particular, our data provide strong support for the pulsed accretion phenomenon, in which enhanced accretion occurs during and after perihelion passage.

We obtained J-, H- and JH-band photometry of known extrasolar planet transiting systems at the 2.1-m Kitt Peak National Observatory Telescope using the FLAMINGOS infrared camera between October 2008 and October 2011. From the derived lightcurves we have extracted the mid-transit times, transit depths and transit durations for these events. The precise mid-transit times obtained help improve the orbital periods and also constrain transit-time variations of the systems. For most cases the published system parameters successfully accounted for our observed lightcurves, but in some instances we derive improved planetary radii and orbital periods. We complemented our 2.1-m infrared observations using CCD z'-band and B-band photometry (plus two Hydrogen Alpha filter observations) obtained with the Kitt Peak Visitor's Center telescope, and with four H-band transits observed in October 2007 with the NSO's 1.6-m McMath-Pierce Solar Telescope. The principal highlights of our results are: 1) our ensemble of J-band planetary radii agree with optical radii, with the best-fit relation being: (Rp/R*)J = 0.0017 + 0.979 (Rp/R*)optical, 2) We observe star spot crossings during the transit of WASP-11/HAT-P-10, 3) we detect star spot crossings by HAT-P-11b (Kepler-3b), thus confirming that the magnetic evolution of the stellar active regions can be monitored even after the Kepler mission has ended, and 4) we confirm a grazing transit for HAT-P-27/WASP-40. In total we present 57 individual transits of 32 known exoplanet systems.

Photometry of the unique pre-main sequence binary system KH 15D is presented, spanning the years 2005-2010. This system has exhibited photometric variations and eclipses over the last 50 years caused by a precessing circumbinary disk. Advancement of the occulting edge across the binary orbit has continued and the photospheres of both stars are now completely obscured at all times. The system is now visible only by scattered light, and yet it continues to show a periodic variation on the orbital cycle with an amplitude exceeding two magnitudes. This variation, which depends only on the binary phase, has likely been present in the data since at least 1995. It can, by itself, account for shoulders on the light curve prior to ingress and following egress, obviating the need for components of extant models such as a scattering halo around star A or forward scattering from a fuzzy disk edge. A plausible source for the variable scattering component is reflected light from the far side of a warped occulting disk. We have detected color changes in V-I of several tenths of a magnitude to both the blue and red that occur during times of minima. These may indicate the presence of a third source of light (faint star) within the system, or a change in the reflectance properties of the disk as the portion being illuminated varies with the orbital motion of the stars. The data support a picture of the circumbinary disk as a geometrically thin, optically thick layer of perhaps mm or cm-sized particles that has been sculpted by the binary stars and possibly other components into a decidedly nonplanar configuration. A simple (infinitely sharp) knife-edge model does a good job of accounting for all of the recent (2005-2010) occultation data.

Na I D lines in the spectrum of the young binary KH 15D have been analyzed in detail. We find an excess absorption component that may be attributed to foreground interstellar absorption, and to gas possibly associated with the solids in the circumbinary disk. The derived column density is log N_NaI = 12.5 cm^-2, centered on a radial velocity that is consistent with the systemic velocity. Subtracting the likely contribution of the ISM leaves log N_NaI ~ 12.3 cm^-2. There is no detectable change in the gas column density across the "knife edge" formed by the opaque grain disk, indicating that the gas and solids have very different scale heights, with the solids being highly settled. Our data support a picture of this circumbinary disk as being composed of a very thin particulate grain layer composed of millimeter-sized or larger objects that are settled within whatever remaining gas may be present. This phase of disk evolution has been hypothesized to exist as a prelude to the formation of planetesimals through gravitational fragmentation, and is expected to be short-lived if much gas were still present in such a disk. Our analysis also reveals the presence of excess Na I emission relative to the comparison spectrum at the radial velocity of the currently visible star that plausibly arises within the magnetosphere of this still-accreting young star.

We present evidence that the spectroscopically identified bipolar jets of the pre-main sequence binary KH 15D are a common product of the whole binary system, rather than being launched from either star individually. They may be launched from the innermost part of the circumbinary disk (CBD) or may result from the merging of two outflows driven by the individual stars. This evidence is based on high-resolution H-alpha and [OI] 6300A line profiles obtained during eclipse phases of this nearly edge-on system. The occultation of star A (the only currently visible star) by the disk strongly suppresses the stellar H-alpha and continuum emission and allows one to study the faint redshifted and blueshifted emission components of the bipolar jets. The strongest evidence for jet production by the whole binary system comes from the observed radial velocity symmetry of the two jet components relative to the systemic velocity of the binary, in combination with current accretion models from the CBD onto a binary system.

We show that grains have grown to ~mm size (sand sized) or larger in the terrestrial zone (within ~3 AU) of the protoplanetary disk surrounding the 3 Myr old binary star KH 15D. We also argue that the reflected light in the system reaches us by back scattering off the far side of the same ring whose near side causes the obscuration.

The unusual flux variations of the pre-main-sequence binary star KH 15D have been attributed to occultations by a circumbinary disk. We test whether or not this theory is compatible with newly available data, including recent radial velocity measurements, CCD photometry over the past decade, and photographic photometry over the past 50 years. We find the model to be successful, after two refinements: a more realistic motion of the occulting feature, and a halo around each star that probably represents scattering by the disk. The occulting feature is exceptionally sharp-edged, raising the possibility that the dust in the disk has settled into a thin layer, and providing a tool for fine-scale mapping of the immediate environment of a T Tauri star. However, the window of opportunity is closing, as the currently visible star may be hidden at all orbital phases by as early as 2008.

We present results from the most recent (2002-2004) observing campaigns of the eclipsing system KH 15D, in addition to re-reduced data obtained at Van Vleck Observatory (VVO) between 1995 and 2000. Phasing nine years of photometric data shows substantial evolution in the width and depth of the eclipses. The most recent data indicate that the eclipses are now approximately 24 days in length, or half the orbital period. These results are interpreted and discussed in the context of the recent models for this system put forward by Winn et al. and Chiang & Murray-Clay. A periodogram of the entire data set yields a highly significant peak at 48.37 +/- 0.01 days, which is in accord with the spectroscopic period of 48.38 +/- 0.01 days determined by Johnson et al. Another significant peak, at 9.6 days, was found in the periodogram of the out-of-eclipse data at two different epochs. We interpret this as the rotation period of the visible star and argue that it may be tidally locked in pseudosynchronism with its orbital motion. If so, application of Hut's theory implies that the eccentricity of the orbit is e = 0.65 +/- 0.01. Analysis of the UVES/VLT spectra obtained by Hamilton et al. shows that the v sin(i) of the visible star in this system is 6.9 +/- 0.3 km/sec. Using this value of v sin(i) and the measured rotation period of the star, we calculate the lower limit on the radius to be R = (1.3 +/- 0.1), R_Sun, which concurs with the value obtained by Hamilton et al. from its luminosity and effective temperature. Here we assume that i = 90 degrees since it is likely that the spin and orbital angular momenta vectors are nearly aligned.

We present the results of a high-resolution spectroscopic monitoring program of the eclipsing pre-main-sequence star KH 15D that reveal it to be a single-line spectroscopic binary. We find that the best-fit Keplerian model has a period P = 48.38 days, which is nearly identical to the photometric period. Thus, we find the best explanation for the periodic dimming of KH 15D is that the binary motion carries the currently visible star alternately above and below the edge of an obscuring cloud. The data are consistent with the models involving an inclined circumstellar disk, as recently proposed by Winn et al. (2004) and Chiang & Murray-Clay (2004). We show that the mass ratio expected from models of PMS evolution, together with the mass constraints for the visible star, restrict the orbital eccentricity to 0.68 < e < 0.80 and the mass function to 0.125 < Fm < 0.5 Msun.

We present high resolution (R $\sim$ 44,000) UVES spectra of the eclipsing pre-main sequence star KH 15D covering the wavelength range 4780 to 6810 {\AA} obtained at three phases: out of eclipse, near minimum light and during egress. The system evidently acts like a natural coronagraph, enhancing the contrast relative to the continuum of hydrogen and forbidden emission lines during eclipse. At maximum light the H$\alpha$ equivalent width was $\sim$2 {\AA} and the profile showed broad wings and a deep central absorption. During egress the equivalent width was much higher ($\sim$70 {\AA}) and the broad wings, which extend to $\pm$ 300 km/s, were prominent. During eclipse totality the equivalent width was less than during egress ($\sim$40 {\AA}) and the high velocity wings were much weaker. H$\beta$ showed a somewhat different behavior, revealing only the blue-shifted portion of the high velocity component during eclipse and egress. [OI] $\lambda\lambda$6300, 6363 lines are easily seen both out of eclipse and when the photosphere is obscured and exhibit little or no flux variation with eclipse phase. Our interpretation is that KH 15D, although clearly a weak-line T Tauri star by the usual criteria, is still accreting matter from a circumstellar disk, and has a well-collimated bipolar jet. As the knife-edge of the occulting matter passes across the close stellar environment it is evidently revealing structure in the magnetosphere of this pre-main sequence star with unprecedented spatial resolution. We also show that there is only a small, perhaps marginally significant, change in the velocity of the K7 star between the maximum light and egress phases probed here.

Results of an international campaign to photometrically monitor the unique pre-main sequence eclipsing object KH 15D are reported. An updated ephemeris for the eclipse is derived that incorporates a slightly revised period of 48.36 d. There is some evidence that the orbital period is actually twice that value, with two eclipses occurring per cycle. The extraordinary depth (~3.5 mag) and duration (~18 days) of the eclipse indicate that it is caused by circumstellar matter, presumably the inner portion of a disk. The eclipse has continued to lengthen with time and the central brightness reversals are not as extreme as they once were. V-R and V-I colors indicate that the system is slightly bluer near minimum light. Ingress and egress are remarkably well modeled by the passage of a knife-edge across a limb-darkened star. Possible models for the system are briefly discussed.

Photometry and spectroscopy of the unique pre-main sequence eclipsing object KH15D in the young cluster NGC 2264 are presented. The orbital period is 48.34 days and both the length (~16 d) and depth (~3 mag) of the eclipse have increased with time. A brightening near the time of central eclipse is confirmed in the recent data but at a much smaller amplitude than was originally seen. Spectra taken when the star is bright show that the primary is a weak T Tauri star of spectral type K7. During eclipse there is no detectable change in spectral type or reddening, indicating that the obscuration is caused by rather large dust grains and/or macroscopic objects. Evidently the star is eclipsed by an extended feature in its circumstellar disk orbiting with a semi-major axis of ~0.2 AU. Continued photometric monitoring should allow us to probe the disk structure with a spatial resolution of ~3 x 10^6 km or better.