Aquí, , y son constantes. Los dos términos del lado derecho constituyen una repulsión y una atracción, porque sus primeras derivadas con respecto a son negativas y positivas, respectivamente.
Como se explica en el artículo original de Buckingham y, por ejemplo, en la sección 2.2.5 del texto de Jensen, [2] la repulsión se debe a la interpenetración de las capas electrónicas cerradas . "Por lo tanto, existe cierta justificación para elegir la parte repulsiva (del potencial) como una función exponencial ". El potencial de Buckingham se ha utilizado ampliamente en simulaciones de dinámica molecular .
Como el término exponencial converge a una constante cuando → , mientras que el término diverge, el potencial de Buckingham se vuelve atractivo cuando se vuelve pequeño. Esto puede ser problemático cuando se trata de una estructura con distancias interatómicas muy cortas, ya que todos los núcleos que cruzan un cierto umbral se unirán fuertemente (y de manera no física) entre sí a una distancia de cero. [2]
Potencial Buckingham modificado (Exp-Six)
El potencial de Buckingham modificado, también llamado potencial "exp-seis", se utiliza para calcular las fuerzas interatómicas de los gases basándose en la teoría de colisión de Chapman y Cowling. [3] El potencial tiene la forma
donde es el potencial interatómico entre el átomo i y el átomo j, es la energía potencial mínima, es la medida de la pendiente de la energía repulsiva que es la relación , es el valor de donde es cero, y es el valor de que puede alcanzar el potencial interatómico mínimo . Esta función de potencial solo es válida cuando , ya que el potencial decaerá hacia . Esto se corrige identificando , que es el valor de en el que se maximiza el potencial; cuando , el potencial se establece en infinito.
Potencial de Coulomb-Buckingham
El potencial de Coulomb-Buckingham es una extensión del potencial de Buckingham para su aplicación en sistemas iónicos (por ejemplo, materiales cerámicos ). La fórmula para la interacción es
La ecuación anterior puede escribirse en su forma alternativa como
donde es la distancia mínima de energía, es un parámetro adimensional libre y es la profundidad de la energía mínima.
Potencial de Beest Kramer van Santen (BKS)
El potencial BKS es un campo de fuerza que se puede utilizar para simular el potencial interatómico entre átomos de vidrio de sílice . [4] En lugar de confiar solo en datos experimentales, el potencial BKS se deriva combinando métodos de química cuántica ab initio en pequeños grupos de sílice para describir la interacción precisa entre los vecinos más cercanos, que es la función del campo de fuerza preciso . Los datos experimentales se aplican para ajustar la información de fuerza a mayor escala más allá de los vecinos más cercanos. Al combinar la información microscópica y macroscópica , la aplicabilidad del potencial BKS se ha extendido tanto a los polimorfos de sílice como a otros sistemas de óxidos de red tetraédrica que tienen la misma estructura de grupo, como aluminofosfatos, carbono y silicio .
La forma de este potencial interatómico es la forma Buckingham habitual, con la adición de un término de fuerza de Coulomb . La fórmula para el potencial BKS se expresa como
donde es el potencial interatómico entre el átomo i y el átomo j, y son las magnitudes de las cargas, es la distancia entre átomos, y , y son parámetros constantes basados en el tipo de átomos. [5]
Los parámetros potenciales BKS para átomos comunes se muestran a continuación: [5]
Parámetros BKS
yo
Un ij (eV)
B ij (Å −1 )
C ij (eV•Å 6 )
O-O
1388.7730
2.76000
175.0000
O-Si
18.003.757
4.87318
133.5381
Sí - Sí
0
0
0
Al-O
16.008,5345
4.79667
130.5659
Al-Al
0
0
0
P-O
9.034,2080
5.19098
19.8793
P-P
0
0
0
Una versión actualizada del potencial BKS introdujo un nuevo término repulsivo para evitar la superposición de átomos. [6] El potencial modificado se toma como
donde los parámetros constantes fueron elegidos para tener los siguientes valores para el vidrio de sílice:
Valor del parámetro para el vidrio de sílice
Sí - Sí
Si-O
O-O
D ij (eV•Å 24 )
3423200
29
113
Referencias
^ Buckingham, RA (1938). "La ecuación clásica de estado del helio, el neón y el argón gaseosos". Actas de la Royal Society A . 168 (933): 264–283. Bibcode :1938RSPSA.168..264B. doi :10.1098/rspa.1938.0173. JSTOR 97239.
^ ab F. Jensen, Introducción a la química computacional , 2.ª ed., Wiley, 2007,
^ Mason, Edward A. (29 de diciembre de 2004). "Propiedades de transporte de gases que obedecen a un potencial Buckingham (Exp-Six) modificado". Revista de Física Química . 22 (2): 169–186. doi :10.1063/1.1740026. ISSN 0021-9606.
^ van Beest, BWH; Kramer, GJ; van Santen, RA (16 de abril de 1990). "Campos de fuerza para sílices y aluminofosfatos basados en cálculos ab initio". Physical Review Letters . 64 (16): 1955–1958. Bibcode :1990PhRvL..64.1955V. doi :10.1103/physrevlett.64.1955. ISSN 0031-9007. PMID 10041537.
^ ab Kramer, GJ; Farragher, NP; van Beest, BWH; van Santen, RA (1991-02-15). "Campos de fuerza interatómicos para sílices, aluminofosfatos y zeolitas: Derivación basada en cálculos ab initio". Physical Review B . 43 (6): 5068–5080. Bibcode :1991PhRvB..43.5068K. doi :10.1103/physrevb.43.5068. ISSN 0163-1829. PMID 9997885.
^ Carré, Antoine; Ispas, Simona; Horbach, Jürgen; Kob, Walter (1 de noviembre de 2016). "Desarrollo de potenciales empíricos a partir de simulaciones ab initio: el caso de la sílice amorfa". Ciencia de materiales computacionales . 124 : 323–334. doi :10.1016/j.commatsci.2016.07.041. ISSN 0927-0256.