Low-mass Star Clusters Research Articles

Violent Starbursts and Quiescence Induced by Far-ultraviolet Radiation Feedback in Metal-poor Galaxies at High Redshift

JWST observations of galaxies at z ≳ 8 suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one galaxy at subparsec resolution with different radiative feedback models. We find that the suppression of cooling by far-UV (FUV) radiation (i.e., H2 dissociating radiation) from Population II stars is the main physical process triggering the formation of compact and massive star clusters and is responsible for the bursty star formation observed in metal-poor galaxies at z ≳ 10. Indeed, artificially suppressing FUV radiation leads to a less intense continuous mode of star formation distributed into numerous but low-mass open star clusters. Due to the intense FUV field, low-metallicity clouds remain warm (∼104 K) until they reach a relatively high density (≳103 cm−3), before becoming self-shielded and transitioning to a colder (∼100 K), partially molecular phase. As a result, star formation is delayed until the clouds accumulate enough mass to become gravitationally unstable. At this point, the clouds undergo rapid star formation, converting gas into stars with high efficiency. We therefore observe exceptionally bright galaxies (10 times brighter than for continuous star formation) and subsequent quenched “dead” galaxies that did not form stars for tens of Myr.

CO in the Draco nebula: the atomic–molecular transition

ABSTRACT This paper presents maps of the J=2-1 transition of CO towards the Draco nebula intermediate-velocity cloud (IVC). The maps cover 8500 square arcmin with a velocity resolution of 0.33 km s$^{-1}$ and angular resolution of 38 ${\rm arcsec}$, or 0.11 pc at the cloud distance of 600 pc. The mapped area includes all the emission detected by the Herschel satellite with 250 $\mu$m intensity >5 MJy sr−1. Previously published observations of the far-infrared (far-IR) emission and the 21 cm line of HI are used to derive the column density distribution of H$_2$ and the abundance ratio CO/H$_2$, as well as the distribution of the molecular fraction of hydrogen, which approaches 90 per cent over much of the brighter parts of the nebula. The CO emission is highly clumpy and closely resembles the structures seen in far-IR images. The kinematics of the CO show supersonic motions between clumps but near-thermal to trans-sonic motions within clumps, consistent with model predictions that the scale length for dissipation of supersonic turbulence should be $\sim 0.1$ pc, mediated by kinematic viscosity and/or ambipolar diffusion. Different parts of the nebula show evidence for a spread of molecular formation time-scales of a few 10$^5$ yr, comparable with the dynamical time-scale of the infalling gas. The IVC will likely merge with the Galactic interstellar medium in $\sim 10^7$ yr, and the densest clumps may form an unbound cluster of low-mass stars.

Compact Binary Formation in Open Star Clusters II: Difficulty of Gaia NS formation in low-mass star clusters

Gaia mission offers opportunities to search for compact binaries not involved in binary interactions (hereafter inert compact binaries), and results in the discoveries of binaries containing one black hole (BH) or one neutron star (NS), called “Gaia BHs’’ and”Gaia NSs’’, respectively. We have assessed if Gaia BHs and NSs can be formed in open clusters through dynamical interactions. In order to obtain a large number of inert compact binaries similar to Gaia BHs and NSs, we have performed gravitational N-body simulations for a large number of open clusters whose total mass is 1.2×108M⊙. These clusters have various masses, metallicities, densities, and binary fractions. We have found that open clusters form Gaia BHs ( 10−6- 10−5M⊙−1) much more efficiently than Gaia NSs ( ≲10−7M⊙−1) for any cluster parameters. This is quite inconsistent with observational results, because the reported numbers of Gaia BHs and NSs are 3 and 21, respectively. Additionally, we have switched off NS natal kicks for 104 open clusters each weighing 103M⊙ in order to retain a large number of NSs in open clusters. Then, open clusters form inert NS binaries originating from primordial binaries rather than formed through dynamical interactions. This means that Gaia NSs are formed dominantly on isolated fields, not in open clusters, if there is no NS natal kick. We have concluded that Gaia BHs can be dominantly formed in open clusters, however Gaia NSs cannot.

