Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Containers”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Confocal Raman imaging of skin sections containing hair follicles using classical least squares regression and multivariate curve resolution – alternating least squares

Confocal Raman microscopy (CRM) is applied ex vivo for imaging of the spatial distribution of different skin components in skin sections containing hair follicles. For multivariate data analysis, different methods are used in order to spectrally decompose the reference spectra of the skin components (dermis, viable epidermis, stratum corneum and hair). Classical least squares regression (CLS) and multivariate curve resolution–alternating least squares (MCR-ALS) are chosen as suitable methods. In comparison to other optical methods, the advantage of CRM is molecular specificity and dispensability of labelling dyes, which is e.g. necessary in fluorescence microscopy. Therefore, a useful future application of CRM in combination with multivariate data analysis lies in the analysis of penetration routes of topically applied substances, such as cosmetic formulations or drugs into the skin, which is particularly interesting in and around hair follicles. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Effect of laser radiation on the gamma activity of aqueous salt solutions containing {sup 152}Eu

We report an experimental study of the effect of laser irradiation of aqueous salt solutions containing the {sup 152}Eu isotope in the presence of nickel nanoparticles on the gamma activity of the isotope. Nickel nanoparticles are generated using laser ablation of the target in water and added to the radionuclide solution before irradiation. Two types of Nd : YAG lasers with a pulse duration of 10 ns and different repetition rates are used to irradiate the solutions. It is found that the {sup 152}Eu activity as a result of repeated irradiation of solutions by a series of pulses decreases with each irradiation cycle. The obtained results are compared with data on the effect of maser radiation with a wavelength of 8 mm on the activity of the same isotope. Possible mechanisms of the effect of electromagnetic radiation on the {sup 152}Eu gamma activity are discussed. (interaction of laser radiation with matter)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structural basis for differential recognition of phosphohistidine-containing peptides by 1-pHis and 3-pHis monoclonal antibodies

In 2015, monoclonal antibodies (mAbs) that selectively recognize the 1-pHis or 3-pHis isoforms of phosphohistidine were developed by immunizing rabbits with degenerate Ala/Gly peptides containing the nonhydrolyzable phosphohistidine (pHis) analog- phosphotriazolylalanine (pTza). Here, we report structures of five rabbit mAbs bound to cognate pTza peptides: SC1-1 and SC50-3 that recognize 1-pHis, and their 3-pHis–specific counterparts, SC39-4, SC44-8, and SC56-2. These cocrystal structures provide insights into the binding modes of the pTza phosphate group that are distinct for the 1- and 3-pHis mAbs with the selectivity arising from specific contacts with the phosphate group and triazolyl ring. The mode of phosphate recognition in the 3-pHis mAbs recapitulates the Walker A motif, as present in kinases. The complementarity-determining regions (CDRs) of four of the Fabs interact with the peptide backbone rather than peptide side chains, thus conferring sequence independence, whereas SC44-8 shows a proclivity for binding a GpHAGA motif mediated by a sterically complementary CDRL3 loop. Specific hydrogen bonding with the triazolyl ring precludes recognition of pTyr and other phosphoamino acids by these mAbs. Kinetic binding experiments reveal that the affinity of pHis mAbs for pHis and pTza peptides is submicromolar. Bound pHis mAbs also shield the pHis peptides from rapid dephosphorylation. The epitope–paratope interactions illustrate how these anti-pHis antibodies are useful for a wide range of research techniques and this structural information can be utilized to improve the specificity and affinity of these antibodies toward a variety of pHis substrates to understand the role of histidine phosphorylation in healthy and diseased states.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Family-wide Characterization of Histone Binding Abilities of Human CW Domain-containing Proteins

