BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ronald E. Taam

Ronald E. Taam appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-phastro-ph.GAastro-ph.COastro-ph.EPastro-ph.IM

Catalog footprint

What is connected

40works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2020arXiv

Common Envelope Shaping of Planetary Nebulae. II. Magnetic Solutions and Self-Collimated Outflows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pulsational pair-instability supernovae in very close binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wind accretion in Cygnus X-1

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Common Envelope Evolution of Massive Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Time-dependent Corotation Resonance in Barred Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-resolution CO observation of the carbon star CIT 6 revealing the spiral structure and a nascent bipolar outflow

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$).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-Gravitational Force Calculation of Second Order Accuracy for Infinitesimally Thin Gaseous Disks in Polar Coordinates

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulations of the symbiotic recurrent nova V407 Cyg. I. Accretion and shock evolutions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Search for surviving companions in type Ia supernova remnants

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spiral Density Waves in M81. I. Stellar Spiral Density Waves

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spiral Density Waves in M81. II. Hydrodynamic Simulations for Gas Response to Stellar Spiral Density Waves

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 $μ$m and 21 cm. In the outer disk, the ring-like structure observed in 8 $μ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on the Viscosity and Magnetic Field in Hot Accretion Flows around Black Holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evidence of a Binary-Induced Spiral from an Incomplete Ring Pattern of CIT 6

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamical Simulations of the Barred Spiral Galaxy NGC 1097

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Method of Determining the Characteristics of Evolved Binary Systems Revealed in the Observed Circumstellar Patterns: Application to AFGL 3068

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An AMR Study of the Common Envelope Phase of Binary Evolution

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disk corona interaction: mechanism for the disk truncation and spectrum change in low luminosity AGN

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 ($α= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Meridional tilt of the stellar velocity ellipsoid during bar buckling instability

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wide Binary Effects on Asymmetries in Asymptotic Giant Branch Circumstellar Envelopes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing Substellar Companions of AGB Stars through Spirals and Arcs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Central Region of the Nearby Seyfert 2 Galaxy NGC 4945 : A Pair of Spirals

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$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theoretical Spectra of the Am CVn Binary System SDSS J0926+3624 : Effects of Irradiation onto the Donor Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical Formation of Millisecond Pulsars in Globular Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evidence for gamma-ray emission from the low-mass X-ray binary system FIRST J102347.6+003841

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of Type Ia Supernova Ejecta on a Helium-star Binary Companion

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, $δM_{\rm ub}$, is related to the initial binary separation, $a$, by a power law of the form $δ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$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Masses and envelope binding energies of primary stars at the onset of a common envelope

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 {λ_\mathrm{env}} from the binding energy and compare its distribution to typical assumptions that are made in population-synthesis codes. We conclude that {λ_\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 {α_\mathrm{CE}}. Finally, we discuss accurate existing analytic fits for the envelope binding energy, which make these oversimplified assumptions for {λ_\mathrm{env}}, and the use of {λ_\mathrm{env}} in general, unnecessary.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Population Synthesis of Common Envelope Mergers: I. Giant Stars with Stellar or Substellar Companions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of bars and transient spirals on the vertical heating in disk galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of coronal radiation on a residual inner disk in the low/hard spectral state of black hole X-ray binary systems

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 $α=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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Interaction of Planets and Brown Dwarfs with AGB Stellar Winds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Application of the Disk Evaporation Model to AGNs

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, α, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Equilibrium Structure of Prolate Magnetized Molecular Cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

On the rarity of double black hole binaries: consequences for gravitational-wave detection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Turn-On of Mass Transfer in AM CVn Binaries: Implications for RX J0806+1527 and RX J1914+2456

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$ ($ν= 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ν \approx \ee{1.0-1.5}{-28}$ Hz s$^{-2}$ and argue that timing observations can probe $\ddotν$ at this level with a total $\approx 20$ yr baseline. We also place constraints on each system's initial parameters given current observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

On the Chandra X-ray Sources in the Galactic Center

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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Pulsar Wind Nebulae and the X-Ray Emission of Non-Accreting Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Thermal Timescale Mass Transfer and the Evolution of White Dwarf Binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Double-Diffusive Mixing-Length Theory, Semiconvection, and Massive Star Evolution

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.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

On the morphology of accretion flows with small non-zero specific angular momentum

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.

0 citations0 reviewsNo code linkedOpen access