Performance of the coupled cluster and DFT methods for through-space magnetic interactions of nitroxide dimer

Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids

Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue

The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S₁ and dark S₂ surfaces near the S₁ minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S₁/S₂ state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.

Assessing the Calculation of Exchange Coupling Constants and Spin Crossover Gaps Using the Approximate Projection Model To Improve Density Functional Calculations

This work evaluates the quality of exchange coupling constant and spin crossover gap calculations using density functional theory corrected by the approximate projection model. Results show that improvements using the approximate projection model range from modest to significant. This study demonstrates that, at least for the class of systems examined here, spin projection generally improves the quality of density functional theory calculations of J-coupling constants and spin crossover gaps. Furthermore, it is shown that spin projection can be important for both geometry optimization and energy evaluations. The approximate projection model provides an affordable and practical approach for effectively correcting spin-contamination errors in such calculations.

The Role of Orbital Symmetries in Enforcing Ferromagnetic Ground State in Mixed Radical Dimers

One of the first steps in designing ferromagnetic (FM) molecular materials of p-block radicals is the suppression of covalent radical-radical interactions that stabilize a diamagnetic ground state. In this contribution, we demonstrate that FM coupling between p-block radicals can be achieved by constructing mixed dimers from different radicals with differing symmetries of their singly occupied molecular orbitals. The applicability of this approach is demonstrated by studying magnetic interactions in organic radical dimers built from different derivatives of the well-known phenalenyl radical. The calculated enthalpies of dimerization for different homo- and heterodimers show that the formation of a mixed dimer with FM coupling is favored compared to the formation of homodimers with antiferromagenetic (AFM) coupling. We argue that cocrystallization of radicals with specifically tuned morphologies of their singly occupied molecular orbitals is a feasible and promising approach in designing new organic magnetic materials.

Synthesis of a labile sulfur-centred ligand, [S(H)C(PPh2S)2](-): structural diversity in lithium(i), zinc(ii) and nickel(ii) complexes

A high-yield synthesis of [Li{S(H)C(PPh2S)2}]2 [Li2·(3)2] was developed and this reagent was used in metathesis with ZnCl2 and NiCl2 to produce homoleptic complexes 4 and 5b in 85 and 93% yields, respectively. The solid-state structure of the octahedral complex [Zn{S(H)C(PPh2S)2}2] (4) reveals notable inequivalence between the Zn-S(C) and Zn-S(P) contacts (2.274(1) Å vs. 2.842(1) and 2.884(1) Å, respectively). Two structural isomers of the homoleptic complex [Ni{S(H)C(PPh2S)2}2] were isolated after prolonged crystallization processes. The octahedral green Ni(ii) isomer 5a exhibits the two monoprotonated ligands bonded in a tridentate (S,S',S'') mode to the Ni(ii) centre with three distinctly different Ni-S bond lengths (2.3487(8), 2.4500(9) and 2.5953(10) Å). By contrast, in the red-brown square-planar complex 5b the two ligands are S,S'-chelated to Ni(ii) (d(Ni-S) = 2.165(2) and 2.195(2) Å) with one pendant PPh2S group. DFT calculations revealed that the energetic difference between singlet and triplet state octahedral and square-planar isomers of the Ni(ii) complex is essentially indistinguishable. Consistently, VT and (31)P CP/MAS NMR spectroscopic investigations indicated that a mixture of isomers exists in solution at room temperature, while the singlet state square-planar isomer 5b becomes favoured at -40 °C.