Source author record

Nanase Harada

Nanase Harada 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

18works
4topics
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

18 published item(s)

preprint2024arXiv

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

preprint2022arXiv

AGN-driven Cold Gas Outflow of NGC 1068 Characterized by Dissociation-Sensitive Molecules

Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16 well-detected molecular line intensity maps convolved to a common 150 pc resolution. In addition, we include the [SIII]/[SII] line ratio and [CI] $^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show remarkable resemblance with that of a kpc-scale biconical outflow from the central AGN. We identify two prominent features: (1) central concentration at the circumnuclear disk (CND) and (2) two peaks across the center that coincide with the biconical outflow peaks. The concentrated molecular lines in the CND are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line emissions from molecules known to be enhanced in irradiated interstellar medium, CN, C$_2$H, and HNC, show similar concentrations and extended components along the bicone, suggesting that molecule dissociation is a dominant chemical effect of the cold molecular outflow of this galaxy. Although further investigation should be made, this scenario is consistent with the faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons and electrons.

preprint2022arXiv

Energizing Star Formation: The Cosmic Ray Ionization Rate in NGC 253 Derived From ALCHEMI Measurements of H$_3$O$^+$ and SO

The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidence for a high CRIR value; $10^3-10^6$ times the average cosmic ray ionization rate within the Milky Way. This is a broad constraint and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations towards the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H$_3$O$^+$ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 due to emission from their transitions. We find that each of the four regions modelled has a CRIR in the range $(1-80)\times10^{-14}$ s$^{-1}$ and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays including C$_2$H, HCO$^+$, HOC$^+$, and CO. From shock and PDR/XDR models, we further find that neither UV/X-ray driven nor shock dominated chemistry are a viable single alternative as none of these processes can adequately fit the abundances of all of these species.

preprint2022arXiv

The Kiloparsec-scale Neutral Atomic Carbon Outflow in the Nearby Type-2 Seyfert Galaxy NGC 1068: Evidence for Negative AGN Feedback

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback.

preprint2021arXiv

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: I. ALMA Observations at $λ= $1.4-0.4 mm and Continuum Analysis

We observed with ALMA three deeply buried nuclei in two galaxies, NGC 4418 and Arp 220, at $\sim$0.2$''$ resolution over a total bandwidth of 67 GHz in $f_{\rm rest}$ = 215 - 697 GHz. Here we (1) introduce our program, (2) describe our data reduction method for wide-band, high-resolution imaging spectroscopy, (3) analyze in visibilities the compact nuclei with line forests, (4) develop a continuum-based estimation method of dust opacity and gas column density in heavily obscured nuclei, which uses the BGN (buried galactic nuclei) model and is sensitive to $\log(N_{\rm H_2}/{\rm cm}^{-2}) \sim $ 25 - 26 at $λ\sim 1$ mm, and (5) present the continuum data and diagnosis of our targets. The three continuum nuclei have major-axis FWHM of $\sim$0.1$''$-0.3$''$ (20-140 pc) aligned to their rotating nuclear disks of molecular gas. However, each nucleus is described better with two or three concentric components than with a single Gaussian. The innermost cores have sizes of 0.05$''$-0.10$''$ (8-40 pc), peak brightness temperatures of ~100-500 K at 350 GHz, and more fractional flux at lower frequencies. The intermediate components correspond to the nuclear disks. They have axial ratios of $\approx$0.5 and hence inclinations $\stackrel{>}{\sim} 60$ deg. The outermost elements include the bipolar outflow from Arp 220W. We estimate 1 mm dust opacity of $τ_{\rm d,1mm} \approx 2.2$, $1.2$, and $\stackrel{<}{\sim} 0.4$ respectively for NGC 4418, Arp 220W, and Arp 220E. The first two correspond to $\log(N_{\rm H}/{\rm cm}^{-2}) \sim 26$ for conventional dust-opacity laws, and hence the nuclei are highly Compton thick.

preprint2021arXiv

Deeply Buried Nuclei in the Infrared-Luminous Galaxies NGC 4418 and Arp 220: II. Line Forests at $λ= $1.4--0.4 mm and Circumnuclear Gas Observed with ALMA

