Showing 1–15 of 15 results for author: Ratcliff, L E

Revealing the bonding nature and electronic structure of early transition metal dihydrides

Authors: Curran Kalha, Laura E. Ratcliff, Giorgio Colombi, Christoph Schlueter, Bernard Dam, Andrei Gloskovskii, Tien-Lin Lee, Pardeep K. Thakur, Prajna Bhatt, Yujiang Zhu, Jürg Osterwalder, Francesco Offi, Giancarlo Panaccione, Anna Regoutz

Abstract: Hydrogen as a fuel plays a crucial role in driving the transition to net zero greenhouse gas emissions. To realise its potential, obtaining a means of efficient storage is paramount. One solution is using metal hydrides, owing to their good thermodynamical absorption properties and effective hydrogen storage. Although metal hydrides appear simple compared to many other energy materials, understand… ▽ More Hydrogen as a fuel plays a crucial role in driving the transition to net zero greenhouse gas emissions. To realise its potential, obtaining a means of efficient storage is paramount. One solution is using metal hydrides, owing to their good thermodynamical absorption properties and effective hydrogen storage. Although metal hydrides appear simple compared to many other energy materials, understanding the electronic structure and chemical environment of hydrogen within them remains a key challenge. This work presents a new analytical pathway to explore these aspects in technologically relevant systems using Hard X-ray Photoelectron Spectroscopy (HAXPES) on thin films of two prototypical metal dihydrides: YH$_{2-δ}$ and TiH$_{2-δ}$. By taking advantage of the tunability of synchrotron radiation, a non-destructive depth profile of the chemical states is obtained using core level spectra. Combining experimental valence band spectra collected at varying photon energies with theoretical insights from density functional theory (DFT) calculations, a description of the bonding nature and the role of d versus sp contributions to states near the Fermi energy are provided. Moreover, a reliable determination of the enthalpy of formation is proposed by using experimental values of the energy position of metal s band features close to the Fermi energy in the HAXPES valence band spectra. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Roadmap on Electronic Structure Codes in the Exascale Era

Authors: Vikram Gavini, Stefano Baroni, Volker Blum, David R. Bowler, Alexander Buccheri, James R. Chelikowsky, Sambit Das, William Dawson, Pietro Delugas, Mehmet Dogan, Claudia Draxl, Giulia Galli, Luigi Genovese, Paolo Giannozzi, Matteo Giantomassi, Xavier Gonze, Marco Govoni, Andris Gulans, François Gygi, John M. Herbert, Sebastian Kokott, Thomas D. Kühne, Kai-Hsin Liou, Tsuyoshi Miyazaki, Phani Motamarri , et al. (16 additional authors not shown)

Abstract: Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources d… ▽ More Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini (vikramg@umich.edu) and Danny Perez (danny_perez@lanl.gov)

Tackling Disorder in $γ$-Ga$_2$O$_3$

Authors: Laura E. Ratcliff, Takayoshi Oshima, Felix Nippert, Benjamin M. Janzen, Elias Kluth, Rüdiger Goldhahn, Martin Feneberg, Piero Mazzolini, Oliver Bierwagen, Charlotte Wouters, Musbah Nofal, Martin Albrecht, Jack E. N. Swallow, Leanne A. H. Jones, Pardeep K. Thakur, Tien-Lin Lee, Curran Kalha, Christoph Schlueter, Tim D. Veal, Joel B. Varley, Markus R. Wagner, Anna Regoutz

Abstract: Ga$_2$O$_3$ and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga$_2$O$_3$ offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. $γ$-Ga$_2$O$_3$ presents a particular challenge across synthesis, characterisation, and theory due to its inherent… ▽ More Ga$_2$O$_3$ and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga$_2$O$_3$ offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. $γ$-Ga$_2$O$_3$ presents a particular challenge across synthesis, characterisation, and theory due to its inherent disorder and resulting complex structure -- electronic structure relationship. Here, density functional theory is used in combination with a machine learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the $γ$-phase. Theoretical results are compared with surface and bulk sensitive soft and hard X-ray photoelectron spectroscopy, X-ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of $γ$-Ga$_2$O$_3$. The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by PLE spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step towards exploring how their electronic structure can be understood in terms of local coordination and overall structure. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Exploring Metastable States in UO$_2$ using Hybrid Functionals and Dynamical Mean Field Theory

