Complex Numbers for Physics

Complex numbers, Euler's formula, modulus, and phase applied to wave physics and harmonic analysis.

What it is

A complex number is a mathematical entity of the form $z = x + i y$, where $x = \text{Re}(z)$ and $y = \text{Im}(z)$ are real numbers and $i$ is the imaginary unit defined by the fundamental property $i^2 = -1$. Geometrically, complex numbers extend the one-dimensional real line into the two-dimensional complex plane (Argand plane), providing points $(x,y)$ with an algebraic structure where addition, subtraction, multiplication, and division are universally defined without algebraic limitations.

By Euler's formula, $e^{i\theta} = \cos\theta + i\sin\theta$, any non-zero complex number can be expressed in polar or exponential form: $z = r e^{i\theta}$, where $r = |z| = \sqrt{x^2 + y^2}$ is the modulus (amplitude) and $\theta = \arg(z) = \arctan(y/x)$ is the argument (phase angle). This exponential representation transforms multiplications and divisions into simple angle additions and modulus scalings, providing an indispensable analytical framework for oscillatory and wave physics.

Why it appears in physics

Complex numbers appear in physics not as artificial contrivances, but as the most concise and powerful tool to package two physically coupled degrees of freedom: amplitude and phase. Instead of manipulating cumbersome trigonometric identities when solving linear differential equations, complex variables linearize the time derivative operator: each temporal derivative of a harmonic state $e^{i\omega t}$ reduces to an algebraic multiplication by $i\omega$.

They emerge ubiquitously in simple and damped harmonic oscillators, electromagnetic theory and wave optics (complex refractive index for absorbing media, wave polarization), alternating current circuits (complex impedance), acoustics, and fundamentally in quantum mechanics, where the wave function $\psi(x,t)$ inherently populates a complex Hilbert space.

Visual or geometric idea

Visualize the Cartesian plane where the horizontal axis represents the real part and the vertical axis the imaginary part. A complex number $z$ is a position vector directed from $(0,0)$ to $(x,y)$. The addition of two complex numbers follows the exact parallelogram law of vector addition.

However, unlike ordinary two-dimensional vectors, complex numbers possess an intrinsic multiplication rule governed by rotation and scaling: multiplying any complex quantity by $i = e^{i\pi/2}$ corresponds to rotating the vector counterclockwise by 90 degrees without altering its magnitude. Multiplying by $r e^{i\theta}$ scales the length by $r$ and rotates it by $\theta$. Consequently, a harmonic oscillation is envisioned as a constant-magnitude vector rotating with angular frequency $\omega$; the physically observed measurable variable is merely the projection of this vector onto the horizontal real axis.

Minimal symbolic language

• Cartesian (algebraic) form: $z = x + i y$, with $x = \text{Re}(z) \in \mathbb{R}$ and $y = \text{Im}(z) \in \mathbb{R}$.
• Imaginary unit: $i = \sqrt{-1}$, satisfying $i^2 = -1$, $i^3 = -i$, $i^4 = 1$.
• Modulus (magnitude): $|z| = r = \sqrt{x^2 + y^2} = \sqrt{z z^*}$.
• Argument (phase): $\theta = \arg(z) = \text{atan2}(y, x)$.
• Complex conjugate: $z^* = x - i y = r e^{-i\theta}$; reflection across the real axis.
• Euler's formula: $e^{i\theta} = \cos\theta + i\sin\theta$, featuring the famous identity $e^{i\pi} + 1 = 0$.
• Multiplicative inverse: $\frac{1}{z} = \frac{z^*}{|z|^2} = \frac{x - i y}{x^2 + y^2}$.

What the result means

When analyzing a physical problem with complex variables, the derived solution is typically expressed as $Z(t) = A e^{i(\omega t + \phi)}$. The physical interpretation is retrieved by taking the real part: $\text{Re}[Z(t)] = A \cos(\omega t + \phi)$, which represents the measurable displacement, electric field, current, or voltage in an experimental setup.

The modulus $A$ determines the maximum amplitude reached by the system, while the argument $\phi$ quantifies the temporal lead or lag relative to a reference signal. If the frequency acquires an imaginary component $\omega = \omega_0 + i\gamma$, the term becomes $e^{-\gamma t} e^{i\omega_0 t}$, indicating that the oscillation frequency is $\omega_0$ while its amplitude undergoes exponential damping with decay rate $\gamma$.

