Showing 1–50 of 106 results for author: Taam, R E

Common Envelope Shaping of Planetary Nebulae. IV. From Proto-planetary to Planetary Nebula

Authors: Guillermo García-Segura, Ronald E. Taam, Paul M. Ricker

Abstract: We present 2D hydrodynamical simulations of the transition of a proto-planetary nebula to a planetary nebula for central stars in binary systems that have undergone a common envelope event. After 1,000 yr of magnetically driven dynamics (proto-planetary nebula phase), a line-driven stellar wind is introduced into the computational domain and the expansion of the nebula is simulated for another 10,… ▽ More We present 2D hydrodynamical simulations of the transition of a proto-planetary nebula to a planetary nebula for central stars in binary systems that have undergone a common envelope event. After 1,000 yr of magnetically driven dynamics (proto-planetary nebula phase), a line-driven stellar wind is introduced into the computational domain and the expansion of the nebula is simulated for another 10,000 yr, including the effects of stellar photoionization. In this study we consider central stars with main sequence (final) masses of 1 (0.569) and 2.5 (0.677) \Mo, together with a 0.6 \Mo ma in sequence companion. Extremely bipolar, narrow-waisted proto-planetary nebulae result in bipolar planetary nebulae, while the rest of the shapes mainly evolve into elliptical planetary nebulae. The initial magnetic field's effects on the collimated structures, such as jets, tend to disappear in most of the cases, leaving behind the remnants of those features in only a few cases. Equatorial zones fragmented mainly by photoionization ( 1 \Mo progenitors), result in ``necklace'' structures made of cometary clumps aligned with the radiation field. On the other hand, fragmentation by photoionization and shocked wind ( 2.5 \Mo progenitors) give rise to the formation of multiple clumps in the latitudinal direction, which remain within the lobes, close to the center, which are immersed and surrounded by hot shocked gas, not necessarily aligned with the radiation field. These results reveal that the fragmentation process has a dependence on the stellar mass progenitor. This fragmentation is made possible by the distribution of gas in the previous post-common envelope proto-planetary nebula as sculpted by the action of the jets. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 12 pages, 7 figures, accepted by MNRAS

Common Envelope Shaping of Planetary Nebulae. III. The Launching of Jets in Proto-Planetary Nebulae

Authors: Guillermo Garcia-Segura, Ronald E. Taam, Paul M. Ricker

Abstract: We compute successfully the launching of two magnetic winds from two circumbinary disks formed after a common envelope event. The launching is produced by the increase of magnetic pressure due to the collapse of the disks. The collapse is due to internal torques produced by a weak poloidal magnetic field. The first wind can be described as a wide jet, with an average mass-loss rate of… ▽ More We compute successfully the launching of two magnetic winds from two circumbinary disks formed after a common envelope event. The launching is produced by the increase of magnetic pressure due to the collapse of the disks. The collapse is due to internal torques produced by a weak poloidal magnetic field. The first wind can be described as a wide jet, with an average mass-loss rate of $\sim 1.3 \times 10^{-7}$ \Moy\ and a maximum radial velocity of $\sim 230$ \kms. The outflow has a half-opening angle of $\sim 20^{\circ}$. Narrow jets are also formed intermittently with velocities up to 3,000 \kms, with mass-loss rates of $\sim 6 \times 10^{-12} $ \Moy\ during short periods of time. The second wind can be described as a wide X-wind, with an average mass-loss rate of $\sim 1.68 \times 10^{-7}$ \Moy\ and a velocity of $\sim 30$ \kms. A narrow jet is also formed with a velocity of 250 \kms, and a mass-loss rates of $\sim 10^{-12} $ \Moy. The computed jets are used to provide inflow boundary conditions for simulations of proto-planetary nebulae. The wide jet evolves into a molecular collimated outflow within a few astronomical units, producing proto-planetary nebulae with bipolar, elongated shapes, whose kinetic energies reach $\sim 4 \times 10^{45}$ erg at 1,000 years. Similarities with observed features in W43A, OH231.8+4.2, and Hen 3-1475 are discussed. The computed wide X-wind produces proto-planetary nebulae with slower expansion velocities, with bipolar and elliptical shapes, and possible starfish type and quadrupolar morphology. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: 35 pages, 17 figures, accepted by the Astrophysical Journal

The role of mass transfer and common envelope evolution in the formation of merging binary black holes

Authors: Pablo Marchant, Kaliroë M. W. Pappas, Monica Gallegos-Garcia, Christopher P. L. Berry, Ronald E. Taam, Vicky Kalogera, Philipp Podsiadlowski