Covalent modifications of histone N-terminal tails play a critical role in regulating chromatin structure and controlling gene expression. These modifications are controlled by histone-modifying enzymes and read out by histone-binding proteins. Numerous proteins have been identified as histone modification readers. Here we report the family-wide characterization of histone binding abilities of human CW domain-containing proteins. We demonstrate that the CW domains in ZCWPW2 and MORC3/4 selectively recognize histone H3 trimethylated at Lys-4, similar to ZCWPW1 reported previously, while the MORC1/2 and LSD2 lack histone H3 Lys-4 binding ability. Our crystal structures of the CW domains of ZCWPW2 and MORC3 in complex with the histone H3 trimethylated at Lys-4 peptide reveal the molecular basis of this interaction. In each complex, two tryptophan residues in the CW domain form the “floor” and “right wall,” respectively, of the methyllysine recognition cage. Our mutation results based on ZCWPW2 reveal that the right wall tryptophan residue is essential for binding, and the floor tryptophan residue enhances binding affinity. Here, our structural and mutational analysis highlights the conserved roles of the cage residues of CW domain across the histone methyllysine binders but also suggests why some CW domains lack histone binding ability.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Wavefunction stability analysis without analytical electronic Hessians: application to orbital-optimised second-order Møller–Plesset theory and VV10-containing density functionals

Wavefunction stability analysis is commonly applied to converged self-consistent field (SCF) solutions to verify whether the electronic energy is a local minimum with respect to second-order variations in the orbitals. By iterative diagonalisation, the procedure calculates the lowest eigenvalue of the stability matrix or electronic Hessian. However, analytical expressions for the electronic Hessian are unavailable for most advanced post-Hartree-Fock (HF) wave function methods and even some Kohn-Sham (KS) density functionals. To address such cases, here we formulate the Hessian-vector product within the iterative diagonalisation procedure as a finite difference of the electronic gradient with respect to orbital perturbations in the direction of the vector. As a model application, following the lowest eigenvalue of the orbital-optimised second-order Møller-Plesset perturbation theory (OOMP2) Hessian during H 2 dissociation reveals the surprising stability of the spin-restricted solution at all separations, with a second independent unrestricted solution. We show that a single stable solution can be recovered by using the regularised OOMP2 method (δ-OOMP2), which contains a level shift. Internal and external stability analyses are also performed for SCF solutions of a recently developed range-separated hybrid density functional, ωB97X-V, for which the analytical Hessian is not yet available due to the complexity of its long-range non-local VV10 correlation functional.

Brueckner orbitals↗
DOE OSTI2013

CHARACTERISTIC LENGTH OF ENERGY-CONTAINING STRUCTURES AT THE BASE OF A CORONAL HOLE

An essential parameter for models of coronal heating and fast solar wind acceleration that rely on the dissipation of MHD turbulence is the characteristic energy-containing length {lambda} of the squared velocity and magnetic field fluctuations (u{sup 2} and b{sup 2}) transverse to the mean magnetic field inside a coronal hole (CH) at the base of the corona. The characteristic length scale directly defines the heating rate. We use a time series analysis of solar granulation and magnetic field measurements inside two CHs obtained with the New Solar Telescope at Big Bear Solar Observatory. A data set for transverse magnetic fields obtained with the Solar Optical Telescope/Spectro-Polarimeter on board the Hinode spacecraft was utilized to analyze the squared transverse magnetic field fluctuations b{sub t}{sup 2}. Local correlation tracking was applied to derive the squared transverse velocity fluctuations u {sup 2}. We find that for u {sup 2} structures, the Batchelor integral scale {lambda} varies in a range of 1800-2100 km, whereas the correlation length sigmav and the e-folding length L vary between 660 and 1460 km. Structures for b{sub t}{sup 2} yield {lambda} Almost-Equal-To 1600 km, sigmav Almost-Equal-To 640 km, and L Almost-Equal-To 620 km. An averaged (over {lambda}, sigmav, and L) value of the characteristic length of u {sup 2} fluctuations is 1260 {+-} 500 km, and that of b{sub t}{sup 2} is 950 {+-} 560 km. The characteristic length scale in the photosphere is approximately 1.5-50 times smaller than that adopted in previous models (3-30 Multiplication-Sign 10{sup 3} km). Our results provide a critical input parameter for current models of coronal heating and should yield an improved understanding of fast solar wind acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Gas-phase reactions of polycyclic aromatic hydrocarbon cations and their nitrogen-containing analogs with H atoms

