Source author record

Jongsoo Kim

Jongsoo Kim 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

25works
6topics
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

25 published item(s)

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2021arXiv

The intrinsic structure of Sagittarius A* at 1.3 cm and 7 mm

Sagittarius A* (Sgr A*), the Galactic Center supermassive black hole (SMBH), is one of the best targets to resolve the innermost region of SMBH with very long baseline interferometry (VLBI). In this study, we have carried out observations toward Sgr A* at 1.349 cm (22.223 GHz) and 6.950 mm (43.135 GHz) with the East Asian VLBI Network, as a part of the multi-wavelength campaign of the Event Horizon Telescope (EHT) in 2017 April. To mitigate scattering effects, the physically motivated scattering kernel model from Psaltis et al. (2018) and the scattering parameters from Johnson et al. (2018) have been applied. As a result, a single, symmetric Gaussian model well describes the intrinsic structure of Sgr A* at both wavelengths. From closure amplitudes, the major-axis sizes are ~704$\pm$102 $μ$as (axial ratio $\sim$1.19$^{+0.24}_{-0.19}$) and $\sim$300$\pm$25 $μ$as (axial ratio $\sim$1.28$\pm$0.2) at 1.349 cm and 6.95 mm respectively. Together with a quasi-simultaneous observation at 3.5 mm (86 GHz) by Issaoun et al. (2019), we show that the intrinsic size scales with observing wavelength as a power-law, with an index $\sim$1.2$\pm$0.2. Our results also provide estimates of the size and compact flux density at 1.3 mm, which can be incorporated into the analysis of the EHT observations. In terms of the origin of radio emission, we have compared the intrinsic structures with the accretion flow scenario, especially the radiatively inefficient accretion flow based on the Keplerian shell model. With this, we show that a nonthermal electron population is necessary to reproduce the source sizes.

preprint2019arXiv

Kinematics of the M87 jet in the collimation zone: gradual acceleration and velocity stratification

We study the kinematics of the M87 jet using the first year data of the KVN and VERA Array (KaVA) large program, which has densely monitored the jet at 22 and 43 GHz since 2016. We find that the apparent jet speeds generally increase from $\approx0.3c$ at $\approx0.5$ mas from the jet base to $\approx2.7c$ at $\approx20$ mas, indicating that the jet is accelerated from subluminal to superluminal speeds on these scales. We perform a complementary jet kinematic analysis by using archival Very Long Baseline Array monitoring data observed in $2005-2009$ at 1.7 GHz and find that the jet is moving at relativistic speeds up to $\approx5.8c$ at distances of $200-410$ mas. We combine the two kinematic results and find that the jet is gradually accelerated over a broad distance range that coincides with the jet collimation zone, implying that conversion of Poynting flux to kinetic energy flux takes place. If the jet emission consists of a single streamline, the observed trend of jet acceleration ($Γ\propto z^{0.16\pm0.01}$) is relatively slow compared to models of a highly magnetized jet. This indicates that Poynting flux conversion through the differential collimation of poloidal magnetic fields may be less efficient than expected. However, we find a non-negligible dispersion in the observed speeds for a given jet distance, making it difficult to describe the jet velocity field with a single power-law acceleration function. We discuss the possibility that the jet emission consists of multiple streamlines following different acceleration profiles, resulting in jet velocity stratification.

preprint2015arXiv

A New Hardware Correlator in Korea: Performance Evaluation using KVN observations

