Source author record

M. W. C. Dharma-wardana

M. W. C. Dharma-wardana appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Comment on "Water sources and kidney function: investigating chronic kidney disease of unknown etiology in a prospective study", by P. Vlahos et al

Vlahos et al., Ref. 1, NPJ Clean water. 4, 50 (2021) have reported the presence of pesticide contamination above safe levels in a "single time-point analysis" of well water in a region in Sri Lanka where chronic kidney disease of unknown etiology (CKDu) is endemic. They conclude "that agrochemical use in paddy and other agricultural practices ... of the Green Revolution in Sri Lanka may now be contributing to ill health, rapid progression of disease, and mortality". The authors also propose "reducing ... agrochemical contaminants in Sri Lanka and other tropical countries to reduce ... CKDu". These conclusions, based on what they call a "single time-point analysis", tantamount to an identification of the etiology of CKDu are unsupported by the evidence presented by Vlahos et al. They do not satisfy, say, even the simplest of Bradford-Hill criteria for causation. In particular, (i) similar but non-persistent pesticide excesses have been detected sporadically in most parts of the country including where there is no CKDu; (ii) the pesticides reported in (1) cause both hepatotoxicity and nephrotoxicity; the latter with glomerular damage while CKDu is associated with tubulo-interstitial damage where no hepatotoxic symptoms have been reported; (iii) the pesticides detected have short half-lives and are used over short periods during farming; so the one time-point analysis is inadequate and misleading; (iv) farming communities that use pesticides in the same way but remain essentially without CKDu are found to exist adjacent to communities with CKDu; (v) the CKDu prevalence seems to correlate with local geomorphology but without correlation to agriculture which is practiced in most parts of the country. .

preprint2022arXiv

Evidence for multiple Liquid-liquid phase transitions in carbon, and the Friedel-ordering of its liquid state

Carbon, the fourth most abundant element in the Universe forms a metallic fluid with transient covalent bonds on melting. Its liquid-liquid phase transitions, intensely sought using simulations had been elusive. Here we use density functional theory (DFT) simulations with up to 108 atoms using molecular dynamics, as well as one-atom DFT as implemented in the neutral pseudo-atom method where multi-atom effects are treated by ion-ion correlation functionals. Both methods use electron-electron exchange correlation functionals for electron many-body effects. Here we show using both methods, that liquid carbon displays multiple liquid-liquid transitions linked to changes in coordination number in the density range 3 g/cm$^3$, to $\sim$ 6 g/cm$^3$ when a coordination number of 12 is reached. The transitions disappear by 4 eV in temperature. The calculated pressures and transition densities are shown to be sensitive to the exchange-correlation functionals used. Significantly, we find that a simple metallic model yields the structure factors and thermodynamics with quantitative accuracy, without invoking any covalent-bonding features. The ion-ion structure factor for these densities and temperatures is found to have a subpeak tied to twice the Fermi wavevector, constraining the fluid in momentum space. The dominant Friedel oscillations forming the pair interactions correlate the ions and drive the multiple liquid-liquid phase transitions. Our results suggest that liquid carbon typifies a class of fluids whose structure is ordered by the long-range Friedel oscillations in the pair-potentials. These results are critical to terrestrial and astrophysical studies, inertial fusion using carbon drivers, refined shock experiments, and in seeking new carbon-based materials.

preprint2022arXiv

Yukawa Friedel-Tail pair potentials for warm dense matter applications

Accurate equations of state (EOS) and plasma transport properties are essential for numerical simulations of warm dense matter encountered in many high-energy-density situations. Molecular dynamics (MD) is a simulation method that generates EOS and transport data using an externally provided potential to dynamically evolve the particles without further reference to the electrons. To minimize computational cost, pair potentials needed in MD may be obtained from the neutral-pseudoatom (NPA) approach, a form of single-ion density functional theory (DFT), where many-ion effects are included via ion-ion correlation functionals. Standard $N$-ion DFT-MD provides pair potentials via the force matching technique but at much greater computational cost. Here we propose a simple analytic model for pair potentials with physically meaningful parameters based on a Yukawa form with a thermally damped Friedel tail (YFT) applicable to systems containing free electrons. The YFT model accurately fits NPA pair potentials or the non-parametric force-matched potentials from $N$-ion DFT-MD, showing excellent agreement for a wide range of conditions. The YFT form provides accurate extrapolations of the NPA or force-matched potentials for small and large particle separations within a physical model. Our method can be adopted to treat plasma mixtures, allowing for large-scale simulations of multi-species warm dense matter.

