Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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.

At least 289 records · Page 16

Reactivity of CO 2 with Utica, Marcellus, Barnett, and Eagle Ford Shales and Impact on Permeability

We report that in order to reduce greenhouse gas emissions while recovering hydrocarbons from unconventional shale formations, processes that make use of carbon dioxide to enhance oil recovery while storing carbon dioxide (CO 2 ) should be considered. Here, we examine samples from three shale basins across the United States (Utica and Marcellus Shales in the Appalachian Basin, Barnett Shale in the Bend Arch-Ft. Worth Basin, and Eagle Ford in the Western Gulf Basin) to address the following questions: (1) do changes from reaction with CO 2 and fluids at the micrometer and nanometer scale alter flow pathways and, in turn, impact hydrocarbon production, CO 2 storage, and seal integrity and (2) can CO 2 or fluid reactivity be predicted based on physical or chemical properties of shale formations? Experiments were conducted at 40 °C and 10.3 MPa to characterize the interaction between CO 2 and shale using X-ray diffraction (XRD), carbon and sulfur analysis, in situ Fourier transform infrared spectroscopy (FT-IR), feature relocation scanning electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), mercury (Hg) intrusion porosimetry, and Brunauer-Emmett-Teller (BET) surface area and pore size analysis coupled with density functional theory (DFT) methods. Changes in mechanical, physical, and flow properties of shale cores due to CO 2 exposure were addressed using a New England Research Autolab 1500 and Xenon X-ray computed tomography (CT) scanning. Results showed that CO 2 did not promote significant reactivity with the shale if water was not present; only shales with swelling clays or residual interstitial pore water reacted with dry CO 2 to promote reactivity in shale. When water was added as a reactant, CO 2 formed carbonic acid and reacted with the shale to dissolve carbonate pockets, etched and pitted the shale matrix surfaces, and increased the microporosity and decreased nanoporosity. Porosity and permeability increased appreciably in core shale samples after exposure to CO 2 saturated fluid due to dissolution of carbonate. Shale mechanical properties were not altered. Trends were not observed that could tie CO 2 or fluid reactivity to physical or chemical properties of the shale formations at the basin scale from the samples we examined. However, if the shale contained significant amounts of carbonate and water was available to react with the CO 2 , pore sizes were altered in the matrix and permeability and porosity increased.

04 OIL SHALES AND TAR SANDS↗

Radiation Effects on Uranyl Tetrachloro Coordination Compounds: Impact of Lattice Water

Nuclear materials, such as uranium-bearing solids, are exposed to high levels of ionizing radiation throughout the nuclear fuel cycle; thus, it is important to develop a molecular-level understanding of how these materials behave and degrade in the presence of gamma (γ) irradiation. In the current study, three U(VI) tetrachloride complexes, M 2 [UO 2 Cl 4 ]·xH 2 O (where M = K + , Rb + , or Cs + and x = 0 or 2), and their respective chloride salts were exposed to 1–50 kGy of γ radiation using a 60 Co source. Irradiated materials were evaluated by using electron paramagnetic resonance (EPR) and Raman spectroscopy and were further explored by using density functional theory (DFT) methods. EPR spectra of the irradiated materials suggest the formation of a Cl-based radical for both the alkali salts and the uranyl tetrachloride compounds, and DFT calculations provide evidence that the Cl 2 –• radical is formed within these materials. The presence of water in the K + and Rb + compounds leads to additional spectroscopic signatures that could be traced back to water radiolysis and the formation of peroxide and superoxide species. DFT results support the formation of HO 2 • in the lattice and potentially the formation of a [UO 2 Cl 3 (O 2 )] 3– species, highlighting the impact of water within the hydrated material to alter U(VI) speciation by radiolysis.

36 MATERIALS SCIENCE↗

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmark Data Set of Crystalline Organic Semiconductors

