Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relativistic effects”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Relativistic Density Functional NMR Tensors Analyzed with Spin–free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the 15 N shielding in a heavy metal complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Relativistic coupled‐cluster and equation‐of‐motion coupled‐cluster methods

Abstract The development of relativistic coupled‐cluster (CC) and equation‐of‐motion coupled‐cluster (EOM‐CC) methods is reviewed. An emphasis is placed on recent efforts to improve the computational efficiency of CC and EOM‐CC calculations with non‐perturbative treatments of spin‐orbit coupling (SO‐CC and EOM‐CC) by partially recovering spin symmetry in the formulations. Example calculations of electronic ground state as well as valence‐excited and core‐excited states for molecules containing heavy elements are presented to demonstrate the applicability and usefulness of the SO‐CC and EOM‐CC methods. Future directions for the development of the SO‐CC and EOM‐CC methods are also discussed. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Gravitational Waves from Nnaturalness

Abstract We study the prospects for probing the Nnaturalness solution to the electroweak hierarchy problem with future gravitational wave observatories. Nnaturalness, in its simplest incarnation, predictsNcopies of the Standard Model with varying Higgs mass parameters. We show that in certain parameter regions the scalar reheaton transfers a substantial energy density to the sector with the smallest positive Higgs squared mass while remaining consistent with bounds on additional effective relativistic species. In this sector, all six quarks are much lighter than the corresponding QCD confinement scale, allowing for the possibility of a first-order chiral symmetry-breaking phase transition and an associated stochastic gravitational wave signal. We consider several scenarios characterizing the strongly-coupled phase transition dynamics and estimate the gravitational wave spectrum for each. Pulsar timing arrays (SKA), spaced-based interferometers (BBO, Ultimate-DECIGO,μAres, asteroid ranging), and astrometric measurements (THEIA) all have the potential to explore new regions of Nnaturalness parameter space, complementing probes from next generation cosmic microwave background radiation experiments.

Physics↗
DOE OSTI2020

A testable hidden-sector model for Dark Matter and neutrino masses

We consider a minimal extension of the Standard Model with a hidden sector charged under a dark local U(1)' gauge group, accounting simultaneously for light neutrino masses and the observed Dark Matter relic abundance. The model contains two copies of right-handed neutrinos which give rise to light neutrino-masses via an extended seesaw mechanism. The presence of a stable Dark-Matter candidate and a massless state naturally arise by requiring the simplest anomaly-free particle content without introducing any extra symmetries. We investigate the phenomenology of the hidden sector considering the U(1)' breaking scale of the order of the electroweak scale. Confronting the thermal history of this hidden-sector model with existing and future constraints from collider, direct and indirect detection experiments provides various possibilities of probing the model in complementary ways as every particle of the dark sector plays a specific cosmological role. Across the identified viable parameter space, a large region predicts a sizable contribution to the effective relativistic degrees-of-freedom in the early Universe that allows to alleviate the recently reported tension between late and early measurements of the Hubble constant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The radial action from probe amplitudes to all orders

We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic effect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe off a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Coupled cluster spectroscopic properties of the coinage metal nitrosyls, M–NO (M = Cu, Ag, Au)

Ab initio near-equilibrium potential energy and dipole moment surfaces for the bent CuNO, AgNO, and AuNO molecules have been calculated under the Feller-Peterson-Dixon (FPD) composite framework at the coupled cluster level of theory including complete basis set extrapolation, outer-core correlation, scalar relativistic effects, and spin-orbit coupling. The Brueckner coupled cluster doubles with perturbative triples method, BCCD(T), was used to greatly improve upon CCSD(T), which was particularly problematic for CuNO. In the latter case the BCCD(T) vibrational frequencies showed significant differences compared to CCSD(T), e.g., nearly 65 cm -1 for the NO stretching frequency, and BCCD(T) also resulted in much better agreement with the available experimental frequencies. A full range of ro-vibrational spectroscopic constants are given for all three molecules of this study using the accurate composite potential energy functions and employing 2nd-order vibrational perturbation theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2021

The Story of 5d Metallocorroles: From Metal–Ligand Misfits to New Building Blocks for Cancer Phototherapeutics

