## 1. Introduction

A wave is the propagation of oscillations in space through the coupling between individual oscillators. There is an obvious experiment to hit the right piano key, which sounds the tuning-fork tuned to the same frequency as the key we hit. Another example is breaking a fragile wine glass by the intense sound of the frequency that corresponds to the resonant frequency of the glass.

The crystal of the substance is formed by the regular arrangement of individual atoms, ions, or molecules that are bound together by electrical forces. Since the crystal keeps its form, it means that there are attractive forces between the particles. Each particle has its equilibrium position, which is given by the minimum potential energy in the force field of the surrounding particles. The particles oscillate around the equilibrium position concerning the homogeneity of the substance. Besides, the frequencies of the self-oscillations of all particles are the same. When one particle oscillates, energy transfers to the adjacent particle by resonance, and then the initial excitation travels through space. This excitation propagating is called the * wave*.

First, we show the oscillations of a system consisting of two or three particles. Then, we show the energy transfer between two particles and finally, the wave propagation in an infinite chain of the same oscillators.

### 1.1 Molecular vibration spectra

Radiation passing through the material medium vibrates the particles of the substance that are in resonance with the radiation. The particles then absorb the energy of the radiation, and thus, the radiation is attenuated. This phenomenon is used in spectroscopy to identify the particles with common properties in a specific medium.

#### 1.1.1 The frequency of vibrations of the diatomic molecule

One of the oscillators of a gaseous environment is a molecule. The molecules consist of the individual atoms bound together. Thus, the system of bonded atoms represents an oscillator. The simplest is a diatomic molecule.

Consider a molecule in which two atoms with masses of _{1}, _{2} are bound by the force * F* = −

*(*k

r - r

_{0}), where the

*is the stiffness coefficient of the molecule, and the*k

r

_{0}is the equilibrium distance of the centers of the atoms. If we denote

r

_{1}and

r

_{2}the coordinates of atoms, in the coordinate system of a straight line passing through the centers of atoms, the equations of motion of both atoms are.

Subtracting the equation, we get.

or

The solution to the equation is a harmonic function of * r* =

r

_{m}sin

*.*ω t

These are vibrations of a molecule with an angular frequency.

where

The molecule, thus, represents a simple oscillating system.

If we apply an external periodic force, then we can achieve resonance when a molecule or array of molecules absorbs the energy. This phenomenon finds its place in spectroscopy. If we irradiate the gas consisting of the diatomic molecules with electromagnetic radiation, the resonance wavelengths are attenuated in the transmitted radiation.

Thus, we can determine the stiffness coefficient of the molecule by measuring the resonance frequency.

Example 1 The vibration frequency of some molecules.

Oxygen O_{2}: * m* = 16 × 1.67 × 10

^{−27}kg,

*=*μ

*/2,*m

*= 1133 N/m, hence*k

*= 46.3 THz, wavelength in vacuum*f

*= 6.47 μm.*λ

Nitrogen N_{2}: * m* = 17 × 1.67 × 10

^{−27}kg,

*=*μ

*/2,*m

*= 2287 N/m, hence*k

*= 63.9 THz, wavelength in vacuum*f

*= 4.70 μm.*λ

From this, we can see that the air in the atmosphere absorbs infrared radiation.

Example 2 Stiffness coefficient of the molecule.

The HCl molecule has the resonance frequency of * f* = 86.7 THz (wavelength 5.33 μm). For

m

_{H}= 1.67 × 10

^{−27}kg and

m

_{Cl}= 35 × 1.67 × 10

^{−27}kg,

*= 1.62 × 10*μ

^{−27}kg, we get the

*≈ 481 N/m.*k

The NO molecule has the resonance frequency of * f* = 56.3 THz (wavelength 3.46 μm). For

m

_{N}= 17 × 1.67 × 10

^{−27}kg and

m

_{O}= 16 × 1.67 × 10

^{−27}kg,

*= 13.76 × 10*μ

^{−27}kg, we get the

