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

Let Me Speak Freely? A Study on the Impact of Format Restrictions on Performance of Large Language Models

Authors: Zhi Rui Tam, Cheng-Kuang Wu, Yi-Lin Tsai, Chieh-Yen Lin, Hung-yi Lee, Yun-Nung Chen

Abstract: Structured generation, the process of producing content in standardized formats like JSON and XML, is widely utilized in real-world applications to extract key output information from large language models (LLMs). This study investigates whether such constraints on generation space impact LLMs abilities, including reasoning and domain knowledge comprehension. Specifically, we evaluate LLMs perform… ▽ More Structured generation, the process of producing content in standardized formats like JSON and XML, is widely utilized in real-world applications to extract key output information from large language models (LLMs). This study investigates whether such constraints on generation space impact LLMs abilities, including reasoning and domain knowledge comprehension. Specifically, we evaluate LLMs performance when restricted to adhere to structured formats versus generating free-form responses across various common tasks. Surprisingly, we observe a significant decline in LLMs reasoning abilities under format restrictions. Furthermore, we find that stricter format constraints generally lead to greater performance degradation in reasoning tasks. △ Less

Submitted 14 October, 2024; v1 submitted 5 August, 2024; originally announced August 2024.

Comments: 18 pages

FPT+: A Parameter and Memory Efficient Transfer Learning Method for High-resolution Medical Image Classification

Authors: Yijin Huang, Pujin Cheng, Roger Tam, Xiaoying Tang

Abstract: The success of large-scale pre-trained models has established fine-tuning as a standard method for achieving significant improvements in downstream tasks. However, fine-tuning the entire parameter set of a pre-trained model is costly. Parameter-efficient transfer learning (PETL) has recently emerged as a cost-effective alternative for adapting pre-trained models to downstream tasks. Despite its ad… ▽ More The success of large-scale pre-trained models has established fine-tuning as a standard method for achieving significant improvements in downstream tasks. However, fine-tuning the entire parameter set of a pre-trained model is costly. Parameter-efficient transfer learning (PETL) has recently emerged as a cost-effective alternative for adapting pre-trained models to downstream tasks. Despite its advantages, the increasing model size and input resolution present challenges for PETL, as the training memory consumption is not reduced as effectively as the parameter usage. In this paper, we introduce Fine-grained Prompt Tuning plus (FPT+), a PETL method designed for high-resolution medical image classification, which significantly reduces memory consumption compared to other PETL methods. FPT+ performs transfer learning by training a lightweight side network and accessing pre-trained knowledge from a large pre-trained model (LPM) through fine-grained prompts and fusion modules. Specifically, we freeze the LPM and construct a learnable lightweight side network. The frozen LPM processes high-resolution images to extract fine-grained features, while the side network employs the corresponding down-sampled low-resolution images to minimize the memory usage. To enable the side network to leverage pre-trained knowledge, we propose fine-grained prompts and fusion modules, which collaborate to summarize information through the LPM's intermediate activations. We evaluate FPT+ on eight medical image datasets of varying sizes, modalities, and complexities. Experimental results demonstrate that FPT+ outperforms other PETL methods, using only 1.03% of the learnable parameters and 3.18% of the memory required for fine-tuning an entire ViT-B model. Our code is available at https://github.com/YijinHuang/FPT. △ Less

Submitted 5 August, 2024; originally announced August 2024.

I Need Help! Evaluating LLM's Ability to Ask for Users' Support: A Case Study on Text-to-SQL Generation

Authors: Cheng-Kuang Wu, Zhi Rui Tam, Chao-Chung Wu, Chieh-Yen Lin, Hung-yi Lee, Yun-Nung Chen

Abstract: This study explores the proactive ability of LLMs to seek user support. We propose metrics to evaluate the trade-off between performance improvements and user burden, and investigate whether LLMs can determine when to request help under varying information availability. Our experiments show that without external feedback, many LLMs struggle to recognize their need for user support. The findings hi… ▽ More This study explores the proactive ability of LLMs to seek user support. We propose metrics to evaluate the trade-off between performance improvements and user burden, and investigate whether LLMs can determine when to request help under varying information availability. Our experiments show that without external feedback, many LLMs struggle to recognize their need for user support. The findings highlight the importance of external signals and provide insights for future research on improving support-seeking strategies. Source code: https://github.com/appier-research/i-need-help △ Less

