Interplay of 3d and 4f Magnetism in Chiral Y 6 FeSi 2 S 14 and Tb 6 FeSi 2 S 14 Chalcogenides
Explore the source record for details and available documents.
Gamage, Eranga H. [Department of Chemistry; Iowa S↗
SEARCH · Engineering Papers
Results for “Occupant”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Structural Phase Transitions in the van der Waals Ferromagnets Fe x Pd y Te 2
Two-dimensional van der Waals (vdW) magnets are attracting significant attention, both as platforms for studying fundamental magnetic interactions and for the exciting possibility of utilizing them as building blocks in devices and heterostructures, which may lead to new physical phenomena and functionalities. Here, we provide a detailed study of the crystal structure and physical properties of the recently discovered vdW ferromagnet FePd 2 Te 2 . We find this compound has a relatively wide width of formation, and grow single crystals with compositions Fe x Pd y Te 2 where x ranges from 0.9 to 1.1 and y from 1.8 to 2.5, respectively. Temperature-dependent X-ray diffraction and transport measurements reveal that a first-order structural transition occurs in the range of T = 360–420 K, where the critical temperature, modulation wave vector, and corresponding room-temperature crystal structures all depend on chemical composition. Above the transition, the compounds with Pd fraction y > 2 adopt a disordered derivative of the tetragonal FeTe structure, with the Fe layer showing mixed Fe/Pd occupancy and the extra Pd atoms partially occupying interstitial sites. Below 370 K, the structure is incommensurately modulated, likely associated with the complex ordering of Pd/Fe atoms in the metal layers or the interstitial Pd in the vdW gaps. For y < 2, the composition Fe 1.1 Pd 1.8 Te 2 has monoclinic symmetry at room temperature that is consistent with the reported structure of FePd 2 Te 2 . This phase undergoes a structural transition at 420 K for which the high temperature structure is yet to be determined; however, based on the similarities with the y > 2 compounds, we speculate that its T > 420 K structure is also tetragonal. Importantly, the high temperature, symmetry-breaking structural transition observed here provides a likely explanation for the origin of the structural domains previously observed in FePd 2 Te 2 . All compounds investigated in the Fe x Pd y Te 2 series show metallic behavior, with magnetic characterization indicating that they are easy-plane, hard, ferromagnets with T C spanning 98–180 K. Both the critical temperature for the structural transition and the Curie temperature are moderately suppressed with increasing Pd fraction y and corresponding decreasing Fe fraction x, indicating that synthetic control over x and y paves way for the further exploration of these compounds.
crystal structure↗
Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal–Organic Framework
In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.
adsorption capacities↗
Understanding the Stabilization and Tunability of Divalent Europium 2.2.2B Cryptates
Lanthanides such as europium with more accessible divalent states are useful for studying redox stability afforded by macrocyclic organic ligands. Substituted cryptands, such as 2.2.2B cryptand, that increase the oxidative stability of divalent europium also provide coordination environments that support synthetic alterations of Eu(II) cryptate complexes. Two single crystal structures were obtained containing nine-coordinate Eu(II) 2.2.2B cryptate complexes that differ by a single coordination site, the occupation of which is dictated by changes in reaction conditions. A crystal structure containing a [Eu(2.2.2B)Cl] + complex is obtained from a methanol–THF solvent mixture, while a methanol–acetonitrile solvent mixture affords a [Eu(2.2.2B)(CH 3 OH)] 2+ complex. While both crystals exhibit the typical blue emission observed in most Eu(II) containing compounds as a result of 4f 6 5d 1 to 4f 7 transitions, computational results show that the substitution of a Cl– anion in the place of a methanol molecule causes mixing of the 5d excited states in the Eu(II) 2.2.2B cryptate complex. Additionally, magnetism studies reveal the identity of the capping ligand in the Eu(II) 2.2.2B cryptate complex may also lead to exchange between Eu(II) metal centers facilitated by π-stacking interactions within the structure, slightly altering the anticipated magnetic moment. Finally, the synthetic control present in these systems makes them interesting candidates for studying less stable divalent lanthanides and the effects of precise modifications of the electronic structures of low valent lanthanide elements.
Anions↗
Main and Satellite Features in the Ni 2p XPS of NiO
The origin and assignment of the complex main and satellite XPS features of the cations in ionic compounds has been the subject of extensive theoretical studies using different methods. There is agreement that within a molecular orbital model one needs to take into account different types of configurations. Specifically, those where a core electron is removed but no other configuration changes are made and those where in addition to ionization there are also shake or charge transfer changes to the ionic configuration. However, there are strong disagreements about the assignment of XPS features to these configurations. The present work is directed toward resolving the origin of main and satellite features for the Ni 2p XPS of NiO based on ab initio molecular orbital wavefunctions for a cluster model of NiO. A major problem in earlier ab initio XPS studies of ionic compounds has been the use of a common set of orbitals that was not able to properly describe all the ionic configurations that contribute to the full XPS spectra. This is resolved in the present work by using orbitals that are optimized for averages of the occupations of the different configurations that contribute to the XPS. The approach of using State Averaged orbitals is validated through comparisons between different averages and through use of higher order excitations in the wavefunctions for the ionic states. It represents a major extension of our earlier work on the main and satellite features of the Fe 2p XPS of Fe 2 O 3 and proves the reliability and the generality of the assignments of the character and origin of the different features of the XPS obtained with orbitals optimized for State Averages. These molecular orbital methods permit the characterization of the ionic states in terms of the importance of shake excitations and of the coupling of ionization of 2p 1/2 and 2p 3/2 spin-orbit split sub shells. Further, the work lays the foundation for definitive assignments of the character of main and satellite XPS features and points to their origin in the electronic structure of the material.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Cause, Consequence, and Control of Ag Vacancies in BaAg 2– x As 2
Explore the source record for details and available documents.
Kyveryga, Victoria↗
Bis( tert -butoxydiphenylsilyl)amide Divalent Lanthanide Complexes
The development of new ligand systems to stabilize “nontraditional/ non-classical” divalent lanthanides is key to tuning the chemical and physical properties of their mixed principal quantum number 4f n 5d 1 ground states. The design and study of novel ligand systems which stabilize occupation of differing orbitals within the 5d manifold for these ions constitutes an area ripe for exploration. Our efforts toward the development of redox-innocent bulky silylamide ligands to stabilize pseudo-octahedral coordination geometries for divalent lanthanides have resulted in the synthesis of the bis( tert -butoxydiphenylsilyl) amide ligand, whose coordination complexes with Sm 2+ , Eu 2+ , and Yb 2+ are reported herein. These systems have been fully characterized by single-crystal X-ray diffraction, elemental analysis, cyclic voltammetry, direct-current magnetometry, and infrared, nuclear magnetic resonance, and electronic absorption spectroscopies. Attempts to extend this system to the more reducing Tm 2+ ion resulted in an inseparable mixture of products from which crystals of the analogous Tm 2+ species and a reduced dinitrogen, bimetallic Tm 3+ -Tm 3+ complex bridged by a η 2 -N 2 3− radical could be identified. Though progress toward six-coordinate complexes of reducing “traditional/ classical” divalent ions is noted for these systems, further work is needed to improve the synthetic utility of this ligand framework for the study of “non-traditional/non-classical” divalent lanthanides with a mixed-principal quantum number 4f n 5d 1 ground state.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structural Investigation of Six Quinary Sulfides Synthesized via the Flux-Assisted Boron-Chalcogen Mixture (BCM) Method: Eu 2+ Containing Members of the RE 3 MTQ 7 (M and T = Transition or Main Group Metals, Q = Chalcogens) Family
For this work, a series of six quinary rare-earth sulfides Ce 4+ 1.85 Eu 2+ 1.15 Na 0.30 SiS 7 , Ce 4+ 1.91 Eu 2+ 1.09 K 0.18 SiS 7 , Ce 4+ 1.96 Eu 2+ 1.04 Rb 0.08 SiS 7 , Ce 4+ 1.98 Eu 2+ 1.02 Cs 0.05 SiS 7 , Ce 4+ 1.97 Eu 2+ 1.03 Ag 0.06 SiS 7 , and Ce 4+ 1.50 Eu 2+ 1.50 CuSiS 7 were obtained in an alkali iodide flux using the boron-chalcogen mixture (BCM) method. Single crystal X-ray diffraction was used to determine the structures of the high quality single crystals that were grown; their elemental compositions were confirmed by energy-dispersive spectroscopy (EDS). The compounds crystallize in the hexagonal crystal system in the noncentrosymmetric space group P63. The crystal structure consists of a three-dimensional network composed of mixed cerium and europium bicapped trigonal prisms, isolated SiS4 tetrahedra, and monovalent metals (Na, K, Rb, Cs, Ag, and Cu) located in cavities created by linked Ce/EuS 8 polyhedra. The structures are charge-balanced when Ce and Eu are in their +4 and +2 oxidation states, respectively. The effective magnetic moment of Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 determined from the temperature dependence of the magnetic susceptibility data is consistent with the presence of Ce 4+ and Eu 2+ . Clear correlations between the alkali ion site occupancy, the ionic radius of the alkali cations, and the average bond length of Ce 4+ /Eu 2+ –S, were established. UV–vis diffuse reflectance data were collected for Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 and a band gap of 1.9(1) eV was established.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
How Does Metal Spin State Affect Electronic Communication in Mixed-Valence Dimers? Insights from Ultrafast Near-Infrared and Soft X-ray Transient Absorption Spectroscopy
Recent advancements in photocatalysis, photovoltaics, and quantum information science take advantage of electron spin, and determining how spin multiplicity affects electron transfer is key to understanding these phenomena. Here, in this study, we examine how metal spin state affects electronic communication in an organometallic mixed-valence dimer, ferrocenyl cobaltocenium ([Fe II Cp 2 Co III Cp 2 ] + ). This complex can be photoexcited from its low-spin singlet Fe II ground state to form intermediate-spin triplet Fe II and high-spin quintet Fe II excited states. Using femtosecond optical transient absorption (OTA) spectroscopy with visible (vis), near-infrared (NIR), and short-wave IR (SWIR) probes, supported by time-dependent density functional theory (TD-DFT) calculations, we measure Fe II Co III →Fe III Co II intervalence charge transfer (IVCT) bands in each of the Fe II spin states. Mulliken–Hush analysis of the excited-state IVCT bands was used to compute the electronic coupling between the metal centers in all three spin states, which increased as quintet < triplet < singlet. Meanwhile, the peak energy of the bands, and thus the ΔG of the IVCT transition, increased as triplet < quintet < singlet. Then, to directly probe the electronic structure at each metal center, we employed picosecond soft X-ray transient absorption (XTA) spectroscopy at the Fe and Co L 3 edges. Our results show that the low-spin and high-spin states of [Fe II Cp 2 Co III Cp 2 ] + are valence-localized, while the intermediate-spin state is partially delocalized. The differences in charge delocalization are attributed to differences in orbital occupation and geometry that affect the free energy and superexchange coupling.
Burke, John H. [Univ. of Illinois at Urbana-Champa↗
Synthesis, Crystal Structure, and Cooperative 3d–5d Magnetism in Rock Salt Type Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6
Two new transition metal oxides with the nominal chemical compositions of Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6 were successfully synthesized. Both compounds crystallize in an ordered rock salt structure type in the monoclinic C2/m space group. The crystal structures were determined using both synchrotron X-ray and time-of-flight neutron, powder diffraction data. In both phases, Ni2+ ions are present while oxidation states of osmium are +6 and +5 in Li 4 NiOsO 6 and Li 3 Ni 2 OsO 6 , respectively. Ni2+ ions in the hypothetical fully ordered phase form a honeycomb arrangement in the ab crystallographic plane and these hexagons are centered by osmium ions. Furthermore, the magnetic layers are separated along the c axis by the octahedra, which are centered by Li+ (or Li+/Ni 2+ , depending on the chemical compositions). Crystal structure refinements reveal that there is some degree of mixed occupancy in cationic positions. Temperature dependent magnetic susceptibility data for both phases show ferrimagnetic transitions with predominant antiferromagnetic (AFM) interactions among 3d electrons of nickel and Sd electrons of osmium. Iso-thermal magnetization loops as a function of the applied magnetic field below the transition temperatures confirm the ferrimagnetic nature in magnetic transitions. Temperature dependent heat capacity data, however, did not exhibit any anomaly in either phase, indicating the absence of long-range magnetic ordering. The lack of long-range order for both Os 5+ and Os 6+ -based compounds was also confirmed by low temperature neutron diffraction data down to 10 K. Temperature dependent AC magnetic susceptibility data in various frequencies for both samples indicate that Li 4 NiOsO 6 exhibits spin-glass-like behavior, while the transition temperature for Li 3 Ni 2 OsO 6 is nearly frequency independent.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results
Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics
Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo
The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective
Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
I ntera C hem : Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics
InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree–Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry
Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters
The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)
In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.
Cluster chemistry↗