Coupled-Cluster Treatment of Complex Open-Shell Systems: The Case of Single-Molecule Magnets

Maristella Alessio,Garrette Pauley Paran,Cansu Utku,Andreas Grüneis,Thomas Christian Jagau

DOI: https://doi.org/10.1039/d4cp01129e

IF: 3.3

2024-05-28

Physical Chemistry Chemical Physics

Abstract:We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

chemistry, physical,physics, atomic, molecular & chemical

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to calculate spin - state energies and spin - related properties efficiently and accurately in multi - configuration electronic systems, especially in transition - metal complexes and single - molecule magnets (SMMs). Specifically, the paper evaluates the performance of two computationally less - expensive coupled - cluster methods - the second - order approximate coupled - cluster singles and doubles excitation (CC2) method and the equation - of - motion coupled - cluster singles and doubles excitation based on projection embedding (EOM - CCSD - in - DFT) method - in describing these complex systems.

### Background of the Paper
- **Research Field**: The application of quantum - chemical methods in studying the electronic states and properties of transition - metal complexes and their related materials, especially in fields such as catalysis, battery design, and molecular magnetism.
- **Challenge**: Describing open - shell electronic states is usually more challenging than closed - shell electronic states, and standard single - reference methods such as density functional theory (DFT) are often qualitatively wrong. Traditionally, multi - reference methods (such as CASPT2 and NEVPT2) are used to approximate the exact multi - configuration wave functions of these systems, and magnetic properties are extracted through phenomenological spin Hamiltonians.
- **Existing Methods**: The equation - of - motion coupled - cluster (EOM - CC) method has been used to deal with transition - metal SMMs, but its computational cost increases sharply with the size of the system and is not suitable for large - molecule magnets.

### Research Objectives
- **Evaluate CC2 and EOM - CCSD - in - DFT**: Evaluate the performance of these two methods in calculating spin - state energies and spin - related properties (such as spin - orbit coupling constants, magnetization intensity, and magnetic susceptibility) of transition - metal complexes.
- **Test Set**: Selected transition - metal complexes with different electronic configurations as a test set, including two ferric hexahydrate complexes (\[Fe(H_2O)_6\]^{3+} and \[Fe(H_2O)_6\]^{2+}), and a cobalt - based single - center SMM (Co(C(SiMe_2ONaph)_3)_2) and its simplified model (Co(C(SiH_3)_3)_2).

### Main Findings
- **CC2 Method**: It is more accurate in describing excited states, and the difference from the EOM - CCSD excitation energy is within 0.05 eV.
- **EOM - CCSD - in - DFT Method**: It performs excellently in describing orbital angular momentum and spin - orbit coupling. For Co(II) molecule magnets, the spin - reversal energy barrier, temperature - dependent and field - dependent magnetization intensity, and magnetic susceptibility calculated using EOM - CCSD - in - DFT are highly consistent with experimental values within spectral accuracy.

### Conclusions
- **Efficiency and Accuracy**: The CC2 method is efficient and accurate in calculating the state energies of multi - configuration, open - shell systems, while the EOM - CCSD - in - DFT method has a higher cost - performance ratio in calculating spin - orbit coupling and magnetic properties.
- **Application Prospects**: These methods provide reliable and efficient tools for studying large and complex molecule magnets, which are helpful for further understanding and designing single - molecule magnets with specific behaviors.

Through these studies, the paper provides important guidance and inspiration for the future development and application of quantum - chemical methods.