Four-component relativistic density functional calculations of heavy diatomic molecules

Toward a correct treatment of core properties with local hybrid functionals

In local hybrid functionals (LHs), a local mixing function (LMF) determines the position-dependent exact-exchange admixture. We report new LHs that focus on an improvement of the LMF in the core region while retaining or partly improving upon the high accuracy in the valence region exhibited by the LH20t functional. The suggested new pt-LMFs are based on a Padé form and modify the previously used ratio between von Weizsäcker and Kohn-Sham local kinetic energies by different powers of the density to enable flexibly improved approximations to the correct high-density and iso-orbital limits relevant for the innermost core region. Using TDDFT calculations for a set of K-shell core excitations of second- and third-period systems including accurate state-of-the-art relativistic orbital corrections, the core part of the LMF is optimized, while the valence part is optimized as previously reported for test sets of atomization energies and reaction barriers (Haasler et al., J Chem Theory Comput 2020, 16, 5645). The LHs are completed by a calibration function that minimizes spurious nondynamical correlation effects caused by the gauge ambiguities of exchange-energy densities, as well as by B95c meta-GGA correlation. The resulting LH23pt functional relates to the previous LH20t functional but specifically improves upon the core region.

Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques

Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.

Interactions of Pt nanoparticles with molecular components in polymer electrolyte membrane fuel cells: multi-scale modeling approach

Three phase model consists of Pt nanoparticles, Nafion, and graphite with oxygen, water, and hydronium.

Potential energy curves for the ground and low-lying excited states of CuAg

The ground and low-lying excited states of heteronuclear diatomic CuAg are examined by multi-reference configuration interaction (MRCI) method. Relativistic effects were treated and probed in two steps. Scalar terms were considered using the spin-free DKH Hamiltonian as a priori and spin-orbit coupling was calculated perturbatively via the spin-orbit terms of the Breit-Pauli Hamiltonian based on MRCI wavefunctions. Potential energy curves of the spin-free states and their corresponding Ω components correlating with the separated atom limits (2)S(Cu) + (2)S(Ag) and (2)D(Cu) + (2)S(Ag) are obtained. The results are in fine agreement with the experimental measurements and tentative conclusions for the ion-pair B0(+) state are confirmed by our theoretical calculations. Illustrative results are presented to reveal the relative importance and magnitude of the scalar and spin-orbit effects on the spectroscopic properties of this molecule. Time dependent density functional theory calculations, using the LDA, BLYP, B3LYP, and SAOP functionals have been carried out for CuAg and the accuracy of TD-DFT has been compared with ab initio results.

Relativistic and electron-correlation effects on the nuclear magnetic resonance shieldings of molecules containing tin and lead atoms

The reference values for NMR magnetic shieldings, σ(ref), are of the highest importance when theoretical analysis of chemical shifts are envisaged. The fact that the nonrelativistically valid relationship among spin-rotation constants and magnetic shieldings is not any longer valid for heavy atoms requires that the search for σ(ref) for such atoms needs new strategies to follow. We present here results of σ(ref) that were obtained by applying our own simple procedure which mixes accurate experimental chemical shifts (δ) and theoretical magnetic shieldings (σ). We calculated σ(Sn) and σ(Pb) in a family of heavy-halogen-containing molecules. We found out that σ(ref)[Sn;Sn(CH3)4] in gas phase should be close to 3864.11 ± 20.05 ppm (0.5%). For Pb atom, σ(ref)[Pb;Pb(CH3)4] should be close to 14475.1 ± 500.7 ppm. Such theoretical values correspond to calculations with the relativistic polarization propagator method, RelPPA, at the RPA level of approach. They are closer to experimental values as compared to those obtained applying few different functionals such as PBE0, B3LYP, BLYP, BP86, KT2, and KT3 of the density functional theory, DFT. We studied tin and lead shieldings of the XY(4-n)Z(n) (X = Sn, Pb; Y, Z = H, F, Cl, Br, I) and PbH(4-n)I(n) (n = 0, 1, 2, 3, 4) family of compounds with four-component functionals as implemented in the DIRAC code. For these systems results of calculations with RelPPA-RPA are more reliable than DFT ones. We argue about why those DFT functionals must be modified in order to obtain more accurate results of NMR magnetic shieldings within the relativistic regime: first, there is a dependence among both electron-correlation and relativistic effects that should be introduced in some way in the functionals; and second, the DIRAC code uses standard nonrelativistic functionals and the functionals B3LYP and PBE0 were parametrized only with data taken from light elements. It can explain why they are not able to properly introduce relativistic effects on nuclear magnetic shieldings. We finally show that in the analysis of magnetic shieldings for the family of compounds mentioned above, one must consider the newest and so-called heavy-atom effect on vicinal heavy atoms, HAVHA. Such effects are among the most important relativistic effects in these kind of compounds.

