-
MICROSIM: A high performance phase-field solver based on CPU and GPU implementations
Authors:
Tanmay Dutta,
Dasari Mohan,
Saurav Shenoy,
Nasir Attar,
Abhikshek Kalokhe,
Ajay Sagar,
Swapnil Bhure,
Swaroop . S. Pradhan,
Jitendriya Praharaj,
Subham Mridha,
Anshika Kushwaha,
Vaishali Shah,
M. P. Gururajan,
V. Venkatesh Shenoi,
Gandham Phanikumar,
Saswata Bhattacharyya,
Abhik Choudhury
Abstract:
The phase-field method has become a useful tool for the simulation of classical metallurgical phase transformations as well as other phenomena related to materials science. The thermodynamic consistency that forms the basis of these formulations lends to its strong predictive capabilities and utility. However, a strong impediment to the usage of the method for typical applied problems of industria…
▽ More
The phase-field method has become a useful tool for the simulation of classical metallurgical phase transformations as well as other phenomena related to materials science. The thermodynamic consistency that forms the basis of these formulations lends to its strong predictive capabilities and utility. However, a strong impediment to the usage of the method for typical applied problems of industrial and academic relevance is the significant overhead with regard to the code development and know-how required for quantitative model formulations. In this paper, we report the development of an open-source phase-field software stack that contains generic formulations for the simulation of multi-phase and multi-component phase transformations. The solvers incorporate thermodynamic coupling that allows the realization of simulations with real alloys in scenarios directly relevant to the materials industry. Further, the solvers utilize parallelization strategies using either multiple CPUs or GPUs to provide cross-platform portability and usability on available supercomputing machines. Finally, the solver stack also contains a graphical user interface to gradually introduce the usage of the software. The user interface also provides a collection of post-processing tools that allow the estimation of useful metrics related to microstructural evolution.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Fluid dynamics alters liquid-liquid phase separation in confined aqueous two-phase systems
Authors:
Eric W. Hester,
Sean P. Carney,
Vishwesh Shah,
Alyssa Arnheim,
Bena Patel,
Dino Di Carlo,
Andrea L. Bertozzi
Abstract:
Liquid-liquid phase separation is key to understanding aqueous two-phase systems (ATPS) arising throughout cell biology, medical science, and the pharmaceutical industry. Controlling the detailed morphology of phase-separating compound droplets leads to new technologies for efficient single-cell analysis, targeted drug delivery, and effective cell scaffolds for wound healing. We present a computat…
▽ More
Liquid-liquid phase separation is key to understanding aqueous two-phase systems (ATPS) arising throughout cell biology, medical science, and the pharmaceutical industry. Controlling the detailed morphology of phase-separating compound droplets leads to new technologies for efficient single-cell analysis, targeted drug delivery, and effective cell scaffolds for wound healing. We present a computational model of liquid-liquid phase separation relevant to recent laboratory experiments with gelatin-polyethylene glycol mixtures. We include buoyancy and surface-tension-driven finite viscosity fluid dynamics with thermally induced phase separation. We show that the fluid dynamics greatly alters the evolution and equilibria of the phase separation problem. Notably, buoyancy plays a critical role in driving the ATPS to energy-minimizing crescent-shaped morphologies and shear flows can generate a tenfold speedup in particle formation. Neglecting fluid dynamics produces incorrect minimum-energy droplet shapes. The model allows for optimization of current manufacturing procedures for structured microparticles and improves understanding of ATPS evolution in confined and flowing settings important in biology and biotechnology.
△ Less
Submitted 26 October, 2023;
originally announced October 2023.
-
Electron-phonon interaction contribution to the total energy of group IV semiconductor polymorphs: evaluation and implications
Authors:
R. Arjun Varma,
Shilpa Paul,
Anup Itale,
Pranav Pable,
Radhika Tibrewala,
Samruddhi Dodal,
Harshal Yerunkar,
Saurav Bhaumik,
Vaishali Shah,
M. P. Gururajan,
T. R. S. Prasanna
Abstract:
In density functional theory (DFT) based total energy studies, the van der Waals (vdW) and zero-point vibrational energy (ZPVE) correction terms are included to obtain energy differences between polymorphs. We propose and compute a new correction term to the total energy, due to electron-phonon interactions (EPI). We rely on Allen's general formalism, which goes beyond the Quasi-Harmonic Approxima…
▽ More
In density functional theory (DFT) based total energy studies, the van der Waals (vdW) and zero-point vibrational energy (ZPVE) correction terms are included to obtain energy differences between polymorphs. We propose and compute a new correction term to the total energy, due to electron-phonon interactions (EPI). We rely on Allen's general formalism, which goes beyond the Quasi-Harmonic Approximation (QHA), to include the free energy contributions due to quasiparticle interactions. We show that, for semiconductors and insulators, the EPI contributions to the free energies of electrons and phonons are the corresponding zero-point energy contributions. Using an approximate version of Allen's formalism in combination with the Allen-Heine theory for EPI corrections, we calculate the zero-point EPI corrections to the total energy for cubic and hexagonal polytypes of Carbon, Silicon and Silicon Carbide. The EPI corrections alter the energy differences between polytypes. In SiC polytypes, the EPI correction term is more sensitive to crystal structure than the vdW and ZPVE terms and is thus essential in determining their energy differences. It clearly establishes that the cubic SiC-3C is metastable and hexagonal SiC-4H is the stable polytype. Our results are consistent with the experimental results of Kleykamp. Our study enables the inclusion of EPI corrections as a separate term in the free energy expression. This opens the way to go beyond the QHA by including the contribution of EPI on all thermodynamic properties.
△ Less
Submitted 16 March, 2023; v1 submitted 14 April, 2022;
originally announced April 2022.
-
AutoMat: Accelerated Computational Electrochemical systems Discovery
Authors:
Emil Annevelink,
Rachel Kurchin,
Eric Muckley,
Lance Kavalsky,
Vinay I. Hegde,
Valentin Sulzer,
Shang Zhu,
Jiankun Pu,
David Farina,
Matthew Johnson,
Dhairya Gandhi,
Adarsh Dave,
Hongyi Lin,
Alan Edelman,
Bharath Ramsundar,
James Saal,
Christopher Rackauckas,
Viral Shah,
Bryce Meredig,
Venkatasubramanian Viswanathan
Abstract:
Large-scale electrification is vital to addressing the climate crisis, but several scientific and technological challenges remain to fully electrify both the chemical industry and transportation. In both of these areas, new electrochemical materials will be critical, but their development currently relies heavily on human-time-intensive experimental trial and error and computationally expensive fi…
▽ More
Large-scale electrification is vital to addressing the climate crisis, but several scientific and technological challenges remain to fully electrify both the chemical industry and transportation. In both of these areas, new electrochemical materials will be critical, but their development currently relies heavily on human-time-intensive experimental trial and error and computationally expensive first-principles, meso-scale and continuum simulations. We present an automated workflow, AutoMat, that accelerates these computational steps by introducing both automated input generation and management of simulations across scales from first principles to continuum device modeling. Furthermore, we show how to seamlessly integrate multi-fidelity predictions such as machine learning surrogates or automated robotic experiments "in-the-loop". The automated framework is implemented with design space search techniques to dramatically accelerate the overall materials discovery pipeline by implicitly learning design features that optimize device performance across several metrics. We discuss the benefits of AutoMat using examples in electrocatalysis and energy storage and highlight lessons learned.
△ Less
Submitted 13 May, 2022; v1 submitted 3 November, 2020;
originally announced November 2020.
-
Physics-aware Deep Generative Models for Creating Synthetic Microstructures
Authors:
Rahul Singh,
Viraj Shah,
Balaji Pokuri,
Soumik Sarkar,
Baskar Ganapathysubramanian,
Chinmay Hegde
Abstract:
A key problem in computational material science deals with understanding the effect of material distribution (i.e., microstructure) on material performance. The challenge is to synthesize microstructures, given a finite number of microstructure images, and/or some physical invariances that the microstructure exhibits. Conventional approaches are based on stochastic optimization and are computation…
▽ More
A key problem in computational material science deals with understanding the effect of material distribution (i.e., microstructure) on material performance. The challenge is to synthesize microstructures, given a finite number of microstructure images, and/or some physical invariances that the microstructure exhibits. Conventional approaches are based on stochastic optimization and are computationally intensive. We introduce three generative models for the fast synthesis of binary microstructure images. The first model is a WGAN model that uses a finite number of training images to synthesize new microstructures that weakly satisfy the physical invariances respected by the original data. The second model explicitly enforces known physical invariances by replacing the traditional discriminator in a GAN with an invariance checker. Our third model combines the first two models to reconstruct microstructures that respect both explicit physics invariances as well as implicit constraints learned from the image data. We illustrate these models by reconstructing two-phase microstructures that exhibit coarsening behavior. The trained models also exhibit interesting latent variable interpolation behavior, and the results indicate considerable promise for enforcing user-defined physics constraints during microstructure synthesis.
△ Less
Submitted 21 November, 2018;
originally announced November 2018.
-
Finite temperature electronic structure of Diamond and Silicon
Authors:
Vaishali Shah,
Bhavik Sanghavi,
Rahul Ramchandani,
M. P. Gururajan,
T. R. S. Prasanna
Abstract:
The electron-phonon interaction contribution to the electronic energies is included in density functional total energy calculations with ab initio pseudopotentials via the Allen formalism [Phys. Rev. B 18, 5217 (1978)] to obtain temperature dependent electronic structure of diamond and silicon. This method allows us to obtain the thermally-averaged ab initio electronic structure in a straightforwa…
▽ More
The electron-phonon interaction contribution to the electronic energies is included in density functional total energy calculations with ab initio pseudopotentials via the Allen formalism [Phys. Rev. B 18, 5217 (1978)] to obtain temperature dependent electronic structure of diamond and silicon. This method allows us to obtain the thermally-averaged ab initio electronic structure in a straightforward and computationally inexpensive way. Our investigations on the finite temperature electronic structure of diamond and silicon lead to a new criterion, temperature transferability, which is required in the ab initio pseudopotentials for temperature dependent studies. The temperature transferability of the Troullier-Martins pseudopotentials used in this work is strongly dependent on the cut-off radius and the inclusion of the unbound 3d$^0$ state. The finite temperature indirect band gaps are highly sensitive to the choice of cut-off radius used in the pseudopotentials. The finite temperature band structures and density of states show that thermal vibrations affect the electron energies throughout the valence and conduction band. We compare our results on the band gap shifts with that due to the Debye-Waller term in the Allen-Heine theory and discuss the observed differences in the zero point and high temperature band gap shifts. Although, the electron energy shifts in the highest occupied valence band and lowest unoccupied conduction band enable to obtain the changes in the indirect and direct band gaps at finite temperatures, the shifts in other electronic levels with temperature enable investigations into the finite temperature valence charge distribution in the bonding region. Thus, we demonstrate that the Allen theory provides a simple and theoretically justified formalism to obtain finite temperature valence electron charge densities that go beyond the rigid pseudo-atom approximation.
△ Less
Submitted 16 February, 2018;
originally announced February 2018.
-
Magnetic properties of 2D nickel nanostrips: structure dependent magnetic anomaly
Authors:
Vikas Kashid,
Vaishali Shah,
H. G. Salunke,
Yuriy Mokrousov,
Stefan Blügel
Abstract:
We have investigated different geometries of two dimensional (2D) infinite length Ni nanowires of increasing width using spin density functional theory calculations. Our simulations demonstrate that the parallelogram motif is the most stable and structures that incorporate the parallelogram motif are more stable as compared to rectangular structures. The wires are conducting and the conductance ch…
▽ More
We have investigated different geometries of two dimensional (2D) infinite length Ni nanowires of increasing width using spin density functional theory calculations. Our simulations demonstrate that the parallelogram motif is the most stable and structures that incorporate the parallelogram motif are more stable as compared to rectangular structures. The wires are conducting and the conductance channels increase with increasing width. The wires have a non-linear behavior in the ballistic anistropic magnetoresistance ratios with respect to the magnetization directions. All 2D nanowires as well as Ni (111) and Ni (100) monolayer investigated are ferromagnetic under the Stoner criterion and exhibit enhanced magnetic moments as compared to bulk Ni and the respective Ni monolayers. The Stoner parameter is seen to depend on the structure and the dimension of the Nws. The easy axis for all nickel nanowires under investigation is observed to be along the wire axis. The double rectangular nanowire exhibits a magnetic anomaly with a smaller magnetic moment when compared to Ni (100) monolayer and is the only structure with an easy axis perpendicular to the wire axis.
△ Less
Submitted 11 February, 2015;
originally announced February 2015.
-
Dzyaloshinskii-Moriya interaction and chiral magnetism in 3$d$-5$d$ zig-zag chains: Tight-binding model and ab initio calculations
Authors:
Vikas Kashid,
Timo Schena,
Bernd Zimmermann,
Yuriy Mokrousov,
Stefan Blügel,
Vaishali Shah,
H. G. Salunke
Abstract:
We investigate the chiral magnetic order in free-standing planar 3$d$-5$d$ bi-atomic metallic chains (3$d$: Fe, Co; 5$d$: Ir, Pt, Au) using first-principles calculations based on density functional theory. We found that the antisymmetric exchange interaction, commonly known as Dzyaloshinskii-Moriya interaction (DMI), contributes significantly to the energetics of the magnetic structure. We used th…
▽ More
We investigate the chiral magnetic order in free-standing planar 3$d$-5$d$ bi-atomic metallic chains (3$d$: Fe, Co; 5$d$: Ir, Pt, Au) using first-principles calculations based on density functional theory. We found that the antisymmetric exchange interaction, commonly known as Dzyaloshinskii-Moriya interaction (DMI), contributes significantly to the energetics of the magnetic structure. We used the full-potential linearized augmented plane wave method and performed self-consistent calculations of homogeneous spin spirals, calculating the DMI by treating the effect of spin-orbit interaction (SOI) in the basis of the spin-spiral states in first-order perturbation theory. To gain insight into the DMI results of our ab initio calculations, we develop a minimal tight-binding model of three atoms and 4 orbitals that contains all essential features: the spin-canting between the magnetic $3d$ atoms, the spin-orbit interaction at the $5d$ atoms, and the structure inversion asymmetry facilitated by the triangular geometry. We found that spin-canting can lead to spin-orbit active eigenstates that split in energy due to the spin-orbit interaction at the $5d$ atom. We show that, the sign and strength of the hybridization, the bonding or antibonding character between $d$-orbitals of the magnetic and non-magnetic sites, the bandwidth and the energy difference between states occupied and unoccupied states of different spin projection determine the sign and strength of the DMI. The key features observed in the trimer model are also found in the first-principles results.
△ Less
Submitted 2 June, 2014;
originally announced June 2014.
-
Superconducting platinum silicide for electron cooling in silicon
Authors:
M J Prest,
J S Richardson-Bullock,
Q T Zhao,
J T Muhonen,
D Gunnarsson,
M Prunnila,
V A Shah,
T E Whall,
E H C Parker,
D R Leadley
Abstract:
We demonstrate electron cooling in silicon using platinum silicide as a superconductor contact to selectively remove the highest energy electrons. The superconducting critical temperature of bulk PtSi is reduced from around 1 K to 0.79 K using a thin film (10 nm) of PtSi, which enhances cooling performance at lower temperatures and enables electron cooling to be demonstrated from 100 mK to 50 mK.
We demonstrate electron cooling in silicon using platinum silicide as a superconductor contact to selectively remove the highest energy electrons. The superconducting critical temperature of bulk PtSi is reduced from around 1 K to 0.79 K using a thin film (10 nm) of PtSi, which enhances cooling performance at lower temperatures and enables electron cooling to be demonstrated from 100 mK to 50 mK.
△ Less
Submitted 16 May, 2014;
originally announced May 2014.
-
Structural and magnetic properties of Pr-alloyed MnBi nanostructures
Authors:
P. Kharel,
V. R. Shah,
X. Z. Li,
W. Y. Zhang,
R. Skomski,
J. E. Shield,
D. J. Sellmyer
Abstract:
The structural and magnetic properties of Pr-alloyed MnBi (short MnBi-Pr) nanostructures with a range of Pr concentrations have been investigated. The nanostructures include thin films having Pr concentrations 0, 2, 3, 5 and 9 atomic percent and melt-spun ribbons having Pr concentrations 0, 2, 4 and 6 percent respectively. Addition of Pr into the MnBi lattice has produced a significant change in t…
▽ More
The structural and magnetic properties of Pr-alloyed MnBi (short MnBi-Pr) nanostructures with a range of Pr concentrations have been investigated. The nanostructures include thin films having Pr concentrations 0, 2, 3, 5 and 9 atomic percent and melt-spun ribbons having Pr concentrations 0, 2, 4 and 6 percent respectively. Addition of Pr into the MnBi lattice has produced a significant change in the magnetic properties of these nanostructures including an increase in coercivity and structural phase transition temperature, and a decrease in saturation magnetization and anisotropy energy. The highest value of coercivity measured in the films is 23 kOe and in the ribbons is 5.6 kOe. The observed magnetic properties are explained as the consequences of competing ferromagnetic and antiferromagnetic interactions.
△ Less
Submitted 28 January, 2013;
originally announced January 2013.
-
Strain enhanced electron cooling in a degenerately doped semiconductor
Authors:
M. J. Prest,
J. T. Muhonen,
M. Prunnila,
D. Gunnarsson,
V. A. Shah,
J. S. Richardson-Bullock,
A. Dobbie,
M. Myronov,
R. J. H. Morris,
T. E. Whall,
E. H. C. Parker,
D. R. Leadley
Abstract:
Enhanced electron cooling is demonstrated in a strained-silicon/superconductor tunnel junction refrigerator of volume 40 um^3. The electron temperature is reduced from 300 mK to 174 mK, with the enhancement over an unstrained silicon control (300 mK to 258 mK) being attributed to the smaller electron-phonon coupling in the strained case. Modeling and the resulting predictions of silicon-based cool…
▽ More
Enhanced electron cooling is demonstrated in a strained-silicon/superconductor tunnel junction refrigerator of volume 40 um^3. The electron temperature is reduced from 300 mK to 174 mK, with the enhancement over an unstrained silicon control (300 mK to 258 mK) being attributed to the smaller electron-phonon coupling in the strained case. Modeling and the resulting predictions of silicon-based cooler performance are presented. Further reductions in the minimum temperature are expected if the junction sub-gap leakage and tunnel resistance can be reduced. However, if only tunnel resistance is reduced, Joule heating is predicted to dominate.
△ Less
Submitted 18 November, 2011; v1 submitted 2 November, 2011;
originally announced November 2011.
-
Electronic structure effects in stability and quantum conductance in 2D gold nanowires
Authors:
Vikas Kashid,
Vaishali Shah,
H. G. Salunke
Abstract:
We have investigated the stability and conductivity of unsupported, two dimensional infinite gold nanowires using ab-initio density functional theory (DFT). Two dimensional ribbon like nanowires, with 1-5 rows of gold atoms in the non-periodic direction and with different possible structures have been considered. The nanowires with > 2 rows of atoms exhibit dimerization, similar to finite wires, a…
▽ More
We have investigated the stability and conductivity of unsupported, two dimensional infinite gold nanowires using ab-initio density functional theory (DFT). Two dimensional ribbon like nanowires, with 1-5 rows of gold atoms in the non-periodic direction and with different possible structures have been considered. The nanowires with > 2 rows of atoms exhibit dimerization, similar to finite wires, along the non-periodic direction. Our results show that in these zero thickness nanowires, the parallelogram motif is the most stable. A comparison between parallelogram and rectangular shaped nanowires of increasing width indicates that zero thickness (111) oriented wires have a higher stability over (100). A detailed analysis of the electronic structure, reveals that the (111) oriented structures show increased delocalization of s and p electrons in addition to a stronger delocalization of the d electrons and hence are the most stable. The density of states show that the nanowires are metallic and conducting except for the double zigzag structure, which, is semiconducting. Conductance calculations show transmission for a wide range of energies in all the stable nanowires with more than two rows of atoms. The conductance channels are not purely s and have strong contributions from the d levels and weak contributions from the p levels.
△ Less
Submitted 28 March, 2011;
originally announced March 2011.
-
Strain control of electron-phonon energy loss rate in many-valley semiconductors
Authors:
J. T. Muhonen,
M. J. Prest,
M. Prunnila,
D. Gunnarsson,
V. A. Shah,
A. Dobbie,
M. Myronov,
R. J. H. Morris,
T. E. Whall,
E. H. C. Parker,
D. R. Leadley
Abstract:
We demonstrate significant modification of the electron-phonon energy loss rate in a many-valley semiconductor system due to lattice mismatch induced strain. We show that the thermal conductance from the electron system to the phonon bath in strained n + Si, at phonon temperatures between 200 mK and 450 mK, is more than an order of magnitude lower than that for a similar unstrained sample.
We demonstrate significant modification of the electron-phonon energy loss rate in a many-valley semiconductor system due to lattice mismatch induced strain. We show that the thermal conductance from the electron system to the phonon bath in strained n + Si, at phonon temperatures between 200 mK and 450 mK, is more than an order of magnitude lower than that for a similar unstrained sample.
△ Less
Submitted 31 January, 2011;
originally announced January 2011.
-
Anisotropic weakly localized transport in nitrogen-doped ultrananocrystalline diamond films
Authors:
Kunjal V Shah,
Dmitry Churochkin,
Zivayi Chiguvare,
Somnath Bhattacharyya
Abstract:
We establish the dominant effect of anisotropic weak localization (WL) in three dimensions associated with a propagative Fermi surface, on the conductivity correction in heavily nitrogen doped ultrananocrystalline diamond (UNCD) films based on magneto-resistance studies at low temperatures. Also, low temperature electrical conductivity can show weakly localized transport in 3D combined with the ef…
▽ More
We establish the dominant effect of anisotropic weak localization (WL) in three dimensions associated with a propagative Fermi surface, on the conductivity correction in heavily nitrogen doped ultrananocrystalline diamond (UNCD) films based on magneto-resistance studies at low temperatures. Also, low temperature electrical conductivity can show weakly localized transport in 3D combined with the effect of electron-electron interactions in these materials, which is remarkably different from the conductivity in 2DWL or strong localization regime. The corresponding dephasing time of electronic wavefunctions in these systems described as ~ T^-p with p < 1, follows a relatively weak temperature dependence compared to the generally expected nature for bulk dirty metals having $p \geq 1$. The temperature dependence of Hall (electron) mobility together with an enhanced electron density has been used to interpret the unusual magneto-transport features and show delocalized electronic transport in these n-type UNCD films, which can be described as low-dimensional superlattice structures.
△ Less
Submitted 24 November, 2010; v1 submitted 22 July, 2010;
originally announced July 2010.
-
Ab initio molecular dynamics using density based energy functionals: application to ground state geometries of some small clusters
Authors:
Dinesh Nehete,
Vaishali Shah,
D. G. Kanhere
Abstract:
The ground state geometries of some small clusters have been obtained via ab initio molecular dynamical simulations by employing density based energy functionals. The approximate kinetic energy functionals that have been employed are the standard Thomas-Fermi $(T_{TF})$ along with the Weizsacker correction $T_W$ and a combination $F(N_e)T_{TF} + T_W$. It is shown that the functional involving…
▽ More
The ground state geometries of some small clusters have been obtained via ab initio molecular dynamical simulations by employing density based energy functionals. The approximate kinetic energy functionals that have been employed are the standard Thomas-Fermi $(T_{TF})$ along with the Weizsacker correction $T_W$ and a combination $F(N_e)T_{TF} + T_W$. It is shown that the functional involving $F(N_e)$ gives superior charge densities and bondlengths over the standard functional. Apart from dimers and trimers of Na, Mg, Al, Li, Si, equilibrium geometries for $Li_nAl, n=1,8$ and $Al_{13}$ clusters have also been reported. For all the clusters investigated, the method yields the ground state geometries with the correct symmetries with bondlengths within 5\% when compared with the corresponding results obtained via full orbital based Kohn-Sham method. The method is fast and a promising one to study the ground state geometries of large clusters.
△ Less
Submitted 18 October, 1995; v1 submitted 14 October, 1995;
originally announced October 1995.
-
Ab initio molecular dynamics via density based energy functionals
Authors:
Vaishali Shah,
Dinesh Nehete,
D. G. Kanhere
Abstract:
The use of energy functionals based on density as the basic variable is advocated for ab initio molecular dynamics. It is demonstrated that the constraint of positivity of density can be incorporated easily by using square root density for minimization of the energy functional. An ad hoc prescription for including nonlocal pseudopotentials for plane wave based calculations is proposed and is sho…
▽ More
The use of energy functionals based on density as the basic variable is advocated for ab initio molecular dynamics. It is demonstrated that the constraint of positivity of density can be incorporated easily by using square root density for minimization of the energy functional. An ad hoc prescription for including nonlocal pseudopotentials for plane wave based calculations is proposed and is shown to yield improved results. Applications are reported for equilibrium geometries of small finite systems, viz. dimers and trimers of simple metal atoms like Na and Mg, which represent a rather stringent test for approximate kinetic energy functionals involved in such calculations.
△ Less
Submitted 8 July, 1994;
originally announced July 1994.