Source author record

Hee-Won Lee

Hee-Won Lee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2022arXiv

Formation of the Asymmetric Accretion Disk from Stellar Wind Accretion in an S-type Symbiotic Star

The accretion process in a typical S-type symbiotic star, targeting~AG Draconis, is investigated through 3D hydrodynamical simulations using the FLASH code. Regardless of the wind velocity of the giant star, an accretion disk surrounding the white dwarf is always formed. In the wind models faster than the orbital velocity of the white dwarf, the disk size and accretion rate are consistent with the predictions under the Bondi-Hoyle-Lyttleton (BHL) condition. In slower wind models, unlike the BHL predictions, the disk size does not grow and the accretion rate increases to a considerably higher level, up to $>20\%$ of the mass-loss rate of the giant star. The accretion disk in our fiducial model is characterized by a flared disk with a radius of 0.16~au and a scale height of 0.03 au. The disk mass of $\sim 5 \times 10^{-8} M_\odot$ is asymmetrically distributed with the density peak toward the giant star, being about $50\%$ higher than the density minimum in the disk. Two inflowing spiral features are clearly identified and their relevance to the azimuthal asymmetry of disk is pointed out. The flow in the accretion disk is found to be sub-Keplerian with about $90\%$ of the Keplerian speed, which indicates the caveat of overestimating the O VI emission region from spectroscopy of Raman-scattered O VI features at 6825 Å and 7082 Å.

preprint2020arXiv

3D Grid-Based Monte Carlo Code for Radiative Transfer through Raman and Rayleigh Scattering with Atomic Hydrogen -- STaRS

Emission features formed through Raman scattering with atomic hydrogen provide unique and crucial information to probe the distribution and kinematics of a thick neutral region illuminated by a strong far UV emission source. We introduce a new 3 dimensional Monte-Carlo code in order to describe the radiative transfer of line photons that are subject to Raman and Rayleigh scattering with atomic hydrogen. In this code entitled "${\bf S}$ejong Radiative ${\bf T}$r${\bf a}$nsfer through ${\bf R}$aman and Rayleigh ${\bf S}$cattering (${\it STaRS}$), each photon is traced until escape with a tag attached carrying information including the position, direction, wavelength, and polarization. The thick neutral scattering region is divided into numerous cells with each cell being characterized by its velocity and density, which ensures huge flexibility of the code in analyzing Raman-scattered features formed in a neutral region with complicated kinematics and density distribution. As a test of the code, we revisit the formation of Balmer wings through Raman scattering of far UV continuum near Ly$β$ and Ly$γ$ in a static neutral region. An additional check is made to investigate Raman scattering of O VI in an expanding neutral medium. We find fairly good agreement of our results with previous works, demonstrating the capability of dealing with radiative transfer modeling that can be applied to spectropolarimetric imaging observations of various objects including symbiotic stars, young planetary nebulae, and active galactic nuclei.

preprint2019arXiv

Line Formation of Raman-Scattered He II $λ$ 4851 in an Expanding Spherical H I Shell in Young Planetary Nebulae

We investigate line formation of Raman-scattered He II at 4851 in an expanding neutral spherical shell that surrounds a point-like He II source located at the center. A new grid-based Monte Carlo code is used to take into consideration the H I density variation along each photon path. In the case of a monochromatic He II emission source, the resultant line profiles are characterized by an asymmetric double peak structure with a tertiary peak and a significant red tail that may extend to line centers of He II$λ$4859 and H$β$. The peak separation corresponds to the expansion velocity, which we consider is in the range $20-40{\rm\ km\ s^{-1}}$ in this work. Tertiary red peaks are formed as a result of multiple Rayleigh reflections at the inner surface of a hollow spherical shell of \ion{H}{1}. Due to a sharp increase of scattering cross section near resonance, the overall Raman conversion efficiency is significantly enhanced as the expansion speed increases. In the case of a He II line source with a Gaussian line profile with a full width at half maximum of $30 - 70{\rm\ km\ s^{-1}}$, we obtain distorted redward profiles due to increasing redward cross section of H I. A simple application to the young planetary nebula IC 5117 is consistent with a neutral shell expanding with a speed $\sim 30{\rm\ km\ s^{-1}}$.