Authors: Laura E. Ratcliff, Luigi Genovese, Hyowon Park, Peter B. Littlewood, Alejandro Lopez-Bezanilla

Abstract: A detailed exploration of the $f$-atomic orbital occupancy space for UO$_2$ is performed using a first principles approach based on density functional theory (DFT), employing a full hybrid functional within a systematic basis set. Specifically, the PBE0 functional is combined with an occupancy biasing scheme implemented in a wavelet-based algorithm which is adapted to large supercells. The results… ▽ More A detailed exploration of the $f$-atomic orbital occupancy space for UO$_2$ is performed using a first principles approach based on density functional theory (DFT), employing a full hybrid functional within a systematic basis set. Specifically, the PBE0 functional is combined with an occupancy biasing scheme implemented in a wavelet-based algorithm which is adapted to large supercells. The results are compared with previous DFT+U calculations reported in the literature, while dynamical mean field theory (DMFT) is also performed to provide a further base for comparison. This work shows that the computational complexity of the energy landscape of a correlated $f$-electron oxide is much richer than has previously been demonstrated. The resulting calculations provide evidence of the existence of multiple previously unexplored metastable electronic states of UO$_2$, including those with energies which are lower than previously reported ground states. △ Less

Submitted 30 September, 2021; originally announced September 2021.

Lifetime effects and satellites in the photoelectron spectrum of tungsten metal

Authors: Curran Kalha, Laura E. Ratcliff, Julio J. Gutiérrez Moreno, Stephan Mohr, Mervi Mantsinen, Nathalie K. Fernando, Pardeep K. Thakur, Tien-Lin Lee, Hsiang-Han Tseng, Tim S. Nunney, Juhan M. Kahk, Johannes Lischner, Anna Regoutz

Abstract: Tungsten is an important and versatile transition metal and has a firm place at the heart of many technologies. A popular experimental technique for the characterisation of tungsten and tungsten-based compounds is X-ray photoelectron spectroscopy (XPS), which enables the assessment of chemical states and electronic structure through the collection of core level and valence band spectra. However, i… ▽ More Tungsten is an important and versatile transition metal and has a firm place at the heart of many technologies. A popular experimental technique for the characterisation of tungsten and tungsten-based compounds is X-ray photoelectron spectroscopy (XPS), which enables the assessment of chemical states and electronic structure through the collection of core level and valence band spectra. However, in the case of metallic tungsten, open questions remain regarding the origin, nature, and position of satellite features that are prominent in the photoelectron spectrum. These satellites are a fingerprint of the electronic structure of the material and have not been thoroughly investigated, at times leading to their misinterpretation. The present work combines high-resolution soft and hard X-ray photoelectron spectroscopy (SXPS and HAXPES) with reflection electron energy loss spectroscopy (REELS) and a multi-tiered ab-initio theoretical approach, including density functional theory (DFT) and many-body perturbation theory (G0W0 and GW+C), to disentangle the complex set of experimentally observed satellite features attributed to the generation of plasmons and interband transitions. This combined experiment-theory strategy is able to uncover previously undocumented satellite features, improving our understanding of their direct relationship to tungsten's electronic structure. Furthermore, it lays the groundwork for future studies into tungsten based mixed-metal systems and holds promise for the re-assessment of the photoelectron spectra of other transition and post-transition metals, where similar questions regarding satellite features remain. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Thermal and Oxidation Stability of Ti$_x$W$_{1-x}$ Diffusion Barriers Investigated by Soft and Hard X-ray Photoelectron Spectroscopy