When not to use it

Do not use complex superposition linear techniques in systems governed by nonlinear equations (such as anharmonic oscillators or nonlinear optics with higher-order susceptibility), because the real part of a product does not equal the product of the real parts: $\text{Re}(z_1 z_2) \neq \text{Re}(z_1) \text{Re}(z_2)$.

Furthermore, take care when computing quadratic energetic observables: instantaneous electric power cannot be calculated by simply multiplying complex voltage and current signals. Instead, use the time-averaged formula $P_{\text{avg}} = \frac{1}{2} \text{Re}[V I^*]$ or compute products strictly after extracting real parts.

Operational summary

Transform real harmonic observables $A \cos(\omega t + \phi)$ into their complex representation $A e^{i\phi} e^{i\omega t}$. Carry out all linear differential or algebraic steps in the complex domain by replacing $\frac{d}{dt} \to i\omega$ and $\int dt \to \frac{1}{i\omega}$. Once the general solution is obtained, take the real part to recover physical trajectories or evaluate $|z|^2$ for probability densities and energy fluxes.

Procedure

Before calculating

Verify that the system is described by linear differential equations or harmonic signals where complex algebraic representation is valid.

Determine whether the given values are in Cartesian ($x + iy$) or polar ($r e^{i\theta}$) form, and establish a uniform time sign convention (typically $e^{i\omega t}$ in classical physics and electrical engineering).

Step 1: identify the mathematical object

Classify all quantities: distinguish between real amplitudes, complex impedances, wave functions, and phasors.

Identify the angular frequency $\omega$, the peak amplitude $A$, and the initial phase shift $\phi$.

Step 2: choose the operation

For addition or subtraction of signals at the same frequency, use Cartesian components: add real parts together and imaginary parts together: $(x_1 + x_2) + i(y_1 + y_2)$.

For multiplication, division, or powers (such as applying generalized Ohm's law $\tilde{V} = \tilde{Z}\tilde{I}$ or De Moivre's formula), convert terms to polar form: $z_1 z_2 = r_1 r_2 e^{i(\theta_1 + \theta_2)}$.

Step 3: preserve units and signs

Ensure that phase angles $\theta$ are evaluated in radians when substituting into exponential or trigonometric expressions.

Remember that taking the complex conjugate strictly reverses the sign of the imaginary component ($i \to -i$), which inverts the relative phase orientation.

Step 4: calculate

Perform complex operations, eliminating imaginary units from denominators by multiplying numerator and denominator by the complex conjugate: $\frac{w}{z} = \frac{w z^*}{|z|^2}$.

When solving linear differential equations with constant coefficients, replace time derivatives $\frac{d}{dt}$ with $i\omega$ and algebraically isolate the unknown complex amplitude.

Step 5: interpret the result

Extract the modulus $|z| = \sqrt{x^2 + y^2}$ to determine the observable physical amplitude.

Compute the net phase $\phi = \text{atan2}(y, x)$ and reconstruct the physical time-domain signal as $s(t) = |z| \cos(\omega t + \phi)$.

Self-check checklist

• Did you correctly account for the Cartesian quadrant when computing $\theta = \text{atan2}(y, x)$?
• Did you multiply numerator and denominator by the conjugate of the denominator when dividing?
• Did you extract the real part at the end to state the physically observable quantity?
• Do the physical dimensions of the modulus match the expected units of the problem?

Guided example

Physical situation

A series RLC circuit driven by an AC voltage $V(t) = V_0 \cos(\omega t)$ with $V_0 = 100\text{ V}$ and angular frequency $\omega = 1000\text{ rad/s}$ consists of a resistor $R = 40\;\Omega$, an inductor $L = 0.05\text{ H}$, and a capacitor $C = 25\;\mu\text{F}$. We wish to find the steady-state current $I(t)$.

Mathematical translation

Represent the input voltage by the complex phasor $\tilde{V} = V_0 e^{i 0} = 100\text{ V}$. The total complex impedance of the series circuit is $\tilde{Z} = R + i\left(\omega L - \frac{1}{\omega C}\right)$. The complex current phasor satisfies generalized Ohm's law: $\tilde{I} = \frac{\tilde{V}}{\tilde{Z}}$.