We report results of the performance evaluation of a new hardware correlator in Korea, the Daejeon correlator, developed by the Korea Astronomy and Space Science Institute (KASI) and the National Astronomical Observatory of Japan (NAOJ). We conducted Very Long Baseline Interferometry (VLBI) observations at 22~GHz with the Korean VLBI Network (KVN) in Korea and the VLBI Exploration of Radio Astrometry (VERA) in Japan, and correlated the aquired data with the Daejeon correlator. For evaluating the performance of the new hardware correlator, we compared the correlation outputs from the Daejeon correlator for KVN observations with those from a software correlator, the Distributed FX (DiFX). We investigated the correlated flux densities and brightness distributions of extragalactic compact radio sources. The comparison of the two correlator outputs show that they are consistent with each other within $<8\%$, which is comparable with the amplitude calibration uncertainties of KVN observations at 22~GHz. We also found that the 8\% difference in flux density is caused mainly by (a) the difference in the way of fringe phase tracking between the DiFX software correlator and the Daejeon hardware correlator, and (b) an unusual pattern (a double-layer pattern) of the amplitude correlation output from the Daejeon correlator. The visibility amplitude loss by the double-layer pattern is as small as 3\%. We conclude that the new hardware correlator produces reasonable correlation outputs for continuum observations, which are consistent with the outputs from the DiFX software correlator.

preprint2015arXiv

Amplitude Correction Factors of KVN Observations

We report results of investigation of amplitude calibration for very long baseline interferometry (VLBI) observations with Korean VLBI Network (KVN). Amplitude correction factors are estimated based on comparison of KVN observations at 22~GHz correlated by Daejeon hardware correlator and DiFX software correlator in Korea Astronomy and Space Science Institute (KASI) with Very Long Baseline Array (VLBA) observations at 22~GHz by DiFX software correlator in National Radio Astronomy Observatory (NRAO). We used the observations for compact radio sources, 3C 454.3, NRAO 512, OJ 287, BL Lac, 3C 279, 1633+382, and 1510-089, which are almost unresolved for baselines in a range of 350-477 km. Visibility data of the sources obtained with similar baselines at KVN and VLBA are selected, fringe-fitted, calibrated, and compared for their amplitudes. We found that visibility amplitudes of KVN observations should be corrected by factors of 1.10 and 1.35 when correlated by DiFX and Daejeon correlators, respectively. These correction factors are attributed to the combination of two steps of 2-bit quantization in KVN observing systems and characteristics of Daejeon correlator.

preprint2015arXiv

Dynamical evolution of supernova remnants breaking through molecular clouds

We carry out three-dimensional hydrodynamic simulations of the supernova remnants (SNRs) produced inside molecular clouds (MCs) near their surface using the HLL code (Harten et al. 1983). We explore the dynamical evolution and the X-ray morphology of SNRs after breaking through the MC surface for ranges of the explosion depths below the surface and the density ratios of the clouds to the intercloud media (ICM). We find that if an SNR breaks out through an MC surface in its Sedov stage, the outermost dense shell of the remnant is divided into several layers. The divided layers are subject to the Rayleigh-Taylor instability and fragmented. On the other hand, if an SNR breaks through an MC after the remnant enters the snowplow phase, the radiative shell is not divided to layers. We also compare the predictions of previous analytic solutions for the expansion of SNRs in stratified media with our onedimensional simulations. Moreover, we produce synthetic X-ray surface brightness in order to research the center-bright X-ray morphology shown in thermal composite SNRs. In the late stages, a breakout SNR shows the center-bright X-ray morphology inside an MC in our results. We apply our model to the observational results of the X-ray morphology of the thermal composite SNR 3C 391.

preprint2015arXiv

High-resolution CO observation of the carbon star CIT 6 revealing the spiral structure and a nascent bipolar outflow

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$).

preprint2015arXiv

Power spectra of outflow-driven turbulence