This work reports a Benchmark Data set of Crystalline Organic Semiconductors to test calculations of the structural and electronic properties of these materials in the solid state. The data set contains 67 crystals consisting of mostly rigid molecules with a single dominant conformer, covering the majority of known structural types. The experimental crystal structure is available for the entire data set, whereas zero-temperature unit cell volume can be reliably estimated for a subset of 28 crystals. Using this subset, we benchmark r 2 SCAN-D3 and PBE-D3 density functionals. Then, for the entire data set, we benchmark approximate density functional theory (DFT) methods, including GFN1-xTB and DFTB3(3ob-3-1), with various dispersion corrections against r 2 SCAN-D3. Our results show that r 2 SCAN-D3 geometries are accurate within a few percent, which is comparable to the statistical uncertainty of experimental data at a fixed temperature, but the unit cell volume is systematically underestimated by 2% on average. The several times faster PBE-D3 provides an unbiased estimate of the volume for all systems except for molecules with highly polar bonds, for which the volume is substantially overestimated in correlation with the underestimation of atomic charges. Considered approximate DFT methods are orders of magnitude faster and provide qualitatively correct but overcompressed crystal structures unless the dispersion corrections are fitted by unit cell volume.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium–Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D e and bond length r e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15–0.19. For the single-reference methods, Kohn–Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D e and r e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core–valence terms; this results in highly accurate dissociation energies and good bond lengths. In conclusion, averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl 2 by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl 2 and 50 mol % MgCl 2 + 50 mol % KCl, and ZnCl 2 , thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl 2 upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl 2 , which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl 5 3– complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl 2 content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unveiling the Redox Noninnocence of Metallocorroles: Exploring K-Edge X-ray Absorption Near-Edge Spectroscopy with a Multiconfigurational Wave Function Approach

X-ray absorption near-edge spectroscopy (XANES) is an advanced technique for probing the local electronic structure of catalysts, effectively identifying the noninnocent nature of ligands in transition-metal complexes. Metallocorroles with noninnocent corrole rings exhibit unusual electronic structures that challenge traditional density functional theory (DFT) methods, necessitating more rigorous approaches to describe electron correlation accurately. We explored K-edge XANES spectra of Fe, Mn, and Co metallocorroles using TDDFT and wave function-based methods. This is the first investigation employing multireference methods, specifically RASSCF, RASPT2, and MC-PDFT, to analyze the redox noninnocent nature of metallocorroles reflected in their XANES spectra. We quantified the noninnocent character of the corrole and the oxidation states of the metals, capturing more than singly excited excitations responsible for the pre-edge peak. Our findings demonstrate the importance of these advanced computational techniques for accurately predicting XANES spectra, providing a reliable understanding of the electronic properties of such complexes. In conclusion, this study offers a new strategy for investigating ligand redox noninnocence via integrated experimental and computational XANES.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling and Optimizing Photoinduced Charge Transfer across Ultrathin Silica Separation Membrane with Embedded Molecular Wires for Artificial Photosynthesis

Ultrathin amorphous silica membranes with embedded organic molecular wires (oligo( p -phenylenevinylene), three aryl units) provide chemical separation of incompatible catalytic environments of CO 2 reduction and H 2 O oxidation while maintaining electronic and protonic coupling between them. For an efficient nanoscale artificial photosystem, important performance criteria are high rate and directionality of charge flow. Here, the visible-light-induced charge flow from an anchored Ru bipyridyl light absorber across the silica nanomembrane to Co 3 O 4 water oxidation catalyst is quantitatively evaluated by photocurrent measurements. Charge transfer rates increase linearly with wire density, with 5 nm -2 identified as an optimal target. Accurate measurement of wire and light absorber densities is accomplished by the polarized FT-IRRAS method. Guided by density functional theory (DFT) calculations, four wire derivatives featuring electron-donating (methoxy) and -withdrawing groups (sulfonate, perfluorophenyl) with highest occupied molecular orbital (HOMO) potentials ranging from 1.48 to 0.64 V vs NHE were synthesized and photocurrents evaluated. Charge transfer rates increase sharply with increasing driving force for hole transfer from the excited light absorber to the embedded wire, followed by a decrease as the HOMO potential of the wire moves beyond the Co 3 O 4 valence band level toward more negative values, pointing to an optimal wire HOMO potential around 1.3 V vs NHE. Comparison with photocurrents of samples without nanomembrane indicates that silica layers with optimized wires are able to approach undiminished electron flux at typical solar intensities. Combined with the established high proton conductivity and small-molecule blocking property, the charge transfer measurements demonstrate that oxidation and reduction catalysis can be efficiently integrated on the nanoscale under separation by an ultrathin silica membrane.

photocurrent measurements↗

Unusual Role of Point Defects in Perovskite Nickelate Electrocatalysts

