Parsec Scales Research Articles

Star formation in infrared dark clouds: Self-gravity and dynamics

The role played by gravity in the transfer of interstellar matter from molecular cloud scales to protostellar scales is still highly debated. Only detailed studies on the kinematics of large samples of star-forming clouds will settle the issue. We present new IRAM 30m observations of a sample of 27 infrared dark clouds covering a large range of sizes, masses, and aspect ratios. Preliminary results suggest that gravity is regulating the dynamical evolution of these clouds on a couple of parsec scales.

FIRST-EPOCH VLBA IMAGING OF 20 NEW TeV BLAZARS

We present Very Long Baseline Array (VLBA) images of 20 TeV blazars (HBLs) not previously well-studied on the parsec scale. Observations were made between August and December 2013, at a frequency of 8.4 GHz. These observations represent the first epoch of a VLBA monitoring campaign on these blazars, and they significantly increase the fraction of TeV HBLs studied with high-resolution imaging. The peak VLBI flux densities of these sources range from about 10 to 100 mJy/beam, and parsec-scale jet structure is detected in all sources. About half of the VLBI cores are resolved, with brightness temperature upper limits of a few times 10^10 K, and we find that a brightness temperature of about 2x10^10 K is consistent with the VLBI data for all but one of the sources. Such brightness temperatures do not require any relativistic beaming to reduce the observed value below commonly invoked intrinsic limits; however, the lack of detection of counter-jets does place a modest limit on the bulk Lorentz factor of greater than about 2. These data are thus consistent with a picture where weak-jet sources like the TeV HBLs develop significant velocity structures on parsec-scales. We also extend consideration to the full sample of TeV HBLs by combining the new VLBI data with VLBI and gamma-ray data from the literature. By comparing measured VLBI and TeV fluxes to samples with intrinsically uncorrelated luminosities generated by Monte Carlo simulations, we find a marginally significant correlation between the VLBI and TeV fluxes for the full TeV HBL sample.

Kiloparsec-scale radio emission in Seyfert and LINER galaxies

Seyfert and LINER galaxies are known to exhibit compact radio emission on $\sim$ 10 to 100 parsec scales, but larger Kiloparsec-Scale Radio structures (KSRs) often remain undetected in sub-arcsec high resolution observations. We investigate the prevalence and nature of KSRs in Seyfert and LINER galaxies using the 1.4 GHz VLA FIRST and NVSS observations. Our sample consists of 2651 sources detected in FIRST and of these 1737 sources also have NVSS counterparts. Considering the ratio of total to peak flux density ($\theta$ $=$ ${\rm (S_{\rm int}/S_{\rm peak})^{1/2}}$) as a parameter to infer the presence of extended radio emission we show that $\geq$ 30$\%$ of FIRST detected sources possess extended radio structures on scales larger than 1.0 kpc. The use of low-resolution NVSS observations help us to recover faint extended KSRs that are resolved out in FIRST observations and results in $\geq$ 42.5$\%$ KSR sources in FIRST-NVSS subsample. This fraction is only a lower limit owing to the combination of projection, resolution and sensitivity effects. Our study demonstrates that KSRs may be more common than previously thought and are found across all redshifts, luminosities and radio-loudness. The extranuclear radio luminosity of KSR sources is found to be positively correlated with the core radio luminosity as well as the [O~III] $\lambda$5007{\AA} line luminosity and this can be interpreted as KSRs being powered by AGN rather than star-formation. The distributions of the FIR-to-radio ratios and mid-IR colors of KSR sources are also consistent with their AGN origin. However, contribution from star-formation cannot be ruled out particularly in sources with low radio luminosities.

DIVERSE PROPERTIES OF INTERSTELLAR MEDIUM EMBEDDING GAMMA-RAY BURSTS AT THE EPOCH OF REIONIZATION

Analysis is performed on ultra-high resolution large-scale cosmological radiation-hydrodynamic simulations to, for the first time, quantify the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high density ($\ge 10^4$cm$^{-3}$) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density ($\sim 10^{-2.5}$cm$^{-3}$) high temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature and metallicity of intervening absorbers, vary strongly from sightline to sightline. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium of the host galaxy based on a very small number of (usually one) sightlines would be precarious. Utilizing high-z GRBs to probe interstellar medium and intergalactic medium should be undertaken properly taking into consideration the physical diversities of the interstellar medium.