Submitted 29 September, 2024; v1 submitted 20 July, 2024; originally announced July 2024.

Comments: Accepted by EMNLP 2024 Main Conference

StreamBench: Towards Benchmarking Continuous Improvement of Language Agents

Authors: Cheng-Kuang Wu, Zhi Rui Tam, Chieh-Yen Lin, Yun-Nung Chen, Hung-yi Lee

Abstract: Recent works have shown that large language model (LLM) agents are able to improve themselves from experience, which is an important ability for continuous enhancement post-deployment. However, existing benchmarks primarily evaluate their innate capabilities and do not assess their ability to improve over time. To address this gap, we introduce StreamBench, a pioneering benchmark designed to evalu… ▽ More Recent works have shown that large language model (LLM) agents are able to improve themselves from experience, which is an important ability for continuous enhancement post-deployment. However, existing benchmarks primarily evaluate their innate capabilities and do not assess their ability to improve over time. To address this gap, we introduce StreamBench, a pioneering benchmark designed to evaluate the continuous improvement of LLM agents over an input-feedback sequence. StreamBench simulates an online learning environment where LLMs receive a continuous flow of feedback stream and iteratively enhance their performance. In addition, we propose several simple yet effective baselines for improving LLMs on StreamBench, and provide a comprehensive analysis to identify critical components that contribute to successful streaming strategies. Our work serves as a stepping stone towards developing effective online learning strategies for LLMs, paving the way for more adaptive AI systems in streaming scenarios. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Fine-grained Prompt Tuning: A Parameter and Memory Efficient Transfer Learning Method for High-resolution Medical Image Classification

Authors: Yijin Huang, Pujin Cheng, Roger Tam, Xiaoying Tang

Abstract: Parameter-efficient transfer learning (PETL) is proposed as a cost-effective way to transfer pre-trained models to downstream tasks, avoiding the high cost of updating entire large-scale pre-trained models (LPMs). In this work, we present Fine-grained Prompt Tuning (FPT), a novel PETL method for medical image classification. FPT significantly reduces memory consumption compared to other PETL metho… ▽ More Parameter-efficient transfer learning (PETL) is proposed as a cost-effective way to transfer pre-trained models to downstream tasks, avoiding the high cost of updating entire large-scale pre-trained models (LPMs). In this work, we present Fine-grained Prompt Tuning (FPT), a novel PETL method for medical image classification. FPT significantly reduces memory consumption compared to other PETL methods, especially in high-resolution input contexts. To achieve this, we first freeze the weights of the LPM and construct a learnable lightweight side network. The frozen LPM takes high-resolution images as input to extract fine-grained features, while the side network is fed low-resolution images to reduce memory usage. To allow the side network to access pre-trained knowledge, we introduce fine-grained prompts that summarize information from the LPM through a fusion module. Important tokens selection and preloading techniques are employed to further reduce training cost and memory requirements. We evaluate FPT on four medical datasets with varying sizes, modalities, and complexities. Experimental results demonstrate that FPT achieves comparable performance to fine-tuning the entire LPM while using only 1.8% of the learnable parameters and 13% of the memory costs of an encoder ViT-B model with a 512 x 512 input resolution. △ Less

Submitted 2 July, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: MICCAI 2024

MProtoNet: A Case-Based Interpretable Model for Brain Tumor Classification with 3D Multi-parametric Magnetic Resonance Imaging

Authors: Yuanyuan Wei, Roger Tam, Xiaoying Tang

