Source author record

Seokho Lee

Seokho 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)

preprint2025arXiv

$\ell_0$-Regularized Item Response Theory Model for Robust Ideal Point Estimation

Ideal point estimation methods face a significant challenge when legislators engage in protest voting -- strategically voting against their party to express dissatisfaction. Such votes introduce attenuation bias, making ideologically extreme legislators appear artificially moderate. We propose a novel statistical framework that extends the fast EM-based estimation approach of \cite{Imai2016} using $\ell_0$ regularization method to handle protest votes. Through simulation studies, we demonstrate that our proposed method maintains estimation accuracy even with high proportions of protest votes, while being substantially faster than MCMC-based methods. Applying our method to the 116th and 117th U.S. House of Representatives, we successfully recover the extreme liberal positions of ``the Squad'', whose protest votes had caused conventional methods to misclassify them as moderates. While conventional methods rank Ocasio-Cortez as more conservative than 69\% of Democrats, our method places her firmly in the progressive wing, aligning with her documented policy positions. This approach provides both robust ideal point estimates and systematic identification of protest votes, facilitating deeper analysis of strategic voting behavior in legislatures.

preprint2022arXiv

A new method for direct measurement of isotopologue ratios in protoplanetary disks: a case study of the $^{12}$CO/$^{13}$CO ratio in the TW Hya disk

Planetary systems are thought to be born in protoplanetary disks. Isotope ratios are a powerful tool for investigating the material origin and evolution from molecular clouds to planetary systems via protoplanetary disks. However, it is challenging to measure the isotope (isotopologue) ratios, especially in protoplanetary disks, because the emission lines of major species are saturated. We developed a new method to overcome these challenges by using optically thin line wings induced by thermal broadening. As a first application of the method, we analyzed two carbon monoxide isotopologue lines, $^{12}$CO $3-2$ and $^{13}$CO $3-2$, from archival observations of a protoplanetary disk around TW Hya with the Atacama Large Millimeter/sub-millimeter Array. The $^{12}$CO/$^{13}$CO ratio was estimated to be ${ 20\pm5}$ at disk radii of ${ 70-110}$ au, which is significantly smaller than the value observed in the local interstellar medium, $\sim69$. It implies that an isotope exchange reaction occurs in a low-temperature environment with $\rm C/O>1$ . In contrast, it is suggested that $^{12}$CO/$^{13}$CO is higher than $\sim{ 84}$ in the outer disk ($r > { 130}$ au), which can be explained by the difference in the binding energy of the isotopologues on dust grains and the CO gas depletion processes. Our results imply that the gas-phase $^{12}$CO/$^{13}$CO can vary by a factor of ${ > 4}$ even inside a protoplanetary disk, and therefore, can be used to trace material evolution in disks.

preprint2022arXiv

Detection of HC18O+ in a protoplanetary disk: exploring oxygen isotope fractionation of CO

The oxygen isotope fractionation scenario, which has been developed to explain the oxygen isotope anomaly in the solar system materials, predicts that CO gas is depleted in 18O in protoplanetary disks, where segregation between solids and gas inside disks had already occurred. Based on ALMA observations, we report the first detection of HC18O+(4-3) in a Class II protoplanetary disk (TW Hya). This detection allows us to explore the oxygen isotope fractionation of CO in the TW Hya disk from optically thin HCO+ isotopologues as a proxy of optically thicker CO isotopologues. Using the H13CO+(4-3) data previously obtained with SMA, we find that the H13CO+/HC18O+ ratio in the central <100 au regions of the disk is 10.3 +- 3.2. We construct a chemical model of the TW Hya disk with carbon and oxygen isotope fractionation chemistry, and estimate the conversion factor from H13CO+/HC18O+ to 13CO/C18O. With the conversion factor (= 0.8), the 13CO/C18O ratio is estimated to be 8.3 +- 2.6, which is consistent with the elemental abundance ratio in the local ISM (8.1 +- 0.8) within error margin. Then there is no clear evidence of 18O depletion in CO gas of the disk, although we could not draw any robust conclusion due to large uncertainties. In conclusion, optically thin lines of HCO+ isotopologues are useful tracers of CO isotopic ratios, which are hardly constrained directly from optically thick lines of CO isotopologues. Future higher sensitivity observations of H13CO+ and HC18O+ would be able to allow us to better constrain the oxygen fractionation in the disk.

preprint2020arXiv