THE BLOB CRASHES INTO THE MIRROR: MODELING THE EXCEPTIONAL γ-RAY FLARING ACTIVITY OF 3C 454.3 IN 2010 NOVEMBER

We focus on the exceptional flaring activity of 3C 454.3 in November 2010 and we discuss a theoretical framework addressing all data in their overall evolution. For two weeks the source has shown a plateau of enhanced GeV emission preceding a sudden major flare lasting about 3 days before decaying. The gamma-ray flare onset is abrupt (about 6 hours), and is characterized by a prominent "Compton dominance" with the GeV flux exceeding the pre-flare values by a factor of 4-5, whereas the optical and X-ray fluxes increased only by a factor 2. We explore two alternatives. Case 1, with high-energy emission originating within the BLR; and Case 2, with most of it produced outside. We show that Case 1 has considerable problems in explaining the whole set of multifrequency data. Case 2, instead, leads to a consistent and interesting interpretation based on the enhanced inverse Compton radiation that is produced as the jet crashes onto a mirror cloud positioned at few parsec from the BH. This model explains the gamma-ray vs. optical/X-ray behavior of 3C 454.3, including the otherwise puzzling phenomena such as the prominent "orphan" optical flare, and the enhanced line emission with no appreciable gamma-ray counterpart that preceded the GeV flare. It also accounts for the delayed onset of the latter on top of the long plateau. Our modelling of the exceptional 3C 454.3 gamma-ray flare shows that, while emission inside the canonical BLR is problematic, major and rapid variations can be produced at parsec scales with moderate bulk Lorentz factors $\Gamma\approx 15$.

FROM GAS TO STARS IN ENERGETIC ENVIRONMENTS: DENSE GAS CLUMPS IN THE 30 DORADUS REGION WITHIN THE LARGE MAGELLANIC CLOUD

We present parsec scale interferometric maps of HCN(1-0) and HCO$^{+}$(1-0) emission from dense gas in the star-forming region 30 Doradus, obtained using the Australia Telescope Compact Array. This extreme star-forming region, located in the Large Magellanic Cloud (LMC), is characterized by a very intense ultraviolet ionizing radiation field and sub-solar metallicity, both of which are expected to impact molecular cloud structure. We detect 13 bright, dense clumps within the 30 Doradus-10 giant molecular cloud. Some of the clumps are aligned along a filamentary structure with a characteristic spacing that is consistent with formation via the varicose fluid instability. Our analysis shows that the filament is gravitationally unstable and collapsing to form stars. There is a good correlation between HCO$^{+}$ emission in the filament and signatures of recent star formation activity including H$_{2}$O masers and young stellar objects (YSOs). YSOs seem to continue along the same direction of the filament toward the massive compact star cluster R136 in the southwest. We present detailed comparisons of clump properities (masses, linewidths, sizes) in 30Dor-10 to those in other star forming regions of the LMC (N159, N113, N105, N44). Our analysis shows that the 30Dor-10 clumps have similar mass but wider linewidths and similar HCN/HCO$^{+}$ (1-0) line ratios as clumps detected in other LMC star-forming regions. Our results suggest that the dense molecular gas clumps in the interior of 30Dor-10 are well-shielded against the intense ionizing field that is present in the 30Doradus region.

Japanese VLBI Network Observations of a Gamma-Ray Narrow-Line Seyfert 1 Galaxy 1H 0323+342

We made simultaneous single-dish and VLBI observations of a gamma-ray narrow-line Seyfert 1 (NLS1) galaxy 1H 0323+342. We found significant flux variation at 8 GHz on a time scale of one month. The total flux density varied by 5.5% in 32 days, corresponding to a variability brightness temperature of 7.0 × 1011 K. We also obtained brightness temperatures of greater than 5.2 × 1010 K from the VLBI images. These high brightness temperatures suggest that the source has nonthermal processes in the central engine. The source structure could be modelled by two elliptical Gaussian components on the parsec scales. The flux of the central component decreases in the same way as the total flux density, showing that the short-term variability is mainly associated with this component.