Authors: Curran Kalha, Sebastian Bichelmaier, Nathalie K. Fernando, Judith. V. Berens, Pardeep K. Thakur, Tien-Lin Lee, Julio. J. Gutiérrez-Moreno, Stephan Mohr, Laura E. Ratcliff, Michael Reisinger, Johannes Zechner, Michael Nelhiebel, Anna Regoutz

Abstract: The binary alloy of titanium-tungsten (TiW) is an established diffusion barrier in high-power semiconductor devices, owing to its ability to suppress the diffusion of copper from the metallisation scheme into the surrounding silicon substructure. However, little is known about the response of TiW to high temperature events or its behaviour when exposed to air. Here, a combined soft and hard X-ray… ▽ More The binary alloy of titanium-tungsten (TiW) is an established diffusion barrier in high-power semiconductor devices, owing to its ability to suppress the diffusion of copper from the metallisation scheme into the surrounding silicon substructure. However, little is known about the response of TiW to high temperature events or its behaviour when exposed to air. Here, a combined soft and hard X-ray photoelectron spectroscopy (XPS) characterisation approach is used to study the influence of post-deposition annealing and titanium concentration on the oxidation behaviour of a 300~nm-thick TiW film. The combination of both XPS techniques allows for the assessment of the chemical state and elemental composition across the surface and bulk of the TiW layer. The findings show that in response to high-temperature annealing, titanium segregates out of the mixed metal system and upwardly migrates, accumulating at the TiW/air interface. Titanium shows remarkably rapid diffusion under relatively short annealing timescales and the extent of titanium surface enrichment is increased through longer annealing periods or by increasing the precursor titanium concentration. Surface titanium enrichment enhances the extent of oxidation both at the surface and in the bulk of the alloy due to the strong gettering ability of titanium. Quantification of the soft X-ray photoelectron spectra highlights the formation of three tungsten oxidation environments, attributed to WO$_2$, WO$_3$ and a WO$_3$ oxide coordinated with a titanium environment. This combinatorial characterisation approach provides valuable insights into the thermal and oxidation stability of TiW alloys from two depth perspectives, aiding the development of future device technologies. △ Less

Submitted 19 February, 2021; originally announced February 2021.

A Combined Density Functional Theory and X-ray Photoelectron Spectroscopy Study of the Aromatic Amino Acids

Authors: Anna Regoutz, Marta S. Wolinska, Nathalie K. Fernando, Laura E. Ratcliff

Abstract: Amino acids are essential to all life. However, our understanding of some aspects of their intrinsic structure, molecular chemistry, and electronic structure is still limited. In particular the nature of amino acids in their crystalline form, often essential to biological and medical processes, faces a lack of knowledge both from experimental and theoretical approaches. An important experimental t… ▽ More Amino acids are essential to all life. However, our understanding of some aspects of their intrinsic structure, molecular chemistry, and electronic structure is still limited. In particular the nature of amino acids in their crystalline form, often essential to biological and medical processes, faces a lack of knowledge both from experimental and theoretical approaches. An important experimental technique that has provided a multitude of crucial insights into the chemistry and electronic structure of materials is X-ray photoelectron spectroscopy. Whilst the interpretation of spectra of simple bulk inorganic materials is often routine, interpreting core level spectra of complex molecular systems is complicated to impossible without the help of theory. We have previously demonstrated the ability of density functional theory to calculate binding energies of simple amino acids, using $Δ$SCF implemented in a systematic basis set for both gas phase (multiwavelets) and solid state (plane waves) calculations. In this study, we use the same approach to successfully predict and rationalise the experimental core level spectra of phenylalanine (Phe), tyrosine (Tyr), tryptophan (Trp), and histidine (His) and gain an in-depth understanding of their chemistry and electronic structure within the broader context of more than 20 related molecular systems. The insights gained from this study provide significant information on the nature of the aromatic amino acids and their conjugated side chains. △ Less

Submitted 30 October, 2020; originally announced October 2020.

Pseudo-Fragment Approach for Extended Systems Derived from Linear-Scaling DFT

Authors: Laura E. Ratcliff, Luigi Genovese