*≈ 1722 N/m.*k

The CO molecule has the resonance frequency of * f* = 66.3 THz (wavelength 4.52 μm). Pre

m

_{C}= 12 × 1.67 × 10

^{−27}kg and

m

_{O}= 16 × 1.67 × 10

^{−27}kg,

*= 11.45 × 10*μ

^{−27}kg, we get the

*≈ 1987 N/m.*k

Example 3 Coupling stiffness coefficient.

The bonds of pairs of atoms in more complex molecules also show resonance character. The binding force and hence the frequency of the oscillations vary for different binding strengths, e.g., carbon bonds: -C ≡ C- (* λ* = 4.73 μm), =C=C = (

*= 5.90 μm), ≡C-C≡ (*λ

*= 9.48 μm), heterogeneous bonds ≡C-H (*λ

*= 3.00 μm), -O-H (*λ

*= 2.80 μm), -C ≡ N (*λ

*= 4.40 μm), =C=N- (*λ

*= 6.30 μm).*λ

#### 1.1.2 Oscillations of multiatomic molecules

In the case of multiatomic molecules, vibrations of individual pairs of atoms, but also other collective stationary modes will occur. The spectrum of such molecules is more complex, and its theoretical calculation is complicated.

As a simple case, consider a CO_{2} molecule that has a simple linear geometry. The O=C=O atoms are on one line, as opposed, e.g., to the two-dimensional H_{2}O molecule, in which the H-O-H atoms form a “V” structure with an opening angle of 104.5° and an O atom at the midpoint. Another example is the 3-dimensional molecule of NH_{3}. A similar linear structure has an N_{2}O molecule having the arrangement of N ≡ N^{+}-O^{−} or N^{−} = N^{+} = O.

Consider the general case of a triple atom in a linear arrangement, see Figure 1.

The equations of motion are.

The stationary mode represents oscillations of the system in which all particles vibrate at the same frequency and constant amplitudes. If we use phasor symbolism, we get.

After substituting into differential equations, we get a system of algebraic equations.

The condition of nontrivial solution is zero value of the system determinant (characteristic equation)

This is the equation of the third degree for ^{2}.

The first solution is ^{2} = 0 and _{1} = _{2} = _{3}. It represents the translation of the molecule by uniform motion.

The other two solutions are the solution of the quadratic equation for ^{2}

which has two solutions.

Example 4 The CO_{2} molecule oscillations.

As an example, we will find the stationary mode frequencies of the linear molecule O=C=O (carbon dioxide), for which the _{1} = _{2} a _{1} = _{3}. If we substitute it to the corresponding relation, we get.

from where

_{1} = 40.92 THz (_{1} = 7.33 μm), _{1} = 1766 N/m, _{2} = 78.35 THz (_{2} = 3.83 μm).

The first oscillating mode with the angular frequency _{1} represents the simultaneous movement of the outer atoms _{1} = _{3} and the movement of the central atom in the opposite direction. At the same time, the mass center remains at rest. In the second mode _{2}, the middle C atom is at rest, and the other O atoms oscillate against each other (_{1} = −_{3}). In the CO_{2} absorption spectrum, these absorption lines correspond to two wavelengths of _{1} and _{2}.

In the case of two or three-dimensional molecules, there are several stationary modes. The complex structure of the spectrum of resonance frequencies is unique to each molecule and thus allows identifying the substance by measuring the resonance spectrum. This method represents spectral analysis. As can be seen from the above examples, the vibration modes have frequencies that belong to the region of infrared radiation (in the wavelength range from one to hundreds of μm).

Simple molecules of CO_{2}, SO_{2}, N_{x}O_{y}, and others, which enter the atmosphere as a result of human activity, absorb thermal radiation and create a * greenhouse effect*. These gases are so-called

*.*greenhouse gases

The gases such as H_{2}, O_{2,} and N_{2} have non-polar molecules. These molecules have zero dipole moment as opposed to the polar molecules of H_{2}O, CO_{2}, and others. Therefore, they interact very weakly with electromagnetic radiation, and they do not absorb the heat. Thus, the elementary atmospheric gases O_{2}, N_{2} do not cause a greenhouse effect.