On detecting oscillations of gamma rays into axion-like particles in turbulent and coherent magnetic fields

Background radiation fields pervade the Universe, and above a certain energy any γ-ray flux emitted by an extragalactic source should be attenuated due to e+e- pair production. The opacity could be alleviated if photons oscillated into hypothetical axion-like particles (ALPs) in ambient magnetic fields, leading to a γ-ray excess especially at high optical depths that could be detected with imaging air Cherenkov telescopes (IACTs). Here, we introduce a method to search for such a signal in γ-ray data and to estimate sensitivities for future observations. Different magnetic fields close to the γ-ray source are taken into account in which photons can convert into ALPs that then propagate unimpeded over cosmological distances until they re-convert in the magnetic field of the Milky Way. Specifically, we consider the coherent field at parsec scales in a blazar jet as well as the turbulent field inside a galaxy cluster. For the latter, we explicitly derive the transversal components of a magnetic field with gaussian turbulence which are responsible for the photon-ALP mixing. To illustrate the method, we apply it to a mock IACT array with characteristics similar to the Cherekov Telescope Array and investigate the dependence of the sensitivity to detect a γ-ray excess on the magnetic-field parameters.

Multifrequency VLBI follow up study of strong γ-ray flares in the blazars 3C273 and 3C279

AbstractWe present multifrequency VLBI observations of the blazars 3C273 and 3C279 after detecting strong γ-ray flares in both of them. 3C273 exhibited a prominent flare in γ-rays in September 2009 which was followed by a strong flare in the 7 mm VLBI core and emergence of a new feature in the parsec scale jet. We have used time delay between flares in different wavebands together with kinematic analysis to determine that the γ-ray emission zone in 3C273 is located 3.6-5.3 pc upstream from the apparent 7 mm core. We have also analyzed frequency dependent core position to measure a deprojected distance between 7 mm core and the true base of the jet: 1-6 pc for 3C273 and 1-3 pc for 3C279, depending on observing epoch. For 3C279 light curve analysis did not give a robust γ-radio delay because there were too many overlapping flares in this source during considered period.

MAGNETIC FIELDS AND MASSIVE STAR FORMATION

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

How to detect supermassive binary black holes at parsec scales

AbstractIt is difficult to find or identifying the binary black holes in parsec scales, since the dual AGN may be merged quickly. It is required to explore more possibilities to identifying binary black holes in parsec scales, we give some discussions, especially with the VLBI methods.

Compact object mergers: observations of supermassive binary black holes and stellar tidal disruption events

AbstractThe capture and disruption of stars by supermassive black holes (SMBHs), and the formation and coalescence of binaries, are inevitable consequences of the presence of SMBHs at the cores of galaxies. Pairs of active galactic nuclei (AGN) and binary SMBHs are important stages in the evolution of galaxy mergers, and an intense search for these systems is currently ongoing. In the early and advanced stages of galaxy merging, observations of the triggering of accretion onto one or both BHs inform us about feedback processes and BH growth. Identification of the compact binary SMBHs at parsec and sub-parsec scales provides us with important constraints on the interaction processes that govern the shrinkage of the binary beyond the “final parsec”. Coalescing binary SMBHs are among the most powerful sources of gravitational waves (GWs) in the universe. Stellar tidal disruption events (TDEs) appear as luminous, transient, accretion flares when part of the stellar material is accreted by the SMBH. About 30 events have been identified by multi-wavelength observations by now, and they will be detected in the thousands in future ground-based or space-based transient surveys. The study of TDEs provides us with a variety of new astrophysical tools and applications, related to fundamental physics or astrophysics. Here, we provide a review of the current status of observations of SMBH pairs and binaries, and TDEs, and discuss astrophysical implications.

COMPARISON OF PRESTELLAR CORE ELONGATIONS AND LARGE-SCALE MOLECULAR CLOUD STRUCTURES IN THE LUPUS I REGION

