https://doi.org/10.22463/2011642X.3122
Recibido: 12 de junio de 2018 - Aprobado: 6 de noviembre de 2018
Como citar:
J.J. Barba-Ortega, M. Rincón-Joya & J.A. Faúndez-Chaura, “Curva voltaje-tiempo en un proceso aniquilación-creación de pares vórtice-antivórtice”, Revista Ingenio, vol. 15(1), pp. 31-37, 2018
En esta contribución se estudió la respuesta eléctrica de la creación-aniquilación de vórtices cinemáticos en una película superconductora en presencia de una corriente continua externa Ja a campo magnético nulo. La película se encuentra en contacto con diferentes tipos de materiales simulados a través del parámetro de extrapolación de deGennes b; a su vez, posee una pequeña región metálica localizada a una distancia D del centro. Se usó la teoría de Ginzburg-Landau dependiente del tiempo para describir la variación voltaje en función del tiempo φ(t) en un proceso de creación-aniquilación de vórtices cinemáticos considerando varias posiciones del defecto metálico y del parámetro de deGennes. Nuestra investigación muestra que las corrientes críticas, amplitud y frecuencia de las oscilaciones en la curva φ(t) dependen fuertemente del parámetro de deGennes y de la posición del defecto en la muestra.
Palabras clave:Agua, Acueducto, Comunicación, Control, Optimización
In this contribution we study the electrical response of the creation-annihilation kinematics vortexs in a superconducting film in the presence of an external dc current Ja and a null magnetic field. The film is in contact with different types of materials simulated via the deGennes extrapolation parameter b; also it has a small metallic region located at a distance D from its center. We use the time-dependent Ginzburg-Landau theory to describe the voltage variation as a function of time φ(t) for a kinematic vortex creation-annihilation process considering various positions of the metallic defect and deGnnes parameter. Our research shows that the critical currents, amplitude and frequency of the oscillations in the φ(t) curve strongly depend on the deGennes parameter and the position of the defect in the sample.
Keywords: DeGennes, Ginzburg-Landau, Mesoscopic, Superconductor, Vortices.
En consecuencia, trabajos experimentales implementan la técnica de corriente pulsada (del orden de los nanosegundos) en películas superconductoras de niobio (Nb), reportan “phase-slip-center” estables en el tiempo cerca de la temperatura crítica y la formación aparente y preferente de “phase-slip-center” antes de la formación de los vórtices Abrikosov [10-11]. Se ha demostrado que la corriente pulsada tiene una serie de ventajas sobre las técnicas convencionales, como las mediciones de voltaje con múltiples sondas [12], imágenes láser [7] y sincronización por radiofrecuencia [13]
Por consiguiente, varios trabajos teóricos se han desarrollado analizando el estado resistivo de películas superconductoras en presencia de corrientes, por ejemplo, Sardella y colaboradores encontraron vórtices cinemáticos en un superconductor mesoscópico en presencia de un campo magnético externo; ellos describieron el pico de la curva característica de corriente-resistencia teniendo en cuenta la influencia del auto-campo magnético inducido por las corrientes y encontraron que las tasas de aniquilación de la super-corriente son altamente dependiente con la velocidad de aniquilación de los vórtices cinemáticos [14-16]. Berdiyorov y colaboradores analizaron el estado de vórtice en una película con centros de anclaje y puntos magnéticos; ellos mostraron que el estado resistivo depende de la corriente externa, el numero de defectos y la intensidad de los dipolos en los puntos magnéticos [17-19].
El efecto del parámetro deGennes sobre la respuesta electromagnética en muestras bidimensionales y tridimensionales lisas y rugosas fue estudiado por Barba-Ortega y colaboradores; ellos encontraron que las propiedades magnéticas del superconductor se modifican cuando los efectos de de-magnetización se tienen en cuenta; a su vez, se observó que la rugosidad de la muestra influye en los campos críticos y en la configuración de vórtices [20-22].
En este trabajo de investigación se estudió la respuesta magnética de una película superconductora con un defecto metálico no centrado en su interior, en presencia de una corriente continua (cc) aplicada y en ausencia de campo magnético, usando la teoría de Ginzburg-Landau dependiente del tiempo. Se encontró que el estado resistivo de las muestras superconductoras se caracteriza por la presencia de vórtices cinemáticos; además, analizamos los “phase-slip-centers” (es decir, vórtices de movimiento rápido) los cuales se caracterizan por un voltaje saturado de salida. El voltaje de respuesta del condensado superconductor al pulso de corriente aplicada se estudia en función del tiempo, las condiciones de contorno y la posición del defeco metálico; por consiguiente, se encontró que el valor de los tiempos en los cuales se da la primera entrada de vórtice cinemático a la muestra y en que ocurre la transición superconductor-normal dependen fuertemente de las condiciones de contorno y de la posición del defecto metálico.
De acuerdo con lo anterior, en la sección 2 describimos el formalismo teórico utilizado para estudiar una película superconductora en presencia de una corriente aplicada; en la sección 3, presentamos los resultados que surgen de la solución numérica de las ecuaciones Ginzburg-Landau para la curva voltaje-tiempo φ(t) a diferentes condiciones de contorno y posición del defecto.
2. Formalismo TeóricoEl sistema estudiado es tratado como un problema bidimensional porque el espesor s de la muestra es mucho menor que la longitud de coherencia ξ y la profundidad de penetración λ; además, despreciamos el apantallamiento del campo magnético, lo cual es válido para las muestras con un ancho efectivo ε<<2λ/s, donde s es el ancho real de la muestra o para muestras infinitas en la dirección z.
El sistema de ecuaciones Ginzburg-Landau, dependiente del tiempo generalizado, para el parámetro de orden ψ, el potencial vectorial A y el potencial escalar φ, en forma adimensional está dado por la Ecuación (1) y Ecuación (2) [22-27].
donde se toma g=1 para simular todas las regiones completamente superconductoras, mientras que g=0 simula superficies con la superconductividad completamente depreciada.
Las ecuaciones Ginzburg-Landau (Ecuación (1) y Ecuación (2)), fueron escalada así: ψ en unidades de ψ∞, las distancias se escalan por la longitud de coherencia ξ, el tiempo t está en unidades del tiempo de Ginzburg-Landau t0=πℏ/8KB Tu~10-12 s, el potencial electrostático φ es dado en unidades de φ0=ℏ/2et0, y el potencial vectorial A se escala por Hc2 ξ (Hc2 es el campo crítico superior en el volumen); a partir de los primeros principios, obtenemos los parámetros u=5.79 y Γ=10 [24].
Las condiciones de contorno de Neumann se seleccionan en los límites de la muestra, excepto en los electrodos, donde e usó ψ=0, y φ.n=-J (J se expresa en unidades de J0=hc2/16eπ2 λ2 ξ). La condición de contorno general para el parámetro de orden, que representa la imposibilidad del transporte de pares de Cooper a través de las componentes normales de las interfaces esta definida por la Ecuación (3).
donde n es el versor perpendicular a las interfaces, b es la longitud de extrapolación de deGennes y simula la condición de contorno [21]; si b→∞ reproduce una interfaz superconductor-vacío; si b>0, simula un superconductor en contacto con un metal; y si b<0, simula una frontera superconductor-superconductor a mayor temperatura crítica en estado normal. El voltaje es calculado por medio de la Ecuación (4).
En la Figura (2) graficamos el voltaje φ en función del tiempo característico para la película rodeada por diferentes tipos de materiales, consideramos interfaces (a) superconductor-vacío, b→∞, (b) superconductor-metal, b=0.5ξ, (c) superconductor-ferromagneto, b=0, y (d) superconductor–superconductor a mayor temperatura critica, b=-0.5ξ, con el defecto ubicado en D1=4.2ξ.
En la Figura (2) se observan dos tiempos críticos t1 y t2 que identifican el tiempo que le lleva al sistema permitir el ingreso del primer par vórtice-antivórtice y en ocurrir la transición al estado normal, respectivamente; como se aprecia, en la Figura (2), al considerar una interfaz con un superconductor a mayor temperatura crítica la súper-corriente y la barrera de energía aumentan, aumentando los tiempos críticos de transición y de ingreso del primer par de vórtices. Por consiguiente, se encuentra t2max≈15 ns, t1max≈6.5 ns para b=-0.5ξ y un t2min≈12 ns, t1max≈4.5 ns para b=0.
En las Figura (3) y Figura (4) se presenta el comportamiento del voltaje φ en función del tiempo característico t0 para b→∞, considerando al defecto ubicado a una distancia D1=4.2ξ, D2=4.2ξ, D3=2.8ξ, y D4=5.6ξ, del centro de la muestra. El voltaje de salida oscila en el tiempo con máximos bien definidos; cada máximo en la curva φ(t) corresponde a la penetración o aniquilación de un par de vórtice-antivórtice y ocurre de manera bien definida de cuantos de grupo de ondas. Además, encontramos los valores de la amplitud media y la frecuencia del grupo de oscilaciones de los casos.D3=2.8ξ, y D4=5.6ξ; llamamos Δti al tiempo medio de vida del i-ésimo paquete de ondas.
De acuerdo con lo anterior, se obtiene Δt1=625 ns, Δt2=375 ns, Δt3=415 ns, Δt4=175 ns y después de este tiempo, los paquetes de onda alcanzan un valor igual Δti>4=175 ns para el caso D3=2.8ξ. Asimismo, para el caso D4=5.6ξ, conseguimos Δt1=315 ns, Δt2=285 ns, y después de este tiempo, los paquetes de onda obtienen un valor igual Δti>2=285 ns; por consiguiente, existe una clara incidencia de la posición del defecto sobre el tiempo de vida de cada paquete de ondas. Analizando más a detalle la naturaleza de las ondas podemos calcular la amplitud no del paquete de ondas si no, de una única onda.
Para los casos D3 y D4, el tiempo de viaje y la amplitud de la onda de aniquilación par vórtice-antivórtice es ΔT y ΔV/2 respectivamente, ΔT=8 ns es el mismo para los casos D4 y D3, pero con un máximo de amplitud 2ΔV=1.05 para D3 y 2ΔV=0.88 para D4. Este resultado muestra la incidencia de la presencia del defecto sobre el tiempo de vida del proceso de aniquilación de los vórtices cinemáticos.
En las Figura (5) y Figura (6) se muestra el comportamiento de la curva de voltaje φ(t) para b=10ξ, y b=5ξ, en varios intervalos temporales.
En consecuencia, en los casos estudiados anteriormente, se escoge una posición constante del defecto D1=4.2ξ y variamos las condiciones de contorno; entonces, se obtiene Δt1=200 ns, Δt2=250 ns, para b=5ξ y Δt1=185 ns, para b=10ξ. Por consiguiente, el tiempo de viaje de la onda del par de vórtices cinemático es Δt=6 ns, 2ΔV=0.52, para b=5ξ, y Δt=4 ns, 2ΔV=0.34 para b=10ξ; Por lo tanto, este resultado muestra la incidencia del material en contacto con la muestra sobre los tiempos de vida de estos procesos.
[1] B. I. Ivlev, N. B. Kopnin, “Electric currents and resistive states in thin superconductors”, Advances in Physics, vol. 33(1), pp. 47-114, June 2006. Doi: https://doi.org/10.1080/00018738400101641
[2] W. W. Webb, R. J. Warburton, “Intrinsic Quantum Fluctuations in Uniform Filamentary Superconductors”, Physical Review Letters, vol.20(9), pp. 461-465, Febraury 1968. Doi: https://doi.org/10.1103/PhysRevLett.20.461
[3] J. D. Meyer, “Voltage steps in the U(T)-transition curves and in the U(I)-characteristic curves of current-carrying whiskers of tin”, Applied Physics, vol. 2, pp. 303-320, December 1973.
[4] W. J. Skocpol, M. R. Beasley, M. Tinkham, “Phase-slip centers and nonequilibrium processes in superconducting tin microbridges”, Journal of Low Temperature Physics, vol. 16, pp. 145-167, July 1974. Doi: https://doi.org/10.1007/BF00655865
[5] A. Andronov, I. Gordion, V. Kurin, I. Nefedov, I. Shereshevsky, “Kinematic vortices and phase slip lines in the dynamics of the resistive state of narrow superconductive thin film channels”, Physica C, vol. 213, pp. 193-199, August 1993. Doi: https://doi.org/10.1016/0921-4534(93)90777-N
[6] G. R. Berdiyorov, M. V. Milosevic, F. M. Peeters, “Kinematic vortex-antivortex lines in strongly driven superconducting stripes”, Physical Review B, vol. 79(184506), pp.1-8, May 2009. Doi: 10.1103/PhysRevB.79.184506
[7] A. G. Sivakov, A. M. Glukhov, A. N. Omelyanchouk, Y. Koval, P. Muller, A. V. Ustinov, “Josephson Behavior of Phase-Slip Lines in Wide Superconducting Strips”, Physical Review Letters, vol. 91(267001), pp. 1-4, December 2003.
[8] J. Barba-Ortega, E. D. V. Niño, E. Sardella, “Fenómeno de transporte en superconductores: vórtices cinemáticos”, Iteckne, vol. 14(1), pp. 11-16, June 2017
[9] J. Barba-Ortega, R. Zadorosny, E. Sardella, “Influence of the deGennes extrapolation parameter on the resistive state of a superconducting strip”, Physics Letters A, vol. 382, pp. 215-219, January 2018.Doi: https://doi.org/10.1016/j.physleta.2017.11.010
[10] F. R. Ladan, K. h. Harrabi, M. Rosticher, C. Villard, P. Mathieu, J. P. Maneval, “Current-temperature diagram of resistive states in long superconducting niobium filaments”, Journal of Low Temperature Physics, vol. 153, pp. 103-122, October 2008. Doi: 10.1007/s10909-008-9821-7
[11] K. Harrabi, F. R. Ladan, Vu Dinh Lam, J. P. Maneval, J. F. Hamet, J.P. Villegier, R. W. Bland, “Current-temperature diagram of resistive states in long superconducting YBa2Cu3O7 strips”, Journal of Low Temperature Physics, vol. 157(XXX), pp. 36-56, September 2009. Doi: 10.1007/s10909-009-9930-y
[12] G. J. Dolan, L. D. Jackel, “Voltage measurements within the nonequilibrium region near phase-slip centers”, Physical Review Letters, vol. 39(25), pp. 1628-1632, December 1977. Doi: https://doi.org/10.1103/PhysRevLett.39.1628
[13] R. B. Laibowitz, A. N. Broers, J. T. C. Yeh,and J. M. Viggiano, “Josephson effect in Nb nanobridges”,Applied Physics Letters, vol. 35, pp. 891-899, August 2008.
[14] E. Sardella, P. N. Lisboa-Filho, C. C. S. Silva, L. R. E. Cabral and W. A. Ortiz, “Vortex-antivortex annihilation dynamics in a square mesoscopic superconducting cylinder”, Physical Review B, vol. 80(012506), pp. 1-4, July 2009.
[15] E. C. S. Duarte, E. Sardella, W. A. Ortiz, R. Zadorosny, “Dynamics and heat diffusion of Abrikosov’s vortex-antivortex pairs during an annihilation process”, Journal of Physics Condensed Matter, vol. 29(405605), pp. 1-9 July 2017.
[16] E. C. S. Duarte, E. Sardella, W. A. Ortiz, R. Zadorosny, “Dynamics and heat diffusion of Abrikosov’s vortex-antivortex pairs during an annihilation process” Journal of Physics: Condensed Matter, vol. 29(40), pp. 1-16, September 2017.
[17] G. Berdiyorov, K. Harrabi, F. Oktasendra, K. Gasmi, A. I. Mansour, J. P. Maneval, F. M. Peeters, “Dynamics of current-driven phase-slip centers in superconducting strips”, Physical Review B, vol. 90(054506), pp. 1-7, August 2014. Doi: https://doi.org/10.1103/PhysRevB.90.054506
[18] G. Berdiyorov, K. Harrabi, J. P. Maneval, F. M. Peeters, “Effect of pinning on the response of superconducting strips to an external pulsed current”, Superconductor Science and Technology, vol. 28(025004), pp. 1-8, December 2015. Doi: 10.1088/0953-2048/28/2/025004
[19] G. R. Berdiyorov, X. H. Chao, F. M. Peeters, H. B. Wang, V. V. Moshchalkov, B. Y. Zhu “Magnetoresistance oscillations in superconducting strips: a Ginzburg-Landau study”, Physical Review B, vol. 86(224504), pp. 1-8, December 2012. Doi: https://doi.org/10.1103/PhysRevB.86.224504
[20] J. Barba-Ortega, E. Sardella, J. A. Aguiar, “Superconducting boundary conditions for mesoscopic circular samples”, Superconductor Science and Technology, vol. 24(015001), pp. 1-7, December 2010. Doi: 10.1088/0953-2048/24/1/015001
[21] J. Barba-Ortega, E. Sardella, “Superconducting properties of a mesoscopic parallelepiped with anisotropic surface conditions”, Physics Letters A, vol. 379(XXX), pp. 3130-3135, October 2015. Doi: 10.1016/j.physleta.2015.10.013
[22] J. Barba-Ortega, E. Sardella, J. A. Aguiar, E. H. Brandt, “Vortex state in a mesoscopic flat disk with rough surface”, Physica C, vol. 479(XXX), pp. 49-52, Febraury 2012.Doi: 10.1016/j.physc.2012.01.021
[23] Q. Du, M. D. Gunzburger, “A model for superconducting thin films having variable thickness”, Physica D, vol. 69, no. XXX, pp. 215-231, December 1993.
[24] L. Kramer, R. J. Watts-Tobin, “Theory of dissipative current-carrying states in superconducting filaments”, Physical Review Letters, vol. 40(XXX), pp. 1041-1044, April 1978.
[25] R. J. Watts-Tobin, Y. Krähenbühl, L. Kramer “Nonequilibrium theory of dirty, current-carrying superconductors: phase-slip oscillators in narrow filaments near Tc”, Journal of Low Temperature Physics, vol. 45(69), pp. 459-501, January 1981.
[26] M. Silaev, “Microscopic derivation of two-component Ginzburg-Landau model and conditions of its applicability in two-band systems”, Physical Review B, vol. 85(134514), pp. 1-8, April 2012. Doi: https://doi.org/10.1103/PhysRevB.85.134514
[27] D. Y. Vodolazov, F. M. Peeters, M. Morelle, V. V. Moshchalkov, “Masking effect of heat dissipation on the current-voltage characteristics of a mesoscopic superconducting sample with leads”, Physical Review B, vol. 71(184502), pp. 1-8, May 2005.Doi: https://doi.org/10.1103/PhysRevB.71.184502
* Doctor.Correo: jjbarbao@unal.edu.co
** Doctor.Correo: mrinconj@unal.edu.co
*** Magister. Correo: julian@fisica.ufmt.br