Abstract: We present a computational approach which is tailored for reducing the complexity of the description of extended systems at the density functional theory level. We define a recipe for generating a set of localized basis functions which are optimized either for the accurate description of pristine, bulk like Wannier functions, or for the \emph{in situ} treatment of deformations induced by defective… ▽ More We present a computational approach which is tailored for reducing the complexity of the description of extended systems at the density functional theory level. We define a recipe for generating a set of localized basis functions which are optimized either for the accurate description of pristine, bulk like Wannier functions, or for the \emph{in situ} treatment of deformations induced by defective constituents such as boundaries or impurities. Our method enables one to identify the regions of an extended system which require dedicated optimization of the Kohn-Sham degrees of freedom, and provides the user with a reliable estimation of the errors -- if any -- induced by the locality of the approach. Such a method facilitates on the one hand an effective reduction of the computational degrees of freedom needed to simulate systems at the nanoscale, while in turn providing a description that can be straightforwardly put in relation to effective models, like tight binding Hamiltonians. We present our methodology with SiC nanotube-like cages as a test bed. Nonetheless, the wavelet-based method employed in this paper makes possible calculation of systems with different dimensionalities, including slabs and fully periodic systems. △ Less

Submitted 12 April, 2019; v1 submitted 16 January, 2019; originally announced January 2019.

Accurate, Large-Scale and Affordable Hybrid-PBE0 Calculations with GPU-Accelerated Supercomputers

Authors: Laura E. Ratcliff, A. Degomme, José A. Flores-Livas, Stefan Goedecker, Luigi Genovese

Abstract: Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding - when not out of reach - if high quality basis sets are used. We present a highly efficient multiple GPU implementation of the exact exchange operator which allows hybrid functional density-functional theory calculations with systemat… ▽ More Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding - when not out of reach - if high quality basis sets are used. We present a highly efficient multiple GPU implementation of the exact exchange operator which allows hybrid functional density-functional theory calculations with systematic basis sets without additional approximations for up to a thousand atoms. This method is implemented in a portable real-space-based algorithm, released as an open-source package. With such a framework hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution of the same order of magnitude as traditional semilocal-GGA functionals. △ Less

Submitted 21 December, 2017; originally announced December 2017.

Zirconia and hafnia polymorphs -- ground state structural properties from diffusion Monte Carlo

Authors: Hyeondeok Shin, Anouar Benali, Ye Luo, Emily Crabb, Alejandro Lopez-Bezanilla, Laura E. Ratcliff, Andrea M. Jokisaari, Olle Heinonen

Abstract: Zirconia (zirconium dioxide) and hafnia (hafnium dioxide) are binary oxides used in a range of applications. Because zirconium and hafnium are chemically equivalent, they have three similar polymorphs, and it is important to understand the properties and energetics of these polymorphs. However, while density functional theory calculations can get the correct energetic ordering, the energy differen… ▽ More Zirconia (zirconium dioxide) and hafnia (hafnium dioxide) are binary oxides used in a range of applications. Because zirconium and hafnium are chemically equivalent, they have three similar polymorphs, and it is important to understand the properties and energetics of these polymorphs. However, while density functional theory calculations can get the correct energetic ordering, the energy differences between polymorphs depend very much on the specific density functional theory approach, as do other quantities such as lattice constants and bulk modulus. We have used highly accurate quantum Monte Carlo simulations to model the three zirconia and hafnia polymorphs. We compare our results for structural parameters, bulk modulus, and cohesive energy with results obtained from density functional theory calculations. We also discuss comparisons of our results with existing experimental data, in particular for structural parameters where extrapolation to zero temperature can be attempted. We hope our results of structural parameters as well as for cohesive energy and bulk modulus can serve as benchmarks for density-functional theory based calculations and as a guidance for future experiments. △ Less

Submitted 27 June, 2018; v1 submitted 30 August, 2017; originally announced August 2017.

Journal ref: Phys. Rev. Materials 2, 075001 (2018)

