Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor ($Z_\odot/10$) single helium stars using the \texttt{MESA} software instrument. For PPISN progenitors with pre-pulse masses $>50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4\;$~yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.

Context: Cygnus X-1 is a black hole X-ray binary system in which the black hole captures and accretes gas from the strong stellar wind emitted by its supergiant O9.7 companion star. The irradiation of the supergiant star essentially determines the flow properties of the stellar wind and the X-ray luminosity from the system. The results of three-dimensional hydrodynamical simulations of wind-fed X-ray binary systems reported in recent work reveal that the ionizing feedback of the X-ray irradiation leads to the existence of two stable states with either a soft or a hard spectrum. Aims: We discuss the observed radiation of Cygnus X-1 in the soft and hard state in the context of mass flow in the corona and disk, as predicted by the recent application of a condensation model. Methods: The rates of gas condensation from the corona to the disk for Cygnus X-1 are determined, and the spectra of the hard and soft radiation are computed. The theoretical results are compared with the MAXI observations of Cygnus X-1 from 2009 to 2018. In particular, we evaluate the hardness-intensity diagrams (HIDs) for its ten episodes of soft and hard states which show that Cygnus X-1 is distinct in its spectral changes as compared to those found in the HIDs of low-mass X-ray binaries. Results: The theoretically derived values of photon counts and hardness are in approximate agreement with the observed data in the HID. However, the scatter in the diagram is not reproduced. Improved agreement could result from variations in the viscosity associated with clumping in the stellar wind and corresponding changes of the magnetic fields in the disk. The observed dipping events in the hard state may also contribute to the scatter and to a harder spectrum than predicted by the model.

Magnetic fields of order $10^1-10^2$ gauss that are present in the envelopes of red giant stars are ejected in common envelope scenarios. These fields could be responsible for the launching of magnetically driven winds in proto-planetary nebulae. Using 2D simulations of magnetized winds interacting with an envelope drawn from a 3D simulation of the common envelope phase, we study the confinement, heating, and magnetic field development of post-common envelope winds. We find that the ejected magnetic field can be enhanced via compression by factors up to $\sim 10^4$ in circumbinary disks during the self-regulated phases. We find values for the kinetic energy of the order of $10^{46}$ erg that explain the large values inferred in proto-planetary nebula outflows. We show that the interaction of the formed circumbinary disk with a spherical, stellar wind produces a "tapered" flow that is almost indistinguishable from an imposed tapered flow. This increases the uncertainty of the origin of proto-planetary nebula winds, which could be either stellar, circumstellar (stellar accretion disk), circumbinary (circumbinary accretion disk), or a combination of all three. Within this framework, a scenario for self-collimation of weakly magnetized winds is discussed, which can explain the two objects where the collimation process is observationally resolved, HD 101584 and Hen 3-1475. An explanation for the equatorial, molecular hydrogen emission in CRL 2688 is also presented.

Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to model the inspiral of a 1.4M$_{\odot}$ neutron star (NS) inside the envelope of a 12$M_{\odot}$ red supergiant star. We self-consistently calculate the drag force experienced by the NS as well as the back-reaction onto the expanding envelope as the NS spirals in. Nearly all of the hydrogen envelope escapes, expanding to large radii ($\sim$10$^2$ AU) where it forms an optically thick envelope with temperatures low enough that dust formation occurs. We simulate the NS orbit until only 0.8M$_{\odot}$ of the hydrogen envelope remains around the giant star's core. Our results suggest that the inspiral will continue until another $\approx$0.3M$_{\odot}$ are removed, at which point the remaining envelope will retract. Upon separation, a phase of dynamically stable mass transfer onto the NS accretor is likely to ensue, which may be observable as an ultraluminous X-ray source. The resulting binary, comprised of a detached 2.6M$_{\odot}$ helium-star and a NS with a separation of 3.3-5.7R$_{\odot}$, is expected to evolve into a merging double neutron-star, analogous to those recently detected by LIGO/Virgo. For our chosen combination of binary parameters, our estimated final separation (including the phase of stable mass transfer) suggests a very high $\alpha_{\rm CE}$-equivalent efficiency of $\approx$5.

The morphological properties of the outflowing circumstellar envelopes surrounding mass-losing stars in eccentric binary systems are presented from a set of three-dimensional hydrodynamical model simulations. Based on four template models of the envelope viewed for a range of inclination angles of the systems, we implement interpretative tools for observations at high spectral/angular resolutions (as illustrated via velocity channel maps as well as position-velocity, radius-velocity, and angle-radius diagrams). Within this framework, the image and kinematical structures can be used to place constraints on the orbital parameters of the system. Specifically, three unique characteristic patterns in the envelopes are found that distinguish these systems from those in binary systems in circular orbits. Bifurcation of the spiral pattern, asymmetry in the interarm density depression, and a concurrent spiral/ring appearance all point to a binary system with an eccentric orbit. The methodology presented in this paper is illustrated in an analysis of recent radio observations of several asymptotic giant branch stars.

The structure and evolution of gas flows within the cavity of a circumbinary disk (CBD) surrounding the stellar components in eccentric binaries are examined via two-dimensional hydrodynamical simulations. The degree to which gas fills the cavity between the circumstellar disks (CSDs) and the CBD is found to be greater for highly eccentric systems, in comparison to low-eccentricity systems, reflecting the spatial extent over which mass enters into the cavity throughout the orbit. The pattern of the gas flow in the cavity differs for eccentric binaries from that of binaries in a circular orbit. In particular, the former reveals tightly wound gas streams and figure-eight-like structures for systems characterized by eccentricities, $e \ge 0.4$, whereas the latter only reveal relatively loosely bent streams from the CBD to the CSDs. Hence, the description of the stream structures can be a probe of sufficient non-circularity of the binary orbital motion. Given that the inner edge of the CBD is not very well defined for highly eccentric systems due to the complex gas structures, it is suggested that the area of the cavity for high-sensitivity imaging observations may prove to be a more useful diagnostic for probing the effectiveness of CBD clearing in the future.

The discovery via gravitational waves of binary black hole systems with total masses greater than $60M_\odot$ has raised interesting questions for stellar evolution theory. Among the most promising formation channels for these systems is one involving a common envelope binary containing a low metallicity, core helium burning star with mass $\sim 80-90M_\odot$ and a black hole with mass $\sim 30-40M_\odot$. For this channel to be viable, the common envelope binary must eject more than half the giant star's mass and reduce its orbital separation by as much as a factor of 80. We discuss issues faced in numerically simulating the common envelope evolution of such systems and present a 3D AMR simulation of the dynamical inspiral of a low-metallicity red supergiant with a massive black hole companion.

Based on a high quality $N$-body simulation of a double bar galaxy model, we investigate the evolution of the bar properties, including their size, strength and instantaneous pattern speed derived by using three distinct methods: the Fourier, Jacobi integral, and moment of inertia methods. The interaction of the two bars, which rotate at distinct speeds, primarily affects the size, strength and pattern speed of the inner bar. When the two bars are perpendicular to each other, the size and the pattern speed of the inner bar decrease and its strength increases. The emergence of a strong Fourier $m=1$ mode increases the oscillation amplitude of the size, strength and pattern speed of the inner bar. On the other hand, the characteristics of the outer bar are substantially influenced by its adjacent spiral structure. When the spiral structure disappears, the size of the outer bar increases and its strength and pattern speed decrease. Consequently, the ratio of the pattern speed of the outer bar with respect to the inner bar is not constant and increases with time. Overall, the double bar and disk system displays substantial high frequency semi-chaotic fluctuations of the pattern strengths and speeds both in space and time, superposed on the slow secular evolution, which invalidates the assumption that the actions of individual stars should be well conserved in barred galaxies, such as the Milky Way.

Among the four black hole binary merger events detected by LIGO, six progenitor black holes have masses greater than 20\,$M_\odot$. The existence of such massive BHs calls for extreme metal-poor stars as the progenitors. An alternative possibility that a pair of stellar mass black holes each with mass $\sim7\,M_\odot$ increases to $>20\,M_\odot$ via accretion from a disk surrounding a super massive black hole in an active galactic nucleus is considered. The growth of mass of the binary and the transfer of orbital angular momentum to the disk accelerates the merger. Based on the recent numerical work of Tang et al. (2017), it is found that, in the disk of a low mass AGN with mass $\sim10^6\,M_\odot$ and Eddington ratio $>0.01$, the mass of an individual BH in the binary can grow to $>20\,M_\odot$ before coalescence provided that accretion takes place at a rate more than 10 times the Eddington value. The mechanism predicts a new class of gravitational wave sources involving the merger of two extreme Kerr black holes associated with active galactic nuclei and a possible electromagnetic wave counterpart.

A new accretion picture based on a small disk surrounding a black hole is developed for the wind-fed source Cyg X-1. The hard and soft spectral states of Cyg X-1 are interpreted in terms of co-spatial two component flows for the innermost region of an accretion disk. The state transitions result from the outward expansion and inward recession of this inner disk for the hard to soft and soft to hard transition respectively. The theoretical framework for state transitions in black hole X-ray binaries with high mass companions involving a change in the inner disk size, thus, differs from systems with low mass companions involving the change in the outer disk size. This fundamental difference stems from the fact that matter captured and supplied to the black hole in wind-fed systems has low specific angular momentum and is hot essentially heated in the bow and spiral shocks, whereas it has high specific angular momentum and is cool in Roche lobe overflow systems. The existence of a weak cool disk around the ISCO region in the hard state allows for the presence of a relativistically broadened Fe K line. The small disk fed by gas condensation forms without an extensive outer disk, precluding thermal instabilities and large outbursts, resulting in the lack of large amplitude outbursts and hysteresis effects in the light curve of high mass black hole X-ray binaries. Their relatively persistent X-ray emission is attributed to their wind-fed nature.

The morphology of planetary nebulae emerging from the common envelope phase of binary star evolution is investigated. Using initial conditions based on the numerical results of hydrodynamical simulations of the common envelope phase it is found that the shapes and sizes of the resulting nebula are very sensitive to the effective temperature of the remnant core, the mass-loss rate at the onset of the common envelope phase, and the mass ratio of the binary system. These parameters are related to the efficiency of the mass ejection after the spiral-in phase, the stellar evolutionary phase (i.e., RG, AGB or TP-AGB), and the degree of departure from spherical symmetry in the stellar wind mass loss process itself respectively. It is found that the shapes are mostly bipolar in the early phase of evolution, but can quickly transition to elliptical and barrel-type shapes. Solutions for nested lobes are found where the outer lobes are usually bipolar and the inner lobes are elliptical, bipolar or barrel-type, a result due to the flow of the photo-evaporated gas from the equatorial region. It is found that the lobes can be produced without the need for two distinct mass ejection events. In all the computations, the bulk of the mass is concentrated in the orbital or equatorial plane, in the form of a large toroid, which can be either neutral (early phases) or photoionized (late phases), depending of the evolutionary state of the system.

The X-ray emission from bright active galactic nuclei (AGNs) is believed to originate in a hot corona lying above a cold, geometrically thin accretion disk. A highly concentrated corona located within $\sim10$ gravitational radii above the black hole is inferred from observations. Based on the accretion of interstellar medium/wind, a disk corona model has been proposed in which the corona is well coupled to the disk by radiation, thermal conduction, as well as by mass exchange \citep{Liu2015, Qiao2017}. Such a model avoids artificial energy input to the corona and has been used to interpret the spectral features observed in AGN. In this work, it is shown that the bulk emission size of the corona is very small for the extended accretion flow in our model. More than 80\% of the hard X-ray power is emitted from a small region confined within 10 Schwarzschild radii around a non-spinning black hole, which is expected to be even smaller accordingly for a spinning black hole. Here, the corona emission is more extended at higher Eddington ratios. The compactness parameter of the corona, $l={L\over R}{\sigma_{\rm T}\over m_{\rm e} c^3}$, is shown to be in the range of 1-33 for Eddington ratios of 0.02 - 0.1. Combined with the electron temperature in the corona, this indicates that electron--positron pair production is not dominant in this regime. A positive relation between the compactness parameter and photon index is also predicted. By comparing the above model predictions with observational features, we find that the model is in agreement with observations.

Preplanetary nebulae and planetary nebulae are evolved, mass-losing stellar objects that show a wide variety of morphologies. Many of these nebulae consist of outer structures that are nearly spherical (spiral/shell/arc/halo) and inner structures that are highly asymmetric (bipolar/multipolar). The coexistence of such geometrically distinct structures is enigmatic because it hints at the simultaneous presence of both wide and close binary interactions, a phenomenon that has been attributed to stellar binary systems with eccentric orbits. Here we report new high-resolution molecular-line observations of the circumstellar spiral-shell pattern of AFGL 3068, an asymptotic giant branch star transitioning to the preplanetary nebula phase. The observations clearly reveal that the dynamics of the mass loss is influenced by the presence of an eccentric-orbit binary. This quintessential object opens a new window on the nature of deeply embedded binary stars through the circumstellar spiral-shell patterns that reside at distances of several thousand Astronomical Units from the stars.

The effective potential neighboring the corotation resonance region in barred galaxies is shown to be strongly time-dependent in any rotating frame because of the competition of nearby perturbations of similar strengths with differing rotation speeds. Contrary to the generally adopted assumption, that in the bar rotating frame the corotation region should possess four stationary equilibrium points (Lagrange points), with high quality N-body simulations we localize the instantaneous equilibrium points and find that they circulate or oscillate broadly in azimuth with respect to the pattern speeds of the inner or outer perturbations. This implies that at the particle level the Jacobi integral is not well conserved around the corotation radius. That is, angular momentum exchanges decouple from energy exchanges, enhancing the chaotic diffusion of stars through the corotation region.

