# Let's speak about Quantum Chemistry! Part 1: Fundamental Classical Mechanics Ideas

in #science4 months ago (edited)

Disclaimer: The reason I wrote this post is that I enjoyed doing so, the challenge to create a good post about one subject I like is enough fuel for me to try... I'm far from being an expert in this field but I enjoyed the research I needed to conduct in order to grasp the basics of the subject. The key concepts and the workflow of the following ideas are based on the first chapter of the book: "Fundamentals of Quantum Chemistry", Molecular Spectroscopy and Modern Electronic Structure Computations, M. Mueller, Kluwer Academic Publishers, 2001.

The foundations of Classical Mechanics are what we could label as everyday experiences, through careful observations scientists made a theoretical model capable of explaining macroscopic phenomena and quantitatively describing it's properties. One key concept needed to be understood is that of a "particle", a particle is a small object which position and speed can be determined with the only limitations of those arising from the uncertainty of the instruments used to measure the particle's properties. The key idea here is that in Classical Mechanics if we can know all the forces acting on the particle, then we can determine precisely where a particle is as well as determine it's motion trajectory.

Now we know this is not true on a microscopic scale, but the results yielded by Classical Mechanics apply so successfully to the bulk scale that the quantum theory has to converge with Classical Mechanics results when we treat macroscopic phenomena. From the Newtonian law of motion, we have that for any inertial reference frame the second derivative of the position times the mass of the particle equals the resultant force acting on the particle:

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\huge&space;\vec{F}=m\cdot&space;\ddot{q}$

Here, F = resultant force, m = mass of the particle, the q with two dots stands for the second derivative of the position in a cartesian, polar or spherical inertial coordinate frame. Most of us are very familiar with the standard approach to solving Classical Mechanics problems. We break the Resultant Force vector into its components in our reference frame, then we solve each component individually and finally add them up using the additive properties of vectors. These types of problems often yield differential equations system which in some cases can be readily integrated.

### Hamiltonian Mechanics

In 1834, the Scottish mathematician William R. Hamilton developed a method under which some complex problems could be solved. In order to grasp the Hamiltonian approach, we have to recall that a conservative system is one in which the forces acting over an object don't vary as time passes, i.e. those forces are just functions of the position of the aftermentioned object. Examples of such conservative forces are the gravity force and the electric force. Since they are functions of a characteristic magnitude of the interacting objects (i.e. mass and electric charge) and the distance between them. A non-conservative force usually is a dissipative one (i.e. they dissipate energy) such as frictional forces.

Another important concept is that of generalized coordinates, a regular coordinate system describes the position of a point in the space by stating the distance to this object from an arbitrary point called the origin. Thus, if we have a particle in a plane, we can describe its position as "x" horizontal units away from the origin and "y" vertical units away from the said origin. The utility of generalized systems is that there are times when the position of an object requires fewer points from what this regular coordinates offer, thus we can describe that object's position in a more simpler and compact set of coordinates called "generalized coordinates".

For instance, let's think about a particle moving on a circle in a 2D plane, we could use x and y coordinates to describe its motion, but we don't need to. The particle is constrained to move on the circular trajectory, hence we can just state the radius from the center of the circle and formulate the angle it covers as a function of time. Thereafter, the particle's movement will be perfectly described using a constant (the radius) and a single variable (the angle) instead of two cartesian variables.

In general, an ordinary n-dimensional coordinate system would require "n" individual variables to describe the position of each particle within it, while the generalized coordinates system would require that for "n" dimensions, "m" particles moving under "k" constraints a total of mb-k generalized coordinates.

With these concepts clearly understood, we can inspect next the definition of the Hamiltonian for a particle (or a lot of particles) under the effects of conservative forces:

Here H is the Hamiltonian, T is obtained from the kinetic energy of the system and V from the potential energy.