Magnetic anisotropy of heteronuclear dimers in the gas phase and supported on graphene: relativistic density-functional calculations

The structural and magnetic properties of mixed PtCo, PtFe, and IrCo dimers in the gas phase and supported on a free-standing graphene layer have been calculated using density-functional theory, both in the scalar-relativistic limit and self-consistently including spin-orbit coupling. The influence of the strong magnetic moments of the 3d atoms on the spin and orbital moments of the 5d atoms, and the influence of the strong spin-orbit coupling contributed by the 5d atom on the orbital moments of the 3d atoms have been studied in detail. The magnetic anisotropy energy is found to depend very sensitively on the nature of the eigenstates in the vicinity of the Fermi level, as determined by band filling, exchange splitting and spin-orbit coupling. The large magnetic anisotropy energy of free PtCo and IrCo dimers relative to the easy direction parallel to the dimer axis is coupled to a strong anisotropy of the orbital magnetic moments of the Co atom for both dimers, and also on the Ir atom in IrCo. In contrast the PtFe dimer shows a weak perpendicular anisotropy and only small spin and orbital anisotropies of opposite sign on the two atoms. For dimers supported on graphene, the strong binding within the dimer and the stronger interaction of the 3d atom with the substrate stabilizes an upright geometry. Spin and orbital moments on the 3d atom are strongly quenched, but due to the weaker binding within the dimer the properties of the 5d atom are more free-atom-like with increased spin and orbital moments. The changes in the magnetic moment are reflected in the structure of the electronic eigenstates near the Fermi level, for all three dimers the easy magnetic direction is now parallel to the dimer axis and perpendicular to the graphene layer. The already very large magnetic anisotropy energy (MAE) of IrCo is further enhanced by the interaction with the support, the MAE of PtFe changes sign, and that of the PtCo dimer is reduced. These changes are discussed in relation to the relativistic electronic structure of free and supported dimers and it is demonstrated that the existence of a partially occupied quasi-degenerate state at the Fermi level favors the formation of a large magnetic anisotropy.

The Beijing Density Functional (BDF) Program Package: Methodologies and Applications

The Beijing Density Functional (BDF) program package is such a code that can perform nonrelativistic, one-, two-, and four-component relativistic density functional calculations on medium-sized molecular systems with various functionals in most compact and yet sufficient basis set expansions. The mergence of different approaches in a single code facilitates direct and systematic comparisons between different Hamiltonians, since they share all the same numerical and technical issues. In this account, the methodologies adopted in the code will be discussed in great detail and some applications of the code will be briefly presented.

Infrared spectra of the (AgCO)2 and AgnCO (n=2-4) molecules in rare-gas matrices

Reactions of laser-ablated silver atoms with carbon monoxide molecules in solid argon and neon have been investigated using matrix-isolation IR spectroscopy. Small silver cluster carbonyls, (AgCO)2 and AgnCO (n=2-4), as well as mononuclear silver carbonyls, Ag(CO)2 and Ag(CO)3, are generated upon sample annealing in the argon experiments and are characterized on the basis of the isotopic substitution, the CO concentration change, and the comparison with theoretical predictions. However, these polynuclear carbonyls are absent from the neon experiments. Density functional theory calculations have been performed on these silver carbonyls and the corresponding ligand-free silver clusters, which support the identification of these silver carbonyls from the matrix IRspectrum. A terminal CO has been found in the most stable structures of (AgCO)2, Ag2CO, Ag3CO, and Ag4CO. Furthermore, a plausible reaction mechanism has been proposed to account for the formation of the (AgCO)2 and AgnCO (n=2-4) molecules.