We extend the work of Yen et al. (2012) and develop 2nd order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a two-dimensional kernel derived for infinitesimally thin disks and is free of artificial boundary conditions. The self-gravitational force calculation is presented in generalized convolution forms for a nested grid configuration. A numerical technique derived from a fast Fourier transform is employed to reduce the computational complexity to be nearly linear. By comparing with analytic potential-density pairs associated with the generalized Maclaurin disks, the extended approach is verified to be of second order accuracy using numerical simulations. The proposed method is accurate, computationally fast and has the potential to be applied to the studies of planetary migration and the gaseous morphology of disk galaxies.

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.

It is commonly believed that the optical/UV and X-ray emissions in luminous AGN are produced in an accretion disk and an embedded hot corona respectively. The inverse Compton scattering of disk photons by hot electrons in the corona can effectively cool the coronal gas if the mass supply is predominantly via a cool disk like flow as in BHXRBs. Thus, the application of such a model to AGNs fails to produce their observed X-ray emission. As a consequence, a fraction of disk accretion energy is usually assumed to be transferred to the corona. To avoid this assumption, we propose that gas in a vertically extended distribution is supplied to a supermassive black hole by the gravitational capture of interstellar medium or stellar wind material. In this picture, the gas partially condenses to an underlying cool disk as it flows toward the black hole, releasing accretion energy as X-ray emission and supplying mass for the disk accretion. Detailed numerical calculations reveal that the X-ray luminosity can reach a few tens of percent of the bolometric luminosity. The value of $\alpha_{\rm ox}$ varies from 0.9 to 1.2 for the mass supply rate ranging from 0.03 to 0.1 times the Eddington value. The corresponding photon index in the 2-10 keV energy band varies from 1.9 to 2.3. Such a picture provides a natural extension of the model for low luminosity AGN where condensation is absent at low mass accretion rates and no optically thick disk exists in the inner region.

The shock interaction and evolution of nova ejecta with a wind from a red giant star in a symbiotic binary system are investigated via three-dimensional hydrodynamics simulations. We specifically model the March 2010 outburst of the symbiotic recurrent nova V407~Cygni from the quiescent phase to its eruption phase. The circumstellar density enhancement due to wind-white dwarf interaction is studied in detail. It is found that the density-enhancement efficiency depends on the ratio of the orbital speed to the red giant wind speed. Unlike another recurrent nova, RS~Ophiuchi, we do not observe a strong disk-like density enhancement, but instead observe an aspherical density distribution with $\sim 20\%$ higher density in the equatorial plane than at the poles. To model the 2010 outburst, we consider several physical parameters, including the red giant mass loss rate, nova eruption energy, and ejecta mass. A detailed study of the shock interaction and evolution reveals that the interaction of shocks with the red giant wind generates strong Rayleigh-Taylor instabilities. In addition, the presence of the companion and circumstellar density enhancement greatly alter the shock evolution during the nova phase. The ejecta speed after sweeping out most of the circumstellar medium decreases to $\sim 100-300$ km-s$^{-1}$, depending on model, which is consistent with the observed extended redward emission in [N~II] lines in April 2011.

Gas response to the underlying stellar spirals is explored for M81 using unmagnetized hydrodynamic simulations. Constrained within the uncertainty of observations, 18 simulations are carried out to study the effects of selfgravity and to cover the parameter space comprising three different sound speeds and three different arm strengths. The results are confronted with those data observed at wavelengths of 8 $\mu$m and 21 cm. In the outer disk, the ring-like structure observed in 8 $\mu$m image is consistent with the response of cold neutral medium with an effective sound speed 7 km s$^{-1}$, while for the inner disk, the presence of spiral shocks can be understood as a result of 4:1 resonances associated with the warm neutral medium with an effective sound speed 19 km s$^{-1}$. Simulations with single effective sound speed alone cannot simultaneously explain the structures in the outer and inner disks. This justifies the coexistence of cold and warm neutral media in M81. The anomalously high streaming motions observed in the northeast arm and the outward shifted turning points in the iso-velocity contours seen along the southwest arm are interpreted as signatures of interactions with companion galaxies. The level of simulated streaming motions narrows down the uncertainty of observed arm strength toward larger amplitudes.

The nature of the progenitor systems of type~Ia supernovae is still unclear. One way to distinguish between the single-degenerate scenario and double-degenerate scenario for their progenitors is to search for the surviving companions. Using a technique that couples the results from multi-dimensional hydrodynamics simulations with calculations of the structure and evolution of main-sequence- and helium-rich surviving companions, the color and magnitude of main-sequence- and helium-rich surviving companions are predicted as functions of time. The surviving companion candidates in Galactic type~Ia supernova remnants and nearby extragalactic type~Ia supernova remnants are discussed. We find that the maximum detectable distance of main-sequence surviving companions (helium-rich surviving companions) is $0.6-4$~Mpc ($0.4-16$~Mpc), if the apparent magnitude limit is 27 in the absence of extinction, suggesting that the Large and Small Magellanic Clouds and the Andromeda Galaxy are excellent environments in which to search for surviving companions. However, only five Ia~SNRs have been searched for surviving companions, showing little support for the standard channels in the singe-degenerate scenario. To better understand the progenitors of type Ia supernovae, we encourage the search for surviving companions in other nearby type Ia supernova remnants.

Aside from the grand-design stellar spirals appearing in the disk of M81, a pair of stellar spiral arms situated well inside the bright bulge of M81 has been recently discovered by Kendall et al. (2008). The seemingly unrelated pairs of spirals pose a challenge to the theory of spiral density waves. To address this problem, we have constructed a three component model for M81, including the contributions from a stellar disk, a bulge, and a dark matter halo subject to observational constraints. Given this basic state for M81, a modal approach is applied to search for the discrete unstable spiral modes that may provide an understanding for the existence of both spiral arms. It is found that the apparently separated inner and outer spirals can be interpreted as a single trailing spiral mode. In particular, these spirals share the same pattern speed 25.5 km s$^{-1}$ kpc$^{-1}$ with a corotation radius of 9.03 kpc. In addition to the good agreement between the calculated and the observed spiral pattern, the variation of the spiral amplitude can also be naturally reproduced.

With the advent of high-resolution high-sensitivity observations, spiral patterns have been revealed around several asymptotic giant branch (AGB) stars. Such patterns can provide possible evidence for the existence of central binary stars embedded in outflowing circumstellar envelopes. Here, we suggest the viability of explaining the previously observed incomplete ring-like patterns with the spiral-shell structure due to the motion of (unknown) binary components viewed at an inclination with respect to the orbital plane. We describe a method of extracting such spiral-shells from an incomplete ring-like pattern to place constraints on the characteristics of the central binary stars. The use of gas kinematics is essential in facilitating a detailed modeling for the three-dimensional structure of the circumstellar pattern. We show that a hydrodynamic radiative transfer model can reproduce the structure of the HC3N molecular line emission of the extreme carbon star, CIT 6. This method can be applied to other sources in the AGB phase and to the outer ring-like patterns of pre-planetary nebulae for probing the existence of embedded binary stars, which are highly anticipated with future observations using the Atacama Large Millimeter/submillimeter Array.

A model for the formation and evolution of binary millisecond radio pulsars in systems with low mass companions (< 0.1 Msun) is investigated using a binary population synthesis technique. Taking into account the non conservative evolution of the system due to mass loss from an accretion disk as a result of propeller action and from the companion via ablation by the pulsar, the transition from the accretion powered to rotation powered phase is investigated. It is shown that the operation of the propeller and ablation mechanisms can be responsible for the formation and evolution of black widow millisecond pulsar systems from the low mass X-ray binary phase at an orbital period of ~0.1 day. For a range of population synthesis input parameters, the results reveal that a population of black widow millisecond pulsars characterized by orbital periods as long as ~0.4 days and companion masses as low as ~0.005 Msun can be produced. The orbital periods and minimum companion mass of this radio millisecond pulsar population critically depend on the thermal bloating of the semi-degenerate hydrogen mass losing component, with longer orbital periods for a greater degree of bloating. Provided that the radius of the companion is increased by about a factor of 2 relative to a fully degenerate, zero temperature configuration, an approximate agreement between observed long orbital periods and theoretical modeling of hydrogen rich donors can be achieved. We find no discrepancy between the estimated birth rates for LMXBs and black widow systems, which on average are $\sim~1.3\times10^{-5} {\rm yr}^{-1}$ and $1.3\times10^{-7} {\rm yr}^{-1}$ respectively.

The magnitude of the viscosity and magnetic field parameters in hot accretion flows is investigated in low luminosity active galactic nuclei (LLAGNs). Theoretical studies show that a geometrically thin, optically thick disk is truncated at mass accretion rates less than a critical value by mass evaporated vertically from the disk to the corona, with the truncated region replaced by an advection dominated accretion flow (ADAF). The critical accretion rate for such a truncation is a function of the viscosity and magnetic field. Observations of X-ray photon indices and spectral fits of a number of LLAGNs published in the literature provide an estimate of the critical rate of mass accretion and the truncation radius respectively. By comparing the observational results with theoretical predictions, the viscosity and magnetic field parameters in the hot accretion flow region are estimated. Specifically, the mass accretion rates inferred in different sources constrain the viscosity parameter, whereas the truncation radii of the disk, as inferred from spectral fits, further constrain the magnetic field parameter. It is found that the value of the viscosity parameter in the corona/ADAF ranges from 0.17 to 0.5, with values clustered about 0.2-0.3. Magnetic pressure is required by the relatively small truncation radii for some LLAGNs and is found to be as high as its equipartition value with the gas pressure. The inferred values of the viscosity parameter are in agreement with those obtained from the observations of non-stationary accretion in stellar mass black hole X-ray transients. This consistency provides support for the paradigm that a geometrically thin disk is truncated by means of a mass evaporation process from the disk to the corona at low mass accretion rates.

The progenitor systems of Type Ia supernovae (SNe Ia) are still under debate. Based on recent hydrodynamics simulations, non-degenerate companions in the single-degenerate scenario (SDS) should survive the supernova impact. One way to distinguish between the SDS and the double-degenerate scenario is to search for the post-impact remnant stars (PIRSs) in SN Ia remnants. Using a technique that combines multi-dimensional hydrodynamics simulations with one-dimensional stellar evolution simulations, we have examined the post-impact evolution of helium-rich binary companions in the SDS. It is found that these helium-rich PIRSs (He PIRSs) dramatically expand and evolve to a luminous phase ($L\sim 10^4 L_\odot$) about 10 years after a supernova explosion. Subsequently, they contract and evolve to become hot blue-subdwarf-like (sdO-like) stars by releasing gravitational energy, persisting as sdO-like stars for several million years before evolving to the helium red-giant phase. We therefore predict that a luminous OB-like star should be detectable within $\sim 30$ years after the SN explosion. Thereafter, it will shrink and become an sdO-like star in the central regions of SN Ia remnants within star-forming regions for SN Ia progenitors evolved via the helium-star channel in the SDS. These He PIRSs are predicted to be rapidly rotating ($v_{\rm rot} \gtrsim 50$ km s$^{-1}$) and to have high spatial velocities ($v_{\rm linear} \gtrsim 500$ km s$^{-1}$). Furthermore, if supernova remnants have diffused away and are not recognizable at a later stage, He PIRSs could be an additional source of single sdO stars and/or hypervelocity stars.

NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in \citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_\sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.

The structure and evolution of the stellar velocity ellipsoid plays an important role in shaping galaxies undergoing bar driven secular evolution and the eventual formation of a boxy/peanut bulge such as present in the Milky Way. Using collisionless N-body simulations, we show that during the formation of such a boxy/peanut bulge, the meridional shear stress of stars, which can be measured by the meridional tilt of the velocity ellipsoid, reaches a characteristic peak in its time evolution. It is shown that the onset of a bar buckling instability is closely connected to the maximum meridional tilt of the stellar velocity ellipsoid. Our findings bring new insight to this complex gravitational instability of the bar which complements the buckling instability studies based on orbital models. We briefly discuss the observed diagnostics of the stellar velocity ellipsoid during such a phenomenon.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

The nature of the progenitor systems of Type Ia supernovae is still uncertain. One way to distinguish between the single-degenerate scenario (SDS) and double-degenerate scenario (DDS) is to search for the post-impact remnant star. To examine the characteristics of the post-impact remnant star, we have carried out three-dimensional hydrodynamic simulations of supernova impacts on main sequence-like stars. We explore the evolution of the post-impact remnants using the stellar evolution code MESA. We find that the luminosity and radius of the remnant star dramatically increase just after the impact. After the explosion, post-impact companions continue to expand on a progenitor-dependent timescale of $\sim 10^{2.5-3} $ yr before contracting. It is found that the time evolution of the remnant star is dependent on not only the amount of energy absorbed, but also on the depth of the energy deposition. We examine the viability of the candidate star Tycho G as the possible remnant companion in Tycho's supernova by comparing it to the evolved post-impact remnant stars in our simulations. The closest model in our simulations has a similar effective temperature, but the luminosity and radius are twice as large. By examining the angular momentum distribution in our simulations, we find that the surface rotational speed could drop to $\sim 10$ km s$^{-1}$ if the specific angular momentum is conserved during the post-impact evolution, implying that Tycho G cannot be completely ruled out because of its low surface rotation speed.

The binary characteristics of asymptotic giant branch (AGB) stars are imprinted in the asymmetric patterns of their circumstellar envelopes. We develop a simple method for constraining the orbital parameters of such binary stars from the characteristics of a spiral-like pattern observed at large distances from the central stars. We place constraints on the properties of AFGL 3068 (i.e., the masses of binary components, the viewing inclination of the orbital plane, as well as the orbital period, velocity, and separation). In particular, the mass of the companion star of AFGL 3068 is estimated to be greater than 2.6 solar mass. This method is applicable to other AGB stars, providing a step toward understanding the role binary stars can play in explaining the diverse patterns in observed circumstellar envelopes.

