Detailed ab initio first-Principles study of the magnetic anisotropy in a family of trigonal pyramidal iron(II) pyrrolide complexes

A theoretical, computational, and conceptual framework for the interpretation and prediction of the magnetic anisotropy of transition metal complexes with orbitally degenerate or orbitally nearly degenerate ground states is explored. The treatment is based on complete active space self-consistent field (CASSCF) wave functions in conjunction with N-electron valence perturbation theory (NEVPT2) and quasidegenerate perturbation theory (QDPT) for treatment of magnetic field- and spin-dependent relativistic effects. The methodology is applied to a series of Fe(II) complexes in ligand fields of almost trigonal pyramidal symmetry as provided by several variants of the tris-pyrrolylmethyl amine ligand (tpa). These systems have recently attracted much attention as mononuclear single-molecule magnet (SMM) complexes. This study aims to establish how the ligand field can be fine tuned in order to maximize the magnetic anisotropy barrier. In trigonal ligand fields high-spin Fe(II) complexes adopt an orbitally degenerate ⁵E ground state with strong in-state spin–orbit coupling (SOC). We study the competing effects of SOC and the ⁵E⊗ε multimode Jahn–Teller effect as a function of the peripheral substituents on the tpa ligand. These subtle distortions were found to have a significant effect on the magnetic anisotropy. Using a rigorous treatment of all spin multiplets arising from the triplet and quintet states in the d⁶ configuration the parameters of the effective spin-Hamiltonian (SH) approach were predicted from first principles. Being based on a nonperturbative approach we investigate under which conditions the SH approach is valid and what terms need to be retained. It is demonstrated that already tiny geometric distortions observed in the crystal structures of four structurally and magnetically well-documented systems, reported recently, i.e., [Fe(tpa®)]⁻ (R = tert-butyl, Tbu (1), mesityl, Mes (2), phenyl, Ph (3), and 2,6-difluorophenyl, Dfp (4), are enough to lead to five lowest and thermally accessible spin sublevels described sufficiently well by S = 2 SH provided that it is extended with one fourth order anisotropy term. Using this most elementary parametrization that is consistent with the actual physics, the reported magnetization data for the target systems were reinterpreted and found to be in good agreement with the ab initio results. The multiplet energies from the ab initio calculations have been fitted with remarkable consistency using a ligand field (angular overlap) model (ab initio ligand field, AILFT). This allows for determination of bonding parameters and quantitatively demonstrates the correlation between increasingly negative D values and changes in the σ-bond strength induced by the peripheral ligands. In fact, the sigma-bonding capacity (and hence the Lewis basicity) of the ligand decreases along the series 1 > 2 > 3 > 4

Zero field splitting of the chalcogen diatomics using relativistic correlated wave-function methods.

International audienceThe spectrum arising from the (π*)(2) configuration of the chalcogen dimers, namely, the X(2)1, a2, and b0(+) states, is calculated using wave-function theory based methods. Two-component (2c) and four-component (4c) multireference configuration interaction (MRCI) and Fock-space coupled cluster (FSCC) methods are used as well as two-step methods spin-orbit complete active space perturbation theory at 2nd order (SO-CASPT2) and spin-orbit difference dedicated configuration interaction (SO-DDCI). The energy of the X(2)1 state corresponds to the zero-field splitting of the ground state spin triplet. It is described with high accuracy by the 2- and 4-component methods in comparison with experiment, whereas the two-step methods give about 80% of the experimental values. The b0(+) state is well described by 4c-MRCI, SO-CASPT2, and SO-DDCI, but FSCC fails to describe this state and an intermediate Hamiltonian FSCC ansatz is required. The results are readily rationalized by a two-parameter model; Δε, the π* spinor splitting by spin-orbit coupling and K, the exchange integral between the π(1)* and the π(-1)* spinors with, respectively, angular momenta 1 and -1. This model holds for all systems under study with the exception of Po(2)

Improving I/O Performance for Exascale Applications through Online Data Layout Reorganization

The applications being developed within the U.S. Exascale Computing Project (ECP) to run on imminent Exascale computers will generate scientific results with unprecedented fidelity and record turn-around time. Many of these codes are based on particle-mesh methods and use advanced algorithms, especially dynamic load-balancing and mesh-refinement, to achieve high performance on Exascale machines. Yet, as such algorithms improve parallel application efficiency, they raise new challenges for I/O logic due to their irregular and dynamic data distributions. Thus, while the enormous data rates of Exascale simulations already challenge existing file system write strategies, the need for efficient read and processing of generated data introduces additional constraints on the data layout strategies that can be used when writing data to secondary storage. We review these I/O challenges and introduce two online data layout reorganization approaches for achieving good tradeoffs between read and write performance. We demonstrate the benefits of using these two approaches for the ECP particle-in-cell simulation WarpX, which serves as a motif for a large class of important Exascale applications. We show that by understanding application I/O patterns and carefully designing data layouts we can increase read performance by more than 80 percent

Study of the ρ\rho, ω\omega, ϕ→ηγ→7γ\phi\to\eta\gamma\to 7\gamma Decays with an SND Detector on a VEPP-2M Collider

The $e^+e^-\to\eta\gamma\to 7\gamma$ process was studied in the energy range $2E=600\div 1060$ MeV with an SND detector on a VEPP-2M $e^+e^-$ collider. The decay branching ratios $B(\phi\to\eta\gamma)=(1.343\pm 0.012\pm 0.055)\cdot 10^{-2}$, $B(\omega\to\eta\gamma)=(4.60\pm 0.72\pm 0.19)\cdot 10^{-4}$, and $B(\rho\to\eta\gamma)=(2.69\pm 0.32\pm 0.16)\cdot 10^{-4}$ were measured.Comment: 5 pages, 4 figure

Search for the radiative decay η→π0γγ\eta \to \pi^0 \gamma \gamma in the SND experiment at VEPP-2M

The $\eta \to \pi^0 \gamma \gamma$ decay was investigated by the SND detector at VEPP-2M $e^+e^-$ collider in the reaction $e^+e^-\to\phi\to \eta\gamma$. Here we present the results and some details of this study. We report an upper limit (90% c.l.) $Br(\eta \to \pi^0 \gamma \gamma)<8.4\times 10^{-4}$ as our final result. Our upper limit does not contradict the earlier measurement by GAMS spectrometer. To facilitate future studies a rather detailed review of the problem is also given.Comment: 24 pages, 6 figures, LaTex. To be published in Nucl. Phys.

Density functional theory

Density functional theory (DFT) finds increasing use in applications related to biological systems. Advancements in methodology and implementations have reached a point where predicted properties of reasonable to high quality can be obtained. Thus, DFT studies can complement experimental investigations, or even venture with some confidence into experimentally unexplored territory. In the present contribution, we provide an overview of the properties that can be calculated with DFT, such as geometries, energies, reaction mechanisms, and spectroscopic properties. A wide range of spectroscopic parameters is nowadays accessible with DFT, including quantities related to infrared and optical spectra, X-ray absorption and Mössbauer, as well as all of the magnetic properties connected with electron paramagnetic resonance spectroscopy except relaxation times. We highlight each of these fields of application with selected examples from the recent literature and comment on the capabilities and limitations of current methods

Antisymmetric Magnetic Interactions in Oxo-Bridged Copper(II) Bimetallic Systems

The antisymmetric magnetic interaction is studied using correlated wave-function-based calculations in oxo-bridged copper bimetallic complexes. All of the anisotropic multispin Hamiltonian parameters are extracted using spin-orbit state interaction and effective Hamiltonian theory. It is shown that the methodology is accurate enough to calculate the antisymmetric terms, while the small symmetric anisotropic interactions require more sophisticated calculations. The origin of the antisymmetric anisotropy is analyzed, and the effect of geometrical deformations is addressed.