Modeling Nitrogen Fractionation in the Protoplanetary Disk around TW Hya: Model Constraints on Grain Population and Carbon-to-Oxygen Elemental Abundance Ratio

Observations conducted using the Atacama Large Millimeter/submillimeter Array on the protoplanetary disk around TW Hya show the nitrogen fractionation of HCN molecules in HC$^{14}$N/HC$^{15}$N $\sim$120 at a radius of $\sim$20 AU. In this study, we investigated the physical and chemical conditions that control this nitrogen fractionation process. To this end, a new disk model was developed, in which the isotope-selective photodissociation of N$_2$ and isotope-exchange chemical reactions have been incorporated. Our model can successfully reproduce the observed HCN column density when the elemental abundances of the gas-phase carbon and oxygen are depleted by two orders of magnitude relative to those in the interstellar medium and carbon is more abundant than oxygen ([C/O]$_{\rm elem}>$ 1). The isotope-selective photodissociation of N$_2$ is the dominant nitrogen fractionation process in our models. The observed HC$^{14}$N/HC$^{15}$N ratio, which increases outwards, can also be reproduced by the model by assuming that the small dust grains in the atmosphere of the outer disk are depleted more than those in the inner disk. This is consistent with grain evolution models, according to which small dust grains are continuously replenished in the inner disk due to fragmentation of the large dust grains that radially drift from the outer disk.

preprint2019arXiv

The circumstellar environment around the embedded protostar EC 53

EC53 is an embedded protostar with quasi-periodic emission in the near-IR and sub-mm. We use ALMA high-resolution observations of continuum and molecular line emission to describe the circumstellar environment of EC 53. The continuum image reveals a disk with a flux that suggests a mass of 0.075 Msun, much less than the estimated mass in the envelope, and an in-band spectral index that indicates grain growth to centimeter sizes. Molecular lines trace the outflow cavity walls, infalling and rotating envelope, and/or the Keplerian disk. The rotation profile of the C17O 3-2 line emission cannot isolate the Keplerian motion clearly although the lower limit of the protostellar mass can be calculated as 0.3 +- 0.1 Msun if the Keplerian motion is adopted. The weak CH3OH emission, which is anti-correlated with the HCO+ 4-3 line emission, indicates that the water snow line is more extended than what expected from the current luminosity, attesting to bygone outburst events. The extended snow line may persist for longer at the disk surface because the lower density increases the freeze-out timescale of methanol and water.

preprint2019arXiv

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

preprint2016arXiv

IGRINS spectroscopy of Class I sources: IRAS 03445+3242 and IRAS 04239+2436

We have detected molecular and atomic line emission from the hot and warm disks of two Class I sources, IRAS 03445+3242 and IRAS 04239+2436 using the high resolution Immersion GRating INfrared Spectrograph (IGRINS). CO overtone band transitions and near-IR lines of Na I and Ca I, all in emission, trace the hot inner disk while CO rovibrational absorption spectra of the first overtone transition trace the warm gas within the inner few AU of the disk. The emission-line profiles for both sources show evidence for Keplerian disks. A thin Keplerian disk with power-law temperature and column density profiles with a projected rotational velocity of $\sim$60--75 km s$^{-1}$ and a gas temperature of $\sim$3500 K at the innermost annulus can reproduce the CO overtone band emission. Na I and Ca I emission lines also arise from this disk, but they show complicated line features possibly affected by photospheric absorption lines. Multi-epoch observations show asymmetric variations of the line profiles on one-year (CO overtone bandhead and atomic lines for IRAS 03445+3242) or on one-day (atomic lines for IRAS 04239+2436) time scales, implying non-axisymmetric features in disks. The narrow CO rovibrational absorption spectra ($v$=0$\rightarrow$2) indicate that both warm ($>$ 150 K) and cold ($\sim$20--30 K) CO gas are present along the line of sight to the inner disk. This study demonstrates the power of IGRINS as a tool for studies of the sub-AU scale hot and AU-scale warm protoplanetary disks with its simultaneous coverage of the full H and K bands with high spectral resolution ($R$= 45,000) allowing many aspects of the sources to be investigated at once.

preprint2015arXiv

"Dust, Ice, and Gas In Time" (DIGIT) Herschel Observations of GSS30-IRS1 in Ophiuchus