Step-by-step calculation

1. Compute the inductive reactance: $X_L = \omega L = 1000 \cdot 0.05 = 50\;\Omega$.
2. Compute the capacitive reactance: $X_C = \frac{1}{\omega C} = \frac{1}{1000 \cdot 25 \times 10^{-6}} = 40\;\Omega$.
3. Net reactance: $X = X_L - X_C = 50 - 40 = 10\;\Omega$.
4. Total complex impedance: $\tilde{Z} = 40 + 10 i\;\Omega$.
5. Impedance magnitude: $|\tilde{Z}| = \sqrt{40^2 + 10^2} = \sqrt{1700} \approx 41.23\;\Omega$.
6. Impedance phase angle: $\theta_Z = \arctan(10/40) = \arctan(0.25) \approx 0.245\text{ rad} \approx 14.04^\circ$.
7. Current amplitude: $I_0 = \frac{V_0}{|\tilde{Z}|} = \frac{100}{41.23} \approx 2.425\text{ A}$.
8. Current phase: $\phi_I = 0 - \theta_Z = -0.245\text{ rad}$.

Reading the result

The physical real current is $I(t) = 2.425 \cos(1000 t - 0.245)\text{ A}$. The negative phase confirms that the current lags the voltage by $14.04^\circ$, characteristic of an inductive circuit where $X_L > X_C$.

Common error avoided

We avoided summing impedances as simple scalar values ($40 + 50 + 40 = 130\;\Omega$), which would fail to account for the $\pm 90^\circ$ relative phase shifts introduced by reactive components.

Inverse example

Given result

In a driven mechanical oscillator at steady state, the measured complex displacement is $\tilde{x} = -i \frac{F_0}{m \gamma \omega_0}$.

What physical question it can answer

This expression reveals whether the driving frequency matches the resonance condition: because $-i = e^{-i\pi/2}$ introduces a phase lag of exactly $90^\circ$ relative to the driving force, the system is in exact velocity resonance, allowing one to deduce the damping factor directly from the measured response amplitude: $\gamma = \frac{F_0}{m \omega_0 |\tilde{x}|}$.

Common errors

Error 1: ignoring the quadrant when evaluating argument via arctan

It seems reasonable to compute $\theta = \arctan(y/x)$ directly using a calculator function for any complex value. It fails because standard arctangent functions restrict outputs to the interval $(-\pi/2, \pi/2)$, ignoring negative real parts and erroneously placing second- and third-quadrant numbers into the fourth or first quadrants. It is detected by plotting $(x,y)$ on the Argand plane and checking the sign of $x$. It is corrected by using $\text{atan2}(y, x)$ or adding $\pi$ radians whenever $x < 0$. Mini-example: for $z = -1 + i$, $\arctan(1/(-1)) = -45^\circ$, whereas the actual physical argument is $-45^\circ + 180^\circ = 135^\circ$.

Error 2: distributing complex conjugation incorrectly across products and quotients

It seems reasonable to assume conjugation only affects isolated terms with explicit $i$, forgetting to conjugate variables within compound fractions or products. It fails because $(z_1 z_2)^* = z_1^* z_2^*$ and $(z_1/z_2)^* = z_1^*/z_2^*$; failing to conjugate either component severely corrupts the resulting phase angle. It is detected by ensuring every complex variable in the algebraic expression was systematically replaced by its conjugate. It is corrected by applying $(A B)^* = A^* B^*$ consistently. Mini-example: if $z = \frac{1}{1 + i\omega}$, its conjugate is $z^* = \frac{1}{1 - i\omega}$, not $\frac{1}{1 + i\omega}$.

Error 3: evaluating nonlinear physical quantities directly with complex phasors

