The Central Molecular Zone (CMZ) of the Galactic Center has to date only been fully mapped at mm wavelengths with singledish telescopes, with resolution about 30$^{\prime\prime}$ (1.2 pc). Using CARMA, we mapped the innermost 0.25 square degrees of the CMZ over the region between -0$.\kern-.25em ^{^\circ}$2$ \leq l \leq $0$.\kern-.25em ^{^\circ}$5 and -0$.\kern-.25em ^{^\circ}$2$ \leq b \leq $0$.\kern-.25em ^{^\circ}$2 (90$\times$50 pc) with spatial and spectral resolution of $\sim$10$^{\prime\prime}$ (0.4 pc) and $\sim$2.5 km/s, respectively. We provide a catalog of 3mm continuum sources as well as spectral line images of SiO(J=2-1), HCO$^{+}$(J=1-0), HCN(J=1-0), N$_{2}$H$^{+}$(J=1-0), and CS(J=2-1), with velocity coverage VLSR= -200 to 200 km/s. To recover the large scale structure resolved out by the interferometer, the continuum-subtracted spectral line images were combined with data from the Mopra 22-m telescope survey, thus providing maps containing all spatial frequencies down to the resolution limit. We find that integrated intensity ratio of I(HCN)/I(HCO$^{+}$) is anti-correlated with the intensity of the 6.4 keV Fe K$\alpha$, which is excited either by high energy photons or low energy cosmic rays, and the gas velocity dispersion as traced by HCO$^{+}$ is correlated with Fe K$\alpha$ intensity. The intensity ratio and velocity dispersion patterns are consistent with variation expected from the interaction of low energy cosmic rays with molecular gas.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

We present radio and infrared observations indicating on-going star formation activity inside the $\sim2-5$ pc circumnuclear ring at the Galactic center. Collectively these measurements suggest a continued disk-based mode of on-going star formation has taken place near Sgr A* over the last few million years. First, VLA observations with spatial resolution 2.17$"\times0.81"$ reveal 13 water masers, several of which have multiple velocity components. The presence of interstellar water masers suggests gas densities that are sufficient for self-gravity to overcome the tidal shear of the 4$\times10^6$ \msol\, black hole. Second, SED modeling of stellar sources indicate massive YSO candidates interior to the molecular ring, supporting in-situ star formation near Sgr A* and appear to show a distribution similar to that of the counter-rotating disks of $\sim$100 OB stars orbiting Sgr A*. Some YSO candidates (e.g., IRS~5) have bow shock structures suggesting that they have have gaseous disks that are phototoevaporated and photoionized by the strong radiation field. Third, we detect clumps of SiO (2-1) and (5-4) line emission in the ring based on CARMA and SMA observations. The FWHM and luminosity of the SiO emission is consistent with shocked protostellar outflows. Fourth, two linear ionized features with an extent of $\sim0.8$ pc show blue and redshifted velocities between $+50$ and $-40$ \kms, suggesting protostellar jet driven outflows with mass loss rates of $\sim5\times10^{-5}$ solar mass yr$^{-1}$. Finally, we present the imprint of radio dark clouds at 44 GHz, representing a reservoir of molecular gas that feeds star formation activity close to Sgr A*.

Phase fluctuations introduced by the atmosphere are the main limiting factor in attaining diffraction limited performance in extended interferometric arrays at millimeter and submillimeter wavelengths. We report the results of C-PACS, the Combined Array for Research in Millimeter-Wave Astronomy Paired Antenna Calibration System. We present a systematic study of several hundred test observations taken during the 2009-2010 winter observing season where we utilize CARMA's eight 3.5-m antennas to monitor an atmospheric calibrator while simultaneously acquiring science observations with 6.1-m and 10.4-m antennas on baselines ranging from a few hundred meters to ~2 km. We find that C-PACS is systematically successful at improving coherence on long baselines under a variety of atmospheric conditions. We find that the angular separation between the atmospheric calibrator and target source is the most important consideration, with consistently successful phase correction at CARMA requiring a suitable calibrator located $\lesssim$6$^\circ$ away from the science target. We show that cloud cover does not affect the success of C-PACS. We demonstrate C-PACS in typical use by applying it to the observations of the nearby very luminous infrared galaxy Arp 193 in $^{12}$CO(2-1) at a linear resolution of ~70 pc (0.12" x 0.18"), 3 times better than previously published molecular maps of this galaxy. We resolve the molecular disk rotation kinematics and the molecular gas distribution and measure the gas surface densities and masses on 90 pc scales. We find that molecular gas constitutes $\sim30\%$ of the dynamical mass in the inner 700 pc of this object with a surface density $\sim10^4 M_\odot$ pc$^{-2}$; we compare these properties to those of the starburst region of NGC 253.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

