Mitsubishi Materials Quantum Simulation for Copper Alloy Design

The Problem

Electric vehicle motors and power electronics demand copper alloys that simultaneously achieve high electrical conductivity and mechanical strength, properties that classical metallurgy trades off against each other. Transition metal dopants (cobalt, nickel, iron) introduce precipitate phases that strengthen the copper matrix, but predicting which dopant concentration and crystal phase is most stable requires accurate modeling of d-orbital electron correlations. Standard density functional theory (DFT) handles weakly correlated systems well but systematically fails for transition metal d-electrons, where the Hubbard U correction (DFT+U) introduces a tunable parameter with no first-principles basis. Dynamical mean-field theory (DMFT) is more accurate but scales poorly with system size.

Mitsubishi Materials’ advanced materials division identified this as a candidate for quantum advantage: the active space consisting of Cu 3d and Co 3d orbitals in a Cu-Co cluster is small enough (roughly 10-14 orbitals) to fit current quantum hardware while being precisely the regime where classical methods are least reliable.

VQE Setup for Cu-Co Electronic Structure

The workflow used PySCF to generate the molecular integrals for a Cu6Co cluster representing a Co dopant in an FCC copper lattice. Qiskit Nature transformed the second-quantized Hamiltonian via Jordan-Wigner encoding into a qubit operator, then ran a UCCSD (Unitary Coupled Cluster Singles and Doubles) ansatz on IBM Heron 133Q hardware. An active space of 10 electrons in 10 orbitals was selected, focusing on the d-orbital manifold of the Co atom and its nearest Cu neighbors.

from pyscf import gto, scf, mcscf
from qiskit_nature.second_q.drivers import PySCFDriver
from qiskit_nature.second_q.mappers import JordanWignerMapper
from qiskit_nature.second_q.algorithms import VQEUCCFactory
from qiskit_algorithms import VQE
from qiskit_algorithms.optimizers import COBYLA
from qiskit.primitives import Estimator

# Define Cu6Co cluster (simplified; real calculation uses periodic boundary conditions)
driver = PySCFDriver(
    atom="""
    Cu  0.000  0.000  0.000;
    Cu  1.278  1.278  0.000;
    Cu  0.000  1.278  1.278;
    Cu  1.278  0.000  1.278;
    Co  2.556  1.278  1.278;
    Cu  1.278  2.556  1.278
    """,
    basis="def2-svp",
    charge=0,
    spin=1,
    unit="Angstrom"
)

problem = driver.run()

# Reduce to active space: 10 electrons in 10 orbitals (d-manifold)
problem.active_space = (10, 10)

mapper = JordanWignerMapper()
qubit_op = mapper.map(problem.second_q_ops()[0])

print(f"Active space qubits required: {qubit_op.num_qubits}")
print(f"Hamiltonian terms: {len(qubit_op)}")

# VQE with UCCSD ansatz
vqe_factory = VQEUCCFactory(
    estimator=Estimator(),
    ansatz_type="UCCSD",
    optimizer=COBYLA(maxiter=500)
)

result = vqe_factory.get_solver(problem, mapper).solve(problem)
print(f"VQE ground state energy: {result.total_energies[0]:.6f} Ha")
print(f"Correlation energy captured: {result.electronic_energies[0] - result.hartree_fock_energy:.6f} Ha")

Phase Stability Comparison

At 0.3 atomic percent cobalt, VQE predicted that the Cu3Co L12-ordered phase is more stable than the disordered FCC solid solution by 0.14 eV per formula unit. Standard DFT with the PBE functional predicted the disordered phase as marginally more stable (+0.03 eV), while DFT+U (U=3.5 eV) gave qualitatively correct ordering but with large uncertainty in the U parameter. The VQE result aligns with experimental observations of Cu-Co precipitation hardening at this composition range, which DFT consistently mischaracterizes.

The conductivity projection follows from the reduced scattering cross-section of ordered Co precipitates versus disordered dopant atoms. Phonon scattering from L12 precipitates is geometrically coherent with the copper matrix, reducing electron-defect scattering. The 12% conductivity improvement and 95% strength retention are projected from the VQE-informed phase diagram combined with classical transport models.

Path to Industrial Application

Mitsubishi Materials’ near-term validation plan involves synthesizing Cu-Co alloys at 0.1, 0.3, and 0.5 at% Co, measuring conductivity and tensile strength, and comparing against both DFT and VQE predictions. The goal is to establish which quantum chemistry level of theory is necessary for reliable alloy screening, a finding that would justify scaling quantum simulation workflows to IBM Heron’s full 133-qubit capacity for larger active spaces. Longer-term, the methodology extends to Cu-Ni and Cu-Fe systems relevant to Mitsubishi’s transformer wire and connector product lines.