We present the line observations in our ALMA imaging spectral scan toward three deeply buried nuclei in NGC 4418 and Arp 220. We cover 67 GHz in $f_{\rm rest}$=215-697 GHz at about 0.2$"$ (30, 80 pc) resolution. All the nuclei show dense line forests; we report our initial line identification using 55 species. The line velocities generally indicate gas rotation around each nucleus, tracing nuclear disks of $\sim$100 pc sizes. We confirmed the counter-rotation of the nuclear disks in Arp 220 and that of the nuclear disk and the galactic disk in NGC 4418. While the brightest lines exceed 100 K, most of the major lines and many $^{13}$C isotopologues show absorption against even brighter continuum cores of the nuclei. The lines with higher upper-level energies, including those from vibrationally-excited molecules, tend to arise from smaller areas, indicating radially varying conditions in these nuclei. The outflows from the two Arp 220 nuclei cause blueshifted line absorption below the continuum level. The absorption mostly has small spatial offsets from the continuum peaks to indicate the outflow orientations. The bipolar outflow from the western nucleus is also imaged in multiple emission lines, showing the extent of $\sim$1$"$ (400 pc). Redshifted line absorption against the nucleus of NGC 4418 indicates either an inward gas motion or a small collimated outflow slanted to the nuclear disk. We also resolved some previous confusions due to line blending and misidentification.

preprint2020arXiv

ALMA Observations of Multiple-CO and C Lines Toward the Active Galactic Nucleus of NGC 7469: X-Ray-dominated Region Caught in the Act

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 2-1), $^{12}$CO($J$ = 3-2), $^{13}$CO($J$ = 2-1), and [CI]($^3P_1$-$^3P_0$) emission lines around the type 1 active galactic nucleus (AGN) of NGC 7469 ($z = 0.0164$) at $\sim 100$ pc resolutions. The CO lines are bright both in the circumnuclear disk (central $\sim 300$ pc) and the surrounding starburst (SB) ring ($\sim 1$ kpc diameter), with two bright peaks on either side of the AGN. By contrast, the [CI]($^3P_1$-$^3P_0$) line is strongly peaked on the AGN. Consequently, the brightness temperature ratio of [CI]($^3P_1$-$^3P_0$) to $^{13}$CO(2-1) is $\sim 20$ at the AGN, as compared to $\sim 2$ in the SB ring. Our local thermodynamic equilibrium (LTE) and non-LTE models indicate that the enhanced line ratios (or CI enhancement) are due to an elevated C$^0$/CO abundance ratio ($\sim 3-10$) and temperature ($\sim 100-500$ K) around the AGN as compared to the SB ring (abundance ratio $\sim 1$, temperature $\lesssim 100$ K), which accords with the picture of the X-ray-dominated Region (XDR). Based on dynamical modelings, we also provide CO(1-0)-to- and [CI]($^3P_1$-$^3P_0$)-to-molecular mass conversion factors at the central $\sim 100$ pc of this AGN as $α_{\rm CO} = 4.1$ and $α_{\rm CI} = 4.4~M_\odot$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. Our results suggest that the CI enhancement is potentially a good marker of AGNs that could be used in a new submillimeter diagnostic method toward dusty environments.

preprint2020arXiv

The MALATANG Survey: Dense Gas and Star Formation from High Transition HCN and HCO+ maps of NGC253

To study the high-transition dense-gas tracers and their relationships to the star formation of the inner $\sim$ 2 kpc circumnuclear region of NGC253, we present HCN $J=4-3$ and HCO$^+ J=4-3$ maps obtained with the James Clerk Maxwell Telescope (JCMT). With the spatially resolved data, we compute the concentration indices $r_{90}/r_{50}$ for the different tracers. HCN and HCO$^+$ 4-3 emission features tend to be centrally concentrated, which is in contrast to the shallower distribution of CO 1-0 and the stellar component. The dense-gas fraction ($f_\text{dense}$, traced by the velocity-integrated-intensity ratios of HCN/CO and HCO$^+$/CO) and the ratio $R_\text{31}$ (CO 3-2/1-0) decline towards larger galactocentric distances, but increase with higher SFR surface density. The radial variation and the large scatter of $f_\text{dense}$ and $R_\text{31}$ imply distinct physical conditions in different regions of the galactic disc. The relationships of $f_\text{dense}$ versus $Σ_\text{stellar}$, and SFE$_\text{dense}$ versus $Σ_\text{stellar}$ are explored. SFE$_\text{dense}$ increases with higher $Σ_\text{stellar}$ in this galaxy, which is inconsistent with previous work that used HCN 1-0 data. This implies that existing stellar components might have different effects on the high-$J$ HCN and HCO$^+$ than their low-$J$ emission. We also find that SFE$_\text{dense}$ seems to be decreasing with higher $f_\text{dense}$, which is consistent with previous works, and it suggests that the ability of the dense gas to form stars diminishes when the average density of the gas increases. This is expected in a scenario where only the regions with high-density contrast collapse and form stars.

