Catálisis de muones en términos de química cuántica. Parte II: electrónica vs. enlace químico muón

Mnogabukaff que la química cuántica piensa en el principio de la catálisis de muones: cómo exactamente el muón reduce la temperatura del plasma deseado. En dos partes (la primera parte se puede leer aquí ).

La esencia de la segunda parte es simple: el muón es más pesado que el electrón, por lo que proporciona un enlace químico más fuerte y un enfoque más cercano de los núcleos, lo que reduce la temperatura del plasma requerida para encender la reacción termonuclear.

Pero aquellos que quieran ver las fórmulas, los gráficos y ver la esencia conceptual de la química cuántica aplicada a las moléculas más simples (cuasi), bienvenidos bajo cat.


Introduccion


En la primera parte (ver aquí ) examinamos la diferencia entre un átomo de hidrógeno  m a t h r m H  c d o t = m a t h r m p + e -  de su contraparte de muones pesados  m a t h r m p + m u -  : en el segundo caso, el muón estará atado con más fuerza y ​​se sentará a una distancia más cercana del protón. Al mismo tiempo, examinamos algunas cosas importantes que necesitaremos aquí (formas de orbitales y el sistema atómico de unidades).

En la segunda parte (es decir, aquí) trataremos de entender por qué, cómo y cuánto disminuirá la temperatura del plasma requerida para encender la reacción termonuclear. Las reacciones que nos interesan son:

 m a t h r mn H +mH rightarrow  textnuevosnúcleos+energía

úí


donde n, m = 1,2,3 corresponden a protón, deuterio y tritio, respectivamente. Naturalmente, estos núcleos tienen una carga positiva, por lo que si intenta acercarlos, comenzarán a repeler de acuerdo con la ley de Coulomb (consulte la parte anterior ), y esta es la barrera que impide la aparición de reacciones de fusión. Por cierto, en el caso de las reacciones de desintegración nuclear, esta repulsión tiene el papel opuesto, porque después de la separación del núcleo común, los fragmentos, que se repelen entre sí, adquieren energía cinética adicional, y es esta energía la que se calienta en las centrales nucleares.

Para superar esta barrera de Coulomb, se requiere un aumento en la temperatura del plasma ( T ), que, como todos recuerdan del curso escolar de MKT , está asociado con la velocidad promedio de las partículas en el plasma ( v ) por la fórmula

mv2=3k mathrmBT

donde m es la masa de partículas, y k mathrmB - Constante de Boltzmann .

Pero, imaginemos que hemos combinado dos núcleos de hidrógeno en una determinada partícula, donde ya están ubicados cerca y, por lo tanto, el resto de la barrera para ellos ya es muy pequeña. Entonces necesitaríamos acelerar significativamente estas partículas (léase: necesitamos temperaturas más bajas) para combinarlas en algo nuevo. Y tal papel debería jugar un ion intermedio ( mathrmnH mu− mathrmmH)+ , un análogo del ion de una molécula de hidrógeno  mathrmH+2=( mathrmHe−H)+ .
Habiendo examinado las diferencias entre estas dos partículas, nos daremos cuenta de cuán efectivo es el muón para reducir la temperatura de ignición de la fusión termonuclear.

Método MILK MO LKAO


Entonces, tenemos nuestro sistema molecular, que consiste en 2 núcleos de hidrógeno con una carga + e (un módulo de carga de electrones) y una partícula (electrón o muón) con una carga - e . Nuestro sistema, hasta que colisiona con otras partículas, está aislado y, por lo tanto, su energía puede descomponerse en sus partes constituyentes:

E=T( mathrmH1)+T( mathrmH2)+ underbraceT( mathrme−/ mu−)+V( mathrmH1 textfromH2)+V( mathrm mathrme−/ mu− text mathrmH1)+V( mathrm mathrme−/ mu− textk mathrmH2)E mathrme