Abstract: As the number of observed merging binary black holes (BHs) grows, accurate models are required to disentangle multiple formation channels. In models with isolated binaries, important uncertainties remain regarding the stability of mass transfer (MT) and common-envelope (CE) evolution. To study some of these uncertainties, we have computed simulations using MESA of a $30M_\odot$, low metallicity (… ▽ More As the number of observed merging binary black holes (BHs) grows, accurate models are required to disentangle multiple formation channels. In models with isolated binaries, important uncertainties remain regarding the stability of mass transfer (MT) and common-envelope (CE) evolution. To study some of these uncertainties, we have computed simulations using MESA of a $30M_\odot$, low metallicity ($Z_\odot/10$) star with a BH companion. We developed a prescription to compute MT rates including possible outflows from outer Lagrangian points, and a method to self-consistently determine the core-envelope boundary in the case of CE evolution. We find that binaries survive a CE only if unstable MT happens after the formation of a deep convective envelope, resulting in a narrow range (0.2 dex) in period for envelope ejection. All cases where interaction is initiated with a radiative envelope have large binding energies ($\sim 10^{50}$ erg), and merge during CE even under the assumption that all the internal and recombination energy of the envelope, as well as the energy from an inspiral, is used for ejection. This is independent of core helium ignition for the donor, a condition under which various rapid-population synthesis calculations assume a successful ejection is possible. Moreover, we find that the critical mass ratio for instability is such that for periods between $\sim 1-1000$ days merging binary BHs can be formed via stable MT. A large fraction of these systems overflow their L$_2$ equipotential, in which case we find stable MT produces merging binary BHs even under extreme assumptions of mass and angular momentum outflows. Our conclusions are limited to the study of one donor star, but suggest that population synthesis calculations overestimate the formation rate of merging binary BHs produced by CE evolution, and that stable MT could dominate the rate from isolated binaries. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 22 pages (including 4 pages of appendices). Abstract abridged for submission to the arXiv. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 650, A107 (2021)

Wind accretion in Cygnus X-1

Authors: Emmi Meyer-Hofmeister, B. F. Liu, E. Qiao, Ronald E. Taam

Abstract: 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-… ▽ More 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. △ Less

Submitted 27 April, 2020; originally announced April 2020.

Comments: 8 pages, 7 figures, to appear in A&A 2020

Journal ref: A&A 637, A66 (2020)

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

Authors: Guillermo García-Segura, Ronald E. Taam, Paul M. Ricker

Abstract: 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,… ▽ More 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. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: 21 pages, 7 figures, accepted by ApJ

The Complete Evolution of a Neutron-Star Binary through a Common Envelope Phase Using 1D Hydrodynamic Simulations

Authors: Tassos Fragos, Jeff J. Andrews, Enrico Ramirez-Ruiz, Georges Meynet, Vicky Kalogera, Ronald E. Taam, Andreas Zezas

Abstract: 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 mode… ▽ More 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 $α_{\rm CE}$-equivalent efficiency of $\approx$5. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages, 5 figures, submitted to The Astrophysical Journal Letters. Comments welcome

Templates of binary-induced spiral-shell patterns around mass-losing post-main sequence stars

Authors: Hyosun Kim, Sheng-Yuan Liu, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 14 June, 2019; originally announced June 2019.

Comments: 58 pages, 49 figures (17 for main text and 32 in appendix), accepted for publication in the ApJS

Gas Flows Within Cavities of Circumbinary Discs in Eccentric Binary Protostellar Systems

Authors: Philipp Mösta, Ronald E. Taam, Paul C. Duffell

Abstract: 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, reflec… ▽ More 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. △ Less

Submitted 22 April, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 8 pages, 5 figures

Journal ref: ApJL, 875, L21 (2019)

Common Envelope Evolution of Massive Stars

Authors: Paul M. Ricker, Frank X. Timmes, Ronald E. Taam, Ronald F. Webbink

Abstract: 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… ▽ More 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. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 6 pages, 3 figures, oral contribution: IAU Symposium 346 "High Mass X-ray Binaries: illuminating the passage from massive binaries to merging compact objects", Vienna, Austria, 27-31 August 2018

Journal ref: Proc. IAU 14 (2018) 449-454

Pulsational pair-instability supernovae in very close binaries

Authors: Pablo Marchant, Mathieu Renzo, Robert Farmer, Kaliroe M. W. Pappas, Ronald E. Taam, Selma de Mink, Vassiliki Kalogera

Abstract: 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 w… ▽ More 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. △ Less

Submitted 12 July, 2020; v1 submitted 31 October, 2018; originally announced October 2018.

Comments: 22 pages, 17 figures, including 4 pages of appendices

Journal ref: The Astrophysical Journal, Volume 882, Issue 1, article id. 36, 17 pp. (2019)

Time-dependent Pattern Speeds in Barred Galaxies

Authors: Yu-Ting Wu, Daniel Pfenniger, Ronald E. Taam

Abstract: 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 pa… ▽ More 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. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: 25 pages, 17 figures, 1 table. ApJ accepted

The growth of stellar mass black hole binaries trapped in an accretion disk of active galactic nuclei

Authors: Shu-Xu Yi, K. S. Cheng, Ronald E. Taam

Abstract: 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 massiv… ▽ More 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. △ Less

Submitted 17 May, 2018; originally announced May 2018.

Comments: Accepted for publication in ApJL

A Model for Spectral States and Their Transition in Cyg X-1

Authors: Ronald E. Taam, Erlin Qiao, B. F. Liu, E. Meyer-Hofmeister

Abstract: 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 har… ▽ More 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. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 13 pages, 2 figures. Accepted for publication in ApJ

Journal ref: ApJ 2018

Common Envelope Shaping of Planetary Nebulae

Authors: Guillermo Garcia-Segura, Paul M. Ricker, Ronald E. Taam

Abstract: 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 t… ▽ More 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. △ Less

Submitted 2 May, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: 25 pages, 9 figures, Accepted by ApJ

Centrally Concentrated X-ray Radiation from an Extended Accreting Corona in Active Galactic Nuclei

Authors: B. F. Liu, Ronald E. Taam, Erlin Qiao, Weimin Yuan

Abstract: 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… ▽ More 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}{σ_{\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. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: 11 pages, 4 figures

Journal ref: ApJ 847:96, 2017 Oct.1

The large-scale nebular pattern of a superwind binary in an eccentric orbit

Authors: Hyosun Kim, Alfonso Trejo, Sheng-Yuan Liu, Raghvendra Sahai, Ronald E. Taam, Mark R. Morris, Naomi Hirano, I-Ta Hsieh

Abstract: 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 simultane… ▽ More 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. △ Less

Submitted 3 April, 2017; originally announced April 2017.

Comments: Published in Nature Astronomy (cover story of the March issue, 2017), Authors' edited version (17 pages, 3 main figures, 4 supplementary figures, and 1 supplementary video)

Journal ref: Nat. Astron. 1, 0060 (2017)

Time-dependent Corotation Resonance in Barred Galaxies

Authors: Yu-Ting Wu, Daniel Pfenniger, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 3 August, 2016; originally announced August 2016.

Comments: 10 pages, 4 figures, 1 table. ApJ accepted

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

Authors: Hsiang-Hsu Wang, Ronald E. Taam, David C. C. Yen

Abstract: 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 genera… ▽ More 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. △ Less

Submitted 9 March, 2016; originally announced March 2016.

Comments: 4 figures, 5 tables, accepted for publication in ApJS

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

Authors: Hyosun Kim, Sheng-Yuan Liu, Naomi Hirano, Ronny Zhao-Geisler, Alfonso Trejo, Hsi-Wei Yen, Ronald E. Taam, Francisca Kemper, Jongsoo Kim, Do-Young Byun, Tie Liu

Abstract: 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… ▽ More 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$). △ Less

Submitted 13 October, 2015; originally announced October 2015.

Comments: 35 pages, 14 figures, 3 tables, accepted for publication in the Astrophysical Journal

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

Authors: Hsiang-Hsu Wang, David C. C. Yen, Ronald E. Taam

Abstract: 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- for… ▽ More 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. △ Less

Submitted 4 September, 2015; originally announced September 2015.

Comments: 35 pages, 7 figures, accepted for publication in ApJS

A Hybrid Two Component Accretion Flow Surrounding Supermassive Black Holes in AGN

Authors: B. F. Liu, R. E. Taam, E. Qiao, W. Yuan

Abstract: 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… ▽ More 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 $α_{\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. △ Less

Submitted 4 May, 2015; originally announced May 2015.

Comments: 21 pages,4 figures, accepted for publication in ApJ

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

Authors: Kuo-Chuan Pan, Paul M. Ricker, Ronald E. Taam

Abstract: 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 studie… ▽ More 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. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Comments: ApJ, In Press. Simulation Animation: https://youtu.be/g5Nu7vDfCOg

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

Authors: Hsiang-Hsu Wang, Wing-Kit Lee, Ronald E. Taam, Chien-Chang Feng, Lien-Hsuan Lin

Abstract: 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 wavele… ▽ More 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. △ Less

Submitted 18 December, 2014; originally announced December 2014.

Comments: 20 pages, 10 figures

Search for surviving companions in type Ia supernova remnants

Authors: Kuo-Chuan Pan, Paul M. Ricker, Ronald E. Taam

Abstract: 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 he… ▽ More 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. △ Less

Submitted 25 July, 2014; originally announced July 2014.

Comments: 25 pages, 5 figures, and 2 tables. Accepted for publication in ApJ

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

Authors: Chien-Chang Feng, Lien-Hsuan Lin, Hsiang-Hsu Wang, Ronald E. Taam

Abstract: 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 contrib… ▽ More 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. △ Less

Submitted 28 February, 2014; originally announced February 2014.

Comments: 25 pages, 6 figures, accepted for publication in ApJ

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

Authors: Hyosun Kim, I-Ta Hsieh, Sheng-Yuan Liu, Ronald E. Taam

Abstract: 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 spira… ▽ More 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. △ Less

Submitted 19 August, 2013; originally announced August 2013.

Comments: 21 pages, 6 figures, 1 table, accepted for publication in the Astrophysical Journal

Accretion, Ablation and Propeller Evolution in Close Millisecond Pulsar Binary Systems

Authors: P. D. Kiel, R. E Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 5 August, 2013; originally announced August 2013.

Comments: Accepted to Astrophysics and Space Science, 21 pages, 12 figures

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

Authors: B. F. Liu, Ronald E. Taam

Abstract: 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 do… ▽ More 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. △ Less

Submitted 25 June, 2013; originally announced June 2013.

Comments: Accepted for publication in ApJS

Hydrodynamical Simulations of the Barred Spiral Galaxy NGC 1097

Authors: Lien-Hsuan Lin, Hsiang-Hsu Wang, Pei-Ying Hsieh, Ronald E. Taam, Chao-Chin Yang, David C. C. Yen

Abstract: 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… ▽ More 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. △ Less

Submitted 8 May, 2013; originally announced May 2013.

Comments: 32 pages, 14 figures, 1 table, accepted for publication in the ApJ

Meridional tilt of the stellar velocity ellipsoid during bar buckling instability

Authors: Kanak Saha, Daniel Pfenniger, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 24 December, 2012; v1 submitted 6 November, 2012; originally announced November 2012.

Comments: 10 pages, 16 Figures, accepted by the ApJ

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

Authors: Hyosun Kim, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 27 September, 2012; originally announced September 2012.

Comments: 14 pages, 5 figures, accepted for publication in the Astrophysical Journal Letters

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

Authors: Ronald E. Taam, B. F. Liu, W. Yuan, E. Qiao

Abstract: 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… ▽ More 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. △ Less

Submitted 22 September, 2012; originally announced September 2012.

Comments: 21 pages, 2 figures. Accepted for publication in ApJ

Wide Binary Effects on Asymmetries in Asymptotic Giant Branch Circumstellar Envelopes

Authors: Hyosun Kim, Ronald E. Taam

Abstract: 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, wh… ▽ More 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. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Comments: 33 pages, 11 figures, 4 tables, accepted for publication in the Astrophysical Journal

X-ray and Gamma-ray Emissions from Different Evolutionary Stage of Rotation Powered Millisecond Pulsars

Authors: J. Takata, K. S. Cheng, R. E. Taam

Abstract: The $Fermi$-LAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the Galactic $γ$-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 proper… ▽ More The $Fermi$-LAT has revealed that rotation powered millisecond pulsars (MSPs) are a major contributor to the Galactic $γ$-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 $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 18 November, 2011; v1 submitted 15 November, 2011; originally announced November 2011.

Comments: 46page, 10 figures, 3 tables, to Appear in ApJ

Probing Substellar Companions of AGB Stars through Spirals and Arcs

Authors: Hyosun Kim, Ronald E. Taam

Abstract: 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 circum… ▽ More 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. △ Less

Submitted 6 October, 2011; originally announced October 2011.

Comments: 23 pages, 7 figures, 1 table, accepted for publication in Astrophysical Journal

An AMR Study of the Common Envelope Phase of Binary Evolution

Authors: Paul M. Ricker, Ronald E. Taam

Abstract: 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 eje… ▽ More 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. △ Less

Submitted 15 January, 2012; v1 submitted 19 July, 2011; originally announced July 2011.

Comments: 18 pages, 11 figures, submitted to ApJ; accepted version

Report number: NSF-KITP-11-085

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

Authors: Sujan Sengupta, Ronald E. Taam

Abstract: 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-c… ▽ More 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. △ Less

Submitted 7 July, 2011; originally announced July 2011.

Comments: 14 pages Latex (AASTEX) including 1 black & white figure and 3 colour figures (all eps), Accepted for publication in The Astrophysical Journal

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

Authors: Lien-Hsuan Lin, Ronald E. Taam, David C. C. Yen, S. Muller, J. Lim

Abstract: 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 resul… ▽ More 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. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: 26 pages, 8 figures, 1 table, accepted for publication in ApJ

The Interaction of Planets and Brown Dwarfs with AGB Stellar Winds

Authors: Hyosun Kim, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2010; originally announced December 2010.

Comments: Part of PlanetsbeyondMS/2010 proceedings http://arxiv.org/html/1011.6606v1, Proc. of the workshop on "Planetary Systems beyond the Main Sequence" (Bamberg, 11-14 August 2010), AIPC in press (eds. S. Schuh, H. Drechsel and U. Heber), 7 pages, 2 figures

Report number: PlanetsbeyondMS/2010/05

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

Authors: B. F. Liu, C. Done, Ronald E. Taam

Abstract: 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 form… ▽ More 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. △ Less

Submitted 11 November, 2010; originally announced November 2010.

Comments: To appear in The Astrophysical Journal

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

Authors: P. H. Thomas Tam, Chung-Yue Hui, Regina H. H. Huang, Albert K. H. Kong, Jumpei Takata, Lupin C. C. Lin, Y. J. Yang, K. S. Cheng, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2010; v1 submitted 20 October, 2010; originally announced October 2010.

Comments: 10 pages, 3 figures, 2 tables, ApJ Letters, in press; corrections made for the gamma flux values

Journal ref: Astrophysical Journal Letters 724 (2010) L207-L211

The Impact of Type Ia Supernova Ejecta on Binary Companions

Authors: P. M. Ricker, K. -C. Pan, R. E. Taam

Abstract: 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 init… ▽ More 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. △ Less

Submitted 11 October, 2010; originally announced October 2010.

Comments: 6 pages, 2 figures, to appear in the proceedings of the conference "Binary Star Evolution: Mass Loss, Accretion, and Mergers" at Mykonos, Greece, June 22-25, 2010

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

Authors: Marc van der Sluys, Michael Politano, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 4 October, 2010; originally announced October 2010.

Comments: 6 pages, 3 figures, 1 table; to be published in the proceedings of the conference "Binary Star Evolution", in Mykonos, Greece, held in June 22-25, 2010

On the Transition from Accretion Powered to Rotation Powered Millisecond Pulsars

Authors: J. Takata, K. S. Cheng, R. E. Taam

Abstract: The heating associated with the deposition of $γ$-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 $γ$-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… ▽ More The heating associated with the deposition of $γ$-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 $γ$-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 $γ$-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. △ Less

Submitted 24 September, 2010; originally announced September 2010.

Comments: 13 pages, 3 figures, accepted in ApJL

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

Authors: Kanak Saha, Yao Huan Tseng, Ronald E. Taam

Abstract: 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 pot… ▽ More 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. △ Less

Submitted 4 August, 2010; originally announced August 2010.

Comments: 14 pages, 20 figures, Accepted for publication in ApJ

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

Authors: Michael Politano, Marc van der Sluys, Ronald E. Taam, Bart Willems

Abstract: 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 evolutio… ▽ More 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. △ Less

Submitted 26 July, 2010; originally announced July 2010.

Comments: 16 pages, 12 figures, accepted for publication in the Astrophysical Journal

Adiabatic Mass Loss and the Outcome of the Common Envelope Phase of Binary Evolution

Authors: Christopher J. Deloye, Ronald E. Taam

Abstract: 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 imp… ▽ More 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. △ Less

Submitted 21 July, 2010; v1 submitted 6 July, 2010; originally announced July 2010.

Comments: 4 pages, 3 figures; corrected typo in equation (1); updated reference data

Journal ref: Astrophysical Journal Letters 719 (2010) L28-L31

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

Authors: Kuo-Chuan Pan, Paul M. Ricker, Ronald E. Taam

Abstract: 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… ▽ More 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$. △ Less

Submitted 15 April, 2010; v1 submitted 2 April, 2010; originally announced April 2010.

Comments: Accepted in ApJ, 9 pages, 9 figures

Dynamical Formation of Millisecond Pulsars in Globular Clusters

Authors: C. Y. Hui, K. S. Cheng, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 22 March, 2010; originally announced March 2010.

Comments: 18 pages, 3 figures, 4 tables, accepted for publication in ApJ

The Equilibrium Structure of Prolate Magnetized Molecular Cores

Authors: Michael J. Cai, Ronald E. Taam

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2009; originally announced December 2009.

Comments: 5 pages, 3 figures. Accepted for publication by ApJL