A new Kilogram? What? How?

A new Kilogram? What? How?

To understand this, let us see how time is defined. We all know the SI unit of time is seconds, and it is defined as the international unit of time, the second, is defined by measuring the electronic transition frequency of caesium atoms. Similarly length (SI unit metre) is defined as the length traveled by light in $\frac{1}{3 \times 10^8}$ seconds or to be more precise $\frac{1}{299792458}$ seconds.
The inspiration to write this blog was derived from Veritasium|How We’re Redefining the kg.
We see that, all these are fundamental quantities are define on a given/fixed standard.
However, the kilogram is set as the weight of a metal cylinder in Paris. That is not that ideal, huh!
So, the NIST is trying to define the unit of mass, i.e., kilogram using some already standardized quantities like the Planck’s constant and the Avogadro’s number. The Planck’s constant takes a dive into advanced applied physics (no, not Quantum physics :p ); therefore, demands my blog’s space and time.
The engineering behind it is pretty simple, to put it in the simplest of the terms;

There is a balancing device, which has a mass unit and a coil
unit. The mass unit is balanced with the magnetic field from
the coil using a motor fixed to it, unit and unless, $F_{magnetic} = F_{gravitational}$. To know more, read Watt Balance or Kibble Balance.

Let’s look into the working now, in this process we will generate some cool equations and in the process will learn some new concepts, as and when required.
First, the watt balance, the principle of operation itself says that,
$$F_{magnetic} = F_{gravitational}$$
Taking into account all the usual scientific nomenclature, we can write,
$$\boxed{m \times g = B \times i \times l} \tag*{(1)}$$
Here, $m$ is the mass, $g$ is the acceleration due to gravity, $B$ is the magnetic flux density (in Tesla($T$)), $i$ is the current in the coil and $l$ length of the conductor in the field.
Equation $(1)$ is for the weighing mode operation of Watt Balance. In which, the weights are matched on both sides.
There is another mode, called the velocity mode in which the mass($m$) is lifted at a height and then the coil is moved back and forth in the magnetic field. This motion induces a voltage, $V$, in the coil. By Faraday’s motional emf expression,
$$\boxed{V = B \times l \times v} \tag*{(2)}$$
Here, $v$ is the velocity of the conductor in the magnetic field.
Now, $(1)$ can be written as,
$$B \times l = \frac{m \times g}{i} \tag*{(3)}$$
and $(2)$ can be written as,
$$B \times l = \frac{V}{v} \tag*{(4)}$$
Equating $(3)$ and $(4)$, we have,
$$\frac{m \times g}{i} = \frac{V}{v}$$
Which can be written as,
$$\boxed{m \times g \times v = V \times i} \tag*{(5)}$$
Interestingly, $m \times g \times v$ is the mechanical power (refer Encyclopedia of Electrochemical Power Sources for more info) and $V \times i$ is the electrical power.
You must all be thinking, how does Planck’s constant come to play. Well, please hold your horses. We are almost there.
To measure $V$ in $(5)$, we go into a concept of superconductivity. Called the Josephson Phenomenon. Please read up on its working here. It is also significant because of being the standard of Voltage.
When DC voltage is applied to a Josephson Junction, the junction experiences an oscillation of frequency(read more),
$$f = \frac{2eV}{h}$$
Here, $f$ is the frequency, $e$ is the elementary charge, and $h$ is the Planck’s constant($6.62607004 \times 10^{-34} m^2 kg \cdot s^{-1}$).
The above equation can be written as,
$$\boxed{V = \frac{h\times f}{2e}} \tag*{(6)}$$
For many Junctions ( say $n$) it is,
$$V = n\frac{h\times f}{2e}$$
The Voltage measure here is accurate to $10^{10}$ parts, refer this.
Now, if we write, $i$ as $\frac{V}{R}$ (where $R$ is the resistance offerered by the junction), then, $(5)$ changes to,
$$\boxed{\frac{h^2\times f^2}{4e^2R} = m \times g \times v} \tag*{(7)}$$
Another question now is how do we measure $R$ ?, for that we will use the idea of Quantum Hall Effect. Quantum Hall effect is the standard for resistance, please refer the paper for more information. But suffice to say, the resistance is defined as,
$$\boxed{R = \frac{1}{k}\frac{h}{e^2}}$$
Here, $k$ is $1, 2, 3, ...$
Please note that, without the integer fraction i.e., $R_k = \frac{h}{e^2}$ is called the von-Klitzing constant, this guy got a Nobel for this. Please read more here.

