This repository contains MATLAB scripts, example files, and license terms for the LDV model.

Calculation workflow

All calculations were conducted using MATLAB R2024a. Please note that, due to our data policy, no third-party data are included in this repository. Therefore, the MATLAB calculations cannot be executed directly using the provided GCAM_ssp_SSP2_RCP2.6_W_PHEV_2025_exNCB_PUBLIC.xlsx file as is.

1. GCAM calculation

GCAM is executed for a selected scenario. The standard calculations in this repository were performed under the SSP2-RCP2.6 scenario, with the LDV vehicle lifetime modified to follow a Weibull distribution and with PHEVs incorporated into to the model.

2. Preparation of input files

The GCAM calculation results are exported to Excel files. Please note that the provided file GCAM_ssp_SSP2_RCP2.6_W_PHEV_2025_exNCB_PUBLIC.xlsx does not include the original GCAM input data in the UCD_base, load_factor_raw, and share-weight sheets, which contain the energy intensities, load factors, and share weights of each technology. The missing data can be derived from transportation_UCD_CORE.xml in GCAM.

3. Conversion to annual data

Mineral_GCAM_read.m reads the input Excel file (.xlsx), converts the data from 5-year intervals to annual data, and outputs a .mat file. The input file must be placed in the same directory as Mineral_GCAM_read.m.

4. LDV model calculation

Mineral_GCAM_calc.m reads the .mat file. After specifying calculation conditions through the GUI, it iteratively calculates critical mineral prices, LIB costs, LDV technology choice, and net mineral demand. The calculation results are appended to base.xlsx and exported as an .xlsx file. Please also note that the LIB mineral-intensity scenarios corresponding to the LIB cathode chemistry scenarios (HiNMC and LFP) are not included in this repository.

Files section

• Mineral_GCAM_read.m
  Reads the Excel input file, converts 5-year data to annual data, and outputs a .mat file.

• Mineral_GCAM_calc.m
  Reads the .mat file and performs the iterative LDV model calculation.

• deflator_conv.m
  Applies an internal price deflator function to convert prices between different dollar-year bases.

• GCAM_ssp_SSP2_RCP2.6_W_PHEV_2025_exNCB_PUBLIC.xlsx
  Example/model input file prepared from GCAM outputs. The original GCAM input data are not included.

• output_base.xlsx
  Template/output file to which the calculation results are appended.

Notes on reproducibility

This repository contains the MATLAB-side calculation workflow and related input/output files, but does not include third-party data. As a result, the full end-to-end workflow cannot be reproduced using this repository alone.

In particular, the following items are not included:

• the GCAM input data referenced in GCAM_ssp_SSP2_RCP2.6_W_PHEV_2025_exNCB_PUBLIC.xlsx The missing data can be obtained from transportation_UCD_CORE.xml in GCAM v7.0 or from GCAM query results under standard conditions. PHEV data were estimated based on the parameters of HEV and BEV. The calculation can be reproduced by supplementing the missing data in the corresponding Excel sheets.
• the LIB mineral-intensity scenarios corresponding to the HiNMC and LFP cathode chemistry scenarios. These scenarios can be obtained from Xu et al. (2020) and Li et al. (2024), respectively. Note that additional modifications to Mineral_GCAM_calc.m are required to incorporate these scenarios into the calculations.