We investigate the power spectra of outflow-driven turbulence through high-resolution three-dimensional isothermal numerical simulations where the turbulence is driven locally in real-space by a simple spherical outflow model. The resulting turbulent flow saturates at an average Mach number of ~2.5 and is analysed through density and velocity power spectra, including an investigation of the evolution of the solenoidal and compressional components. We obtain a shallow density power spectrum with a slope of ~-1.2 attributed to the presence of a network of localised dense filamentary structures formed by strong shock interactions. The total velocity power spectrum slope is found to be ~-2.0, representative of Burgers shock dominated turbulence model. The density weighted velocity power spectrum slope is measured as ~-1.6, slightly less than the expected Kolmogorov scaling value (slope of -5/3) found in previous works. The discrepancy may be caused by the nature of our real space driving model and we suggest there is no universal scaling law for supersonic compressible turbulence. We find that on average, solenoidal modes slightly dominate in our turbulence model as the interaction between strong curved compressible shocks generates solenoidal modes, and compressible modes decay faster.

preprint2015arXiv

The East-Asian VLBI Network

The East-Asian VLBI Network (EAVN) is the international VLBI facility in East Asia and is conducted in collaboration with China, Japan, and Korea. The EAVN consists of VLBI arrays operated in each East Asian country, containing 21 radio telescopes and three correlators. The EAVN will be mainly operated at 6.7 (C-band), 8 (X-band), 22 (K-band), and 43 GHz (Q-band), although the EAVN has an ability to conduct observations at 1.6 - 129 GHz. We have conducted fringe test observations eight times to date at 8 and 22 GHz and fringes have been successfully detected at both frequencies. We have also conducted science commissioning observations of 6.7 GHz methanol masers in massive star-forming regions. The EAVN will be operational from the second half of 2017, providing complementary results with the FAST on AGNs, massive star-forming regions, and evolved stars with high angular resolution at cm- to mm-wavelengths.

preprint2014arXiv

Astrometrically Registered Simultaneous Observations of the 22 GHz H$_2$O and the 43GHz SiO masers towards R Leonis Minoris using KVN and Source/Frequency Phase Referencing

Oxygen-rich Asymptotic Giant Branch (AGB) stars can be intense emitters of SiO ($v$=1 and 2, J=1$\rightarrow$0) and H$_2$O maser lines at 43 and 22 GHz, respectively. VLBI observations of the maser emission provide a unique tool to probe the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned \water\ and $v$=1 and $v$=2 SiO maser maps have traditionally limited the physical constraints that can be placed on the SiO maser pumping mechanism. We present phase referenced simultaneous spectral-line VLBI images for the SiO $v$=1 and $v$=2, J=1$\rightarrow$0, and H$_2$O maser emission around the AGB star R\,LMi, obtained from the Korean VLBI Network (KVN). The simultaneous multi-channel receivers of the KVN offer great possibilities for astrometry in the frequency domain. With this facility we have produced images with bona-fide absolute astrometric registration between high frequency maser transitions of different species to provide the positions of the \water\ maser emission, and the centre of the SiO maser emission, and hence reducing the uncertainty in the proper motion for R\,LMi by an order of magnitude over that from Hipparcos. This is the first successful demonstration of source frequency phase referencing for mm-VLBI spectral-line observations and also where the ratio between the frequencies is not an integer.

preprint2014arXiv

Magnetic Field and Faraday Rotation Measure in the Turbulent Warm Ionized Medium

Wu et al. (2009) suggested an empirical relation between the magnetic field strength along the line of sight (LOS) and the dispersion of Faraday rotation measure (RM) distribution in turbulent media with root-mean-square sonic Mach number $M_s \simeq 1$. In this paper, we extend the work by incorporating the Mach number dependence. Media with $\sim 0.5 < M_s < \sim 2$ are considered to cover the Mach number range of the warm ionized medium (WIM) of our Galaxy. Three-dimensional, magnetohydrodynamic isothermal turbulence simulations with solenoidal forcing are used. We suggest a new relation among the LOS magnetic field strength, the dispersion of RM distribution, and the Mach number, which approximately represents the relation for Alfvénic perturbations. In addition, we suggest a relation between the Mach number and the dispersion of log-normal distribution of emission measure (EM), which is basically the relation for the Mach number and the density dispersion. The relations could be used for a quick and rough estimation of the LOS magnetic field strength in the turbulent WIM.

