Astronomy
Colisión de galaxias
Explora un encuentro de galaxias: inclina los discos, sigue una estrella y compara rotaciones opuestas. Los puentes de marea y las colas emergen de las trayectorias de las estrellas de prueba. Los halos oscuros opcionales y la fricción dinámica ilustran cómo los núcleos pueden fusionarse.
Qué ocurre
Las estrellas rara vez chocan directamente durante un encuentro galáctico. La gravedad diferencial estira los discos estelares en puentes y colas. Las estrellas son trazadores sin masa, sin autogravedad del disco o gas. Con los halos oscuros habilitados, los potenciales rígidos de Plummer profundizan el campo; la fricción dinámica opcional elimina la energía orbital. Los núcleos se unen en una pequeña separación prescrita. Esto ilustra un mecanismo de fusión, no un tiempo de fusión calibrado o una predicción para la Vía Láctea y Andrómeda. La energía estelar positiva es solo un diagnóstico instantáneo en el potencial cambiante.
Ecuaciones del módulo
Ecuaciones unidades: G = 1, Masa primaria 1, radio del disco 1; masa del intruso q Radiador y disco de radio sqrt(q); Centro de masa en reposo Núcleos: potencial suavizado, epsilon = 0.3; mu = q/(1+q); E = mu*v²/2 - q/sqrt(r²+epsilon²) comienzo: D0 = 8; v_escape = sqrt(2*(1+q)/sqrt(D0²+epsilon²)); v = f*v_escape; E0 = q*(f-1)*(f+1)/sqrt(D0²+epsilon²) Sin halos: f < 1 significa ligado; f = 1 es el límite de escape; f > 1 significa no ligado. Las órbitas suavizadas no son cónicas exactas de Kepler; una órbita ligada por sí sola no garantiza la coalescencia. Punto de giro: h²/(2*r_p²) - (1+q)/sqrt(r_p²+epsilon²) = E/mu; h = b*v estrellas: 2100 partículas sin masa por disco a resolución base, escaladas por el multiplicador de partículas; 20 anillos; las velocidades iniciales sqrt(M/r) utilizan una aproximación de punto-masa y la suavización 0.05 se utiliza en las fuerzas estelares. integración: kick-drift-kick leapfrog dt = 0.005; el pronóstico utiliza el mismo integrador nuclear y muestra su trayectoria para mostrar Fracción de cola: distancia > 2 en unidades de radio de disco primario desde el núcleo de nacimiento; esto es un diagnóstico geométrico, no una prueba de escape Sin halos, energía estelar: v²/2 - 1/sqrt(r1²+0.05²) - q/sqrt(r2²+0.05²), en el marco baricéntrico.La energía positiva no garantiza el escape permanente de un potencial dependiente del tiempo. Encuentros estelares (estimación geométrica sin enfoque gravitatorio): n = 10^11/(pi*15000²*300) pc^(-3); radio solar = 2.3e-8 pc; sigma = pi*(2*2.3e-8)²; trayectoria libre media = 1/(n*sigma), sobre 3.19e14 pc Fricción: a = −4π ln Λ M ρ(r) [erf X − 2X/√π e^(−X²)] v/|v|³, X = |v|/(√2 σ) ; ρ = 3M_h/(4π a³)(1 + r²/a²)^(−5/2), σ² = M_h/(6√(r² + a²)) Con halos: E_pair = mu v²/2 - q/sqrt(r²+epsilon²) - 2 eta q/sqrt(r²+a²) - eta² q/sqrt(r²+2a²); mu = q(1+eta)/(1+q) El potencial estelar también contiene -eta/sqrt(r1²+a²) - eta q/sqrt(r2²+a²). eta = M_h/M_nucleus; a = 1; R_merge = 0.3