preprint2016arXiv

A Monte Carlo Study of Flux Ratios of Raman Scattered O~VI Features at 6825 Å and 7082 Å in Symbiotic Stars

Symbiotic stars are regarded as wide binary systems consisting of a hot white dwarf and a mass losing giant. They exhibit unique spectral features at 6825 Å and 7082 Å, which are formed via Raman scattering of \ion{O}{6}$λλ$ 1032 and 1038 with atomic hydrogen. We adopt a Monte Carlo technique to generate the same number of \ion{O}{6}$λ$1032 and $λ$1038 line photons and compute the flux ratio $F(6825)/F(7082)$ of these Raman scattered \ion{O}{6} features formed in neutral regions with a simple geometric shape as a function of \ion{H}{1} column density $N_{HI}$. In cylindrical and spherical neutral regions with the \ion{O}{6} source embedded inside, the flux ratio $F(6825)/F(7082)$ shows an overall decrease from 3 to 1 as $N_{HI}$ increases in the range $10^{22-24}{\rm\ cm^{-2}}$. In the cases of a slab geometry and other geometries with the \ion{O}{6} source outside the \ion{H}{1} region, Rayleigh escape operates to lower the flux ratio considerably. For moderate values of $N_{HI}\sim 10^{23}{\rm\ cm^{-2}}$ the flux ratio behaves in a complicated way to exhibit a broad bump with a peak value of 3.5 in the case of a sphere geometry. We find that the ratio of Raman conversion efficiencies of \ion{O}{6}$λλ$1032, 1038 ranges from 0.8 to 3.5. Our high resolution spectra of 'D' type HM~Sge and 'S' type AG~Dra obtained with the {\it Canada-France-Hawaii-Telescope} show that the flux ratio $F(6825)/F(7082)$ of AG~Dra is significantly smaller than that of HM~Sge, implying that 'S' type symbiotics are characterized by higher $N_{HI}$ than 'D' type symbiotics.

preprint2016arXiv

A Profile Analysis of Raman-scattered O VI Bands at 6825 Å and 7082 Å in Sanduleak's Star

We present a detailed modeling of the two broad bands observed at 6825 Å and 7082 Å in Sanduleak's star, a controversial object in the Large Magellanic Cloud. These bands are known to originate from Raman-scattering of O VI $λλ$ 1032 and 1038 photons with atomic hydrogen and are only observed in bona fide symbiotic stars. Our high-resolution spectrum obtained with the Magellan Inamori Kyocera Echelle (MIKE) spectrograph at the Magellan-Clay Telescope reveals, quite surprisingly, that the profiles of the two bands look very different: while the Raman 6825 Å band shows a single broad profile with a redward extended bump, the Raman 7082 Å band exhibits a distinct triple-peak profile. Our model suggests that the O VI emission nebula can be decomposed into a red, blue and central emission regions from an accretion disk, a bipolar outflow and a further compact, optically thick region. We also perform Monte Carlo simulations with the aim of fitting the observed flux ratio $F(6825)/F(7082) \sim 4.5$, which indicate that the neutral region in Sanduleak's star is characterized by the column density $N_{HI} \sim 1 \times 10^{23} {\rm\ cm^{-2}}$.

preprint2016arXiv

Polarization of Rayleigh scattered Lyα in active galactic nuclei