preprint2014arXiv

Testing assumptions and predictions of star-formation theories

(Abridged). We present numerical simulations of isothermal, MHD, supersonic turbulence, designed to test various hypotheses frequently assumed in star formation(SF) theories. We consider three simulations, each with a different combination of physical size, rms sonic Mach number, and Jeans parameter, but chosen as to give the same value of the virial parameter and to conform with Larson's scaling relations. As in the non-magnetic case: we find no simultaneously subsonic and super-Jeans structures in our MHD simulations. We find that the fraction of small-scale super-Jeans structures increases when self gravity is turned on, and that the production of gravitationally unstable dense cores by turbulence alone is very low. This implies that self-gravity is in general necessary not only to induce the collapse of Jeans-unstable cores, but also to form them. We find that denser regions tend to have more negative values of the velocity divergence, implying a net inwards flow towards the regions' centers. We compare the results from our simulations with the predictions from the recent SF theories by Krumholz & McKee, Padoan & Nordlund, and Hennebelle & Chabrier, using the expressions recently provided by Federrath & Klessen. We find that none of these theories reproduces the dependence of the SFEff with Ms observed in our simulations in the MHD case. The SFEff predicted by the theories ranges between half and one order of magnitude larger than what we observe in the simulations in both the HD and the MHD cases. We conclude that the type of flow used in simulations like the ones presented here and assumed in recent SF theories, may not correctly represent the flow within actual clouds, and that theories that assume it does may be missing a fundamental aspect of the flow. We suggest that a more realistic regime may be that of hierarchical gravitational collapse.

preprint2014arXiv

The First Very Long Baseline Interferometry Image of 44 GHz Methanol Maser with the KVN and VERA Array (KaVA)

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

preprint2014arXiv

Unbiased water and methanol maser surveys of NGC 1333

We present the results of unbiased 22 GHz H2O water and 44 GHz class I CH3OH methanol maser surveys in the central 7x10 arcmin area of NGC 1333 and two additional mapping observations of a 22 GHz water maser in a ~3x3arcmin area of the IRAS4A region. In the 22 GHz water maser survey of NGC 1333 with sensitivity of sigma~0.3Jy, we confirmed masers toward H2O(B) in the region of HH 7-11 and IRAS4B. We also detected new water masers at ~20arcsec away in the western direction of IRAS4B or ~25arcsec away in the southern direction of IRAS4A. We could not however find young stellar objects or molecular outflows associated with them. They showed two different velocity components of ~0 and ~16 km/s, which are blue- and red-shifted relative to the adopted systemic velocity of ~7 km/s for NGC 1333. They also showed time variabilities in both intensity and velocity from multi-epoch observations and an anti-correlation between the intensities of the blue- and the red-shifted velocity components. We suggest that the unidentified powering source of these masers might be in the earliest evolutionary stage of star formation before the onset of molecular outflows. Finding this kind of water masers is only possible by an unbiased blind survey. In the 44 GHz methanol maser survey with sensitivity of sigma~0.5 Jy, we confirmed masers toward the IRAS4A2 and the eastern shock region of the IRAS2A. Both sources are also detected in 95 and 132 GHz methanol maser lines. In addition, we had new detections of methanol masers at 95 and 132 GHz toward IRAS4B. In terms of the isotropic luminosity, we detected the methanol maser sources brighter than ~5x1025 erg/s from our unbiased survey.

preprint2014arXiv

Verification of the Astrometric Performance of the Korean VLBI Network, using comparative SFPR studies with the VLBA at 14/7 mm