As a part of the "Dust, Ice, and Gas In Time" (DIGIT) key program on Herschel, we observed GSS30-IRS1, a Class I protostar located in Ophiuchus (d = 120 pc), with Herschel/Photodetector Array Camera and Spectrometer (PACS). More than 70 lines were detected within a wavelength range from 50 micron to 200 micron, including CO, H2O, OH, and two atomic [O I] lines at 63 and 145 micron. The [C II] line, known as a tracer of externally heated gas by the interstellar radiation field, is also detected at 158 micron. All lines, except [O I] and [C II], are detected only at the central spaxel of 9.4" X 9.4". The [O I] emissions are extended along a NE-SW orientation, and the [C II] line is detected over all spaxels, indicative of external PDR. The total [C II] intensity around GSS30 reveals that the far-ultraviolet radiation field is in the range of 3 to 20 G0, where G0 is in units of the Habing Field, 1.6 X 10^{-3} erg cm^{-2} s^{-1}. This enhanced external radiation field heats the envelope of GSS30-IRS1, causing the continuum emission to be extended, unlike the molecular emission. The best-fit continuum model of GSS30-IRS1 with the physical structure including flared disk, envelope, and outflow shows that the internal luminosity is 10 Lsun, and the region is externally heated by a radiation field enhanced by a factor of 130 compared to the standard local interstellar radiation field.

preprint2015arXiv

Simultaneous observation of water and class I methanol masers toward class II methanol maser sources

We present a simultaneous single-dish survey of 22 GHz water maser and 44 GHz and 95 GHz class I methanol masers toward 77 6.7 GHz class II methanol maser sources, which were selected from the Arecibo methanol maser Galactic plane survey (AMGPS) catalog.Water maser emission is detected in 39 (51%) sources, of which 15 are new detections. Methanol maser emission at 44 GHz and 95 GHz is found in 25 (32%) and 19 (25%) sources, of which 21 and 13 sources are newly detected, respectively. We find 4 high-velocity (> 30 km/s) water maser sources, including 3 dominant blue- or redshifted outflows.The 95 GHz masers always appear with the 44 GHz maser emission. They are strongly correlated with 44 GHz masers in velocity, flux density, and luminosity, while they are not correlated with either water or 6.7 GHz class II methanol masers. The average peak flux density ratio of 95 GHz to 44 GHz masers is close to unity, which is two times higher than previous estimates. The flux densities of class I methanol masers are more closely correlated with the associated BGPS core mass than those of water or class II methanol masers. Using the large velocity gradient (LVG) model and assuming unsaturated class I methanol maser emission, we derive the fractional abundance of methanol to be in a range of 4.2*10^-8 to 2.3*10^-6, with a median value of 3.3\pm2.7*10^-7.

preprint2015arXiv

The warm CO gas along the UV-heated outflow cavity walls: a possible interpretation for the Herschel/PACS CO spectra of embedded YSOs

A fraction of the mid-$J$ ($J$= 14--13 to $J$= 24--23) CO emission detected by the \textit{Herschel}/PACS observations of embedded young stellar objects (YSOs) has been attributed to the UV-heated outflow cavity walls. We have applied our newly developed self-consistent models of Photon-Dominated Region (PDR) and non-local thermal equilibrium line Radiative transfer In general Grid (RIG) to the \textit{Herschel} FIR observations of 27 low mass YSOs and one intermediate mass YSO, NGC7129-FIRS2. When the contribution of the hot component (traced by transitions of $J> 24$) is removed, the rotational temperature of the warm component is nearly constant with $\sim250$ K. This can be reproduced by the outflow cavity wall ($n \geq 10^6\, \mathrm{cm}^{-3}$, $\log G_{0}/n \geq-4.5$, $\mathrm{log} G_0\ge 3$, $T_{\rm gas} \ge 300 $K, and X(CO)$ \ge 10^{-5}$) heated by a UV radiation field with a black body temperature of 15,000 K or 10,000 K. However, a shock model combined with an internal PDR will be required to determine the quantitative contribution of a PDR relative to a shock to the mid-$J$ CO emission.

preprint2014arXiv

A PDR model for the FIR mid-$J$ CO ladder with universal rotational temperature in star forming region