preprint2020arXiv

Discussion on a "Dynamic model to conceptualize multiple environmental pathways to the epidemic of Chronic Kidney Disease of unknown etiology (CKDu)"

Jayasinghe et al. [Science of the Total Environment, 705 (2020) 135766] propose a 'dynamical' model of Chronic Kidney Disease of Unknown etiology (CKDu) wherein CKDu arises as an emergent property of a complex system where they claim that weak multiple environmental factors contribute. They criticize the usual approaches as being "reductionist". We use their model as a basis of a discussion on the possibility of treating CKDu as an emergent property resulting from the interaction of multiple weak environmental factors with the organism. The model does not reveal anything beyond what is already known from simple considerations of well-known feed-back loops, but has the merit of re-stating those issues in a different format. If a proper weighting of the possible environmental factors is included, most proposed environmental factors drop out and what Jayasinghe et al. call the "reductionist" approach naturally arises due to the weight of evidence. The theory that the consumption of water containing fluoride and magnesium ions as found in water from regolith aquifers drawn via house-hold wells is found to clearly hold within this model when proper weighting is included. However, we show by examples that such models can be easily misused, leading to completely misleading conclusions. A response formalism useful in the theory of complex systems and emergent modes is presented in the context of the current problem. In addition to there being a lack of adequate data to fully implement such a theory, it is seen that such elaborations are unnecessary and misleading in the present context.

preprint2020arXiv

Liquid-liquid Phase transitions in silicon