Adsorption of superheavy elements on metal surfaces

Fully relativistic four-component density functional theory with the general gradient approximation calculations have been performed to determine the adsorption energy and position of the superheavy element 112 on a Au surface. Extended cluster as well as embedded cluster calculations were used to simulate the surface which allow for the top, bridge, and hollow adsorption sites without losing the advantage of symmetry considerations. Comparison with analogous calculations of the adsorption of the homologue element Hg allows to predict the adsorption of element 112 at a bridge site with a binding energy of 0.67 eV.

Relativistic density functional calculations using two-spinor minimax finite-element method and linear combination of atomic orbitals for ZnO, CdO, HgO, UubO and Cu2, Ag2, Au2, Rg2

In previous work the authors have presented a highly accurate two-spinor fully relativistic solution of the two-center Coulomb problem utilizing the finite-element method (FEM) and furthermore developed a relativistic minimax two-spinor linear combination of atomic orbitals (LCAO). In the present paper the authors present Dirac-Fock-Slater (DFS-) density functional calculations for two-atomic molecules up to super heavy systems using the fully nonlinear minimax FEM and the minimax LCAO in its linearized approximation (linear approximation to relativistic minimax). The FEM gives highly accurate benchmark results for the DFS functional. Especially considering molecules with up to super heavy atoms such as UubO and Rg2, the authors found that LCAO fails to give the correct systematic trends. The accurate FEM results shed a new light on the quality of the DFS-density functional.

Exploring the Hydration of Pb2+: Ab Initio Studies and First-Principles Molecular Dynamics

Even though lead is a well-known toxicant widely scattered throughout the world since antiquity, its chemistry is poorly documented at the molecular level. Here we investigate the hydration of the Pb(2+) ion by means of first-principles molecular dynamics (Car-Parrinello molecular dynamics, CPMD). We found that the hydrated cation is heptacoordinated in a dynamically holodirected arrangement roughly corresponding to a fluxional distorted pentagonal bipyramid. The time-averaged Pb-O bond length is especially large and amounts to 2.70 A with an associated root-mean-square deviation of 0.26 A. This results from a dynamic exchange between short (<2.6 A), intermediate (2.6-3.0 A) and long (>3.0 A) Pb-O bonds. The latter very long Pb-O distance implies that the determination of the coordination number n(c) from experimental work may not necessarily yield values directly comparable to the theoretical value of n(c)=7, since not all experimental techniques would recognize such a long distance as a bond to the metal cation. Pronounced disorders are evidenced in the second shell, characteristic of a chaotropic cation, and exchanges between the first and second shells cannot be excluded on a timescale of a few tens of picoseconds.

[Pb(H2O)]2+ and [Pb(OH)]+: Four-component density functional theory calculations, correlated scalar relativistic constrained-space orbital variation energy decompositions, and topological analysis

Within the scope of studying the molecular implications of the Pb(2+) cation in environmental and polluting processes, this paper reports Hartree-Fock and density functional theory (B3LYP) four-component relativistic calculations using an all-electron basis set applied to [Pb(H(2)O)](2+) and [Pb(OH)](+), two complexes expected to be found in the terrestrial atmosphere. It is shown that full-relativistic calculations validate the use of scalar relativistic approaches within the framework of density functional theory. [Pb(H(2)O)](2+) is found C(2v) at any level of calculations whereas [Pb(OH)](+) can be found bent or linear depending of the computational methodology used. When C(s) is found the barrier to inversion through the C(infinityv) structure is very low, and can be overcome at high enough temperature, making the molecule floppy. In order to get a better understanding of the bonding occurring between the Pb(2+) cation and the H(2)O and OH(-) ligands, natural bond orbital and atoms-in-molecule calculations have been performed. These approaches are supplemented by a topological analysis of the electron localization function. Finally, the description of these complexes is refined using constrained-space orbital variation complexation energy decompositions.