Star-forming environments in smoothed particle magnetohydrodynamics simulations II: re-simulating isolated clumps to determine equivalence of extracted clumps and parent simulations

ABSTRACT What is the numerical reproducibility of a stellar system (including its discs) when evolving only a subset of (partially-evolved) smoothed particle hydrodynamics (SPH) particles? To investigate this, we modelled the evolution of 29 star-forming clumps that were extracted from our previous simulations that investigated the formation and early evolution of low-mass star clusters. These clumps were evolved using a three-dimensional smoothed particle radiation magnetohydrodynamics code, where we included or excluded non-ideal magnetohydrodynamics to match the cluster simulation. While star formation proceeded as expected, we were unable to identically reproduce any of the systems present at the end of the cluster simulations. However, the final distributions of stellar mass, stellar system mass, disc mass, and disc radii were reproduced statistically; unfortunately, the distribution of average magnetic field strengths in the discs was not reproduced statistically, but this may be a result of our updated algorithms governing the evolution of the magnetic field. Therefore, given that our clumps yield stellar masses that are statistically similar to those in the original low-mass star clusters, we have demonstrated that we can statistically reproduce systems (aside from their magnetic field strength) by evolving a subset of SPH particles. Therefore, clumps such as these can be used as initial conditions to investigate the formation of isolated stars from less-contrived initial environments.

PECULIAR SPECTRUM OF THE WATER MASER IN THE DARK NEBULA MSXDCG24.33+011(G24.33+014)

Giant molecular clouds (GMC) in our and other galaxies and small dense molecular clouds inside the Galaxy (IRDC) form cores due to gravitational instability, in which massive stars and clusters of low-mass stars arise. The high background of infrared radiation inside the Galaxy creates advantages in favor of IRDCs in the study of star formation processes and accompanying phenomena such as accretion, the appearance of HII regions, bipolar outflows, and others that cause various responses in their molecular composition. As a part of studying the evolutionary state in the cloud IRDC MSXDCG24.33+011 (the alternative name is G24.33+014), observations of the water vapor maser were made. On November 28 (2022) during the observations using the RT-22 of the Pushchino Radio Astronomy Observatory, the H2O maser detail at the velocity of \({{V}_{{{\text{LSR}}}}} = 103.15\) km/s with the linewidth of 0.52 km/s was detected. Peak flux of 49.5(\( \pm 6\)) Jy was recorded. This detail was not detected on RT-22 in the PRAO on July 5 (2022) and has not been seen before by other researchers.

The growth of intermediate mass black holes through tidal captures and tidal disruption events

ABSTRACTWe present N-body simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 M⊙. The models are evolved with the N-body code bifrost to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruption events (TDEs). We model star–BH tidal interactions using a velocity-dependent drag force, which causes orbital energy and angular momentum loss near the BH. About ∼20–30 per cent of the stars within the spheres of influence of the black holes form Bahcall–Wolf cusps and prevent the systems from core collapse. Within the first 40 Myr of evolution, the systems experience 500–1300 TDEs, depending on the initial cluster structure. Most (>95 per cent) of the TDEs originate from stars in the Bahcall–Wolf cusp. We derive an analytical formula for the TDE rate as a function of the central BH mass, density, and velocity dispersion of the clusters ($\dot{N}_{\mathrm{TDE}} \propto M\mathrm{_{BH}}\rho \sigma ^{-3}$). We find that TDEs can lead a 300 M⊙ BH to reach $\sim 7000 \, \mathrm{{M}_{\odot }}$ within a Gyr. This indicates that TDEs can drive the formation and growth of massive BHs in sufficiently dense environments, which might be present in the central regions of nuclear star clusters.

Compact object mergers in hierarchical triples from low-mass young star clusters