The truncation of an optically thick, geometrically thin accretion disk is investigated in the context of low luminosity AGN (LLAGN). We generalize the disk evaporation model used in the interpretative framework of black hole X-ray binaries by including the effect of a magnetic field in accretion disks surrounding supermassive black holes. The critical transition mass accretion rate for which the disk is truncated is found to be insensitive to magnetic effects, but its inclusion leads to a smaller truncation radius in comparison to a model without its consideration. That is, a thin viscous disk is truncated for LLAGN at an Eddington ratio less than 0.03 for a standard viscosity parameter ($\alpha = 0.3$). An increase of the viscosity parameter results in a higher critical transition mass accretion rate and a correspondingly smaller truncation distance, the latter accentuated by greater magnetic energy densities in the disk. Based on these results, the truncation radii inferred from spectral fits of LLAGN published in the literature are consistent with the disk evaporation model. The infrared emission arising from the truncated geometrically thin accretion disks may be responsible for the red bump seen in such LLAGN.

Observations of increasingly higher spatial resolution reveal the existence of asymmetries in the circumstellar envelopes of a small fraction of asymptotic giant branch (AGB) stars. Although there is no general consensus for their origin, a binary companion star may be responsible. Within this framework, we investigate the gravitational effects associated with a sufficiently wide binary system, where Roche lobe overflow is unimportant, on the outflowing envelopes of AGB stars using three dimensional hydrodynamic simulations. The effects due to individual binary components are separately studied, enabling investigation of the stellar and circumstellar characteristics in detail. The reflex motion of the AGB star alters the wind velocity distribution, thereby, determining the overall shape of the outflowing envelope. On the other hand, the interaction of the companion with the envelope produces a gravitational wake, which exhibits a vertically thinner shape. The two patterns overlap and form clumpy structures. To illustrate the diversity of shapes, we present the numerical results as a function of inclination angle. Not only is spiral structure produced by the binary interaction, but arc patterns are also found that represent the former structure when viewed at different inclinations. The arcs reveal a systematic shift of their centers of curvature for cases when the orbital speed of the AGB star is comparable to its wind speed. They take on the shape of a peanut for inclinations nearly edge-on. In the limit of slow orbital motion of the AGB star relative to the wind speed, the arc pattern becomes nearly spherically symmetric. We find that the aspect ratio of the overall oblate shape of the pattern is an important diagnostic probe of the binary as it can be used to constrain the orbital velocity of the AGB star, and moreover the binary mass ratio.

Type Ia supernovae are thought to be caused by thermonuclear explosions of a carbon-oxygen white dwarf in close binary systems. In the single-degenerate scenario (SDS), the companion star is non-degenerate and can be significantly affected by the explosion. We explore this interaction by means of multi-dimensional adaptive mesh refinement simulations using the FLASH code. We consider several different companion types, including main-sequence-like stars (MS), red giants (RG), and helium stars (He). In addition, we include the symmetry-breaking effects of orbital motion, rotation of the non-degenerate star, and Roche-lobe overflow. A detailed study of a sub-grid model for Type Ia supernovae is also presented. We find that the dependence of the unbound stellar mass and kick velocity on the initial binary separation can be fitted by power-law relations. By using the tracer particles in FLASH, the process leading to the unbinding of matter is dominated by ablation, which has usually been neglected in past analytical studies. The level of Ni/Fe contamination of the companion that results from the passage of the supernova ejecta is found to be a $\sim 10^{-5} M_{\odot}$ for the MS star, $\sim 10^{-4} M_{\odot}$ for the He star, and $\sim 10^{-8} M_{\odot}$ for the RG. The spinning MS companion star loses about half of its initial angular momentum during the impact, causing the rotational velocity to drop to a quarter of the original rotational velocity, suggesting that the Tycho G star is a promising progenitor candidate in the SDS.

The hydrodynamic evolution of the common envelope phase of a low mass binary composed of a 1.05 Msun red giant and a 0.6 Msun companion has been followed for five orbits of the system using a high resolution method in three spatial dimensions. During the rapid inspiral phase, the interaction of the companion with the red giant's extended atmosphere causes about 25% of the common envelope to be ejected from the system, with mass continuing to be lost at the end of the simulation at a rate ~ 2 Msun/yr. In the process the resulting loss of angular momentum and energy reduces the orbital separation by a factor of seven. After this inspiral phase the eccentricity of the orbit rapidly decreases with time. The gravitational drag dominates hydrodynamic drag at all times in the evolution, and the commonly-used Bondi-Hoyle-Lyttleton prescription for estimating the accretion rate onto the companion significantly overestimates the true rate. On scales comparable to the orbital separation, the gas flow in the orbital plane in the vicinity of the two cores is subsonic with the gas nearly corotating with the red giant core and circulating about the red giant companion. On larger scales, 90% of the outflow is contained within 30 degrees of the orbital plane, and the spiral shocks in this material leave an imprint on the density and velocity structure. Of the energy released by the inspiral of the cores, only about 25% goes toward ejection of the envelope.

The $Fermi$-LAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the Galactic $\gamma$-ray source population. Such pulsars may also be important in modeling the quiescent state of several low mass X-ray binaries (LMXBs), where optical observations of the companion star suggest the possible existence of rotation powered MSPs. To understand the observational properties of the different evolutionary stages of MSPs, the X-ray and $\gamma$-ray emission associated with the outer gap model is investigated. For rotation powered MSPs, the size of the outer gap and the properties of the high-energy emission are controlled by either the photon-photon pair-creation process or magnetic pair-creation process near the surface. For these pulsars, we find that the outer gap model controlled by the magnetic pair-creation process is preferable in explaining the possible correlations between the $\gamma$-ray luminosity or non-thermal X-ray luminosity versus the spin down power. For the accreting MSPs in quiescent LMXBs, the thermal X-ray emission at the neutron star surface resulting from deep crustal heating can control the conditions in the outer gap. We argue that the optical modulation observed in the quiescent state of several LMXBs originates from the irradiation of the donor star by $\gamma$-rays from the outer gap. In these systems, the irradiation luminosity required for the optical modulation of the source such as SAX J1808.4-3658 can be achieved for a neutron star of high mass. Finally, we discuss the high-energy emission associated with an intra-binary shock in black widow systems, e.g. PSR B1957+20.

Recent observations of strikingly well-defined spirals in the circumstellar envelopes of asymptotic giant branch (AGB) stars point to the existence of binary companions in these objects. In the case of planet or brown dwarf mass companions, we investigate the observational properties of the spiral-onion shell wakes due to the gravitational interaction of these companions with the outflowing circumstellar matter. Three dimensional hydrodynamical simulations at high resolution show that the substellar mass objects produce detectable signatures at 100 AU distance, for the wake induced by a Jupiter to brown dwarf mass object orbiting a solar mass AGB star. In particular, the arm pattern propagates with a speed depending on the local wind and sound speeds, implying possible variations in the arm separation in the wind acceleration region and/or in a slow wind with significant temperature variation. The pattern propagation speeds of the inner and outer boundaries differ by twice the sound speed, leading to the overlap of high density boundaries in slow winds and producing a subpattern of the spiral arm feature. Vertically, the wake forms concentric arcs with angular sizes anticorrelated to the wind Mach number. We provide an empirical formula for the peak density enhancement as a function of the mass, orbital distance, and velocity of the object as well as the wind and local sound speeds. In typical condition of AGB envelopes, the arm-interarm density contrast can be greater than 30 % of the background density within a distance of ~10(M_p/M_J) AU for the object mass M_p in units of Jupiter mass M_J. These results suggest that such features may probe unseen substellar mass objects embedded in the winds of AGB stars and may be useful in planning future high sensitivity/resolution observations with ALMA.

Taking into account a range of parameters determined from the evolutionary models and available observational data, the detailed non-LTE spectra for the primary star and the irradiated donor star in the only known eclipsing AM CVn system to date SDSS J0926+3624 are constructed based on the TLUSTY stellar atmosphere code. The combined spectrum of the primary and the donor stars along with a multi-color blackbody spectrum of the accretion disk that reproduces a detailed numerical model is compared to the SDSS optical spectrum of the system. The photometric flux of the primary star inferred from eclipse observations is compared with the synthetic spectrum. The model fit of the two independent observations provides an upper limit on the distance of the system for different effective temperatures of the primary. In addition, an upper limit on the combined flux of the disk and the donor in the infrared region wherein the contribution of the primary is negligible is also determined. It is shown that the spectrum of a sufficiently cool donor can exhibit emission lines due to irradiation from a hot primary and the emission features should be detectable in the infrared even though the contribution of the flux from the disk dominates. Thus, it is pointed out that infrared observations of the system would provide important information on the thermal state of the donor as well as useful insight on the thermal properties of the primary star and the accretion disk.