Binary Clusters AuPt and Au6Pt: Structure and Reactivity within Density Functional Theory

Within density functional theory with the general gradient approximation for the exchange and correlation, the bimetallic clusters AuPt and Au(6)Pt have been studied for their structure and reactivity. The bond strength of AuPt lies between those of Au(2) and Pt(2), and it is closer to that of Au(2). The Pt atom is the reactive center in both AuPt and AuPt(+) according to electronic structure analysis. AuPt(+) is more stable than AuPt. Au(6)Pt prefers electronic states with low multiplicity. The most stable conformation of Au(6)Pt is a singlet and has quasi-planar hexagonal frame with Pt lying at the hexagonal center. The doping of Pt in Au cluster enhances the chemical regioselectivity of the Au cluster. The Pt atom essentially serves as electron donor and the Au atoms bonded to the Pt atom acts as electron acceptor in Au(6)Pt. The lowest triplet of edge-capped rhombus Au(6)Pt clusters is readily accessible with very small singlet-triplet energy gap (0.32 eV). O(2) prefers to adsorb on Au and CO prefers to adsorb on Pt. O(2) and CO have stronger adsorption on AuPt than they do on Au(6)Pt. CO has a much stronger adsorption on AuPt bimetallic cluster than O(2) does. The adsorption of CO on Pt modifies the geometry of AuPt bimetallic clusters.

An atlas of forecasted molecular data. 2. Vibration frequencies of main-group and transition-metal neutral gas-phase diatomic molecules in the ground state

This atlas of diatomic-molecular vibration frequencies parallels the previously offered Atlas of Internuclear Separations. The Atlas was produced by mining the data from Huber and Herzberg and training neural network software to forecast new data. New protocols were employed with the powerful software, which was originally designed for forecasting the financial markets. The Atlas presents 1920 additional vibration frequencies for use until critical tables are available to fill the needs more precisely. The precision of the predictions is characterized by the average fractional 1% confidence limit, that is, 10.66%. The accuracies of the predictions are determined in two ways. First, 221 of the 224 Huber and Herzberg data values used for training and validation fall within the prediction confidence limits or fall outside by less than 10% of the Huber and Herzberg values, and 181 values agree (within the limits). Second, 87 of 101 comparison data values, consisting of literature data and some additional Huber and Herzberg values, fall within the prediction confidence limits or fall outside by less than half the prediction values, and 44 of the 101 values agree (within the limits).

Computational strategies for a four-component Dirac–Kohn–Sham program: Implementation and first applications

An implementation of the generalized gradient approximation within the four-component formulation of relativistic density-functional theory using G-spinor basis sets is presented. This approach is based on the direct evaluation of the relativistic density and its gradient from the G-spinor amplitudes and gradients without explicit reference to the total density matrix. This proves to be a particularly efficient scheme, with an intrinsic computational cost that scales linearly with the number of G-spinor basis functions. In order to validate this new implementation, incorporated in the parallel version of the program BERTHA, a detailed study of the diatomic system CsAu is also reported. The spectroscopic constants D(e),r(e),omega(e), and x(e)omega(e) and the dipole moment mu have been calculated and compared with the best available theoretical and experimental data. The sensitivity of our results to the details of the numerical schemes used to evaluate the matrix elements is analyzed in detail. Also presented is a comparative study of molecular properties in the alkali auride series which have been obtained using several standard non-relativistic density functionals.

Analysis of the Effect of Spin−Orbit Coupling on the Electronic Structure and Excitation Spectrum of the Bi22- Anion in (K-crypt)2Bi2 on the Basis of Relativistic Electronic Structure Calculations