Low-cost transition-metal oxide is regarded as a promising electrocatalyst family for an oxygen evolution reaction (OER). The classic design principle for an oxide electrocatalyst believes that point defect engineering, such as oxygen vacancies (VO .. ) or heteroatom doping, offers the opportunities to manipulate the electronic structure of material toward optimal OER activity. Oppositely, in this work, we discover a counterintuitive phenomenon that both VO .. and an aliovalent dopant (i.e., proton (H + )) in perovskite nickelate (i.e., NdNiO 3 (NNO)) have a considerably detrimental effect on intrinsic OER performance. Detailed characterizations unveil that the introduction of these point defects leads to a decrease in the oxidative state of Ni and weakens Ni-O orbital hybridization, which triggers the local electron-electron correlation and a more insulating state. Evidenced by first-principles calculation using the density functional theory (DFT) method, the OER on nickelate electrocatalysts follows the lattice oxygen mechanism (LOM). The incorporation of point defect increases the energy barrier of transformation from OO*(V O ) to OH*(V O ) intermediates, which is regarded as the rate-determining step (RDS). Furthermore, this work offers a new and significant perspective of the role that lattice defects play in the OER process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemistry and Interfacial Structure Promoting Quasi-van der Waals Epitaxial Growth of WS 2 Nanosheets on Sapphire for Prospective Application in Field-Effect Transistors

How do chemical and structural modifications to the supporting crystal surface affect the subsequent van der Waals (vdW) or quasi(Q)-vdW epitaxial growth of 2D nanocrystals? Developing an atomic-scale picture of such an interfacial system is crucial for understanding its impact on the physical and chemical properties of the supported 2D materials. The elucidation of the interfacial structure and chemistry needed to promote the Q-vdW epitaxial growth of 2D tungsten disulfide (WS 2 ) nanocrystals contributes to the growth mechanism understanding, thus pushing forward the integration of such atomically thin semiconductors toward real field-effect transistor applications. In addition to an atomic-force microscopy top view, we showcase a combination of X-ray techniques for a top-to-bottom investigation of the complexities of the buried interface structures. Furthermore, this approach uses X-ray photoelectron spectroscopy, X-ray standing wave excited X-ray fluorescence, and crystal truncation rod scattering to produce a highly resolved chemical-state-specific 3D atomic map for the extended interface structure of WS 2 /α-Al 2 O 3 (001). Employing these detailed analysis methods, along with density functional theory to further refine the picoscale structure, we demonstrate how two different types of interface engineering during the pregrowth stage lead to significant differences in the chemical and structural modifications to the terminal surface of c-face sapphire, which in turn leads to substantial differences in the submonolayer growth of supported WS 2 2D nanocrystals in terms of lateral domain sizes, epitaxial registry, vdW gaps, and stability.

36 MATERIALS SCIENCE↗

Revealing Fast Cu-Ion Transport and Enhanced Conductivity at the CuInP 2 S 6 –In 4/3 P 2 S 6 Heterointerface

Van der Waals layered ferroelectrics, such as CuInP 2 S 6 (CIPS), offer a versatile platform for miniaturization of ferroelectric device technologies. Control of the targeted composition and kinetics of CIPS synthesis enables the formation of stable self-assembled heterostructures of ferroelectric CIPS and nonferroelectric In 4/3 P 2 S 6 (IPS). Here, we use quantitative scanning probe microscopy methods combined with density functional theory (DFT) to explore in detail the nanoscale variability in dynamic functional properties of the CIPS-IPS heterostructure. Within, we report evidence of fast ionic transport which mediates an appreciable out-of-plane electromechanical response of the CIPS surface in the paraelectric phase. Further, we map the nanoscale dielectric and ionic conductivity properties as we thermally stimulate the ferroelectric-paraelectric phase transition, recovering the local dielectric behavior during this phase transition. Finally, aided by DFT, we reveal a substantial and tunable conductivity enhancement at the CIPS/IPS interface, indicating the possibility of engineering its interfacial properties for next generation device applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

74 ATOMIC AND MOLECULAR PHYSICS↗

Kinetics and mechanisms of polycondensation reactions between aryl halides and bisphenol A

Aryl chlorides (ArCl) or aryl fluorides (ArF) were used in polycondensation reactions to form poly(arylene ether sulfone)s (PAES). Interestingly, the kinetics of the ArF reaction fit a third-order rate law, which is attributed to the activation of the carbon–fluorine bond by two potassium cations (at least one bound to phenolate), which form a three-body complex. Here, the ArCl monomer follows a second-order rate law, where a two-body complex forms at the initial state of the aromatic nucleophilic substitution (S N Ar) pathway. These metal cation-activated complexes act as intermediates during the attack by the nucleophile. This finding was reproduced with either the potassium or the sodium counterion (introduced via potassium carbonate or sodium carbonate). Through a combination of experimental analysis of reaction kinetics and computational calculations with density functional theory (DFT) methods, the present work extends the fundamental understanding of polycondensation mechanisms for two aryl halides and highlights the importance of the CX–metal interaction(s) in the S N Ar reaction, which is translational to other ion-activated substitution reactions.

36 MATERIALS SCIENCE↗