We have studied the reactions of polycyclic aromatic hydrocarbon cations and their nitrogen-containing analogs with H atoms. Reaction rate constants are measured at 300 K using a flowing afterglow-selected ion flow tube. We have implemented the laser induced acoustic desorption technique to allow the study of large, non-volatile species in the gas phase. The extension of this work from previous studies shows that the reactivity of polycyclic aromatic hydrocarbon cations with H atoms reaches a constant value for large cations. There is a small difference in reactivity when comparing molecules of different size and geometry; however, no difference in reactivity was found when nitrogen was incorporated into the ring.

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

FORMATION OF S-BEARING SPECIES BY VUV/EUV IRRADIATION OF H{sub 2}S-CONTAINING ICE MIXTURES: PHOTON ENERGY AND CARBON SOURCE EFFECTS

Carbonyl sulfide (OCS) is a key molecule in astrobiology that acts as a catalyst in peptide synthesis by coupling amino acids. Experimental studies suggest that hydrogen sulfide (H{sub 2}S), a precursor of OCS, could be present in astrophysical environments. In the present study, we used a microwave-discharge hydrogen-flow lamp, simulating the interstellar UV field, and a monochromatic synchrotron light beam to irradiate CO:H{sub 2}S and CO{sub 2}:H{sub 2}S ice mixtures at 14 K with vacuum ultraviolet (VUV) or extreme ultraviolet (EUV) photons in order to study the effect of the photon energy and carbon source on the formation mechanisms and production yields of S-containing products (CS{sub 2}, OCS, SO{sub 2}, etc.). Results show that (1) the photo-induced OCS production efficiency in CO:H{sub 2}S ice mixtures is higher than that of CO{sub 2}:H{sub 2}S ice mixtures; (2) a lower concentration of H{sub 2}S enhances the production efficiency of OCS in both ice mixtures; and (3) the formation pathways of CS{sub 2} differ significantly upon VUV and EUV irradiations. Furthermore, CS{sub 2} was produced only after VUV photoprocessing of CO:H{sub 2}S ices, while the VUV-induced production of SO{sub 2} occurred only in CO{sub 2}:H{sub 2}S ice mixtures. More generally, the production yields of OCS, H{sub 2}S{sub 2}, and CS{sub 2} were studied as a function of the irradiation photon energy. Heavy S-bearing compounds were also observed using mass spectrometry during the warm-up of VUV/EUV-irradiated CO:H{sub 2}S ice mixtures. The presence of S-polymers in dust grains may account for the missing sulfur in dense clouds and circumstellar environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

N- AND O-HETEROCYCLES PRODUCED FROM THE IRRADIATION OF BENZENE AND NAPHTHALENE IN H{sub 2}O/NH{sub 3}-CONTAINING ICES

Aromatic heterocyclic molecules are an important class of molecules of astrophysical and biological significance that include pyridine, pyrimidine, and their derivatives. Such compounds are believed to exist in interstellar and circumstellar environments, though they have never been observed in the gas phase. Regardless of their presence in the gas phase in space, numerous heterocycles have been reported in carbonaceous meteorites, which indicates that they are formed under astrophysical conditions. The experimental work described here shows that N- and O-heterocyclic molecules can form from the ultraviolet (UV) irradiation of the homocyclic aromatic molecules benzene (C{sub 6}H{sub 6}) or naphthalene (C{sub 10}H{sub 8}) mixed in ices containing H{sub 2}O and NH{sub 3}. This represents an alternative way to generate aromatic heterocycles to those considered before and may have important implications for astrochemistry and astrobiology.

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

A unified energy-reservoir model containing contributions from {sup 56}Ni and neutron stars and its implication for luminous type Ic supernovae