preprint2019arXiv

The Nuclear Filaments inside the Circumnuclear Disk in the Central 0.5 pc of the Galactic center

We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*.

preprint2015arXiv

ALMA Observations of the Submillimeter Dense Molecular Gas Tracers in the Luminous Type-1 Active Nucleus of NGC 7469

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $μ$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

preprint2015arXiv

Submillimeter-HCN Diagram for an Energy Diagnostics in the Centers of Galaxies

Compiling data from literature and the ALMA archive, we show enhanced HCN(4-3)/HCO$^+$(4-3) and/or HCN(4-3)/CS(7-6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN-enhancement). The number of sample galaxies is significantly increased from our previous work. We expect this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN-to-HCO$^+$ and HCN-to-CS in AGNs as compared to SB galaxies by a factor of a few to even $>$ 10 is a plausible explanation for the submillimeter HCN-enhancement. However, a counter argument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we could not fully discriminate these two scenarios at this moment due to insufficient amount of multi-transition, multi-species data, the former scenario equivalently claims for abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest it is difficult with conventional gas phase X-ray dominated region (XDR) ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral-neutral reactions that are efficiently activated at high temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.

preprint2014arXiv

A Multi-Transition Study of Molecules toward NGC 1068 based on High-Resolution Imaging Observations with ALMA

We present 0.8-mm band molecular images and spectra obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) toward one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068. Distributions of CO isotopic species ($^{13}$CO and C$^{18}$O) $\it{J}$ = 3--2, CN $\it{N}$ = 3--2 and CS $\it{J}$ = 7--6 are observed toward the circumnuclear disk (CND) and a part of the starburst ring with an angular resolution of $\sim$1.$^{\prime\prime}$3 $\times$ 1.$^{\prime\prime}$2. The physical properties of these molecules and shock-related molecules such as HNCO, CH$_{3}$CN, SO, and CH$_{3}$OH detected in the 3-mm band were estimated using rotation diagrams under the assumption of local thermodynamic equilibrium. The rotational temperatures of the CO isotopic species and the shock-related molecules in the CND are, respectively, 14--22 K and upper limits of 20--40 K. Although the column densities of the CO isotopic species in the CND are only from one-fifth to one-third of that in the starburst ring, those of the shock-related molecules are enhanced by a factor of 3--10 in the CND. We also discuss the chemistry of each species, and compare the fractional abundances in the CND and starburst ring with those of Galactic sources such as cold cores, hot cores, and shocked molecular clouds in order to study the overall characteristics. We find that the abundances of shock-related molecules are more similar to abundances in hot cores and/or shocked clouds than to cold cores. The CND hosts relatively complex molecules, which are often associated with shocked molecular clouds or hot cores. Because a high X-ray flux can dissociate these molecules, they must also reside in regions shielded from X-rays.

preprint2014arXiv

Distributions of molecules in the circumnuclear disk and surrounding starburst ring in the Seyfert galaxy NGC 1068 observed with ALMA