The unification scheme of active galactic nuclei (AGNs) invokes an optically thick molecular torus component hiding the broad emission line region. Assuming the presence of a thick neutral component in the molecular torus characterized by a \ion{H}{I} column density > $10^{22}{\rm\ cm^{-2}}$, we propose that far UV radiation around Ly$α$ can be significantly polarized through Rayleigh scattering. Adopting a Monte Carlo technique we compute polarization of Rayleigh scattered radiation near Ly$α$ in a thick neutral region in the shape of a slab and a cylindrical shell. It is found that radiation near Ly$α$ Rayleigh reflected from a very thick slab can be significantly polarized in a fairly large range of wavelength $Δλ\sim 50$ Å exhibiting a flux profile similar to the incident one. Rayleigh transmitted radiation in a slab is characterized by the central dip with a complicated polarization behavior. The optically thick part near Ly$α$ center is polarized in the direction perpendicular to the slab normal, which is in contrast to weakly polarized wing parts in the direction parallel to the slab normal. A similar polarization flip phenomenon is also found in the case of a tall cylindrical shell, in which the spatial diffusion along the vertical direction near the inner cylinder wall for core photons leads to a tendency of the electric field aligned to the direction perpendicular to the vertical axis. Observational implications are briefly discussed including spectropolarimetry of the quasar PG~1630+377 by Koratkar et al. in 1990 where Ly$α$ is strongly polarized with no other emission lines polarized.

preprint2015arXiv

Formation of Raman Scattering Wings around H alpha, H beta and Pa alpha in Active Galactic Nuclei

Powered by a supermassive black hole with an accretion disk, the spectra of active galactic nuclei (AGNs) are characterized by prominent emission lines including Balmer lines. The unification schemes of AGNs require the existence of a thick molecular torus that may hide the broad emission line region from the view of observers near the equatorial direction. In this configuration, one may expect that the far UV radiation from the central engine can be Raman scattered by neutral hydrogen to reappear around Balmer and Paschen emission lines which can be identified with broad wings. We produce H$α$, H$β$ and Pa$α$ wings using a Monte Carlo technique to investigate their properties. The neutral scattering region is assumed to be a cylindrical torus specified by the inner and outer radii and the height. While the covering factor of the scattering region affects the overall strengths of the wings, the wing widths are primarily dependent on the neutral hydrogen column density $N_{\rm HI}$ being roughly proportional to $N_{\rm HI}^{1/2}$. In particular, with $N_{\rm HI}=10^{23}{\rm\ cm^{-2}}$ the H$α$ wings typically show a width $\sim 2\times 10^4{\rm\ km\ s^{-1}}$. We also find that H$α$ and Pa$α$ wing profiles are asymmetric with the red part stronger than the blue part and an opposite behavior is seen for H$β$ wings.

preprint2014arXiv

Raman Scattered Ne VII$λ$973 at 4881 Å in the Symbiotic Star V1016 Cygni

We present the high resolution spectra of the symbiotic star V1016 Cygni obtained with the Bohyunsan Optical Echelle Spectrograph in 2003 and 2005, from which we find a broad emission feature at 4881 Å. We propose that this broad feature is formed from Raman scattering of Ne VII$\lambda973$ by atomic hydrogen. Thus far, the detection of Raman scattered lines by atomic hydrogen is limited to O VI$λλ$1032, 1038 and He II$λλ$940, 972 and 1025. With the adoption of the center wavelength 973.302 Å of Ne VII$λ$973 and consideration of the air refractive index of $n_{air}=1.000279348$, the atomic line center of the Raman scattered Ne VII feature is determined to be 4880.53 Å. The total cross section at the line center of Ne VII$\lambda973$ is computed to be $2.62\times 10^{-22}{\rm\ cm^2}$ with the branching ratio of 0.17. We perform Monte Carlo simulations to fit the Raman scattered Ne VII$\lambda973$. Assuming that the Ne VII and He II emission regions share the same kinematics with respect to the neutral scattering region, we find that the Raman scattered He II$λ$972 at 4850 Å and Ne VII$λ$973 at 4881 Å are excellently fitted. We also propose that the He II and Ne VII emission regions were stationary with respect to the H I region in 2003 but that they were receding from it with a velocity $\sim 20{\rm\ km\ s^{-1}}$ in 2005.

preprint2013arXiv