Turbulence and magnetic fields are expected to be important for regulating molecular cloud formation and evolution. However, their effects on subparsec to 100 parsec scales, leading to the formation of starless cores, is not well understood. We investigate the prestellar core structure morphologies obtained from analysis of the Herschel-SPIRE 350 $\mu$m maps of the Lupus I cloud. This distribution is first compared on a statistical basis to the large scale shape of the main filament. We find the distribution of the elongation position angle of the cores to be consistent with a random distribution, which means no specific orientation of the morphology of the cores is observed with respect to a large-scale filament shape model for Lupus I, or relative to a large-scale bent filament model. This distribution is also compared to the mean orientation of the large-scale magnetic fields probed at 350 $\mu$m with the Balloon-borne Large Aperture Telescope for Polarimetry (BLASTPol) during its 2010 campaign. Here again we do not find any correlation between the core morphology distribution and the average orientation of the magnetic fields on parsec scales. Our main conclusion is that the local filament dynamics - including secondary filaments that often run orthogonally to the primary filament - and possibly small-scale variations in the local magnetic field direction, could be the dominant factors for explaining the final orientation of each core.

On the spine-layer scenario for the very high-energy emission of NGC 1275

We discuss the $\gamma$-ray emission of the radiogalaxy NGC 1275 (the central galaxy of the Perseus Cluster), detected by Fermi-LAT and MAGIC, in the framework of the "spine-layer" scenario, in which the jet is assumed to be characterized by a velocity structure, with a fast spine surrounded by a slower layer. The existence of such a structure in the parsec scale jet of NGC 1275 has been recently proved through VLBI observations. We discuss the constraints that the observed spectral energy distribution imposes to the parameters and we present three alternative models, corresponding to three different choices of the angles between the jet and the line of sight ($\theta_{\rm v}=6^{\circ}, 18^{\circ}$ and 25$^{\circ}$). While for the the case with $\theta_{\rm v}=6^{\circ}$ we obtain an excellent fit, we consider this solution unlikely, since such small angles seems to be excluded by radio observations of the large-scale jet. For $\theta_{\rm v}=25^{\circ}$ the required large intrinsic luminosity of the soft (IR--optical) component of the spine determines a large optical depth for $\gamma$-rays through the pair production scattering $\gamma \gamma\rightarrow e^+ e^-$, implying a narrow cut--off at $\sim50$~GeV. We conclude that intermediate angles are required. In this case the low frequency and the high--energy emissions are produced by two separate regions and, in principle, a full variety of correlations is expected. The correlation observed between the optical and the $\gamma$-ray flux, close to linearity, is likely linked to variations of the emissivity of the spine.

Chemo-Kinematic Survey of z ∼ 1 Star Forming Galaxies using Keck OSIRIS LGS-AO

AbstractWe present first results from the Intermediate Redshift OSIRIS Chemo-Kinematic Survey (IROCKS) of z ∼ 1 star forming galaxies (Mieda et al. in prep). We have targeted Hα and [NII] emission lines in J-band and have spatially resolved the galaxies at sub-kilo parsec scale. We have combined our sample with deep HST continuum images, and are able to reveal the dynamics, morphologies, metallicity distribution, emission-line diagnostics, and star formation rates of galaxies spanning this crucial z ∼ 1 epoch.

The structure function of Galactic H i opacity fluctuations on au scales based on MERLIN, VLA and VLBA data