Abstract: Recent applications of deep convolutional neural networks in medical imaging raise concerns about their interpretability. While most explainable deep learning applications use post hoc methods (such as GradCAM) to generate feature attribution maps, there is a new type of case-based reasoning models, namely ProtoPNet and its variants, which identify prototypes during training and compare input imag… ▽ More Recent applications of deep convolutional neural networks in medical imaging raise concerns about their interpretability. While most explainable deep learning applications use post hoc methods (such as GradCAM) to generate feature attribution maps, there is a new type of case-based reasoning models, namely ProtoPNet and its variants, which identify prototypes during training and compare input image patches with those prototypes. We propose the first medical prototype network (MProtoNet) to extend ProtoPNet to brain tumor classification with 3D multi-parametric magnetic resonance imaging (mpMRI) data. To address different requirements between 2D natural images and 3D mpMRIs especially in terms of localizing attention regions, a new attention module with soft masking and online-CAM loss is introduced. Soft masking helps sharpen attention maps, while online-CAM loss directly utilizes image-level labels when training the attention module. MProtoNet achieves statistically significant improvements in interpretability metrics of both correctness and localization coherence (with a best activation precision of $0.713\pm0.058$) without human-annotated labels during training, when compared with GradCAM and several ProtoPNet variants. The source code is available at https://github.com/aywi/mprotonet. △ Less

Submitted 14 April, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 15 pages, 5 figures, 1 table; accepted for oral presentation at MIDL 2023 (https://openreview.net/forum?id=6Wbj3QCo4U4 ); camera-ready version

SSiT: Saliency-guided Self-supervised Image Transformer for Diabetic Retinopathy Grading

Authors: Yijin Huang, Junyan Lyu, Pujin Cheng, Roger Tam, Xiaoying Tang

Abstract: Self-supervised Learning (SSL) has been widely applied to learn image representations through exploiting unlabeled images. However, it has not been fully explored in the medical image analysis field. In this work, Saliency-guided Self-Supervised image Transformer (SSiT) is proposed for Diabetic Retinopathy (DR) grading from fundus images. We novelly introduce saliency maps into SSL, with a goal of… ▽ More Self-supervised Learning (SSL) has been widely applied to learn image representations through exploiting unlabeled images. However, it has not been fully explored in the medical image analysis field. In this work, Saliency-guided Self-Supervised image Transformer (SSiT) is proposed for Diabetic Retinopathy (DR) grading from fundus images. We novelly introduce saliency maps into SSL, with a goal of guiding self-supervised pre-training with domain-specific prior knowledge. Specifically, two saliency-guided learning tasks are employed in SSiT: (1) Saliency-guided contrastive learning is conducted based on the momentum contrast, wherein fundus images' saliency maps are utilized to remove trivial patches from the input sequences of the momentum-updated key encoder. Thus, the key encoder is constrained to provide target representations focusing on salient regions, guiding the query encoder to capture salient features. (2) The query encoder is trained to predict the saliency segmentation, encouraging the preservation of fine-grained information in the learned representations. To assess our proposed method, four publicly-accessible fundus image datasets are adopted. One dataset is employed for pre-training, while the three others are used to evaluate the pre-trained models' performance on downstream DR grading. The proposed SSiT significantly outperforms other representative state-of-the-art SSL methods on all downstream datasets and under various evaluation settings. For example, SSiT achieves a Kappa score of 81.88% on the DDR dataset under fine-tuning evaluation, outperforming all other ViT-based SSL methods by at least 9.48%. △ Less

Submitted 12 March, 2024; v1 submitted 19 October, 2022; originally announced October 2022.

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

Identifying the key components in ResNet-50 for diabetic retinopathy grading from fundus images: a systematic investigation

Authors: Yijin Huang, Li Lin, Pujin Cheng, Junyan Lyu, Roger Tam, Xiaoying Tang

Abstract: Although deep learning based diabetic retinopathy (DR) classification methods typically benefit from well-designed architectures of convolutional neural networks, the training setting also has a non-negligible impact on the prediction performance. The training setting includes various interdependent components, such as objective function, data sampling strategy and data augmentation approach. To i… ▽ More Although deep learning based diabetic retinopathy (DR) classification methods typically benefit from well-designed architectures of convolutional neural networks, the training setting also has a non-negligible impact on the prediction performance. The training setting includes various interdependent components, such as objective function, data sampling strategy and data augmentation approach. To identify the key components in a standard deep learning framework (ResNet-50) for DR grading, we systematically analyze the impact of several major components. Extensive experiments are conducted on a publicly-available dataset EyePACS. We demonstrate that (1) the DR grading framework is sensitive to input resolution, objective function, and composition of data augmentation, (2) using mean square error as the loss function can effectively improve the performance with respect to a task-specific evaluation metric, namely the quadratically-weighted Kappa, (3) utilizing eye pairs boosts the performance of DR grading and (4) using data resampling to address the problem of imbalanced data distribution in EyePACS hurts the performance. Based on these observations and an optimal combination of the investigated components, our framework, without any specialized network design, achieves the state-of-the-art result (0.8631 for Kappa) on the EyePACS test set (a total of 42670 fundus images) with only image-level labels. We also examine the proposed training practices on other fundus datasets and other network architectures to evaluate their generalizability. Our codes and pre-trained model are available at https://github.com/YijinHuang/pytorch-classification. △ Less

Submitted 17 October, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Best Practices for Data Publication in the Astronomical Literature

Authors: Tracy X. Chen, Marion Schmitz, Joseph M. Mazzarella, Xiuqin Wu, Julian C. van Eyken, Alberto Accomazzi, Rachel L. Akeson, Mark Allen, Rachael Beaton, G. Bruce Berriman, Andrew W. Boyle, Marianne Brouty, Ben Chan, Jessie L. Christiansen, David R. Ciardi, David Cook, Raffaele D'Abrusco, Rick Ebert, Cren Frayer, Benjamin J. Fulton, Christopher Gelino, George Helou, Calen B. Henderson, Justin Howell, Joyce Kim , et al. (20 additional authors not shown)

Abstract: We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to… ▽ More We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multi-dimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union (IAU) for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors. △ Less

Submitted 16 April, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 19 pages, 1 figure, 3 tables, accepted for publication in ApJS

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)