It seems reasonable to multiply complex phasors $\tilde{A} e^{i\omega t}$ and $\tilde{B} e^{i\omega t}$ to obtain instantaneous physical products such as power $p(t) = v(t) i(t)$. It fails because $\text{Re}(\tilde{A}) \cdot \text{Re}(\tilde{B}) \neq \text{Re}(\tilde{A} \tilde{B})$; the product of raw phasors yields double-frequency complex terms that do not reflect the physical real product. It is detected when calculated energy or power expressions retain an imaginary component. It is corrected by taking real parts before multiplying or using the time-averaged phasor relation $\frac{1}{2} \text{Re}(\tilde{A} \tilde{B}^*)$. Mini-example: for $v(t) = \cos(\omega t)$ and $i(t) = \cos(\omega t)$, the true product is $\cos^2(\omega t) = \frac{1}{2}(1 + \cos(2\omega t))$, whereas $\text{Re}(e^{i\omega t} e^{i\omega t}) = \cos(2\omega t)$, omitting the essential constant $1/2$ dissipation term.

Error 4: confusing temporal sign conventions across disciplines

It seems reasonable to interchange $e^{i\omega t}$ and $e^{-i\omega t}$ arbitrarily, assuming signs are trivial aesthetic choices. It fails because altering the temporal exponential sign reverses the physical meaning of phase lead vs. lag, and turns a rightward traveling wave ($e^{i(kx - \omega t)}$) into a leftward wave unless the spatial wavevector is concurrently inverted. It is detected by verifying whether a delayed signal $t - t_0$ produces the expected negative phase shift. It is corrected by stating the chosen convention at the outset and strictly preserving it throughout. Mini-example: under $e^{i\omega t}$, an inductive reactance is $+i\omega L$; under $e^{-i\omega t}$, the identical physical inductor requires $-i\omega L$.

Error 5: dividing complex numbers without rationalizing via the conjugate

It seems reasonable to attempt dividing $\frac{x_1 + i y_1}{x_2 + i y_2}$ by dividing components term by term: $\frac{x_1}{x_2} + i \frac{y_1}{y_2}$. It fails because division over the field of complex numbers requires converting the denominator to a purely real scalar, which mixes Cartesian parts in cross products. It is detected whenever an imaginary unit lingers in the denominator of a final result. It is corrected by multiplying numerator and denominator by the conjugate $x_2 - i y_2$, producing the real denominator $x_2^2 + y_2^2$. Mini-example: $\frac{1}{1+i} = \frac{1-i}{(1+i)(1-i)} = \frac{1-i}{2} = 0.5 - 0.5i$, which is completely different from $1/1 + i/1 = 1 + i$.

Quick detection rule

If a denominator retains an imaginary unit $i$, if an energy or power dissipation calculation yields an imaginary component, or if a polar angle does not match the Cartesian quadrant of the Argand plane, there is an algebraic error in the complex formulation.

Physical applications

1. Alternating current circuits and complex impedance

It appears when analyzing electrical networks driven by sinusoidal power sources containing resistors, inductors, and capacitors. It is calculated by assigning each element its complex impedance: $Z_R = R$, $Z_L = i\omega L$, and $Z_C = \frac{1}{i\omega C}$, converting integro-differential loop equations into linear algebraic Ohm and Kirchhoff laws. The result provides the amplitude ratio and phase difference between driving voltage and resulting current. The dangerous error is adding reactive terms algebraically as real scalars rather than maintaining orthogonal phase quadrature. As a final check, the student must explain why an ideal inductor advances voltage by $90^\circ$ while a capacitor delays it.

2. Damped and driven harmonic mechanical oscillators

It appears when modeling mechanical vibrations subjected to viscous friction and external periodic driving forces $F(t) = F_0 e^{i\omega t}$. It is calculated by postulating a steady-state trial solution $x(t) = \tilde{x}_0 e^{i\omega t}$, transforming the second-order differential equation into an algebraic equation for the complex amplitude $\tilde{x}_0$. The result yields the resonance response curve and the frequency-dependent phase lag between driver and displacement. The dangerous error is discarding the imaginary part of complex displacement prior to differentiating to find velocity. As a final check, the student must verify that at velocity resonance the mechanical velocity is precisely in phase with the applied force.

3. Electromagnetic waves and absorbing dielectric media

It appears when examining the propagation of light and radiofrequency waves through lossy conductors or dissipative dielectrics. It is calculated by introducing a complex refractive index $\tilde{n} = n + i\kappa$, where $n$ dictates phase velocity and $\kappa$ is the extinction coefficient. The result produces a plane wave factor $e^{i\tilde{k} x} = e^{-\kappa \frac{\omega}{c} x} e^{i n \frac{\omega}{c} x}$, capturing both spatial phase propagation and exponential attenuation according to the Beer-Lambert law. The dangerous error is confusing the field amplitude attenuation factor with the optical power attenuation factor, which is twice as large due to squaring. As a final check, the student must relate the electrical conductivity of a metal to the imaginary part of its refractive index.