The Kinetic Energy of a particle is usually defined as the product of the mass times the half of the square of the velocity, but we shall define it as the dot product of the momentum vector divided by two times the mass of each particle. The potential energy only depends on the positions of each particle and it's usually the trickiest part of writing the Hamiltonian.
$\dpi{100}&space;\bg_white&space;\fn_phv&space;\huge&space;T=\sum&space;\frac{\vec{p}\cdot&space;\vec{p}}{2m}=\sum&space;\frac{p^{2}}{2m}$
Equation 2

For these types of systems, Hamilton proposed that for a generalized coordinate frame, the motion equations can be obtained from the Hamiltonian and the following identities:

### The Hamiltonian in action

Let's consider a vertical thrown particle of mass m:

Gif made by the author using picasion.com and Microsoft PowerPoint

Also, let's first solve the problem using the Newtonian mechanics:

Since the motion takes place alongside a line, we can consider it unidimensional, thus we can take vectorial magnitudes as scalar ones.

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;F_{y}=m\ddot{y}=-mg$

After integrating two times with respect to the variable "t" (the time) we get the familiar result:

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;y(t)=y_{o}+(V(t))t-\frac{1}{2}gt^{2}$

Let's see how we can handle this using the Hamiltonian, we have to make an expression for the kinetic energy of the particle in terms of its momentum plus its potential energy.

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;H_{(p,y)}=T+V=(\frac{p^{2}}{2m})+mgy$

Afterwards, we take the partial derivative of the Hamiltonian with respect to the particle's coordinates and consider it's momentum constant, we get:

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;(\frac{\partial&space;H}{\partial&space;q})_{p}=mg=-\dot{p}_{y}$

Then we pay attention to the partial derivate of the Hamiltonian with respect to the particle's momentum and considered its coordinate variables constant.

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;(\frac{\partial&space;H}{\partial&space;p})_{q}=\frac{p}{m}=\dot{y}$