The Bi2(2-) anions that have been characterized in (K-crypt)2Bi2 are isoelectronic with O2 but are diamagnetic and EPR-silent, unlike O2. The UV-vis spectrum measured for (K-crypt)2Bi2 shows two broad absorption peaks located at 2.05 and 2.85 eV, but no absorption at lower energies down to 0.62 eV. To account for these observations, the electronic structures of the isoelectronic diatomic dianions Q2(2-) (Q = N, P, As, Sb, Bi) were compared on the basis of relativistic density functional theory calculations, and the electronic excitations of Bi2(2-) were analyzed on the basis of relativistic configuration interaction calculations. The extent of spin-orbit coupling, brought about by the relativistic effect, increases steadily in the order N < P < As < Sb < Bi such that the "closed-shell" state is more stable than the "open-shell" state for Bi2(2-), while the opposite is the case for N2(2-), P2(2-), As2(2-), and Sb2(2-). The nature of the electronic excitations of Bi2(2-) was assigned and discussed from the viewpoint of molecular orbitals in the absence of spin-orbit coupling.

Reactions of Gold Atoms and Small Clusters with CO: Infrared Spectroscopic and Theoretical Characterization of AunCO (n = 1−5) and Aun(CO)2 (n = 1, 2) in Solid Argon

Laser-ablated Au atoms have been co-deposited with CO molecules in solid argon to produce gold carbonyls. In addition to the previously reported Au(CO)n (n = 1, 2) and Au2(CO)2 molecules, small gold cluster monocarbonyls Au(n)CO (n = 2-5) are formed on sample annealing and characterized using infrared spectroscopy on the basis of the results of the isotopic substitution and CO concentration change and comparison with theoretical predictions. Of particular interest is that the mononuclear gold carbonyls, Au(CO)n (n = 1, 2), are favored under the experimental conditions of higher CO concentration and lower laser energy, whereas the yields of the gold cluster carbonyls, Au(n)CO (n = 2-5) and Au2(CO)2, remarkably increase with lower CO concentration and higher laser power. Density functional theory (DFT) calculations have been performed on these molecules and the corresponding small naked gold clusters. The identities of these gold carbonyls Au(n)CO (n = 1-5) and Au(n)(CO)2 (n = 1, 2) are confirmed by the good agreement between the experimental and calculated vibrational frequencies, relative absorption intensities, and isotopic shifts.

Time-dependent four-component relativistic density functional theory for excitation energies

Time-dependent four-component relativistic density functional theory within the linear response regime is developed for calculating excitation energies of heavy element containing systems. Since spin is no longer a good quantum number in this context, we resort to time-reversal adapted Kramers basis when deriving the coupled Dirac-Kohn-Sham equation. The particular implementation of the formalism into the Beijing density functional program package utilizes the multipolar expansion of the induced density to facilitate the construction of the induced Coulomb potential. As the first application, pilot calculations on the valence excitation energies and fine structures of the rare gas (Ne to Rn) and Group 12 (Zn to Hg) atoms are reported. To the best of our knowledge, it is the first time to be able to account for spin-orbit coupling within time-dependent density functional theory for excitation energies.

Numerical examination of performance of some exchange-correlation functionals for molecules containing heavy elements

The performance of 17 exchange-correlation functionals for molecules containing heavy elements are numerically examined through four-component relativistic density DFT calculations. The examined functionals show the similar accuracy as they do for the molecules containing light elements only except for bond lengths. LDA and OP86 produce good results for bond lengths and frequencies but bad bond energies. Different functionals do not show much different performance for bond energies except LDA. BP86 and GP86 produce results with average accuracy while LYP does not perform well. Although encouraging results are obtained with functional B97GGA-1, other heavily parameterized and meta-GGA functionals do not produce impressive results.

Relativistic density functional calculations for Pt2

First full-relativistic density functional calculations with the extension of the spin-polarization functional for the relativistic density functional theory in their collinear and noncollinear form are presented here for the molecular system Pt2. The agreement with experiment is very good.

Analytical energy gradient evaluation in relativistic and nonrelativistic density functional calculations