NGC 4945 is a Seyfert 2 galaxy at a distance of 3.82 Mpc. Its relative proximity has permitted a detailed SMA study of the circumnuclear molecular gas in a galaxy exhibiting an AGN. Based on an analysis of the high-resolution velocity field of the central region (20" X 20", 1" = 19 pc), we demonstrate that the S-shaped structure of the isovelocity contours is well reproduced by the numerical results of a two dimensional hydrodynamical simulation. In particular, the velocity structure is represented by the bending produced by a shock along the spiral density waves, which are excited at the outer-inner Lindblad resonance by a fast rotating bar. The simulated density map reveals a pair of tightly wound spirals in the center which pass through most of the ring-like (claimed to be a circumnuclear starburst ring by other authors) high intensity region in the observations as well as intersect several Pa$\alpha$ emission line knots located outside the ring-like region. The calculated mass inflow rate at a scale of 50 pc is about three times the inferred mass accretion rate of the AGN of NGC 4945. We find that self-gravity of the gas is important and should be included in our model for NGC 4945. The model is compared with the gas orbit model discussed in Lim et al. (2009), and it is shown that the hydrodynamic model provides a better match to the observed position-velocity diagram and, hence, provides a more reliable prediction of the outer inner Lindblad resonance position.

Beyond the main sequence solar type stars undergo extensive mass loss, providing an environment where planet and brown dwarf companions interact with the surrounding material. To examine the interaction of substellar mass objects embedded in the stellar wind of an asymptotic giant branch (AGB) star, three dimensional hydrodynamical simulations at high resolution have been calculated utilizing the FLASH adaptive mesh refinement code. Attention is focused on the perturbation of the substellar mass objects on the morphology of the outflowing circumstellar matter. In particular, we determine the properties of the resulting spiral density wake as a function of the mass, orbital distance, and velocity of the object as well as the wind velocity and its sound velocity. Our results suggest that future observations of the spiral pattern may place important constraints on the properties of the unseen low mass companion in the outflowing stellar wind.

Thermal conduction between a cool accretion disk and a hot inner corona can result in either evaporation of the disk or condensation of the hot corona. At low mass accretion rates, evaporation dominates and can completely remove the inner disk. At higher mass accretion rates, condensation becomes more efficient in the very inner regions, so that part of the mass accretes via a weak (initially formed) inner disk which is separated from the outer disk by a fully evaporated region at mid radii. At still higher mass accretion rates, condensation dominates everywhere, so there is a continuous cool disk extending to the innermost stable circular orbit. We extend these calculations by including the effect of irradiation by the hot corona on the disk structure. The flux which is not reflected is reprocessed in the disk, adding to the intrinsic thermal emission from gravitational energy release. This increases the seed photons for Compton cooling of the hot corona, enhancing condensation of the hot flow and re-inforcing the residual inner disk rather than evaporating it. Our calculations confirm that a residual inner disk can co-exist with a hard, coronally dominated, spectrum over a range of $0.006<\dot m<0.016$ (for $\alpha=0.2$). This provides an explanation for the weak thermal component seen recently in the low/hard state of black hole X-ray binary systems.

The low-mass X-ray binary (LMXB) system FIRST J102347.6+003841 hosts a newly born millisecond pulsar (MSP) PSR J1023+0038 that was revealed as the first and only known rotation-powered MSP in a quiescent LMXB. While the system is shown to have an accretion disk before 2002, it remains unclear how the accretion disk has been removed in order to reveal the radio pulsation in 2007. In this Letter, we report the discovery of gamma-rays spatially consistent with FIRST J102347.6+003841, at a significance of 7 standard deviations, using data obtained by the Fermi Gamma-ray Space Telescope. The gamma-ray spectrum can be described by a power law (PL) with a photon index of 2.9+-0.2, resulting in an energy flux above 200 MeV of (5.5+-0.9)x10^{-12} erg cm^{-2}s^{-1}. The gamma-rays likely originate from the MSP PSR J1023+0038, but also possibly from an intrabinary shock between the pulsar and its companion star. To complement the gamma-ray study, we also re-investigate the XMM-Newton data taken in 2004 and 2008. Our X-ray spectral analysis suggests that a broken PL with two distinct photon indices describes the X-ray data significantly better than a single PL. This indicates that there exists two components and that both components appear to vary with the orbital phase. The evidence for gamma-ray emission conforms with a recent suggestion that gamma-rays from PSR J1023+0038 may be responsible for ejecting the disk material out of the system.

We present adaptive mesh refinement (AMR) hydrodynamical simulations of the interaction between Type Ia supernovae and their companion stars within the context of the single-degenerate model. Results for 3D red-giant companions without binary evolution agree with previous 2D results by Marietta et al. We also consider evolved helium-star companions in 2D. For a range of helium-star masses and initial binary separations, we examine the mass unbound by the interaction and the kick velocity delivered to the companion star. We find that unbound mass versus separation obeys a power law with index between -3.1 and -4.0, consistent with previous results for hydrogen-rich companions. Kick velocity also obeys a power-law relationship with binary separation, but the slope differs from those found for hydrogen-rich companions. Assuming accretion via Roche-lobe overflow, we find that the unbound helium mass is consistent with observational limits. Ablation (shock heating) appears to be more important in removing gas from helium-star companions than from hydrogen-rich ones, though stripping (momentum transfer) dominates in both cases.

We present basic properties of primary stars that initiate a common envelope (CE) in a binary, while on the giant branch. We use the population-synthesis code described in Politano et al. (2010) and follow the evolution of a population of binary stars up to the point where the primary fills its Roche lobe and initiates a CE. We then collect the properties of each system, in particular the donor mass and the binding energy of the donor's envelope, which are important for the treatment of a CE. We find that for most CEs, the donor mass is sufficiently low to define the core-envelope boundary reasonably well. We compute the envelope-structure parameter {\lambda_\mathrm{env}} from the binding energy and compare its distribution to typical assumptions that are made in population-synthesis codes. We conclude that {\lambda_\mathrm{env}} varies appreciably and that the assumption of a constant value for this parameter results in typical errors of 20--50%. In addition, such an assumption may well result in the implicit assumption of unintended and/or unphysical values for the CE parameter {\alpha_\mathrm{CE}}. Finally, we discuss accurate existing analytic fits for the envelope binding energy, which make these oversimplified assumptions for {\lambda_\mathrm{env}}, and the use of {\lambda_\mathrm{env}} in general, unnecessary.

The heating associated with the deposition of $\gamma$-rays in an accretion disk is proposed as a mechanism to facilitate the transformation of a low mass X-ray binary to the radio millisecond pulsar phase. The $\gamma$-ray emission produced in the outer gap accelerator in the pulsar magnetosphere likely irradiates the surrounding disk, resulting in its heating and to the possible escape of matter from the system. We apply the model to PSR J1023+0038, which has recently been discovered as a newly born rotation powered millisecond pulsar. The predicted $\gamma$-ray luminosity $\sim 6 \times 10^{34}~\mathrm{erg~s^{-1}}$ can be sufficient to explain the disappearance of the truncated disk existing during the 8~month$\sim 2$~yr period prior to the 2002 observations of J1023+0038 and the energy input required for the anomalously bright optical emission of its companion star.

The nature of vertical heating of disk stars in the inner as well as the outer region of disk galaxies is studied. The galactic bar (which is the strongest non-axisymmetric pattern in the disk) is shown to be a potential source of vertical heating of the disk stars in the inner region. Using a nearly self-consistent high-resolution N-body simulation of disk galaxies, the growth rate of the bar potential is found to be positively correlated with the vertical heating exponent in the inner region of galaxies. We also characterize the vertical heating in the outer region where the disk dynamics is often dominated by the presence of transient spiral waves and mild bending waves. Our simulation results suggest that the non-axisymmetric structures are capable of producing the anisotropic heating of the disk stars.

Using a population synthesis technique, we have calculated detailed models of the present-day field population of objects that have resulted from the merger of a giant primary and a main-sequence or brown dwarf secondary during common-envelope evolution. We used a grid of 116 stellar and 32 low-mass/brown dwarf models, a crude model of the merger process, and followed the angular momentum evolution of the binary orbit and the primary's rotation prior to merger, as well as the merged object's rotation after the merger. We find that present-day merged objects that are observable as giant stars or core-helium burning stars in our model population constitute between 0.24% and 0.33% of the initial population of ZAMS binaries, depending upon the input parameters chosen. The median projected rotational velocity of these merged objects is ~16 km/sec, an order of magnitude higher than the median projected rotational velocity in a model population of normal single stars calculated using the same stellar models and initial mass function. The masses of the merged objects are typically less than ~2 solar masses, with a median mass of 1.28 solar masses, which is slightly more than, but not significantly different from, their normal single star counterparts. The luminosities in our merged object population range from ~10-100 solar luminosities, with a strong peak in the luminosity distribution at ~60 solar luminosities, since the majority of the merged objects (57%) lie on the horizontal branch at the present epoch. The results of our population synthesis study are discussed in terms of possible observational counterparts either directly involving the high rotational velocity of the merger product or indirectly, via the effect of rotation on envelope abundances and on the amount and distribution of circumstellar matter.

We have developed a new method for calculating common envelope (CE) events based on explicit consideration of the donor star's structural response to adiabatic mass loss. In contrast to existing CE prescriptions, which specify a priori the donor's remnant mass, we determine this quantity self-consistently and find it depends on binary and CE parameters. This aspect of our model is particularly important to realistic modeling for upper main sequence star donors without strongly degenerate cores (and hence without a clear core/envelope boundary). We illustrate the central features of our method by considering CE events involving 10 solar mass donors on or before their red giant branch. For such donors, the remnant core mass can be as much as 30% larger than the star's He-core mass. Applied across a population of such binaries, our methodology results in a significantly broader remnant mass and final orbital separation distribution and a 20% increase in CE survival rates as compared to previous prescriptions for the CE phase.

The impact of Type Ia supernova ejecta on a helium-star companion is investigated via high-resolution, two-dimensional hydrodynamic simulations. For a range of helium-star models and initial binary separations it is found that the mass unbound in the interaction, $\delta M_{\rm ub}$, is related to the initial binary separation, $a$, by a power law of the form $\delta M_{\rm ub} \propto a^{m}$. This power-law index is found to vary from -3.1 to -4.0, depending on the mass of the helium star. The small range of this index brackets values found previously for hydrogen-rich companions, suggesting that the dependence of the unbound mass on orbital separation is not strongly sensitive to the nature of the binary companion. The kick velocity is also related to the initial binary separation by a power law with an index in a range from -2.7 to -3.3, but the power-law index differs from those found in previous studies for hydrogen-rich companions. The space motion of the companion after the supernova is dominated by its orbital velocity in the pre-supernova binary system. The level of Ni/Fe contamination of the companion resulting from the passage of the supernova ejecta is difficult to estimate, but an upper limit on the mass of bound nickel is found to be $\sim 5\times 10^{-4}\ M_\odot$.

The cumulative luminosity distribution functions (CLFs) of radio millisecond pulsars (MSPs) in globular clusters (GCs) and in the Galactic field at a frequency of 1.4 GHz have been examined. Assuming a functional form, $N \propto L^q$ where $N$ is the number of MSPs and $L$ is the luminosity at 1.4 GHz, it is found that the CLFs significantly differ with a steeper slope, $q=-0.83 \pm 0.05$, in GCs than in the Galactic field ($q=-0.48 \pm 0.04$), suggesting a different formation or evolutionary history of MSPs in these two regions of the Galaxy. To probe the production mechanism of MSPs in clusters, a search of the possible relationships between the MSP population and cluster properties was carried out. The results of an investigation of 9 GCs indicate positive correlations between the MSP population and the stellar encounter rate and metallicity. This provides additional evidence suggesting that stellar dynamical interactions are important in the formation of the MSP population in GCs.

The structure of molecular cloud cores supported by thermal pressure and a poloidal magnetic field is reinvestigated in the magnetohydrostatic and axisymmetric approximation. In addition to oblate configurations found in earlier work, solutions yielding prolate spheroidal shapes have also been obtained for a reference state described by a uniform sphere threaded by a uniform background magnetic field. The solutions for prolate configurations are found to be relevant for lower masses than for their oblate counterparts. Of particular importance is the result that the prolate cloud cores have radii less than a maximum given by $0.25 pc (a/0.2km/s)^2 (P_{\rm ext}/10^{-12} dyne/cm^2)^{-1/2}$, where $a$ is the sound speed and $P_{\rm ext}$ is the external pressure of the background medium. The existence of such solutions obviates the presence of toroidal fields in such modeled structures.

The disk corona evaporation model extensively developed for the interpretation of observational features of black hole X-ray binaries (BHXRBs) is applied to AGNs. Since the evaporation of gas in the disk can lead to its truncation for accretion rates less than a maximal evaporation rate, the model can naturally account for the soft spectrum in high luminosity AGNs and the hard spectrum in low luminosity AGNs. The existence of two different luminosity levels describing transitions from the soft to hard state and from the hard to soft state in BHXRBs, when applied to AGNs, suggests that AGNs can be in either spectral state within a range of luminosities. For example, at a viscosity parameter, \alpha, equal to 0.3, the Eddington ratio from the hard to soft transition and from the soft to hard transition occurs at 0.027 and 0.005 respectively. When the Eddington ratio of the AGN lies below the critical value corresponding to its evolutionary state, the disk is truncated. With decreasing Eddington ratios, the inner edge of the disk increases to greater distances from the black hole with a concomitant increase in the inner radius of the broad line region, $R_{BLR}$. The absence of an optically thick inner disk at low luminosities gives rise to region in the size of borad line-luminosity plane for which the relation $R_{BLR} \propto L^{1/2}$ inferred at high luminosities is excluded. As a result, a lower limit to the accretion rate is predicted for the observability of broad emission lines, if the broad line region is associated with an optically thick accretion disk. Thus, true Seyfert 2 galaxies may exist at very low accretion rates/luminosities. The differences between BHXRBs and AGNs in the framework of the disk corona model are discussed and possible modifications to the model are briefly suggested.

The spectral energy distribution from the X-ray to the very high energy regime ($>100$ GeV) has been investigated for the $\gamma$-ray binary system PSR B1259-63/SS2883 as a function of orbital phase within the framework of a simple model of a pulsar wind nebula. The emission model is based on the synchrotron radiation process for the X-ray regime and the inverse Compton scattering process boosting stellar photons from the Be star companion to the very high energy (100GeV-TeV) regime. With this model, the observed temporal behavior can, in principle, be used to probe the pulsar wind properties at the shock as a function of the orbital phase. Due to theoretical uncertainties in the detailed microphysics of the acceleration process and the conversion of magnetic energy into particle kinetic energy, the observed X-ray data for the entire orbit are fit using two different methods.

The unresolved X-ray emission in the cores of 10 globular clusters hosting millisecond pulsars is investigated. Subtraction of the known resolved point sources leads to detectable levels of unresolved emission in the core region of M28, NGC 6440, M62, and NGC 6752. The X-ray luminosities in the 0.3-8 keV energy band of this emission component were found to lie in the range $\sim 1.5 \times 10^{31}$ ergs s$^{-1}$ (NGC 6752) to $\sim 2.2 \times 10^{32}$ ergs s$^{-1}$ (M28). The lowest limiting luminosity for X-ray source detections amongst these four clusters was $1.1 \times 10^{30}$ ergs s$^{-1}$ for NGC 6752. The spectrum of the unresolved emission can be fit equally well by a power-law, a thermal bremsstrahlung model, a black body plus power-law, or a thermal bremsstrahlung model plus black body component. The unresolved emission is considered to arise from the cumulative contribution of active binaries, cataclysmic variables, and faint millisecond pulsars with their associated pulsar wind nebulae. In examining the available X-ray data, no evidence for any pulsar wind nebular emission in globular clusters is found. It is shown that the X-ray luminosity contribution of a faint source population based on an extrapolation of the luminosity function of detected point sources is compatible with the unresolved X-ray emission in the cores of NGC 6440 and NGC 6752. Adopting the same slope for the luminosity function for M62 as for NGC 6440 and NGC 6752 leads to a similar result for M62. For M28, the contribution from faint sources in the core can attain a level comparable with the observed value if a steeper slope is adopted. The characteristics on the faint source population as constrained by the properties of the unresolved X-ray emission are briefly discussed.

We observed SAX J1808.4-3658 (1808), the first accreting millisecond pulsar, in deep quiescence with XMM-Newton and (near-simultaneously) Gemini-South. The X-ray spectrum of 1808 is similar to that observed in quiescence in 2001 and 2006, describable by an absorbed power-law with photon index 1.74+-0.11 and unabsorbed X-ray luminosity L_X=7.9+-0.7*10^{31} ergs/s, for N_H=1.3*10^{21} cm^{-2}. Fitting all the quiescent XMM-Newton X-ray spectra with a power-law, we constrain any thermally emitting neutron star with a hydrogen atmosphere to have a temperature less than 30 eV and L_{NS}(0.01-10 keV)<6.2*10^{30} ergs/s. A thermal plasma model also gives an acceptable fit to the continuum. Adding a neutron star component to the plasma model produces less stringent constraints on the neutron star; a temperature of 36^{+4}_{-8} eV and L_{NS}(0.01-10 keV)=1.3^{+0.6}_{-0.8}*10^{31} ergs/s. In the framework of the current theory of neutron star heating and cooling, the constraints on the thermal luminosity of 1808 and 1H 1905+000 require strongly enhanced cooling in the cores of these neutron stars. We compile data from the literature on the mass transfer rates and quiescent thermal flux of the largest possible sample of transient neutron star LMXBs. We identify a thermal component in the quiescent spectrum of the accreting millisecond pulsar IGR J00291+5934, which is consistent with the standard cooling model. The contrast between the cooling rates of IGR J00291+5934 and 1808 suggests that 1808 may have a significantly larger mass. This can be interpreted as arising from differences in the binary evolution history or initial neutron star mass in these otherwise similar systems.

We observed the accreting millisecond pulsar SAX J1808.4-3658 with Gemini-South in g' and i' bands, nearly simultaneous with XMM-Newton observations. A clear periodic flux modulation on the system's orbital period is present, consistent with the varying aspect of the donor star's heated face. We model the contributions of a disk and donor star to these optical bands. To produce the observed modulation amplitudes, we conclude that the donor must be irradiated by an external flux 2 orders of magnitude greater than provided by the measured X-ray luminosity. A possible explanation for this irradiation is that the radio pulsar mechanism becomes active during the quiescent state as suggested by Burderi et al., with relativistic particles heating the donor's day-side face. Our modelling constrains the binary inclination to be between 36 and 67 degrees. We obtain estimates for the pulsar mass of >2.2 solar masses (although this limit is sensitive to the source's distance), consistent with the accelerated NS cooling in this system indicated by X-ray observations. We also estimate the donor mass to be in the range of 0.07-0.11 solar masses, providing further indications that the system underwent non-standard binary evolution to reach its current state.