The Korean VLBI Network (KVN) is a new mm-VLBI dedicated array with capability for simultaneous observations at multiple frequencies, up to 129 GHz. The innovative multi-channel receivers present significant benefits for astrometric measurements in the frequency domain. The aim of this work is to verify the astrometric performance of the KVN using a comparative study with the VLBA, a well established instrument. For that purpose, we carried out nearly contemporaneous observations with the KVN and the VLBA, at 14/7 mm, in April 2013. The KVN observations consisted of simultaneous dual frequency observations, while the VLBA used fast frequency switching observations. We used the Source Frequency Phase Referencing technique for the observational and analysis strategy. We find that having simultaneous observations results in a superior performance for compensation of all atmospheric terms in the observables, in addition to offering other significant benefits for astrometric analysis. We have compared the KVN astrometry measurements to those from the VLBA. We find that the structure blending effects introduce dominant systematic astrometric shifts and these need to be taken into account. We have tested multiple analytical routes to characterize the impact of the low resolution effects for extended sources in the astrometric measurements. The results from the analysis of KVN and full VLBA datasets agree within 2-$σ$ of the thermal error estimate. We interpret the discrepancy as arising from the different resolutions. We find that the KVN provides astrometric results with excellent agreement, within 1-$σ$, when compared to a VLBA configuration which has a similar resolution. Therefore this comparative study verifies the astrometric performance of KVN using SFPR at 14/7 mm, and validates the KVN as an astrometric instrument.

preprint2014arXiv

VLBI observations of bright AGN jets with KVN and VERA Array (KaVA): Evaluation of Imaging Capability

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

preprint2013arXiv

Density distributions of outflow driven turbulence

Protostellar jets and outflows are signatures of star formation and promising mechanisms for driving supersonic turbulence in molecular clouds. We quantify outflow-driven turbulence through three-dimensional numerical simulations using an isothermal version of the robust total variation diminishing code. We drive turbulence in real-space using a simplified spherical outflow model, analyse the data through density probability distribution functions (PDF), and investigate the Core Formation Rate per free-fall time (CFR_ff). The real-space turbulence driving method produces a negatively skewed density PDF possessing an enhanced tail on the low-density side. It deviates from the log-normal distributions typically obtained from Fourier-space turbulence driving at low densities, but can provide a good fit at high-densities, particularly in terms of mass weighted rather than volume weighted density PDF. Due to this fact, we suggest that the CFR_ff determined from a Fourier-driven turbulence model could be comparable to that of our particular real-space driving model, which has a ratio of solenoidal to compressional components from the resulting turbulence velocity fields of ~0.6.

preprint2013arXiv

Simulated Faraday Rotation Measures toward High Galactic Latitudes

We study the Faraday rotation measure (RM) due to the Galactic magnetic field (GMF) toward high Galactic latitudes. The RM arises from the global, regular component as well as from the turbulent, random component of the GMF. We model the former based on observations and the latter using the data of magnetohydrodynamic turbulence simulations. For a large number of different GMF models, we produce mock RM maps around the Galactic poles and calculate various statistical quantities with the RM maps. We find that the observed medians of RMs toward the north and south Galactic poles, $\sim 0.0\pm 0.5 {\rm rad m^{-2}}$ and $\sim +6.3\pm 0.5 {\rm rad m^{-2}}$, are difficult to explain with any of our many alternate GMF models. The standard deviation of observed RMs, $\sim 9 {\rm rad m^{-2}}$, is clearly larger than that of simulated RMs. The second-order structure function of observed RMs is substantially larger than that of simulated RMs, especially at small angular scales. We discuss other possible contributions to RM toward high Galactic latitudes. Besides observational errors and the intrinsic RM of background radio sources against which RM is observed, we suggest that the RM due to the intergalactic magnetic field may account for a substantial fraction of the observed RM. Finally we note that reproducing the observed medians may require additional components or/and structures of the GMF that are not present in our models.

preprint2013arXiv

The Density Distribution in Turbulent Bi-stable Flows

