1. 丝状分子云中Bonnor-Ebert球体塌缩的物理基础Bonnor-Ebert球体作为描述星际介质中自引力平衡的经典模型其稳定性判据最早由Ebert(1955)和Bonnor(1956)独立建立。这个理论框架描述了在等温条件下由气体压力与自引力维持平衡的球对称云核。当考虑丝状分子云环境时情况变得更为复杂——丝状结构的轴向对称性和外部压力梯度会显著改变传统Bonnor-Ebert球体的稳定性条件。在典型的星际介质条件下Bonnor-Ebert球体的临界质量可以表示为M_BE ≈ 1.18 c_s^4 / (G^(3/2) P_ext^(1/2))其中c_s为等温声速G为引力常数P_ext为外部压力。这个经典公式假设球体处于均匀外部压力环境中。然而当球体嵌入丝状结构时压力梯度会沿径向变化形成各向异性的力学环境。我们的模拟显示这种情况下需要考虑三个关键压力分量(1)核心内部热压力(2)丝状结构产生的径向压力梯度(3)轴向扰动导致的压力不均匀性。2. 数值模拟方法与初始条件设置我们采用基于RAMSES代码的自适应网格细化(AMR)方法进行三维流体动力学模拟耦合了自引力、气体动力学和简化的辐射冷却过程。计算域设置为边长5pc的立方体最高分辨率达到0.005pc(约1000AU)足以解析原恒星核心的形成过程。初始条件构建遵循以下步骤建立轴向均匀的丝状结构密度剖面采用Plummer-like分布 ρ(r) ρ_c / [1 (r/R_flat)^2]^p 其中ρ_c5×10^−20g/cm^3R_flat0.1pcp2.5与观测到的典型分子云丝参数一致沿轴向引入线质量扰动扰动幅度ϵ满足 f_cyl (1 ϵ) × f_cyl0 其中f_cyl0为背景线质量扰动波长λ在0.3-3pc范围内变化设置等温状态方程温度T10K对应声速c_s0.19km/s添加1%级别的随机密度扰动模拟湍流涨落特别值得注意的是我们采用了两种不同的扰动实现方式对于小扰动(ϵ0.3)使用传统的密度扰动而对于大扰动(ϵ≥0.3)则直接修改线质量分布后者更符合实际分子云中观察到的强烈非均匀性。3. 塌缩阈值与非线性演化特征模拟结果揭示了丝状环境中Bonnor-Ebert球体塌缩的三个关键特征尺度主导波长(λ_dom)线性扰动理论预测的最快增长模式在我们的参数下约为0.6pc临界波长(λ_crit)约为λ_dom/2低于此尺度的扰动理论上应该被抑制BE修正尺度(λ_BEλ_dom/2.61)我们发现的实际塌缩最小间距当线质量参数f_cyl超过0.66时系统表现出显著的非线性行为。图7展示了核心质量(实线)与可支撑的BE质量(虚线)随f_cyl的变化关系。值得注意的是在f_cyl0.66时系统行为与经典线性理论预测一致当0.66≤f_cyl≤0.89时局域质量超过BE极限即可触发塌缩即使间距小于λ_crit对于f_cyl0.89的情况整个丝状结构变得径向不稳定关键发现塌缩由局部BE质量判据主导而非整体线质量稳定性。这意味着高线质量丝状结构中可以形成比传统预测更密集的恒星团。4. 塌缩动力学的时域特征通过追踪最大线质量随时间演化(图10)我们识别出两种不同的塌缩模式直接塌缩当初始扰动波长λ_init满足λ_BEλ_initλ_crit时核心在约0.2Myr内快速坍缩典型表现为轴向速度场出现明显汇聚流(v_z0.5km/s)中心密度呈指数增长(ρ∝e^(t/τ)τ≈0.05Myr)形成明确的激波前沿振荡延迟塌缩当λ_initλ_BE时系统会经历以下阶段 (1) 初始声波振荡阶段(持续≈0.8Myr) (2) 主导模式增长阶段(≈0.5Myr) (3) 最终在λ_dom尺度上塌缩这种差异源于压力平衡的重建时间尺度。当扰动过于紧凑时核心需要更长时间调整密度剖面以满足BE平衡条件。5. 天体物理意义与观测验证我们的研究为解释以下观测现象提供了理论依据密集恒星形成区如Orion A分子云中观测到的0.1-0.2pc核心间距远小于经典预测的λ_dom≈0.6pc分层碎裂结构在SDC13等丝状云中观察到的丝中丝现象可能反映不同尺度的BE塌缩过程年轻星团的高空间密度如NGC 6334中每立方秒差距数十颗原恒星的极端情况观测检验建议通过ALMA等高分辨率设备测量核心间距分布寻找λ_dom/2.61的特征尺度结合HCO(3-2)等示踪分子分析速度场识别非线性塌缩的动力学特征利用Herschel等远红外数据重建质量分布验证f_cyl阈值效应6. 模型局限性与未来方向当前研究存在几个需要改进的方面磁场效应未考虑磁场对塌缩的抑制作用实际临界f_cyl可能更高辐射反馈新生恒星的辐射场会改变局部热力学状态更真实的初始条件需要从大尺度湍流模拟自洽地产生丝状结构未来工作将聚焦于引入非理想磁流体动力学(MHD)效应耦合简化的辐射传输计算发展新型示踪算法更好区分不同形成机制的恒星这些改进将帮助我们更准确地预测恒星形成效率及初始质量函数为星系演化模型提供更可靠的物理基础。