Since this example is extremely simple, here we get a trivial relationship because the momentum is just the product of the velocity times the mass of the particle, but it is this capability to consider the momentum of the particle independent of its coordinates the feature which makes the Hamiltonian so useful in quantum mechanics (later we'll handle fairly complex examples where the usefulness of the Hamiltonian will be clearly appreciated).

### Ok, but this seems like a Physics class, what's the relationship with Chemistry?

Well, it turns out that variables governing the Harmonic oscillator can be readily described using the Hamiltonian, this Harmonic Oscillator model is a fundamental tool to describe the vibrational energy alongside the chemical bond between the atoms in a molecule.

Let's consider first the simple model where one of the ends of the spring is attached to a fixed location and the other end is attached to a particle of mass m (This model yields good results when the mass of one atom in a diatomic molecule which greatly exceeds the mass of the other, such as Hydrogen Bromide). The mass is bound to move only on the x axis and it will have potential energy according to Hook's law.

Gif made by the author using gifmaker.me and Microsoft PowerPoint

The Hamiltonian for such system is of the form:

$\inline&space;\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;H=\frac{p_{x}^{2}}{2m}+\frac{1}{2}k(x-x_{o})^{2}$

When we take the partial derivative of the Hamiltonian on equation 4 we get:

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;\left&space;(&space;\frac{\partial&space;H}{\partial&space;x}&space;\right&space;)_{p}=k(x-x_{o})=-\dot{p}$

Then we can notice that the derivative on equation 5 produces a trivial result since the momentum is just the velocity times the mass:

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;\left&space;(&space;\frac{\partial&space;H}{\partial&space;p}&space;\right&space;)_{q}=\frac{p}{m}=\dot{x}$

The first derivative allows us to find the trajectory of the particle since the first derivative of the momentum is just the resultant force, or the mass times the acceleration, thus we have:

$\dpi{150}&space;\bg_white&space;\fn_phv&space;\large&space;\frac{dp}{dt}=m\frac{d^{2}x}{dt}=-k(x-x_o)\Rightarrow&space;\ddot{x}=-\left&space;(&space;\frac{k}{m}&space;\right&space;)(x-x_o)$
Equation 5

The solution to this differential equation is well known

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;x_{\left&space;(&space;t&space;\right&space;)}=x_{o}+a\sin&space;\omega&space;t+b\cos&space;\omega&space;t$
Equation 6

It also can be written using imaginary numbers as

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;x_{\left&space;(&space;t&space;\right&space;)}=x_{o}+Ae^{i\omega&space;t}-Be^{-i\omega&space;t}$
Equation 7

Taking a closer look to the above equations, we can know that xo is the equilibrium point of the spring (the initial position), also we notice that omega must have units of the inverse of time since the exponential, sine, and cosine functions should have dimensionless arguments. Furthermore, we can state that this omega is related to the oscillation frequency of the spring. To have an explicit expression for omega we can take the second derivative of both the equation 6 or equation 7 to get:

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;\ddot{x}=-\omega&space;^{2}(x(t)-x_{o})$
Equation 8

From equation 5 and equation 8 we get an expression for omega (the oscillation frequency) in terms of the variables of our model:

$\dpi{100}&space;\bg_white&space;\fn_phv&space;\LARGE&space;\omega&space;=\sqrt{\frac{k}{m}}$

From the fact that the cosine and sine functions oscillate from -1 to 1, the constants (a, b, A, B) on equation 6 and equation 7 are related to the amplitude and the phase of the oscillation. Since we haven't quantized this model yet, there are no constraints for the values of these constants... but the fact that in molecular bonds these constants have a limited set of allowed values, because of the spacial constraints caused by the other bonds and atoms in the surroundings, and because of the nature of the bond itself (for instance simple bonds and double bonds vibration is not the same alongside two carbon atoms), we can identify the identity of the bonds present in a sample via infrared spectroscopy (On this technique, a device called infrared spectrometer is used to apply infrared radiation -the radiation between 0.7-1000 μm, the resulting spectra allow us to identify which bonds are present in the sample since the wavelength absorption range is a unique characteristic of each bond).

### Let's see how the scope of the Harmonic Oscillator model can be broadened

Let's remove the limitation of one mas being greater than the other to describe better diatomic molecules.

Gif made by the author using gifmaker.me and Microsoft PowerPoint

Now, we will consider two masses m1 and m2 separated by a string with a force constant k and equilibrium length xo. We shall see here how taking generalized coordinates instead of ordinary coordinates allows us to simplify the equations.

Gif made by the author using gifmaker.me and Microsoft PowerPoint

The Hamiltonian for such configuration of masses is as follows:

$\dpi{150}&space;\bg_white&space;\fn_phv&space;\large&space;H\left&space;(&space;x_{1},x_{2},p_1,p_2&space;\right&space;)=\frac{p_1^{2}}{2m_{1}}+\frac{p_2^{2}}{2m_2}+\frac{1}{2}k(x_2-x_0-x_1)^{2}$
Equation 9

This Hamiltonian is not as easy to solve as the previous ones, but if we can introduce some convenient generalized coordinates we can simplify it. Let's define r as the length of the spring outside from its equilibrium length and let s be the position of the mass center of the system.

$\dpi{150}&space;\bg_white&space;\fn_phv&space;\large&space;r\equiv&space;x_{2}-x_{1}-x_o$
$\dpi{150}&space;\bg_white&space;\fn_phv&space;\large&space;s\equiv&space;\frac{(m_1x_1+m_2x_2)}{(m_1+m_2)}$

After introducing these simplifications, the potential energy of the system becomes into the product of half the elastic constant times r square (the new coordinate for the displacement of the string from its equilibrium point). After this, we have to transform the momentum using the new coordinates taking the derivative of r and s with respect to the time:

From equation 10 and equation 11 we can find expressions for the time derivative of X1 and X2 with respect to the time (the velocities).

From the last two equations, we can create expressions for the momentum of the system in terms of the generalized coordinates:

$\dpi{200}&space;\bg_white&space;\fn_phv&space;\large&space;p_1=m_1\dot{x}_1=m_1\dot{s}-(\frac{m_1m_2}{m_1+m_2})\dot{r}$
$\dpi{200}&space;\bg_white&space;\fn_phv&space;\large&space;p_2=m_2\dot{x}_2=m_2\dot{s}-(\frac{m_1m_2}{m_1+m_2})\dot{r}$

Further simplification can be made if we define the reduced mass "μ" of this system as the term in parenthesis in the above couple of equations.

Now we can merge all these expressions into our Hamiltonian to get:

$\dpi{300}&space;\bg_white&space;\fn_phv&space;\tiny&space;H(r,s,p_r,p_s)=\frac{p_1^{2}}{2m_1}+\frac{p_2^{2}}{2m_2}+\frac{1}{2}Kr^{2}=\frac{1}{2}\left&space;[&space;\left&space;(&space;m_1+m_2&space;\right&space;)\dot{s}^{2}+\mu&space;\dot{r}^{2}&space;\right&space;]+\frac{1}{2}kr^{2}=\frac{1}{2}\left&space;[&space;M\dot{s}^{2}+\mu&space;\dot{r}^{2}&space;\right&space;]+\frac{1}{2}kr^{2}$
$\dpi{200}&space;\bg_white&space;\fn_phv&space;\small&space;H(r,s,p_r,p_s)=\frac{p_r^{2}}{2\mu&space;}+\frac{p_s^{2}}{2M}+\frac{1}{2}Kr^{2}$

Since "s" is the generalized coordinate for the center of mass, we can assume that the momentum of this coordinate corresponds to the translational kinetic energy. Because we are only interested in the vibrational energy, our final Hamiltonian has the form:

$\dpi{200}&space;\bg_white&space;\fn_phv&space;\small&space;H(r,p_r)=\frac{p_r^{2}}{2\mu&space;}+\frac{1}{2}Kr^{2}$

Our generalized coordinates have proven useful, now we can recognize the shape of this last Hamiltonian as the one we solved previously for the Harmonic Oscillator, thus the solution for this one is also of the form:

$\dpi{200}&space;\bg_white&space;\fn_phv&space;\small&space;r(t)=a\sin(&space;\omega&space;t)+b\cos&space;(\omega&space;t)=Ae^{i\omega&space;t}+Be^{-i\omega&space;t}$
$\dpi{200}&space;\bg_white&space;\fn_phv&space;\small&space;\omega&space;=\sqrt{\frac{k}{\mu&space;}}$

This solution grants us a glimpse of how Classical Mechanics concepts lead us to the doorstep of solving atomic problems, the main goal of this article was to emphasize the usefulness of the Hamiltonian and the generalized coordinates to solve otherwise complicate equations and to lay the foundations to comprehend quantum phenomena. In the next part, we shall discuss a little bit of the fundamental ideas of Quantum Mechanics and how these concepts are applied to describe better properties of atoms and molecules from a chemist's perspective.

All equations were made on codecogs.com latex online editor.
Sort:

I guess it will take me few hours to grind through your post. For now, I want this curation! Now! @tipu curate

Posted using Partiko Android

Don't worry, take your time. I made this post in order to try to explain this subject to people who at least has some basic understanding of integral calculus and some basics on elemental physics. I can break it down further if needed.

I usually work around this topics, but since none of my usual contact list seemed interested in this specific content I never posted any of it before. Since you said that you like this kind of stuff I gave it a try😊

Correct tag is #stemgeek, I hope it's not too late to change the tag 🙂

Posted using Partiko Android

#stemgeek and #steemstem tag left out

Posted using Partiko Android

The #steemstem community has really tough quality criteria, I don't want to spam their tag since they emphasize that users shall post original content.

The ideas on my post are not of my own (99.99% most of the scientific knowledge is from somebody else), I just digested the best I could a very complicated subjected and tried to come up with some understandable ideas for people who has none specialized knowledge about quantum chemistry but with some academic background.

I dunno if they would approve my content or not, even though I followed the copyright guidelines they have for their posts. Anyhow, this is just the first one, I will ask some of the SteemSTEM guys if they like it, and then I will add their tag in the next issues of the series.

You can use it provides you're using it as an educational purpose. If you have a new discovery of a quantum calculation, please don't bring it to #steem because that doesn't get you noble prize, other than a few extra upvote.

Posted using Partiko Android

Hahaha trust me, if I ever (most unlikely since most of my time is expended trying to get enough resources to sustain myself, not the best scenario to develop a scientist career but we don't get to choose the place where we are born 😃. I don't complain much about it, I just keep moving forward with a smile on my face) make such a discovery I would post it as Satoshi Nakamoto did. I really enjoy having a simple life.

You know the math better than any of us. Digital currency or fiat, they need to have a limit for supply and demand in order to function as a value holder(asset). You can have 5000 steem today, one day when steem price is $1000, one of the whale decided to dump his/her 5000000steem, and your 5000steem will instantly get back to worth only$1000 for the all of them. Can you imagine how many of them really managed to sell their BTC at 20k?

But again, no right or wrong. Different measure for different situation. And congratulations on the curation 🙂

Posted using Partiko Android

4 months ago (edited)

I know, if this structure is meant to be successful then having great influence on it would prove more worthy than the net value of your so holdings in $terms. I didn't post this article because of the curation, of course it is nice to get rather big upvotes... But I posted this stuff cuz you said you liked the subject, and it would be super interesting to have a debate about this and the future topics (at least for me here is where true value lies on). But giving it a second thought, if a nobel price would ever start posting his/her content on Steem and make the coin value rice to 1,000$. My current SP would be as worthy as the novel prize award lol

Thanks for using eSteem!
Your post has been voted as a part of eSteem encouragement program. Keep up the good work! Install Android, iOS Mobile app or Windows, Mac, Linux Surfer app, if you haven't already!

Finish reading. Took me 5 hours just to grind through once, and another 2 hours to reread. Still, like most of the formula I have absolutely no idea, and I did not even try to solve them.

What intriguing to me, is the bounce and the radius. Coincidentally, this is what exactly a warp speed concept about. Create an explosion and bounce us to the "front". In space, there's no front, so we bend the time in front, so we warp over the radius 😂 I don't know how to explain the picture in my mind now after reading your article, but I hope sometime I will be able to write them down with my own word, with my own drawing.

By the way, I failed my maths, and I did not take up my chemistry nor physics. Instead I went for accounting, but that failed too…eventually 😂 and I graduated in computer science, supposed to be a programmer, and I'm a marketing guy in an insurance firm. I am still very much a sci-fi guy, but the theory is driving me nuts.

Posted using Partiko Android

How cool! I would like to learn how to make codes, we could help each other in this issue. You give me a hand to learn how to code and I'll give you a hand to comprehend what is the science part of current sci-fi ideas... but keep in mind that what once has been a dream for mankind then we have made it a routinary activity in our daily lifestyles.

The important issue here is to actually know how to start from the bulk observations and the models that work well at the macroscopic scale to start describing the atomic world. General and Special Relativity is another issue (I will post about that as well since it is a mandatory subject in my Modern Physics course), but have some patience this is a subject worthy of the effort to be understood. I will make a comment breaking down the math of each equation step by step, I didn't write this on the post because it would be too lengthy.

4 months ago (edited)

This post was manually curated by @upmewhale, in combined efforts with the eSteem curation team to bring further support to our valued Steem community! Post curated, courtesy of @davidke20

~ eSteem Curation Team

Use #steemstem tag for science blogs

Hola, aqui en mi tierra para llegar a esto, tenemos que fumarnos minimo un lumpia. De verdad estoy boqui abierto. jejejejejjeje

Lo felicito, a tu universidad también. Como le dije una vez, estas en la Simón Bolívar, son grandes ligas. Y que bueno que ven física moderna y química cuantica. Cuanda yo estudie se veia química descriptiva, cualitativa y cuantitativa, osea la ciencia del siglo antes pasado......

;-)