### 1.2 Energy transfer in the system of the coupled oscillators

As stated earlier, the nature of the wave is the propagation of oscillations in space. Let us consider a pair of coupled oscillators and then an extension to an infinite line chain of oscillators.

#### 1.2.1 Two coupled oscillation systems

Consider a mechanical model of a pair of particles of an equal mass of the * m*bound to each other and stationary walls by springs with a stiffness of the

*, Figure 2.*k

Suppose that particles can deviate from equilibrium positions in the longitudinal * x*-direction. Then, denote the displacements as

x

_{1}and

x

_{2}. The following equations describe the motion of the particles:

or

where

Consequently, we rewrite the equations to the form.

which has a harmonic solution

where

The * A*and

_{ij}

*are constants, which values are determined by the initial conditions such as displacements of the*α

_{ij}

x

_{10},

x

_{20}and velocities of the

v

_{10},

v

_{20}at the time of

*= 0 and Eq. (8).*t

#### 1.2.1.1 Stationary modes

Stationary oscillations are oscillations with a single frequency and a constant amplitude. There are two stationary modes in the system - one with a frequency of _{1} and the other with a frequency of _{2}.

For the first mode, we have.

and after substituting to the (8), we get.

The first stationary mode is characterized in that both particles vibrate with the same amplitude and the same phase, i.e., together.

The second stationary mode is then.

and after substituting into the (8), we get.

The second stationary mode characterizes two particles oscillating with the same amplitude but the opposite phase, which means opposite direction to each other.

_{2}, which is represented by a significant minimum in the spectral transmission characteristic of the substance.

#### 1.2.1.2 Stationary mode interference

Any movement of the described system is a superposition of stationary modes.

The oscillation character depends on the initial conditions such as the displacement of _{10}, _{20}, _{10} a _{20} and the particle velocities at the time * t* = 0.

As an example, we assume the case, when the first particle is shifted from the equilibrium and the second one is quiescent. The initial conditions are _{10} = _{m}, _{20} = 0, _{10} = _{20} = 0. After substituting the initial conditions, we get.

from where, we obtain amplitudes of _{1} = _{2} = * A* =

x

_{m}/2 and phases

α

_{1}=

α

_{2}= π/2 rad.

Time courses of displacements are in Figure 3.

It is seen from relations (13) and Figure 3 that oscillations are periodically transmitted from one body to another and backward. In the figure, the envelope curves of both oscillators are indicated by a dashed line. As a result of the coupling, there is a periodic exchange of energy between oscillators. At the time of * t* = 0 s, the first particle has the entire energy of the oscillations. On the other hand, in time of.

the whole energy is transferred to the second particle.

Example 2.5 Oscillation energy transfer rate.

Let us assume the aluminum ions from Example 1.5 have the angular frequency of _{0} ≈ 6.1 × 10^{13} rad/s and a mutual distance * a* ≈ 2.8 × 10

^{−10}m. If we suppose the time of transfer

t

_{1}≈ 1.4 × 10

^{−13}s, the transfer rate is

*=*v

