|
Métodos de Montecarlo |
La clase que describe el sistema físico
Los materiales dieléctricos están formados por dipolos eléctricos. Un dipolo eléctrico, véase la figura, es un sistema formado por dos cargas iguales y de signo contrario, separadas una distancia d. Se define el momento dipolar, como un vector cuyo módulo es el producto de la carga por la separación entre las mismas, cuya dirección es la recta que las une, y cuyo sentido es la flecha que va de la carga negativa a la positiva.
Si colocamos el dipolo en un campo eléctrico E, actúan sobre él un par de fuerzas, cuyo momento tiende a orientar al dipolo en la dirección del campo. Sin embargo, esta tendencia estará contrarrestada por la agitación térmica de las moléculas. Para cada campo y temperatura tendremos una orientación media resultado del compromiso entre ambas tendencias contrapuestas. La energía del dipolo en el campo eléctrico es
Como el ángulo q puede tener cualquier valor entre 0 y p, las energías de los estados disponibles no son discretas sino que varían de forma continua. La probabilidad P(q) de que el vector momento dipolar p apunte según un ángulo comprendido entre q y q +dq es proporcional a
La constante de proporcionalidad se determina a partir de la condición de que los dipolos con dicha energía apuntan con sus vectores momento dipolar en todas las direcciones del espacio cuyos ángulos están comprendidos entre q y q +dq. En la última figura, se representa el ángulo sólido dW correspondiente a dicha apertura angular.
La polarización de la sustancia dieléctrica será
Siendo N el número de moléculas y el resto de la fórmula, la componente media del momento dipolar en la dirección del campo. La integración conduce a la conocida ecuación de Langevin
En el caso límite u<<<1, es decir, para campos eléctricos débiles o temperaturas elevadas dicha función se puede aproximar a la recta
la polarización es una función lineal del cociente campo/temperatura, esta fórmula es análoga a la expresión de la ley de Curie para los materiales paramagnéticos.
Cuando u es grande la ecuación de Langevin tiende hacia un valor límite constante. Para campos eléctricos intensos o bajas temperaturas los dipolos están alineados, es decir, predomina el efecto orientador del campo.
Para aplicar el método de Montecarlo hemos de transformar la variable aleatoria angular continua q, distribuida con probabilidad P(q) en el intervalo [0, p) en una variable aleatoria uniformemente distribuida g en el intervalo [0, 1), véase el apartado variable aleatoria continua. La fórmula de transformación es
Despejamos de la fórmula anterior cos(q).
La clase Dielectrico es similar a la clase Paramagnetico descrita en la sección anterior, la única diferencia estriba en que hemos de transformar una variable continua en vez de una variable discreta, por lo que no tenemos que usar un array para guardar probabilidades sino tan solo aplicar la fórmula final del apartado anterior.
Los miembros dato de la clase son el momento dipolar po de las moléculas del material dieléctrico polar, y el número de partículas nParticulas que forman la sustancia dieléctrica. Se pasan al constructor estos dos datos, para crear un objeto de la clase Dielectrico que describa una muestra de dicha sustancia.
public class Dielectrico {
private double po; //momento dipolar
public int nParticulas;
//números aleatorios
private Random rnd=new Random();
private static final double COEFICIENTE=0.7246;
public Dielectrico(int nParticulas, double momento) {
this.nParticulas=nParticulas;
po=momento;
}
//funciones miembro ...
}
|
La constante COEFICIENTE, es el producto de 10-30 (unidad del momento dipolar) por 107 (unidad del campo eléctrico) dividido por la constante de Boltzmann 1.38 10-23.
La función miembro sorteo devuelve el cos(q) de cada molécula particular. Dicho valor se calcula mediante la última fórmula del apartado anterior.
private double sorteo(double campo, double temperatura){
double aleatorio=rnd.nextDouble();
double u=COEFICIENTE*po*campo/temperatura;
double coseno=Math.log((1-aleatorio)*Math.exp(u)+aleatorio*Math.exp(-u))/u;
return coseno;
}
|
La función momentoDipolarMedio calcula y devuelve el momento dipolar medio de una muestra significativa de moléculas dipolares para cada campo eléctrico y temperatura absoluta. Se suma todos los cosenos devueltos por la función sorteo, los multiplica por el momento dipolar po de las moléculas y los divide entre el número de moléculas. Como caso particular, consideramos aquél en el que la temperatura es cero, todos los momentos dipolares apuntan en la dirección del campo.
public double momentoDipolarMedio(double campo, double temperatura){
if(temperatura<=0.0){
return po;
}
double media=0.0;
for(int i=0; i<nParticulas; i++){
media+=sorteo(campo, temperatura);
}
return (media*po)/nParticulas;
}
|
El uso de la clase Dielectrico es completamente similar al uso de la clase Paramagnetico. Fijaremos el momento dipolar de las moléculas de la sustancia y la intensidad del campo eléctrico que actúa sobre la muestra de material dieléctrico y variaremos la temperatura.
Creamos un objeto di de la clase Dielectrico y le pasamos en su constructor el número de moléculas de la muestra y el tipo de material dieléctrico por ejemplo, el agua que tiene un momento dipolar de 6.2 (x 10-30 C m). Fijamos el campo eléctrico cuya intensidad vale 40 (x 107 N/C). El objeto di llama a la función miembro momentoDipolarMedio y le pasa la intensidad del campo y la temperatura absoluta, devolviendo el momento dipolar medio de la muestra.
Dielectrico di=new Dielectrico(1000, 6.2);
double campo=40.0; //en unidades 10 elevado a 7 N/C
System.out.println("temperatura \t m. dipolar medio");
for(double T=0.0; T<100; T+=10.0){
System.out.println(" "+T+ " \t \t"+
Math.floor(100*di.momentoDipolarMedio(campo, T))/100);
}
En la siguiente tabla se proporcionan los momentos dipolares de algunas sustancias dieléctricas
| Sustancia | Momento dipolar (x 10-30 Cm) |
| Agua | 6.2 |
| Nitrobenceno | 13.2 |
| Fenol | 5.2 |
| Bromhídrico | 2.6 |
| Clorhídrico | 3.5 |
| Iodhídrico | 1.3 |
· Completar tablas en las que fijada la intensidad del campo eléctrico y el tipo de material, se obtenga el momento dipolar medio en función de la temperatura.
· Que la polarización es una función lineal del campo temperatura, para altas temperaturas o baja intensidad del campo eléctrico (un comportamiento similar a la ley de Curie)
