Search

In the Milky Way and other main-sequence galaxies, stars form exclusively in molecular gas, which is traced by CO emission. However, low metallicity dwarf galaxies are often `CO-dark' in the sense that CO emission is not observable even at the high resolution and sensitivities of modern observing facilities. In this work we use ultra high-resolution simulations of four low-metalicity dwarf galaxies (which resolve star formation down to the scale of star-forming cores, 0.01 pc) combined with a time-dependent treatment of the chemistry of the interstellar medium, to investigate the star formation environment in this previously hidden regime. By generating synthetic observations of our models we show that the galaxies have high to extremely high dark gas fractions (0.13 to 1.00 dependent on beam size and conditions), yet despite this form stars. However, when examined on smaller scales, we find that the stars still form in regions dominated by molecular gas, it is simply that these are far smaller than the scale of the beam (1.5"). Thus, while stars in CO-dark dwarf galaxies form in small molecular cores like larger galaxies, their cloud-scale environment is very different., Comment: Re-submitted to MNRAS after referee report, comments and questions welcome, 23 pages, 16 figures

Achieving low-threshold infrared stimulated emission in solution-processed quantum dots is critical to enable real-life application including photonic integrated circuits (PICs), LIDAR application and optical telecommunication. However, realization of low threshold infrared gain is fundamentally challenging due to high degeneracy of the first emissive state (e.g., 8-fold) and fast Auger recombination. In this letter, we demonstrate ultralow-threshold infrared stimulated emission with an onset of 110 uJ.cm-2 employing cascade charge transfer (CT) in Pb-chalcogenide colloidal quantum dot (CQD) solids. In doing so, we investigate this idea in two different architectures including a mixture of multiband gap CQDs and layer-by-layer (LBL) configuration. Using transient absorption spectroscopy, we show ultrafast cascade CT from large band-gap PbS CQD to small band-gap PbS/PbSSe core/shell CQDs in LBL (~ 2 ps) and mixture (~ 9 ps) configuration. These results indicate the feasibility of using cascade CT as an efficient method to reduce optical gain threshold in CQD solid films.

We have investigated the possibility of molecular cloud formation via the Collision-induced Magnetic Reconnection (CMR) mechanism of the cold neutral medium (CNM). Two atomic gas clouds with conditions typical of the CNM were set to collide at the interface of reverse magnetic fields. The cloud-cloud collision triggered magnetic reconnection and produced a giant 20pc filamentary structure which was not seen in the control models without CMR. The cloud, with rich fiber-like sub-structures, developed a fully molecular spine at 5Myr. Radiative transfer modeling of dust emission at far infrared wavelengths showed that the middle part of the filament contained dense cores over a span of 5pc. Some of the cores were actively forming stars and typically exhibited both connecting fibers in dust emission and high-velocity gas in CO line emission, indicative of active accretion through streamers. Supersonic turbulence was present in and around the CMR-filament due to inflowing gas moving at supersonic velocities in the collision mid-plane. The shocked gas was condensed and transported to the main filament piece by piece by reconnected fields, making the filament and star formation a bottom-up process. Instead of forming a gravitationally bounded cloud which then fragments hierarchically (top-down) and forms stars, the CMR process creates dense gas pieces and magnetically transports them to the central axis to constitute the filament. Since no turbulence is manually driven, our results suggest that CMR is capable of self-generating turbulence. Finally, the resulting helical field should show field-reversal on both sides of the filament from most viewing angles., Comment: accepted by ApJ

The relation between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance to many areas of astrophysics, from protostellar disks to galaxy evolution. We present and compare Bayesian analyses of the B - n relation for a comprehensive observational data set, as well as a large body of numerical MHD simulations. We extend the original Zeeman relation of Crutcher et al. (2010) with a large body of magnetic data that includes 700 observations with the Davis-Chandrasekhar-Fermi method. By using a new multiparameter Bayesian analysis we present a new, more general, time-averaged observational relation: B \propto n^{0.27 \pm 0.017} for n \leq n_0 and B \propto n^{0.54 \pm 0.18} for n \geq n_0, with n_0 = 924^(+145-144) cm^-3. We perform a separate analysis on 19 numerical magnetohydrodynamics simulations that cover a wide range of scales, resolutions, initial conditions, and completed with a variety of codes: arepo, flash, pencil, and ramses. The power law exponents derived from the simulations depend on several physical factors including: dynamo effects, time scales, turbulence, and the initial seed field strength. In particular, early-time simulations where the density, velocity and magnetic fields are unevolved do not match the observational scalings. The simulations that trace the observed density range best, the evolved dwarf galaxy and Milky Way like galaxy simulations, settle into a near consistent exponent of = 0.5 in the dense gas, with variability in the diffuse gas exponent., Comment: 23 figures, 28 pages, submitted to MNRAS, Comments welcome

A crucial aspect in addressing the challenge of measuring the core mass function, that is pivotal for comprehending the origin of the initial mass function, lies in constraining the temperatures of the cores. We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. High angular resolution mapping is required to capture the properties of protostellar and pre-stellar cores and to effectively separate them from larger features, such as dusty filaments. We employed the point process mapping (PPMAP) technique, enabling us to perform spectral energy distribution fitting of far-infrared and submillimeter observations across the 15 ALMA-IMF fields, at an unmatched 2.5" angular resolution. By combining the modified blackbody model with near-infrared data, we derived bolometric luminosity maps. We estimated the errors impacting values of each pixel in the temperature, column density, and luminosity maps. Subsequently, we employed the extraction algorithm getsf on the luminosity maps in order to detect luminosity peaks and measure their associated masses. We obtained high-resolution constraints on the luminosity, dust temperature, and mass of protoclusters, that are in agreement with previously reported measurements made at a coarser angular resolution. We find that the luminosity-to-mass ratio correlates with the evolutionary stage of the studied regions, albeit with intra-region variability. We compiled a PPMAP source catalog of 313 luminosity peaks using getsf on the derived bolometric luminosity maps. The PPMAP source catalog provides constraints on the mass and luminosity of protostars and cores, although one source may encompass several objects. Finally, we compare the estimated luminosity-to-mass ratio of PPMAP sources with evolutionary tracks and discuss the limitations imposed by the 2.5" beam., Comment: 35 pages, 19 figures, 5 tables Accepted by A&A

We use the H$41\alpha$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41\alpha$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41\alpha$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $\Delta V_\mathrm{dyn}$ with $D$ and $Q_0$. $\Delta V_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind., Comment: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $\sigma_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $\sigma_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $\mu$G. Finally, we discuss the implications of our results within the observational context in the CMZ.

High-latitude ($|b|>30^{\circ}$) molecular clouds have virial parameters that exceed 1, but whether these clouds can form stars has not been studied systematically. Using JCMT SCUBA-2 archival data, we surveyed 70 fields that target high-latitude Planck galactic cold clumps (HLPCs) to find dense cores with density of $10^{5}$-$10^{6}$ cm$^{-3}$ and size of $<0.1$ pc. The sample benefits from both the representativeness of the parent sample and covering densest clumps at the high column density end ($>1\times10^{21}$ cm$^{-2}$). At an average noise rms of 15 mJy/beam, we detected Galactic dense cores in only one field, G6.04+36.77 (L183), while also identifying 12 extragalactic objects and two young stellar objects. Compared to the low-latitude clumps, dense cores are scarce in HLPCs. With synthetic observations, the densities of cores are constrained to be $n_c\lesssim10^5$ cm$^{-3}$, should they exist in HLPCs. Low-latitude clumps, Taurus clumps, and HLPCs form a sequence where a higher virial parameter corresponds to a lower dense core detection rate. If HLPCs were affected by the Local Bubble, the scarcity should favor turbulence-inhibited rather than supernova-driven star formation. Studies of the formation mechanism of the L183 molecular cloud are warranted., Comment: 9 pages for the main text. 4 figures, 1 table. Published in Astrophysical Journal Letter

We report 850 $\mu$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity., Comment: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,000 years, yielding median mass, momentum, and energy rates of $\dot{M}$ = 10$^{-4.4}$ M$_{\odot}$ yr$^{-1}$, $\dot{P}$ = 10$^{-3.2}$ M$_{\odot}$ km/s yr$^{-1}$, and $\dot{E}$ = 1 L$_{\odot}$. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (~3$-$5$\sigma$), and no correlations above 3$\sigma$ with clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass - which has been previously described in the literature - may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15-30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force-bolometric luminosity relationship in the literature, and our data extend this relationship up to L $\geq$ 10$^6$ L$_{\odot}$ and $\dot{P}$ $\geq$ 1 M$_{\odot}$ km/s yr$^{-1}$. Our lack of correlation with clump L/M is inconsistent with models of protocluster formation in which all protostars start forming at the same time., Comment: 46 pages, 14 figures, 10 tables. This publication has an associated Zenodo entry, which can be found here: https://zenodo.org/records/8350595

The use of colloidal quantum dots (CQDs) as a gain medium in infrared laser devices has been underpinned by the need for high pumping intensities, very short gain lifetimes and low gain coefficients.

We present observations of polarized dust emission at 850 $\mu$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $\mu$G in the main starless core and up to $\sim$90$\pm$40 $\mu$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfv\'enic. We also present a new method of data reduction for these denser but fainter objects like starless cores., Comment: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

Filamentary structures are often identified in column density maps of molecular clouds, and appear to be important for both low- and high-mass star formation. Theoretically, these structures are expected to form in regions where the supersonic cloud-scale turbulent velocity field converges. While this model of filament formation successfully reproduces several of their properties derived from column densities, it is unclear whether it can also reproduce their kinematic features. We use a combination of hydrodynamical, chemical and radiative transfer modelling to predict the emission properties of these dynamically-forming filaments in the $^{13}$CO, HCN and N$_2$H$^+$ $J=1-0$ rotational lines. The results are largely in agreement with observations; in particular, line widths are typically subsonic to transonic, even for filaments which have formed from highly supersonic inflows. If the observed filaments are formed dynamically, as our results suggest, no equilibrium analysis is possible, and simulations which presuppose the existence of a filament are likely to produce unrealistic results., Comment: 9 pages, 9 figures. MNRAS accepted

We present BISTRO Survey 850 {\mu}m dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965)., Comment: 14 pages, 5 figures. ApJ accepted

The Cold Neutral Medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star forming gas, yet its distribution is still not fully understood. In this work we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total HI due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii.However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM doesn't extend out radially as far as HI, and the vertical scale height is smaller in the outer galaxy compared to HI with no flaring. The CNM column density scales with total midplane pressure and disappears from the gas phase below values of PT/kB =1000 K/cm3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt-Schmidt law. In the outer galaxy we produce realistic vertical velocity dispersions in the HI purely from galactic dynamics but our models do not predict CNM at the extremely large radii observed in HI absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii., Comment: 13 pages, 20 figures, accepted by MNRAS

The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $\alpha=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $\alpha=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores., Comment: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work has discovered a tight, often close-to-linear relationship between the integrated intensity of molecular lines and the total column density of the cloud material. We combine magnetohydrodynamical simulations, time-dependent chemistry, and radiative transfer to produce synthetic molecular line observations of model clouds. We find similarly tight correlations between line intensity and column density to those observed, although the linear behaviour is only seen in isolated (as opposed to colliding) model clouds. This linear relationship is not due to optically thin emission; all lines investigated have high optical depths, and the increase in integrated intensity with column density is due to higher velocity dispersion along the line of sight. Overall, the idealised models commonly used in the literature appear to be reasonably accurate representations of real molecular clouds., Comment: 10 pages, 10 figures. MNRAS accepted. Data publicly available at http://cloudzoo.astro.cf.ac.uk/Downloads/202211/

We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines., Comment: Accepted for publication in Astrophysical Journal (ApJ)

We present 850 $\mu$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments., Comment: 26 pages, 13 figures, ApJ accepted

Synthetic data generation has recently gained widespread attention as a more reliable alternative to traditional data anonymization. The involved methods are originally developed for image synthesis. Hence, their application to the typically tabular and relational datasets from healthcare, finance and other industries is non-trivial. While substantial research has been devoted to the generation of realistic tabular datasets, the study of synthetic relational databases is still in its infancy. In this paper, we combine the variational autoencoder framework with graph neural networks to generate realistic synthetic relational databases. We then apply the obtained method to two publicly available databases in computational experiments. The results indicate that real databases' structures are accurately preserved in the resulting synthetic datasets, even for large datasets with advanced data types., Comment: 8 pages, 2 figures, 2 tables, Synthetic Data 4 ML workshop of the Thirty-sixth Conference on Neural Information Processing Systems (NeurIPS 2022)

Models of prestellar cores often assume that the cores are isolated from their environment - material outside the core boundary plays no role in the subsequent evolution. This is unlikely to be the case in reality, where cores are located within hierarchically substructured molecular clouds. We investigate the dynamical and chemical evolution of prestellar cores, modelled as Bonnor-Ebert spheres, and show that the density of the ambient medium has a large impact on the resulting chemical properties of the cores. Models embedded in high-density, low-temperature surroundings have greatly enhanced abundances of several molecules, such as CO and CS, compared to models with more diffuse surroundings, corresponding to relatively isolated cores. The predicted intensities and profile shapes of molecular lines are also affected. The density of the ambient medium has a stronger effect on the chemical evolution than whether the cores are initially in or out of equilibrium. This suggests that the impact of environment cannot be neglected when modelling chemistry in prestellar cores; the results of these models are highly sensitive to the assumptions made about the core surroundings., Comment: 7 pages, 5 figures. MNRAS accepted

We develop a simple analytic model for what happens when an O star (or compact cluster of OB stars) forms in a shock compressed layer and carves out an approximately circular hole in the layer, at the waist of a bipolar HII Region (HIIR). The model is characterised by three parameters: the half-thickness of the undisturbed layer, Zlay, the mean number-density of hydrogen molecules in the undisturbed layer, nlay, and the (collective) ionising output of the star(s), NdotLyC. The radius of the circular hole is given by WIF ~ 3.8 pc [Zlay/0.1pc]^{-1/6} [nlay/10^4cm^{-3}]^{-1/3} [NdotLyC/10^{49} s^{-1}]^{1/6} [t/Myr]^{2/3}. Similar power-law expressions are obtained for the rate at which ionised gas is fed into the bipolar lobes; the rate at which molecular gas is swept up into a dense ring by the shock front (SF) that precedes the ionisation front (IF); and the density in this dense ring. We suggest that our model might be a useful zeroth-order representation of many observed HIIRs. From viewing directions close to the midplane of the layer, the HIIR will appear bipolar. From viewing directions approximately normal to the layer it will appear to be a limb-brightened shell but too faint through the centre to be a spherically symmetric bubble. From intermediate viewing angles more complicated morphologies can be expected., Comment: To appear in MNRAS, 10 pages, 3 figures, 1 table