Sensitive observations with ALMA allow astronomers to observe the detailed distributions of molecules with relatively weak intensity in nearby galaxies. In particular, we report distributions of several molecular transitions including shock and dust related species ($^{13}$CO $J$ = 1--0, C$^{18}$O $J$ = 1--0, $^{13}$CN $N$ = 1--0, CS $J$ = 2--1, SO $J_N$ = 3$_2$--2$_1$, HNCO $J_{Ka,Kc}$ = 5$_{0,5}$--4$_{0,4}$, HC$_3$N $J$ = 11--10, 12--11, CH$_3$OH $J_K$ = 2$_K$--1$_K$, and CH$_3$CN $J_K$ = 6$_K$--5$_K$) in the nearby Seyfert 2 galaxy NGC 1068 observed with the ALMA early science program. The central $\sim$1 arcmin ($\sim$4.3 kpc) of this galaxy was observed in the 100 GHz region covering $\sim$96--100 GHz and $\sim$108--111 GHz with an angular resolution of $\sim4"\times2"$ (290 pc$\times$140 pc) to study the effects of an active galactic nucleus and its surrounding starburst ring on molecular abundances. Here, we present images and report a classification of molecular distributions into three main categories: (1) Molecules concentrated in the circumnuclear disk (CND) (SO $J_N$ = 3$_2$--2$_1$, HC$_3$N $J$ = 11--10, 12--11, and CH$_3$CN $J_K$ = 6$_K$--5$_K$), (2) Molecules distributed both in the CND and the starburst ring (CS $J$ = 2--1 and CH$_3$OH $J_K$ = 2$_K$--1$_K$), (3) Molecules distributed mainly in the starburst ring ($^{13}$CO $J$ = 1--0 and C$^{18}$O $J$ = 1--0). Since most of the molecules such as HC$_3$N observed in the CND are easily dissociated by UV photons and X-rays, our results indicate that these molecules must be effectively shielded. In the starburst ring, the relative intensity of methanol at each clumpy region is not consistent with those of $^{13}$CO, C$^{18}$O, and CS. This difference is probably caused by the unique formation and destruction mechanisms of CH$_3$OH.

preprint2013arXiv

Modeling The Molecular Composition in an AGN Disk

We use a high-temperature chemical network to derive the molecular abundances in axisymmetric accretion disk models around active galactic nuclei (AGNs) within 100 pc using simple radial and vertical density and temperature distributions motivated by more detailed physical models. We explore the effects of X-ray irradiation and cosmic ray ionization on the spatial distribution of the molecular abundances of CO, CN, CS, HCN, HCO+, HC3N, C2H, and c-C3H2 using a variety of plausible disk structures. These simple models have molecular regions with a layer of X-ray dominated regions, a midplane without the strong influence of X-rays, and a high-temperature region in the inner portion with moderate X-ray flux where families of polyynes (C$_{\rm n}$H$_{2}$) and cyanopolyynes can be enhanced. For the high midplane density disks we explore, we find that cosmic rays produced by supernovae do not significantly affect the regions unless the star formation efficiency significantly exceeds that of the Milky Way. We highlight molecular abundance observations and ratios that may distinguish among theoretical models of the density distribution in AGN disks. Finally, we assess the importance of the shock crossing time and the accretion time relative to the formation time for various chemical species. Vertical column densities are tabulated for a number of molecular species at both the characteristic shock crossing time and steady state. Although we do not attempt to fit any particular system or set of observations, we discuss our models and results in the context of the nearby AGN NGC 1068.

preprint2013arXiv

Molecules in the Circumnuclear Disk of the Galactic Center

Within a few parsecs around the central Black Hole Sgr A*, chemistry in the dense molecular cloud material of the circumnuclear disk (CND) can be affected by many energetic phenomena such as high UV-flux from the massive central star cluster, X-rays from Sgr A*, shock waves, and an enhanced cosmic-ray flux. Recently, spectroscopic surveys with the IRAM 30 meter and the APEX 12 meter telescopes of substantial parts of the 80--500 GHz frequency range were made toward selected positions in and near the CND. These datasets contain lines from the molecules HCN, HCO$^+$, HNC, CS, SO, SiO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$, H$_3$O$^+$ and others. We conduct Large Velocity Gradient analyses to obtain column densities and total hydrogen densities, $n$, for each species in molecular clouds located in the southwest lobe of CND. The data for the above mentioned molecules indicate 10$^5\,$cm$^{-3} \lesssim n <10^6\,$cm$^{-3}$, which shows that the CND is tidally unstable. The derived chemical composition is compared with a chemical model calculated using the UCL\_CHEM code that includes gas and grain reactions, and the effects of shock waves. Models are run for varying shock velocities, cosmic-ray ionization rates, and number densities. The resulting chemical composition is fitted best to an extremely high value of cosmic-ray ionization rate $ζ\sim 10^{-14}\,$s$^{-1}$, 3 orders of magnitude higher than the value in regular Galactic molecular clouds, if the pre-shock density is $n=10^5\,$cm$^{-3}$.

preprint2013arXiv

Submillimeter ALMA Observations of the Dense Gas in the Low-Luminosity Type-1 Active Nucleus of NGC 1097