We use MERLIN, VLA and VLBA observations of Galactic \HI absorption towards 3C~138 to estimate the structure function of the \HI opacity fluctuations at AU scales. Using Monte Carlo simulations, we show that there is likely to be a significant bias in the estimated structure function at signal-to-noise ratios characteristic of our observations, if the structure function is constructed in the manner most commonly used in the literature. We develop a new estimator that is free from this bias and use it to estimate the true underlying structure function slope on length scales ranging $5$ to $40$~AU. From a power law fit to the structure function, we derive a slope of $0.81^{+0.14}_{-0.13}$, i.e. similar to the value observed at parsec scales. The estimated upper limit for the amplitude of the structure function is also consistent with the measurements carried out at parsec scales. Our measurements are hence consistent with the \HI opacity fluctuation in the Galaxy being characterized by a power law structure function over length scales that span six orders of magnitude. This result implies that the dissipation scale has to be smaller than a few AU if the fluctuations are produced by turbulence. This inferred smaller dissipation scale implies that the dissipation occurs either in (i) regions with densities $\gtrsim 10^3 $cm$^-3$ (i.e. similar to that inferred for "tiny scale" atomic clouds or (ii) regions with a mix of ionized and atomic gas (i.e. the observed structure in the atomic gas has a magneto-hydrodynamic origin).

The MOG weak field approximation – II. Observational test of Chandra X-ray clusters

We apply the weak field approximation limit of the covariant Scalar-Tensor-Vector Gravity (STVG) theory, so-called MOdified gravity (MOG), to the dynamics of clusters of galaxies by using only baryonic matter. The MOG effective gravitational potential in the weak field approximation is composed of an attractive Newtonian term and a repulsive Yukawa term with two parameters $\alpha$ and $\mu$. The numerical values of these parameters have been obtained by fitting the predicted rotation curves of galaxies to observational data, yielding the best fit result: $\alpha = 8.89 \pm 0.34$ and $\mu= 0.042\pm 0.004$ kpc$^{-1}$~\cite{rah13}. We extend the observational test of this theory to clusters of galaxies, using data for the ionized gas and the temperature profile of nearby clusters obtained by the Chandra X-ray telescope. Using the MOG virial theorem for clusters, we compare the mass profiles of clusters from observation and theory for eleven clusters. The theoretical mass profiles for the inner parts of clusters exceed the observational data. However, the observational data for the inner parts of clusters (i.e., $r<0.1 r_{500}$) is scattered, but at distances larger than $\sim 300$ kpc, the observed and predicted mass profiles converge. Our results indicate that MOG as a theory of modified gravity is compatible with the observational data from the the solar system to Mega parsec scales without invoking dark matter.

HIERARCHICAL STAR FORMATION IN NEARBY LEGUS GALAXIES

Hierarchical structure in ultraviolet images of 12 late-type LEGUS galaxies is studied by determining the numbers and fluxes of nested regions as a function of size from ~1 to ~200 pc, and the number as a function of flux. Two starburst dwarfs, NGC 1705 and NGC 5253, have steeper number-size and flux-size distributions than the others, indicating high fractions of the projected areas filled with star formation. Nine subregions in 7 galaxies have similarly steep number-size slopes, even when the whole galaxies have shallower slopes. The results suggest that hierarchically structured star-forming regions several hundred parsecs or larger represent common unit structures. Small galaxies dominated by only a few of these units tend to be starbursts. The self-similarity of young stellar structures down to parsec scales suggests that star clusters form in the densest parts of a turbulent medium that also forms loose stellar groupings on larger scales. The presence of super star clusters in two of our starburst dwarfs would follow from the observed structure if cloud and stellar subregions more readily coalesce when self-gravity in the unit cell contributes more to the total gravitational potential.

On the deceleration of Fanaroff–Riley Class I jets: mass loading by stellar winds

Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.

The dusty torus in the Circinus galaxy: a dense disk and the torus funnel

(Abridged) With infrared interferometry it is possible to resolve the nuclear dust distributions that are commonly associated with the dusty torus in active galactic nuclei (AGN). The Circinus galaxy hosts the closest Seyfert 2 nucleus and previous interferometric observations have shown that its nuclear dust emission is well resolved. To better constrain the dust morphology in this active nucleus, extensive new observations were carried out with MIDI at the Very Large Telescope Interferometer. The emission is distributed in two distinct components: a disk-like emission component with a size of ~ 0.2 $\times$ 1.1 pc and an extended component with a size of ~ 0.8 $\times$ 1.9 pc. The disk-like component is elongated along PA ~ 46{\deg} and oriented perpendicular to the ionisation cone and outflow. The extended component is elongated along PA ~ 107{\deg}, roughly perpendicular to the disk component and thus in polar direction. It is interpreted as emission from the inner funnel of an extended dust distribution and shows a strong increase in the extinction towards the south-east. We find no evidence of an increase in the temperature of the dust towards the centre. From this we infer that most of the near-infrared emission probably comes from parsec scales as well. We further argue that the disk component alone is not sufficient to provide the necessary obscuration and collimation of the ionising radiation and outflow. The material responsible for this must instead be located on scales of ~ 1 pc, surrounding the disk. The clear separation of the dust emission into a disk-like emitter and a polar elongated source will require an adaptation of our current understanding of the dust emission in AGN. The lack of any evidence of an increase in the dust temperature towards the centre poses a challenge for the picture of a centrally heated dust distribution.