4. Wave optics and light polarization states

It appears when characterizing the polarization state (linear, circular, or elliptical) of transverse electromagnetic waves via complex Jones vectors $\begin{pmatrix} E_{0x} e^{i\phi_x} \\ E_{0y} e^{i\phi_y} \end{pmatrix}$. It is calculated by extracting the phase difference $\Delta\phi = \phi_y - \phi_x$ between orthogonal electric field vectors. The result defines the geometric trajectory traced by the electric field vector over one optical period. The dangerous error is assuming that equal orthogonal field amplitudes guarantee circular polarization without confirming that the relative phase shift is strictly $\pm 90^\circ$. As a final check, the student must demonstrate why zero phase shift yields an inclined linear polarization state.

5. Quantum mechanical wavefunctions and state evolution

It appears in the time-dependent Schrödinger equation $i\hbar \frac{\partial \psi}{\partial t} = \hat{H}\psi$. It is calculated by solving for stationary eigenstates $\psi(x,t) = \phi(x) e^{-i E t/\hbar}$, where the complex phase precesses at angular frequency $\omega = E/\hbar$. The result yields the observable probability density via the complex conjugate product $\rho(x) = \psi^* \psi = |\psi|^2$. The dangerous error is treating $\psi$ itself as a directly measurable classical voltage or position, rather than recognizing that only conjugate bilinear products $\psi^* \psi$ and operator expectation values represent measurable quantities. As a final check, the student must prove that total probability $\int |\psi|^2 dx = 1$ is conserved over time due to the hermiticity of the Hamiltonian operator.

Practice

Level 1: recognition

Convert the following complex numbers into polar form $r e^{i\theta}$ and sketch them on the complex Argand plane: a) $z_1 = 1 + i$ b) $z_2 = -3$ c) $z_3 = -2i$ d) $z_4 = -\sqrt{3} + i$

Exercise: For $z_4$, confirm that the angle lies in the second quadrant by computing its principal argument in radians.

Level 2: direct application

Evaluate the following complex algebraic expressions, reducing each to standard Cartesian form $x + i y$: a) $(3 + 4i)(1 - 2i)$ b) $\frac{5 + 5i}{3 - 4i}$ c) $(1 + i)^8$ applying Euler's formula and De Moivre's theorem.

Exercise: Determine the modulus and complex conjugate of $\frac{2 + i}{1 - 3i}$.

Level 3: choosing the tool

Two coherent light waves of identical frequency arrive at a photodetector with amplitudes $A_1$ and $A_2$ and relative phase difference $\Delta\phi$.

Exercise: Formulate their superposition as the sum of two complex phasors $\tilde{E}_1 = A_1$ and $\tilde{E}_2 = A_2 e^{i\Delta\phi}$. Prove algebraically that the resultant intensity $|\tilde{E}_1 + \tilde{E}_2|^2$ yields $A_1^2 + A_2^2 + 2 A_1 A_2 \cos(\Delta\phi)$ without tedious trigonometric expansions.

Level 4: physical interpretation

In a lossy plasma or conducting medium, the electromagnetic dispersion relation yields a complex wave vector $k = \beta + i \alpha$.

Exercise: Insert $k$ into the spatial factor of a plane wave $e^{i k x}$ and explain physically how $\beta$ and $\alpha$ dictate the phase velocity and skin depth (penetration length) in the medium.

Level 5: intentional error

A student attempts to determine the average power dissipated in an AC load with voltage $\tilde{V} = 100 e^{i 0}$ and current $\tilde{I} = 5 e^{-i \pi/3}$ by computing: $P = \text{Re}(\tilde{V} \cdot \tilde{I}) = 100 \cdot 5 \cos(-\pi/3) = 250\text{ W}$.

Exercise: Identify why the rigorous complex power formula requires the complex conjugate of the current $\frac{1}{2} \text{Re}(\tilde{V} \tilde{I}^*)$, correct the missing $1/2$ factor for peak values, and explain what would happen if the voltage had a non-zero initial phase.