*/*a

t

_{1}≈ 2.0 × 10

^{3}m/s, which corresponds to the order of the ultrasound velocity in aluminum (5.1 × 10

^{3}m·s

^{−1}).

As we see from a simple example, the energy of oscillations in a substance is transmitted due to mutual bonds, which are manifested externally as a coherent propagation of mechanical waves, e.g., sound and ultrasound, or incoherent heat propagation.

#### 1.2.2 Propagation of oscillations in a long particle chain

Substances consist of a great number of less or more ordered fundamental particles, which can be molecules, atoms, or ions. Solids and liquids are formed by particles whose attraction to each other is sufficiently large to form a compact body. Each particle of solid or liquid substance has its equilibrium position, which is caused by the force of the surrounding particles. These forces are mainly electrical. The equilibrium position corresponds to the minimum potential energy, and around the minimum of the potential energy, each particle performs an oscillating movement. As a simplified crystal model, we will use an infinite chain of bound particles arranged in a single line.

#### 1.2.2.1 Propagation of longitudinal oscillations along the chain

We see the line chain model in Figure 4 suppose that at the beginning, only the leftmost particle is put into oscillating motion, and the remaining particles of the chain are quiescent. Let us further consider that an external time-dependent force keeps the extreme left particle in a steady oscillating motion with time dependence.

The vibrations of the extreme left particle are gradually transferred to the other particles of the chain via bonds. After some time, all the particles of the system oscillate in harmonic oscillating motion with the same angular frequency of * ω*, but different amplitudes and different phases.

Consider three adjacent particles with indexes * i*- 1,

*,*i

*+ 1, indicated by a dotted fill in the figure.*i

The equation of motion of the * i*-th particle has the form.

which can be expressed as.

where

If the particle chain is very long (theoretically infinite), then, the relationship between neighboring particles does not depend on * i*index.

We can use phasor annotation for the harmonic solution of the equation:

Subsequently, we express the relationship between complex particle displacements by the constant using complex exponential form.

where =

*+ j*β

*.*α

Substituting to (15) we get.

and from there.

Let us consider cosh * γ*=cosh

*cos*β

*+ jsin*α

*sinh*α

*. As right side of the (17) is real, the imaginary part of cosh*β

**is zero. Thus, we get sin**γ

*= 0 or sinh*α

*= 0.*β

In the case that sinh * β* = 0, then

*= 0, cosh*β

*= 1 and cosh*β

**= cos**γ

*:*α

From −1 ≤ cos * α* ≤ 1, we get the corresponding interval of angular frequencies:

where

From (16), we get _{i + 1} = * x*.

_{i}

All particles oscillate with the same amplitude, but the neighboring particles have different phase shift * α*of their oscillations. Thus, there is propagating

*along the chain.*the undamped traveling harmonic wave

The phase shift * α*corresponds to a time shift of Δ

*=*t

*/*α

*. If the distance between the adjacent particles is*ω

*, the phase velocity is*a

For small angular frequencies * ω* <<

ω

_{m}, we can use an approximation for

*<< 1, and thus,*x

The low-frequency waves propagate along the chain having a constant phase velocity of _{0}.

In the case of * ω* >

ω

_{m}, then cosh

**< −1. As comes from cosh**γ

*≥ 1, we can fulfill (17), if sin*x

*= 0 and cos*a

*= −1.*a

The (17) gets then.

The relationship between the displacements of two adjacent particles has the form.

The adjacent particles oscillate with opposite phases and the amplitude of the oscillation exponentially decreases along the chain. Thus, the * exponentially damped standing wave*is arising in the chain.

#### 1.2.2.2 Wave energy transfer in chain of particles

To investigate oscillation propagation along the particle chain, it is important to observe the transfer of energy in this bound system. Let us assume that an external source generates oscillations according to (14). We denote these particle oscillations as _{0}.

For subcritical frequencies * ω* <

ω

_{k}, the power of the source is.

The power has time-dependent parts and constant parts. The constant part represents a steady energy flow.

At angular frequencies of * ω* <

ω

_{m}, the chain can transfer energy without any losses.

In the case of supercritical angular frequencies of * ω* >

ω

_{k}, the power is.

We can see that the power mean value is * ω* >

ω

_{m}, the chain cannot transfer energy, and thus, the power of the source is reactive.

#### 1.2.2.3 Heat transfer in substances

