47 results on '"REAL-SPACE"'
Search Results
2. SPARC: Simulation Package for Ab-initio Real-space Calculations
- Author
-
Qimen Xu, Abhiraj Sharma, Benjamin Comer, Hua Huang, Edmond Chow, Andrew J. Medford, John E. Pask, and Phanish Suryanarayana
- Subjects
Kohn–Sham ,Density functional theory ,Electronic structure ,Real-space ,Finite-differences ,Computer software ,QA76.75-76.765 - Abstract
We present SPARC: Simulation Package for Ab-initio Real-space Calculations. SPARC can perform Kohn–Sham density functional theory calculations for isolated systems such as molecules as well as extended systems such as crystals and surfaces, in both static and dynamic settings. It is straightforward to install/use and highly competitive with state-of-the-art planewave codes, demonstrating comparable performance on a small number of processors and increasing advantages as the number of processors grows. Notably, SPARC brings solution times down to a few seconds for systems with O(100–500)atoms on large-scale parallel computers, outperforming planewave counterparts by an order of magnitude and more.
- Published
- 2021
- Full Text
- View/download PDF
3. Implementation of Laplace Transformed MP2 for Periodic Systems With Numerical Atomic Orbitals
- Author
-
Honghui Shang and Jinlong Yang
- Subjects
MP2 ,NAO ,real-space ,Hartree–Fock ,periodic system ,Chemistry ,QD1-999 - Abstract
We present an implementation of the canonical and Laplace-transformed formulation of the second-order Møller–Plesset perturbation theory under periodic boundary conditions using numerical atomic orbitals. To validate our approach, we show that our results of the Laplace-transformed MP2 correlation correction for the total energy and the band gap are in excellent agreement with the results of the canonical MP2 formulation. We have calculated the binding energy curve for the stacked trans-polyacetylene at the Hartree–Fock + MP2 level as a preliminary application.
- Published
- 2020
- Full Text
- View/download PDF
4. M-SPARC: Matlab-Simulation Package for Ab-initio Real-space Calculations
- Author
-
Qimen Xu, Abhiraj Sharma, and Phanish Suryanarayana
- Subjects
Kohn–Sham ,Density Functional Theory ,Electronic structure ,Real-space ,Matlab ,Computer software ,QA76.75-76.765 - Abstract
We present M-SPARC: Matlab-Simulation Package for Ab-initio Real-space Calculations. It can perform pseudopotential spin-polarized and unpolarized Kohn–Sham Density Functional Theory (DFT) simulations for isolated systems such as molecules as well as extended systems such as crystals, surfaces, and nanowires. M-SPARC provides a rapid prototyping platform for the development and testing of new algorithms and methods in real-space DFT, with the potential to significantly accelerate the rate of advancements in the field. It also provides a convenient avenue for the accurate first principles study of small to moderate sized systems.
- Published
- 2020
- Full Text
- View/download PDF
5. Crystallographic Structure Refinement in a Nutshell
- Author
-
Afonine, Pavel V., Adams, Paul D., Read, Randy, editor, Urzhumtsev, Alexandre G., editor, and Lunin, Vladimir Y., editor
- Published
- 2013
- Full Text
- View/download PDF
6. On real-space Density Functional Theory for non-orthogonal crystal systems: Kronecker product formulation of the kinetic energy operator.
- Author
-
Sharma, Abhiraj and Suryanarayana, Phanish
- Subjects
- *
DENSITY functional theory , *KRONECKER products , *KINETIC energy , *ELECTROSTATICS , *MULTIPLICATION - Abstract
We present an accurate and efficient real-space Density Functional Theory (DFT) framework for the ab initio study of non-orthogonal crystal systems. Specifically, employing a local reformulation of the electrostatics, we develop a novel Kronecker product formulation of the real-space kinetic energy operator that significantly reduces the number of operations associated with the Laplacian-vector multiplication, the dominant cost in practical computations. In particular, we reduce the scaling with respect to finite-difference order from quadratic to linear, thereby significantly bridging the gap in computational cost between non-orthogonal and orthogonal systems. We verify the accuracy and efficiency of the proposed methodology through selected examples. [ABSTRACT FROM AUTHOR]
- Published
- 2018
- Full Text
- View/download PDF
7. Quantitative real-space analysis of colloidal structures and dynamics with confocal scanning light microscopy
- Author
-
van Blaaderen, A., Kremer, F., editor, Lagaly, G., editor, Palberg, T., editor, and Ballauff, M., editor
- Published
- 1997
- Full Text
- View/download PDF
8. Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids
- Author
-
N. Hamaus, M. Aubert, A. Pisani, S. Contarini, G. Verza, M.-C. Cousinou, S. Escoffier, A. Hawken, G. Lavaux, G. Pollina, B. D. Wandelt, J. Weller, M. Bonici, C. Carbone, L. Guzzo, A. Kovacs, F. Marulli, E. Massara, L. Moscardini, P. Ntelis, W. J. Percival, S. Radinović, M. Sahlén, Z. Sakr, A. G. Sánchez, H. A. Winther, N. Auricchio, S. Awan, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou, G. Congedo, L. Conversi, Y. Copin, L. Corcione, M. Cropper, A. Da Silva, H. Degaudenzi, M. Douspis, F. Dubath, C. A. J. Duncan, X. Dupac, S. Dusini, A. Ealet, S. Ferriol, P. Fosalba, M. Frailis, E. Franceschi, P. Franzetti, M. Fumana, B. Garilli, B. Gillis, C. Giocoli, A. Grazian, F. Grupp, S. V. H. Haugan, W. Holmes, F. Hormuth, K. Jahnke, S. Kermiche, A. Kiessling, M. Kilbinger, T. Kitching, M. Kümmel, M. Kunz, H. Kurki-Suonio, S. Ligori, P. B. Lilje, I. Lloro, E. Maiorano, O. Marggraf, K. Markovic, R. Massey, S. Maurogordato, M. Melchior, M. Meneghetti, G. Meylan, M. Moresco, E. Munari, S. M. Niemi, C. Padilla, S. Paltani, F. Pasian, K. Pedersen, V. Pettorino, S. Pires, M. Poncet, L. Popa, L. Pozzetti, R. Rebolo, J. Rhodes, H. Rix, M. Roncarelli, E. Rossetti, R. Saglia, P. Schneider, A. Secroun, G. Seidel, S. Serrano, C. Sirignano, G. Sirri, J.-L. Starck, P. Tallada-Crespí, D. Tavagnacco, A. N. Taylor, I. Tereno, R. Toledo-Moreo, F. Torradeflot, E. A. Valentijn, L. Valenziano, Y. Wang, N. Welikala, G. Zamorani, J. Zoubian, S. Andreon, M. Baldi, S. Camera, S. Mei, C. Neissner, E. Romelli, Hamaus, N., Aubert, M., Pisani, A., Contarini, S., Verza, G., Cousinou, M. -C., Escoffier, S., Hawken, A., Lavaux, G., Pollina, G., Wandelt, B. D., Weller, J., Bonici, M., Carbone, C., Guzzo, L., Kovacs, A., Marulli, F., Massara, E., Moscardini, L., Ntelis, P., Percival, W. J., Radinovic, S., Sahlen, M., Sakr, Z., Sanchez, A. G., Winther, H. A., Auricchio, N., Awan, S., Bender, R., Bodendorf, C., Bonino, D., Branchini, E., Brescia, M., Brinchmann, J., Capobianco, V., Carretero, J., Castander, F. J., Castellano, M., Cavuoti, S., Cimatti, A., Cledassou, R., Congedo, G., Conversi, L., Copin, Y., Corcione, L., Cropper, M., Da Silva, A., Degaudenzi, H., Douspis, M., Dubath, F., Duncan, C. A. J., Dupac, X., Dusini, S., Ealet, A., Ferriol, S., Fosalba, P., Frailis, M., Franceschi, E., Franzetti, P., Fumana, M., Garilli, B., Gillis, B., Giocoli, C., Grazian, A., Grupp, F., Haugan, S. V. H., Holmes, W., Hormuth, F., Jahnke, K., Kermiche, S., Kiessling, A., Kilbinger, M., Kitching, T., Kummel, M., Kunz, M., Kurki-Suonio, H., Ligori, S., Lilje, P. B., Lloro, I., Maiorano, E., Marggraf, O., Markovic, K., Massey, R., Maurogordato, S., Melchior, M., Meneghetti, M., Meylan, G., Moresco, M., Munari, E., Niemi, S. M., Padilla, C., Paltani, S., Pasian, F., Pedersen, K., Pettorino, V., Pires, S., Poncet, M., Popa, L., Pozzetti, L., Rebolo, R., Rhodes, J., Rix, H., Roncarelli, M., Rossetti, E., Saglia, R., Schneider, P., Secroun, A., Seidel, G., Serrano, S., Sirignano, C., Sirri, G., Starck, J. -L., Tallada-Crespi, P., Tavagnacco, D., Taylor, A. N., Tereno, I., Toledo-Moreo, R., Torradeflot, F., Valentijn, E. A., Valenziano, L., Wang, Y., Welikala, N., Zamorani, G., Zoubian, J., Andreon, S., Baldi, M., Camera, S., Mei, S., Neissner, C., Romelli, E., Department of Physics, Helsinki Institute of Physics, Universitats-Sternwarte [München], Ludwig-Maximilians-Universität München (LMU), Institut de Physique des 2 Infinis de Lyon (IP2I Lyon), Université Claude Bernard Lyon 1 (UCBL), Université de Lyon-Université de Lyon-Institut National de Physique Nucléaire et de Physique des Particules du CNRS (IN2P3)-Centre National de la Recherche Scientifique (CNRS), Centre de Physique des Particules de Marseille (CPPM), Aix Marseille Université (AMU)-Institut National de Physique Nucléaire et de Physique des Particules du CNRS (IN2P3)-Centre National de la Recherche Scientifique (CNRS), Department of Astrophysical Sciences [Princeton], Princeton University, Dipartimento di Fisica e Astronomia [Bologna], Alma Mater Studiorum Università di Bologna [Bologna] (UNIBO), Istituto Nazionale di Fisica Nucleare, Sezione di Bologna (INFN, Sezione di Bologna), Istituto Nazionale di Fisica Nucleare (INFN), INAF - Osservatorio Astronomico di Bologna (OABO), Istituto Nazionale di Astrofisica (INAF), Istituto Nazionale di Fisica Nucleare, Sezione di Padova (INFN, Sezione di Padova), Dipartimento di Fisica e Astronomia 'Galileo Galilei', Università degli Studi di Padova = University of Padua (Unipd), Institut d'Astrophysique de Paris (IAP), Institut national des sciences de l'Univers (INSU - CNRS)-Sorbonne Université (SU)-Centre National de la Recherche Scientifique (CNRS), Max Planck Institute for Extraterrestrial Physics (MPE), Max-Planck-Gesellschaft, Istituto Nazionale di Fisica Nucleare, Sezione di Genova (INFN, Sezione di Genova), Università degli studi di Genova = University of Genoa (UniGe), Istituto Nazionale di Fisica Nucleare, Sezione di Milano (INFN), INAF-IASF Milano, INAF - Osservatorio Astronomico di Brera (OAB), Università degli Studi di Milano = University of Milan (UNIMI), Universidad de La Laguna [Tenerife - SP] (ULL), Instituto de Astrofisica de Canarias (IAC), University of Bologna/Università di Bologna, University of Waterloo [Waterloo], Department of Physics and Astronomy [Waterloo], Perimeter Institute for Theoretical Physics [Waterloo], Institute of Theoretical Astrophysics [Oslo], University of Oslo (UiO), Swedish Collegium for Advanced Study [Uppsala], Uppsala University, Institut de recherche en astrophysique et planétologie (IRAP), Université Toulouse III - Paul Sabatier (UT3), Université de Toulouse (UT)-Université de Toulouse (UT)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire Midi-Pyrénées (OMP), Institut de Recherche pour le Développement (IRD)-Université Toulouse III - Paul Sabatier (UT3), Université de Toulouse (UT)-Université de Toulouse (UT)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Météo-France -Institut de Recherche pour le Développement (IRD)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Météo-France -Centre National de la Recherche Scientifique (CNRS), Unité de Recherche Environnement, Génomique Fonctionnelle et Études Mathématiques [Beyrouth] (UR-EGFEM), Université Saint-Joseph de Beyrouth (USJ), Institut National de Physique Nucléaire et de Physique des Particules du CNRS (IN2P3), Centre National d'Études Spatiales [Toulouse] (CNES), Institut d'astrophysique spatiale (IAS), Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Centre National d’Études Spatiales [Paris] (CNES), Astrophysique Interprétation Modélisation (AIM (UMR_7158 / UMR_E_9005 / UM_112)), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Université Paris Cité (UPCité), Observatoire de la Côte d'Azur (OCA), Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS), AstroParticule et Cosmologie (APC (UMR_7164)), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Institut National de Physique Nucléaire et de Physique des Particules du CNRS (IN2P3)-Observatoire de Paris, Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Centre National de la Recherche Scientifique (CNRS)-Université Paris Cité (UPCité), ANR-16-CE23-0002,BIG4,Grosses données, Grosses simulations, Big Bang et Grands problèmes: Algorithes de reconstruction bayésiennes contraintes par la physique et application à l'analyse de données cosmologiques(2016), Agence Nationale de la Recherche (France), German Research Foundation, European Space Agency, National Aeronautics and Space Administration (US), Agenzia Spaziale Italiana, Institut National de Physique Nucléaire et de Physique des Particules du CNRS (IN2P3)-Université Claude Bernard Lyon 1 (UCBL), Université de Lyon-Université de Lyon-Centre National de la Recherche Scientifique (CNRS), Universita degli Studi di Padova, Universita degli studi di Genova, Università degli Studi di Milano [Milano] (UNIMI), Università di Bologna Dipartimento di Fisca e Astronomia, INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, University of Bologna, Institut national des sciences de l'Univers (INSU - CNRS)-Université Toulouse III - Paul Sabatier (UT3), Université Fédérale Toulouse Midi-Pyrénées-Université Fédérale Toulouse Midi-Pyrénées-Observatoire Midi-Pyrénées (OMP), Météo France-Centre National d'Études Spatiales [Toulouse] (CNES)-Université Fédérale Toulouse Midi-Pyrénées-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD)-Météo France-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD)-Centre National de la Recherche Scientifique (CNRS), Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS), Commissariat à l'énergie atomique et aux énergies alternatives (CEA)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris-Saclay-Centre National de la Recherche Scientifique (CNRS)-Université de Paris (UP), Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Centre National de la Recherche Scientifique (CNRS)-Université de Paris (UP), and Astronomy
- Subjects
Void (astronomy) ,Methods: data analysis / surveys ,Cosmological parameter ,Cosmology and Nongalactic Astrophysics (astro-ph.CO) ,Large-scale structure of Universe ,[SDU.ASTR.CO]Sciences of the Universe [physics]/Astrophysics [astro-ph]/Cosmology and Extra-Galactic Astrophysics [astro-ph.CO] ,Cosmological parameters ,FOS: Physical sciences ,Astrophysics ,Astrophysics::Cosmology and Extragalactic Astrophysics ,Surveys ,01 natural sciences ,Cosmology: observations ,Dark energy ,Methods: data analysis ,114 Physical sciences ,Cosmology: observation ,Cosmology ,Redshift-space distortions ,real-space ,0103 physical sciences ,Dark matter ,Large-scale structure of the Universe ,luminosity function ,observations [Cosmology] ,data analysis [Methods] ,010303 astronomy & astrophysics ,dark energy survey ,Physics ,survey cosmological implications ,galaxy troughs ,density ,COSMIC cancer database ,010308 nuclear & particles physics ,Angular diameter distance ,Astronomy and Astrophysics ,oscillation spectroscopic survey ,115 Astronomy, Space science ,Redshift ,Galaxy ,matter ,gravity ,Space and Planetary Science ,gravitational-instability ,Astrophysics - Cosmology and Nongalactic Astrophysics ,Methods: data analysi - Abstract
Euclid Consortium: N. Hamaus et al., Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the three-dimensional large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmological studies. In particular, we study the imprints of dynamic (redshift-space) and geometric (Alcock–Paczynski) distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11 000 voids and we estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on f/b and DMH, where f is the linear growth rate of density fluctuations, b the galaxy bias, DM the comoving angular diameter distance, and H the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects in the analysis. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of f/b and 0.5% on DMH in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter, w, for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts., NH, GP and JW are supported by the Excellence Cluster ORIGINS, which is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311. MA, MCC and SE are supported by the eBOSS ANR grant (under contract ANR-16-CE31-0021) of the French National Research Agency, the OCEVU LABEX (Grant No. ANR-11-LABX-0060) and the A*MIDEX project (Grant No. ANR-11-IDEX-0001-02) funded by the Investissements d’Avenir French government program, and by CNES, the French National Space Agency. AP is supported by NASA ROSES grant 12-EUCLID12-0004, and NASA grant 15-WFIRST15-0008 to the Nancy Grace Roman Space Telescope Science Investigation Team “Cosmology with the High Latitude Survey”. GL is supported by the ANR BIG4 project, grant ANR-16-CE23-0002 of the French Agence Nationale de la Recherche. PN is funded by the Centre National d’Etudes Spatiales (CNES). We acknowledge use of the Python libraries NumPy (Harris et al. 2020), SciPy (Virtanen et al. 2020), Matplotlib (Hunter 2007), Astropy (Astropy Collaboration 2013, 2018), emcee (Foreman-Mackey et al. 2019), GetDist (Lewis 2019), healpy (Górski et al. 2005; Zonca et al. 2019), and PyAbel (Hickstein et al. 2019). This work has made use of CosmoHub (Carretero et al. 2017; Tallada et al. 2020). CosmoHub has been developed by the Port d’Informació Científica (PIC), maintained through a collaboration of the Institut de Física d’Altes Energies (IFAE) and the Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT) and the Institute of Space Sciences (CSIC & IEEC), and was partially funded by the “Plan Estatal de Investigación Científica y Técnica y de Innovación” program of the Spanish government. The Euclid Consortium acknowledges the European Space Agency and a number of agencies and institutes that have supported the development of Euclid, in particular the Academy of Finland, the Agenzia Spaziale Italiana, the Belgian Science Policy, the Canadian Euclid Consortium, the Centre National d’Etudes Spatiales, the Deutsches Zentrum für Luft- und Raumfahrt, the Danish Space Research Institute, the Fundação para a Ciência e a Tecnologia, the Ministerio de Economia y Competitividad, the National Aeronautics and Space Administration, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Romanian Space Agency, the State Secretariat for Education, Research and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency.
- Published
- 2022
- Full Text
- View/download PDF
9. SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems.
- Author
-
Ghosh, Swarnava and Suryanarayana, Phanish
- Subjects
- *
DENSITY functional theory , *FINITE differences , *ELECTROSTATICS , *CHEBYSHEV polynomials , *NUCLEAR forces (Physics) , *ELECTRONIC structure , *PARALLEL processing - Abstract
As the second component of SPARC (Simulation Package for Ab-initio Real-space Calculations), we present an accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory (DFT) for extended systems. Specifically, employing a local formulation of the electrostatics, the Chebyshev polynomial filtered self-consistent field iteration, and a reformulation of the non-local force component, we develop a finite-difference framework wherein both the energy and atomic forces can be efficiently calculated to within desired accuracies in DFT. We demonstrate using a wide variety of materials systems that SPARC achieves high convergence rates in energy and forces with respect to spatial discretization to reference plane-wave result; exponential convergence in energies and forces with respect to vacuum size for slabs and wires; energies and forces that are consistent and display negligible ‘egg-box’ effect; accurate properties of crystals, slabs, and wires; and negligible drift in molecular dynamics simulations. We also demonstrate that the weak and strong scaling behavior of SPARC is similar to well-established and optimized plane-wave implementations for systems consisting up to thousands of electrons, but with a significantly reduced prefactor. Overall, SPARC represents an attractive alternative to plane-wave codes for performing DFT simulations of extended systems. Program summary Program title: SPARC Catalogue identifier: AFBR_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AFBR_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU GPL v3 No. of lines in distributed program, including test data, etc.: 93822 No. of bytes in distributed program, including test data, etc.: 1386659 Distribution format: tar.gz Programming language: C/C++. Computer: Any system with C/C++ compiler. Operating system: Linux. RAM: Problem dependent. Ranges from 80 GB to 800 GB for a system with 2500 electrons. Classification: 7.3. External routines: PETSc 3.5.3 ( http://www.mcs.anl.gov/petsc ), MKL 11.2 ( https://software.intel.com/en-us/intel-mkl ), and MVAPICH2 2.1 ( http://mvapich.cse.ohio-state.edu/ ). Does the new version supersede the previous version?: Yes Nature of problem: Calculation of the static and dynamic properties of isolated and extended systems in the framework of Kohn–Sham Density Functional Theory (DFT). Solution method: High-order finite-difference discretization. Local reformulation of the electrostatics in terms of the electrostatic potential and pseudocharge densities. Application of Bloch-periodic and zero-Dirichlet boundary conditions on the orbitals in the direction of periodicity and vacuum, respectively. Application of periodic and Dirichlet boundary conditions on the electrostatic potential in the direction of periodicity and vacuum, respectively. Integration over the Brillouin zone for extended systems using the Monkhorst–Pack grid. Calculation of the electronic ground-state using the Chebyshev polynomial filtered self-consistent field iteration in conjunction with Anderson based extrapolation/mixing schemes. Reformulation of the non-local component of the force. Geometry optimization using the Polak–Ribiere variant of non-linear conjugate gradients with secant line search. NVE molecular dynamics using the leapfrog method. Parallelization via domain decomposition and over Brillouin zone integration. Reasons for new version: To enable the study of extended systems like crystals, slabs, and wires using SPARC. Summary of revisions: Incorporated the ability to study the static and dynamic properties of crystals, slabs, and wires. Restrictions: System size less than ∼ 4000 electrons. Local Density Approximation (LDA). Troullier–Martins pseudopotentials without relativistic or non-linear core corrections. Domain has to be cuboidal. Running time: Problem dependent. Timing results for selected examples provided in the paper. [ABSTRACT FROM AUTHOR]
- Published
- 2017
- Full Text
- View/download PDF
10. SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters.
- Author
-
Ghosh, Swarnava and Suryanarayana, Phanish
- Subjects
- *
CHEBYSHEV polynomials , *NUCLEAR forces (Physics) , *ELECTROSTATICS , *FINITE difference method , *STOCHASTIC convergence , *VIBRATIONAL spectra , *DENSITY functional theory - Abstract
As the first component of SPARC (Simulation Package for Ab-initio Real-space Calculations), we present an accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory (DFT) for isolated clusters. Specifically, utilizing a local reformulation of the electrostatics, the Chebyshev polynomial filtered self-consistent field iteration, and a reformulation of the non-local component of the force, we develop a framework using the finite-difference representation that enables the efficient evaluation of energies and atomic forces to within the desired accuracies in DFT. Through selected examples consisting of a variety of elements, we demonstrate that SPARC obtains exponential convergence in energy and forces with domain size; systematic convergence in the energy and forces with mesh-size to reference plane-wave result at comparably high rates; forces that are consistent with the energy, both free from any noticeable ‘egg-box’ effect; and accurate ground-state properties including equilibrium geometries and vibrational spectra. In addition, for systems consisting up to thousands of electrons, SPARC displays weak and strong parallel scaling behavior that is similar to well-established and optimized plane-wave implementations, but with a significantly reduced prefactor. Overall, SPARC represents an attractive alternative to plane-wave codes for practical DFT simulations of isolated clusters. Program summary Program title: SPARC Catalogue identifier: AFBL_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AFBL_v1_0.html Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions: GNU GPL v3 No. of lines in distributed program, including test data, etc.: 47525 No. of bytes in distributed program, including test data, etc.: 826436 Distribution format: tar.gz Programming language: C/C++. Computer: Any system with C/C++ compiler. Operating system: Linux. RAM: Problem dependent. Ranges from 80 GB to 800 GB for a system with 2500 electrons. Classification: 7.3. External routines: PETSc 3.5.3 ( http://www.mcs.anl.gov/petsc ), MKL 11.2 ( https://software.intel.com/en-us/intel-mkl ), and MVAPICH2 2.1 ( http://mvapich.cse.ohio-state.edu/ ). Nature of problem: Calculation of the electronic and structural ground-states for isolated clusters in the framework of Kohn–Sham Density Functional Theory (DFT). Solution method: High-order finite-difference discretization. Local reformulation of the electrostatics in terms of the electrostatic potential and pseudocharge densities. Calculation of the electronic ground-state using the Chebyshev polynomial filtered Self-Consistent Field (SCF) iteration in conjunction with Anderson extrapolation/mixing. Evaluation of boundary conditions for the electrostatic potential through a truncated multipole expansion. Reformulation of the non-local component of the force. Geometry optimization using the Polak–Ribiere variant of non-linear conjugate gradients with secant line search. Restrictions: System size less than ∼4000 electrons. Local Density Approximation (LDA). Troullier–Martins pseudopotentials without relativistic or non-linear core corrections. Running time: Problem dependent. Timing results for selected examples provided in the paper. [ABSTRACT FROM AUTHOR]
- Published
- 2017
- Full Text
- View/download PDF
11. Monte Carlo Simulation of a Model Cuprate
- Author
-
Panov, Yu. D., Moskvin, A. S., Chikov, A. A., Ulitko, V. A., Panov, Yu. D., Moskvin, A. S., Chikov, A. A., and Ulitko, V. A.
- Abstract
We develop a classical Monte Carlo algorithm based on a quasi-classical approximation for a pseudospin S = 1 Hamiltonian in real space to construct a phase diagram of a model cuprate with a high Tc. A model description takes into account both local and nonlocal correlations, Heisenberg spin-exchange interaction, correlated single-particle, and two-particle transport. We formulate a state selection algorithm for a given parameterization of the wave function in order to ensure a uniform distribution of states in the phase space. The simulation results show a qualitative agreement with the experimental phase diagrams. © 2021 Institute of Physics Publishing. All rights reserved.
- Published
- 2021
12. Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory.
- Author
-
Ghosh, Swarnava and Suryanarayana, Phanish
- Subjects
- *
FINITE difference method , *DENSITY functional theory , *SELF-consistent field theory , *KINETIC energy , *STOCHASTIC convergence - Abstract
We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state and forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms. [ABSTRACT FROM AUTHOR]
- Published
- 2016
- Full Text
- View/download PDF
13. Multislice algorithms revisited: Solving the Schrödinger equation numerically for imaging with electrons.
- Author
-
Wacker, C. and Schröder, R.R.
- Subjects
- *
ALGORITHMS , *APPROXIMATION theory , *SIMULATION methods & models , *ELECTRON microscopy , *NUMERICAL analysis - Abstract
For a long time, the high-energy approximation was sufficient for any image simulation in electron microscopy. This changed with the advent of aberration correctors that allow high-resolution imaging at low electron energies. To deal with this fact, we present a numerical solution of the exact Schrödinger equation that is novel in the field of electron microscopy. Furthermore, we investigate systematically the advantages and problems of several multislice algorithms, especially the real-space algorithms. [ABSTRACT FROM AUTHOR]
- Published
- 2015
- Full Text
- View/download PDF
14. Augmented Lagrangian formulation of orbital-free density functional theory.
- Author
-
Suryanarayana, Phanish and Phanish, Deepa
- Subjects
- *
LAGRANGIAN functions , *DENSITY functional theory , *CONSTRAINTS (Physics) , *PLANE wavefronts , *LAGRANGE multiplier , *FINITE differences - Abstract
We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic Orbital-Free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a sequence of minimization problems without any constraint, thereby making it amenable to powerful unconstrained optimization algorithms. Further, we develop a parallel implementation of this approach for the Thomas–Fermi–von Weizsacker (TFW) kinetic energy functional in the framework of higher-order finite-differences and the conjugate gradient method. With this implementation, we establish that the Augmented Lagrangian approach is highly competitive compared to the penalty and Lagrange multiplier methods. Additionally, we show that higher-order finite-differences represent a computationally efficient discretization for performing OF-DFT simulations. Overall, we demonstrate that the proposed formulation and implementation are both efficient and robust by studying selected examples, including systems consisting of thousands of atoms. We validate the accuracy of the computed energies and forces by comparing them with those obtained by existing plane-wave methods. [ABSTRACT FROM AUTHOR]
- Published
- 2014
- Full Text
- View/download PDF
15. TTDFT: A GPU accelerated Tucker tensor DFT code for large-scale Kohn-Sham DFT calculations.
- Author
-
Lin, Chih-Chuen and Gavini, Vikram
- Subjects
- *
DENSITY functional theory , *GRAPHICS processing units , *CODING theory , *SYMMETRIC matrices , *QUANTUM dots , *PROGRAMMING languages - Abstract
We present the Tucker tensor DFT (TTDFT) code which uses a tensor-structured algorithm with graphic processing unit (GPU) acceleration for conducting ground-state DFT calculations on large-scale systems. The Tucker tensor DFT algorithm uses a localized Tucker tensor basis computed from an additive separable approximation to the Kohn-Sham Hamiltonian. The discrete Kohn-Sham problem is solved using Chebyshev filtered subspace iteration method that relies on matrix-matrix multiplications of a sparse symmetric Hamiltonian matrix and a dense wavefunction matrix, expressed in the localized Tucker tensor basis. These matrix-matrix multiplication operations, which constitute the most computationally intensive step of the solution procedure, are GPU accelerated providing ∼8-fold GPU-CPU speedup for these operations on the largest systems studied. The computational performance of the TTDFT code is presented using benchmark studies on aluminum nano-particles and silicon quantum dots with system sizes ranging up to ∼7,000 atoms. Program Title: TTDFT: Tucker tensor density functional theory code CPC Library link to program files: https://doi.org/10.17632/8dgmcs8ys2.1 Licensing provisions: LGPL Programming language: C/C++ External routines/libraries: TuckerMPI (https://gitlab.com/tensors/TuckerMPI), cuBLAS (https://docs.nvidia.com/cuda/cublas/index.html), cuSparse (https://docs.nvidia.com/cuda/cusparse/index.html), ALGLIB (http://www.alglib.net/), Boost (https://www.boost.org/), BLAS (http://www.netlib.org/blas/), LAPACK (http://www.netlib.org/lapack/), PETSc (https://www.mcs.anl.gov/petsc), SLEPc (http://slepc.upv.es) Nature of problem: Real-space Kohn-Sham density functional theory calculations using localized Tucker tensor basis. Solution method: We present a real-space Kohn-Sham density functional theory code based on tensor-structured techniques with GPU acceleration. Tensor-structured techniques are adopted for computing a Tucker tensor basis, representing the eigenfunctions of an additive separable approximation to the Kohn-Sham Hamiltonian. The Tucker tensor basis is further localized using L 1 regularization to improve the sparsity of the Kohn-Sham Hamiltonian matrix, and improve the computational efficiency and parallel scalability of the proposed algorithm. The solution to the Kohn-Sham problem in the localized Tucker tensor basis is computed using the Chebyshev filtered subspace iteration (ChFSI) method. Additional comments including restrictions and unusual features: The code works with Troullier-Martin (TM) pseudopotentials in Kleinman-Bylander form. The current release supports only non-periodic DFT calculations with the local density approximation (LDA) for exchange-correlation functional. This TTDFT project uses GitHub via Git, a free distributed version control software. The archived version at the time of submission of this work can be found on the CPC program library through program files DOI provided above. The GitHub repository of this project can be found on https://github.com/ttdftdev/ttdft_public. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
16. Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite system
- Author
-
Física de materiales, Materialen fisika, Tancogne Dejean, Nicolas, Oliveira, Micael J. T., Andrade, Xavier, Appel, Heiko, Borca, Carlos H., Le Breton, Guillaume, Buchholz, Florian, Castro, Alberto, Corni, Stefano, Correa, Alfredo A., De Giovannini, Umberto, Delgado, Alain, Eich, Florian G., Flick, Johannes, Gil, Gabriel, Gómez, Adrián, Helbig, Nicole, Hübener, Hannes, Jestädt, René, Jornet Somoza, Joaquim, Larsen, Ask H., Lebedeva, Irina V., Lüders, Martin, Lopes Marques, Miguel A., Ohlmann, Sebastian T., Pipolo, Silvio, Rampp, Markus, Rozzi, Carlo A., Strubbe, David A., Sato, Shunsuke A., Schäfer, Christian, Theophilou, Iris, Welden, Alicia, Rubio Secades, Angel, Física de materiales, Materialen fisika, Tancogne Dejean, Nicolas, Oliveira, Micael J. T., Andrade, Xavier, Appel, Heiko, Borca, Carlos H., Le Breton, Guillaume, Buchholz, Florian, Castro, Alberto, Corni, Stefano, Correa, Alfredo A., De Giovannini, Umberto, Delgado, Alain, Eich, Florian G., Flick, Johannes, Gil, Gabriel, Gómez, Adrián, Helbig, Nicole, Hübener, Hannes, Jestädt, René, Jornet Somoza, Joaquim, Larsen, Ask H., Lebedeva, Irina V., Lüders, Martin, Lopes Marques, Miguel A., Ohlmann, Sebastian T., Pipolo, Silvio, Rampp, Markus, Rozzi, Carlo A., Strubbe, David A., Sato, Shunsuke A., Schäfer, Christian, Theophilou, Iris, Welden, Alicia, and Rubio Secades, Angel
- Abstract
Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).
- Published
- 2020
17. Electrophoresis of concentrated colloidal dispersions in low-polar solvents
- Author
-
Vissers, Teun, Imhof, Arnout, Carrique, Félix, Delgado, Ángel V., and van Blaaderen, Alfons
- Subjects
- *
ELECTROPHORESIS , *DISPERSION (Chemistry) , *SOLVENTS , *POLARITY (Chemistry) , *COLLOIDS , *CONFOCAL microscopy , *ZETA potential , *DIFFUSION - Abstract
Abstract: We present a method to accurately measure the electrophoretic mobility of spherical colloids at high volume fractions in real space using confocal laser scanning microscopy (CLSM) and particle tracking. We show that for polymethylmethacrylate (PMMA) particles in a low-polar, density- and refractive-index-matched mixture of cyclohexylbromide and cis-decahydronaphthalene, the electrophoretic mobility decreases nonlinearly with increasing volume fraction. From the electrophoretic mobilities, we calculate the ζ-potential and the particle charge with and without correcting for volume fraction effects. For both cases, we find a decreasing particle charge as a function of volume fraction. This is in accordance with the fact that the charges originate from chemical equilibria that represent so-called weak association and/or dissociation reactions. Finally, as our methodology also provides data on particle self-diffusion in the presence of an electric field, we also analyze the diffusion at different volume fractions and identify a nonlinear decreasing trend for increasing volume fraction. [Copyright &y& Elsevier]
- Published
- 2011
- Full Text
- View/download PDF
18. A massively-parallel electronic-structure calculations based on real-space density functional theory
- Author
-
Iwata, Jun-Ichi, Takahashi, Daisuke, Oshiyama, Atsushi, Boku, Taisuke, Shiraishi, Kenji, Okada, Susumu, and Yabana, Kazuhiro
- Subjects
- *
ELECTRONIC structure , *DENSITY functionals , *FINITE differences , *PARALLEL computers , *COMPUTER software , *PROBABILITY measures - Abstract
Abstract: Based on the real-space finite-difference method, we have developed a first-principles density functional program that efficiently performs large-scale calculations on massively-parallel computers. In addition to efficient parallel implementation, we also implemented several computational improvements, substantially reducing the computational costs of operations such as the Gram–Schmidt procedure and subspace diagonalization. Using the program on a massively-parallel computer cluster with a theoretical peak performance of several TFLOPS, we perform electronic-structure calculations for a system consisting of over 10,000 Si atoms, and obtain a self-consistent electronic-structure in a few hundred hours. We analyze in detail the costs of the program in terms of computation and of inter-node communications to clarify the efficiency, the applicability, and the possibility for further improvements. [Copyright &y& Elsevier]
- Published
- 2010
- Full Text
- View/download PDF
19. Monte Carlo simulation of a model cuprate
- Author
-
V. A. Ulitko, A. S. Moskvin, Yu. D. Panov, and A. A. Chikov
- Subjects
History ,PSEUDOSPIN ,Uniform distribution (continuous) ,CLASSICAL APPROXIMATION ,NONLOCAL CORRELATIONS ,Monte Carlo method ,FOS: Physical sciences ,MONTE CARLO ALGORITHMS ,Space (mathematics) ,Education ,Condensed Matter - Strongly Correlated Electrons ,symbols.namesake ,82M31 ,Statistical physics ,Monte Carlo algorithm ,Physics ,MODEL DESCRIPTION ,MONTE CARLO'S SIMULATION ,Strongly Correlated Electrons (cond-mat.str-el) ,LOCAL CORRELATIONS ,PHASE SPACE METHODS ,CUPRATES ,MONTE CARLO METHODS ,Function (mathematics) ,Computational Physics (physics.comp-ph) ,PHASE DIAGRAMS ,Computer Science Applications ,WAVE FUNCTIONS ,Model description ,APPROXIMATION ALGORITHMS ,HEISENBERG SPINS ,Phase space ,symbols ,I.6.8 ,REAL-SPACE ,Hamiltonian (quantum mechanics) ,Physics - Computational Physics ,INTELLIGENT SYSTEMS ,COPPER COMPOUNDS - Abstract
We develop a classical Monte Carlo algorithm based on a quasi-classical approximation for a pseudospin S=1 Hamiltonian in real space to construct a phase diagram of a model cuprate with a high Tc. A model description takes into account both local and nonlocal correlations, Heisenberg spin-exchange interaction, correlated single-particle, and two-particle transport. We formulate a state selection algorithm for a given parameterization of the wave function in order to ensure a uniform distribution of states in the phase space. The simulation results show a qualitative agreement with the experimental phase diagrams., Comment: 9 pages, 3 figures
- Published
- 2021
- Full Text
- View/download PDF
20. Deformation-Space Method for the Design of Biplanar Transverse Gradient Coils in Open MRI Systems.
- Author
-
Minhua Zhu, Ling Xia, Feng Liu2, and Crozier, Stuart
- Subjects
- *
MAGNETIC resonance imaging , *MAGNETIC fields , *MAGNETICS , *ELECTROMAGNETIC induction , *MAGNETIC resonance , *ELECTROMAGNETIC fields - Abstract
We propose an efficient real-space algorithm for the design of biplanar transverse gradient coils for use in open magnetic resonance imaging (MRI) systems. In our method, each wire arc is represented by a closed contour (Limaçon). Using parametric equations, we deform/reshape an ensemble of closed contours in a simple manner, controllable by just a few parameters. These parameters are used to define system rearrangements in the design procedure. We use an iterative optimization procedure to adjust the control parameters in order to minimize cost functions such as gradient homogeneity and inductance. Here, we comapare the coil pattern designed by our deformation-space method with a pattern designed by the conventional stream function approach, and we discuss the merit of the new method. [ABSTRACT FROM AUTHOR]
- Published
- 2008
- Full Text
- View/download PDF
21. Real-space density-functional calculations for Si divacancies with large size supercell models
- Author
-
Iwata, J.-I., Oshiyama, A., and Shiraishi, K.
- Subjects
- *
SEMICONDUCTORS , *SILICON , *PSEUDOPOTENTIAL method , *FUNCTIONAL analysis - Abstract
Abstract: First-principles density-functional calculations of divacancies in crystalline silicon are performed with large size super-cell models which correspond to the unit cell of 64–1000 Si atoms. The calculations are performed by the newly developed real-space finite-difference pseudopotential code in parallel computation. It is found that the model size is large enough to provide a converged divacancy structures. A structure which is predicted by the early experiment and the recent cluster model calculation, is not found. Our own cluster calculation also corroborate the finding. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
22. DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization.
- Author
-
Das, Sambit, Motamarri, Phani, Subramanian, Vishal, Rogers, David M., and Gavini, Vikram
- Subjects
- *
DENSITY functional theory , *CODING theory , *CHEBYSHEV polynomials , *PSEUDOPOTENTIAL method , *ELECTROSTATIC interaction , *SUPERCOMPUTERS , *GRAPHICS processing units - Abstract
We present DFT-FE 1.0, building on DFT-FE 0.6 Motamarri et al. (2020) [28] , to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ∼ 100 , 000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our GPU acceleration, which yields ∼20× speed-up on hybrid CPU-GPU nodes of the Summit supercomputer. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80 − 140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ∼ 6 , 000 − 15 , 000 electrons using 64 − 224 nodes of the Summit supercomputer. Program Title: DFT-FE CPC Library link to program files: https://doi.org/10.17632/c5ghfc6ctn.1 Developer's repository link: https://github.com/dftfeDevelopers/dftfe Licensing provisions: LGPL v3 Programming language: C/C++ External routines/libraries: p4est (http://www.p4est.org/), deal.II (https://www.dealii.org/), BLAS (http://www.netlib.org/blas/), LAPACK (http://www.netlib.org/lapack/), ELPA (https://elpa.mpcdf.mpg.de/), ScaLAPACK (http://www.netlib.org/scalapack/), Spglib (https://atztogo.github.io/spglib/), ALGLIB (http://www.alglib.net/), LIBXC (http://www.tddft.org/programs/libxc/), PETSc (https://www.mcs.anl.gov/petsc), SLEPc (http://slepc.upv.es), NCCL (optional- https://github.com/NVIDIA/nccl). Nature of problem: Density functional theory calculations. Solution method: We employ a local real-space variational formulation of Kohn-Sham density functional theory that is applicable for both pseudopotential and all-electron calculations on periodic, semi-periodic and non-periodic geometries. Higher-order adaptive spectral finite-element basis is used to discretize the Kohn-Sham equations. Chebyshev polynomial filtered subspace iteration procedure (ChFSI) is employed to solve the nonlinear Kohn-Sham eigenvalue problem self-consistently. ChFSI in DFT-FE employs Cholesky factorization based orthonormalization, and spectrum splitting based Rayleigh-Ritz procedure in conjunction with mixed precision arithmetic. Configurational force approach is used to compute ionic forces and periodic cell stresses for geometry optimization. Additional comments including restrictions and unusual features: Exchange correlation functionals are restricted to Local Density Approximation (LDA) and Generalized Gradient Approximation (GGA), with and without spin. The pseudopotentials available are optimized norm conserving Vanderbilt (ONCV) pseudopotentials and Troullier–Martins (TM) pseudopotentials. Calculations are non-relativistic. DFT-FE handles all-electron and pseudopotential calculations in the same framework, while accommodating periodic, non-periodic and semi-periodic boundary conditions. [ABSTRACT FROM AUTHOR]
- Published
- 2022
- Full Text
- View/download PDF
23. Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems
- Author
-
René Jestädt, Alain Delgado, Christian Schäfer, Andrea Castro, Guillaume Le Breton, M. Lüders, Gabriel Gil, Hannes Hübener, Florian Buchholz, Adrián Gomez, Nicolas Tancogne-Dejean, Alfredo A. Correa, Sebastian T. Ohlmann, Micael J. T. Oliveira, Miguel A. L. Marques, Joaquim Jornet-Somoza, Carlos H. Borca, Markus Rampp, Angel Rubio, David A. Strubbe, Alicia Rae Welden, Iris Theophilou, F. G. Eich, Ask Hjorth Larsen, Carlo Andrea Rozzi, Umberto De Giovannini, Shunsuke A. Sato, Nicole Helbig, Johannes Flick, Heiko Appel, Irina V. Lebedeva, Silvio Pipolo, Stefano Corni, Xavier Andrade, European Commission, European Research Council, Simons Foundation, Department of Energy (US), German Research Foundation, University of California, Tancogne-Dejean N., Oliveira M.J.T., Andrade X., Appel H., Borca C.H., Le Breton G., Buchholz F., Castro A., Corni S., Correa A.A., De Giovannini U., Delgado A., Eich F.G., Flick J., Gil G., Gomez A., Helbig N., Hubener H., Jestadt R., Jornet-Somoza J., Larsen A.H., Lebedeva I.V., Luders M., Marques M.A.L., Ohlmann S.T., Pipolo S., Rampp M., Rozzi C.A., Strubbe D.A., Sato S.A., Schafer C., Theophilou I., Welden A., and Rubio A.
- Subjects
spectroscopy ,Photon ,electronic-structure calculations ,Computer science ,spectra ,Quantum dynamics ,molecular-dynamics ,Complex system ,General Physics and Astronomy ,FOS: Physical sciences ,010402 general chemistry ,spin ,01 natural sciences ,Settore FIS/03 - Fisica Della Materia ,Engineering ,TDDFT ,real-space ,0103 physical sciences ,octopus ,generalized gradient approximation ,Physical and Theoretical Chemistry ,density-functional theory ,Massively parallel ,Quantum ,Chemical Physics ,real time ,010304 chemical physics ,Computational Physics (physics.comp-ph) ,scientific software ,0104 chemical sciences ,total-energy calculations ,physics.comp-ph ,Physical Sciences ,Chemical Sciences ,polarizable continuum model ,State of matter ,Systems engineering ,Light driven ,Density functional theory ,Physics - Computational Physics - Abstract
Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials)., This work was supported by the European Research Council (Grant No. ERC-2015-AdG694097), the Cluster of Excellence “Advanced Imaging of Matter” (AIM), Grupos Consolidados (IT1249-19), and SFB925. The Flatiron Institute is a division of the Simons Foundation. X.A., A.W., and A.C. acknowledge that part of this work was performed under the auspices of the U.S. Department of Energy at Lawrence Livermore National Laboratory under Contract No. DE-AC52-07A27344. J.J.-S. gratefully acknowledges the funding from the European Union Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie Grant Agreement No. 795246-StrongLights. J.F. acknowledges financial support from the Deutsche Forschungsgemeinschaft (DFG Forschungsstipendium FL 997/1-1). D.A.S. acknowledges University of California, Merced start-up funding.
- Published
- 2020
- Full Text
- View/download PDF
24. SPARC: Simulation Package for Ab-initio Real-space Calculations
- Author
-
Andrew J. Medford, Benjamin M. Comer, Edmond Chow, Abhiraj Sharma, John E. Pask, Qimen Xu, Phanish Suryanarayana, and Hua Huang
- Subjects
Electronic structure ,Kohn–Sham ,Real-space ,Plane wave ,Ab initio ,Kohn–Sham equations ,FOS: Physical sciences ,Space (mathematics) ,01 natural sciences ,QA76.75-76.765 ,03 medical and health sciences ,Finite-differences ,0103 physical sciences ,Computer software ,010306 general physics ,030304 developmental biology ,Physics ,0303 health sciences ,Condensed Matter - Materials Science ,Finite difference ,Materials Science (cond-mat.mtrl-sci) ,Computational Physics (physics.comp-ph) ,Computer Science Applications ,Computational physics ,Density functional theory ,Physics - Computational Physics ,Software ,Order of magnitude - Abstract
We present SPARC: Simulation Package for Ab-initio Real-space Calculations. SPARC can perform Kohn-Sham density functional theory calculations for isolated systems such as molecules as well as extended systems such as crystals and surfaces, in both static and dynamic settings. It is straightforward to install/use and highly competitive with state-of-the-art planewave codes, demonstrating comparable performance on a small number of processors and increasing advantages as the number of processors grows. Notably, SPARC brings solution times down to a few seconds for systems with $\mathcal{O}(100-500)$ atoms on large-scale parallel computers, outperforming planewave counterparts by an order of magnitude and more., Comment: 17 pages, 3 figures
- Published
- 2020
- Full Text
- View/download PDF
25. Real-time in society by nonreal-time in network.
- Author
-
Ikegami, T.
- Abstract
The last century started with quantum mechanics, and at the end of the century we are enjoying our life by fruits originated and cultivated by quantum electronics. The technological impart on telecommunication has been great and marvellous and presently, in our perception about information, the meaning of distance, far or near, has been almost removed. The cyber-space and the real-space encounter each other by click and we feel the pace and the speed of our life is passing in real-time. Agility is the bottom line in businesses and alliance becomes nothing special by the Internet. Under the circumstance, the bit-greedy society suddenly comes out and optical communication businesses are dramatically expanding day by day beyond our forecast. Let us go back to ten years ago. We had cultivated the fruits of quantum electronics in laboratories and succeeded in ultra long-distance broad-band transmission experiments. The question in front of the demonstration was “who would use this technology?” We answered, “We know it might be too much for telephone service. Please wait for the video era.” “When?” and we were struck for a word. Back again now, we can put the different question, “Why does the bit-greedy society appear now without video services?” The paper addresses a part of the answer. On the contrary to our perception of real-time life, the telecommunication network has changed from “real-time” to “nonreal-time” in the procedure, in which messages are endlessly and infinitely stored in the network. This special feature of the Internet forces us to push into the bit-greedy world without video services, or only with data services [ABSTRACT FROM PUBLISHER]
- Published
- 2000
- Full Text
- View/download PDF
26. Quantitative 3D real-space studies of arrested colloidal structures and processes
- Author
-
van der Wee, Ernest Benjamin, Afd Soft Condensed Matter and Biophysics, Soft Condensed Matter and Biophysics, and van Blaaderen, Alfons
- Subjects
dissertation ,Confocal microscopy ,stimulated emission depletion microscopy (STED) ,homogeneous nucleation ,electron microscopy ,supraparticles ,Laves phases ,proefschrift ,real-space ,Colloids ,FIB-SEM tomography - Abstract
In this thesis colloidal dispersions were quantitatively studied by means of both light and electron microscopy. These colloidal particles, typically with a size somewhere between 1 and 1000 nm, have proven to be a good model system for physical phenomena in condensed matter physics, such as crystallization and the glass transition. As the colloidal particles are much larger than single atoms or molecules, they are also much slower, making it possible to study the phenomena in detail. Moreover, the larger size of the particles allows to characterize them with microscopy. We developed a method to effectively arrest colloidal dispersions for 3D study using light microscopy. We show that this method is applicable to a wide variety of colloidal dispersions. In addition, this method can also be used to preserve structures induced by an external field, such as an electric field, which would otherwise collapse when the external field is removed. This enables the development of new exciting materials. We also demonstrate that with the arrest method the nucleation of crystals of colloids can be studied in great detail. While previous microscopy studies were limited in their sample size, we show that we can now study volumes containing up to one hundred times more particles. In addition, we studied spherical assemblies of mixtures of binary colloids which form one of the Laves phases in bulk. We found that depending on the number ratio between the large and small colloids, the spherical geometry of the assembly plays no role, or surprisingly, leads to the formation of an icosahedral quasicrystal.
- Published
- 2019
27. Quantitative 3D real-space studies of arrested colloidal structures and processes
- Subjects
dissertation ,Confocal microscopy ,stimulated emission depletion microscopy (STED) ,homogeneous nucleation ,electron microscopy ,supraparticles ,Laves phases ,proefschrift ,real-space ,digestive, oral, and skin physiology ,Colloids ,FIB-SEM tomography - Abstract
In this thesis colloidal dispersions were quantitatively studied by means of both light and electron microscopy. These colloidal particles, typically with a size somewhere between 1 and 1000 nm, have proven to be a good model system for physical phenomena in condensed matter physics, such as crystallization and the glass transition. As the colloidal particles are much larger than single atoms or molecules, they are also much slower, making it possible to study the phenomena in detail. Moreover, the larger size of the particles allows to characterize them with microscopy. We developed a method to effectively arrest colloidal dispersions for 3D study using light microscopy. We show that this method is applicable to a wide variety of colloidal dispersions. In addition, this method can also be used to preserve structures induced by an external field, such as an electric field, which would otherwise collapse when the external field is removed. This enables the development of new exciting materials. We also demonstrate that with the arrest method the nucleation of crystals of colloids can be studied in great detail. While previous microscopy studies were limited in their sample size, we show that we can now study volumes containing up to one hundred times more particles. In addition, we studied spherical assemblies of mixtures of binary colloids which form one of the Laves phases in bulk. We found that depending on the number ratio between the large and small colloids, the spherical geometry of the assembly plays no role, or surprisingly, leads to the formation of an icosahedral quasicrystal.
- Published
- 2019
28. (e,2e) impact ionization processes for surface science
- Author
-
Didier Sébilleau, Junqing Xu, Calogero R. Natoli, Rakesh Choubisa, Institut de Physique de Rennes (IPR), Université de Rennes (UR)-Centre National de la Recherche Scientifique (CNRS), Université de Rennes 1 (UR1), Université de Rennes (UNIV-RENNES)-Université de Rennes (UNIV-RENNES)-Centre National de la Recherche Scientifique (CNRS), Institut de Physique de Rennes ( IPR ), Université de Rennes 1 ( UR1 ), and Université de Rennes ( UNIV-RENNES ) -Université de Rennes ( UNIV-RENNES ) -Centre National de la Recherche Scientifique ( CNRS )
- Subjects
Physics ,Surface (mathematics) ,[PHYS]Physics [physics] ,Electron-electron interactions ,[ PHYS ] Physics [physics] ,Scattering ,Real-space ,Reflection geometry ,Ionization process ,Electron ,Space (mathematics) ,Kinetic energy ,01 natural sciences ,010305 fluids & plasmas ,Impact ionization ,Surface science ,0103 physical sciences ,Surface structure ,Core-level electrons ,Atomic physics ,010306 general physics ,Spectroscopy - Abstract
International audience; We present a scattering theoretic approach to the calculation of the cross-section of (e,2e) impact spectroscopy where all the electrons involved are treated within the real space multiple scattering framework. This approach is particularly suited to the reflection geometry at low kinetic energies, with the ejection of a core-level electron. In this case, we expect (e,2e) spectroscopy can be turned into an extremely sensitive surface structure probe. © 2018, Springer International Publishing AG, part of Springer Nature.
- Published
- 2018
- Full Text
- View/download PDF
29. Maximum probability domains for Hubbard models
- Author
-
Patrick Bultinck, Dimitri Van Neck, Ward Poelmans, Guillaume Acke, Stijn De Baerdemacker, Mario Van Raemdonck, and Pieter W. Claeys
- Subjects
Hubbard model ,Biophysics ,Hubbard benzene ,FOS: Physical sciences ,010402 general chemistry ,01 natural sciences ,Projection (linear algebra) ,5-hexatriene ,Condensed Matter - Strongly Correlated Electrons ,MOLECULAR WAVE-FUNCTIONS ,PAIR ,Physics - Chemical Physics ,0103 physical sciences ,Statistical physics ,Physical and Theoretical Chemistry ,Molecular Biology ,Eigenvalues and eigenvectors ,Mathematics ,Chemical Physics (physics.chem-ph) ,Condensed Matter::Quantum Gases ,010304 chemical physics ,Strongly Correlated Electrons (cond-mat.str-el) ,Hubbard 1 ,LOCALIZATION ,CHEMICAL-BONDS ,Condensed Matter Physics ,0104 chemical sciences ,ELECTRONS ,REPRESENTATIONS ,Physics and Astronomy ,generating function ,Line (geometry) ,Maximum probability domains ,LOCALIZABILITY ,Slater determinant ,Valence bond theory ,Condensed Matter::Strongly Correlated Electrons ,REAL-SPACE ,QUANTUM-THEORY ,POPULATION ANALYSIS - Abstract
The theory of maximum probability domains (MPDs) is formulated for the Hubbard model in terms of projection operators and generating functions for both exact eigenstates as well as Slater determinants. A fast MPD analysis procedure is proposed, which is subsequently used to analyse numerical results for the Hubbard model. It is shown that the essential physics behind the considered Hubbard models can be exposed using MPDs. Furthermore, the MPDs appear to be in line with what is expected from Valence Bond (VB) Theory-based knowledge.
- Published
- 2016
30. Evaluating Hellmann–Feynman forces within non-local pseudopotentials.
- Author
-
Novák, Matyáš, Vackář, Jiří, and Cimrman, Robert
- Subjects
- *
PSEUDOPOTENTIAL method , *DENSITY functional theory , *FINITE element method - Abstract
A new approach for evaluating Hellmann–Feynman forces within a non-local potential is introduced. Particularly, the case of Hellmann–Feynman theorem applied within density functional theory in combination with nonlocal ab-initio pseudopotentials, discretized by the finite-element method, is discussed in detail. The validity of the new approach is verified using test calculations on simple molecules and the convergence properties (w.r.t. the DFT loop) are analyzed. A comparison to other previously published approaches to Hellmann–Feynman forces calculations is shown to document that the new approach mitigates, for l -dependent as well as for separable forms of nonlocal pseudopotentials, the efficiency and/or accuracy problems arising in the methods published so far. [ABSTRACT FROM AUTHOR]
- Published
- 2020
- Full Text
- View/download PDF
31. On preconditioning the self-consistent field iteration in real-space Density Functional Theory.
- Author
-
Kumar, Shashikant, Xu, Qimen, and Suryanarayana, Phanish
- Subjects
- *
DENSITY functional theory , *FUNCTION spaces , *LINEAR systems , *MEAN field theory - Abstract
• Real-space formulation for isotropic Fourier-space SCF preconditioners in DFT. • Rational function approximation of the preconditioner function in Fourier space. • Real-space application in terms of the solution of sparse Helmholtz-type systems. • Solitary linear system required for truncated-Kerker and Resta preconditioners. • Proposed approach accelerates SCF to similar extent as exact counterparts. We present a real-space formulation for isotropic Fourier-space preconditioners used to accelerate the self-consistent field iteration in Density Functional Theory calculations. Specifically, after approximating the preconditioner in Fourier space using a rational function, we express its real-space application in terms of the solution of sparse Helmholtz-type systems. Using the truncated-Kerker and Resta preconditioners as representative examples, we show that the proposed real-space method is both accurate and efficient, requiring the solution of a single linear system, while accelerating self-consistency to the same extent as its exact Fourier-space counterpart. [ABSTRACT FROM AUTHOR]
- Published
- 2020
- Full Text
- View/download PDF
32. Application of time-dependent current-density-functional theory to nonlocal exchange-correlation effects in polymers
- Author
-
Johannes Berger, Jaap G. Snijders, R. van Leeuwen, P. L. de Boeij, M. van Faassen, Zernike Institute for Advanced Materials, and Theoretical Chemistry
- Subjects
CORRELATION POTENTIALS ,NONLINEAR POLARIZABILITIES ,POLARIZATION ,Chemistry ,STATIC LONGITUDINAL POLARIZABILITY ,General Physics and Astronomy ,MOLECULAR CHAINS ,Polarization (waves) ,ELECTRIC-FIELD ,Bond length ,DIELECTRICS ,Polarizability ,SYSTEMS ,Quantum mechanics ,Density functional theory ,REAL-SPACE ,Physical and Theoretical Chemistry ,Local-density approximation ,LINEAR-RESPONSE ,Adiabatic process ,Wave function ,Current density - Abstract
We provide a successful approach towards the solution of the longstanding problem of the large overestimation of the static polarizability of conjugated oligomers obtained using the local density approximation within density-functional theory. The local approximation is unable to describe the highly nonlocal exchange and correlation effects found in these quasi-one-dimensional systems. Time-dependent current-density-functional theory enables us to describe ultranonlocal exchange-correlation effects within a local current description. Recently a brief account was given of the application of the Vignale-Kohn current-functional [G. Vignale and W. Kohn, Phys. Rev. Lett. 77, 2037 (1996)] to the axial polarizability of oligomer chains [M. van Faassen, P. L. de Boeij, R. van Leeuwen, J. A. Berger, and J. G. Snijders, Phys. Rev. Lett. 88, 186401 (2002)]. With the exception of the model hydrogen chain, our results were in excellent agreement with best available wavefunction methods. In the present work we further outline the underlying theory and describe how the Vignale-Kohn functional was implemented. We elaborate on earlier results and present new results for the oligomers of polyethylene, polysilane, polysilene, polymethineimine, and polybutatriene. The adiabatic local density approximation gave good results for polyethylene, which were slightly modified by the Vignale-Kohn functional. In all other cases the Vignale-Kohn functional gave large improvements upon the adiabatic local density approximation. The Vignale-Kohn results were in agreement with best available data from wave function methods. We further analyze the hydrogen chain model for different bond length alternations. In all these cases the Vignale-Kohn correction upon the adiabatic local density approximation was too small. Arguments are given that further improvements of the functional are needed. (C) 2003 American Institute of Physics.
- Published
- 2003
- Full Text
- View/download PDF
33. Sample Corrugation Affects the Apparent Bond Lengths in Atomic Force Microscopy
- Author
-
Boneschanscher, Mark P., Hamalainen, Sampsa K., Liljeroth, Peter, Swart, Ingmar, Sub Condensed Matter and Interfaces, Sub Practicum, Condensed Matter and Interfaces, Sub Condensed Matter and Interfaces, Sub Practicum, and Condensed Matter and Interfaces
- Subjects
Electrostatic force microscope ,corrugation ,Analytical chemistry ,SCANNING PROBE MICROSCOPY ,General Physics and Astronomy ,Atomic force acoustic microscopy ,Molecular physics ,law.invention ,Scanning probe microscopy ,law ,General Materials Science ,Spectroscopy ,SINGLE-MOLECULE ,atomic force microscopy ,SPECTROSCOPY ,IDENTIFICATION ,Graphene ,Chemistry ,graphene ,General Engineering ,Conductive atomic force microscopy ,Bond length ,GRAPHENE MOIRE ,tip relaxation ,bond length ,REAL-SPACE ,Non-contact atomic force microscopy - Abstract
Frequency modulation atomic force microscopy (AFM) allows the chemical structure of planar molecules to be determined with atomic resolution. Typically, these measurements are carried out in constant-height mode using carbon monoxide (CO) terminated tips. Such tips exhibit considerable flexibility, i.e., the CO molecule can bend laterally due to the tip-sample interaction. Using epitaxial graphene as a model system, we demonstrate experimentally that the apparent atomic positions measured by AFM depend on the sample corrugation. Using molecular mechanics simulations, we explain these observations by the interplay of the CO bending and the nonlinear background signal arising from the neighboring atoms. These effects depend nontrivially on the tip-sample distance and limit the achievable accuracy on the bond length determination based on AFM experiments.
- Published
- 2014
- Full Text
- View/download PDF
34. Large Scale Electronic Structure Studies on the Energetics of Dislocations in Al-Mg Materials System and Its Connection to Mesoscale Models
- Author
-
Das, Sambit
- Subjects
- Electronic structure, Real-space, Finite-elements, Dislocation core, Discrete dislocation dynamics
- Abstract
Computational modeling of dislocation behavior is vital for designing new lightweight metallic alloys. However, extraordinary challenges are posed by the multiscale physics ranging over a vast span of interacting length-scales from electronic-structure and atomic-scale effects at the dislocation core ($< 10^{-9} {rm m}$) to long-ranged elastic interactions at the continuum scale ($sim 10 upmu$). In particular, quantification of the energetics associated with electronic-structure effects inside the dislocation core and its interaction with the external macroscopic elastic fields have not been explored due to limitations of current electronic-structure methods based on the widely used plane-wave based discretization. This thesis seeks to address the above challenges by developing computational methodologies to conduct large-scale real-space electronic-structure studies of energetics of dislocations in Aluminum and Magnesium, and use these results to develop phenomenological connections to mesoscale models of plasticity like discrete dislocation dynamics (DDD), which study the collective behavior of the dislocations at longer length scales ($sim$ 1--15 $upmu$). First, a local real-space formulation of orbital-free Density Functional Theory is developed based on prior work, and implemented using finite-element discretization. The local real-space formulation coupled with bulk Dirichlet boundary conditions enables a direct computation of the isolated dislocation core energy. Studies on dislocations in Aluminum and Magnesium suggest that the core-size---region with significant contribution of electronic effects to dislocation energetics---is around seven to eleven times the magnitude of the Burgers vector. This is in stark contrast to prior displacement field based core size estimates of one to three times the magnitude of the Burgers vector. Interestingly, our study further indicates that the core-energy of the dislocations in both Aluminum and Magnesium is strongly dependent on external macroscopic strains with a non-zero slope at zero external strain. Next, the computed dislocation core energetics is used to develop a continuum model for an arbitrary aggregate of dislocations in an infinite isotropic elastic continua. This model, which accounts for the core energy dependence on macroscopic deformation provides a phenomenological approach to incorporate the electronic structure effects into mesoscale DDD simulations. Application of this model to derive nodal forces in a discrete dislocation network, leads to additional configurational forces beyond those considered in existing DDD models. Using case studies, we show that even up to distances of $10-15$ nm between the dislocations, these additional configurational forces are non-trivial in relation to the elastic Peach-Koehler force. Furthermore, the core force model is incorporated into a DDD implementation, where significant influence of core force on elementary dislocation mechanisms in Aluminum such as critical stress of a Frank-Read source and structure of a dislocation binary junction are demonstrated. To enable the above electronic-structure studies of dislocations in generic material systems, calculations using the more accurate and transferable Kohn-Sham Density Functional Theory (KS-DFT) are required. The final part of this thesis extends previous work on real-space adaptive spectral finite-element discretization of KS-DFT to develop numerical strategies and implementation innovations, which significantly reduce the computational pre-factor, while increasing the arithmetic intensity and lowering the data movement costs on both many-core and heterogeneous architectures. This has enabled systematically convergent and massively parallel (demonstrated up to 192,000 MPI tasks) KS-DFT calculations on material systems up to $sim 100,000$ electrons. Using GPUs, an unprecedented sustained performance of 46 PFLOPS (27.8% peak FP64 performance) is demonstrated on a large-scale benchmark dislocation system in Magnesium containing 105,080 electrons.
- Published
- 2019
35. DFT-FE – A massively parallel adaptive finite-element code for large-scale density functional theory calculations.
- Author
-
Motamarri, Phani, Das, Sambit, Rudraraju, Shiva, Ghosh, Krishnendu, Davydov, Denis, and Gavini, Vikram
- Subjects
- *
DENSITY functional theory , *AB-initio calculations , *FUNCTIONALS , *SPECTRAL element method , *CHEBYSHEV polynomials , *PROGRAMMING languages , *INTERIOR-point methods - Abstract
We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching ∼ 100 , 000 electrons) using Kohn–Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn–Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn–Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes—both in CPU-times and wall-times, and on both non-periodic and periodic systems—at systems sizes beyond a few thousand electrons, with over 5 − 10 fold speedups in systems with more than 10,000 electrons. The benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated up to 192,000 MPI tasks. Program Title: DFT-FE Program Files doi: http://dx.doi.org/10.17632/tgdmgvmfft.1 Licensing provisions: LGPL v3 Programming language: C/C++ External routines/libraries: p4est (http://www.p4est.org/), deal.II (https://www.dealii.org/), BLAS (http://www.netlib.org/blas/), LAPACK (http://www.netlib.org/lapack/), ELPA (https://elpa.mpcdf.mpg.de/), ScaLAPACK (http://www.netlib.org/scalapack/), Spglib (https://atztogo.github.io/spglib/), ALGLIB (http://www.alglib.net/), LIBXC (http://www.tddft.org/programs/libxc/), PETSc (https://www.mcs.anl.gov/petsc), SLEPc (http://slepc.upv.es). Nature of problem: Density functional theory calculations. Solution method: We employ a local real-space variational formulation of Kohn–Sham density functional theory that is applicable for both pseudopotential and all-electron calculations with arbitrary boundary conditions. Higher-order adaptive spectral finite-element basis is used to discretize the Kohn–Sham equations. Chebyshev polynomial filtered subspace iteration procedure (ChFSI) is employed to solve the nonlinear Kohn–Sham eigenvalue problem self-consistently. ChFSI in DFT-FE employs Cholesky factorization based orthonormalization, and spectrum splitting based Rayleigh–Ritz procedure in conjunction with mixed precision arithmetic. Configurational force approach is used to compute ionic forces and periodic cell stresses for geometry optimization. Restrictions: Exchange–correlation functionals are restricted to Local Density Approximation (LDA) and Generalized Gradient Approximation (GGA), with and without spin. The pseudopotentials available are optimized norm conserving Vanderbilt (ONCV) pseudopotentials and Troullier–Martins (TM) pseudopotentials. Calculations are non-relativistic. Unusual features: DFT-FE handles all-electron and pseudopotential calculations in the same framework, while accommodating periodic, non-periodic and semi-periodic boundary conditions. Additional comments: DFT-FE github https://github.com/dftfeDevelopers/dftfe [ABSTRACT FROM AUTHOR]
- Published
- 2020
- Full Text
- View/download PDF
36. Electronic and magnetic properties of bulk Cr tips for scanning tunneling spectroscopy
- Author
-
Guido Fratesi, Matteo Passoni, Lixin Ning, Alberto Brambilla, Carlo Spartaco Casari, Mario Italo Trioni, A. Li Bassi, Fabio Donati, Donati, F, Fratesi, G, Ning, L, Brambilla, A, Trioni, M, Bassi, A, Casari, C, and Passoni, M
- Subjects
density ,Materials science ,Condensed matter physics ,Scanning tunneling spectroscopy ,Condensed Matter Physics ,AU(111) surface ,Electronic, Optical and Magnetic Materials ,law.invention ,generalized-gradient-approximation ,law ,real-space ,Condensed Matter::Superconductivity ,antiferromagnetism ,Microscopy ,augmented-wave method ,microscopy ,Antiferromagnetism ,Density functional theory ,surface-state ,chromium ,Scanning tunneling microscope ,Quantum tunnelling ,Surface states ,Electronic properties - Abstract
The electronic and magnetic properties of bulk Cr tips for scanning tunneling microscopy have been investigated by means of density functional theory (DFT) calculations and scanning tunneling spectroscopy (STS) measurements. Spin-resolved densities of states (DOS) were calculated for model tips, i.e., Cr adatoms and clusters on ideal Cr surfaces. STS measurements on Au(111) and Si(111)-7×7 have been interpreted by modeling the tunneling process in a Wentzel-Kramers-Brillouin approximation in order to ascertain the role of tip DOS features. Calculated spin-resolved electronic properties of Cr tips have been used to revisit spin-polarized (SP) STS measurements from the literature. The agreement between experimental findings and DFT calculations confirms the relevance of the knowledge of tip electronic properties and of an appropriate model for the tunneling process in the interpretation of STS and SP-STS data.
- Published
- 2013
- Full Text
- View/download PDF
37. Packing confined hard spheres denser with adaptive prism phases
- Author
-
Francisco Meseguer, Matthieu Marechal, I. Rodriguez, Fernando Ramiro-Manzano, Hartmut Löwen, René Messina, and Erdal C. Oğuz
- Subjects
Hard spheres ,Materials science ,Monte Carlo method ,Real-space ,General Physics and Astronomy ,FOS: Physical sciences ,Geometry ,02 engineering and technology ,Condensed Matter - Soft Condensed Matter ,01 natural sciences ,Square (algebra) ,Optics ,Close packing ,0103 physical sciences ,Freezing ,Mesoscale and Nanoscale Physics (cond-mat.mes-hall) ,Computer Simulation ,Particle Size ,Internal structure ,Prismatic structures ,010306 general physics ,Plate separation ,Condensed Matter - Statistical Mechanics ,Range (particle radiation) ,Statistical Mechanics (cond-mat.stat-mech) ,Condensed Matter - Mesoscale and Nanoscale Physics ,Stabilized colloids ,business.industry ,Close-packing of equal spheres ,Monte Carlo Simulation ,Models, Theoretical ,021001 nanoscience & nanotechnology ,Soft Condensed Matter (cond-mat.soft) ,SPHERES ,Prism ,0210 nano-technology ,business ,Crystallization - Abstract
[EN] We show that hard spheres confined between two parallel hard plates pack denser with periodic adaptive prismatic structures which are composed of alternating prisms of spheres. The internal structure of the prisms adapts to the slit height which results in close packings for a range of plate separations, just above the distance where three intersecting square layers fit exactly between the plates. The adaptive prism phases are also observed in real-space experiments on confined sterically stabilized colloids and in Monte Carlo simulations at finite pressure. © 2012 American Physical Society., We thank Elvira Bonet, Moises Garin, and Kevin Mutch for helpful discussions. This work was partially supported by the DFG within the SFB TR6 (Project D1), and by the Spanish CICyT Projects FIS2009-07812 and PROMETEO/2010/043. F. R.-M. acknowledges the support from the EU Marie Curie Project APPCOPTOR-275150 (FP7-PEOPLE-2010-IEF).
- Published
- 2012
- Full Text
- View/download PDF
38. The zCOSMOS Survey. The dependence of clustering on luminosity and stellar mass at z=0.2-1
- Author
-
C. Knobel, Graziano Coppa, Anton M. Koekemoer, M. Fumana, E. Zucca, Alvio Renzini, Lidia Tasca, A. Iovino, J. F. Le Borgne, Cristiano Porciani, M. Scodeggio, B. Meneux, P. Memeo, John D. Silverman, Luigi Guzzo, U. Abbas, Y. Peng, C. M. Carollo, Daniela Vergani, Pascal Oesch, Roberto Scaramella, V. Le Brun, V. Mainieri, Simon J. Lilly, E. Ricciardelli, Christian Maier, Thierry Contini, Christian Marinoni, Roser Pello, Olga Cucciati, L. de Ravel, Laurence Tresse, G. Zamorani, K. Kovac, Karina Caputi, P. Kampczyk, H. J. McCracken, J. P. Kneib, A. Cappi, Micol Bolzonella, Alexie Leauthaud, S. Bardelli, Bianca Garilli, Angela Bongiorno, E. Perez Montero, Alessandro Cimatti, Masayuki Tanaka, S. de la Torre, P. Franzetti, O. Le Fèvre, D. Maccagni, F. Lamareille, Paolo Cassata, Lucia Pozzetti, D. Bottini, INAF - Osservatorio Astronomico di Bologna (OABO), Istituto Nazionale di Astrofisica (INAF), Max-Planck-Institut für Extraterrestrische Physik (MPE), Universitats-Sternwarte [München], Ludwig-Maximilians-Universität München (LMU), INAF - Osservatorio Astronomico di Brera (OAB), Laboratoire d'Astrophysique de Marseille (LAM), Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Aix Marseille Université (AMU)-Centre National d'Études Spatiales [Toulouse] (CNES), AUTRES, INAF- Milano, University of Chicago, Laboratoire Astrophysique de Toulouse-Tarbes (LATT), Université Toulouse III - Paul Sabatier (UT3), Université Fédérale Toulouse Midi-Pyrénées-Université Fédérale Toulouse Midi-Pyrénées-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire Midi-Pyrénées (OMP), Météo France-Centre National d'Études Spatiales [Toulouse] (CNES)-Université Fédérale Toulouse Midi-Pyrénées-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD)-Météo France-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD)-Centre National de la Recherche Scientifique (CNRS), Laboratoire d'Astrophysique de l'Observatoire Midi-Pyrénées (LATT), Institute of Astronomy [ETH Zürich], Department of Physics [ETH Zürich] (D-PHYS), Eidgenössische Technische Hochschule - Swiss Federal Institute of Technology [Zürich] (ETH Zürich)- Eidgenössische Technische Hochschule - Swiss Federal Institute of Technology [Zürich] (ETH Zürich), Human Genetics Center, The University of Texas Health Science Center at Houston (UTHealth), School of Biotechnology and Biomolecular Sciences, University of New South Wales [Sydney] (UNSW), Space Telescope Science Institute (STSci), Centre de Physique Théorique - UMR 6207 (CPT), Centre National de la Recherche Scientifique (CNRS)-Université de Toulon (UTLN)-Université de Provence - Aix-Marseille 1-Université de la Méditerranée - Aix-Marseille 2, Centre de Physique Théorique - UMR 7332 (CPT), Aix Marseille Université (AMU)-Université de Toulon (UTLN)-Centre National de la Recherche Scientifique (CNRS), Institut d'Astrophysique de Paris (IAP), Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Pierre et Marie Curie - Paris 6 (UPMC), Istituto di Astrofisica Spaziale e Fisica cosmica - Bologna (IASF-Bo), Meneux B., Guzzo L., de La Torre S., Porciani C., Zamorani G., Abbas U., Bolzonella M., Garilli B., Iovino A., Pozzetti L., Zucca E., Lilly S. J., Le Fèvre O., Kneib J.-P., Carollo C. M., Contini T., Mainieri V., Renzini A., Scodeggio M., Bardelli S., Bongiorno A., Caputi K., Coppa G., Cucciati O., de Ravel L., Franzetti P., Kampczyk P., Knobel C., Kovač, K., Lamareille F., Le Borgne J.-F., Le Brun V., Maier C., Pellò R., Peng Y., Perez Montero E., Ricciardelli E., Silverman J. D., Tanaka M., Tasca L., Tresse L., Vergani D., Bottini D., Cappi A., Cimatti A., Cassata P., Fumana M., Koekemoer A. M., Leauthaud A., Maccagni D., Marinoni C., McCracken H. J., Memeo P., Oesch P., Scaramella R., Astronomy, Aix Marseille Université (AMU)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS), Université de Toulouse (UT)-Université de Toulouse (UT)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire Midi-Pyrénées (OMP), Institut de Recherche pour le Développement (IRD)-Université Toulouse III - Paul Sabatier (UT3), Université de Toulouse (UT)-Université de Toulouse (UT)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Météo-France -Institut de Recherche pour le Développement (IRD)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National d'Études Spatiales [Toulouse] (CNES)-Centre National de la Recherche Scientifique (CNRS)-Météo-France -Centre National de la Recherche Scientifique (CNRS), Université de la Méditerranée - Aix-Marseille 2-Université de Provence - Aix-Marseille 1-Université de Toulon (UTLN)-Centre National de la Recherche Scientifique (CNRS), and Université Pierre et Marie Curie - Paris 6 (UPMC)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS)
- Subjects
Cosmology and Nongalactic Astrophysics (astro-ph.CO) ,Stellar mass ,Large-scale structure of Universe ,media_common.quotation_subject ,SPACE CORRELATION-FUNCTIONS ,FOS: Physical sciences ,Astrophysics ,Astrophysics::Cosmology and Extragalactic Astrophysics ,Surveys ,01 natural sciences ,Luminosity ,[PHYS.ASTR.CO]Physics [physics]/Astrophysics [astro-ph]/Cosmology and Extra-Galactic Astrophysics [astro-ph.CO] ,GALAXY REDSHIFT SURVEY ,1ST EPOCH DATA ,0103 physical sciences ,LARGE-SCALE STRUCTURE ,Cluster analysis ,010303 astronomy & astrophysics ,PHOTOMETRIC REDSHIFTS ,Astrophysics::Galaxy Astrophysics ,media_common ,Physics ,Galaxy: evolution ,010308 nuclear & particles physics ,Cosmology: observations ,Astronomy and Astrophysics ,VLT DEEP SURVEY ,Covariance ,ENVIRONMENTAL DEPENDENCE ,Astrophysics - Astrophysics of Galaxies ,Galaxy ,Redshift ,Universe ,HALO OCCUPATION DISTRIBUTION ,Amplitude ,Space and Planetary Science ,Astrophysics of Galaxies (astro-ph.GA) ,DIGITAL SKY SURVEY ,REAL-SPACE ,Astrophysics - Cosmology and Nongalactic Astrophysics - Abstract
We study the dependence of galaxy clustering on luminosity and stellar mass at redshifts z ~ [0.2-1] using the first zCOSMOS 10K sample. We measure the redshift-space correlation functions xi(rp,pi) and its projection wp(rp) for sub-samples covering different luminosity, mass and redshift ranges. We quantify in detail the observational selection biases and we check our covariance and error estimate techniques using ensembles of semi-analytic mock catalogues. We finally compare our measurements to the cosmological model predictions from the mock surveys. At odds with other measurements, we find a weak dependence of galaxy clustering on luminosity in all redshift bins explored. A mild dependence on stellar mass is instead observed. At z~0.7, wp(rp) shows strong excess power on large scales. We interpret this as produced by large-scale structure dominating the survey volume and extending preferentially in direction perpendicular to the line-of-sight. We do not see any significant evolution with redshift of the amplitude of clustering for bright and/or massive galaxies. The clustering measured in the zCOSMOS data at 0.5=10 is only marginally consistent with predictions from the mock surveys. On scales larger than ~2 h^-1 Mpc, the observed clustering amplitude is compatible only with ~1% of the mocks. Thus, if the power spectrum of matter is LCDM with standard normalization and the bias has no unnatural scale-dependence, this result indicates that COSMOS has picked up a particularly rare, ~2-3 sigma positive fluctuation in a volume of ~10^6 h^-1 Mpc^3. These findings underline the need for larger surveys of the z~1 Universe to appropriately characterize the level of structure at this epoch., Comment: 18 pages, 21 figures, accepted for publication in Astronomy and Astrophysics
- Published
- 2009
- Full Text
- View/download PDF
39. Electronic and magnetic properties of bulk Cr tips for scanning tunneling spectroscopy
- Author
-
Donati, F, Fratesi, G, Ning, L, Brambilla, A, Trioni, M, Bassi, A, Casari, C, Passoni, M, FRATESI, GUIDO, Passoni, M., Donati, F, Fratesi, G, Ning, L, Brambilla, A, Trioni, M, Bassi, A, Casari, C, Passoni, M, FRATESI, GUIDO, and Passoni, M.
- Abstract
The electronic and magnetic properties of bulk Cr tips for scanning tunneling microscopy have been investigated by means of density functional theory (DFT) calculations and scanning tunneling spectroscopy (STS) measurements. Spin-resolved densities of states (DOS) were calculated for model tips, i.e., Cr adatoms and clusters on ideal Cr surfaces. STS measurements on Au(111) and Si(111)-7x7 have been interpreted by modeling the tunneling process in a Wentzel-Kramers-Brillouin approximation in order to ascertain the role of tip DOS features. Calculated spin-resolved electronic properties of Cr tips have been used to revisit spin-polarized (SP) STS measurements from the literature. The agreement between experimental findings and DFT calculations confirms the relevance of the knowledge of tip electronic properties and of an appropriate model for the tunneling process in the interpretation of STS and SP-STS data.
- Published
- 2013
40. Packing confined hard spheres denser with adaptive prism phases
- Author
-
Universitat Politècnica de València. Centro de Tecnologías Físicas: Acústica, Materiales y Astrofísica - Centre de Tecnologies Físiques: Acústica, Materials i Astrofísica, European Commission, Generalitat Valenciana, Deutsche Forschungsgemeinschaft, Ministerio de Ciencia e Innovación, Oguz, Erdal C., Merechal, Matthieu, Ramiro Manzano, Fernando, Rodríguez, Marie-Isabelle, Messina, Rene, Meseguer Rico, Francisco Javier, Loewen, Hartmut, Universitat Politècnica de València. Centro de Tecnologías Físicas: Acústica, Materiales y Astrofísica - Centre de Tecnologies Físiques: Acústica, Materials i Astrofísica, European Commission, Generalitat Valenciana, Deutsche Forschungsgemeinschaft, Ministerio de Ciencia e Innovación, Oguz, Erdal C., Merechal, Matthieu, Ramiro Manzano, Fernando, Rodríguez, Marie-Isabelle, Messina, Rene, Meseguer Rico, Francisco Javier, and Loewen, Hartmut
- Abstract
[EN] We show that hard spheres confined between two parallel hard plates pack denser with periodic adaptive prismatic structures which are composed of alternating prisms of spheres. The internal structure of the prisms adapts to the slit height which results in close packings for a range of plate separations, just above the distance where three intersecting square layers fit exactly between the plates. The adaptive prism phases are also observed in real-space experiments on confined sterically stabilized colloids and in Monte Carlo simulations at finite pressure. © 2012 American Physical Society.
- Published
- 2012
41. On complex solutions of the eikonal equation
- Author
-
E. Hasanov, Işık Üniversitesi, Fen Edebiyat Fakültesi, Matematik Bölümü, Işık University, Faculty of Arts and Sciences, Department of Mathematics, and Hasanoğlu, Elman
- Subjects
Differential equations ,New approaches ,Advanced applications ,Equations ,Differential equation ,Astrophysics::High Energy Astrophysical Phenomena ,Real-space ,Refractive index ,Geometry ,Rays ,Space (mathematics) ,Gaussian beams ,Transient analysis ,Electromagnetism ,Inhomogeneous media ,Minkowski space ,Inhomogeneous medium ,Physics ,Geometical optics ,Wave propagation ,Geometrical optics ,Eikonal equation ,Magnetism ,Electro magnetics ,Complex solutions ,Ketones ,Eikonal equations ,Eikonal approximation ,Minkowski geometries ,Classical mechanics ,Optical variables control ,Optical sensors ,Characteristic methods ,Optical refraction ,Extraterrestrial measurements ,Complex Rays ,International conferences ,Mathematics ,Gaussian optics - Abstract
In this paper a new approach of complex rays in an inhomogeneous medium is presented. Complex rays are complex solutions of the eikonal equation, the main equation of the geometical optics. It is shown that solving the eikonal equation by using the characteristic method naturally leads to the pseudoriemann and Minkowski geometries. In framework of these geometries complex rays , like the real ones, can be drawn in real space and they may have caustics, and caustics also can be drawn in real space. Publisher's Version
- Published
- 2007
42. Time-dependent current-density-functional theory for the metallic response of solids
- Author
-
P. L. de Boeij, Pina Romaniello, Theoretical Chemistry, and Zernike Institute for Advanced Materials
- Subjects
NOBLE-METALS ,Physics ,ELECTRON-ELECTRON SCATTERING ,III-V-SEMICONDUCTORS ,Condensed matter physics ,BRILLOUIN-ZONE ,Fermi surface ,NONMETALLIC CRYSTALS ,Dielectric ,Condensed Matter Physics ,DIELECTRIC-CONSTANT ,Drude model ,OPTICAL 2ND-HARMONIC GENERATION ,Electronic, Optical and Magnetic Materials ,Adiabatic theorem ,LOCAL-FIELD ,Condensed Matter::Strongly Correlated Electrons ,Density functional theory ,QUADRATIC INTEGRATION ,REAL-SPACE ,Local-density approximation ,Adiabatic process ,Local field - Abstract
We extend the formulation of time-dependent current-density-functional theory for the linear response properties of dielectric and semi-metallic solids [Kootstra , J. Chem. Phys. 112, 6517 (2000)] to treat metals as well. To achieve this, the Kohn-Sham response functions have to include both interband and intraband transitions with an accurate treatment of the Fermi surface in the Brillouin-zone integrations. The intraband contributions in particular have to be evaluated using a wave-vector-dependent description. To test the method we calculate the optical properties of the two noble metals Cu and Ag. The dielectric and energy loss functions are compared with experiments and with the classical Drude theory. In general we find a good agreement with the experiments for the calculated results obtained within the adiabatic local density approximation. In order to describe the Drude-like absorption below the interband onset and the sharp plasma feature in silver exchange-correlation, effects beyond the adiabatic local density approximation are needed, which may be included in a natural way in the present current-density-functional approach.
- Published
- 2005
- Full Text
- View/download PDF
43. The Keldysh formalism applied to time-dependent current-density-functional theory
- Author
-
Van Leeuwen, Robert, Gidopoulos, NI, Wilson, S, and Theoretical Chemistry
- Subjects
NONEQUILIBRIUM PROCESSES ,DIELECTRIC FUNCTION ,GAS ,Physics::Atomic and Molecular Clusters ,NONMETALLIC CRYSTALS ,REAL-SPACE ,FIELD ,ELECTRON-SYSTEMS - Abstract
In this work we demonstrate how to derive the Kohn-Sham equations of time-dependent current-density functional theory from a generating action functional defined on a Keldysh time contour. These Kohn-Sham equations contain an exchange-correlation contribution to the vector potential. For this quantity we derive an integral equation. We further derive an integral equation for its functional derivative, the exchange-correlation kernel, which plays an essential role in response theory. The exchange-only limits of the latter equation is studied in detail for the electron gas and future applications are discussed.
- Published
- 2003
44. The Keldysh formalism applied to time-dependent current-density-functional theory
- Subjects
NONEQUILIBRIUM PROCESSES ,DIELECTRIC FUNCTION ,GAS ,NONMETALLIC CRYSTALS ,REAL-SPACE ,FIELD ,ELECTRON-SYSTEMS - Abstract
In this work we demonstrate how to derive the Kohn-Sham equations of time-dependent current-density functional theory from a generating action functional defined on a Keldysh time contour. These Kohn-Sham equations contain an exchange-correlation contribution to the vector potential. For this quantity we derive an integral equation. We further derive an integral equation for its functional derivative, the exchange-correlation kernel, which plays an essential role in response theory. The exchange-only limits of the latter equation is studied in detail for the electron gas and future applications are discussed.
- Published
- 2003
45. The Keldysh formalism applied to time-dependent current-density-functional theory
- Subjects
NONEQUILIBRIUM PROCESSES ,DIELECTRIC FUNCTION ,GAS ,Physics::Atomic and Molecular Clusters ,NONMETALLIC CRYSTALS ,REAL-SPACE ,FIELD ,ELECTRON-SYSTEMS - Abstract
In this work we demonstrate how to derive the Kohn-Sham equations of time-dependent current-density functional theory from a generating action functional defined on a Keldysh time contour. These Kohn-Sham equations contain an exchange-correlation contribution to the vector potential. For this quantity we derive an integral equation. We further derive an integral equation for its functional derivative, the exchange-correlation kernel, which plays an essential role in response theory. The exchange-only limits of the latter equation is studied in detail for the electron gas and future applications are discussed.
- Published
- 2003
46. The 2dF Galaxy Redshift Survey: The dependence of galaxy clustering on luminosity and spectral type
- Author
-
Carlton M. Baugh, Matthew Colless, Warrick J. Couch, Ed Hawkins, Chris A. Collins, Carole Jackson, Shaun Cole, Ivan K. Baldry, Ian Lewis, John A. Peacock, Richard S. Ellis, Keith Taylor, Gavin Dalton, Peder Norberg, T. Bridges, Steve Maddox, Bruce A. Peterson, dFGRS Team, R. D. Cannon, Joss Bland-Hawthorn, George Efstathiou, O. Lahav, Stuart Lumsden, Darren Madgwick, Carlos S. Frenk, Simon P. Driver, William J. Sutherland, and Karl Glazebrook
- Subjects
Galaxies formation ,Statistical methods ,Evolution ,Population ,Catalog ,Real-space ,Velocity ,Large-scale structure of universe ,FOS: Physical sciences ,Formation ,Astrophysics ,Astrophysics::Cosmology and Extragalactic Astrophysics ,Correlation function (astronomy) ,Stellar classification ,Luminosity ,Emission ,Bias ,Large-scale structure of Universe ,Cluster analysis ,education ,Space distortions ,Astrophysics::Galaxy Astrophysics ,2dF Galaxy Redshift Survey ,Physics ,education.field_of_study ,Astrophysics (astro-ph) ,Astronomy and Astrophysics ,Statistical ,Galaxies ,Galaxy ,Amplitude ,Space and Planetary Science ,Numerical methods ,Morphological segregation - Abstract
We investigate the dependence of galaxy clustering on luminosity and spectral type using the 2dF Galaxy Redshift Survey (2dFGRS). Spectral types are assigned using the principal component analysis of Madgwick et al. We divide the sample into two broad spectral classes: galaxies with strong emission lines (`late-types'), and more quiescent galaxies (`early-types'). We measure the clustering in real space, free from any distortion of the clustering pattern due to peculiar velocities, for a series of volume-limited samples. The projected correlation functions of both spectral types are well described by a power law for transverse separations in the range 2 < (sigma/Mpc/h) < 15, with a marginally steeper slope for early-types than late-types. Both early and late types have approximately the same dependence of clustering strength on luminosity, with the clustering amplitude increasing by a factor of ~2.5 between L* and 4 L*. At all luminosities, however, the correlation function amplitude for the early-types is ~50% higher than that of the late-types. These results support the view that luminosity, and not type, is the dominant factor in determining how the clustering strength of the whole galaxy population varies with luminosity., accepted by MNRAS after minor revision. 13 pages, 10 figures
- Published
- 2001
47. Computing accurate solutions to the Kohn-Sham problem quickly in real space
- Author
-
Schofield, Grady Lynn
- Subjects
- Hermitian eigenproblem, Kohn-Sham equation, Density functional theory, Quantum forces, Spectrum slicing, Real-space, Pseudopotential
- Abstract
Matter on a length scale comparable to that of a chemical bond is governed by the theory of quantum mechanics, but quantum mechanics is a many body theory, hence for the sake of chemistry or solid state physics, finding solutions to the governing equation, Schrodinger's equation, is hopeless for all but the smallest of systems. As the number of electrons increases, the complexity of solving the equations grows rapidly without bound. One way to make progress is to treat the electrons in a system as independent particles and to attempt to capture the many-body effects in a functional of the electrons' density distribution. When this approximation is made, the resulting equation is called the Kohn-Sham equation, and instead of requiring solving for one function of many variables, it requires solving for many functions of the three spatial variables. This problem turns out to be easier than the many body problem, but it still scales cubically in the number of electrons. In this work we will explore ways of obtaining the solutions to the Kohn-Sham equation in the framework of real-space pseudopotential density functional theory. The Kohn-Sham equation itself is an eigenvalue problem, just as Schrodinger's equation. For each electron in the system, there is a corresponding eigenvector. So the task of solving the equation is to compute many eigenpairs of a large Hermitian matrix. In order to mitigate the problem of cubic scaling, we develop an algorithm to slice the spectrum into disjoint segments. This allows a smaller eigenproblem to be solved in each segment where a post-processing step combines the results from each segment and prevents double counting of the eigenpairs. The efficacy of this method depends on the use of high order polynomial filters that enhance only a segment of the spectrum. The order of the filter is the number of matrix-vector multiplication operations that must be done with the Hamiltonian. Therefore the performance of these operations is critical. We develop a scalable algorithm for computing these multiplications and introduce a new density functional theory code implementing the algorithm.
- Published
- 2014
Catalog
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.