Introduction to a novel T2 relaxation analysis method SAME-ECOS: Spectrum Analysis for Multiple Exponentials via Experimental Condition Oriented Simulation

Authors: Hanwen Liu, Qing-San Xiang, Roger Tam, Piotr Kozlowski, David K. B. Li, Alex L. Mackay, John K. Kramer, Cornelia Laule

Abstract: We propose a novel T2 relaxation data analysis method which we have named spectrum analysis for multiple exponentials via experimental condition oriented simulation (SAME-ECOS). SAME-ECOS, which was developed based on a combination of information theory and machine learning neural network algorithms, is tailored for different MR experimental conditions, decomposing multi-exponential decay data int… ▽ More We propose a novel T2 relaxation data analysis method which we have named spectrum analysis for multiple exponentials via experimental condition oriented simulation (SAME-ECOS). SAME-ECOS, which was developed based on a combination of information theory and machine learning neural network algorithms, is tailored for different MR experimental conditions, decomposing multi-exponential decay data into T2 spectra, which had been considered an ill-posed problem using conventional fitting algorithms, including the commonly used non-negative least squares (NNLS) method. Our results demonstrated that, compared with NNLS, the simulation-derived SAME-ECOS model yields much more reliable T2 spectra in a dramatically shorter time, increasing the feasibility of multi-component T2 decay analysis in clinical settings. △ Less

Submitted 14 September, 2020; originally announced September 2020.

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

Scanner Invariant Multiple Sclerosis Lesion Segmentation from MRI

Authors: Shahab Aslani, Vittorio Murino, Michael Dayan, Roger Tam, Diego Sona, Ghassan Hamarneh