Abstract A binary star orbited by an outer companion constitutes a hierarchical triple system. The outer body may excite the eccentricity of the inner binary through the von Zeipel-Lidov-Kozai (ZLK) mechanism, triggering the gravitational wave (GW) coalescence of the inner binary when its members are compact objects. Here, we study a sample of hierarchical triples with an inner black hole (BH) – BH binary, BH – neutron star (NS) binary, and BH – white dwarf (WD) binary, formed via dynamical interactions in low-mass young star clusters. Our sample of triples was obtained self-consistently from direct N-body simulations of star clusters which included up-to-date stellar evolution. We find that the inner binaries in our triples cannot merge via GW radiation alone, and the ZLK mechanism is essential to trigger their coalescence. Contrary to binaries assembled dynamically in young star clusters, binary BHs merging in triples have preferentially low mass ratios (q ≃ 0.3) and higher primary masses ($m_{\rm p} \gtrsim 40 {\, \rm M_\odot }$). We derive a local merger rate density of 0.60, 0.11 and $0.5 {-1} \, \rm Gpc^{-3}$ for BH-BH, BH-NS and BH-WD binaries, respectively. Additionally, we find that merging binaries have high eccentricities across the GW spectrum, including the LIGO-Virgo-KAGRA (LVK), LISA, and DECIGO frequencies. About 7 per cent of BH-BH and 60 per cent of BH-NS binaries will have detectable eccentricities in the LVK band. Our results indicate that the eccentricity and the mass spectrum of merging binaries are the strongest features for the identification of GW mergers from triples.

Synthetic photometry of OB star clusters with stochastically sampled IMFs: analysis of models and HST observations

ABSTRACT We present a pilot library of synthetic NUV, U, B, V, and I photometry of star clusters with stochastically sampled IMFs and ionized gas for initial masses, Mi = 103, 104, and 105 M⊙; t = 1, 3, 4, and 8 Myr; Z = 0.014 and Z = 0.002; and log(US) = −2 and −3. We compare the library with predictions from deterministic models and observations of isolated low-mass (<104 M⊙) star clusters with co-spatial compact H ii regions. The clusters are located in NGC 7793, one of the nearest galaxies observed as part of the HST LEGUS and Hα-LEGUS surveys. (1) For model magnitudes that only account for the stars: (a) the residual |deterministic mag - median stochastic mag| can be ≥0.5 mag, even for Mi = 105 M⊙; and (b) the largest spread in stochastic magnitudes occurs when Wolf–Rayet stars are present. (2) For Mi = 105 M⊙: (a) the median stochastic mag with gas can be >1.0 mag more luminous than the median stochastic magnitude without gas; and (b) nebular emission lines can contribute with $\gt 50{{\ \rm per\ cent}}$ and $\gt 30{{\ \rm per\ cent}}$ to the total emission in the V and I bands, respectively. (3) Age-dating OB-star clusters via deterministic tracks in the U-B versus V-I plane is highly uncertain at Z = 0.014 for Mi ∼ 103 M⊙ and Z = 0.002 for Mi ∼ 103–105 M⊙. (4) For low-mass clusters, the V-band extinction derived with stochastic models significantly depends on the value of log(US). (5) The youngest clusters tend to have higher extinction. (6) The majority of clusters have multi-peaked age PDFs. (7) Finally, we discuss the importance of characterizing the true variance in the number of stars per mass bin in nature.

Dynamics of binary black holes in low-mass young star clusters

ABSTRACT Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters (∼300–103 M⊙) are more numerous than massive systems and are characterized by a two-body relaxation time-scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact objects. Here, we explore the properties of black holes (BHs) and binary black holes (BBHs) formed in low-mass young star clusters, by means of a suite of 105 direct N-body simulations with a high original binary fraction (100 per cent for stars with mass >5 M⊙). Most BHs are ejected in the first ∼20 Myr by dynamical interactions. Dynamical exchanges are the main formation channel of BBHs, accounting for ∼40–80 per cent of all the systems. Most BBH mergers in low-mass young star clusters involve primary BHs with mass <40 M⊙ and low-mass ratios are extremely more common than in the field. Comparing our data with those of more massive star clusters (103 − 3 × 104 M⊙), we find a strong dependence of the percentage of exchanged BBHs on the mass of the host star cluster. In contrast, our results show just a mild correlation between the mass of the host star cluster and the efficiency of BBH mergers.

