Recibido: 24 de noviembre de 2019
Aceptado: 21 de Abril de 2020
En este trabajo calculamos la dependencia de la sensibilidad de cinco células cancerígenas con el ángulo incidente de la radiación en un cristal fotónico unidimensional. El estudio se realizó para dos tipos de polarización: Transversal Eléctrica (TE) y Transversal Magnética (TM). Las muestras de las células se infiltran mediante biopsia líquida dentro de una nano cavidad que funciona como recipiente, y es revestida por nanocompuestos formados por la mezcla de la sílice con nano partículas de plata de tal forma que rompe la periodicidad del cristal fotónico compuesto por capas de SiO2 y aire. Encontramos picos de resonancia en el espectro de transmisión que caracteriza a cada célula infiltrada, por ende, este mecanismo de detección permite que el cristal fotónico funcione como un biosensor. Mediante el método de la matriz de trasferencia, para la polarización TE se observó un incremento de la sensibilidad al aumentar el ángulo incidente. Sin embargo, en la polarización TM no existieron aumentos significativos. Adicionalmente, se varía el espesor de la nano cavidad y el factor de relleno para obtener una mayor optimización. Los resultados revelan un incremento en la sensibilidad al aumentar el espesor de la nanocavidad, mientras que al aumentar el factor de relleno la sensibilidad decrece.
Palabras clave: cristal fotónico, método de la matriz de transferencia, célula cancerígena, banda fotónica prohibida, biosensor.
In geotechnical engineering, bored-pile wall stability is evaluated using deterministic design methods based on safety factors to establish a margin against failure. In recent years, reliability-based design methods have been adopted to include uncertainty in the assessment of bored-pile wall stability as well as in the calculation of the feasible embedment depth of the walls. In this study, an expanded reliability-based design approach, along with finite element analysis, was applied to conduct parametric analyses of bored-pile wall stability. In serviceability limit state design framework, the results indicate that cohesion and groundwater level are factors that significantly affect bored-pile wall stability. Moreover, high variability in the cohesion range causes great uncertainty to determine the embedment depth of bored-pile wall. The feasible embedment depth can reach 4 times the free height considering the maximum coefficient of variation (50 %) of the cohesion. In turn, when the groundwater level is located at the retained ground surface, the horizontal displacement of the upper end of the wall reaches 15.2 mm, i.e., 0.0038 times the free height of the wall, for which the soil mobilizes active earth pressures. It was also found that the resolution of probabilistic results is highly influenced by the number of iterations in Monte Carlo simulations.
Keywords: Photonic crystal, transfer matrix method, cancer cell, photonic bandgap, biosensor.
Los cristales fotónicos (CF) son estructuras periódicas compuestas por capas alternadas de distinto nanomaterial, que afectan la propagación de la luz, dando origen a características ópticas únicas como las bandas fotónicas prohibidas (BFP). La luz dentro del BFP en determinados rangos de longitudes onda no se le tiene permitido propagarse, debido a la interferencia destructiva en las interfaces de las capas [
Los cristales fotónicos unidimensionales (CF-1D), también conocidos como espejos de Bragg, han recibido considerable atención debido a su alta reflectividad, simplicidad en la fabricación en comparación con cristales fotónicos 2D y 3D, así como también a una gran cantidad de aplicaciones especialmente como reflectores [
Recientemente, Francis et all., propusieron un CF-1D compuestos de capas alternadas de sílice fundida y agua marina, para la simultanea detección de salinidad y temperatura de determinadas muestras de agua. Ellos encontraron que, a medida que la temperatura aumentaba dejando la salinidad constante, el modo defectivo se corre a longitudes de ondas más cortas, y al mantener constante la temperatura aumentando la salinidad, el corrimiento se produce a longitudes de onda más largas [
El cáncer es actualmente una de las enfermedades con una alta tasa de mortalidad, debido a que la mayoría de los casos son diagnosticados en su etapa más avanzada (metástasis) [
Normalmente las células cancerígenas se han caracterizado por tener un índice de refracción mayor que las células normales que oscilan en un rango de 1.35-1.37, la gran acumulación de proteína en el núcleo es lo que origina la división incontrolable de estas células que destruyen los tejidos del cuerpo [
Ramanujan et al. Propusieron un CF-1D para la detección de células cancerígenas, introduciendo una cavidad revestida de nanocompuestos con incidencia normal de la radiación. Ellos optimizaron el factor de relleno de las nanopartículas y el espesor de los nanocompuestos obteniendo una sensibilidad máxima de 43 nm/RIU [
En este trabajo optimizaremos la sensibilidad del dispositivo para la detección de cinco células cancerígenas diferentes, variando el ángulo de incidencia de la radiación y el espesor de la nano cavidad para la polarización TE y TM. El método de la matriz de transferencia (MMT) será implementado para el estudio de la propagación de esta radiación.
En la Fig.1 se muestra la estructura del cristal fotónico unidimensional compuesto por la alternación de los materiales SiO2 (capa A) y aire (capa B) de espesores d1 y d2, respectivamente, con periodicidad N en el eje x. La simetría de translación se rompe al introducir la nano cavidad en el cual las muestras de las células son infiltradas (capa D) revestido por el nano compuesto de Ag/SiO2 (capa C) que se forma al mezclar las nanopartículas de plata con la sílice. El espesor para estos últimos materiales es denotado por d3 y d4, respectivamente.
Por medio del método de la matriz de transferencia [
Se asume que los materiales del CF son isotrópicos y no magnetizados.
Para la polarización TE, la matriz para la j-ésima capa (=A, B, C y D) de espesores dj (d1, d2, d3 y d4) está representada por (1):
Mj=Dj Pj Dj-1 (1)
Donde Pj es la matriz de propagación y está dada en (2):
La fase φj se expresa como (3):
En (3) εj es la constante dieléctrica, dj el espesor y θj el ángulo de incidencia para la j-ésima capa (4).
Dj es la matriz dinámica y está dada por (4):
Por lo tanto, la matriz de transferencia total para la estructura (AB)N(CDC)(AB)N puede ser calculada como se muestra en (5).
Donde Da es la matriz dinámica del aire que rodea al CF-1D defectivo,
La transmitancia es calculada con el elemento m11 de la matriz M como (6).
De forma similar, para obtener la matriz para la j-ésima capa con polarización TM en (1), la matriz de propagación se mantiene intacta y el cambio se hace únicamente en la matriz dinámica que se encuentra dada en (7).
Por consiguiente, se calcula la matriz de transferencia total y la transmitancia utilizando (5) y (6).
La constante dieléctrica efectiva εc(ω) del nanocompuesto es determinada usando (8) de Maxwell-Garnett [
Donde εd es la constante dieléctrica de la sílice, εm(ω) es la constante dieléctrica del metal en este caso plata, f es el factor de relleno de las nanopartículas y ω la frecuencia angular.
Ahora despejando (8) se obtiene εc(ω) explícitamente como se muestra en (9).
La constante dieléctrica de las nanopartículas de metal εm(ω) en este caso plata está dada por el modelo de Drude (10).
Donde εo=5 es la constante dieléctrica relativa para la plata, εp = 13.64 x 1015rad/s es la frecuencia del plasma de gases de electrones libres en un volumen infinito, y γ =3.03 x1013rad/s es la constante de amortiguamiento.
Finalmente, la constante dieléctrica de sílice εd se obtiene por medio de la ecuación de Sellmeier [
Donde B1=0.6961663, B2= 0.4079426, B3=0.8974794, C1=0.00467914826µm2, C2 =0.0135120631µm2, C3= 97.9340025 µm2 son los coeficientes de Sellmeier determinados experimentalmente.
Se considera que el medio de entrada y salida de la radiación es aire. Los espesores de las capas que componen el CF son: para la capa de sílice d1 = 268nm, capa de aire d2 = 388nm, capa del nano compuesto d3 = 90nm y el de la capa de la nano cavidad d4 = 90nm. Los espesores de las capas se eligieron un cuarto de longitud de onda cumpliendo con la relación nHdH=nL dL= λo/4 , donde λo es la longitud de onda con el cual se quiera trabajar. La periodicidad de las capas es N = 10 con un factor de relleno de las nano partículas de plata de 0.1.
El mecanismo de detección se debe a la aparición de picos de resonancia dentro de la banda fotónica prohibida que caracteriza a cada tipo de célula infiltrada por medio de biopsia líquida [
Tipo de célula | Constante dieléctrica |
Célula Normal | 1,8225 |
Jurkat | 1,932100 |
HeLa | 1,937660 |
PC12 | 1,946025 |
MDA-MB-231 | 1,957201 |
MCF-7 | 1,962801 |
En este trabajo se busca optimizar la sensibilidad del biosensor variando parámetros como el ángulo de incidencia, espesor de la nano cavidad y el factor de relleno para la polarización TE y TM.
Es decir, obtener la mayor respuesta óptica o separación de los picos de resonancia a pequeños cambios de índice de refracción.
La Fig. 2, muestra el espectro de transmitancia con incidencia normal de la radiación, al infiltrar distintos tipos de células en la nano cavidad. Se logró observar corrimientos de los picos de resonancia de las células cancerígenas hacia longitudes de onda más largas respecto a la célula normal. Este corrimiento se produce debido al incremento de la constante dieléctrica.
En la Fig. 3, se muestra los efectos de incidencia oblicua de la radiación con polarización TE, consiguiendo corrimientos de los picos de resonancia hacia longitudes de ondas más largas, cuando se cambiaba de célula. Sin embargo, a medida que se incrementaba el ángulo de incidencia, los picos de resonancia se corrían a longitudes de ondas más cortas, con un decrecimiento en la transmitancia; esto se debe al aumento del camino geométrico que debe recorrer el haz de luz. Además, la separación entre los picos de resonancia de las células cancerígenas respecto a la célula normal se hacía mayor; esto se ve reflejado en la Tabla 2 con el cálculo de la sensibilidad. Para el caso de polarización TM, en la Fig. 4, de forma similar, se obtiene corrimientos a longitudes de onda más largas al cambiar la célula.
Células Cancerígenas | Sensibilidad (nm/RIU) | |||
Ángulo 0° | Ángulo 25° | Ángulo 50° | Ángulo 75° | |
Jurkat | 43,5 | 44 | 46 | 49,5 |
HeLa | 43,57 | 44,04 | 46,95 | 49,52 |
PC12 | 43,55 | 44,22 | 46 | 49,55 |
MDA-MB-231 | 43,67 | 44,28 | 46,12 | 49,59 |
MCF-7 | 43,52 | 44,31 | 46,07 | 49,6 |
No obstante, el biosensor es mucho más sensible en el sentido de que los picos de resonancia se corren a longitudes de onda mucho más cortas que para el caso de la polarización TE. Adicionalmente, al aumentar el ángulo de incidencia la transmitancia aumenta y la separación de los picos de resonancia se mantiene casi constante. Esto se puede observar en la Tabla 3, cuando se calculó su sensibilidad en el cual solo existe un pequeño aumento a un ángulo de 50°.
Células Cancerígenas | Sensibilidad (nm/RIU) | |||
Ángulo 0° | Ángulo 25° | Ángulo 50° | Ángulo 75° | |
Jurkat | 43,5 | 43,5 | 43,75 | 43,5 |
HeLa | 43,57 | 43,57 | 43,81 | 43,57 |
PC12 | 43,55 | 43,55 | 43,77 | 43,55 |
MDA-MB-231 | 43,67 | 43,67 | 43,87 | 43,67 |
MCF-7 | 43,52 | 43,72 | 43,92 | 43,52 |
Ahora, se hace una optimización de la sensibilidad variando los espesores de la nano cavidad. Para ello se mantiene constante el ángulo de incidencia con el cual se obtuvo mayor sensibilidad, en este caso para polarización TE con ángulo de incidencia de 75° y para TM con ángulo de incidencia de 50°. La Fig. 5 muestra que para la polarización TE a medida que se incrementa el espesor de la nano cavidad los picos de resonancia se corren a longitudes de onda más largas para cada tipo de célula. Además, la separación de los picos entre la célula normal y las células cancerígenas aumenta. Así, obteniendo una mayor sensibilidad como se evidencia en la Tabla 4. De igual forma, ocurre el mismo efecto para la polarización TM donde existe un corrimiento de los picos de resonancia a longitudes de onda largas como se observa en la Fig. 6 y un aumento en la sensibilidad que se muestra en la Tabla 5.
Células Cancerígenas | Sensibilidad (nm/RIU) para distintos espesores | ||||
70 nm | 80 nm | 90 nm | 100 nm | 110 nm | |
Jurkat | 39,5 | 44,5 | 49,5 | 54,25 | 59 |
HeLa | 39,52 | 44,52 | 49,52 | 54,28 | 59,04 |
PC12 | 39,55 | 44,66 | 49,55 | 54,22 | 58,88 |
MDA-MB-231 | 40,61 | 44,69 | 49,59 | 54,28 | 59,97 |
MCF-7 | 39,6 | 44,71 | 49,61 | 54,31 | 59,01 |
Células Cancerígenas | Sensibilidad (nm/RIU) para distintos espesores | ||||
70 nm | 80 nm | 90 nm | 100 nm | 110 nm | |
Jurkat | 34,5 | 39,25 | 43,75 | 48 | 52,25 |
HeLa | 34,76 | 39,28 | 43,81 | 48,09 | 52,38 |
PC12 | 34,88 | 39,33 | 43,77 | 48,22 | 52,22 |
MDA-MB-231 | 34,89 | 39,38 | 43,87 | 48,16 | 52,24 |
MCF-7 | 34,9 | 39,41 | 43,92 | 48,23 | 50,19 |
Finalmente, se varía el factor de relleno de las nanopartículas de plata, manteniendo constante el espesor de la nano cavidad d4=90nm, y los ángulos de incidencia para polarización TE 75° y para la polarización TM 50°. Como se puede observar tanto en la Fig. 7 (a) para polarización TE como en la Fig. 7 (b) polarización TM para una célula normal, al aumentar el factor de relleno de las nanopartículas de plata, los picos de resonancia dentro del espectro de transmisión se acortan y hay corrimientos a longitudes de onda más largas. En el caso de infiltrar una célula cancerígena Jurkat ocurre el mismo fenómeno, como se muestra en la Fig. 8. El mismo efecto ocurre si se realiza la simulación para las células HeLa, PC-12, MDA-MB-231 y MCF-7.
La Tabla 6 muestra los valores de la sensibilidad para la polarización TE con ángulo de incidencia de 75° para la célula normal y las células cancerígenas. Se pudo notar que existe una reducción de la sensibilidad al incrementar el factor de relleno de las nanopartículas de plata.
Células Cancerígenas | Sensibilidad (nm/RIU) para distintos factores de relleno | ||||
0 | 0,1 | 0,2 | 0,3 | 0,4 | |
Jurkat | 62 | 49,15 | 37,25 | 26 | 16,5 |
HeLa | 61,90 | 49,52 | 37,14 | 25,95 | 16,66 |
PC12 | 62 | 49,55 | 37,11 | 26 | 16,66 |
MDA-MB-231 | 62,04 | 49,59 | 37,14 | 25,91 | 16,53 |
MCF-7 | 61,15 | 49,60 | 37,25 | 26,07 | 16,66 |
De forma similar, la Tabla 7 muestra el mismo comportamiento para la polarización TM con ángulo de incidencia de 50°. En ambas, la mayor sensibilidad se consigue a un factor de relleno f=0. Es decir, cuando la sílice se encuentra en su forma pura. Sin embargo, los analitos preservan suscaracterísticas, debido a la mezcla de las nanopartículas de plata con la sílice. Por consiguiente, la mayor sensibilidad del biosensor se consigue a un factor de relleno de f = 0,1.
Células Cancerígenas | Sensibilidad (nm/RIU) para distintos factores de relleno | ||||
0 | 0,1 | 0,2 | 0,3 | 0, 4 | |
Jurkat | 52,25 | 43,75 | 34,75 | 26,5 | 19 |
HeLa | 52,38 | 43,81 | 34,76 | 26,43 | 19,04 |
PC12 | 52,44 | 43,77 | 34,88 | 26,44 | 19,1 |
MDA-MB-231 | 52,45 | 43,87 | 34,89 | 26,53 | 18,97 |
MCF-7 | 52,35 | 43,92 | 34,90 | 26,47 | 19,01 |
Los corrimientos de los picos de resonancia ocurren debido al cumplimiento de la siguiente relación para ondas estacionarias (14) [
δ =Κλ=n Δd (14)
Donde δ es la diferencia de camino óptico, K es un entero, λ es la longitud de onda de la radiación incidente, n el índice de refracción de cada célula y Δd es la diferencia de camino geométrico (espesor de la nanocavidad). Por tanto, al incrementar n debido al cambio de célula y manteniendo constante Δd y K, debe existir un incremento de λ.
Mediante el método de la matriz de transferencia se calculó el espectro de transmisión en un cristal fotónico unidimensional con la presencia de las células cancerígenas Jurkat, HeLa, PC-12, MDA-MB-231 y MCF-7. Se logró obtener la aparición de picos de resonancia que caracteriza a cada una de estas células, en el cual se corrían a longitudes de onda más largas cuando la constante dieléctrica de las células aumentaba. Se observó que para la polarización TE, al aumentar el ángulo de incidencia de la radiación, se obtuvo corrimientos de longitudes de ondas más cortas acompañado con el decrecimiento de la transmitancia y un aumento en la sensibilidad. Por otro lado, para la polarización TM se obtuvo corrimientos de longitudes de onda más cortas que en la polarización TE con un incremento en la transmitancia. Sin embargo, la sensibilidad se mantiene constante, excepto para el ángulo de 50° que tiene un pequeño aumento.
Finalmente, para una mayor optimización de la sensibilidad se variaron parámetros como: el espesor de la nano cavidad y factor de relleno de las nanopartículas de plata. Al aumentar el espesor de la nano cavidad encontramos corrimientos de los picos de resonancia a longitudes de más largas con un incremento en la sensibilidad, mientras que al aumentar el factor de relleno también se apreciaba corrimientos de los picos de resonancia a longitudes de onda más largas. Sin embargo, la sensibilidad decrecía. De esta forma, se logró obtener las condiciones necesarias para que la separación entre los picos de resonancia que caracterizan a cada célula cancerígena sea más distante respecto a la de una célula normal, de tal forma que permita una detección mucho más clara.
A la Universidad Surcolombiana por su apoyo financiero del presente proyecto de investigación en modalidad de Trabajo de grado.
1 Se encargó de la conceptualización, metodología, cálculos numéricos, análisis de la investigación y escritura del manuscrito.
2 Se encargó de la conceptualización, recursos de la investigación, supervisión y validación de los resultados, y de la revisión final del manuscrito.