Source author record

Seok-Jun Chang

Seok-Jun Chang 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

8works
3topics
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

8 published item(s)

preprint2026arXiv

Rapid emergence of overmassive black holes in the early Universe

The origin of supermassive black holes (SMBHs) remains a long-standing problem in astrophysics. Recent JWST observations reveal an unexpectedly abundant population of overmassive black holes at z>4-6, where the BH masses lie far above local scaling relations and not reproduced by current cosmological models. How such overmassive black holes form and rapidly grow within young galaxies has remained unclear. Here we present fully cosmological radiation-hydrodynamic simulations that, for the first time, self-consistently follow the birth, early growth, and emergent observable signatures of SMBHs in proto-cluster environments. We find that heavy seeds of order $10^6 M_\text{sun}$ naturally form, exceeding typical theoretical expectations by an order of magnitude. These seeds rapidly develop dense, optically thick disks whose strong electron scattering produces broad H$α$ emission comparable to that seen in little red dots (LRDs). Sustained super-Eddington accretion then drives fast growth to $\sim 3 \times 10^7 ~M_\text{sun}$ by $z \sim 8$. These results provide a unified physical scenario in which LRDs correspond to a short-lived, enshrouded phase of heavy-seed formation, naturally evolving into the overmassive quasars detected by JWST and ultimately the progenitors of today's SMBHs.

preprint2025arXiv

Impact of Resonance, Raman, and Thomson Scattering on Hydrogen Line Formation in Little Red Dots

Little Red Dots (LRDs) are compact sources at $z>5$ discovered through JWST spectroscopy. Their spectra exhibit broad Balmer emission lines ($\gtrsim1000\rm~km~s^{-1}$), alongside absorption features and a pronounced Balmer break -- evidence for a dense, neutral hydrogen medium with the $n=2$ state. When interpreted as arising from AGN broad-line regions, inferred black hole masses from local scaling relations exceed expectations given their stellar masses, challenging models of early black hole-galaxy co-evolution. However, radiative transfer effects in dense media may also impact the formation of hydrogen emission lines. We model three scattering processes shaping hydrogen line profiles: resonance scattering by hydrogen in the $n=2$ state, Raman scattering of UV radiation by ground-state hydrogen, and Thomson scattering by free electrons. Using 3D Monte Carlo radiative transfer simulations with multi-branching resonance transitions, we examine their imprint on line shapes and ratios. Resonance scattering produces strong deviations from Case B flux ratios, clear differences between H$α$ and H$β$, and encodes gas kinematics in line profiles but cannot broaden H$β$ due to conversion to Pa$α$. While Raman scattering can yield broad wings, scattering of UV continuum is disfavored given the absence of strong FWHM variations across transitions. Raman scattering of higher Lyman-series emission can produce H$α$/H$β$ wing width ratios of $\gtrsim1.28$, agreeing with observations. Thomson scattering can reproduce the observed $\gtrsim1000~\rm km\, s^{-1}$ wings under plausible conditions, e.g., $T_{\rm e} \sim 10^4\rm \, K$ and $N_{\rm e}\sim10^{24}\rm~cm^{-2}$ -- and lead to black hole mass overestimates by factors $\gtrsim10$. Our results provide a framework for interpreting hydrogen lines in LRDs and similar systems.

preprint2021arXiv

Ly-alpha Radiative Transfer: A Stokes Vector Approach to Ly-alpha Polarization

Ly-alpha emitting galaxies and giant Ly-alpha blobs (LABs) have been extensively observed to study the formation history of galaxies. However, the origin of their extended Ly-alpha emission, especially of LABs, remains controversial. Polarization signals from some LABs have been discovered, and this is commonly interpreted as strong evidence supporting that the extended Ly-alpha emission originates from the resonance scattering. The Monte Carlo Ly-alpha radiative transfer code LaRT is updated to investigate the polarization of Ly-alpha using the Stokes vector formalism. We apply LaRT to a few models to explore the fundamental polarization properties of Ly-alpha. Interestingly, individual Ly-alpha photon packets are found to be almost completely polarized by a sufficient number of scatterings (N_scatt > 10^4-10^5 in a static medium) or Doppler shifts induced by gas motion, even starting from unpolarized light. It is also found that the polarization pattern can exhibit a non-monotonically increasing pattern in some cases, besides the commonly-known trend that the polarization monotonically increases with radius. The polarization properties are primarily determined by the degree of polarization of individual photon packets and the anisotropy of the Ly-alpha radiation field, which are eventually controlled by the medium's optical depth and velocity field. If once Ly-alpha photon packets achieve ~100% polarization, the radial profile of polarization appears to correlate with the surface brightness profile. A steep surface brightness profile tends to yield a rapid increase of the linear polarization near the Ly-alpha source location. In contrast, a shallow surface brightness profile gives rise to a slowly increasing polarization pattern.

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

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.