The Binary INformation from Open Clusters Using SEDs (BINOCS) Project: Reliable Photometric Mass Determinations of Binary Star Systems in Clusters

Abstract We introduce a new binary detection technique, Binary INformation from Open Clusters using SEDs (binocs), which we show is able to determine reliable stellar multiplicity and masses over a much larger mass range than current approaches. This new technique determines accurate component masses of binary and single systems of the open clusters’ main sequence by comparing observed magnitudes from multiple photometric filters to synthetic star spectral energy distributions (SEDs), allowing us to systematically probe the binary population for low-mass stars in clusters for eight well-studied open clusters. We provide new deep, infrared photometric catalogs (1.2–8.0 μm) for the key open clusters NGC 1960 (M36), NGC 2099 (M37), NGC 2420, and NGC 2682 (M67), using observations from NOAO/NEWFIRM and Spitzer/IRAC. Using these deep multiwavelength catalogs, the binocs method is applied to these clusters to determine accurate component masses for unresolved cluster binaries. We explore binary fractions as a function of cluster age, Galactic location, and metallicity.

The possible role of stellar mergers for the formation of multiple stellar populations in globular clusters

ABSTRACTMany possible scenarios for the formation of multiple stellar populations (MSPs) in globular clusters (GCs) have been discussed so far, including the involvement of asymptotic giant branch stars, fast-rotating main-sequence stars, very massive main-sequence stars and mass-transferring massive binaries based on stellar evolution modelling. But self-consistent, dynamical simulations of very young GCs are usually not considered. In this work, we perform direct N-body modelling of such systems with total masses up to 3.2 × 105 M⊙, taking into account the observationally constrained primordial binary properties, and discuss the stellar mergers driven both by binary stellar evolution and dynamical evolution of GCs. The occurrence of stellar mergers is enhanced significantly in binary-rich clusters such that stars forming from the gas polluted by merger-driven ejection/winds would appear as MSPs. We thus emphasize that stellar mergers can be an important process that connects MSP formation with star cluster dynamics, and that multiple MSP formation channels can naturally work together. The scenario studied here, also in view of a possible top-heavy initial mass function, may be particularly relevant for explaining the high mass fraction of MSPs (the mass budget problem) and the absence of MSPs in young and low-mass star clusters.

Black Hole dynamics in Young Star Clusters

AbstractYoung star clusters are a promising environment for forming binary black holes. Such binaries may form dynamically or via binary star evolution or through the interplay of these two channels. To study these formation pathways, we have performed high precision direct N-body simulations of low-mass (M < 1000 M⊙) young star clusters. The simulations were carried out with the code Nbody6++GPU coupled with the population synthesis code MOBSE. Our results highlight the importance of dynamics to form massive black hole binaries even in low-mass young star clusters.

Deriving star cluster parameters with convolutional neural networks

Context. Convolutional neural networks (CNNs) have been proven to perform fast classification and detection on natural images and have the potential to infer astrophysical parameters on the exponentially increasing amount of sky-survey imaging data. The inference pipeline can be trained either from real human-annotated data or simulated mock observations. Until now, star cluster analysis was based on integral or individual resolved stellar photometry. This limits the amount of information that can be extracted from cluster images. Aims. We aim to develop a CNN-based algorithm capable of simultaneously deriving ages, masses, and sizes of star clusters directly from multi-band images. We also aim to demonstrate CNN capabilities on low-mass semi-resolved star clusters in a low-signal-to-noise-ratio regime. Methods. A CNN was constructed based on the deep residual network (ResNet) architecture and trained on simulated images of star clusters with various ages, masses, and sizes. To provide realistic backgrounds, M 31 star fields taken from The Panchromatic Hubble Andromeda Treasury (PHAT) survey were added to the mock cluster images. Results. The proposed CNN was verified on mock images of artificial clusters and has demonstrated high precision and no significant bias for clusters of ages ≲3 Gyr and masses between 250 and 4000 M⊙. The pipeline is end-to-end, starting from input images all the way to the inferred parameters; no hand-coded steps have to be performed: estimates of parameters are provided by the neural network in one inferential step from raw images.