We quantify an evolutionary channel for single sdB stars based on mergers of binaries containing a red giant star and a lower mass main sequence or brown dwarf companion in our Galaxy. Population synthesis calculations that follow mergers during the common envelope phase of evolution of such systems reveal a population of rapidly rotating horizontal branch stars with a distribution of core masses between 0.32 Mo - 0.7 Mo that is strongly peaked between 0.47 Mo - 0.54 Mo. The high rotation rates in these stars are a natural consequence of the orbital angular momentum deposition during the merger and the subsequent stellar contraction of the merged object from the tip of the red giant branch. We suggest that centrifugally enhanced mass loss facilitated by the rapid rotation of these stars may lead to the formation of single sdB stars for some of these objects.

The formation of a cool disk in the innermost regions of black hole X-ray transient systems in the low hard state is investigated. Taking into account the combined cooling associated with the Compton and conductive energy transport processes in a corona, the radial structure of a disk is described for a range of mass accretion rates. The mass flow in an optically thick inner region can be maintained by the condensation of matter from a corona with the disk temperature and luminosity varying continuously as a function of the accretion rate. Although such a disk component can be present, the contribution of the optically thick disk component to the total luminosity can be small since the mass flow due to condensation in the optically thick disk underlying the corona can be significantly less than the mass flow rate in the corona. The model is applied to the observations of the low quiescent state of the black hole source GX 339-4 at luminosities of around $0.01 L_{Edd}$ and is able to explain the temperature of the thermal component at the observed luminosities. Since conductive cooling dominates Compton cooling at low mass accretion rates, the luminosity corresponding to the critical mass accretion rate above which a weak thermal disk component can be present in the low hard state is estimated to be as low as $0.001 L_{Edd}$.

The evolution leading to the formation of a neutron star in the very young Westerlund 1 star cluster is investigated. The turnoff mass has been estimated to be 35 Msun, indicating a cluster age ~ 3-5 Myr. The brightest X-ray source in the cluster, CXO J164710.2-455216, is a slowly spinning (10 s) single neutron star and potentially a magnetar. Since this source was argued to be a member of the cluster, the neutron star progenitor must have been very massive (M_zams > 40 Msun) as noted by Muno et al. (2006). Since such massive stars are generally believed to form black holes (rather than neutron stars), the existence of this object poses a challenge for understanding massive star evolution. We point out while single star progenitors below M_zams < 20 Msun form neutron stars, binary evolution completely changes the progenitor mass range. In particular, we demonstrate that mass loss in Roche lobe overflow enables stars as massive as 50-80 Msun, under favorable conditions, to form neutron stars. If the very high observed binary fraction of massive stars in Westerlund 1 (> 70 percent) is considered, it is natural that CXO J164710.2-455216 was formed in a binary which was disrupted in a supernova explosion such that it is now found as a single neutron star. Hence, the existence of a neutron star in a given stellar population does not necessarily place stringent constraints on progenitor mass when binary interactions are considered. It is concluded that the existence of a neutron star in Westerlund 1 cluster is fully consistent with the generally accepted framework of stellar evolution.

Mergers of double neutron stars are considered the most likely progenitors for short gamma-ray bursts. Indeed such a merger can produce a black hole with a transient accreting torus of nuclear matter (Lee & Ramirez-Ruiz 2007, Oechslin & Janka 2006), and the conversion of a fraction of the torus mass-energy to radiation can power a gamma-ray burst (Nakar 2006). Using available binary pulsar observations supported by our extensive evolutionary calculations of double neutron star formation, we demonstrate that the fraction of mergers that can form a black hole -- torus system depends very sensitively on the (largely unknown) maximum neutron star mass. We show that the available observations and models put a very stringent constraint on this maximum mass under the assumption that a black hole formation is required to produce a short gamma-ray burst in a double neutron star merger. Specifically, we find that the maximum neutron star mass must be within 2 - 2.5 Msun. Moreover, a single unambiguous measurement of a neutron star mass above 2.5 Msun would exclude a black hole -- torus central engine model of short gamma-ray bursts in double neutron star mergers. Such an observation would also indicate that if in fact short gamma-ray bursts are connected to neutron star mergers, the gamma-ray burst engine is best explained by the lesser known model invoking a highly magnetized massive neutron star (e.g., Usov 1992; Kluzniak & Ruderman 1998; Dai et al. 2006; Metzger, Quataert & Thompson 2007).

The evolution of the spin and tilt of black holes in compact black hole - neutron star and black hole - black hole binary systems is investigated within the framework of the coalescing compact star binary model for short gamma ray bursts via the population synthesis method. Based on recent results on accretion at super critical rates in slim disk models, estimates of natal kicks, and the results regarding fallback in supernova models, we obtain the black hole spin and misalignment. It is found that the spin parameter, a_spin}, is less than 0.5 for initially non rotating black holes and the tilt angle, i_tilt, is less than 45 deg for 50% of the systems in black hole - neutron star binaries. Upon comparison with the results of black hole - neutron star merger calculations we estimate that only a small fraction (~ 0.01) of these systems can lead to the formation of a torus surrounding the coalesced binary potentially producing a short-hard gamma ray burst. On the other hand, for high initial black hole spin parameters (a_spin>0.6) this fraction can be significant (~ 0.4). It is found that the predicted gravitational radiation signal for our simulated population does not significantly differ from that for non rotating black holes. Due to the (i) insensitivity of signal detection techniques to the black hole spin and the (ii) predicted overall low contribution of black hole binaries to the signal we find that the detection of gravitational waves are not greatly inhibited by current searches with non spinning templates. It is pointed out that the detection of a black hole - black hole binary inspiral system with LIGO or VIRGO may provide a direct measurement of the initial spin of a black hole.

We use high-resolution, three-dimensional hydrodynamic simulations to study the hydrodynamic and gravitational interaction between stellar companions embedded within a differentially rotating common envelope. Specifically, we evaluate the contributions of the nonaxisymmetric gravitational tides and ram pressure forces to the drag force and, hence, to the dissipation rate and the mass accumulated onto the stellar companion. We find that the gravitational drag dominates the hydrodynamic drag during the inspiral phase, leading to the result that a simple prescription based on a gravitational capture radius formalism significantly underestimates the dissipation rate and overestimates the inspiral decay timescale. Although the rate of mass accretion fluctuates significantly, we observe a secular trend leading to an effective rate of mass accretion which is significantly less than the rate based on a gravitational capture radius. The implications of these results are discussed within the context of accretion of compact objects in the common-envelope phase.

The accreting millisecond pulsar SAX J1808.4-3658 may be a transition object between accreting X-ray binaries and millisecond radio pulsars. We have constrained the thermal radiation from its surface through XMM-Newton X-ray observations, providing strong evidence for neutrino cooling processes from the neutron star core. We have also undertaken simultaneous X-ray and optical (Gemini) observations, shedding light on whether the strong heating of the companion star in quiescence may be due to X-ray irradiation, or to a radio pulsar turning on when accretion stops.

The condensation of matter from a corona to a cool, optically thick inner disk is investigated for black hole X-ray transient systems in the low hard state. A description of a simple model for the exchange of energy and mass between corona and disk originating from thermal conduction is presented, taking into account the effect of Compton cooling of the corona by photons from the underlying disk. It is found that a weak, condensation-fed inner disk can be present in the low hard state of black hole transient systems for a range of luminosities which depend on the magnitude of the viscosity parameter. For $\alpha \sim 0.1-0.4$ an inner disk can exist for luminosities in the range $\sim 0.001- 0.02$ Eddington value. The model is applied to the X-ray observations of the black hole candidate sources GX 339-4 and Swift J1753.5-0127 in their low hard state. It is found that Compton cooling is important in the condensation process, leading to the maintenance of cool inner disks in both systems. As the results of the evaporation/condensation model are independent of the black hole mass, it is suggested that such inner cool disks may contribute to the optical and ultraviolet emission of low luminosity active galactic nuclei.

(Abridged) We calculate the full stellar-structural evolution of donors in AM CVn systems formed through the WD channel coupled to the binary's evolution. Contrary to assumptions made in prior modelling, these donors are not fully convective over much of the AM CVn phase and do not evolve adiabatically under mass loss indefinitely. Instead, we identify three distinct phases of evolution: a mass transfer turn-on phase (during which the orbital period continues to decrease even after contact, the donor contracts, and the mass transfer rate accelerates to its maximum), a phase in which the donor expands adiabatically in response to mass loss, and a cooling phase beginning at orbital periods of approximately 45--55 minutes during which the donor contracts. The physics that determines the behaviour in the first and third phases, both of which are new outcomes of this study, are discussed in some detail. We find the overall duration of the turn-on phase to be between $\sim 10^4$-$10^6$ yrs, significantly longer than prior estimates. We predict the donor's luminosity and effective temperature. During the adiabatic expansion phase (ignoring irradiation effects), the luminosity is approximately $10^{-6}$--$10^{-4} L_\odot$ and the effective temperature is approximately 1000--1800 K. However, the flux generated in the accretion flow dominates the donor's intrinsic light at all times. The impact of irradiation on the donor extends the phase of adiabatic expansion to longer orbital periods and alters the donor's observational characteristics. Irradiated donors during the adiabatic phase can attain a surface luminosity of up to $\approx10^{-2} L_\odot$. We argue that the turn-on and cooling phases both will leave significant imprints on the AM CVn population's orbital period distribution.

We present a comprehensive description of the population synthesis code StarTrack. The original code has been significantly modified and updated. Special emphasis is placed here on processes leading to the formation and further evolution of compact objects (white dwarfs, neutron stars, and black holes). Both single and binary star populations are considered. The code now incorporates detailed calculations of all mass-transfer phases, a full implementation of orbital evolution due to tides, as well as the most recent estimates of magnetic braking. This updated version of StarTrack can be used for a wide variety of problems, with relevance to many current and planned observatories, e.g., studies of X-ray binaries (Chandra, XMM-Newton), gravitational radiation sources (LIGO, LISA), and gamma-ray burst progenitors (HETE-II, Swift). The code has already been used in studies of Galactic and extra-galactic X-ray binary populations, black holes in young star clusters, Type Ia supernova progenitors, and double compact object populations. Here we describe in detail the input physics, we present the code calibration and tests, and we outline our current studies in the context of X-ray binary populations.

The magnetic cataclysmic variable AE Aquarii hosts a rapidly rotating white dwarf which is thought to expel most of the material streaming onto it. Observations of AE Aqr have been obtained in the wavelength range of 5 - 70 microns with the IRS, IRAC, and MIPS instruments on board the Spitzer Space Telescope. The spectral energy distribution reveals a significant excess above the K4V spectrum of the donor star with the flux increasing with wavelength above 12.5 microns. Superposed on the energy distribution are several hydrogen emission lines, identified as Pf alpha and Hu alpha, beta, gamma. The infrared spectrum above 12.5 microns can be interpreted as synchrotron emission from electrons accelerated to a power-law distribution dN=E^{-2.4}dE in expanding clouds with an initial evolution timescale in seconds. However, too many components must then be superposed to explain satisfactorily both the mid-infrared continuum and the observed radio variability. Thermal emission from cold circumbinary material can contribute, but it requires a disk temperature profile intermediate between that produced by local viscous dissipation in the disk and that characteristic of a passively irradiated disk. Future high-time resolution observations spanning the optical to radio regime could shed light on the acceleration process and the subsequent particle evolution.

Double black hole binaries are among the most important sources of gravitational radiation for ground-based detectors such as LIGO or VIRGO. Even if formed with lower efficiency than double neutron star binaries, they could dominate the predicted detection rates, since black holes are more massive than neutron stars and therefore could be detected at greater distances. Here we discuss an evolutionary process that can very significantly limit the formation of close double black hole binaries: the vast majority of their potential progenitors undergo a common envelope (CE) phase while the donor, one of the massive binary components, is evolving through the Hertzsprung gap. Our latest theoretical understanding of the CE process suggests that this will probably lead to a merger, inhibiting double black hole formation. Barring uncertainties in the physics of CE evolution, we use population synthesis calculations, and find that the corresponding reduction in the merger rate of double black holes formed in galactic fields is so great (by ~500) that their contribution to inspiral detection rates for ground-based detectors could become relatively small (1 in 10) compared to double neutron star binaries. A similar process also reduces the merger rates for double neutron stars, by factor of ~5, eliminating most of the previously predicted ultracompact NS-NS systems. Our predicted detection rates for Advanced LIGO are now much lower for double black holes (~2/yr), but are still quite high for double neutron stars (~20/yr). If double black holes were found to be dominant in the detected inspiral signals, this could indicate that they mainly originate from dense star clusters (not included here) or that our theoretical understanding of the CE phase requires significant revision.

Thermal X-ray radiation from neutron star soft X-ray transients in quiescence provides the strongest constraints on the cooling rates of neutron stars, and thus on the interior composition and properties of matter in the cores of neutron stars. We analyze new (2006) and archival (2001) XMM-Newton observations of the accreting millisecond pulsar SAX J1808.4-3658 in quiescence, which provide the most stringent constraints to date. The X-ray spectrum of SAX J1808.4-3658 in the 2006 observation is consistent with a power-law of photon index 1.83\pm0.17, without requiring the presence of a blackbody-like component from a neutron star atmosphere. Our 2006 observation shows a slightly lower 0.5-10 keV X-ray luminosity, at a level of 68^{+15}_{-13}% that inferred from the 2001 observation. Simultaneous fitting of all available XMM data allows a constraint on the quiescent neutron star (0.01-10 keV) luminosity of L_{NS}<1.1*10^{31} erg/s. This limit excludes some current models of neutrino emission mediated by pion condensates, and provides further evidence for additional cooling processes, such as neutrino emission via direct Urca processes involving nucleons and/or hyperons, in the cores of massive neutron stars.

