Variational approach to dequantization

We present a dequantization procedure based on a variational approach whereby quantum fluctuations latent in the quantum momentum are suppressed. This is done by adding generic local deformations to the quantum momentum operator which give rise to a deformed kinetic term quantifying the amount of ``fuzzyness'' caused by such fluctuations. Considered as a functional of such deformations, the deformed kinetic term is shown to possess a unique minimum which is seen to be the classical kinetic energy. Furthermore, we show that extremization of the associated deformed action functional introduces an essential nonlinearity to the resulting field equations which are seen to be the classical Hamilton-Jacobi and continuity equations. Thus, a variational procedure determines the particular deformation that has the effect of suppressing the quantum fluctuations, resulting in dequantization of the system.Comment: 6 pages, 1 figure. v2: changes in presentation and conten

Classical-path integral adaptive resolution in molecular simulation: towards a smooth quantum-classical coupling

Simulations that couple different classical molecular models in an adaptive way by changing the number of degrees of freedom on the fly, are available within reasonably consistent theoretical frameworks. The same does not occur when it comes to classical-quantum adaptivity. The main reason for this is the difficulty in describing a continuous transition between the two different kind of physical principles: probabilistic for the quantum and deterministic for the classical. Here we report the basic principles of an algorithm that allows for a continuous and smooth transition by employing the path integral description of atoms.Comment: 8 pages 4 figure

An analytic expression for the electronic correlation term of the kinetic functional

We propose an analytic formula for the non-local Fisher information functional, or electronic kinetic correlation term, appearing in the expression of the kinetic density functional. Such an explicit formula is constructed on the basis of well-founded physical arguments and a rigorous mathematical prescription

Levy-Lieb-Based Monte Carlo Study of the Dimensionality Behaviour of the Electronic Kinetic Functional

We consider a gas of interacting electrons in the limit of nearly uniform density and treat the one dimensional (1D), two dimensional (2D) and three dimensional (3D) cases. We focus on the determination of the correlation part of the kinetic functional by employing a Monte Carlo sampling technique of electrons in space based on an analytic derivation via the Levy-Lieb constrained search principle. Of particular interest is the question of the behaviour of the functional as one passes from 1D to 3D; according to the basic principles of Density Functional Theory (DFT) the form of the universal functional should be independent of the dimensionality. However, in practice the straightforward use of current approximate functionals in different dimensions is problematic. Here, we show that going from the 3D to the 2D case the functional form is consistent (concave function) but in 1D becomes convex; such a drastic difference is peculiar of 1D electron systems as it is for other quantities. Given the interesting behaviour of the functional, this study represents a basic first-principle approach to the problem and suggests further investigations using highly accurate (though expensive) many-electron computational techniques, such as Quantum Monte Carlo

Polymers near Metal Surfaces: Selective Adsorption and Global Conformations

We study the properties of a polycarbonate melt near a nickel surface as a model system for the interaction of polymers with metal surfaces by employing a multiscale modeling approach. For bulk properties a suitably coarse grained bead spring model is simulated by molecular dynamics (MD) methods with model parameters directly derived from quantum chemical calculations. The surface interactions are parameterized and incorporated by extensive quantum mechanical density functional calculations using the Car-Parrinello method. We find strong chemisorption of chain ends, resulting in significant modifications of the melt composition when compared to an inert wall.Comment: 8 pages, 3 figures (2 color), 1 tabl

Adaptive Resolution Molecular Dynamics Simulation: Changing the Degrees of Freedom on the Fly

We present a new adaptive resolution technique for efficient particle-based multiscale molecular dynamics (MD) simulations. The presented approach is tailor-made for molecular systems where atomistic resolution is required only in spatially localized domains whereas a lower mesoscopic level of detail is sufficient for the rest of the system. Our method allows an on-the-fly interchange between a given molecule's atomic and coarse-grained level of description, enabling us to reach large length and time scales while spatially retaining atomistic details of the system. The new approach is tested on a model system of a liquid of tetrahedral molecules. The simulation box is divided into two regions: one containing only atomistically resolved tetrahedral molecules, the other containing only one particle coarse-grained spherical molecules. The molecules can freely move between the two regions while changing their level of resolution accordingly. The coarse-grained and the atomistically resolved systems have the same statistical properties at the same physical conditions.Comment: 17 pages, 11 figures, 5 table

Adaptive Resolution Molecular Dynamics Technique: Down to the Essential

We investigate the role of the thermodynamic (TD) force, as an essential and sufficient technical ingredient for an efficient and accurate adaptive resolution algorithm. Such a force applied in the coupling region of an adaptive resolution Molecular Dynamics (MD) set-up, assures thermodynamic equilibrium between atomistically resolved and coarse-grained regions, allowing the proper exchange of molecules. We numerically prove that indeed for systems as relevant as liquid water and 1,3-dimethylimidazolium chloride ionic liquid, the combined action of the TD force and thermostat allows for computationally efficient and numerically accurate simulations, beyond the current capabilities of adaptive resolution set-ups, which employ switching functions in the coupling region

Adaptive Resolution Simulation of Liquid Water

We present a multiscale simulation of liquid water where a spatially adaptive molecular resolution procedure allows for changing on-the-fly from a coarse-grained to an all-atom representation. We show that this approach leads to the correct description of all essential thermodynamic and structural properties of liquid water.Comment: 4 pages, 3 figures; changed figure

On the upper bound of the electronic kinetic energy in terms of density functionals

We propose a simple density functional expression for the upper bound of the kinetic energy for electronic systems. Such a functional is valid in the limit of slowly varying density, its validity outside this regime is discussed by making a comparison with upper bounds obtained in previous work. The advantages of the functional proposed for applications to realistic systems is briefly discussed.Comment: 10 pages, no figure

Classical kinetic energy, quantum fluctuation terms and kinetic-energy functionals

We employ a recently formulated dequantization procedure to obtain an exact expression for the kinetic energy which is applicable to all kinetic-energy functionals. We express the kinetic energy of an N-electron system as the sum of an N-electron classical kinetic energy and an N-electron purely quantum kinetic energy arising from the quantum fluctuations that turn the classical momentum into the quantum momentum. This leads to an interesting analogy with Nelson's stochastic approach to quantum mechanics, which we use to conceptually clarify the physical nature of part of the kinetic-energy functional in terms of statistical fluctuations and in direct correspondence with Fisher Information Theory. We show that the N-electron purely quantum kinetic energy can be written as the sum of the (one-electron) Weizsacker term and an (N-1)-electron kinetic correlation term. We further show that the Weizsacker term results from local fluctuations while the kinetic correlation term results from the nonlocal fluctuations. For one-electron orbitals (where kinetic correlation is neglected) we obtain an exact (albeit impractical) expression for the noninteracting kinetic energy as the sum of the classical kinetic energy and the Weizsacker term. The classical kinetic energy is seen to be explicitly dependent on the electron phase and this has implications for the development of accurate orbital-free kinetic-energy functionals. Also, there is a direct connection between the classical kinetic energy and the angular momentum and, across a row of the periodic table, the classical kinetic energy component of the noninteracting kinetic energy generally increases as Z increases.Comment: 10 pages, 1 figure. To appear in Theor Chem Ac