A substance consists of many particles (atoms, ions, molecules) interacting with each other. If any of them oscillate, the oscillations spread to the surrounding medium, and energy is transmitted. If the surface particles of the body vibrate coherently (together with the same phase), the wave propagates through the medium as described in paragraph 1.2.2.1. In this way, a coherent wave arises in the medium. However, if the particles vibrate randomly (uncoordinated), there is also energy transfer, but the disturbance propagation process is not a coherent wave. Such transfer of energy from higher energy particles to lower energy particles is called thermal transfer. Since the energy of the particles is related to the temperature at a given location of the substance, we observe the transfer of energy from the higher temperature locations to the lower temperature locations. This phenomenon is called heat conduction in the substance. The interaction of particles in different substances varies, and therefore some substances conduct heat better (thermal conductors), others worse (thermal insulators). The power density is proportional to the temperature gradient.

where the * λ*is the specific thermal conductivity of the substance.

High thermal conductivity is typical for substances consisting of an ordered structure of small particles, especially metals. In contrast, disordered structures and substances consisting of large molecules (glass, plastics, etc.) are mostly thermal insulators. E.g., pure metal is a good conductor of heat, but steel, which contains many different impurities and structure failures, has a significantly lower thermal conductivity. As an example, compare an aluminum spoon and a stainless-steel spoon immersed in a hot tea - the aluminum end is soon hot, but the stainless-steel end remains cool.

#### 1.2.3 Oscillation’s propagation in a long LC chain

The electrical system analogous to the long chain of bound particles is a chain consisting of a homogeneous series of LC segments. The harmonic voltage _{0} = _{0} sin * ω t*is connected to the input, Figure 5.

If a chain is very long (theoretically infinite), its properties do not change when we add another segment to it. If the input impedance is ** Z**, it does not change by adding an LC segment to the beginning, Figure 6.

Then we get

Complex voltage transfer of the ** Z**loaded segment is then.

If we define a cut-off angular frequency of _{m} and characteristic resistance of _{0} as,

then, we can express the impedance of the chain and complex voltage transfer as.

If * ω* <

ω

_{m}(low-frequency oscillations), the absolute value of voltage transfer and voltage phase shift on one segment are.

The voltage of the * n-*th segment in chain is then.

It follows from this result that the voltage amplitude along the chain remains constant and the phase shift gradually increases: * φ* =

_{n}

*. Assuming harmonic excitation at the input,*n φ

_{k}

*propagates along the chain.*the harmonic undamped wave

If we know the segment length * d*and we express the phase difference as

*=*φ

_{k}

*Δ*ω

*, then, we can determine the velocity of the phase propagation (phase velocity of the wave):*t

For very low frequencies of excitation * ω* <<

ω

_{m}(taken an approximation arctan

*≈*x

*for*x

*<< 1), the velocity is constant and equals*x

For frequencies * ω* →

ω

_{m}the velocity

*→*v

Normalized dependency of velocity * v*versus angular frequency

*is in Figure 7.*ω

Low frequency oscillations having the angular frequency of * ω* <<

ω

_{m}propagate along the chain as an

*with constant velocity of*undamped harmonic wave

v

_{0}, for

*→*ω

ω

_{m}the velocity decreases to the

v

_{m}.

At * ω* >

ω

_{m}we get.

Absolute value of voltage transfer is

The impedance has inductive character (reactive character). Voltage transfer is less than one, which means that the voltage along the chain exponentially decreases.

Since the complex voltage transfer is a negative real number, adjacent segments oscillate with the opposite phase, but the phase constant of the segment does not change. It is, therefore, a * damped standing wave*.

When a harmonic voltage source is connected to the chain, its complex power supplied to the chain is.

For frequencies * ω* <

ω

_{m}the complex power of source consists of the active power part of

*and the reactive part of*P

Q.

Thus, the chain transfers the power (permanently draws energy from the source).

For frequencies * ω* >

ω

_{m}.

The chain is unable to transmit energy. The power source is only reactive, i.e., the mean value of the energy supplied by the source into the chain is zero.

We can see from these two cases that the oscillating events of mechanical and electrical are analogous, and thus there is electro-mechanical duality. One of the practical possibilities of using this feature is to model mechanical systems using electrical circuits. In the field of biomedicine, we can consider an example of blood flow modeling in blood vessels using electrical circuits [1].