Star clusters and young populations in the dwarf irregular galaxy Leo A

AbstractWe have studied young stellar populations and star clusters in the dwarf irregular galaxy Leo A using multicolor (B, V, R, I, Hα) photometry data obtained with the Subaru Suprime-Cam and two-color photometry results measured on archival HST/ACS F475W & F814W frames. The analysis of the main sequence (MS) and blue supergiant (BSG – “blue loop”) stars enabled us to study the star formation history in the Leo A galaxy during the last ~200 Myr. Also, we have discovered 5 low-mass (≲ 400 M⊙) star clusters within the ACS field. This finding, taking into account a low metallicity environment and a yet-undetected molecular gas in Leo A, constrains star formation efficiency estimates and scenarios. Inside the well known “hole” in the H i column density map (Hunter et al. 2012) we found a shock front (prominent in Hα), implying an unseen progenitor and reminding the “hole” problems widely discussed by Warren et al. (2011).

The effects of magnetic fields and protostellar feedback on low-mass cluster formation

We present a large suite of simulations of the formation of low-mass star clusters. Our simulations include an extensive set of physical processes -- magnetohydrodynamics, radiative transfer, and protostellar outflows -- and span a wide range of virial parameters and magnetic field strengths. Comparing the outcomes of our simulations to observations, we find that simulations remaining close to virial balance throughout their history produce star formation efficiencies and initial mass function (IMF) peaks that are stable in time and in reasonable agreement with observations. Our results indicate that small-scale dissipation effects near the protostellar surface provide a feedback loop for stabilizing the star formation efficiency. This is true regardless of whether the balance is maintained by input of energy from large scale forcing or by strong magnetic fields that inhibit collapse. In contrast, simulations that leave virial balance and undergo runaway collapse form stars too efficiently and produce an IMF that becomes increasingly top-heavy with time. In all cases we find that the competition between magnetic flux advection toward the protostar and outward advection due to magnetic interchange instabilities, and the competition between turbulent amplification and reconnection close to newly-formed protostars renders the local magnetic field structure insensitive to the strength of the large-scale field, ensuring that radiation is always more important than magnetic support in setting the fragmentation scale and thus the IMF peak mass. The statistics of multiple stellar systems are similarly insensitive to variations in the initial conditions and generally agree with observations within the range of statistical uncertainty.

Stability of multiplanetary systems in star clusters

Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $\Delta e \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability.

AChandraX-ray study of the young star cluster NGC 6231: low-mass population and initial mass function

NGC6231 is a massive young star cluster, near the center of the Sco OB1 association. While its OB members are well studied, its low-mass population has received little attention. We present high-spatial resolution Chandra ACIS-I X-ray data, where we detect 1613 point X-ray sources. Our main aim is to clarify global properties of NGC6231 down to low masses through a detailed membership assessment, and to study the cluster stars' spatial distribution, the origin of their X-ray emission, the cluster age and formation history, and initial mass function. We use X-ray data, complemented by optical/IR data, to establish cluster membership. The spatial distribution of different stellar subgroups also provides highly significant constraints on cluster membership, as does the distribution of X-ray hardness. We perform spectral modeling of group-stacked X-ray source spectra. We find a large cluster population down to ~0.3 Msun (complete to ~1 Msun), with minimal non-member contamination, with a definite age spread (1-8 Myrs) for the low-mass PMS stars. We argue that low-mass cluster stars also constitute the majority of the few hundreds unidentified X-ray sources. We find mass segregation for the most massive stars. The fraction of circumstellar-disk bearing members is found to be ~5%. Photoevaporation of disks under the action of massive stars is suggested by the spatial distribution of the IR-excess stars. We also find strong Halpha emission in 9% of cluster PMS stars. The dependence of X-ray properties on mass, stellar structure, and age agrees with extrapolations based on other young clusters. The cluster initial mass function, computed over ~2 dex in mass, has a slope Gamma~-1.14. The total mass of cluster members above 1 Msun is 2280 Msun, and the inferred total mass is 4380 Msun. We also study the peculiar, hard X-ray spectrum of the Wolf-Rayet star WR79.

