Quantum-Classical Fragmentation with the Effective Fragment Molecular Orbital Method

We present the quantum effective fragment molecular orbital (Q-EFMO) method, a hybrid framework that evaluates selected molecular-fragment correlation energies with the variational quantum eigensolver and assembles them through the size-consistent EFMO energy expression. A Hartree-Fock EFMO calculation supplies monomer references, many-body polarization, and distant-pair EFP interactions; independent VQE/UCCSD calculations provide correlation increments for monomers and selected near-field dimers. Consequently, the maximum quantum register is determined by the largest active fragment or fragment pair rather than by total system size. For a three-layer LiH benchmark in STO-3G, a full 2.5 to 3.5 Angstrom separation scan shows that two orbital-reduction schemes approach the full CCSD reference monotonically. The more aggressive scheme remains below 1 kcal/mol for separations of 2.6 Angstrom and greater and reaches 0.10 kcal/mol at 3.5 Angstrom, using 6 qubits for monomers and 12 for the largest dimers, compared with 14 qubits for the frozen-core dimer baseline. An implementation-specific cost proxy decreases from 1000 to 40, a factor of 25. These noise-free results establish the equations, software path, and resource scaling; broader chemical validation, larger bases, and hardware tests remain necessary.