We present 74 MHz radio continuum observations of the Galactic center region. These measurements show nonthermal radio emission arising from molecular clouds that is unaffected by free-free absorption along the line of sight. We focus on one cloud, G0.13--0.13, representative of the population of molecular clouds that are spatially correlated with steep spectrum (alpha^{74MHz}_{327MHz}=1.3\pm0.3) nonthermal emission from the Galactic center region. This cloud lies adjacent to the nonthermal radio filaments of the Arc near l~0.2^0 and is a strong source of 74 MHz continuum, SiO (2-1) and FeI Kalpha 6.4 keV line emission. This three-way correlation provides the most compelling evidence yet that relativistic electrons, here traced by 74 MHz emission, are physically associated with the G0.13--0.13 molecular cloud and that low energy cosmic ray electrons are responsible for the FeI Kalpha line emission. The high cosmic ray ionization rate ~10-13 s-1 H-1 is responsible for heating the molecular gas to high temperatures and allows the disturbed gas to maintain a high velocity dispersion. LVG modeling of multi-transition SiO observations of this cloud implies H2 densities ~104-5 cm-3 and high temperatures. The lower limit to the temperature of G0.13-0.13 is ~100K, whereas the upper limit is as high as 1000K. Lastly, we used a time-dependent chemical model in which cosmic rays drive the chemistry of the gas to investigate for molecular line diagnostics of cosmic ray heating. When the cloud reaches chemical equilibrium, the abundance ratios of HCN/HNC and N2H+/HCO+ are consistent with measured values. In addition, significant abundance of SiO is predicted in the cosmic ray dominated region of the Galactic center. We discuss different possibilities to account for the origin of widespread SiO emission detected from Galactic center molecular clouds.

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

