Orbital-Free Density Functional Theory: Kinetic Potentials and Ab-Initio Local Pseudopotentials

In the density functional (DF) theory of Kohn and Sham, the kinetic energy of the ground state of a system of noninteracting electrons in a general external field is calculated using a set of orbitals. Orbital free methods attempt to calculate this directly from the electron density by approximating the universal but unknown kinetic energy density functional. However simple local approximations are inaccurate and it has proved very difficult to devise generally accurate nonlocal approximations. We focus instead on the kinetic potential, the functional derivative of the kinetic energy DF, which appears in the Euler equation for the electron density. We argue that the kinetic potential is more local and more amenable to simple physically motivated approximations in many relevant cases, and describe two pathways by which the value of the kinetic energy can be efficiently calculated. We propose two nonlocal orbital free kinetic potentials that reduce to known exact forms for both slowly varying and rapidly varying perturbations and also reproduce exact results for the linear response of the density of the homogeneous system to small perturbations. A simple and systematic approach for generating accurate and weak ab-initio local pseudopotentials which produce a smooth slowly varying valence component of the electron density is proposed for use in orbital free DF calculations of molecules and solids. The use of these local pseudopotentials further minimizes the possible errors from the kinetic potentials. Our theory yields results for the total energies and ionization energies of atoms, and for the shell structure in the atomic radial density profiles that are in very good agreement with calculations using the full Kohn-Sham theory.Comment: To be published in Phys. Rev.

Explicit characterization of the identity configuration in an Abelian Sandpile Model

Since the work of Creutz, identifying the group identities for the Abelian Sandpile Model (ASM) on a given lattice is a puzzling issue: on rectangular portions of Z^2 complex quasi-self-similar structures arise. We study the ASM on the square lattice, in different geometries, and a variant with directed edges. Cylinders, through their extra symmetry, allow an easy determination of the identity, which is a homogeneous function. The directed variant on square geometry shows a remarkable exact structure, asymptotically self-similar.Comment: 11 pages, 8 figure

Abelian Sandpile Model on the Honeycomb Lattice

We check the universality properties of the two-dimensional Abelian sandpile model by computing some of its properties on the honeycomb lattice. Exact expressions for unit height correlation functions in presence of boundaries and for different boundary conditions are derived. Also, we study the statistics of the boundaries of avalanche waves by using the theory of SLE and suggest that these curves are conformally invariant and described by SLE2.Comment: 24 pages, 5 figure

Three-leg correlations in the two component spanning tree on the upper half-plane

We present a detailed asymptotic analysis of correlation functions for the two component spanning tree on the two-dimensional lattice when one component contains three paths connecting vicinities of two fixed lattice sites at large distance $s$ apart. We extend the known result for correlations on the plane to the case of the upper half-plane with closed and open boundary conditions. We found asymptotics of correlations for distance $r$ from the boundary to one of the fixed lattice sites for the cases $r\gg s \gg 1$ and $s \gg r \gg 1$.Comment: 16 pages, 5 figure

When Models Interact with their Subjects: The Dynamics of Model Aware Systems

A scientific model need not be a passive and static descriptor of its subject. If the subject is affected by the model, the model must be updated to explain its affected subject. In this study, two models regarding the dynamics of model aware systems are presented. The first explores the behavior of "prediction seeking" (PSP) and "prediction avoiding" (PAP) populations under the influence of a model that describes them. The second explores the publishing behavior of a group of experimentalists coupled to a model by means of confirmation bias. It is found that model aware systems can exhibit convergent random or oscillatory behavior and display universal 1/f noise. A numerical simulation of the physical experimentalists is compared with actual publications of neutron life time and {\Lambda} mass measurements and is in good quantitative agreement.Comment: Accepted for publication in PLoS-ON

Aharonov-Bohm Oscillations with Spin: Evidence for Berry's Phase

We report a study of the Aharonov-Bohm effect, the oscillations of the resistance of a mesoscopic ring as a function of a perpendicular magnetic field, in a GaAs two-dimensional hole system with a strong spin-orbit interaction. The Fourier spectra of the oscillations reveal extra structure near the main peak whose frequency corresponds to the magnetic flux enclosed by the ring. A comparison of the experimental data with results of simulations demonstrates that the origin of the extra structure is the geometric (Berry) phase acquired by the carrier spin as it travels around the ring.Comment: To be published in Physical Review Letter

From boundary to bulk in logarithmic CFT

The analogue of the charge-conjugation modular invariant for rational logarithmic conformal field theories is constructed. This is done by reconstructing the bulk spectrum from a simple boundary condition (the analogue of the Cardy `identity brane'). We apply the general method to the c_1,p triplet models and reproduce the previously known bulk theory for p=2 at c=-2. For general p we verify that the resulting partition functions are modular invariant. We also construct the complete set of 2p boundary states, and confirm that the identity brane from which we started indeed exists. As a by-product we obtain a logarithmic version of the Verlinde formula for the c_1,p triplet models.Comment: 35 pages, 2 figures; v2: minor corrections, version to appear in J.Phys.

Logarithmic two-point correlators in the Abelian sandpile model

We present the detailed calculations of the asymptotics of two-site correlation functions for height variables in the two-dimensional Abelian sandpile model. By using combinatorial methods for the enumeration of spanning trees, we extend the well-known result for the correlation $\sigma_{1,1} \simeq 1/r^4$ of minimal heights $h_1=h_2=1$ to $\sigma_{1,h} = P_{1,h}-P_1P_h$ for height values $h=2,3,4$. These results confirm the dominant logarithmic behaviour $\sigma_{1,h} \simeq (c_h\log r + d_h)/r^4 + {\cal O}(r^{-5})$ for large $r$, predicted by logarithmic conformal field theory based on field identifications obtained previously. We obtain, from our lattice calculations, the explicit values for the coefficients $c_h$ and $d_h$ (the latter are new).Comment: 28 page

Structural and doping effects in the half-metallic double perovskite A2A_2CrWO6_6

he structural, transport, magnetic and optical properties of the double perovskite $A_2$CrWO$_6$ with $A=\text{Sr, Ba, Ca}$ have been studied. By varying the alkaline earth ion on the $A$ site, the influence of steric effects on the Curie temperature $T_C$ and the saturation magnetization has been determined. A maximum $T_C=458$ K was found for Sr$_2$CrWO$_6$ having an almost undistorted perovskite structure with a tolerance factor $f\simeq 1$. For Ca$_2$CrWO$_6$ and Ba$_2$CrWO$_6$ structural changes result in a strong reduction of $T_C$. Our study strongly suggests that for the double perovskites in general an optimum $T_C$ is achieved only for $f \simeq 1$, that is, for an undistorted perovskite structure. Electron doping in Sr$_2$CrWO$_6$ by a partial substitution of Sr$^{2+}$ by La$^{3+}$ was found to reduce both $T_C$ and the saturation magnetization $M_s$. The reduction of $M_s$ could be attributed both to band structure effects and the Cr/W antisites induced by doping. Band structure calculations for Sr$_2$CrWO$_6$ predict an energy gap in the spin-up band, but a finite density of states for the spin-down band. The predictions of the band structure calculation are consistent with our optical measurements. Our experimental results support the presence of a kinetic energy driven mechanism in $A_2$CrWO$_6$, where ferromagnetism is stabilized by a hybridization of states of the nonmagnetic W-site positioned in between the high spin Cr-sites.Comment: 14 pages, 10 figure