We present the first 100 pc scale view of the dense molecular gas in the central ~ 1.3 kpc region of the type-1 Seyfert NGC 1097 traced by HCN (J=4-3) and HCO+ (J=4-3) lines afforded with ALMA band 7. This galaxy shows significant HCN enhancement with respect to HCO+ and CO in the low-J transitions, which seems to be a common characteristic in AGN environments. Using the ALMA data, we study the characteristics of the dense gas around this AGN and search for the mechanism of HCN enhancement. We find a high HCN (J=4-3) to HCO+ (J=4-3) line ratio in the nucleus. The upper limit of the brightness temperature ratio of HCN (v2=1^{1f}, J=4-3) to HCN (J=4-3) is 0.08, which indicates that IR pumping does not significantly affect the pure rotational population in this nucleus. We also find a higher HCN (J=4-3) to CS (J=7-6) line ratio in NGC 1097 than in starburst galaxies, which is more than 12.7 on the brightness temperature scale. Combined from similar observations from other galaxies, we tentatively suggest that this ratio appears to be higher in AGN-host galaxies than in pure starburst ones similar to the widely used HCN to HCO+ ratio. LTE and non-LTE modeling of the observed HCN and HCO+ lines using J=4-3 and 1-0 data from ALMA, and J=3-2 data from SMA, reveals a high HCN to HCO+ abundance ratio (5 < [HCN]/[HCO+] < 20: non-LTE analysis) in the nucleus, and that the high-J lines (J=4-3 and 3-2) are emitted from dense (10^{4.5} < n_H2 [/cc] < 10^6), hot (70 < Tkin [K] < 550) regions. Finally we propose that the high temperature chemistry is more plausible to explain the observed enhanced HCN emission in NGC 1097 than the pure gas phase PDR/XDR chemistry.

preprint2011arXiv

Carbon-Chain Species in Warm-up Models

In previous warm-up chemical models of the low-mass star-forming region L1527, we investigated the evolution of carbon-chain unsaturated hydrocarbon species when the envelope temperature is slightly elevated to $T\approx 30$ K. These models demonstrated that enhanced abundances of such species can be explained by gas-phase ion-molecule chemistry following the partial sublimation of methane from grain surfaces. We also concluded that the abundances of hydrocarbon radicals such as the C$_{\rm n}$H family should be further enhanced as the temperatures increase to higher values, but this conclusion stood in contrast with the lack of unambiguous detection of these species toward hot core and corino sources. Meanwhile, observational surveys have identified C$_2$H, C$_4$H, CH$_3$CCH, and CH$_3$OH toward hot corinos (especially IRAS 16293-2422) as well as towards L1527, with lower abundances for the carbon chain radicals and higher abundances for the other two species toward the hot corinos. In addition, the {\it Herschel Space Telescope} has detected the bare linear chain C$_3$ in 50 K material surrounding young high-mass stellar objects. To understand these new results, we revisit previous warm-up models with an augmented gas-grain network that incorporated reactions from a gas-phase network that was constructed for use with increased temperature up to 800 K. Some of the newly adopted reactions between carbon-chain species and abundant H$_2$ possess chemical activation energy barriers. The revised model results now better reproduce the observed abundances of unsaturated carbon chains under hot-corino (100 K) conditions and make predictions for the abundances of bare carbon chains in the 50 K regions observed by Herschel HIFI.

preprint2008arXiv

Modeling Carbon Chain Anions in L1527

The low-mass protostellar region L1527 is unusual because it contains observable abundances of unsaturated carbon-chain molecules including CnH radicals, H2Cn carbenes, cyanopolyynes, and the negative ions C4H- and C6H-, all of which are more associated with cold cores than with protostellar regions. Sakai et al. suggested that these molecules are formed in L1527 from the chemical precursor methane, which evaporates from the grains during the heat-up of the region. With the gas-phase osu.03.2008 network extended to include negative ions of the families Cn-, and CnH-, as well as the newly detected C3N-, we modeled the chemistry that occurs following methane evaporation at T~ 25-30 K. We are able to reproduce most of the observed molecular abundances in L1527 at a time of ~5000 yr. At later times, the overall abundance of anions become greater than that of electrons, which has an impact on many organic species and ions. The anion-to-neutral ratio in our calculation is in good agreement with observation for C6H- but exceeds the observed ratio by more than three orders of magnitude for C4H-. In order to explain this difference, further investigation is needed on the rate coefficients for electron attachment and other reactions regarding anions.