The expressions of analytical energy gradients in density functional theory and their implementation in programs are reported. The evaluation of analytical energy gradients can be carried out in the fully 4-component relativistic, approximate relativistic, and nonrelativistic density functional calculations under local density approximation or general gradient approximation with or without frozen core approximation using different basis sets in our programs. The translational invariance condition and the fact that the one-center terms do not contribute to the energy gradients are utilized to improve the calculation accuracy and to reduce the computational effort. The calculated results of energy gradients and optimized geometry as well as atomization energies of some molecules by the analytical gradient method are in very good agreement with results obtained by the numerical derivative method.

Four-component relativistic Kohn-Sham theory

A four-component relativistic implementation of Kohn-Sham theory for molecular systems is presented. The implementation is based on a nonredundant exponential parametrization of the Kohn-Sham energy, well suited to studies of molecular static and dynamic properties as well as of total electronic energies. Calculations are presented of the bond lengths and the harmonic and anharmonic vibrational frequencies of Au(2), Hg(2+)(2), HgAu(+), HgPt, and AuH. All calculations are based on the full four-component Dirac-Coulomb Hamiltonian, employing nonrelativistic local, gradient-corrected, and hybrid density functionals. The relevance of the Coulomb and Breit operators for the construction of relativistic functionals is discussed; it is argued that, at the relativistic level of density-functional theory and in the absence of a vector potential, the neglect of current functionals follows from the neglect of the Breit operator.

Relativistic electronic structure theory

The theoretical and technical foundations are presented for the efficient relativistic electronic structure theories to treat heavy-atomic molecular systems. This review contains two surveys of four-component and two-component quasi-relativistic approaches. First, we review our highly efficient computational scheme for four-component relativistic ab initio molecular orbital (MO) methods over generally contracted spherical harmonic Gaussian-type spinors (GTSs). Illustrative calculations, which are performed with a new four-component relativistic ab initio molecular orbital program package REL4D, clearly show the efficiency of our computational scheme by the Dirac-Hartree-Fock (DHF) and Dirac-Hartree-Fock (DKS) methods. Next, in the two-component quasi-relativistic framework, two relativistic Hamiltonians, RESC and higher order Douglas-Kroll (DK) Hamiltonians, are introduced, and several illustrative calculations are shown. Numerical results for several systems show that good accuracy can be obtained with our third-order DK (DK3) Hamiltonian.

Spin densities in two-component relativistic density functional calculations: noncollinear versus collinear approach

With present day exchange-correlation functionals, accurate results in nonrelativistic open shell density functional calculations can only be obtained if one uses functionals that do not only depend on the electron density but also on the spin density. We consider the common case where such functionals are applied in relativistic density functional calculations. In scalar-relativistic calculations, the spin density can be defined conventionally, but if spin-orbit coupling is taken into account, spin is no longer a good quantum number and it is not clear what the "spin density" is. In many applications, a fixed quantization axis is used to define the spin density ("collinear approach"), but one can also use the length of the local spin magnetization vector without any reference to an external axis ("noncollinear approach"). These two possibilities are compared in this work both by formal analysis and numerical experiments. It is shown that the (nonrelativistic) exchange-correlation functional should be invariant with respect to rotations in spin space, and this only holds for the noncollinear approach. Total energies of open shell species are higher in the collinear approach because less exchange energy is assigned to a given Kohn-Sham reference function. More importantly, the collinear approach breaks rotational symmetry, that is, in molecular calculations one may find different energies for different orientations of the molecule. Data for the first ionization potentials of Tl, Pb, element 113, and element 114, and for the orientation dependence of the total energy of I+2 and PbF indicate that the error introduced by the collinear approximation is approximately 0.1 eV for valence ionization potentials, but can be much larger if highly ionized open shell states are considered. Rotational invariance is broken by the same amount. This clearly indicates that the collinear approach should not be used, as the full treatment is easily implemented and does not introduce much more computational effort.

The Dirac equation in quantum chemistry: strategies to overcome the current computational problems

A perspective on the use of the relativistic Dirac equation in quantum chemistry is given. It is demonstrated that many of the computational problems that plague the current implementations of the different electronic structure methods can be overcome by utilizing the locality of the small component wave function and density. Possible applications of such new and more efficient formulations are discussed.