Porphyrin chemistry is Shakespearean: over a century of study has not withered the field’s apparently infinite variety. Heme proteins continually astonish us with novel molecular mechanisms, while new porphyrin analogues bowl us over with unprecedented optical, electronic, and metal-binding properties. Within the latter domain, corroles occupy a special place, exhibiting a unique and rich coordination chemistry. The 5d metallocorroles are arguably the icing on that cake. New Zealand chemist Penny Brothers has used the word “misfit” to describe the interactions of boron, a small atom with a predilection for tetrahedral coordination, and porphyrins, classic square-planar ligands. Steve Jobs lionized misfits as those who see things differently and push humanity forward. Both perspectives have inspired us. The 5d metallocorroles are misfits in that they encapsulate a large 5d transition metal ion within the tight cavity of a contracted porphyrin ligand. Given the steric mismatch inherent in their structures, the syntheses of some 5d metallocorroles are understandably capricious, proceeding under highly specific conditions and affording poor yields. Three broad approaches may be distinguished. (a) In the metal–alkyl approach, a free-base corrole is exposed to an alkyllithium and the resulting lithio-corrole is treated with an early transition metal chloride; a variant of the method eschews alkyllithium and deploys a transition metal–alkyl instead, resulting in elimination of the alkyl group as an alkane and insertion of the metal into the corrole. This approach is useful for inserting transition metals from groups 4, 5, and, to some extent, 6, as well as lanthanides and actinides. (b) In our laboratory, we have often deployed a low-valent organometallic approach for the middle transition elements (groups 6, 7, 8, and 9). The reagents are low-valent metal–carbonyl or −olefin complexes, which lose one or more carbon ligands at high temperature, affording coordinatively unsaturated, sticky metal fragments that are trapped by the corrole nitrogens. (c) Finally, a metal acetate approach provides the method of choice for gold and platinum insertion (groups 10 and 11). This Account provides a first-hand perspective of the three approaches, focusing on the last two, which were largely developed in our laboratory. In general, the products were characterized with X-ray crystallography, electrochemistry, and a variety of spectroscopic methods. The physicochemical data, supplemented by relativistic DFT calculations, have provided fascinating insights into periodic trends and relativistic effects. An unexpected feature of many 5d metallocorroles, given their misfit character, is their remarkable stability under thermal, chemical, and photochemical stimulation. Many of them also exhibit long triplet lifetimes on the order of 100 μs and effectively sensitize singlet oxygen formation. Many exhibit phosphorescence in the near-infrared under ambient conditions. Furthermore, water-soluble ReO and Au corroles exhibit impressive photocytotoxicity against multiple cancer cell lines, promising potential applications as cancer phototherapeutics. We thus envision a bright future for the compounds as rugged building blocks for new generations of therapeutic and diagnostic (theranostic) agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Unsaturated Sulfur Crown Ethers Can Extract Mercury(II) and Show Promise for Future Copernicium(II) Studies: A Combined Experimental and Computational Study

The unsaturated hexathia-18-crown-6 (UHT18C6) molecule was investigated for the extraction of Hg(II) in HCl and HNO3 media. This extractant can be directly compared to the recently studied saturated hexathia-18-crown-6 (HT18C6). The default conformation of the S lone pairs in UHT18C6 is endodentate, where the pocket of the charge density, according to the crystal structures, is oriented toward the center of the ring, which should allow better extraction for Hg(II) compared to the exodentate HT18C6. Batch study experiments showed that Hg(II) had better extraction at low acid molarity (ca. 99% in HCl and ca. 95% in HNO 3 ), while almost no extraction was observed above 0.4 M HCl and 4 M HNO 3 (<5%). Speciation studies were conducted with the goal of delineating a plausible extraction mechanism. Density functional theory calculations including relativistic effects were carried out on both Hg(II)-encapsulated HT18C6 and UHT18C6 complexes to shed light on the binding strength and the nature of bonding. Our calculations offer insights into the extraction mechanism. In addition to Hg(II), calculations were performed on the hypothetical divalent Cn(II) ion, and showed that HT18C6 and UHT18C6 could extract Cn(II). Finally, the extraction kinetics were explored to assess whether this crown can extract the short-lived Cn(II) species in a future online experiment.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Linear versus Bent Uranium(II) Metallocenes—A Local Vibrational Mode Study

Uranium metallocenes have recently attracted attention driven by their use as catalysts in organometallic synthesis. In addition to bent U(IV) and U(III), an U(II) metallocene [(η 5 -C 5 i Pr 5 ) 2 U] was synthesized with an unusual linear Cp–U–Cp angle. Here, in this work, we investigated 22 U(II) metallocenes, (i) assessing the intrinsic strength of the U-ring interactions in these complexes with a novel bond strength measure based on our local vibrational mode analysis and (ii) systematically exploring what makes these U(II) metallocenes bent. We included relativistic effects through the NESCau Hamiltonian and complemented the local mode analysis with natural bonding orbital (NBO) and quantum theory of atoms in molecules (QTAIM) data. Our study led to the following results: (i) reduction of bulky U-ring ligand substituents does not lead to bent complexes for alkyl substituents ( i Pr and i Bu) in contrast to SiMe 3 ring substituents, which are all bent. (ii) The most bent complexes are [(η 5 -C 5 H 4 SiMe 3 ) 2 U] (130°) and [η 5 -P 5 H 5 ) 2 U] (143°). (iii) Linear complexes showed one hybridized NBO with s/d character, while bent structures were characterized by s/d/f mixing. (iv) We did not observe a correlation between the strength of the U-ring interaction and the amount of the ring-U-ring bend; the strongest interaction was found for [η 5 -Cp) 2 U] and the weakest for [η 5 -P 5 H 5 ) 2 U]. In conclusion, our results provide a foundation for the design of U(II) metallocenes with specific physicochemical properties and increased reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra

X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structure, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra is complicated by both relativistic effects, including spin-orbit coupling, and by the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space configuration interaction (X2C-MRRASCI) scheme that uses an exact two-component state average complete active space self-consistent-field method (X2C-SA-CASSCF), which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a restricted active space configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron containing complexes, with results showing key features of the spectrum being reproduced, including ligand to metal charge transfer and shake-up excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mapping the Electronic Structure of the Uranium(VI) Dinitride Molecule, UN 2

A combined anion photoelectron spectroscopic and relativistic electronic-structure computational study of the electronic structure of the UN 2 molecule is reported. Since the photoelectron spectrum of uranium dinitride’s negative ion, UN 2 - , (UN 2 plus an additional electron) directly reflects the electronic structure of neutral UN 2 , we have measured and relied upon the photoelectron spectrum of the UN 2 – anion as a means of mapping the electronic structure of neutral UN 2 . In addition to the electron affinity of the UN 2 ground state (the energy released when attaching an electron to UN 2 molecule), energy levels of the UN 2 excited states (electronic states with higher energies than the ground state) were well characterized by the close interplay between experiment and high-level theory. We found that both electron attachment and electronic excitation significantly bend the UN 2 molecule and elongate its U≡N bond. Implications for the activation of UN 2 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2020

Prediction of Structures and Atomization Energies of Coinage Metals, (M)n, n < 20: Extrapolation of Normalized Clustering Energies to Predict the Cohesive Energy

The geometries of the group 11 coinage metals (n = 2-20) were optimized to determine the lowest energy isomers for each cluster size, singlets for even numbers and doublets for odd numbers. For copper and silver, 2-D (planar) geometries were favored up to n = 6. For gold, 2D (planar) geometries were favored up to n = 13. Normalized clustering energies were plotted as a function of cluster size (n-1/3, for n = 4-20) with various DFT functionals and the CCSD(T)-F12b method and were extrapolated to predict the bulk cohesive energy. In the case of copper and silver, there is excellent agreement between the cohesive energies predicted at the CCSD(T)-F12b level of theory and the experimental values. For gold, the CCSD(T)-F12b values needed to be corrected for spin-orbit relativistic effects to obtain good agreement with experiment. Electronic properties including the HOMO-LUMO gaps for the even clusters and the spin densities for the odd clusters were calculated. The lowest gap is predicted to occur for n = 16 where the HOMO and LUMO are very similar in shape.

Persuad, Rudradatt R.↗
DOE OSTI2023

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Creation of an unexpected plane of enhanced covalency in cerium(III) and berkelium(III) terpyridyl complexes

Abstract Controlling the properties of heavy element complexes, such as those containing berkelium, is challenging because relativistic effects, spin-orbit and ligand-field splitting, and complex metal-ligand bonding, all dictate the final electronic states of the molecules. While the first two of these are currently beyond experimental control, covalent M‒L interactions could theoretically be boosted through the employment of chelators with large polarizabilities that substantially shift the electron density in the molecules. This theory is tested by ligating Bk III with 4’-(4-nitrophenyl)-2,2’:6’,2”-terpyridine (terpy*), a ligand with a large dipole. The resultant complex, Bk(terpy*)(NO 3 ) 3 (H 2 O)·THF, is benchmarked with its closest electrochemical analog, Ce(terpy*)(NO 3 ) 3 (H 2 O)·THF. Here, we show that enhanced Bk‒N interactions with terpy* are observed as predicted. Unexpectedly, induced polarization by terpy* also creates a plane in the molecules wherein the M‒L bonds trans to terpy* are shorter than anticipated. Moreover, these molecules are highly anisotropic and rhombic EPR spectra for the Ce III complex are reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Deep-learning atomistic semi-empirical pseudopotential model for nanomaterials

The semi-empirical pseudopotential method (SEPM) has been widely applied to provide computational insights into the electronic structure, photophysics, and charge carrier dynamics of nanoscale materials. We present “DeepPseudopot”, a machine-learned atomistic pseudopotential model that extends the SEPM framework by combining a flexible neural network representation of the local pseudopotential with parameterized non-local and spin-orbit coupling terms. Trained on bulk quasiparticle band structures and deformation potentials from GW calculations, the model captures many-body and relativistic effects with very high accuracy across diverse semiconducting materials, as illustrated for silicon and group III-V semiconductors. DeepPseudopot’s accuracy, efficiency, and transferability make it well-suited for data-driven in silico design and discovery of novel optoelectronic nanomaterials.

Lin, Kailai [University of California, Berkeley, C↗
DOE OSTI2022

The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