The population of non magnetic cataclysmic variables evolving under the influence of a circumbinary disk is investigated for systems above the upper edge of the period gap at orbital periods greater than 2.75hr. For a fractional mass input rate into the disk, corresponding to 3e-4 of the mass transfer rate, the model systems exhibit a bounce at orbital periods greater than 2.75hr. The simulations reveal that (1) some systems can exist as dwarf nova type systems throughout their lifetime, (2) dwarf nova type systems can evolve into nova-like systems as their mass transfer rate increases with increasing circumbinary disk mass, and (3) nova-like systems can evolve back into dwarf nova systems during their postbounce evolution to longer orbital periods. Among these subclasses, nova-like cataclysmic variables would be the best candidates to search for circumbinary disks at wavelengths greater than 10 micron. The theoretical orbital period distribution of our population synthesis model is in reasonable accord with the combined population of dwarf novae and nova-like systems above the period gap, suggesting the possibility that systems with unevolved donors need not detach and evolve below the period gap as in the disrupted magnetic braking model. The resulting population furthermore reveals the possible presence of systems with small mass ratios (a property of systems exhibiting superhump phenomena at long orbital periods) and a preference of O/Ne/Mg white dwarfs in dwarf nova systems in comparison to nova-like systems. The importance of observational bias in accounting for the differing populations is examined, and it is shown that an understanding of these effects is necessary in order to confront the theoretical distributions with the observed ones in a meaningful manner. (abridged)

The common envelope phase of binary star evolution plays a central role in many evolutionary pathways leading to the formation of compact objects in short period systems. Using three dimensional hydrodynamical computations, we review the major features of this evolutionary phase, focusing on the conditions that lead to the successful ejection of the envelope and, hence, survival of the system as a post common envelope binary. Future hydrodynamical calculations at high spatial resolution are required to delineate the regime in parameter space for which systems survive as compact binary systems from those for which the two components of the system merge into a single rapidly rotating star. Recent algorithmic developments will facilitate the attainment of this goal.

We discuss results of our study on AM CVn binaries formed with donors that never ignited He before contact. For the first time, we treat the donor's in these systems in the context of a full stellar structure evolution theory and find that the binary's evolution can described in terms of 3 phases: contact, adiabatic donor expansion, and late-time donor cooling. Details of the first and third phase are new results from this study and we focus on generally characterizing these two phases. Finally, we present our predictions for the donor's light in these systems.

The Low Mach Number Approximation (LMNA) is applied to 2D hydrodynamical modeling of Type I X-ray bursts on a rectangular patch on the surface of a non-rotating neutron star. Because such phenomena involve decidedly subsonic flows, the timestep increase offered by the LMNA makes routine simulations of these deflagrations feasible in an environment where strong gravity produces significant stratification, while allowing for potentially significant lateral differences in temperature and density. The model is employed to simulate the heating, peak, and initial cooling stages in the deep envelope layers of a burst. During the deflagration, Benard-like cells naturally fill up a vertically expanding convective layer. The Mach number is always less than 0.15 throughout the simulation, thus justifying the low Mach number approximation. While the convective layer is superadiabatic on average, significant fluctuations in adiabaticity occur within it on subconvective timescales. Due to convective layer expansion, significant compositional mixing naturally occurs, but tracer particle penetration through the convective layer boundaries on convective timescales is temporary and spatially limited. Thus, mixing occurs on the relatively slow burst timescale through thermal expansion of the convective layer rather than from mass penetration of the convective layer boundary through particle convection. At the convective layer boundaries where mixing is less efficient, the actual temperature gradient more closely follows the Ledoux criteria.

We report on evolutionary calculations of the onset of mass transfer in AM CVn binaries, treating the donor's evolution in detail. We show that during the early contact phase, while the mass transfer rate, $\Mdot$, is increasing, gravity wave (GW) emission continues to drive the binary to shorter orbital period, $\Porb$. We argue that the phase where $\Mdot > 0$ and $\nudot > 0$ ($\nu = 1/\Porb$) can last between $10^3$ and $10^6$ yrs, significantly longer than previously estimated. These results are applied to RX J0806+1527 ($\Porb = 321$ s) and RX J914+2456 ($\Porb=569$ s), both of which have measured $\nudot > 0$. \emph{Thus, a $\nudot > 0$ does not select between the unipolar inductor and accretion driven models proposed as the source of X-rays in these systems}. For the accretion model, we predict for RX J0806 that $\ddot{\nu} \approx \ee{1.0-1.5}{-28}$ Hz s$^{-2}$ and argue that timing observations can probe $\ddot{\nu}$ at this level with a total $\approx 20$ yr baseline. We also place constraints on each system's initial parameters given current observational data.

In this paper, the first of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We examine the processes leading to the creation of X-ray binaries such as cataclysmic variables and AM CVn systems. Using numerical simulations, we identify the dominant formation channels and we predict the expected numbers and characteristics of detectable systems, emphasizing how the cluster sources differ from the field population. We explore the dependence of formation rates on cluster properties and we explain in particular why the distribution of cataclysmic variables has only a weak dependence on cluster density. We also discuss the frequency of dwarf nova outbursts in globular clusters and their connection with moderately strong white dwarf magnetic fields. We examine the rate of Type Ia supernovae via both single and double degenerate channels in clusters and we argue that those rates may contribute to the total SN Ia rate in elliptical galaxies. Considering coalescing white dwarf binaries we discuss possible constraints on the common envelope evolution of their progenitors and we derive theoretical expectations for gravitational wave detection by LISA.

We report on the emission properties of PSR B1929+10 and its putative X-ray trail from a multi-wavelength study performed with XMM-Newton, the ESO NTT, the HST, the Effelsberg 100m Radio Telescope and the Jodrell Bank Radio Observatory. The XMM-Newton observations confirm the existence of the diffuse emission with a trail morphology lying in a direction opposite to the transverse motion of the pulsar. The trail has a length of ~15 arcmin. Its spectrum is non-thermal and produced by electron-synchrotron emission in the shock between the pulsar wind and the surrounding medium. Assuming that the electron lifetime against synchrotron cooling is comparable to the source transit time over the X-ray trail length, the magnetic field strength in the trail emitting region is inferred to be ~5 uG. Inspecting data from the Effelsberg 11cm radio continuum survey of the Galactic plane we discovered an elongated feature apparently coincident with the X-ray trail. The emission properties observed from PSR 1929+10 are found to finally challenge the picture of the emission properties of old non-recycled rotation powered pulsars. Both, the temporal and spectral X-ray emission properties of PSR 1929+10 are in excellent agreement with a non-thermal and, thus, magnetospheric radiation dominated emission scenario. A flux contribution from the thermal emission of heated polar caps of ~7% is inferred from a best fitting composite Planckian and power law spectral model. The X-ray pulse profile is found to be markedly different from the broad sinusoidal pulse profile seen in the low statistic ROSAT data. Simulations in the framework of an outer-gap emission model are able to reproduce the observed X-ray pulse profile and its phase relation relative to the radio pulse.

The non-thermal, non-pulsed X-ray emission of MSPs is investigated. As in young pulsars, MSPs emit a relativistic wind, which interacting with the interstellar medium and/or a binary companion can significantly contribute to the non-pulsed emission of these pulsars. An application and extension of a simple model developed for young pulsars is applied to the old recycled MSP B1957+20. It is found that the pulsar wind can, indeed, contribute to both the resolved and unresolved X-ray emission. For other MSP in the Galactic field where the spectral index of the non-pulsed component has been measured (i.e., PSR B1937+21, J0218+4232) the contribution of the pulsar wind to the non-pulsed X-ray luminosity is estimated. For the MSPs in the core regions of globular clusters, the pulsar wind nebula is likely affected by its interaction with the dense stellar environment, possibly leading to a diminished contribution to the total X-ray emission. In this case, the existence of non-thermal non-pulsed X-ray emission is more likely for binary rather than isolated MSPs with the emission arising from the interaction of the relativistic pulsar wind with a binary companion. Our study suggests that the magnetization parameter in the pulsar wind nebulae of MSPs is significantly larger than that of the Crab nebula by about a factor of 10. The emission from MSPs moving at high velocities may appear spatially extended with a tail-like morphology, which may contribute to the faint filamentary X-ray source subpopulation in the Galaxy.

The population synthesis of cataclysmic variables below the period is investigated. A grid of detailed binary evolutionary sequences has been calculated and included in the simulations to take account of additional angular momentum losses beyond that associated with gravitational radiation and mass loss, due to nova outbursts, from the system. As a specific example, we consider the effect of a circumbinary disk to gain insight into the ingredients necessary to reproduce the observed orbital period distribution. The resulting distributions show that the period minimum lies at about 80 minutes with the number of systems monotonically increasing with increasing orbital period to a maximum near 90 minutes. There is no evidence for an accumulation of systems at the period minimum which is a common feature of simulations in which only gravitational radiation losses are considered. The period distribution is found to be fairly flat for orbital periods ranging from about 85 to 120 minutes. The steepness of the lower edge of the period gap can be reproduced, for example, by an input of systems at periods near 2.25 hrs due to a flow of cataclysmic variable binary systems from orbital periods longer than 2.75 hrs. The good agreement with the cumulated distribution function of observed systems within the framework of our model indicates that the angular momentum loss by a circumbinary disk or a mechanism which mimics its features coupled with a weighting factor to account for selection effects in the discovery of such systems and a flow of systems from above the period gap to below the period gap are important ingredients for understanding the overall period distribution of cataclysmic variable binary systems.

The effectiveness of the thermal coupling of ions and electrons in the context of optically thin, hot accretion flows is investigated. In the limit of complete coupling, we focus on the one-temperature accretion flows. Based on a global analysis, the results are compared with two-temperature accretion flow models and with the observations of black hole sources. Many features are quite similar. That is, hot one-temperature solutions are found to exist for mass flow rates less than a critical value; i.e., $\dot{M}\la 10\alpha^2\dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}= L_{\rm Edd}/c^2$ is the Eddington accretion rate. At low mass flow rates, $\dot{M}\la 10^{-3}\alpha^2 \dot{M}_{\rm Edd}$, the solution is in the advection-dominated accretion flow (ADAF) regime. But at higher rates, radiative cooling is effective and is mainly balanced by advective {\em heating}, placing the solution in the regime of luminous hot accretion flow (LHAF). To test the viability of the one-temperature models, we have fitted the spectra of the two black hole sources, Sgr A* and XTE J1118+480, which have been examined successfully with two-temperature models. It is found that the one-temperature models do not provide acceptable fits to the multi-wavelength spectra of Sgr A* nor to XTE J1118+480 as a result of the higher temperatures characteristic of the one-temperature models. It is concluded that the thermal coupling of ions and electrons cannot be fully effective and that a two-temperature description is required in hot accretion flow solutions.

Recent deep Chandra surveys of the Galactic center region have revealed the existence of a faint, hard X-ray source population. While the nature of this population is unknown, it is likely that several types of stellar objects contribute. For sources involving binary systems, accreting white dwarfs and accreting neutron stars with main sequence companions have been proposed. Among the accreting neutron star systems, previous studies have focused on stellar wind-fed sources. In this paper, we point out that binary systems in which mass transfer occurs via Roche lobe overflow (RLOF) can also contribute to this X-ray source population. A binary population synthesis study of the Galactic center region has been carried out, and it is found that evolutionary channels for neutron star formation involving the accretion induced collapse of a massive ONeMg white dwarf, in addition to the core collapse of massive stars, can contribute to this population. The RLOF systems would appear as transients with quiescent luminosities, above 2 keV, in the range from 10^31-10^32 ergs/s. The results reveal that RLOF systems primarily contribute to the faint X-ray source population in the Muno et al. (2003) survey and wind-fed systems can contribute to the less sensitive Wang et al. (2002) survey. However, our results suggest that accreting neutron star systems are not likely to be the major contributor to the faint X-ray source population in the Galactic center.

The general properties of the non-thermal non-pulsed X-ray emission of rotation powered pulsars are investigated in the context of a pulsar wind nebula model. An examination of the observed X-ray emission from a sample of 23 pulsars from ASCA in the energy range between 2-10 keV reveals that the relation of X-ray luminosity, $L_x$, to the pulsar spin down power, $\dot{E}$, is steeper for the non-pulsed component than for the pulsed component. Specifically, $L_{x}^{npul} \propto \dot{E}^{1.4 \pm 0.1}$ for the non-pulsed component, whereas $L_{x}^{pul} \propto \dot{E}^{1.2 \pm 0.08}$ for the pulsed component. The former relation is consistent with emission from a pulsar wind nebula model in which $L_{x}^{npul} \propto \dot{E}^{p/2}$ where $p$ is the power law index of the electron energy distribution. The relation for the pulsed component, on the other hand, is consistent with a magnetospheric emission model.

Recent RXTE and Chandra discoveries of low mass X-ray binaries with ultra-short orbital periods have initiated theoretical work on the origins of these peculiar systems. Using the StarTrack population synthesis code the formation and evolution of X-ray ultracompact binaries (UCBs) in the Galactic field are analyzed. The relative number of UCBs with a neutron star or a black hole accretor populating our Galaxy is predicted. Our results demonstrate that standard evolutionary scenarios involving primordial binaries can be sufficient to produce the UCBs in the Galactic field without requiring additional processes associated with the dense stellar environments in the cores of globular clusters. In contrast to previous studies we find that the majority of the immediate progenitors of these systems consist of a hydrogen exhausted donor with an ONeMg white dwarf. The evolution of these systems leads to the accretion induced collapse of the white dwarf to a neutron star, which can play an important role in the formation of a majority of Galactic UCBs. We predict that with an increase in the number of X-ray active UCBs hosting neutron stars by an order of magnitude, a system with a black hole accretor may be found.