Asymmetric Absorption Profiles of Lyα and Lyβ in Damped Lyman Alpha Systems

Damped Lyα systems (DLAs) observed in the quasar spectra are characterized by a high neutral hydrogen column density $N_{HI}>2\times 10^{20}{\rm\ cm^{-2}}$. The absorption wing profiles are often fitted using the Voigt function due to the fact that the scattering cross section near resonant line center is approximately described by the Lorentzian function. Since a hydrogen atom has infinitely many $p$ states that participate in the electric dipole interaction, the cross section starts to deviate from the Lorentzian in an asymmetric way in the line wing regions. We investigate this asymmetry in the absorption line profiles around Lyα and Lyβ as a function of the neutral hydrogen column density $N_{HI}$. In terms of $Δλ\equivλ-λ_α$ we expand the Kramers-Heisenberg formula around Lyα to find $σ(λ)\simeq (0.5 f_{12})^2σ_{T}(Δλ/λ_α)^{-2}[1+3.792 (Δλ/λ_α)]$, where $f_{12}$ and $σ_T$ are the oscillator strength of Ly$α$ and the Thomson scattering cross section, respectively. In terms of $Δλ_2\equivλ-λ_β$ in the vicinity of Lyβ, the total scattering cross section, given as the sum of cross sections for Rayleigh and Raman scattering, is shown to be $σ(λ)\simeq σ_T(0.5f_{13})^2(1+R_0) (Δλ_2/λ_β)^{-2}[1-24.68(Δλ_2/λ_β)]$ with $f_{13}$ and the factor $R_0=0.1342$ being the oscillator strength for Lyβ and the ratio of Raman cross section to Rayleigh cross section, respectively. A redward asymmetry develops around Lyα whereas a blue asymmetry is obtained for Lyβ. The absorption center shifts are found to be almost proportional to the neutral hydrogen column density.

preprint2009arXiv

Raman Scattered HeII 6545 in the Young and Compact Planetary Nebula NGC 6790

We present the high resolution spectra of the youn and compact planerary nebula NGC 6790 obtained with the echelle spectrograph at Bohyunsan Optical Astronomy Observatory and report the discovery of Raman scattered He II 6545 in this object. This line feature is formed in a thick neutral region surrounding the hot central star, where He II 1025 line photons are scattered inelastically by hydrogen atoms. A Monte Carlo technique is adopted to compute the line profile with a simple geometric model, in which the neutral region is in the form of a cylindrical shell that is expanding from the central star. From our line profile analysis, the expansion velocity of the HI region lies in the range V_exp = 15- 19 Km/s. Less stringent constraints are put on the HI column density N_HI and covering factor C, where the total flux of Raman He II 6545 is consistent with the product of N_HI and C being CN_HI ~ 0.5 X 10^20 cm^-2. The Monte Carlo profiles from stationary emission models exhibit deficit in the wing parts. A much better fit is obtained when the He II emission region is assumed to take the form of a ring that slowly rotates with a rotation speed ~ 18Km/s. Brief discussions are presented regarding the mass loss processes and future observations.

preprint2000arXiv

Lyman alpha Transfer in a thick, dusty, and static medium

We developed a Monte Carlo code that describes the resonant Lyman alpha line transfer in an optically thick, dusty, and static medium. The code was tested against the analytic formula derived by Neufeld (1990). We explain the line transfer mechanism for a wide range of line center optical depths by tracing histories of photons in the medium. We find that photons escape from the medium by a series of wing scatterings, during which polarization may develop. We applied our code to examine the amount of dust extinction around the Lyman alpha in primeval galaxies. Brief discussions on the astrophysical application of our work are presented.

preprint1998arXiv

Lyman alpha line formation in an expanding H I superbubble in the primeval galaxy DLA 2233+131