donde los dos primeros términos ( T( mathrmH1) y T( mathrmH2) ) Es la energía cinética de los núcleos de hidrógeno, el tercer término ( T( mathrme−/ mu−) ) Es la energía cinética de una partícula negativa (electrón o muón), el cuarto término V( mathrmH1 textfromH2) Es la energía de la repulsión de Coulomb de los hidrógenos entre sí, y los dos restantes son la atracción de Coulomb del electrón / muón a cada uno de los protones. En el caso general, este es un problema de 3 cuerpos, solo un problema cuántico. Naturalmente, resolverlo en la frente es muy difícil. Pero, afortunadamente, los núcleos son al menos 1800 veces más pesados ​​que el electrón y 10 más pesados ​​que el muón, por lo que se moverán claramente más lentamente que las pequeñas partículas negativas. Debido a esto, primero puede resolver el problema a su vez: primero, encuentre la energía de los movimientos que no están relacionados con el movimiento de los núcleos, es decir, E mathrme y luego plena energía. Se ve así.

  1. Se selecciona la disposición de los núcleos de hidrógeno entre sí, y esto determina las interacciones de Coulomb entre ellos y con el electrón / muón. Potencial de Coulomb V(R)=k fracq1q2R depende solo de las cargas de partículas qi y la distancia entre ellos, por lo que para todos los isótopos de hidrógeno este valor será el mismo. Además, se resuelve el problema del movimiento de un electrón / muón en el campo de estos núcleos. Esta es la tarea de un solo cuerpo.
  2. Estas energías E mathrme se calculan para todos los arreglos posibles de los núcleos entre sí, y esta será la energía potencial de movimiento efectiva de los núcleos. En nuestro caso, necesitamos calcular las energías a diferentes distancias entre sí, por lo que el potencial para un par de núcleos es siempre unidimensional. Bueno, entonces solo necesitamos resolver el problema de dos cuerpos del movimiento de dos isótopos de hidrógeno uno con respecto al otro.

Obviamente, la raíz del problema con nosotros es el cálculo de la energía de electrones / muones en el campo de los núcleos. E mathrme . De hecho, este es el enlace químico: un cierto potencial que mantiene unidos los núcleos en ciertos lugares. Y esta misma tarea de encontrar la energía del enlace químico es la principal en la química cuántica.

Desafortunadamente, tanto el muón como el electrón son partículas cuánticas, por lo tanto, para encontrar esta energía, tenemos que recurrir a los métodos de la mecánica cuántica. De hecho, nuestro problema del movimiento de un electrón / muón en el campo de dos núcleos idénticos se resuelve explícitamente (ver aquí ), pero esta solución es muy complicada y el resultado no es tan claro como en el caso de un átomo similar al hidrógeno. Por lo tanto, intentaremos desmontar un enfoque diferente y aproximado, que sea aplicable a cualquier sistema. Este es el llamado método de orbitales moleculares como combinaciones lineales de orbitales atómicos, o MO LKAO.

Echemos un vistazo más de cerca a la ecuación de Schrödinger para el movimiento de un electrón / muón en el campo de los núcleos de hidrógeno:

 hatH psi= underbrace left( overbrace− frac12m( frac partial2 partialx2+ frac partial2 partialy2+ frac partial2 partialz2) hatT+ overbrace− frac1R1 hatV1+ overbrace− frac1R2 hatV2+ overbrace frac1R hatV mathrmHH right) hatH psi=E psi


Esta ecuación se escribió en el sistema atómico de unidades (ver PS en la parte anterior ), por lo tanto, la carga del núcleo de hidrógeno y el electrón / muón es +1, --1, respectivamente, la masa del electrón es m = 1, y para el muón m ≈207.

Y si observa más de cerca, puede ver que en el Hamiltoniano puede seleccionar una pieza conectada puramente con el movimiento de una partícula negativa alrededor de solo uno de los núcleos, que es solo el Hamiltoniano del átomo de hidrógeno, y esto se puede hacer de 2 maneras:

 hatH=( overbrace hatT+ hatV1 hatH1+ hatV2+ hatV mathrmHH)=( overbrace hatT+ hatV2 hatH2+ hatV1+ hatV mathrmHH)


Fuera del hamiltoniano de un átomo similar al hidrógeno (  hatHi, i=$1. ) siempre tenemos 2 piezas: la energía de interacción de un electrón / muón con otro núcleo (  hatVj ) y energía de repulsión nuclear (  hatV mathrmHH ) El segundo de ellos no afecta en absoluto el movimiento de los electrones: es solo un cambio de energía en una cierta cantidad, pero la interacción de un electrón con otro núcleo es algo importante.
Podemos imaginar que en cualquier momento nuestra partícula gira solo alrededor de uno de los núcleos, y la interacción con el segundo es solo una corrección. Como método de rotación alrededor de uno de los núcleos, podemos suponer que el electrón / muón está en el estado fundamental (1s), cuya función de onda nos es bien conocida en la parte anterior:

|1s rangle= frac1 sqrt pi exp left(− fracRR1 right)


donde R1 Es el radio de Bohr para una partícula. En el caso de un electrón R1=1 Boro (que es el radio de Bohr para un electrón, igual a aproximadamente 0.5 angstroms), y en el caso de un muón R1= frac1m mu aprox frac1207 .
Para aproximar de alguna manera la función de onda de electrones / muones en el campo de 2 núcleos, podemos intentar tomar la siguiente representación:

 psi aproxc1|1s1 rangle+c2|1s2 rangle


y luego el problema de resolver una ecuación diferencial parcial compleja con nosotros se reduce a la búsqueda de 2 coeficientes desconocidos c 1 y c 2 . Este es el orbital muy molecular presentado como la suma de los coeficientes (una combinación lineal de los científicos) orbitales atómicos 1s.

Naturalmente, necesitamos una ecuación para estos parámetros. Y obtenerlo es bastante simple si sustituyes esta aproximación en la ecuación de Schrödinger  hatH psi=E psi :

 hatH(c1|1s1 rangle+c2|1s2 rangle)=c1 hatH|1s1 rangle+c2 hatH|1s2 rangle=E(c1|1s1 rangle+c2|1s2 rangle)=c1E|1s1 rangle+c2E|1s2 rangle


En realidad, queremos que esta proporción se satisfaga en todas partes, por lo que de alguna manera podemos calcular los valores promedio de todo esto. Multiplicamos esta ecuación a la izquierda por <1s1| y < 1 s 2 | e integrarse sobre todas las coordenadas. Como resultado, obtenemos un sistema de 2 ecuaciones lineales, donde es necesario encontrar los coeficientes c 1 , c 2 y energía E :

\ begin {pmatrix} \ langle 1s_1 | \ hat {H} | 1s_1 \ rangle & \ langle 1s_1 | \ hat {H} | 1s_2 \ rangle \\ \ langle 1s_2 | \ hat {H} | 1s_1 \ rangle & \ langle 1s_2 | \ hat {H} | 1s_2 \ rangle \\ \ end {pmatrix} \ begin {pmatrix} c_1 \\ c_2 \ end {pmatrix} = E \ begin {pmatrix} \ langle 1s_1 | 1s_1 \ rangle & \ langle 1s_1 | 1s_2 \ rangle \\ \ langle 1s_2 | 1s_1 \ rangle & \ langle 1s_2 | 1s_2 \ rangle \\ \ end {pmatrix} \ begin {pmatrix} c_1 \\ c_2 \ end {pmatrix}