We numerically study the volume density probability distribution function (n-PDF) and the column density probability distribution function (Sigma-PDF) resulting from thermally bistable turbulent flows. We analyze three-dimensional hydrodynamic models in periodic boxes of 100pc by side, where turbulence is driven in the Fourier space at a wavenumber corresponding to 50pc. At low densities (n <= 0.6cm^-3) the n-PDF, is well described by a lognormal distribution for average local Mach number ranging from ~0.2 to ~5.5. As a consequence of the non linear development of thermal instability (TI), the logarithmic variance of the distribution for the diffuse gas increases with M faster than in the well known isothermal case. The average local Mach number for the dense gas (n >= 7.1cm^-3) goes from ~1.1 to ~16.9 and the shape of the high density zone of the n-PDF changes from a power-law at low Mach numbers to a lognormal at high M values. In the latter case the width of the distribution is smaller than in the isothermal case and grows slower with M. At high column densities the Sigma-PDF is well described by a lognormal for all the Mach numbers we consider and, due to the presence of TI, the width of the distribution is systematically larger than in the isothermal case but follows a qualitatively similar behavior as M increases. Although a relationship between the width of the distribution and M can be found for each one of the cases mentioned above, these relations are different form those of the isothermal case.

preprint2012arXiv

Early science with Korean VLBI network: the QCAL-1 43GHz calibrator survey

This paper presents the catalog of correlated flux densities in three ranges of baseline projection lengths of 637 sources from a 43 GHz (Q-band) survey observed with the Korean VLBI Network. Of them, 623 sources have not been observed before at Q-band with VLBI. The goal of this work in the early science phase of the new VLBI array is twofold: to evaluate the performance of the new instrument that operates in a frequency range of 22-129 GHz and to build a list of objects that can be used as targets and as calibrators. We have observed the list of 799 target sources with declinations down to -40 degrees. Among them, 724 were observed before with VLBI at 22 GHz and had correlated flux densities greater than 200 mJy. The overall detection rate is 78%. The detection limit, defined as the minimum flux density for a source to be detected with 90% probability in a single observation, was in a range of 115-180 mJy depending on declination. However, some sources as weak as 70 mJy have been detected. Of 623 detected sources, 33 objects are detected for the first time in VLBI mode. We determined their coordinates with the median formal uncertainty 20 mas. The results of this work set the basis for future efforts to build the complete flux-limited sample of extragalactic sources at frequencies 22 GHz and higher at 3/4 of the celestial sphere.

preprint2010arXiv

Density Power Spectrum in Turbulent Thermally Bi-stable Flows

In this paper we numerically study the behavior of the density power spectrum in turbulent thermally bistable flows. We analyze a set of five three-dimensional simulations where turbulence is randomly driven in Fourier space at a fixed wave-number and with different Mach numbers M (with respect to the warm medium) ranging from 0.2 to 4.5. The density power spectrum becomes shallower as M increases and the same is true for the column density power spectrum. This trend is interpreted as a consequence of the simultaneous turbulent compressions, thermal instability generated density fluctuations, and the weakening of thermal pressure force in diffuse gas. This behavior is consistent with the fact that observationally determined spectra exhibit different slopes in different regions. The values of the spectral indexes resulting from our simulations are consistent with observational values. We do also explore the behavior of the velocity power spectrum, which becomes steeper as M increases. The spectral index goes from a value much shallower than the Kolmogorov one for M=0.2 to a value steeper than the Kolmogorov one for M=4.5.

preprint2010arXiv

Enhanced Core Formation Rate in a Turbulent Cloud by Self-gravity

We performed a numerical experiment designed for core formation in a self-gravitating, magnetically supercritical, supersonically turbulent, isothermal cloud. A density probability distribution function (PDF) averaged over a converged turbulent state before turning self-gravity on is well-fitted with a lognormal distribution. However, after turning self-gravity on, the volume fractions of density PDFs at a high density tail, compared with the lognormal distribution, increase as time goes on. In order to see the effect of self-gravity on core formation rates, we compared the core formation rate per free-fall time (CFR$_{\rm ff}$) from the theory based on the lognormal distribution and the one from our numerical experiment. For our fiducial value of a critical density, 100, normalised with an initial value, the latter CFR$_{\rm ff}$ is about 30 times larger the former one. Therefore, self-gravity plays an important role in significantly increasing CFR$_{\rm ff}$. This result implies that core (star) formation rates or core (stellar) mass functions predicted from theories based on the lognormal density PDF need some modifications. Our result of the increased volume fraction of density PDFs after turning self-gravity on is consistent with power-law like tails commonly observed at higher ends of visual extinction PDFs of active star-forming clouds.

