Researcher profile

Peter Koval

Peter Koval contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2022arXiv

Inelastic scattering of electrons in water from first principles: cross sections and inelastic mean free path for use in Monte Carlo track-structure simulations of biological damage

Modelling the inelastic scattering of electrons in water is fundamental, given their crucial role in biological damage. In Monte Carlo track-structure codes used to assess biological damage, the energy loss function, from which cross sections are extracted, is derived from different semi-empirical optical models. Only recently, first ab-initio results for the energy loss function and cross-sections in water became available. For benchmarking purpose, in this work, we present ab-initio linear-response time-dependent density functional theory calculations of the energy loss function of liquid water. We calculated the inelastic scattering cross sections, inelastic mean free paths, and electronic stopping powers and compared our results with recent calculations and experimental data showing a good agreement. In addition, we provide an in-depth analysis of the contributions of different molecular orbitals, species, and orbital angular momenta to the total energy loss function. Moreover, we present single-differential cross sections computed for each molecular orbital channel, which should prove useful for Monte-Carlo track-structure simulations.

preprint2020arXiv

SIESTA: recent developments and applications

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as Wannier90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low level libraries, as well as data standards and support for them, in particular the PSML definition and library for transferable pseudopotentials, and the interface to the ELSI library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

preprint2011arXiv

An O(N^3) implementation of Hedin's GW approximation

Organic electronics is a rapidly developing technology. Typically, the molecules involved in organic electronics are made up of hundreds of atoms, prohibiting a theoretical description by wavefunction-based ab-initio methods. Density-functional and Green's function type of methods scale less steeply with the number of atoms. Therefore, they provide a suitable framework for the theory of such large systems. In this contribution, we describe an implementation, for molecules, of Hedin's GW approximation. The latter is the lowest order solution of a set of coupled integral equations for electronic Green's and vertex functions that was found by Lars Hedin half a century ago. Our implementation of Hedin's GW approximation has two distinctive features: i) it uses sets of localized functions to describe the spatial dependence of correlation functions, and ii) it uses spectral functions to treat their frequency dependence. Using these features, we were able to achieve a favorable computational complexity of this approximation. In our implementation, the number of operations grows as N^3 with the number of atoms N.

preprint2010arXiv

A Parallel Iterative Method for Computing Molecular Absorption Spectra

We describe a fast parallel iterative method for computing molecular absorption spectra within TDDFT linear response and using the LCAO method. We use a local basis of "dominant products" to parametrize the space of orbital products that occur in the LCAO approach. In this basis, the dynamical polarizability is computed iteratively within an appropriate Krylov subspace. The iterative procedure uses a a matrix-free GMRES method to determine the (interacting) density response. The resulting code is about one order of magnitude faster than our previous full-matrix method. This acceleration makes the speed of our TDDFT code comparable with codes based on Casida's equation. The implementation of our method uses hybrid MPI and OpenMP parallelization in which load balancing and memory access are optimized. To validate our approach and to establish benchmarks, we compute spectra of large molecules on various types of parallel machines. The methods developed here are fairly general and we believe they will find useful applications in molecular physics/chemistry, even for problems that are beyond TDDFT, such as organic semiconductors, particularly in photovoltaics.

preprint2010arXiv

Fast construction of the Kohn--Sham response function for molecules

The use of the LCAO (Linear Combination of Atomic Orbitals) method for excited states involves products of orbitals that are known to be linearly dependent. We identify a basis in the space of orbital products that is local for orbitals of finite support and with a residual error that vanishes exponentially with its dimension. As an application of our previously reported technique we compute the Kohn--Sham density response function $χ_{0}$ for a molecule consisting of $N$ atoms in $N^{2}N_ω$ operations, with $N_ω$ the number of frequency points. We test our construction of $χ_{0}$ by computing molecular spectra directly from the equations of Petersilka--Gossmann--Gross in $N^{2}N_ω$ operations rather than from Casida's equations which takes $N^{3}$ operations. We consider the good agreement with previously calculated molecular spectra as a validation of our construction of $χ_{0}$. Ongoing work indicates that our method is well suited for the computation of the GW self-energy $Σ=\mathrm{i}GW$ and we expect it to be useful in the analysis of exitonic effects in molecules.