We use computationally simple neutral pseudo-atom (`average atom') one-center density functional theory (DFT) and standard N-center DFT-Molecular Dynamics simulations to elucidate liquid-liquid phase transitions (LPTs) in supercooled liquid silicon at 1200K ,and also in silicon under warm-dense matter conditions up to 11604k (1 eV). An ionization-driven transition and three LPTs including the known LPT near 2.5 g/cm$^3$ are found. They are robust even to 1 eV. The pair distributions functions, pair potentials, electrical conductivities, and compressibilites are reported. The LPTs are elucidated within a Fermi liquid picture of electron scattering at the Fermi energy that complements the usual transient covalent bonding picture.

preprint2016arXiv

Current issues in finite-$T$ density-functional theory and Warm-Correlated Matter

Finite-temperature DFT has become of topical interest, partly due to the increasing ability to create novel states of warm-correlated matter (WCM). Subclasses of WCM are Warm-dense matter (WDM), ultra-fast matter (UFM), and high-energy density matter (HEDM), containing electyrons (e) and ions (i). Strong e-e, i-i and e-i correlation effects and partial degeneracies are found in these systems where the electron temperature $T_e $ is comparable to the electron Fermi energy. The ion subsystem may be solid, liquid or plasma, with many states of ionization with ionic charge $Z_j$. Quasi-equilibria with the ion temperature $T_i\ne T_e$ are common. The ion subsystem in WCM can no longer be treated as a passive "external potential", as is customary in $T=0$ density functional theory (DFT) dominated by solid-state theory or quantum chemistry. Hohenberg-Kohn-Mermin theory can be used for WCMs if finite-$T$ exchange-correlation (XC) functionals are available. They are functionals of both the one-body electron density $n_e$ and the one-body ion densities $ρ_j$. A method of approximately but accurately mapping the quantum electrons to a classical Coulomb gas enables one to treat electron-ion systems entirely classically at any temperature and arbitrary spin polarization, using exchange-correlation effects calculated {\it in situ}, directly from the pair-distribution functions. This eliminates the need for any XC-functionals, or the use of a Born-Oppenheimer approximation. This classical map has been used to calculate the equation of state of WDM systems, and construct a finite-$T$ XC functional that is found to be in close agreement with recent quantum path-integral simulation data. In this review current developments and concerns in finite-$T$ DFT, especially in the context of non-relativistic warm-dense matter and ultra-fast matter will be presented.

preprint2016arXiv

The dynamic conductivity and the plasmon profile of Aluminum in the ultra-fast-matter regime; an analysis of recent X-ray scattering data from the LCLS

We use an explicitly isochoric two-temperature theory where the ion temperature is not assumed to be equal to the electron temperature, to analyze recent X-ray laser scattering data for Aluminum in the ultra-fast-matter (UFM) regime up to 6 eV, obtained at the Standford Linac Coherent Light Source (LCLS). The observed surprisingly low conductivities are explained by including strong electron-ion scattering effects using the phase shifts calculated via the density-functional neutral-pseudo-atom model (NPA). The applicability of the Mermin model to UFM is questioned. The static and dynamic conductivity, complex collision frequency and the plasmon line-shape are evaluated within a Born approximation and without the use of the Mermin response. This approach, where all needed quantities are obtained from the NPA model, provides good agreement with the reported experimental results.

preprint2014arXiv

Strongly-Coupled Coulomb Systems using finite-$T$ Density Functional Theory: A review of studies on Strongly-Coupled Coulomb Systems since the rise of DFT and SCCS-1977

The conferences on "Strongly Coupled Coulomb Systems" (SCCS) arose from the "Strongly Coupled Plasmas" meetings, inaugurated in 1977. The progress in SCCS theory is reviewed in an `author-centered' frame to limit its scope. Our efforts, i.e., with François Perrot, sought to apply density functional theory (DFT) to SCCS calculations. DFT was then poised to become the major computational scheme for condensed matter physics. The ion-sphere models of Salpeter and others evolved into useful average-atom models for finite-$T$ Coulomb systems, as in Lieberman's Inferno code. We replaced these by correlation-sphere models that exploit the description of matter via density functionals linked to pair-distributions. These methods provided practical computational means for studying strongly interacting electron-ion Coulomb systems like warm-dense matter (WDM). The staples of SCCS are wide-ranged, viz., equation of state, plasma spectroscopy, opacity (absorption, emission), scattering, level shifts, transport properties, e.g., electrical and heat conductivity, laser- and shock- created plasmas, their energy relaxation and transient properties etc. These calculations need pseudopotentials and exchange-correlation functionals applicable to finite-$T$ Coulomb systems that may be used in ab initio codes, molecular dynamics, etc. The search for simpler computational schemes has proceeded via proposals for orbital-free DFT, statistical potentials, classical maps of quantum systems using classical schemes like HNC to include strong coupling effects (CHNC). Laughlin's classical plasma map for the fractional quantum Hall effect (FQHE) is a seminal example where we report new results for graphene.

preprint2013arXiv

Density-Functional theory, finite-temperature classical maps, and their implications for foundational studies of quantum systems

The advent of the Hohenberg-Kohn theorem in 1964, its extension to finite-T, Kohn-Sham theory, and relativistic extensions provide the well-established formalism of density-functional theory (DFT). This theory enables the calculation of all static properties of quantum systems {\it without} the need for an n-body wavefunction ψ. DFT uses the one-body density distribution instead of ψ. The more recent time-dependent formulations of DFT attempt to describe the time evolution of quantum systems without using the time-dependent wavefunction. Although DFT has become the standard tool of condensed-matter computational quantum mechanics, its foundational implications have remained largely unexplored. While all systems require quantum mechanics (QM) at T=0, the pair-distribution functions (PDFs) of such quantum systems have been accurately mapped into classical models at effective finite-T, and using suitable non-local quantum potentials (e.g., to mimic Pauli exclusion effects). These approaches shed light on the quantum \to hybrid \to classical models, and provide a new way of looking at the existence of non- local correlations without appealing to Bell's theorem. They also provide insights regarding Bohmian mechanics. Furthermore, macroscopic systems even at 1 Kelvin have de Broglie wavelengths in the micro-femtometer range, thereby eliminating macroscopic cat states, and avoiding the need for {\it ad hoc} decoherence models.

preprint2012arXiv

Electron-ion and ion-ion potentials for modeling warm-dense-matter: applications to laser-heated or shock-compressed Al and Si

The pair-interactions U_{ij}(r) determine the thermodynamics and linear transport properties of matter via the pair-distribution functions (PDFs), i.e., g_{ij}(r). Great simplicity is achieved if U_{ij}(r) could be directly used to predict material properties via classical simulations, avoiding many-body wavefunctions. Warm dense matter (WDM) is encountered in quasi-equilibria where the electron temperature $T_e$ differs from the ion temperature T_i, as in laser-heated or in shock-compressed matter. The electron PDFs g_{ee}(r) as perturbed by the ions are used to evaluate fully non-local exchange-correlation corrections to the free energy, using Hydrogen as an example. Electron-ion potentials for ions with a bound core are discussed with Al and Si as examples, for WDM with T_e \ne T_i, and valid for times shorter than the electron-ion relaxation time. In some cases the potentials develop attractive regions, and then become repulsive and `Yukawa-like' for higher $T_e$. These results clarify the origin of initial phonon-hardening and rapid release. Pair-potentials for shock-heated WDM show that phonon hardening would not occur in most such systems. Defining meaningful quasi-equilibrium static transport coefficients consistent with the dynamic values is addressed. There seems to be no meaningful `static conductivity' obtainable by extrapolating experimental or theoretical σ(ω, T_i, T_e) to ω\to 0, unless T_i \to T_e as well. Illustrative calculations of quasi-static resistivities R(T_i,T_e) of laser-heated as well as shock-heated Aluminum and Silicon are presented using our pseudopotentials, pair-potentials and classical integral equations. The quasi-static resistivities display clear differences in their temperature evolutions, but are not the strict ω\to 0 limits of the dynamic values.

preprint2011arXiv

The Classical-Map Hyper-Netted-Chain (CHNC) method and associated novel density-functional methods for Warm Dense Matter

The advent of short-pulse lasers, nanotechnology, as well as shock-wave techniques have created new states of matter (e.g., warm dense matter) that call for new theoretical tools. Ion correlations, electron correlations as well as bound states, continuum states, partial degeneracies and quasi-equilibrium systems need to be addressed. Bogoliubov's ideas of timescales can be used to discuss the quasi-thermodynamics of non-equilibrium systems. A rigorous approach to the associated many-body problem turns out to be the computation of the underlying pair-distribution functions g_ee, g_ei and g_ii, that directly yield non-local exchange-correlation potentials, free energies etc., valid within the timescales of each evolving system. An accurate classical map of the strongly-quantum uniform electron-gas problem given by Dharma-wardana and Perrot is reviewed. This replaces the quantum electrons at T=0 by an equivalent classical fluid at a finite temperature T_q, and having the same correlation energy. It has been shown, but not proven, that the classical fluid g_ij are excellent approximations to the quantum g_ij. The classical map is used with classical molecular dynamics (CMMD) or hyper-netted-chain integral equations (CHNC) to determine the pair-distribution functions (PDFs), and hence their thermodynamic and linear transport properties. The CHNC is very efficient for calculating the PDFs of uniform systems, while CMMD is more adapted to non-uniform systems. Applications to 2D and 3D quantum fluids, Si metal-oxide-field-effect transistors, Al plasmas, shock-compressed deuterium, two-temperature plasmas, pseudopotentials, as well as calculations for parabolic quantum dots are reviewed.

preprint2011arXiv

The electron effective mass in the strongly correlated 2D-uniform electron fluid from finite-temperature calculations

The very-low temperature thermal effective mass m* of paramagnetic and ferromagnetic electrons in a uniform electron fluid in two dimensions is studied. Analytical and numerical evaluations are used to meaningfully define an m*, even in the the Hartree-Fock approximation. The Hartree-Fock m* decreases linearly with the electron-disk radius r_s. Correlation effects lead to strong cancellations between exchange and correlation. Thus the effective mass is enhanced with increasing r_s for the unpolarized fluid, while m* decreases with the r_s of the polarized fluid. The effective mass is calculated from the coefficient of the quadratic temperature dependence of exchange-correlation free energy F_{xc}. This is calculated in a physically transparent manner using a new formula for the effective mass. This uses the T=0 pair-distribution functions of Gori-Giorgi et al., and the temperature derivative of a quantum analogue of the potential of mean-force well known in the statistical mechanics of classical fluids. The results are compared with recent quantum Monte-Carlo simulations at T=0, as well as with other available experimental and theoretical data for the effective mass.

preprint2010arXiv

Surface-directed spinodal decomposition in the pseudobinary alloy (HfO_2)_x(SiO_2)_{1-x}

Hf silicate films (HfO_2)_{0.25}(SiO_2)_{0.75} with thicknesses in the range 4-20 nm were grown on silicon substrate by atomic layer deposition at 350 deg.C.The Hf distributions in as-grown and 800 deg.C annealed films were investigated by high resolution transmission electron microscopy (HRTEM), angle-resolved x-ray photoelectron spectroscopy (ARXPS) and medium energy ion scattering (MEIS). HRTEM images show a layered structure in films thinner than 8 nm. The ARXPS data also reveal a non-uniform distribution of Hf throughout the film depth. Diffusion of SiO_2 to the film surface after a longer time anneal was observed by MEIS. All these observations provide evidence for surface-directed spinodal decomposition in the pseudobinary (HfO_2)_x(SiO_2)_{1-x} alloy system.

preprint2010arXiv

The Coulomb bridge function and the Pair-distribution functions of the 2-dimensional electron liquid in the quantum regime

The electron-electron pair distribution functions (PDF) of the 2-D electron fluid (2DEF) in the quantum regime (at T=0) are calculated using a classical-map-hyper-netted-chain (CHNC) scheme and compared with currently available Quantum Monte-Carlo (QMC) simulations in the coupling range r_s=1 to 50. We iteratively extract the bridge function of the "equivalent" classical 2-D liquid in the quantum regime. These bridge functions B(r) are relatively insensitive to spin-polarization effects. The structure of the bridge functions changes significantly for r_s>6, suggesting the onset of strongly correlated clusters. The new B(r), appropriate for the long-range Coulomb potential, can be used to replace the hard-sphere B(r) previously used in these calculations. They provide accurate classical representations of the QMC-PDFs even at very strong coupling, and probably at finite-T near T=0.

preprint2009arXiv

The Classical-Map Hyper-Netted-Chain (CHNC) technique for inhomogeneous electron systems. Application to quantum dots

The Classical-map Hyper-Netted-Chain (CHNC) technique is a simple method of calculating quantum pair-distribution functions, spin-dependent energies, etc., of strongly-interacting {\it uniform} systems. We present CHNC calculations of charge densities and energies of {\it non-uniform} systems, viz., quantum dots, and compare with quantum Monte Carlo and density -functional results. Results for up to 210 electrons are reported.

preprint2007arXiv

Temperature Relaxation in Hot Dense Hydrogen

Temperature equilibration of hydrogen is studied for conditions relevant to inertial confinement fusion. New molecular-dynamics simulations and results from quantum many-body theory are compared with Landau-Spitzer (LS) predictions for temperatures T from 50 eV to 5000 eV, and densities with Wigner-Seitz radii r_s = 1.0 and 0.5. The relaxation is slower than the LS result, even for temperatures in the keV range, but converges to agreement in the high-T limit.

preprint2004arXiv

Local environment of Nitrogen in GaN{y}As{1-y} epilayers on GaAs (001) studied using X-ray absorption near edge spectroscopy

X-ray absorption near-edge spectroscopy (XANES) is used to study the N environment in bulk GaN and in GaN{y}As{1-y} epilayers on GaAs (001), for y \~5%. Density-functional optimized structures were used to predict XANES via multiple-scattering theory. We obtain striking agreement for pure GaN. An alloy model with nitrogen pairs on Ga accurately predicts the threshold energy, the width of the XANES ``white line'', and features above threshold, for the given X-ray polarization. The presence of N-pairs may point to a role for molecular N_2 in epitaxial growth kinetics.

preprint2004arXiv

The effective mass and g-factor of the strongly correlated 2-D electron fluid. Evidence for a coupled-valley condensate in the Si system

The effective mass m*, and the Lande g-factor of the uniform 2-D electron fluid (2DEF) are calculated as a function of the spin polarization zeta, and the density parameter r_s, using a non-perturbative analytic approach. Our theory is in good accord with the m*g* data of Zhu et al. for zeta=0 for the GaAs-2DEF, and striking agreement with the data of Shashkin et al for the Si-2DEF. While g* is enhanced in GaAs, m* is enhanced in Si. The latter arises from singlet-pair excitations in the two valleys forming a coupled-valley state occurring at the critical density of ~1.10^{11}$ e/cm^2.