preprint2009arXiv

Estimation of Magnetic Field Strength in the Turbulent Warm Ionized Medium

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

preprint2007arXiv

The Virial Balance of Clumps and Cores in Molecular Clouds

We study the instantaneous virial balance of clumps and cores (CCs) in 3D simulations of driven, MHD, isothermal molecular clouds (MCs). The models represent a range of magnetic field strengths in MCs from subcritical to non-magnetic regimes. We identify CCs at different density thresholds, and for each object, we calculate all the terms that enter the Eulerian form of the virial theorem (EVT). A CC is considered gravitationally bound when the gravitational term in the EVT is larger than the amount for the system to be virialized, which is more stringent than the condition that it be large enough to make the total volume energy negative. We also calculate, quantities commonly used in the observations to indicate the state of gravitational boundedness of CCs such as the Jeans number J_c, the mass-to magnetic flux ratio mu_c, and the virial parameter alpha_vir. Our results show that: a) CCs are dynamical out-of-equilibrium structures. b) The surface energies are of the same order than their volume counterparts c) CCs are either in the process of being compressed or dispersed by the velocity field. Yet, not all CCs that have a compressive net kinetic energy are gravitationally bound. d) There is no 1-to-1 correspondence between the state of gravitational boundedness of a CC as described by the virial analysis or as implied by the classical indicators. In general, in the virial analysis, we observe that only the inner regions of the objects are gravitationally bound, whereas J_c, alpha_vir, and mu_c estimates tend to show that they are more gravitationally bound at the lowest threshold levels and more magnetically supercritical. g) We observe, in the non-magnetic simulation, the existence of a bound core with structural and dynamical properties that resemble those of the Bok globule Barnard 68 (B68).

preprint2004arXiv

The Lifetimes and Evolution of Molecular Cloud Cores

We discuss the lifetimes and evolution of clumps and cores formed as turbulent density fluctuations in nearly isothermal molecular clouds. In the non-magnetic case, clumps are unlikely to reach a hydrostatic state, and instead are expected to either proceed directly to collapse, or else ``rebound'' towards the mean pressure and density of the parent cloud. Rebounding clumps are delayed in their re-expansion by their self-gravity. From a simple virial calculation, we find re-expansion times of a few free-fall times. In the magnetic case, we present a series of driven-turbulence, ideal-MHD isothermal numerical simulations in which we follow the evolution of clumps and cores in relation to the magnetic criticality of their ``parent clouds'' (the numerical boxes). In subcritical boxes, magnetostatic clumps do not form. A few moderately-gravitationally bound clumps form which however are dispersed by the turbulence in < 1.3 Myr. An estimate of the ambipolar diffusion (AD) time scale t_AD in these cores gives t_AD > 1.3 Myr, only slightly longer than the dynamical times. In supercritical boxes, some cores become locally supercritical and collapse in typical times ~ 1 Myr. We also observe longer-lived supercritical cores that however do not collapse because they are smaller than the local Jeans length. Fewer clumps and cores form in these simulations than in their non-magnetic counterpart. Our results suggest that a) A fraction of the cores may not form stars, and may correspond to some of the observed starless cores. b) Cores may be out-of-equilibrium structures, rather than quasi-magnetostatic ones. c) The magnetic field may help reduce the star formation efficiency by reducing the probability of core formation, rather than by significantly delaying the collapse of individual cores.