Using the above equation in $(7)$ we get,
$$\frac{h^2\times f^2}{4e^2 \frac{1}{k}\frac{h}{e^2}} = m \times g \times v$$
Which comes to,
$$h^2\times f^2 = m \times g \times v \times 4e^2 \frac{1}{k}\frac{h}{e^2}$$
$$\implies h\times f^2 \times k = m \times g \times v \times 4$$
$$\implies h = \frac{m \times g \times v \times 4}{f^2 \times k} \tag*{(8)}$$
This is for a single Josephson Junction, for $n$ junctions, it becomes,
$$\implies h = \frac{m \times g \times v \times 4}{f^2 \times k \times n^2}$$
Which can look more elegant if we write it as,
$$\boxed{h = \frac{4}{k \times n^2}\frac{g \times v}{f}m}$$
We have seen that $V$ and $R$ were measured very accurately. Similarly, there is a need that we measure the factors very accurately as well. The $v$ is measured using a Laser Interferometer. The $g$ was measured using a Gravimeter.
So, the scientists in NIST, are just putting in some mass $m$ and get the $h$, and keep tuning the value of $m$ till we get a very accurate $h$.
Cheers!

Study of Brachistochrone Problem.

This blog will deal with one of the most elegant topics in mathematics. The famous Brachistochrone curve. Vsuace did a great job at explaining it in this video. This looks good to me, but it doesn’t feel good; because of course the video does not contain tons of fancy mathematical vocabulary and millions of lines of $\LaTeX$.

To get started, we will just put forward a small fact that, this topic is from a larger field of mathematics called the Calculus of Variations.

Before we get into the Analysis of Brachistochrone Curve (which happens to be the heart of this post), let us brush-up on some concepts.

---

Partial Differentiation (advanced preliminary):