Various types of galaxies observed in the cosmological scales show P-Cygni type profiles in the Lyman alpha emission lines. A Monte Carlo code is developed to investigate the Lyman alpha line transfer in an optically thick and moving medium with a careful consideration of the scattering in the damping wing. The main features in emergent line profiles include a primary emission peak in the red part and a much weaker secondary peak in the blue part. The primary peak recedes to the red and the width of the feature increases as $N_{HI}$ increases. The P-Cygni type profile in the Lyman alpha emission line of DLA 2233+131 ($z=3.15$) is noted to be similar to those found in the nearby galaxies and distant ($3<z<3.5$) H II galaxies. Our numerical results are applied to show that the DLA system may possess an expanding H I supershell with bulk flow of $\sim 200 km/s$ and that the H I column density $N_{HI}$ is approximately $10^{20} cm^{-2}$. From the observed Lyman alpha flux and adopting a typical size of the emission region $\sim 1 \kpc$, we estimate the electron density of the H II region to be $\sim 1 \cm^{-3} and the mass of H II region $\sim 10^8 \msun$. We also conclude that it requires $\sim 10^3$ O5 stars for photoionization, which is comparable to first-ranked H II regions found in nearby spiral and irregular galaxies.

preprint1998arXiv

On the Polarization of H-alpha Lines Scattered by Neutral Hydrogen in Active Galactic Nuclei

Raman scattering by atomic hydrogen converts the UV continuum around Ly$β$ into optical continuum around H$α$, and the basic atomic physics has been discussed in several works on symbiotic stars. We propose that the same process may operate in active galactic nuclei (AGN) and calculate the linear polarization of the broad emission lines Raman-scattered by a high column neutral hydrogen compnent. The conversion efficiency of the Raman scattering process is discussed and the expected scattered flux is computed using the spectral energy distribution of an AGN given by a typical power law. The high column H {\sc i} component in AGN is suggested by many observations encompassing radio through UV and X-ray ranges. When the neutral hydrogen component with a column density $\sim 10^{22} cm^{-2}$ is present around the active nucleus, it is found that the scattered H$α$ is characterized by a very broad width $\sim 20,000 km/s$ and that the strength of the polarized flux is comparable to that of the electron-scattered flux expected from a conventional unified model of narrow line AGN. The width of the scattered flux is mainly determined by the column density of the neutral scatterers where the total scattering optical depth becomes of order unity. The asymmetry in the Raman scattering cross section around Ly$β$ introduces red asymmetric polarized profiles around H$α$. The effects of the blended Ly$β$ and O {\sc vi} 1034 doublet are also investigated. We briefly discuss the spectropolarimetric observations performed on the Seyfert galaxy IRAS 110548-1131 and the narrow line radio galaxy Cyg A.

preprint1997arXiv

On the Profiles and Polarization of Raman Scattered Emission Lines in Symbiotic Stars:II. Numerical Simulations

A Monte Carlo method is used to calculate the profiles and the polarization of the Raman scattered O VI lines(lambda lambda 6827,7088) in symbiotic stars. A point-like isotropic UV radiation source is assumed and a simple spherical wind model is adopted for the kinematics of the scattering material from the cool giant. We first investigate the case where the incident line photons are described by a Gaussian profile having a width of 10^4 K. We subsequently investigate the effects of the extended ionized region and non-spherical wind models including a disk-type wind and a bipolar wind. The cases where the emissin source is described by non-Gaussian profiles are briefly studied. Finally, as an additional component for the kinematics of symbiotic stars the orbital motion of the hot component around the cool giant is included and the effect on the spectropolarimetry is investigated. In this case the polarization direction changes around the red part of the Raman-scattered emission lines, when the observer's line of sight is perpendicular to the orbital plane, and no such effect is seen when the line of sight lies in the orbital plane. Furthermore, complex peak structures are seen in the degree of polarization and polarized flux, which have often been observed in several symbiotic systems including RR Tel. Brief observational consequences and preditions are discussed in relation to the present and future spectropolarimetry for symbiotic stars. It is concluded that spectropolarimetry may provide a powerful diagnostic of the physical conditions of symbiotic stars.