A photon dominated region (PDR) is one of the leading candidate mechanisms for the origin of the warm CO gas with near universal $\sim$300~K rotational temperature inferred from the CO emission detected towards embedded protostars by Herschel/PACS. We have developed a PDR model in general coordinates, where we can use the most adequate coordinate system for an embedded protostar having outflow cavity walls, to solve chemistry and gas energetics self-consistently for given UV radiation fields with different spectral shapes. Simple 1D tests and applications show that FIR mid-$J$ ($14\leq J\leq 24$) CO lines are emitted from near the surface of a dense region exposed to high UV fluxes. We apply our model to HH46 and find the UV-heated outflow cavity wall can reproduce the mid-$J$, CO transitions observed by Herschel/PACS. A model with UV radiation corresponding to a blackbody of 10,000~K results in the rotational temperature lower than 300~K, while models with the Draine interstellar radiation field and the 15,000 K blackbody radiation field predict the rotational temperature similar to the observed one.

preprint2014arXiv

Herschel Key Program, "Dust, Ice, and Gas In Time" (DIGIT): the origin of molecular and atomic emission in low-mass protostars in Taurus

Six low-mass embedded sources (L1489, L1551-IRS5, TMR1, TMC1-A, L1527, and TMC1) in Taurus have been observed with Herschel-PACS to cover the full spectrum from 50 to 210 $μ$m as part of the Herschel key program, "Dust, Ice, and Gas In Time (DIGIT)". The relatively low intensity of the interstellar radiation field surrounding Taurus minimizes contamination of the [C II] emission associated with the sources by diffuse emission from the cloud surface, allowing study of the [C II] emission from the source. In several sources, the [C II] emission is distributed along the outflow, as is the [O I] emission. The atomic line luminosities correlate well with each other, as do the molecular lines, but the atomic and molecular lines correlate poorly. The relative contribution of CO to the total gas cooling is constant at $\sim$30 %, while the cooling fraction by H$_2$O varies from source to source, suggesting different shock properties resulting in different photodissociation levels of H$_2$O. The gas with a power-law temperature distribution with a moderately high density can reproduce the observed CO fluxes, indicative of CO close to LTE. However, H$_2$O is mostly subthermally excited. L1551-IRS5 is the most luminous source ($L_{\rm bol}=24.5$ L$_\odot$) and the [O I] 63.1 $μ$m line accounts for more than 70 % of its FIR line luminosity, suggesting complete photodissociation of H$_2$O by a J-shock. In L1551-IRS5, the central velocity shifts of the [O I] line, which exceed the wavelength calibration uncertainty ($\sim$70 km s$^{-1}$) of PACS, are consistent with the known red- and blue-shifted outflow direction.

preprint2013arXiv

L1448-MM observations by the Herschel Key program, "Dust, Ice, and Gas In Time" (DIGIT)

We present Herschel/PACS observations of L1448-MM, a Class 0 protostar with a prominent outflow. Numerous emission lines are detected at 55 < lambda < 210 micrometer including CO, OH, H2O, and [O I]. We investigate the spatial distribution of each transition to find that lines from low energy levels tend to distribute along the outflow direction while lines from high energy levels peak at the central spatial pixel. Spatial maps reveal that OH emission lines are formed in a relatively small area, while [O I] emission is extended. According to the rotational diagram analysis, the CO emission can be fitted by two (warm and hot) temperature components. For H2O, the ortho-to-para ratio is close to 3. The non-LTE LVG calculations suggest that CO and H2O lines could instead be formed in a high kinetic temperature (T > 1000 K) environment, indicative of a shock origin. For OH, IR-pumping processes play an important role in the level population. The molecular emission in L1448-MM is better explained with a C-shock model, but the atomic emission of PACS [O I] and Spitzer/IRS [Si II] emission is not consistent with C-shocks, suggesting multiple shocks in this region. Water is the major line coolant of L1448-MM in the PACS wavelength range, and the best-fit LVG models predict that H2O and CO emit (50-80)% of their line luminosity in the PACS wavelength range.

preprint2010arXiv

Sparse logistic principal components analysis for binary data

We develop a new principal components analysis (PCA) type dimension reduction method for binary data. Different from the standard PCA which is defined on the observed data, the proposed PCA is defined on the logit transform of the success probabilities of the binary observations. Sparsity is introduced to the principal component (PC) loading vectors for enhanced interpretability and more stable extraction of the principal components. Our sparse PCA is formulated as solving an optimization problem with a criterion function motivated from a penalized Bernoulli likelihood. A Majorization--Minimization algorithm is developed to efficiently solve the optimization problem. The effectiveness of the proposed sparse logistic PCA method is illustrated by application to a single nucleotide polymorphism data set and a simulation study.