Partitioning scheme for density functional calculations of extended systems

Area of Science:

• Computational Chemistry
• Quantum Chemistry
• Materials Science

Background:

• Accurate electronic structure calculations are crucial for understanding molecular systems.
• Current methods for large, nonperiodic systems are computationally expensive due to unfavorable scaling with system size.
• Treating extended systems often involves approximations at the boundaries, which can affect accuracy.

Purpose of the Study:

• To develop a more computationally efficient method for calculating the electronic structure of large, nonperiodic molecular systems.
• To demonstrate that calculations on smaller, overlapping subsystems can accurately represent the entire system.
• To reduce the computational cost associated with modern electronic structure methods.

Main Methods:

• Decomposition of the extended system into overlapping subsystems.
• Termination of subsystems with link (or pseudo) atoms to mitigate boundary effects.
• Validation of the method by comparing results to calculations on the entire system.
• Application to serine-water, lysine-water, and lysine dimer systems.

Main Results:

• The proposed subsystem method accurately reproduces the total energy, electronic density, and atomic forces of the entire system.
• The computational efficiency is significantly improved due to favorable scaling compared to traditional methods.
• The effects of link atoms are largely compensated for when subsystems are sufficiently large and appropriately chosen.
• The method shows potential for approximating quantum-classical calculations, minimizing quantum cluster boundary issues.

Conclusions:

• The subsystem-based approach offers a computationally efficient and accurate alternative for electronic structure calculations of extended systems.
• This method significantly reduces the computational burden without sacrificing accuracy for ground electronic states.
• Potential applications extend to quantum-classical simulations, offering a more seamless treatment of quantum-classical boundaries.