ass2: Finalize assignment 2, hand it in

.gitignore

@@ -44,3 +44,6 @@
 # Assignment 1
 sc_elements.pdf
 superconductivity_assignment1_kvkempen.pdf
+
+# Assignment 2
+superconductivity_assignment2_kvkempen.pdf

makefile

@@ -3,3 +3,6 @@
 ass1:
 	./sc_elements.py
 	latexmk -xelatex superconductivity_assignment1_kvkempen
+
+ass2:
+	latexmk -xelatex superconductivity_assignment2_kvkempen

references.bib

@@ -59,3 +59,52 @@ Electrical conductivity or specific conductance is the reciprocal of electrical
 	editor = {Lide, David R.},
 	year = {2003},
 }
+
+@report{abrikosov,
+	title = {Type {II} superconductors and the vortex lattice},
+	url = {https://www.nobelprize.org/prizes/physics/2003/abrikosov/lecture/},
+	abstract = {The Nobel Prize in Physics 2003 was awarded jointly to Alexei A. Abrikosov, Vitaly L. Ginzburg and Anthony J. Leggett "for pioneering contributions to the theory of superconductors and superfluids".},
+	language = {en-US},
+	urldate = {2022-02-23},
+	institution = {The Nobel Foundation},
+	author = {Abrikosov, Alexei A.},
+	year = {2003},
+	pages = {29--67},
+}
+
+@article{dai-synthesis-1995,
+	title = {Synthesis and neutron powder diffraction study of the superconductor {HgBa2Ca2Cu3O8} + δ by {Tl} substitution},
+	volume = {243},
+	issn = {0921-4534},
+	url = {https://www.sciencedirect.com/science/article/pii/0921453494024618},
+	doi = {10.1016/0921-4534(94)02461-8},
+	abstract = {Substitution of Tl for Hg was performed in the Hg based 1223 phase HgBa2Ca2Cu3O8 + δ (Tc = 135 K), resulting in an increase of the superconducting transition temperature to 138 K for samples with a nominal composition of Hg0.8Tl0.2Ba2Ca2Cu3O8 + δ. The crystal structure of this solid solution has been investigated by neutron powder diffraction techniques at room temperature and at 10 K. The compound has the same crystal as Hg-1223 with the space group symmetry P4/mmm and lattice parameters a = 3.8489(1), c = 15.816(1) Å. Rietveld analysis results indicate that Hg is partially replaced by Tl, and the oxygen content, δ, is 0.33. The lattice-parameter changes resulting from the Tl substitution are too small to account for the Tc change by mimicking the effect of pressure. No phase transition occurs down to 10 K.},
+	language = {en},
+	number = {3},
+	urldate = {2022-03-07},
+	journal = {Physica C: Superconductivity},
+	author = {Dai, P. and Chakoumakos, B. C. and Sun, G. F. and Wong, K. W. and Xin, Y. and Lu, D. F.},
+	month = mar,
+	year = {1995},
+	pages = {201--206},
+}
+
+@article{chubukov,
+author = {Chubukov, Andrey and Pines, David and Schmalian, Jörg},
+year = {2002},
+month = {02},
+pages = {51},
+title = {A Spin Fluctuation Model for D-wave Superconductivity}
+}
+
+@misc{ray-2016, title={Master's thesis: Structural investigation of La(2-x)Sr(x)CuO(4+y) - Following staging as a function of temperature}, url={https://figshare.com/articles/thesis/Structural_investigation_of_La_2_x_Sr_x_CuO_4_y_Following_staging_as_a_function_of_temperature/2075680/2}, DOI={10.6084/m9.figshare.2075680.v2}, abstractNote={A thesis submitted to the Niels Bohr Institute at the Faculty of Science at the University of Copenhagen, Denmark, in partial fulfilment of the requirements for the degree of Master of Science in physics. Submission date was November 19, 2015, and the defence was held on November 30, 2015, where the degree was also awarded.
+
+The cuprate La2-xSrxCuO4+y – a high-temperature superconductor – was discovered almost three decades ago. However the mechanisms behind the superconductivity in the material for different doping values x and y are still not fully understood. A small part of this large puzzle is added to the pile with this thesis, where results on the structure for several different samples are presented.
+
+The emphasis in this thesis is on a specific superstructure thought to be connected to the ordering of interstitial oxygen, known from the isostructural compound La2NiO4+y as staging. Four single crystal samples with different co-doping values are investigated by the use of both X-rays and neutrons.
+
+Staging is observed for all four samples at low temperatures with X-ray measurements. The sample with strontium doping x = 0.00 shows several coexisting staging levels with staging numbers between 2 and 8, with the highest contribution from a staging level between 4 and 5. The co-doped samples show increasing staging number with increasing x. It is found that the staging belongs to a structural phase assumed in space group Fmmm, while the unstaged fraction of the samples are in the Bmab space group. These two structural phases are found to have significantly different lengths of the long crystal axis for the two low x samples, in the order of a fraction of a percent, while the two higher x samples had a difference of only a small fraction of a permille.
+
+The temperature dependent phase transitions for both the Bmab structure and the staging reflections are investigated between 5 and 300 K. The critical exponents for the Bmab reflections are found to be significantly lower than results from similar materials in literature, although with transition temperatures consistent with literature for comparable sample compositions. It is found that the critical exponents for the staging reflections increase for increasing doping while the transition temperatures decrease, both consistent with results on the isostructural La2NiO4+y.
+
+Results from previous neutron measurements are found to be consistent with the X-ray measurements in this work, and measured reciprocal space maps from this work show a large variety of other superstructure reflections which will be interesting to investigate in the future.}, publisher={figshare}, author={Ray, Pia Jensen}, year={2016}, month={Feb} } 

superconductivity_assignment2_kvkempen.tex

@@ -0,0 +1,319 @@
+\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 2}
+\author{
+	Kees van Kempen (s4853628)\\
+	\texttt{k.vankempen@student.science.ru.nl}
+}
+\AtBeginDocument{\maketitle}
+
+% Start from 4
+\setcounter{section}{3}
+
+\begin{document}
+
+\section{Temperature dependence in Landau model}
+In the Landau model, free energy is given as function of order parameter $\psi$ and temperature $T$ as
+\[
+	\mathcal{F} = a(T - T_c) \psi^2 + \frac{\beta}{2}\psi^4.
+\]
+
+The equilibrium state as function of temperature $T$ is the state of minimal free energy with respect to the order parameter $\psi(T)$.
+This point we call $F_0(T)$ with order parameter $\psi_0(T)$.
+For this, we will take the derivative of $F$ with respect to $\psi$ and equate it to zero.
+
+\[
+	0 = \pfrac{\mathcal{F}}{\psi} = \pfrac{}{\psi} \left[ a(T-T_c)\psi^2 + \frac{\beta}{2}\psi^4 \right] = 2a(T-T_c)\psi + 2\beta\psi^3
+\]
+
+Extreme points are found at $\psi = 0$ and $\psi = \pm\sqrt{\frac{-a}{\beta}(T-T_c)}$.
+
+For $T \geq T_c$, $\psi_0(T \geq T_c) = 0$ gives the minimum, i.e. $\mathcal{F}_0(T \geq T_c) = 0$.
+
+For $T \leq T_c$, $\psi_0(T \leq T_c) = \pm\sqrt{\frac{-a}{\beta}(T-T_c)}$ is the minimum,
+giving free energy
+\[
+	\mathcal{F}_0(T \leq T_c) = \frac{-a^2}{\beta}(T-T_c)^2 + \frac{a^2}{2\beta}(T-T_c)^2 = \frac{-a^2}{2\beta}(T-T_c)^2 \leq \mathcal{F}_0(T \geq T_c).
+\]
+
+For the specific heat, we find
+\[
+	C(T) = -T\pfrac{^2\mathcal{F}}{T^2} =
+	\begin{cases}
+		0 & T > T_c \\
+		\frac{a^2}{\beta}T & T < T_c
+	\end{cases}.
+\]
+There is thus a discontinuity in $C(T)$ at $T = T_c$ with size $\Delta C(T) = \frac{a^2}{\beta}T_c$.
+See the following sketches of the $T$-dependence of the derived quantities.
+
+\begin{figure}[H]
+	\centering
+	\includegraphics[width=.9\textwidth]{ass2-landau-theory-T.pdf}
+	\caption{For the Landau theory, we find the drawn temperature dependences for equilibrium values of $\mathcal{F}_0$, $\psi_0$ and $C$. Note the minus sign for the free energy $\mathcal{F}_0$. At $T = 0$, there are  $y$ axis intersections for all three quantities, namely a minimum $\mathcal{F}_0(0) = \frac{-a^2}{2\beta}T_c^2$, $\psi_0(0) = \pm\sqrt{\frac{-a}{b}T_c}$, and $C(0) = 0$, which I forgot to indicate in the sketches. Do also note that there thus is an intersection in the $\mathcal{F}_0(T)$ curve at $T = 0$, although the drawing may look asymptotic.}
+\end{figure}
+
+\section{Type-I superconducting foil}
+\begin{enumerate}
+\item
+The screening equation is given as
+\[
+	\nabla^2\vec{B} = \frac{\vec{B}}{\lambda^2}.
+\]
+For easy of calculation, we will use cartesian coordinates,
+and put the external magnetic field $B_E$ along the $x$ axis:
+$\vec{B_E} = B_E \hat{x}$.
+A foil with thickness $a$ we put parallel to the $xy$ plane with the middle of the thickness at $z = 0$ such that the foil fills $-\frac{a}{2} < z < \frac{a}{2}$.
+Due to symmetry in the $xy$ plane of the system, the field inside the foil can only depend on $z$ coordinates.
+So we define the magnitude of the field $|\vec{B}| = B(z)$. 
+
+Using the screening equation, we look for a solution.
+\[
+	\nabla^2\vec{B} = \left[ \pfrac{^2}{x^2} + \pfrac{^2}{y^2} + \pfrac{^2}{z^2} \right] \vec{B}
+\]
+we realize that we only have $z$ dependence, and $B_y = 0 = B_z$.
+\[
+	\nabla^2\vec{B} = \pfrac{^2B_x}{z^2} \hat{x}
+\]
+Rewriting yields
+\[
+	\vec{B} = \lambda^2 \pfrac{^2B_x}{z^2}\hat{x}.
+\]
+For this we know the general solution:
+\[
+	\vec{B} = B_0 \left[ C \cdot e^{z/\lambda} + D \cdot e^{-z/\lambda} \right],
+\]
+with constants $B_0$, $C$, and $D$.
+
+Now we can apply two boundary conditions to find the solution inside the material.
+
+First, due to mirror symmetry in $z$, we require $B(z) = B(-z)$, giving that $C = D$,
+thus we contract the constants as $B'_0 = CB_0 = DB_0$.
+This allows us to write the exponents into $cosh$ form.
+\[
+	B(z) = B'_0 \left[ e^{z/\lambda} + e^{-z/\lambda} \right] = B'_0 \cosh{\frac{z}{\lambda}}
+\]
+
+Second, just outside the foil, at $z = \pm \frac{a}{2}$, the field must be $B_E$, and the field should be continuous across the boundary:
+\[
+	B_E = B(\frac{a}{2}) = B'_0 \cosh{\frac{a}{2\lambda}} \iff B'_0 = \frac{B_E}{\cosh{\frac{a}{2\lambda}}}
+\]
+
+This gives us our final expression for $B(z)$:
+\[
+	B(z) = 
+	\begin{cases}
+		B_E\frac{1}{\cosh{\frac{a}{2\lambda}}}\cosh{\frac{z}{\lambda}} & |z| \leq \frac{a}{2} \\
+		B_E & |z| \geq \frac{a}{2}
+	\end{cases}.
+\]
+
+The supercurrent follows from the Maxwell-Amp\`ere law, considering that there are no other currents, and we look at a current steady over time ($\pfrac{\vec{E}}{t} = 0$):
+\[
+	\nabla\times\vec{B}(z) = \mu_0\vec{J_s}
+\]
+Reordering and calculating the curl gives:
+\[
+	\vec{J_s} = \frac{1}{\mu_0} \nabla \times (\pfrac{B(z)}{z}\hat{x}) = \frac{B_E}{\mu_0 \lambda \cosh{\frac{a}{2\lambda}}} \sinh{\frac{z}{\lambda}} \hat{y}
+\]
+
+\item
+From the derivation of the Ginzburg-Landau theory, we get the following expression for the supercurrent $\vec{J_s}$:
+\[
+	\vec{J_s} = -\frac{2e\hbar n_s}{m}(\nabla\theta + \frac{2e\vec{A}}{\hbar})
+\]
+Using the rigid gauge, we set $\theta = 0$.
+Next, we can equate the previously found supercurrent for our foil to the Ginzburg-Landau found one and reorder to find $\vec{A}$:
+\[
+	\vec{A} = \frac{-B_E m \sinh{\frac{z}{\lambda}}}{4\lambda\mu_0 e^2 n_s \cosh{\frac{a}{2\lambda}}} \hat{y}
+\]
+This we can rewrite using $\lambda = \sqrt{\frac{m}{\mu_0 n_s e^2}}$ for the London penetration depth as
+\[
+	\vec{A} = \frac{-B_E \lambda \sinh{\frac{z}{\lambda}}}{4 \cosh{\frac{a}{2\lambda}}} \hat{y}.
+\]
+
+\item
+\[
+	\nabla \cdot \vec{A} = \pfrac{A_x}{x} + \pfrac{A_y}{y} + \pfrac{A_z}{z} = \pfrac{0}{x} + \pfrac{A_y}{y} + \pfrac{0}{z} = 0,
+\]
+as $A_y$ is independent of $y$.
+
+In our case, indeed the rigid gauge choice gives the criterium for the London gauge ($\nabla \cdot \vec{A} = 0$).
+
+In the rigid gauge, the order parameter $\psi$ is constant in space and time.
+To then also have that $\nabla \cdot \vec{A} = 0$, follows from the expression for the supercurrent as we saw earlier.
+
+Reversely, assume that $\nabla \cdot \vec{A} = 0$, and look at what conditions need to be met in order to imply rigid gauge.
+Again, we look at the expression for the supercurrent as function of $\theta$ and $\vec{A}$,
+\begin{align*}
+	\vec{J_s} &= -\frac{2e\hbar n_s}{m}(\nabla\theta + \frac{2e\vec{A}}{\hbar}) \\
+	\iff \frac{2e}{\hbar}\vec{A} &= -\frac{m}{2e\hbar n_s} \vec{J_s} - \nabla\theta,
+\end{align*}
+and take the divergence,
+\[
+	\frac{2e}{\hbar}\nabla \cdot \vec{A} = -\frac{m}{2e\hbar n_s} \nabla \cdot \vec{J_s} - \Delta\theta = 0
+	\implies \Delta\theta = -\frac{m}{2e\hbar n_s} \nabla \cdot \vec{J_s}.
+\]
+This is what the London gauge implies.
+Now the question is under what circumstances the rigid gauge follows from the London gauge.
+This is the case for $\nabla \cdot \vec{J_s} = 0$, or in words, when there is conservation of superelectrons.
+If this is not the case (if the divergence is non-zero), there is conversion between normal electrons and superelectrons.
+This would take place if the temperature is lowered, as more superelectrons allow for a larger supercurrent, thus a larger critical magnetic field.
+This result seems to agree with Waldram's conclusion in \cite[p. 24--26]{waldram}.
+
+\item
+We apply a gauge transformation as follows.
+\begin{align}
+	\chi(\vec{r}, t) &= \frac{-\hbar}{2e}(\omega t - \vec{k} \cdot \vec{r}) \\
+	\vec{A} &\to \vec{A} + \nabla\chi = \vec{A} + \frac{\hbar}{2e} \vec{k} \\
+	\phi &\to \phi - \pfrac{\chi}{t} = \phi + \frac{\hbar}{2e} \omega
+\end{align}
+
+\end{enumerate}
+
+\section{Type II superconductors and the vortex lattice}
+In 2003, Alexei Abrikosov was one of the winners of the Nobel Prize in Physics ``for pioneering contributions to the theory of superconductors and superfluids''.
+For this occasion, he gave a lecture called ``Type II superconductors and the vortex lattice''\cite{abrikosov}
+explaining the discoveries that led to the understanding of conventional superconductors.
+To get started, let me first explain what superconductors are.
+
+% Begin copy from philosophy
+Superconductors are characterized by perfect diamagnetism and zero resistance.
+Perfect diamagnetism is the ability by superconductors to have a net zero magnetic field inside.
+If you apply an external magnetic field, this thus means that a superconductor will let a current flow on its inside to generate a field to counteract this external field $\vec{H}$.
+This generated current is called a supercurrent.
+Superconductivity is, however, a phase of the material.
+Superconductors only have these properties below a certain temperature, its critical temperature $T_c$, and can only expel a maximum external magnetic field, its critical magnetic field $B_c(T)$, which is a function of the temperature.
+The zero resistance property follows from the perfect diagmagnetism.
+It is impossible for the material to let these supercurrents flow indefinitely with resistance, as heat would be generated.
+
+The class of superconductors we have a model for, is the class of conventional superconductors.
+In this class, there are two types, called type-I and type-II superconductors.
+% End copy from philosophy
+
+In type-I superconductors, there is only one phase in which the superconductor material exhibits perfect diamagnetism:
+when the externally applied magnetic field $H < B_c(T)$ and $T < T_c$.
+
+In type-II superconductors, there are two phases distinct from the normal conducting state.
+One is the superconducting state which behaves as in type-I superconductors, with critical field $B_{c1}(T)$.
+This state is reached for $T < T_c$ and $B_E < B_{c1}(T)$.
+The other state is a mixed state that allows some flux to pass through the material.
+This passing through is done by creating normally conducting channels throughout the material where a fixed amount of flux can pass through.
+This fixed amount is a multiple of the flux quantum $\Phi_0$.
+The material generates current around these channels in accordance to the Maxwell-Amp\`ere law, conforming to the let through magnetic field inside the vortex and cancelling the field on the outside the vortex.
+
+There are lots of applications for both the perfect diagmagnetism and the zero resistivity.
+There is even a Wiki about them: \url{https://en.wikipedia.org/wiki/Technological\_applications\_of\_superconductivity}.
+What is most notable about these applications, is that maintaining a temperature below the critical temperature is the biggest challenge.
+A real breakthrough for superconductivity would be the discovery of room-temperature superconductors at atmospheric pressure, or materials close to that.
+Currently, the highest $T_c$ material we know is carbonaceous sulfur hydride (\ce{CH8S}) with $T_c = \SI{15}{\degreeCelsius}$ but at a pressure of a whopping $\SI{267}{\giga\pascal}$.
+At atmospheric pressure, the highest $T_c$ material known is a cuprate \cite{dai-synthesis-1995} \ce{HgBa2Ca2Cu3O_{8+\delta}} at $T_c = \SI{135}{\kelvin}$.
+The quest for this breakthrough is intensely researched, although most is experimental.
+The clue is that most of the high $T_c$ materials that are being discovered, are unconventional superconductors.
+As there is no theory for them (yet), the search is mostly educational guessing.
+By trying to find patterns in the previously high $T_c$ materials, similar materials are studied to see if they also exhibit superconductivity.
+One of the patterns is that superconductivity in cuprates is high $T_c$.
+We'll focus on these materials in the following.
+
+Currently, most hopeful candidates are cuprates.
+These materials are made of layers of copper oxides (\ce{CuO2}) alternated with layers of other metal oxides.
+The copper oxide layers are the superconductive layers, and the other metal oxides are used as charge reservoirs, doping electrons (or holes) into the copper oxide layers.
+Due to the geometry of these materials, there is anisotropy in the resistivity of the material.
+Parallel to the layers, superconductivity takes place in the copper oxide layers.
+Perpendicular to the layers, this is not the case.
+
+The behaviour of the material can be tuned by tuning the doping, thus the other metal oxides as mentioned before.
+A typical phase diagram as function of the doping can be seen in figure \ref{fig:cuprate-phase}.
+The material can be steered from being antiferrimagnetic to superconductive by increasing doping.
+
+\begin{figure}
+	\centering
+	\includegraphics{cuprate-phase.pdf}
+	\caption{For high $T_c$ superconducting cuprates, a typical phase diagram as function of doping looks like this.\cite{chubukov} }
+	\label{fig:cuprate-phase}
+\end{figure}
+
+As can be seen, there is an optimal doping fraction for achieving highest $T_c$.
+Aiming for this doping yields the desired material.
+
+Now the question is what direction to search for.
+The timeline in figure \ref{fig:timeline} might give a direction for the most promising types of cuprates to look into.
+It could be, however, that other types have higher $T_c$.
+A lot of creativity is therefore needed to find them.
+
+\begin{figure}
+	\centering
+	\includegraphics[width=\textwidth]{Timeline-of-Superconductivity-from-1900-to-2015.pdf}
+	\caption{The last century, a lot of research has been done in the direction of cuprate superconductivity. Pia Jensen Ray made this overview for his master thesis.\cite{ray-2016} The different paths are different types of cuprates. Please see his thesis for the meaning of the labels. On the right side, an idea of the temperature is givin by comparing it to common cooling agents.}
+	\label{fig:timeline}
+\end{figure}
+
+---
+
+The nobel prize lecture by Abriskosov \cite{abrikosov} was really interesting.
+The start was a good recap of the breakthroughs relevant to conventional superconductivity,\footnote{Why is it that every story on superconductivity includes KGB captivity?}
+but in pages 61--63, the theory is worked through a little quickly.
+I might reread it some times.
+
+\section{Currents inside type-II superconducting cylinder}
+For $B_{c1} < B_E < B_{c2}$, the cylinder of type-II superconductor material is in the mixed state.
+In the mixed or vortex state, superconductors let through a number of finite flux quanta $\Phi_0$.
+Some small regions of the material are not superconducting, but in the normal state.
+Flux passes through these regions in multiples of $\Phi_0$, but usually just one $\Phi_0$ per region,
+and a supercurrent is generated to expel the field from the rest of the material.
+These flux allowing regions are called vortices, due to their shape and direction of current flow.
+Vortices look like channels (or tubes), and supercurrents move around these channels in a spiraling fashion.
+One can visualize this as current through a coil such that on the inside of the coil, the field is in one direction perpendicular to it, and on the outside it is the opposite direction.
+The current direction is governed by the Maxwell-Amp\`ere equation.
+In this case, the current is such that the field inside the cylinder but outside these channels is counteracted.
+
+Please see the figure below for a beautiful drawing.
+It was not specified what the direction of $\vec{B_E}$ was with respect to the cylinder orientation, so I chose what I thought was most reasonable as an example.
+
+\begin{figure}
+	\centering
+	\includegraphics[width=.6\textwidth]{SchermafbeeldingKees-vortex-by-fleur-ahlers.png}
+	\caption{The direction of $\vec{J_s}$ is such that a magnetic field is generated to counteract and even expel the external field outside the vortices inside the material. Around the vortices, that means that the supercurrents run anti-clockwise. The field is then along $\vec{B_E}$ inside the vortices, but along $-\vec{B_E}$ outside the vortices but inside the material. Around the outside border of the cylinder, however, $\vec{J_s}$ runs clockwise and again cancels $\vec{B_E}$ on the inside of the material.}
+\end{figure}
+
+\bibliographystyle{vancouver}
+\bibliography{references.bib}
+
+%\appendix
+
+\end{document}