Fragment Approach to Constrained Density Functional Theory Calculations using Daubechies Wavelets

Authors: Laura E. Ratcliff, Luigi Genovese, Stephan Mohr, Thierry Deutsch

Abstract: In a recent paper we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions is optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contr… ▽ More In a recent paper we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions is optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e. without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Accurate and efficient linear scaling DFT calculations with universal applicability

Authors: Stephan Mohr, Laura E. Ratcliff, Luigi Genovese, Damien Caliste, Paul Boulanger, Stefan Goedecker, Thierry Deutsch

Abstract: Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens o… ▽ More Density Functional Theory calculations traditionally suffer from an inherent cubic scaling with respect to the size of the system, making big calculations extremely expensive. This cubic scaling can be avoided by the use of so-called linear scaling algorithms, which have been developed during the last few decades. In this way it becomes possible to perform ab-initio calculations for several tens of thousands of atoms or even more within a reasonable time frame. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis -- which offers ideal properties for accurate linear scaling calculations -- we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large systems with only a moderate demand of computing resources. △ Less

Submitted 23 January, 2015; originally announced January 2015.

Daubechies Wavelets for Linear Scaling Density Functional Theory

Authors: Stephan Mohr, Laura E. Ratcliff, Paul Boulanger, Luigi Genovese, Damien Caliste, Thierry Deutsch, Stefan Goedecker

Abstract: We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelet… ▽ More We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems. △ Less

Submitted 3 June, 2014; v1 submitted 29 January, 2014; originally announced January 2014.

Ab initio calculations of the optical absorption spectra of C60-conjugated polymer hybrids

Authors: Laura E. Ratcliff, Peter D. Haynes

Abstract: A recently developed linear-scaling density-functional theory (LS-DFT) formalism is used to calculate optical absorption spectra of hybrids of C60 and the conjugated polymers poly(para-phenylene) (PPP) and poly(para-phenylene vinylene) (PPV). The use of a LS formalism allows calculations on large systems with realistic proportions of C60, which has been of interest for the use of such materials in… ▽ More A recently developed linear-scaling density-functional theory (LS-DFT) formalism is used to calculate optical absorption spectra of hybrids of C60 and the conjugated polymers poly(para-phenylene) (PPP) and poly(para-phenylene vinylene) (PPV). The use of a LS formalism allows calculations on large systems with realistic proportions of C60, which has been of interest for the use of such materials in photovoltaics. Two different bonding structures are tested for the hybrid PPP and for both systems additional peaks are present in the absorption spectra below the original onset of absorption. By identifying the eigenstates involved in the relevant transitions, a weighted density difference is formed, demonstrating the transfer of charge between the polymer chain and the C60, in agreement with experiment. For the hybrid PPV, no additional peaks are observed in the absorption spectrum. △ Less

Submitted 30 March, 2013; v1 submitted 28 March, 2013; originally announced March 2013.

Calculating optical absorption spectra for large systems using linear-scaling density-functional theory

Authors: Laura E. Ratcliff, Nicholas D. M. Hine, Peter D. Haynes

Abstract: A new method for calculating optical absorption spectra within linear-scaling density-functional theory (LS-DFT) is presented, incorporating a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states. Three different schemes are compared and the most promising of these, based on the use of a projection operator, has been implemented in a fully-functiona… ▽ More A new method for calculating optical absorption spectra within linear-scaling density-functional theory (LS-DFT) is presented, incorporating a scheme for optimizing a set of localized orbitals to accurately represent unoccupied Kohn-Sham states. Three different schemes are compared and the most promising of these, based on the use of a projection operator, has been implemented in a fully-functional LS-DFT code. The method has been applied to the calculation of optical absorption spectra for the metal-free phthalocyanine molecule and the conjugated polymer poly(para-phenylene). Excellent agreement with results from a traditional DFT code is obtained. △ Less

Submitted 11 October, 2011; v1 submitted 15 September, 2011; originally announced September 2011.

Comments: Accepted for publication in Physical Review B