Most type Ic core-collapse supernovae produce {sup 56}Ni and neutron stars (NSs) or black holes. The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitudes M{sub peak} in any band ≳−19.5 mag, while the {sup 56}Ni can be neglected in fitting most type Ic superluminous supernovae whose M{sub peak} in any band is ≲−21 mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of {sup 56}Ni. For luminous SNe Ic with −21≲M{sub peak}≲−19.5 mag, however, both contributions from {sup 56}Ni and NSs cannot be neglected without serious modeling since they are not SLSNe and the {sup 56}Ni mass could be up to ∼0.5 M{sub ⊙}. In this paper we propose a unified model that contains contributions from both {sup 56}Ni and a nascent NS. We select three luminous SNe Ic-broad-lined, SN 2010ay, SN 2006nx, and SN 14475, and show that, if these SNe are powered by {sup 56}Ni, the ratio of M{sub Ni} to M{sub ej} is unrealistic. Alternatively, we invoke the magnetar model and the hybrid ({sup 56}Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on {sup 56}Ni yields and the parameters of putative magnetars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Plasma diagnostics and modeling of lithium-containing plasmas

Thin-film deposition from chemically reactive multi-component plasmas is complex, and the lack of electron collision cross-sections for even the most common metalorganic precursors and their fragments complicates their modeling based on fundamental plasma physics. This study focuses on understanding the plasma physics and chemistry in argon (Ar) plasmas containing lithium bis (trimethylsilyl) amide used to deposit Li x Si y thin films. These films are emerging as potential solid electrolytes for lithium-ion batteries, and the Li-to-Si ratio is a crucial parameter to enhance their ionic conductivity. We deposited Li x Si y films in an axial flow-through plasma reactor and studied the factors that determine the variation of the Li-to-Si ratio in films deposited at various points on a substrate spanning the entire reactor axis. While the Li-to-Si ratio is 1:2 in the precursor, the Li-to-Si ratio is as high as 3:1 in films deposited near the plasma entrance and decreases to 1:1 for films deposited downstream. Optical emission from the plasma is dominated by Li emission near the entrance, but Li emission disappears downstream, which we attribute to the complete consumption of the precursor. We hypothesized that the axially decreasing precursor concentration affects the electron energy distribution function in a way that causes different dissociation efficiencies for the production of Li and Si. We used Li line intensities to estimate the local precursor concentration and Ar line ratios to estimate the local reduced electric field to test this hypothesis. This analysis suggests that the mean electron energy increases along the reactor axis with decreasing precursor concentration. As a result, the decreasing Li-to-Si ratio with axially decreasing precursor concentration may be explained by Li release from the precursor having lower threshold energy than Si release.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The role of isotope mass and transport for H-mode access in tritium containing plasmas at JET with ITER-like wall

The required heating power, $P_{\mathrm{LH}}$, to access the high confinement regime (H-mode) in tritium containing plasmas is investigated in JET with ITER-like wall at a toroidal magnetic field of $B_t = 1.8$ T and a plasma current of $I_p = 1.7$ MA. $P_{\mathrm{LH}}$, also referred to as the L-H power threshold, is determined in plasmas of pure tritium as well as mixtures of hydrogen with tritium (H-T) and mixtures of deuterium with tritium (D-T), and is compared to the L-H power threshold in plasmas of pure hydrogen and pure deuterium. It is found that, for otherwise constant parameters, $P_{\mathrm{LH}}$ is not the same in plasmas with the same effective isotope mass, $A_{\mathrm{eff}}$, when they differ in their isotope composition. Thus, $A_{\mathrm{eff}}$ is not sufficient to describe the isotope effect of $P_{\mathrm{LH}}$ in a consistent manner for all considered isotopes and isotope mixtures. The electron temperature profiles measured at the L-H transition in the outer half of the radius are very similar for all isotopes and isotope mixtures, despite the fact that the L-H power threshold varies by a factor of about six. This finding, together with the observation of an offset linear relation between the L-H power threshold, $P_{\mathrm{LH}}$, and an effective heat diffusivity, $\chi_{\mathrm{eff}}$, indicates that the composition-dependent heat transport in the low confinement mode (L-mode) determines, how much power is needed to reach the necessary electron temperatures at the edge, and hence $P_{\mathrm{LH}}$.

tritium plasmas↗
DOE OSTI2021

Vertex-finding and reconstruction of contained two-track neutrino events in the MicroBooNE detector

In this work, we describe algorithms developed to isolate and accurately reconstruct two-track events that are contained within the MicroBooNE detector. This method is optimized to reconstruct two tracks of lengths longer than 5cm. This code has applications to searches for neutrino oscillations and measurements of cross sections using quasi-elastic-like charged current events. The algorithms we discuss will be applicable to all detectors running in Fermilab's Short Baseline Neutrino program (SBN), and to any future liquid argon time projection chamber (LArTPC) experiment with beam energies ~ 1 GeV. The algorithms are publicly available on a GITHUB repository. This reconstruction offers a complementary and independent alternative to the Pandora reconstruction package currently in use in LArTPC experiments, and provides similar reconstruction performance for two-track events.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Stable isotopes contain substantial additive information about terrestrial carbon and water cycling

Abstract Stable isotope ratios of H ( δ 2 H ), O ( δ 18 O), and C ( δ 13 C) are linked to key biogeochemical processes of the water and carbon cycles; however, the degree to which isotope-associated processes are reflected in macroscale ecosystem flux observations remains unquantified. Here through formal information assessment, new measurements of δ 13 C of net ecosystem exchange ( NEE ) as well as δ 2 H and δ 18 O of latent heat ( LH ) fluxes across the United States National Ecological Observation Network (NEON) are used to determine conditions under which isotope measurements are informative of environmental exchanges. We find all three isotopic datasets individually contain comparable amounts of information about NEE and LH fluxes as wind speed observations. Such information from isotope measurements, however, is largely unique. Generally, δ 13 C provides more information about LH as aridity increases or mean annual precipitation decreases. δ 2 H provides more information about LH as temperatures or mean annual precipitation decreases, and also provides more information about NEE as temperatures decrease. Overall, we show that the stable isotope datasets collected by NEON contribute non-trivial amounts of new information about bulk environmental fluxes useful for interpreting biogeochemical and ecohydrological processes at landscape scales. However, the utility of this new information varies with environmental conditions at continental scales. This study provides an approach for quantifying the value adding non-traditional sensing approaches to environmental monitoring sites and the patterns identified here are expected to aid in modeling and data interpretation efforts focused on constraining carbon and water cycles’ mechanisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Crystal structure of a DNA/Ba 2+ G-quadruplex containing a water-mediated C-tetrad

We have determined the 1.50 Å crystal structure of the DNA decamer, d(CCA CNV KGCGTGG) ( CNV K, 3-cyanovinylcarbazole), which forms a G-quadruplex structure in the presence of Ba 2+ . The structure contains several unique features including a bulged nucleotide and the first crystal structure observation of a C-tetrad. The structure reveals that water molecules mediate contacts between the divalent cations and the C-tetrad, allowing Ba 2+ ions to occupy adjacent steps in the central ion channel. One ordered Mg 2+ facilitates 3'-3' stacking of two quadruplexes in the asymmetric unit, while the bulged nucleotide mediates crystal contacts. Despite the high diffraction limit, the first four nucleotides including the CNV K nucleoside are disordered though they are still involved in crystal packing. This work suggests that the bulky hydrophobic groups may locally influence the formation of non-Watson–Crick structures from otherwise complementary sequences. These observations lead to the intriguing possibility that certain types of DNA damage may act as modulators of G-quadruplex formation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

hLARP7 C-terminal domain contains an xRRM that binds the 3' hairpin of 7SK RNA

The 7SK small nuclear ribonucleoprotein (snRNP) sequesters and inactivates the positive transcription elongation factor b (P-TEFb), an essential eukaryotic mRNA transcription factor. The human La-related protein group 7 (hLARP7) is a constitutive component of the 7SK snRNP and localizes to the 3' terminus of the 7SK long noncoding RNA. hLARP7, and in particular its C-terminal domain (CTD), is essential for 7SK RNA stability and assembly with P-TEFb. The hLARP7 N-terminal Lamodule binds and protects the 3' end from degradation, but the structural and functional role of its CTD is unclear.We report the solution NMR structure of the hLARP7 CTD and show that this domain contains an xRRM, a class of atypical RRM first identified in the Tetrahymena thermophila telomerase LARP7 protein p65. The xRRM binds the 3' end of 7SK RNA at the top of stem-loop 4 (SL4) and interacts with both unpaired and base-paired nucleotides. This study thus confirms that the xRRM is general to the LARP7 family of proteins and defines the binding site for hLARP7 on the 7SK RNA, providing insight into function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Structural basis for preferential binding of human TCF4 to DNA containing 5-carboxylcytosine

The psychiatric risk-associated transcription factor 4 (TCF4) is linked to schizophrenia. Rare TCF4 coding variants are found in individuals with Pitt-Hopkins syndrome—an intellectual disability and autism spectrum disorder. TCF4 contains a C-terminal basic-helix-loop-helix (bHLH) DNA binding domain which recognizes the enhancer-box (E-box) element 5'-CANNTG-3' (where N = any nucleotide). A subset of the TCF4-occupancy sites have the expanded consensus binding specificity 5'-C(A/G)-CANNTG-3', with an added outer Cp(A/G) dinucleotide; for example in the promoter for CNIH3, a gene involved in opioid dependence. In mammalian genomes, particularly brain, the CpG and CpA dinucleotides can be methylated at the 5-position of cytosine (5mC), and then may undergo successive oxidations to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), and 5-carboxyl (5caC) forms. We find that, in the context of 5'- 0 CG- 1 CA- 2 CG- 3 TG-3'(where the numbers indicate successive dinucleotides), modification of the central E-box 2 CG has very little effect on TCF4 binding, E-box 1 CA modification has a negative influence on binding, while modification of the flanking 0 CG, particularly carboxylation, has a strong positive impact on TCF4 binding to DNA. Crystallization of TCF4 in complex with unmodified or 5caC-modified oligonucleotides revealed that the basic region of bHLH domain adopts multiple conformations, including an extended loop going through the DNA minor groove, or the N-terminal portion of a long helix binding in the DNA major groove. The different protein conformations enable arginine 576 (R576) to interact, respectively, with a thymine in the minor groove, a phosphate group of DNA backbone, or 5caC in the major groove. The Pitt-Hopkins syndrome mutations affect five arginine residues in the basic region, two of them (R569 and R576) involved in 5caC recognition. Our analyses indicate, and suggest a structural basis for, the preferential recognition of 5caC by a transcription factor centrally important in brain development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Dirac's equation and its implications for density functional theory based calculations of materials containing heavy elements