Let us look at how we define the derivative of a given function, say we have $g(x)$ and we have,
$$\boxed{\frac{\mathrm{d}}{\mathrm{d}x}g(x)=g'(x)=\lim_{c \to 0}\frac{f(x+ c)-f(x)}{c}}$$
Provided this limit exists.

Now, for a function of many variables it is not easy to compute the total derivative (the usual derivative). Therefore, we calculate the Partial Derivative. You might have already guessed it by now, we define partial derivative as, say we have a function $h(x,y,z)$, then,
$$\boxed{\frac {\partial}{\partial x}h(x,y,z) = \lim_{c \to 0} \frac{h(x+c, y,z)-f(x,y,z)}{c}}$$
Let us look into a small example now, say we have,
$$f(x,y,z) = x^2y^3+x+y+z+\sin(x) \cos(y)$$
Now,
$$\frac{\partial}{\partial x}f(x,y,z) = 2xy^3+1+0+0+\cos(x) \cos(y)$$
similarly,
$$\frac{\partial}{\partial y}f(x,y,z) = 3x^2y^2+0+1+0-\sin(x) \sin(y)$$
and,
$$\frac{\partial}{\partial z}f(x,y,z) = 0 +0+0+1+0$$

---

Idea of Speed and Time (basic preliminary):

Let us consider, the following scenario, (Made using Geogebra)

We have, in the figure we have denoted the path ACB as $\pi_1$ and the path ADB as $\pi_2$.
Consider that, The time taken for a body to move from A to B using path $\pi_1$ is $t_{AB_{\pi_1}}$ similarly, the time taken for a body to move from A to B using path $\pi_2$ is $t_{AB_{\pi_2}}$.
Also, consider that the distance for $\pi_1$ be $l_1$ and the distance for $\pi_2$ be $l_2$.

Then,
Speed ($\pi_1$) =$v_1$= $\frac{l_1}{t_{AB_{\pi_1}}}$ and Speed ($\pi_2$) =$v_2$ =$\frac{l_2}{t_{AB_{\pi_2}}}$.

3. Idea of Distance(basic preliminary):

Distance between two points $(a,b)$ and $(c,d)$ is given by,
$$l = \sqrt{(a-c)^2+(b-d)^2}$$
In case, the distance is measured from origin ($O(0,0)$) to some point $A(x,y)$, we get,
$$l = \sqrt{(x-0)^2+(y-0)^2} = \sqrt{x^2 + y^2}$$

---

Now that the notion is clear, we shall proceed into the analysis.
[NOTE : Some topics which are new to me (or you) will be explained as and when it is required].

---

Time needed for a body to travel from $A$ to $B$ is given by,
$t_{AB}=\frac{l}{v}$ for a linear path.

For a curve or a non-linear path, we consider piece-wise linear distance ($\mathrm{d}l$) and speed ($v$). We can define, $t_{AB}$ as,
$$\boxed{t_{AB} = \int_{A}^{B}\frac{\mathrm{d}l}{v}} \tag*{(1)}$$
Now, we must understand a fact that, all distance is composed of $x$ and $y$ coordinates. We consider the point $A$ as the origin ($O(0,0)$) and $B$ as $(x,y)$. Hence,
$$l = \sqrt{x^2 + y^2}\\ \implies \mathrm{d}l = \sqrt{\mathrm{d}x^2 + \mathrm{d}y^2}\\ \implies \mathrm{d}l = \sqrt{(1 + \frac{\mathrm{d}y^2}{\mathrm{d}x^2})\mathrm{d}x^2}$$
$$\boxed{\implies \mathrm{d}l = \sqrt{(1 + \frac{\mathrm{d}y^2}{\mathrm{d}x^2})}\mathrm{d}x} \tag*{(2)}$$

---

Also, since the translation is both in the $x$ and $y$ axes, we can say that the speed is gained by $KE$ and $PE$ equality,
$$KE = PE\\ \implies \frac{1}{2}mv^2 = mgy \\ $$
$$\implies \boxed{v = \sqrt{2gy}} \tag*{(3)}$$

---

Use $(2)$ and $(3)$ in $(1)$, we have,
$$\boxed{t_{AB} = \int_{A}^{B}\frac{\sqrt{(1 + \frac{\mathrm{d}y^2}{\mathrm{d}x^2})}\mathrm{d}x}{\sqrt{2gy}}} \tag*{(4)}$$
We write
$$\frac{\mathrm{d}y}{\mathrm{d}x}=y'$$
Now, $(4)$ becomes,
$$t_{AB} = \int_{A}^{B}\frac{\sqrt{(1 +y'^2 )}}{\sqrt{2gy}}\mathrm{d}x \tag*{(4.1)}$$
Therefore, the function is,
$$\theta = \frac{\sqrt{(1 +y'^2 )}}{\sqrt{2gy}}\\ \boxed{\theta = (1 +y'^2 )^{1/2}(2gy)^{-1/2}} \tag*{(5)}$$

---

Please notice that, the integral $t_{AB}$ can be written as,
$$t_{AB} = \int_{A}^{B} f(y, y', x) \mathrm{d}x \tag*{(6)}$$
In our case, $f(y, y', x)$ is $\theta$.

Equation$(6)$ closely follows the The Euler-Lagrange equation.

It is pretty simple to define,

---

The Euler-Lagrange equation

For any $I$ such that,
$$I = \int_{a}^{b}f(y,y',x)\mathrm{d}x$$
Where, $y' = \frac{\mathrm{d}y}{\mathrm{d}x}$, then $I$ has an stationary point, if the Euler-Lagrange Equation, given by,
$$\frac{\partial f}{\partial y}-\frac{\mathrm{d}}{\mathrm{d}x}\frac{\partial f}{\partial y'}=0$$
is satisfied.

Therefore, for our analysis,
$$\boxed{\frac{\partial \theta}{\partial y} - \frac{\mathrm{d}}{\mathrm{d}x}\frac{\partial \theta}{\partial y'}=0} \tag*{(7)}$$

---

Stationary value

This is a value at the stationary point.
A stationary point $a$ is the point at which the first derivative of a function becomes zero,
$$f'(x)|_{x=a} = 0 ; a = \text{Stationary point}$$
Side note :
Find the stationary points of $f(x) = x^3-6x^2+9x$.

---

Beltrami Identity

Our $\theta$ (from $(5)$) is such a cool expression, because $x$ does not appear in that; therefore, of course $\boxed{\frac{\partial \theta}{\partial x} = 0}$; this again leads us to another beautiful form, Beltrami Identity, which is given by,
$$\boxed{\theta - y' \frac{\partial \theta}{\partial y'}=C ; C \text{ is a constant}}.\tag*{(8)}$$

---

We now find, $\frac{\partial \theta}{\partial y'}$,
$$\frac{\partial \theta}{\partial y'}=\frac{\partial}{\partial y'}[(1 +y'^2 )^{1/2}(2gy)^{-1/2}] \\ \implies \frac{\partial \theta}{\partial y'} =(2gy)^{-1/2} \frac{1}{2}(1 +y'^2 )^{-1/2}\cdot 2y'$$
$$\therefore \boxed{\frac{\partial \theta}{\partial y'} = (2gy)^{-1/2} (1 +y'^2 )^{-1/2}\cdot y'} \tag*{(9)}$$

---

Use $\theta$ and $\frac{\partial \theta}{\partial y'}$ in equation $(8)$, we get,
$$(1 +y'^2 )^{1/2}(2gy)^{-1/2}-y' \cdot (2gy)^{-1/2} (1 +y'^2 )^{-1/2}\cdot y' = C$$
$$\implies \frac{\sqrt{(1 +y'^2 )}}{\sqrt{2gy}}-\frac{y'^2}{\sqrt{2gy}\sqrt{1+y'^2}} = C$$
$$\implies \frac{1+y'^2-y'^2}{\sqrt{2gy}\sqrt{1+y'^2}}=C$$
$$\implies \boxed{\frac{1}{\sqrt{2gy}\sqrt{1+y'^2}}=C} \tag*{(10)}$$
Rearranging a little gives,
$$\sqrt{1+y'^2}=\frac{1}{\sqrt{2gy}C}$$
$$\implies \boxed{(1+\frac{\mathrm{d}y}{\mathrm{d}x})y = \frac{1}{2gC^2} = M^2} \tag*{(11)}$$
This is the equation of the cycloid as per [4], I am yet to figure out how this happened.
The solution of $(11)$ can be found using a parametric equation,
$$\boxed{x(t) = \frac{1}{2}M^2(t-\sin t)\\ y(t)=\frac{1}{2}M^2(1-\cos t)} \tag*{(12)}$$
To derive the above equation, please refer Math-stackexchange|Derive the parametric equations of a cycloid.

I have even plotted the solutions at $(12)$, we can see that it is a cycloid,

If you wish to play with the visualisation please go to Cycloid | Parametric @ Desmos by Pragyaditya

Cheers!

---

References:
1. Stationary Points
2. Euler Lagrange Equation
3. Derivation of Beltrami Identity
4. Introduction to calculus of variations
5. Brachistochrone @ Wolfram