Cualquiera que haya estudiado álgebra lineal reconocerá un problema generalizado de autovector-autovalor. Antes de resolverlo, analizaremos a qué son iguales los elementos de las 2 matrices, la tuta existente (y al mismo tiempo presentamos su breve designación con una letra).

  • Comencemos con el más simple:  langle1s1|1s1 rangle= langle1s2|1s2 rangle=1 - esta es la normalización de las funciones de onda, y como recordamos, la probabilidad total de encontrar un electrón / muón es al menos 1 en alguna parte.
  •  langle1s1|1s2 rangle= langle1s2|1s1 rangle=S - este es el llamado solapamiento integral, que muestra cómo se superponen las nubes de electrones para cada uno de los átomos.
  •  langle1s1| hatH|1s1 rangle= langle1s2| hatH|1s2 rangle= alpha . Esta integral consta de varias partes:

     langle1s1| hatH|1s1 rangle= underbrace langle1s1| hatH1|1s1 rangle− fracm2+ langle1s1| hatV2|1s1 rangle+ frac1R


  •  langle1s1| hatH|1s2 rangle= langle1s2| hatH|1s1 rangle= beta . Aquí es similar:

    \ langle 1s_2 | \ hat {H} | 1s_1 \ rangle = \ underbrace {\ langle 1s_2 | \ overbrace {\ hat {H} _1 | 1s_1 \ rangle} ^ {- \ frac {m} {2} | 1s_1 \ rangle}} {{\ frac {m} {2} S} + \ langle 1s_2 | \ sombrero {V} _2 | 1s_1 \ rangle + \ frac {S} {R}


    es decir la energía de un átomo similar al hidrógeno y la repulsión internuclear, escalada por la superposición integral (primer y último término) y, por así decir, la energía del salto de electrones / muones de un átomo a otro.

Encontremos las expresiones para las energías de nuestro ion similar al hidrógeno de la ecuación reescrita como

\ begin {pmatrix} \ alpha & \ beta \\ \ beta & \ alpha \ end {pmatrix} \ begin {pmatrix} c_1 \\ c_2 \ end {pmatrix} = E \ begin {pmatrix} 1 & S \\ S & 1 \ end {pmatrix} \ begin {pmatrix} c_1 \\ c_2 \ end {pmatrix}


Para encontrar la energía que necesitas para resolver la ecuación:

\ det \ begin {pmatrix} \ alpha -E & \ beta -ES \\ \ beta - ES & \ alpha -E \ end {pmatrix} = (\ alpha -E) ^ 2 - (\ beta - ES) ^ 2 = 0


donde "det" denota el determinante (determinante de una matriz, en ruso).

Las soluciones de esta ecuación cuadrática con respecto a E son:

E pm= frac alpha pm beta1 pmS=− fracm2+ frac1R+ frac langle1s1| hatV2|1s1 rangle pm langle1s2| hatV2|1s1 rangle1 pmS


La primera pieza es obviamente la energía del átomo, la segunda es la repulsión internuclear, la misma barrera de Coulomb que evita la ignición de la reacción termonuclear, y la última estructura compleja debe ser tratada.

Si descartamos la repulsión internuclear, que es solo un punto de referencia para la energía de electrones / muones, obtenemos que tenemos dos estados con energía

 epsilon pm=− fracm2+ frac langle1s1| hatV2|1s1 rangle pm langle1s2| hatV2|1s1 rangle1 pmS


Dado que ambas funciones de onda |1s1 rangle y |1s2 rangle - positivo, y  hatVi<0 (porque la partícula negativa siempre se atrae a la positiva), entonces  epsilon+<− fracm2 (energía de un solo átomo), y  epsilon−>− fracm2 es decir obtenemos una imagen estándar de los orbitales moleculares:

imagen

Orbital inferior con energía E+ llamado enlace, y la parte superior (con energía E− ) - anti-vinculante, o aflojamiento. Como resultado, si un electrón / muón se encuentra en el orbital molecular más bajo, entonces se beneficia al volar alrededor de 2 núcleos que alrededor de uno, y con su movimiento disminuye la energía total del sistema. Y este es el enlace químico muy mágico que detecta la repulsión internuclear, permitiendo que los núcleos estén uno al lado del otro durante bastante tiempo.

Y aquí las integrales del enlace químico deben calcularse para comprender cuán cerca se permite que estén los núcleos de hidrógeno. De hecho, las tres integrales buscadas se calculan analíticamente, pero es terriblemente hemorroidal y complicada (cualquiera que esté interesado, vea el Capítulo 9 en el libro de Química Cuántica de Flary ). Por lo tanto, iremos de una manera diferente, más simple, y calcularemos estas integrales numéricamente usando el método de Monte Carlo.

Método Metrópolis


Lo veo muy lógico en el texto sobre la energía termonuclear para rendir homenaje a su abuelo: el átomo militar y, más específicamente, el proyecto Manhattan . Fue de él que creció el método Monte Carlo, y en particular el algoritmo Metropolis , uno de los autores, Edward Teller, es el "padre de la bomba de hidrógeno" (es decir, la persona que lanzó la fusión termonuclear en el atolón Eniwetok ).


En general, analizaremos la esencia del método. Está destinado a las tareas de mecánica estadística. La distribución principal es la distribución de Boltzmann: la probabilidad de detectar un sistema en un determinado estado es  exp(− betaE) ,  beta−1=k mathrmBT . Y el valor observado de algún parámetro A para el sistema en equilibrio termodinámico es igual a la integral

 langleA rangle= frac1Z intA(q) exp(− betaE(q))dq


donde q son las coordenadas que parametrizan el estado del sistema (por ejemplo, las coordenadas / momentos de las partículas), y Z es el factor de normalización llamado función de partición:

Z= int exp(− betaE(q))dq


Si hay demasiadas partículas en el sistema, entonces no contar ninguna de las integrales en la frente es completamente poco realista. El ingenuo método de Monte Carlo, en el que simplemente seleccionamos un montón de coordenadas q aleatorias, tampoco dará nada significativo si realmente hay estados posibles del sistema para los cuales la probabilidad  exp(− betaE) notablemente no cero, muy pocos. Y es precisamente para tales casos que necesitamos una muestra por significación, en la que permitimos que el algoritmo muestree solo lugares suficientemente probables en el espacio de estado.

El algoritmo Metropolis se ve de la siguiente manera.

  1. Al iniciar la simulación, seleccionamos alguna aproximación inicial en el espacio de configuración  mathbfq(0) y algún vector del incremento máximo posible  delta mathbfq . En el punto de partida, calculamos la energía del sistema. E(0)=E( mathbfq(0)) (leer - probabilidad p= exp(− betaE(0)) )
  2. La nueva configuración en el enésimo paso es la siguiente.
    1. Calcule la energía de la configuración de prueba E mathrmtrial=E( mathbfq mathrmtrial) (es decir, probabilidad p mathrmtrial= exp(− betaE mathrmtrial) )
    2. Y luego comparamos la vieja probabilidad p(n) con juicio p mathrmtrial

      • si la nueva configuración tiene una probabilidad mayor o igual (  fracp mathrmpruebap(n) geq1 ), o, de manera equivalente, la energía del nuevo punto es menor o igual que en el antiguo ( E mathrmprueba leqE(n) ), se acepta el nuevo punto y el sistema entra en él ( q(n+1)=q mathrmprueba ),
      • si la configuración de prueba es más alta en energía ( E mathrmprueba>E(n) ), que es equivalente  fracp mathrmpruebap(n)<1 , entonces en este caso generamos un número aleatorio P en[0;1) de una distribución uniforme, y compárelo con la razón de las probabilidades, que son las probabilidades de transición. Si P< fracp mathrmpruebap(n) , entonces aceptamos un nuevo punto, y si no ( P geq fracp mathrmpruebap(n) ), luego rechazamos y el sistema permanece en la configuración anterior ( q(n+1)=q(n) ) ...


  3. Tomando muchos pasos de acuerdo con el algoritmo anterior, tomamos muestras de una parte significativa (es decir, realmente importante) del posible espacio de configuraciones del sistema. La integral de interés para nosotros se calcula mediante la fórmula:
     langleA rangle= frac1Z intA( mathbfq) exp(− betaE( mathbfq))d mathbfq= frac1N sumNn=0A( mathbfq(n))

Así funciona el algoritmo Metropolis.

Y ahora sería necesario adaptarlo al cálculo de las 3 integrales que nos interesan. Echemos un vistazo a ellos con más detalle.

  • S(R)= langle1s2|1s1 rangle= int limits+ infty− infty int limits+ infty− infty int limits+ infty− infty overbrace underbrace frac1 sqrt pi exp(−m underbrace| mathbfr− mathbfr2|R2)1s2A( mathbfr) cdot overbrace underbrace frac1 sqrt pi exp(−m underbrace| mathbfr− mathbfr1|R1)1s1p( mathbfr)dxdydz donde  mathbfr=(x,y,z) mathbfT - coordenadas del electrón / muón,  mathbfri=(xi,yi,zi) mathbfT Son las coordenadas de los núcleos de hidrógeno, y Ri=| mathbfr− mathbfri|= sqrt(x−xi)2+(y−yi)2+(z−zi)2 - distancias entre partículas positivas y negativas,
  •  langle1s1| hatV2|1s1 rangle=− int limits+ infty− infty int limits+ infty− infty int límites+ infty− infty overbrace underbrace frac1 sqrt pi exp(−mR1)1s1 frac1R2A( mathbfr) cdot overbrace underbrace frac1 sqrt pi exp(−mR1)1s1p( mathbfr)dxdydz
  •  langle1s2| hatV2|1s1 rangle=− int limits+ infty− infty int limits+ infty− infty int límites+ infty− infty overbrace underbrace frac1 sqrt pi exp(−mR2)1s2 frac1R1A( mathbfr) cdot overbrace underbrace frac1 sqrt pi exp(−mR1)1s1p( mathbfr)dxdydz

Se ve que, si calculamos la función 1s de uno de los átomos para la probabilidad p ,

hacerlo, por supuesto, no es muy bueno,
porque la densidad de probabilidad es el módulo de la función de onda cuadrada El | p s i | 2  , no la función de onda en sí  p s i .

entonces todo lo demás bajo el signo de la integral (la función de segunda onda y en 2 de 3 casos el potencial de atracción del electrón / muón al núcleo) será una función cuyo valor promedio se calcula. Lo único que tendrá que hacerse, en contraste con el cálculo habitual por el método de Metrópolis, es enderezar la normalización de las integrales. El hecho es que la normalización estándar estará activada

Z= int limits+ infty− infty int limits+ infty− infty int limits+ infty− infty exp(−mR)dxdydz=4 pi int limits+ infty0 exp(−mR)R2dR= frac8 pim3


Y necesitamos normalización para  sqrt langle1s1|1s1 rangle donde

 langle1s1|1s1 rangle= int limits+ infty− infty int limits+ infty− infty int limits+ infty− infty exp(−2mR)dxdydz=4 pi int limits+ infty0 exp(−2mR)R2dR= frac pim3$


Esto significa que cada integral calculada de acuerdo con la metrópoli deberá multiplicarse por un factor

 fracZ sqrt langle1s1|1s1 rangle=8 sqrt frac pim3



Esto ya se puede organizar en forma de una secuencia de comandos determinada, por ejemplo, en Python (por ejemplo, el código está debajo).

Por ejemplo, así.
import numpy as np from math import * # r = 0...+infty # phi = 0...2pi # theta = 0...pi # function to convert spherical coordinates into Cartesian def sph2cart(r, phi, theta): xyz=np.zeros(3) xyz[0]=r*cos(phi)*sin(theta) xyz[1]=r*sin(phi)*sin(theta) xyz[2]=r*cos(theta) return xyz # Distance between vectors r1 and r2 def dist(r1, r2): return sqrt(np.dot(r1-r2, r1-r2)) # re -- Cartesian coordinates of electron # rn -- Cartesian coordinates of nucleus # psi_1s returns a value of 1s wavefunction def psi_1s(re, rn, scale=1.0): ren=dist(re,rn) return scale**(3/2)*exp(-scale*ren)/(sqrt(pi)) ######################################## ############ Settings ################## ######################################## NumPtsPerIntegral=100000 # Number of points per integral... duuuh #mass=1.0 # mass of the particle (electron) mass=207.0 # mass of the particle (muon) AllRab=[] #AllRab+=[1.0*(0.1)**n for n in range(0,10) ] AllRab+=[ (1.4+0.25*n)/mass for n in range(0,10)] print(AllRab) ######################################## ######################################## ######################################## dumpster=open("res.dat", "w") # output file to store results of the simulation dx=2.0/mass # maximal increment for the coordinate rna=np.array([0.0, 0.0, 0.0]) # position of nucleus "a" renorm=8.0*sqrt(pi/mass**3) # factor to readjust result from incorrect norm of Metropolis weighting to a correct 1s wavefunction norm # loop for the potential energy calculation at the chosen internuclear distances for npt,Rab in enumerate(AllRab): Norm=0.0 # <1s_a | 1s_a > for check Sab=0.0 # <1s_a | 1s_b > Vaa=0.0 # <1s_a | |r - R_b|**(-1) | 1s_a > Vab=0.0 # <1s_a | |r - R_b|**(-1) | 1s_b > re=np.array([1.0/mass, 0.0, 0.0]) # initial position of the electron rnb=np.array([Rab, 0.0, 0.0]) # position of nucleus "b" NumAcc=0.0 # Number of accepted points for i in range(0,NumPtsPerIntegral): # loop for Metropolis algorithm newre=re+np.random.uniform(low=-dx, high=dx, size=3) # trial position of electron pnew=psi_1s(newre, rna, scale=mass) ## trial probability pold=psi_1s(re, rna, scale=mass) ## previous probability due to dumb and ineffective realization if pnew/pold >= np.random.random(): ## importance sampling step re=newre NumAcc+=1. Norm+=psi_1s(re, rna, scale=mass) Sab+=psi_1s(re, rnb, scale=mass) Vaa+=psi_1s(re, np.zeros(3), scale=mass)/dist(re, rnb) Vab+=psi_1s(re, rnb, scale=mass)/dist(re, rnb) Norm*=renorm/NumPtsPerIntegral Sab*=renorm/NumPtsPerIntegral Vaa*=renorm/NumPtsPerIntegral Vab*=renorm/NumPtsPerIntegral def s_test(x,scale=1.0): # this is an analytical expression for overlap integral S in case of 1s hydrogen wavefunctions return exp(-x*scale)*(1.+x*scale+(1./3.)*(x*scale)**2) #E=-0.5*mass**2 + 1./sqrt(np.dot(rnb,rnb)) - (Vaa + Vab)/(1.0 + Sab) # full energy E= 1./sqrt(np.dot(rnb,rnb)) - (Vaa + Vab)/(1.0 + Sab) # energy adjusted to energy of a single atom as the dissociational limit dumpster.write(" %10.3e %40.10f %15.10f %15.10f %15.5f %15.5f\n" % (Rab, E, Sab, s_test(Rab,scale=mass), 100.0*NumAcc/NumPtsPerIntegral, Norm )) dumpster.flush() dumpster.close() 



Usando tales cálculos, finalmente podemos comparar las energías potenciales en el ion hidrógeno  mathrmH+2 y su homólogo muón.

 mathrmH+2=p+e−p+ vs.  mathrmp+ mu−p+


Entonces, armados con un script, podemos calcular la superficie de la energía potencial del enfoque de los núcleos de hidrógeno unidos por un electrón y un muón. Como punto de referencia de energía, tomamos átomos infinitamente diluidos entre sí (es decir, −m/2 , que es igual al potencial a la distancia entre los núcleos R=+ infty )

En el caso de un electrón, el potencial cercano al mínimo se ve así:



El mínimo se produce a una distancia de aproximadamente 2 boro (es decir, aproximadamente la suma de 2 radios atómicos), y la energía de disociación de la molécula en fragmentos es de aproximadamente 0,06 Hartree, que corresponde al calentamiento a unos 20,000 grados Kelvin (o Celsius, no importa aquí). Para convertir energías, recomiendo usar recursos en línea, como este .

Una situación similar con un ión de hidrógeno unido muonalmente:



Como el radio de Bohr para el muón hidrógeno es más pequeño (ver la parte anterior ), los núcleos de hidrógeno también se sientan unas 200 veces más cerca de la energía potencial mínima. La energía de descomposición de esta molécula ya es más de 10 Hartree, que corresponde a una temperatura de más de tres grados de lyam (  aprox(3.2 cdot106) circ )

Para la ignición, las reacciones generalmente requieren una temperatura del orden de 10 8 K , que es aproximadamente 320 Hartree. Veamos a qué distancias se logra una energía similar en el caso del ion divodorón ordinario y en el caso de su versión muón:



En el caso de los primeros, esto corresponde a una distancia de aproximadamente 0.0058 Boro (línea vertical).

Se alcanza una distancia similar en el muón hidrógeno a una energía de aproximadamente 190 Ha, es decir aproximadamente una vez y media menos. Y esta es la estimación más simple de la temperatura de la catálisis de muones.

Pero, de hecho, todo será aún más genial. El hecho es que si se forma una partícula estable  m a t h r m (m H ( m u - ) n H ) + , entonces estos núcleos, mientras el muón está vivo, oscilarán entre sí. Y aquí puede producirse un túnel desde el estado de "dos átomos de hidrógeno" al estado de "núcleo más pesado", y la probabilidad de túnel depende de la longitud de túnel requerida d aproximadamente como p - d , de modo que al acercar los dos núcleos mediante el muón, aumentaremos en gran medida la probabilidad del curso de túnel de esta reacción. Desafortunadamente, las estimaciones de este efecto ya no requieren química cuántica, sino física nuclear, por lo que esta parte de la discusión está más allá del alcance de esta publicación. Entonces en esto nos detendremos.

PD: ¿Por qué no es tan simple?


De hecho, formar estas partículas no es tan simple en condiciones de plasma. El hecho es que si colisionamos dos partículas, entonces su energía total obviamente excede la energía de disociación (o ionización, en el caso de un núcleo + electrón / muón), por lo que cuando colisionan, no forman una partícula estable (átomo, ion, molécula), sino que pasan volando uno al lado del otro. Para mantenerse unidos, necesitan arrojar energía excedente en algún lugar, y para esto necesitamos un tercer extra que tomará esta energía. Puede ser un fotón, o algún tipo de partícula izquierda volando cerca, pero lo principal es que las condiciones deberían contribuir a este arrastre de exceso de energía.

PPS


Si tiene comentarios / aclaraciones / preguntas, escriba comentarios o PM. Corregiré todo, responderé y explicaré todo.

Source: https://habr.com/ru/post/445672/


All Articles