Abstract: This paper presents a simple and effective generalization method for magnetic resonance imaging (MRI) segmentation when data is collected from multiple MRI scanning sites and as a consequence is affected by (site-)domain shifts. We propose to integrate a traditional encoder-decoder network with a regularization network. This added network includes an auxiliary loss term which is responsible for th… ▽ More This paper presents a simple and effective generalization method for magnetic resonance imaging (MRI) segmentation when data is collected from multiple MRI scanning sites and as a consequence is affected by (site-)domain shifts. We propose to integrate a traditional encoder-decoder network with a regularization network. This added network includes an auxiliary loss term which is responsible for the reduction of the domain shift problem and for the resulting improved generalization. The proposed method was evaluated on multiple sclerosis lesion segmentation from MRI data. We tested the proposed model on an in-house clinical dataset including 117 patients from 56 different scanning sites. In the experiments, our method showed better generalization performance than other baseline networks. △ Less

Submitted 22 October, 2019; originally announced October 2019.

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

Evolution Of Post-Impact Remnant Helium Stars In Type Ia Supernova Remnants Within The Single-Degenerate Scenario

Authors: Kuo-Chuan Pan, Paul Ricker, Ronald Taam

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

Submitted 23 June, 2013; v1 submitted 5 March, 2013; originally announced March 2013.

Comments: 37 pages, 12 figures, accepted for publication in 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

Self-gravitational force calculation of infinitesimally thin gaseous disks

Authors: C. C. Yen, R. E. Taam, Ken H. C. Yeh, K. C. Chang

Abstract: A thin gaseous disk has often been investigated in the context of various phenomena in galaxies, which point to the existence of starburst rings and dense circumnuclear molecular disks. The effect of self-gravity of the gas in the 2D disk can be important in confronting observations and numerical simulations in detail. For use in such applications, a new method for the calculation of the gravitati… ▽ More A thin gaseous disk has often been investigated in the context of various phenomena in galaxies, which point to the existence of starburst rings and dense circumnuclear molecular disks. The effect of self-gravity of the gas in the 2D disk can be important in confronting observations and numerical simulations in detail. For use in such applications, a new method for the calculation of the gravitational force of a 2D disk is presented. Instead of solving the complete potential function problem, we calculate the force in infinite planes in Cartesian and polar coordinates by a reproducing kernel method. Under the limitation of a 2D disk, we specifically represent the force as a double summation of a convolution of the surface density and a fundamental kernel and employ a fast Fourier transform technique. In this method, the entire computational complexity can be reduced from $O(N^2\times N^2)$ to $O((N\times \log_2(N)^2)$, where $N$ is the number of zones in one dimension. This approach does not require softening. The proposed method is similar to a spectral method, but without the necessity of imposing a periodic boundary condition. We further show this approach is of near second order accuracy for a smooth surface density in a Cartesian coordinate system. △ Less

Submitted 6 October, 2012; originally announced October 2012.

Comments: 16 pages

Journal ref: Journal of Computational Physics, 231(2012) 8246-8261

Evolution of Post-Impact Companion Stars in SN Ia Remnants within the Single-Degenerate Scenario

Authors: Kuo-Chuan Pan, Paul Ricker, Ronald Taam

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

Submitted 15 October, 2012; v1 submitted 30 September, 2012; originally announced October 2012.

Comments: 36 pages, 11 figures, Accepted to ApJ; add a footnote for discussing the progenitor of SN1006

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

Common Envelope Evolution: Where we stand and how we can move forward

Authors: N. Ivanova, S. Justham, X. Chen, O. De Marco, C. L. Fryer, E. Gaburov, H. Ge, E. Glebbeek, Z. Han, X. -D. Li, G. Lu, T. Marsh, Ph. Podsiadlowski, A. Potter, N. Soker, R. Taam, T. M. Tauris, E. P. J. van den Heuvel, R. F. Webbink

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

Submitted 6 December, 2012; v1 submitted 19 September, 2012; originally announced September 2012.

Comments: 83 pages (11 chapters, 2 appendixes), 10 figures, Astronomy & Astrophysics Review format. This version contains revisions as a response to community comments

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

Impact of type ia supernova ejecta on binary companions in the single-degenerate scenario

Authors: Kuo-Chuan Pan, Paul Ricker, Ronald Taam

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

Submitted 29 March, 2012; v1 submitted 7 March, 2012; originally announced March 2012.

Comments: Accepted to ApJ

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