ALMA observations of the Galactic center with spatial resolution $2.61"\times0.97"$ resulted in the detection of 11 SiO (5-4) clumps of molecular gas within 0.6pc (15$"$) of Sgr A*, interior to the 2-pc circumnuclear molecular ring. The three SiO (5-4) clumps closest to Sgr A* show the largest central velocities, $\sim150$ \kms, and broadest asymmetric linewidths with full width zero intensity (FWZI) $\sim110-147$ \kms. The remaining clumps, distributed mainly to the NE of the ionized mini-spiral, have narrow FWZI ($\sim18-56$ \kms). Using CARMA SiO (2-1) data, LVG modeling of the the SiO line ratios for the broad velocity clumps, constrains the column density N(SiO) $\sim10^{14}$ cm$^{-2}$, and the H$_2$ gas density n$_{\rm H_2}=(3-9)\times10^5$ cm$^{-3}$ for an assumed kinetic temperature 100-200K. The SiO clumps are interpreted as highly embedded protostellar outflows, signifying an early stage of massive star formation near Sgr A* in the last $10^4-10^5$ years. Support for this interpretation is provided by the SiO (5-4) line luminosities and velocity widths which lie in the range measured for protostellar outflows in star forming regions in the Galaxy. Furthermore, SED modeling of stellar sources shows two YSO candidates near SiO clumps, supporting in-situ star formation near Sgr A*. We discuss the nature of star formation where the gravitational potential of the black hole dominates. In particular, we suggest that external radiative pressure exerted on self-shielded molecular clouds enhances the gas density, before the gas cloud become gravitationally unstable near Sgr A*. Alternatively, collisions between clumps in the ring may trigger gravitational collapse.

The 229 GHz (lambda 1.3mm) radio emission from Orion-KL was mapped with up to 0.14" angular resolution with CARMA, allowing measurements of the flux densities of Source I ('SrcI') and the Becklin-Neugebauer Object (BN), the 2 most massive stars in this region. We find integrated flux densities of 310 +/- 45 mJy for SrcI and 240 +/- 35 mJy for BN. SrcI is optically thick even at 229 GHz. No trace of the H30alpha recombination line is seen in its spectrum, although the v_2=1, 5(5,0)-6(4,3) transition of H2O, 3450 K above the ground state, is prominent. SrcI is elongated at position angle 140 degrees, as in 43 GHz images. These results are most easily reconciled with models in which the radio emission from SrcI arises via the H- free-free opacity in a T < 4500 K disk, as considered by Reid et al. (2007). By contrast, the radio spectrum of BN is consistent with p+/e- free-free emission from a dense (n_e ~ 5x10^7 cm^{-3}), but otherwise conventional, hypercompact HII region. The source is becoming optically thin at 229 GHz, and the H30alpha recombination line, at VLSR = 23.2 +/- 0.5 km/sec, is prominent in its spectrum. A Lyman continuum flux of 5x10^{45} photons/sec, consistent with that expected from a B star, is required to maintain the ionization. Supplementary 90 GHz observations were made to measure the H41alpha and H42alpha recombination lines toward BN. Published 43 and 86 GHz data suggest that SrcI brightened with respect to BN over the 15 year period from 1994 to 2009.

The compact radio source Sgr A* is coincident with a 4 million solar mass black hole at the dynamical center of the Galaxy and is surrounded by dense orbiting ionized and molecular gas. We present high resolution radio continuum images of the central 3' and report a faint continuous linear structure centered on Sgr A* with a PA~60 degrees. The extension of this feature appears to be terminated symmetrically by two linearly polarized structures at 8.4 GHz, ~75" from Sgr A*. A number of weak blobs of radio emission with X-ray counterparts are detected along the axis of the linear structure. The linear structure is best characterized by a mildly relativistic jet from Sgr A* with an outflow rate 10^-6 solar mass per year. The near and far-sides of the jet are interacting with orbiting ionized and molecular gas over the last 1-3 hundred years and are responsible for a 2" hole, the "minicavity", characterized by disturbed kinematics, enhanced FeII/III line emission, and diffuse X-ray gas. The estimated kinetic luminosity of the outflow is ~1.2x10^{41} erg/s, so the interaction with the bar may be responsible for the Galactic center X-ray flash inferred to be responsible for much of the fluorescent Fe Kalpha line emission from the inner 100pc of the Galaxy.

Building on the "drPACS" contribution at ADASS XX of a simple Unix pipeline infrastructure, we implemented a pipeline toolkit using the package MIRIAD to combine Interferometric and Single Dish data (MIS). This was prompted by our observations made with the Combined Array For Research in Millimeter-wave Astronomy (CARMA) interferometer of the star-forming region NGC 1333, a large survey highlighting the new 23-element and singledish observing modes. The project consists of 20 CARMA datasets each containing interferometric as well as simultaneously obtained single dish data, for 3 molecular spectral lines and continuum, in 527 different pointings, covering an area of about 8 by 11 arcminutes. A small group of collaborators then shared this toolkit and their parameters via CVS, and scripts were developed to ensure uniform data reduction across the group. The pipeline was run end-to-end each night as new observations were obtained, producing maps that contained all the data to date. We will show examples of the scripts and data products. This approach could serve as a model for repeated calibration and mapping of large mixed-mode correlation datasets from ALMA.

Using CARMA, we imaged the 87 GHz SiO v=0 J=2-1 line toward Orion-KL with 0.45 arcsec angular resolution. The maps indicate that radio source I drives a bipolar outflow into the surrounding molecular cloud along a NE--SW axis, in agreement with the model of Greenhill et al. (2004). The extended high velocity outflow from Orion-KL appears to be a continuation of this compact outflow. High velocity gas extends farthest along a NW--SE axis, suggesting that the outflow direction changes on time scales of a few hundred years.

We report the highest-fidelity observations of the spiral galaxy M51 in CO emission, revealing the evolution of giant molecular clouds (GMCs) vis-a-vis the large-scale galactic structure and dynamics. The most massive GMCs (so-called GMAs) are first assembled and then broken up as the gas flow through the spiral arms. The GMAs and their H2 molecules are not fully dissociated into atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large-scale galactic dynamics --their coagulation into GMAs is due to spiral arm streaming motions upon entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream side. In M51, the majority of the gas remains molecular from arm entry through the inter-arm region and into the next spiral arm passage.

We report on Combined Array for Research in Millimeter-Wave Astronomy (CARMA) and James Clerk Maxwell Telescope (JCMT) observations toward the Einstein source 1E 1740.7-2942, a LMXB commonly known as the "Great Annihilator." The Great Annihilator is known to be near a small, bright molecular cloud on the sky in a region largely devoid of emission in 12-CO surveys of the Galactic Center. The region is of interest because it is interior to the dust lanes which may be the shock zones where atomic gas from HI nuclear disk is converted into molecular gas. We find that the region is populated with a number of dense (n ~ 10^5 cm^-3) regions of excited gas with small filling factors, and estimate that up to 1-3 x 10^5 solar masses of gas can be seen in our maps. The detection suggests that a significant amount of mass is transported from the shock zones to the GC star-forming regions in the form of small, dense bundles.

We present resolved images of the dust continuum emission from the debris disk around the young (80-200 Myr) solar-type star HD 107146 with CARMA at $\lambda$1.3 mm and the CSO at $\lambda$350 $\mu$m. Both images show that the dust emission extends over an $\sim$10\arcsec diameter region. The high resolution (3\arcsec) CARMA image further reveals that the dust is distributed in a partial ring with significant decrease in flux inward of 97 AU. Two prominent emission peaks appear within the ring separated by $\sim$140 degrees in position angle. The morphology of the dust emission is suggestive of dust captured into a mean motion resonance, which would imply the presence of a planet at an orbital radius of $\sim$45-75 AU.

We report the detection of CO(2-1) emission from the z=1.2 ultraluminous infrared galaxy (ULIRG) GOODS J123634.53+621241.3 (also known as the sub-millimeter galaxy GN26). These observations represent the first discovery of high-redshift CO emission using the new Combined Array for Research in Millimeter-Wave Astronomy (CARMA). Of all high-redshift (z>1) galaxies within the GOODS-North field, this source has the largest far-infrared (FIR) flux observed in the Spitzer 70um and 160um bands. The CO redshift confirms the optical identification of the source, and the bright CO(2-1) line suggests the presence of a large molecular gas reservoir of about 7x10^10 M(sun). The infrared-to-CO luminosity ratio of L(IR)/L'(CO) = 80+/-30 L(sun) (K Km/s pc^2)^-1 is slightly smaller than the average ratio found in local ULIRGs and high-redshift sub-millimeter galaxies. The short star-formation time scale of about 70 Myr is consistent with a starburst associated with the merger event and is much shorter than the time scales for spiral galaxies and estimates made for high-redshift galaxies selected on the basis of their B-z and z-K colors.

We investigate numerically the hydrodynamic instability of an ionization front (IF) accelerating into a molecular cloud, with imposed initial perturbations of different amplitudes. When the initial amplitude is small, the imposed perturbation is completely stabilized and does not grow. When the initial perturbation amplitude is large enough, roughly the ratio of the initial amplitude to wavelength is greater than 0.02, portions of the IF temporarily separate from the molecular cloud surface, locally decreasing the ablation pressure. This causes the appearance of a large, warm HI region and triggers nonlinear dynamics of the IF. The local difference of the ablation pressure and acceleration enhances the appearance and growth of a multimode perturbation. The stabilization usually seen at the IF in the linear regimes does not work due to the mismatch of the modes of the perturbations at the cloud surface and in density in HII region above the cloud surface. Molecular pillars are observed in the late stages of the large amplitude perturbation case. The velocity gradient in the pillars is in reasonably good agreement with that observed in the Eagle Nebula. The initial perturbation is imposed in three different ways: in density, in incident photon number flux, and in the surface shape. All cases show both stabilization for a small initial perturbation and large growth of the second harmonic by increasing amplitude of the initial perturbation above a critical value.

Recent Chandra and XMM-{\it Newton} observations reported evidence of two X-ray filaments G359.88-0.08 (SgrA-E) and G359.54+0.18 (the ripple filament) near the Galactic center. The X-ray emission from these filaments has a nonthermal spectrum and coincides with synchrotron emitting radio sources. Here, we report the detection of a new X-ray feature coincident with a radio filament G359.90-0.06 (SgrA-F) and show more detailed VLA, Chandra and BIMA observations of the radio and X-ray filaments. In particular, we show that radio emission from the nonthermal filaments G359.90-0.06 (SgrA-F) and G359.54+0.18 (the ripple) has a steep spectrum whereas G359.88-0.08 (SgrA-E) has a flat spectrum. The X-ray emission from both these sources could be due to synchrotron radiation. However, given that the 20 \kms molecular cloud, with its intense 1.2mm dust emission, lies in the vicinity of SgrA-F, it is possible that the X-rays could be produced by inverse Compton scattering of far-infrared photons from dust by the relativistic electrons responsible for the radio synchrotron emission. The production of X-ray emission from ICS allows an estimate of the magnetic field strength of ~0.08 mG within the nonthermal filament. This should be an important parameter for any models of the Galactic center nonthermal filaments.

We investigate hydrodynamic instability of accelerating ionization fronts using two dimensional hydrodynamic simulations that include detailed energy deposition and release due to the absorption of UV radiation, recombination of hydrogen, radiative molecular cooling, and magnetic pressure. We consider linear perturbation growth and find that the stabilization mechanism associated with non-accelerated fronts remains a significant factor even when acceleration is present. Conversely, if recombination in the ionized region is turned off,Rayleigh-Taylor (RT) instability becomes effective, and the classical RT growth rate recovered.