Electronic structure calculations based on density functional theory (DFT) give quantitatively accurate predictions of properties of most materials containing light elements. For heavy materials, and in particular for f -electron systems, DFT based methods can fail both qualitatively and quantitatively for two distinct reasons: their failure to describe confinement effects arising from localized f -electron behavior and their incomplete or approximate treatment of relativity. In addition, different methods for incorporating relativistic effects, which give identical results in most light materials, can give different predictions in heavy elements. In order to develop a quantitative capability for calculating the properties of these materials, it is essential to separate the predictions of the underlying equations from the uncertainty introduced in approximations used in computation. Working toward that goal, here we have developed a code, called dirac-fp , which is based directly on solving the Dirac-Kohn-Sham equations and uses the full-potential linear muffin-tin orbital (FP-LMTO) approach to electronic structure. In order to assess the performance of dirac-fp , we perform calculations on three different face-centered cubic materials using different approximate treatments of relativity: the scalar relativistic (SR) approach commonly used in most solid-state DFT codes, the scalar relativistic plus spin-orbit coupling corrections (SR+SO) approach which includes spin-orbit coupling self-consistently using the SR states inside the muffin tins, and the Dirac-Kohn-Sham (Dirac) approach implemented in dirac-fp . Performing calculations on thorium, in which relativistic effects should be strong, aluminum, in which relativistic effects should be negligible, and gold, in which relativistic effects play an intermediate role, we find that the Dirac approach is able to provide theoretically consistent results in the electronic structure and ground-state properties across all three materials.

3-dimensional systems↗