Valence-to-core X-ray emission spectroscopy of vanadium oxide and lithiated vanadyl phosphate materials

In this work, we report valence-to-core (VTC) X-ray emission spectroscopy (XES) measurements and theoretical calculations of the electrochemical sequence ε-VOPO 4 , ε-LiVOPO 4 , ε-Li 2 VOPO 4 and the reference oxides V 2 O 3 , VO 2 , and V 2 O 5 . In our analysis of these results, we establish a framework for interrogating chemical bonding that is generally applicable to a wide range of systems, including complex, extended inorganic compounds. While this latter regime has garnered less focused application than, e.g. , metalloenzymes in many excellent catalysis studies, we show that the technique provides high utility in materials-focused energy storage research. Here, sensitivities to the local atomic structure and hybridization schemes are discussed in detail. Similarly, the effect of lithiation on oxidation, delocalization, and shifts in ligand valence energy levels are all readily apparent in the analyzed results. Finally, the TDDFT projections clearly reveal the directional dependencies of the valence band at each of the vanadium sites. Our results demonstrate laboratory-based X-ray spectroscopy instrumentation is a viable route for attaining well-resolved VTC-XES features for inorganic compounds of 3d transition metals, even for samples of limited quantity or suffering from sensitivity to the atmosphere. The experimental results are in good agreement with results produced by real-space Green's function and time-dependent density functional theory (TDDFT) methods, respectively. Hence, we propose that VTC-XES, when equipped with appropriate theoretical support, can be a valuable complement to X-ray absorption pre-edge features for more detailed characterization of a compound's electronic structure. We expect similar analyses will find application in a broad range of materials chemistry research and provide both fundamental and applied insights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A substituent-dependent deep-blue pyrene-based chemosensor for trace nitroaniline sensing

Nitroaniline (NA) is a kind of high toxicity and pathogenic amine derivative, and the trace monitoring of nitroaniline has become an urgent problem in terms of both the environment and security. Herein, two K-region multiple substituted pyrene-based luminogens with reasonable photoluminescence efficiency and blue emission (408–430 nm) have been synthesized and characterized. Their photophysical properties were carefully investigated in solution and in the solid state by experimental and density functional theory (DFT) methods. It is worth mentioning that both of these π-electron rich and extended π-conjugated materials were used for the detection of electron deficient nitroaromatics (NACs) based on a photo-induced electron transfer (PET) mechanism, and a sensitive response towards NA, especially for p-NA with a limit of detection (LOD) as low as 8.58 × 10 -9 M, was observed, with a fluorescence quenching efficiency up to 99.5%. In conclusion, these research results provided a low-cost and simple strategy for trace nitroaniline sensing.

36 MATERIALS SCIENCE↗

Advanced thermal/environmental barrier coatings of high-entropy rare earth disilicates tuned by strong anharmonicity of Eu 2 Si 2 O 7

Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.

36 MATERIALS SCIENCE↗

Hyperfine and quadrupole interactions for Dy isotopes in DyPc 2 molecules

Nuclear spin levels play a vital role in understanding magnetization dynamics and implementation and control of quantum bits in lanthanide-based single-molecule magnets. We investigate the hyperfine and nuclear quadrupole interactions for 161 Dy and 163 Dy nuclei in anionic DyPc 2 (Pc=phthalocyanine) single-molecule magnets, using multiconfigurational ab initio methods (beyond density-functional theory) including spin–orbit interaction. The two isotopes of Dy are chosen because the others have zero nuclear spin. Both isotopes have the nuclear spin I =5/2, although the magnitude and sign of the nuclear magnetic moment differ from each other. The large energy gap between the electronic ground and first-excited Kramers doublets, allows us to map the microscopic hyperfine and quadrupole interaction Hamiltonian onto an effective Hamiltonian with an electronic pseudo-spin that corresponds to the ground Kramers doublet. Our ab initio calculations illustrate that the coupling between the nuclear spin and electronic orbital angular momentum contributes the most to the hyperfine interaction and that both the hyperfine and nuclear quadrupole interactions for 161 Dy and 163 Dy nuclei are much smaller than those for the 159 Tb nucleus in TbPc 2 single-molecule magnets. The calculated separations of the electronic-nuclear levels are comparable to experimental data reported for 163 DyPc 2 . We demonstrate that hyperfine interaction for the Dy Kramers ion leads to tunnel splitting (or quantum tunneling of magnetization) at zero field. This effect does not occur for TbPc 2 single-molecule magnets. The magnetic field values of the avoided level crossings for 161 DyPc 2 and 163 DyPc 2 are found to be noticeably different, which can be observed from the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