M A T L A B T E C H

The Engineering Guide to Cell Chemistry and the SOC-OCV Curve

Decoding the relationship between state of charge, open-circuit voltage, and crystalline phase transitions for advanced powertrain modeling.

In the realm of electric vehicle engineering and battery management system (BMS) design, the State of Charge to Open-Circuit Voltage (SOC-OCV) curve is the foundational blueprint of a battery cell. It is the mathematical bridge that connects the abstract concept of “remaining capacity” to the physically measurable metric of “voltage.”

However, treating the OCV curve merely as a static 2D lookup table is a critical engineering oversight. The shape, slope, and hysteresis of this curve are directly dictated by the complex thermodynamic reactions and crystalline structural changes occurring at the atomic level within the cathode and anode. For engineers utilizing algorithmic environments to simulate Equivalent Circuit Models (ECM) or design Kalman Filter estimators, a deep understanding of why the curve behaves the way it does is paramount.

1. The Thermodynamics of Voltage: Why OCV Exists

Open-Circuit Voltage (OCV) is defined as the electrical potential difference between the positive and negative terminals of a battery when no external load is connected (i.e., zero current flow), and the cell has reached thermodynamic equilibrium. It is the battery in its absolute resting state.

The voltage of a cell is fundamentally a measure of chemical potential energy. It is governed by the change in Gibbs Free Energy ($\Delta G$) during the electrochemical intercalation process. The relationship between the free energy of the reaction and the resulting open-circuit voltage is expressed by the fundamental thermodynamic equation:

$$ \Delta G = -n \cdot F \cdot E_{OCV} $$
Where:
• $\Delta G$ = Change in Gibbs Free Energy (Joules/mol)
• $n$ = Number of electrons transferred per reaction (for Li-ion, $n = 1$)
• $F$ = Faraday’s Constant (96,485 Coulombs/mol)
• $E_{OCV}$ = Open-Circuit Voltage (Volts)

As the battery discharges, Lithium ions ($Li^+$) migrate from the anode to the cathode. The chemical potential of both electrodes changes continuously as the concentration of Lithium varies. The OCV at any given State of Charge (SOC) is directly derived from the Nernst Equation, which adapts the standard potential based on the activity (concentration) of the reactants and products.

2. Dissecting the SOC-OCV Curve

An SOC-OCV curve maps the resting voltage (Y-axis) against the percentage of remaining capacity (X-axis). While every battery chemistry has a unique curve, they all generally follow a specific macroscopic trend: voltage drops as SOC decreases. However, the rate of this drop—the slope—reveals the underlying phase mechanics of the cell.

Macroscopic OCV Comparison: Solid Solution vs. Two-Phase Coexistence
4.2V 3.8V 3.4V 3.0V 2.6V 0% 25% 50% 75% 100% State of Charge (SOC) Open Circuit Voltage (V) NMC Solid Solution Region LFP Two-Phase Coexistence Plateau

Figure 1: The distinct voltage responses of NMC (Solid Solution) and LFP (Two-Phase) chemistries. Notice the extreme flatness of the LFP curve between 15% and 85% SOC.

2.1 Solid Solution Chemistries (The Sloping Curve)

Cathode chemistries like Nickel Manganese Cobalt (NMC) and Nickel Cobalt Aluminum (NCA) exhibit a distinct, sloping OCV curve. As Lithium ions intercalate into the layered metal-oxide lattice of an NMC cathode, the crystal structure expands and contracts, but it generally maintains a single, continuous phase—a Solid Solution.

Because the concentration of Lithium in this single phase changes continuously, the chemical potential (and thus the voltage) changes continuously. This steady slope is highly advantageous for BMS engineering. A simple voltage measurement (after the cell has rested) can reliably pinpoint the exact SOC through a basic 1D lookup table.

2.2 Two-Phase Coexistence (The Flat Curve)

Lithium Iron Phosphate (LFP) presents a massive algorithmic challenge because its OCV curve is remarkably flat between 15% and 85% SOC. This flatness is a direct result of the Gibbs Phase Rule.

When Lithium intercalates into an LFP cathode, the material does not form a continuous solid solution. Instead, it transitions between two distinct crystalline phases:

  1. Lithium-poor phase: $FePO_4$ (Heterosite)
  2. Lithium-rich phase: $LiFePO_4$ (Triphylite)

According to thermodynamic principles, when two phases coexist in equilibrium at a constant temperature and pressure, the chemical potential remains constant. Therefore, as long as both $FePO_4$ and $LiFePO_4$ are present in the cathode, the voltage is pinned at approximately 3.2V to 3.3V. The voltage only drops sharply when the material is almost entirely converted to the Lithium-poor phase (near 0% SOC), and rises sharply when fully converted to the Lithium-rich phase (near 100% SOC).

The Algorithmic Challenge of LFP

For engineers developing control logic, the flat OCV curve of LFP means that a resting voltage of 3.28V might correspond to 40% SOC, or it might correspond to 60% SOC. Standard voltage-to-SOC lookup algorithms fail entirely. Accurate estimation requires advanced techniques like Coulomb Counting (Current Integration) tightly coupled with Extended Kalman Filters (EKF) to continuously correct the SOC based on dynamic internal resistance rather than static voltage.

3. Mathematical Modeling of the OCV Curve

To simulate a battery pack in software, the physical OCV data must be parameterized into a mathematical function. This allows the simulation environment to calculate the OCV continuously as the SOC state variable changes.

3.1 Polynomial Curve Fitting

For sloping chemistries like NMC, a high-order polynomial equation is often sufficient to map the curve. A 5th or 6th-order polynomial provides a good balance between computational efficiency and accuracy.

% MATLAB Pseudocode for OCV Polynomial Fitting
% OCV(SOC) = p1*SOC^6 + p2*SOC^5 + p3*SOC^4 + … + p7

soc_data = [0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0];
ocv_data = [3.0, 3.4, 3.55, 3.65, 3.75, 3.82, 3.9, 3.98, 4.05, 4.12, 4.2];

% Fit a 6th degree polynomial
p_coeffs = polyfit(soc_data, ocv_data, 6);

% Calculate OCV for a given SOC
current_soc = 0.65;
estimated_ocv = polyval(p_coeffs, current_soc);

3.2 Combined Model (The Plett Model)

Dr. Gregory Plett introduced a highly robust mathematical model that combines polynomial terms with exponential and logarithmic terms. This model is exceptionally good at capturing the sharp voltage elbows near 0% and 100% SOC, which simple polynomials often over-smooth.

$$ OCV(SOC) = K_0 – \frac{K_1}{SOC} – K_2 SOC + K_3 \ln(SOC) + K_4 \ln(1 – SOC) $$

In this equation, the $\ln(SOC)$ term forces the voltage to plunge near complete discharge, while the $\ln(1 – SOC)$ term captures the sharp rise during complete polarization at full charge. Engineers use parameter estimation toolboxes to find the optimal $K$ coefficients based on empirical pulse-discharge testing data.

Conclusion

The SOC-OCV curve is the definitive translation matrix between the thermodynamic reality of cell chemistry and the digital realm of powertrain simulation. By recognizing that the slope of NMC is dictated by solid solution intercalation, and the plateau of LFP is governed by two-phase coexistence, engineers can move beyond basic lookup tables. Implementing advanced polynomial fitting, hysteresis tracking, and temperature compensation within simulation models is the key to unlocking highly accurate range estimations and robust battery management systems.