No Arabic abstract
Charge Distributions in Metallic Alloys: a Charge Excess Functional theory approach
The distribution of local charge excesses (DLC) in metallic alloys, previously obtained as a result of the analysis of order N electronic structure calculations, is derived from a variational principle. A phenomenological Charge Excess Functional (CEF) theory is obtained which is determined by three concentration dependent, material specific, parameters that can be obtained from {it ab initio} calculations. The theory requires modest computational efforts and reproduces with an excellent accuracy the DLC and the electrostatic energies of ordered, substitutionally disordered or segregating metallic alloys and, hence, can be considered an efficient approach alternative to conventional electronic structure calculations. The substantial reduction of computing time opens new perspectives for the understanding of metallic systems and their mechanical properties.
The class of the Generalized Coherent Potential Approximations (GCPA) to the Density Functional Theory (DFT) is introduced within the Multiple Scattering Theory formalism for dealing with, ordered or disordered, metallic alloys. All GCPA theories are based on a common ansatz for the kinetic part of the Hohenberg-Kohn functional and each theory of the class is specified by an external model concerning the potential reconstruction. The GCPA density functional consists of marginally coupled local contributions, does not depend on the details of the charge density and can be exactly rewritten as a function of the appropriate charge multipole moments associated with each lattice site. A general procedure based on the integration of the qV laws is described that allows for the explicit construction the same function. The coarse grained nature of the GCPA density functional implies great computational advantages and is connected with the O(N) scalability of GCPA algorithms. Moreover, it is shown that a convenient truncated series expansion of the GCPA functional leads to the Charge Excess Functional (CEF) theory [E. Bruno, L. Zingales and Y. Wang, Phys. Rev. Lett. {bf 91}, 166401 (2003)] which here is offered in a generalized version that includes multipolar interactions. CEF and the GCPA numerical results are compared with status of art LAPW full-potential density functional calculations for 62, bcc- and fcc-based, ordered CuZn alloys, in all the range of concentrations. These extensive tests show that the discrepancies between GCPA and CEF are always within the numerical accuracy of the calculations, both for the site charges and the total energies. Furthermore, GCPA and CEF very carefully reproduce the LAPW site charges and the total energy trends.
Electronic structure calculations performed on very large supercells have shown that the local charge excesses in metallic alloys are related through simple linear relations to the local electrostatic field resulting from distribution of charges in the whole crystal. By including local external fields in the single site Coherent Potential Approximation theory, we develop a novel theoretical scheme in which the local charge excesses for random alloys can be obtained as the responses to local external fields. Our model maintains all the computational advantages of a single site theory but allows for full charge relaxation at the impurity sites. Through applications to CuPd and CuZn alloys, we find that, as a general rule, non linear charge rearrangements occur at the impurity site as a consequence of the complex phenomena related with the electronic screening of the external potential. This nothwithstanding, we observe that linear relations hold between charge excesses and external potentials, in quantitative agreement with the mentioned supercell calculations, and well beyond the limits of linearity for any other site property.
Electronic structure calculations performed on very large supercells have shown that the local charge excesses in metallic alloys are related through simple linear relations to the local electrostatic field resulting from distribution of charges in the whole crystal. By including local external fields in the single site Coherent Potential Approximation theory, we develop a novel theoretical scheme in which the local charge excesses for random alloys can be obtained as the responses to local external fields. Our model maintains all the computational advantages of a single site theory but allows for full charge relaxation at the impurity sites. Through applications to CuPd and CuZn alloys, we find that, as a general rule, non linear charge rearrangements occur at the impurity site as a consequence of the complex phenomena related with the electronic screening of the external potential. This nothwithstanding, we observe that linear relations hold between charge excesses and external potentials, in quantitative agreement with the mentioned supercell calculations, and well beyond the limits of linearity for any other site property.
A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and separable pairing interaction are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses neutron shell closures at N = 28, 50, 82 and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the differential charge radii. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.