We present a search for excess mid-IR emission due to circumbinary material in the orbital plane of cataclysmic variables (CVs). Our motivation stems from the fact that the strong braking exerted by a circumbinary (CB) disc on the binary system could explain several puzzles in our current understanding of CV evolution. Since theoretical estimates predict that the emission from a CB disc can dominate the spectral energy distribution (SED) of the system at wavelengths > 5 microns, we obtained simultaneous visible to mid-IR SEDs for eight systems. We report detections of SS Cyg at 11.7 microns and AE Aqr at 17.6 microns, both in excess of the contribution from the secondary star. In AE Aqr, the IR likely originates from synchrotron-emitting clouds propelled by the white dwarf. In SS Cyg, we argue that the observed mid-IR variability is difficult to reconcile with simple models of CB discs and we consider free-free emission from a wind. In the other systems, our mid-IR upper limits place strong constraints on the maximum temperature of a putative CB disc. The results show that if any sizeable CB disc are present in these systems, they must be self-shadowed or perhaps dust-free, with the peak thermal emission shifted to far-IR wavelengths.

The evolution of binaries consisting of evolved main sequence stars (1 < M_d/Msun < 3.5) with white dwarf companions (0.7 < M_wd/Msun < 1.2) is investigated through the thermal mass transfer phase. Taking into account the stabilizing effect of a strong, optically thick wind from the accreting white dwarf surface, we have explored the formation of several evolutionary groups of systems for progenitors with initial orbital periods of 1 and 2 days. The numerical results show that CO white dwarfs can accrete sufficient mass to evolve to a Type Ia supernova and ONeMg white dwarfs can be built up to undergo accretion induced collapse for donors more massive than about 2 Msun. For donors less massive than ~2 Msun the system can evolve to form a He and CO or ONeMg white dwarf pair. In addition, sufficient helium can be accumulated (~0.1 Msun) in systems characterized by 1.6 < M_d/Msun < 1.9 and 0.8 < M_wd/Msun < 1 such that sub Chandrasekhar mass models for Type Ia supernovae, involving off center helium ignition, are possible for progenitor systems evolving via the Case A mass transfer phase. For systems characterized by mass ratios > 3 the system likely merges as a result of the occurrence of a delayed dynamical mass transfer instability. A semi-analytical model is developed to delineate these phases which can be easily incorporated in population synthesis studies of these systems.

We present a description for the angular momentum loss rate due to magnetic braking for late type stars taking into account recent observational data on the relationship between stellar activity and rotation. The analysis is based on an idealized two component coronal model subject to constraints imposed on the variation of the coronal gas density with rotation period inferred from the observed variation of X-ray luminosity, L_x, with rotation rate, Omega, (L_x \propto Omega^2) for single rotating dwarfs. An application of the model to high rotation rates leads to a gradual turnover of the X-ray luminosity which is similar to the saturation recently observed in rapidly rotating dwarfs. The resulting angular momentum loss rate, \dot J, depends on Omega in the form \dot J \propto Omega^beta where beta ~ 3 for slow rotators and ~ 1.3 for fast rotators. The relation at high rotation rates significantly differs from the power law exponent for slowly rotating stars, depressing the angular momentum loss rate without necessarily requiring the saturation of the magnetic field. The application of this model to the evolution of cataclysmic variable binary systems leads to mass transfer rates that are in approximate accord with those observed in comparison to rates based on either a Skumanich law or an empirical law based on beta=1.

The observed spectra and X-ray luminosities of millisecond pulsars in 47 Tuc can be interpreted in the context of theoretical models based on strong, small scale multipole fields on the neutron star surface. For multipole fields that are relatively strong as compared to the large scale dipole field, the emitted X-rays are thermal and likely result from polar cap heating associated with the return current from the polar gap. On the other hand, for weak multipole fields, the emission is nonthermal and results from synchrotron radiation of $e^{\pm}$ pairs created by curvature radiation. The X-ray luminosity, $L_x$, is related to the spin down power, $L_{sd}$, expressed in the form $L_x \propto L^{\beta}_{sd}$ with $\beta \sim 0.5$ and $\sim 1$ for strong and weak multipole fields respectively. If the polar cap size is of the order of the length scale of the multipole field, $s$ and $\beta \sim 0.5$, the polar cap temperature is $\sim 3 \times 10^6 K (\frac{L_{sd}}{10^{34}erg s^{-1}})^{1/8} (\frac{s}{3\times 10^4 cm})^{-1/2}$. A comparison of the X-ray properties of millisecond pulsars in globular clusters and in the Galactic field suggests that the emergence of relatively strong small scale multipole fields from the neutron star interior may be correlated with the age and evolutionary history of the underlying neutron star.

We have calculated the evolution of 60 model binary systems consisting of helium stars in the mass range of M_He= 2.5-6Msun with a 1.4Msun neutron star companion to investigate the formation of double neutron star systems.Orbital periods ranging from 0.09 to 2 days are considered, corresponding to Roche lobe overflow starting from the helium main sequence to after the ignition of carbon burning in the core. We have also examined the evolution into a common envelope phase via secular instability, delayed dynamical instability, and the consequence of matter filling the neutron star's Roche lobe. The survival of some close He-star neutron-star binaries through the last mass transfer episode (either dynamically stable or unstable mass transfer phase) leads to the formation of extremely short-period double neutron star systems (with P<~0.1days). In addition, we find that systems throughout the entire calculated mass range can evolve into a common envelope phase, depending on the orbital period at the onset of mass transfer. The critical orbital period below which common envelope evolution occurs generally increases with M_He. In addition, a common envelope phase may occur during a short time for systems characterized by orbital periods of 0.1-0.5 days at low He-star masses (~ 2.6-3.3Msun). The existence of a short-period population of double neutron stars increases the predicted detection rate of inspiral events by ground-based gravitational-wave detectors and impacts their merger location in host galaxies and their possible role as gamma-ray burst progenitors. We use a set of population synthesis calculations and investigate the implications of the mass-transfer results for the orbital properties of DNS populations.

The vertical structure of optically thick accretion disks is investigated in the two-zone approximation. The disk is divided into an underlying disk and a corona, where the latter is defined as the upper surface layers for which the effective optical depth is unity. It is found that a significant part of the accretion flow (or dissipation rate) can take place in the corona if the scale height of the magnetic field is larger than that of the disk. The presence of such a dissipative corona leads to a modification in the topology of local disk solutions. For example, these solutions are found from local stability analysis to be both secularly and thermally stable, for accretion rates which are a factor $\approx$ four higher than those inferred from the stability of standard disk solutions. Thus, the applicability of optically thick disks with dissipative coronas are not as restrictive as disks without such coronas and can provide an attractive explanation for the origin of the soft spectral component observed in black hole X-ray binary systems.

We investigate the structure and evolution of a geometrically thin viscous Keplerian circumbinary (CB) disk, using detailed models of their radiative/convective vertical structure. We use a simplified description for the evolution of the cataclysmic binary and focus on cases where the circumbinary disk causes accelerated mass transfer (> 1e-8 Msun/yr). The inner edge of the disk is assumed to be determined by the tidal truncation radius and the mass input rate into the disk is assumed to be a small fraction (1e-5-0.01) of the mass transfer rate. Under the action of the viscous stresses in the disk the matter drifts outward with the optically thick region extending to several AU. The inner part of the disk is cool with maximum effective temperatures < 3,000 K while the outermost parts of the disk are < 30 K and optically thin. We calculate the effects of thermal instability on a sufficiently massive CB disk. It leads to outbursts reminiscent of those in thermally unstable accretion disks, with the instability remaining confined to the inner regions of the CB disk. However, for most of the evolutionary sequences the surface densities required to trigger instability are not reached. The spectral energy distributions from circumbinary disks are calculated, and the prospects for the detection of such disks in the infrared and submm wavelength regions are discussed.

Observations of Cataclysmic Variables show a number of phenomena that do not fit easily into the standard magnetic braking scenario. These include the large spread in mass transfer rates, the low surface temperatures of many of the companion stars, and evidence for material at low velocities. We propose that these anomalies have a common cause: the presence of a circumbinary (CB) disk. This may be a significant component of mass transferring binaries in general. Direct detection of such CB disks may be possible but difficult, because of their low luminosity Disks surrounding Cataclysmic Binaries and spectral energy distribution peaking in the mid-IR.

An empirical Compton up-scattering model is described which reproduces both the fractional amplitude (RMS) vs. energy and the soft time lags in the 830 Hz QPO observed in 4U1608-52 on Mar. 3, 1996. A combination of two coherent variations in the coronal and soft photon temperatures (with their relative contributions determined by enforcing energy conservation) gives rise to the QPO's energy dependent characteristics. All input parameters to the model, save a characteristic plasma size and the fraction of Comptonized photons impinging on the soft photon source, are derived from the time-averaged photon energy spectrum of the same observation. Fits to the fractional RMS and phase lag data for this kilohertz QPO imply that the spatial extent of the plasma is in the range from 4 to 15 km.

The effect of a possible hydrodynamical wind on the nature of hot accretion disk solutions is studied. It is found that the advection dominated branch, in the presence of a wind, maintains the self-similar solution for the disk structure with the temperature, $\theta \propto 1/r$, optical depth $\tau \propto r^{P-1/2}$ and accretion rate $\dot m \propto r^P$. Based on global solutions, cooling due to wind energy loss and advection are found to be equally important. For a wide range of viscosity and wind parameters, the temperature is about one-tenth of the virial value and $P \approx 0.9$, independent of the mass accretion rate and radius. In the context of cooling by unsaturated Comptonization of soft photons, solutions also exist where radiative cooling, advection, and wind cooling are important. In this case, wind regulated solutions are possible. Here, the radial dependence of the critical mass accretion rate above which solutions do not exist is unchanged from those solutions without winds. The wind/advection dominated solutions are locally unstable to a new type of instability called ``wind driven insbility'' where the presence of a wind causes the disk to be unstable to long wavelength perturbations of the surface density. The growth rate of this instability is inversely proportional to the ratio of the radiative cooling to the gravitational energy dissipation rates and it can grow on a timescale much longer than the viscous timescale in the disk for sufficiently small radiative cooling efficiencies.

Hydrodynamical winds from a spherical two-temperature plasma surrounding a compact object are constructed. The mass-loss rate is computed as a function of electron temperature, optical depth and luminosity of the sphere, the values of which can be constrained by the fitting of the spectral energy distributions for known X-ray binary systems. The sensitive dependence of the mass loss rate with these parameters leads to the identification of two distinct regions in the parameter space separating wind-dominated from non wind dominated systems. A critical optical depth, tau_c, as a function of luminosity and electron temperature, is defined which differentiates these two regions. Systems with optical depths significantly smaller than tau_c are wind-dominated. The results are applied to black hole candidate X-ray binary systems in the hard spectral state (Cyg X-1, GX 339-4 and Nova Muscae), and it is found that the inferred optical depth (tau) is similar to tau_c suggesting that they are wind regulated systems. On the other hand, for X-ray binary systems containing a neutron star (e.g., Cyg X-2) tau is much larger than tau_c indicating the absence of significant hydrodynamical winds.

The influence of a circumbinary (CB) disk on the evolution of cataclysmic variable (CV) binary systems is investigated. We show that CB mass surface densities sufficient to influence the evolution rate are plausibly provided by the outflows observed in CVs, if the net effect of these winds is to deliver $10^{-4}$--$10^{-3}$ of the mass transfer rate to the CB disk. The torque exerted by the CB disk provides a positive feedback between mass transfer rate and CB disk mass which can lead to mass transfer rates of $\sim 10^{-8} -10^{-7} \mpy$. This mechanism may be responsible for causing the range of variation of mass transfer rates in CV's. In particular, it may explain rates inferred for the novalike variables and the supersoft X-ray binary systems observed near the upper edge of the period gap ($P \sim 3 - 4$ hr), as well as the spread in mass transfer rates above and below the period gap. Consquences and the possible observability of such disks are discussed.

We reexamine evolutionary channels for the formation of binary millisecond pulsars in order to understand their observed orbital period distribution. The available paths provide a natural division into systems characterized by long orbital periods (> 60 d) and short orbital periods (< 30 d). Systems with initial periods 1 - 2 d ultimately produce low mass He white dwarfs with short orbital periods (< 1 d). For longer initial periods (> few days), early massive Case B evolution produces CO white dwarfs with orbital periods < 20 d. Common envelope evolution result in short period systems (P < 1 d) from unstable low mass Case B evolution producing He white dwarfs, and from unstable Case C evolution leading to CO white dwarfs. On the other hand, the long orbital period group arises from stable low mass Case B evolution with initial periods > few days producing low mass He white dwarfs and periods > 30 d, and from stable Case C evolution producing CO white dwarfs. The lack of observed systems between 23 and 56 days probably reflects the fact that for comparable initial orbital periods (< few days) low mass Case B and early massive Case B evolution lead to very discrepant final periods. We show in particular that the lower limit (~ 23 d) cannot result from common-envelope evolution.

We consider possible evolutionary models for SS 433. We assume that common-envelope evolution is avoided if radiation pressure is able to expel most of a super-Eddington accretion flow from a region smaller than the accretor's Roche lobe. This condition is satisfied, at least initially, for largely radiative donors with masses in the range 4-12 solar masses. For donors more massive than about 5 solar masses, moderate mass ratios q = M_2/M_1 > 1 are indicated, thus tending to favor black-hole accretors. For lower mass donors, evolutionary considerations do not distinguish between a neutron star or black hole accretor. In all cases the mass transfer (and mass loss) rates are much larger than the likely mass-loss rate in the precessing jets. Almost all of the transferred mass is expelled at radii considerably larger than the jet acceleration region, producing the "stationary" H-alpha line, the infrared luminosity, and accounting for the low X-ray luminosity.