Supernova enrichment of planetary systems in low-mass star clusters

The presence and abundance of short lived radioisotopes (SLRs) $^{26}$Al and $^{60}$Fe in chondritic meteorites implies that the Sun formed in the vicinity of one or more massive stars that exploded as supernovae (SNe). Massive stars are more likely to form in massive star clusters ($>$1000 M$_{\odot}$) than lower mass clusters. However, photoevaporation of protoplanetary discs from massive stars and dynamical interactions with passing stars can inhibit planet formation in clusters with radii of $\sim$1 pc. We investigate whether low-mass (50 - 200 M$_{\odot}$) star clusters containing one or two massive stars are a more likely avenue for early Solar system enrichment as they are more dynamically quiescent. We analyse $N$-body simulations of the evolution of these low-mass clusters and find that a similar fraction of stars experience supernova enrichment than in high mass clusters, despite their lower densities. This is due to two-body relaxation, which causes a significant expansion before the first supernova even in clusters with relatively low (100 stars pc$^{-3}$) initial densities. However, because of the high number of low mass clusters containing one or two massive stars, the absolute number of enriched stars is the same, if not higher than for more populous clusters. Our results show that direct enrichment of protoplanetary discs from supernovae occurs as frequently in low mass clusters containing one or two massive stars (>20 M$_{\odot}$) as in more populous star clusters (1000 M$_\odot$). This relaxes the constraints on the direct enrichment scenario and therefore the birth environment of the Solar System.

Initial Dynamical Evolution of Star Clusters with Tidal Field

AbstractObservations have been suggested that star clusters could form from the rapid collapse and violent relaxation of substructured distributions. We investigate the collapse of fractal stellar distributions in no, weak, and very strong tidal fields. We find that the rapid collapse of substructure into spherical clusters happens quickly with no or a weak tidal field, but very strong tidal fields prevent a cluster forming. However, we also find that dense Plummer spheres are also rapidly destroyed in strong tidal fields. We suggest that this is why the low-mass star clusters cannot survive near the galactic centre which has strong tidal field.

AN OPTICAL SURVEY OF THE PARTIALLY EMBEDDED YOUNG CLUSTER IN NGC 7129

NGC 7129 is a bright reflection nebula located in the molecular cloud complex near l=105.4, b=+9.9, about 1.15 kpc distant. Embedded within the reflection nebula is a young cluster dominated by a compact grouping of four early-type stars: BD+65 1638 (B3V), BD+65 1637 (B3e), SVS 13 (B5e), and LkH-alpha 234 (B8e). About 80 H-alpha emission sources brighter than V~23 are identified in the region, many of which are presumably T Tauri star members of the cluster. We also present deep (V~23), optical (VRI) photometry of a field centered on the reflection nebula and spectral types for more than 130 sources determined from low dispersion, optical spectroscopy. The narrow pre-main sequence evident in the color-magnitude diagram suggests that star formation was rapid and coeval. A median age of about 1.8 Myr is inferred for the H-alpha and literature-identified X-ray emission sources having established spectral types, using pre-main sequence evolutionary models. Our interpretation of the structure of the molecular cloud and the distribution of young stellar objects is that BD+65 1638 is primarily responsible for evacuating the blister-like cavity within the molecular cloud. LkH-alpha 234 and several embedded sources evident in near infrared adaptive optics imaging have formed recently within the ridge of compressed molecular gas. The compact cluster of low-mass stars formed concurrently with the early-type members, concentrated within a central radius of ~0.7 pc. Star formation is simultaneously occurring in a semi-circular arc some ~3 pc in radius that outlines remaining dense regions of molecular gas. High dispersion, optical spectra are presented for BD+65 1638, BD+65 1637, SVS 13, LkH-alpha 234, and V350 Cep. These spectra are discussed in the context of the circumstellar environments inferred for these stars.