\documentclass[a4paper, 11pt]{article}
\usepackage[utf8]{inputenc}

\usepackage[
	a4paper,
	headheight = 20pt,
	margin = 1in,
	tmargin = \dimexpr 1in - 10pt \relax
]{geometry}

\usepackage{fancyhdr} % for headers and footers
\usepackage{graphicx} % for including figures
\usepackage{booktabs} % for professional tables
\setlength{\headheight}{14pt}


\fancypagestyle{plain}{
	\fancyhf{}
	\fancyhead[L]{\sffamily Radboud University Nijmegen}
	\fancyhead[R]{\sffamily Superconductivity (NWI-NM117), Q3+Q4}
	\fancyfoot[R]{\sffamily\bfseries\thepage}
	\renewcommand{\headrulewidth}{0.5pt}
	\renewcommand{\footrulewidth}{0.5pt}
}
\pagestyle{fancy}

\usepackage{siunitx}
\usepackage{hyperref}
\usepackage{float}
\usepackage{mathtools}
\usepackage{amsmath}
\usepackage{todonotes}
\setuptodonotes{inline}
\usepackage{mhchem}

\newcommand{\pfrac}[2]{\frac{\partial #1}{\partial #2}}

\title{Superconductivity - Assignment 3}
\author{
	Kees van Kempen (s4853628)\\
	\texttt{k.vankempen@student.science.ru.nl}
}
\AtBeginDocument{\maketitle}

% Start from 8
\setcounter{section}{7}

\begin{document}

\section{\ce{Nb3Sn} cylinder}
Consider a cylinder of \ce{Nb3Sb}.
From lecture 4, we have the following properties for \ce{Nb3Sn}:
$T_c = \SI{18.2}{\kelvin}$,
$\xi = \SI{3.6}{\nano\meter}$,
$\lambda = \SI{124}{\nano\meter}$,
$\kappa = \frac{\lambda}{\xi} = 34 > \frac{1}{\sqrt{2}}$,
which means we are indeed dealing with a type-II superconductor.
As $B_{c1} < B_E < B_{c2}$, the cylinder is in the vortex state.
From the previous set of assignments, we know what the currents in the cylinder look like.
From free energy considerations, we have found in lecture 4 that for type-II superconductors, it is favorable to allow flux quanta inside the superconductor in this vortex state.
In this derivation, the contribution of one flux quantum is considered, but the consideration holds for many vortices, until they start to interact and repel eachother.
At that point, the vortex-vortex interaction orders the vortices in a lattice.
When the vortex cores start to overlap, there are no superconducting regions left, thus the material enters the normal conducting state.\footnote{I wanted to paint a complete picture althought it is not needed to answer the question.}
Minimizing the free energy over the flux shows the energy is lowered for determined thresholds $B_{c1} < B_E < B_{c2}$.

Let's start with the result from said free energy considerations.
The average field inside the cylinder is given by the following self-consisting equation as
\[
	B = B_E - \frac{\phi_0}{8\pi\lambda^2}\ln{\frac{\phi_0}{4\exp{(1)}\xi^2B}}.
\]
Plugging in the values for \ce{Nb3Sn}, $B_E = \SI{1}{\tesla}$, and $\phi_0 = \SI{2.0678}{\weber}$, $B$ is found as $B = \SI{.986}{\tesla} \approx B_E$ by intersection.
%https://www.wolframalpha.com/input?i=B+%3D+1+-+%282.0678*10%5E-15%29%2F%288*pi*%28124*10%5E-9%29%5E2%29+*+ln%282.0678*10%5E-15%2F%284*exp%281%29*%283.65*10%5E-9%29%5E2+*+B%29%29
This is in the range as provided in the assignment ($B = \SI[separate-uncertainty]{.981\pm.019}{\tesla}$).

To investigate the inhomogenity of the field inside the cylinder, we look at the gradient $\nabla B$ inside the material.
As we assume a vortex lattice that fully fills the cross section of the cylinder,
and we assume that the fields due to each vortex die out quickly enough to not overlap,
it suffices to calculate the gradient over just one vortex.
These assumptions coincide with slide 15 of lecture 4, from which I took figure \ref{fig:lec4-vortexlattice}.

\begin{figure}[H]
	\centering
	\label{fig:lec4-vortexlattice}
	\includegraphics[width=.4\textwidth]{lec4-vortexlattice.png}
	\caption{The vortices are arranged in a lattice to maximize their distance, as this lowers their repulsive interaction and thus the energy.}
\end{figure}

On slide 19 from the same week, we find an expression $B(r)$ for the field at distance $r$ from the vortex core as
\[
	B = \frac{\phi_0}{2\pi\lambda^2} K_0(r/\lambda) = B_0 K_0(r/\lambda),
\]
where $K_0$ is the modified Bessel function of the second kind.
For small $r$ (i.e. $r << \lambda$), we can approximate this and find that
\[
	K_0 \propto - \ln{(r/\lambda)},
\]
and notice a singularity at $r = 0$.
For the gradient we thus find
\[
	\nabla B \propto \pfrac{K_0}{r}(r/\lambda) \propto \pfrac{-\ln{(r/\lambda)}}{r} = -\lambda/r.
\]
The size of the supercurrent density has the same relation, $J_S \propto 1/r$.

\section{Superconducting wire}
The voltage $U = \SI{1.5e-5}{\volt}$ across the wire of length $\ell = \SI{.08}{\meter}$ induces a current $J_t$. % through the resistive wire with unknown resistivity $\rho$ according to Ohm's law.
Due to the presence of the magnetic field $B = \SI{5}{\tesla}$, if the vortices move with velocity $v_L$, a Lorentz force $f_L$ per vortex acts on the vortices.
This results in a power input $P_L = f_Lv_L = J_tBv_L$ per vortex.
%$\epsilon = Bv_L$
This power should come from the current induced by the voltage, thus $P_L = \epsilon J_t = \frac{U}{\ell}J_t$.
Equating these expressions and rewriting yields
\[
	v_L = \frac{U}{B\ell} = \SI{3.75e5}{\meter\per\second}.
	% https://www.wolframalpha.com/input?i=1.5*10%5E-5+%2F+%285*+.08%29
\]

\section{Fine type-II superconducting wire}

\section{Critical currents}

\section{A weak junction}


\bibliographystyle{vancouver}
\bibliography{references.bib}

%\appendix

\end{document}