The evolution of low mass (M < 2.5 Msun) binaries through the common envelope phase has been studied for systems in which one member is on its first ascent of the red giant branch. Three-dimensional hydrodynamical simulations have been carried out for a range of red giant masses (1 - 2 Msun) with degenerate helium cores (0.28 - 0.45 Msun) and companions (0.1 - 0.45 Msun) for initial orbital periods ranging from about 15 - 1000 days. The results suggest that these low mass binary systems can survive the common envelope phase provided that the helium degenerate core is more massive than about 0.2 - 0.25 Msun and that the mass of the red giant progenitor is less than about 2 Msun. For the observed short period double degenerate systems, it is found that evolutionary scenarios involving two phases of common envelope evolution are not likely and that a scenario involving an Algol-like phase of mass transfer followed by a common envelope phase is viable, suggesting that the first-formed white dwarf is often reheated by nuclear burning on its surface. A formation mechanism for two subdwarf B stars observed in eclipsing short period binaries with low mass main sequence stars is also described.

We show that the galactic microquasar GRS 1915+105 exhibits quasi-periodic oscillations (QPOs) whose frequency varies continuously from 1-15 Hz, during spectrally hard dips when the source is in a flaring state. We report here analyses of simultaneous energy spectra and power density spectra at 4 s intervals. The energy spectrum is well fit at each time step by an optically thick accretion disk plus power law model, while the power density spectrum consists of a varying red noise component plus the variable frequency QPO. The features of both spectra are strongly correlated with one another. The 1-15 Hz QPOs appear when the power law component becomes hard and intense, and themselves have an energy spectrum consistent with the power law component (with root mean square amplitudes as high as 10%). The frequency of the oscillations, however, is most strikingly correlated with the parameters of the thermal disk component. The tightest correlation is between QPO frequency and the disk X-ray flux. This fact indicates that the properties of the QPO are not determined by solely a disk or solely a corona.

The evolution of a giant planet within the stellar envelope of a main-sequence star is investigated as a possible mechanism for enhancing the stellar metallicities of the parent stars of extrasolar planetary systems. Three-dimensional hydrodynamical simulations of a planet subject to impacting stellar matter indicate that the envelope of a Jupiter-like giant planet can be completely stripped in the outer stellar convection zone of a solar-mass star. In contrast, Jupiter-like and less massive Saturn-like giant planets are able to survive through the base of the convection zone of a 1.22 solar-mass star. Although strongly dependent on details of planetary interior models, partial or total dissolution of giant planets can result in significant enhancements in the metallicity of host stars with masses between about 1.0 and 1.3 solar masses. The implications of these results with regard to planetary orbital migration are briefly discussed.

The evolution of binary systems consisting of an asymptotic giant branch star of mass equal to 3 M_sun or 5 M_sun, and a main sequence star of mass equal to 0.4 M_sun or 0.6 M_sun with orbital periods > 200 days has been followed from the onset through the late stages of the common envelope phase. Using a nested grid technique, the three-dimensional hydrodynamical simulations of an asymptotic giant branch star with radii approximately 1 A.U. indicate that a significant fraction of the envelope gas is unbound (about 31% and 23% for binaries of 3 M_sun and 0.4 M_sun, and 5 M_sun and 0.6 M_sun respectively) by the ends of the simulations, and that the efficiency of the mass ejection process is about 40%. While the original volume of the giant is virtually evacuated in the late stages, most of the envelope gas remains marginally bound on the grid. At the ends of our simulations, when the orbital decay timescale exceeds about 5 years, the giant core and companion orbit one another with a period of about 1 day (2.4 days for a binary involving a more evolved giant), although this is an upper limit to the final orbital period. For a binary of 5 M_sun and 0.4 M_sun, the common envelope may not be completely ejected.

We present time-resolved GRS 1915+105 energy and power spectra observed by RXTE, during an episode where the X-ray intensity makes an extreme dip. If the spectra are modeled in terms of disk and power law components, both have large variations. When the inner disk is disrupted, the power law dominates, exhibiting quasi-periodic oscillations with varying frequency until the inner disk returns.

The light curves of GRS 1915+105 observed with RXTE on October 15, 1996 reveal a wide range of transient activity including regular bursts with a recurrence time of about 1 minute, irregular bursts, and dips. In contrast to bursts from other sources, a secondary (and a tertiary) weaker burst immediately following the primary burst are observed. Detailed energy spectra indicate that the source softens during the main outburst and successively hardens during the secondary and tertiary bursts. This may imply that the accretion flow has a corona-disk configuration and that the relative contribution of the hot corona decreases during the primary bursts and increases during the secondary and tertiary bursts. The primary burst profile resembles the bursts produced in the time dependent evolutions of accretion disks which are thermally and viscously unstable. The secondary burst may reflect an inward shift of the inner edge of the disk which results in a greater release of gravitational binding energy.

From data obtained from the PCA in the 2-11 keV and 11-30.5 keV energy range, GRS 1915+105 is seen during RXTE observations between 1996 May and October on two separate branches in a hardness intensity diagram. On the hard branch, GRS 1915+105 exhibits narrow quasi-periodic oscillations ranging from 0.5 to 6 Hz with ${\Delta \nu \over \nu} \sim 0.2$. The QPOs are observed over intensities ranging from about 6,000 to 20,000 counts s$^{-1}$ in the 2 - 12.5 keV energy band, indicating a strong dependence on source intensity. Strong harmonics are seen, especially, at lower frequencies. As the QPO frequency increases, the harmonic feature weakens and disappears. On the soft branch, narrow QPOs are absent and the low frequency component of the power density spectrum is approximated by a power-law, with index $\sim -1.25$ for low count rates and $\sim -1.5$ for high count rates ($\gta 18000$ cts/s). Occasionally, a broad peaked feature in the 1-6 Hz frequency range is also observed on this branch. The source was probably in the very high state similar to those of other black hole candidates. Thermal-viscous instabilities in accretion disk models do not predict the correlation of the narrow QPO frequency and luminosity unless the fraction of luminosity from the disk decreases with the total luminosity.

The morphology of adiabatic accretion flows with small non-zero specific angular momentum has been investigated in the axisymmetric and nonviscous limit. For an initial state characterized by a Bondi flow with the specific angular momentum distributed with respect to polar angle, a travelling shock wave forms which propagates more rapidly in the equatorial plane than in the plane perpendicular to it resulting in the formation of a hot torus. In cases where the incoming flow is restricted to lie near the equatorial plane, a strong wind forms directed away from this plane with the tendency for the formation of a nonsteady shock structure. As the height of the incoming flow is increased it is found that the resultant wind weakens. The parameter regime which delineates accretion flows characterized by a travelling shock and a nearly standing shock structure is presented.

We present the results of hydrodynamical simulations of nonaxisymmetric gas flow past a gravitating compact object in two dimensions. Calculations were performed with uniform flow as well as with transverse velocity and density gradients. We find that the flow is highly nonsteady, exhibiting the ``flip-flop'' behavior seen in previous studies in which accretion disks form with alternating directions of rotation. We investigate the periodicity of the flip-flop behavior, and study the effects of spatial resolution on the results. We find that the flip-flop motion creates accretion torques which, in some cases, may be large enough to explain the erratic spin behavior observed in some massive X-ray pulsars.

Double-diffusive convection refers to mixing where the effects of thermal and composition gradients compete to determine the stability of a fluid. In addition to the familiar fast convective instability, such fluids exhibit the slow, direct salt finger instability and the slow, overstable semiconvective instability. Previous approaches to this subject usually have been based on linear stability analyses. We develop here the nonlinear mixing-length theory (MLT) of double-diffusive convection, in analogy to the more familiar MLT for a fluid of homogeneous composition. We present approximate solutions for the mixing rate in the various regimes, and show that the familiar Schwarzschild and Ledoux stability criteria are good approximations to the precise criteria in stellar interiors. We have implemented the self-consistent computation of the temperature gradient and turbulent mixing rate in a stellar evolution code and solved a diffusion equation to mix composition at the appropriate rate. We have evolved $15\msun$ and $30\msun$ stars from the zero-age main sequence to the end of core He-burning. Semiconvective mixing is fast enough to alter stellar composition profiles on relevant time scales, but not so fast that instantaneous readjustment is appropriate.

A framework for the interpretation of the spectral states of black hole X-ray transients based on the diversity of accretion disk models is introduced. Depending on the mass accretion rate, it is proposed that the accretion disk is described by one or a combination of the following structures: optically thick disk, advection-dominated disk, corona-disk, and non-steady state disk. In particular, it is suggested that the very high, high, low, and off states are characterized by mass accretion rates of decreasing magnitude. The very high state corresponds to mass accretion rates near the Eddington limit in which an optically thin non steady inner region is surrounded by an optically thick structure. In the high state, the inner region is optically thin and advection-dominated or optically thick. The low hard state is interpreted in terms of a disk-corona system and the off state in terms of an optically thin disk dominated by advective energy transport into the black hole. The possible observational consequences of such a paradigm are discussed.

We argue that the bursting, transient X-ray source GRO J1744-28 is a binary consisting of a very low-mass ($M \approx 0.2 M_\sun$), highly evolved giant star that is transferring mass by Roche lobe overflow onto a high-mass ($M \approx 1.8 M_\sun$) neutron star. We explore a picture in which the bursts are due to thermonuclear flashes in matter that has accreted onto the neutron star. We attribute the unusually hard spectra of the bursts and the high burst rate compared to normal Type I X-ray burst sources -- indeed the existence of unstable nuclear burning itself -- to the unusually strong surface magnetic field, $ B_s \approx 1 \times 10^{13}$ G, of the neutron star.

The global nonlinear time dependent evolution of accretion disk-corona systems in X-ray binary sources has been investigated to provide an understanding of the low frequency ($\sim 0.04$~Hz) quasi-periodic oscillations (QPOs) observed recently in the Rapid Burster MXB 1730-335 and in some black hole candidates sources (Cyg X--1 and GRO~J0422+32). We consider $\alpha$-viscosity models in which the viscous stress is proportional to the total pressure. In contrast to previous time dependent studies it is assumed that all mass accretion and angular momentum transport take place in an optically thick disk, but that a fraction of the gravitational energy that is released is dissipated in a corona. It is found that the coronal energy dissipation can effectively reduce the amplitudes of the mass flow variations generated from the thermal and viscous instabilities (in comparison to models without a corona) to a state in which the disk exhibits mild oscillatory nonsteady behavior. These oscillations are globally coherent in the unstable regions of the disk. A model for the high and low states of black hole candidate systems is also proposed. It is suggested that the low state, which is characterized by a hard X-ray spectrum, corresponds to a disk configuration in which the inner disk is in an advection dominated hot optically thin state whereas the high state corresponds to a configuration in which the inner disk is in an optically thick state surrounded by a corona. In this model, the mass accretion rate in the system is higher in the low state than in the high state. The hard X-ray spectrum of QPOs observed in the low state can be naturally explained by such a model.

The global nonlinear time-dependent evolution of the inertial-acoustic mode instability in accretion disks surrounding black holes has been investigated. The viscous stress is assumed to be proportional to the gas pressure only, i.e.\,, $\tau = - \alpha p_g$. It is found that an oscillatory instability exists in the inner regions of disks ($r < 10 r_g$, where $r_g$ is the Schwarzschild radius) for sufficiently large $\alpha$ ($\gapprox 0.2$), and for mass accretion rates less than about 0.3 times the Eddington value. The variations of the integrated bolometric luminosity from the disk, $\Delta L/L$, are less than 3\%. A power spectrum analysis of these variations reveals a power spectrum which can be fit to a power law function of the frequency $P \propto f^{-\gamma}$, with index $\gamma \sim 1.4-2.3$ and a low frequency feature at about 4 Hz in one case. In addition, a narrow peak centered at a frequency corresponding to the maximum epicyclic frequency of the disk at $\sim 100-130$ Hz and its first harmonic is also seen. The low frequency modulations are remarkably similar to those observed in black hole candidate systems. The possible existence of a scattering corona in the inner region of the disk and/or to other processes contributing to the power at high frequencies in the inner region of the accretion disk may make the detection of the high frequency component difficult.

We consider the inner regions of accretion disks surrounding black holes and neutron stars and investigate the nonlinear time dependent evolution of thermal-viscous instabilities. The viscous stress is assumed to be proportional to the gas pressure with the viscosity parameter formulated as $\alpha =\min[\alpha_0 (h/r)^n, \alpha_{\max}]$, where $h$ is the local scale height, $r$ is the distance from the central compact object, and $n$, $\alpha_0$ and $\alpha_{\max}$ are constants. It is found that the disk is unstable for $\alpha$ sufficiently sensitive to $h$ ($n \gapprox 1.2$). The instabilities are globally coherent in the entire unstable region of the disk and, depending on the viscosity parameters, the time variability of the mass accretion rates are manifested as periodic or quasi-periodic oscillations. We show that, the low frequency ($\sim 0.04$~Hz) quasi-periodic oscillations (QPOs) discovered recently in some of the black hole candidates (Cyg~X-1 and GRO~J0422+32) and a low mass X-ray binary (Rapid Burster MXB~1730--335) may be explicable by the thermal-viscous instabilities in accretion disks. The observations of QPOs place constraints on the viscosity parameters and suggest that $(n,\alpha_0)$ $\sim (1.6,30)$ for the Rapid Burster with a $1.4\,\msun$ neutron star. In the case of black hole candidates, the dependence of $\alpha$ on $h/r$ is less steep corresponding to $n \sim 1.2-1.3$ for black holes less than $10\,\msun$.