Quasi-relativistic approach to analytical gradients of parity violating potentials.

An analytic gradient approach for the computation of derivatives of parity-violating (PV) potentials with respect to displacements of the nuclei in chiral molecules is described and implemented within a quasirelativistic mean-field framework. Calculated PV potential gradients are utilized for estimating PV frequency splittings between enantiomers in rotational and vibrational spectra of four chiral polyhalomethanes, i.e., CHBrClF, CHClFI, CHBrFI, and CHAtFI. Values calculated within the single-mode approximation for frequency shifts agree well with previously reported theoretical values. The influence of non-separable anharmonic effects (multi-mode effects) on vibrational frequency shifts, which are readily accessible with the present analytic derivative approach, is estimated for the C-F stretching fundamental of all four molecules and computed for each of the fundamentals in CHBrClF and CHAtFI. Multi-mode effects are found to be significant, in particular, for C-F stretching modes, being for some modes and cases of similar size as the single-mode contribution.

Microcanonical RT-TDDFT simulations of realistically extended devices.

In this paper, real-time time-dependent density functional theory (RT-TDDFT) calculations of realistically sized nanodevices are presented. These microcanonical simulations rely on a closed boundary approach based on recent advances in the software package CP2K. The obtained results are compared to those derived from the open-boundary Non-equilibrium Green's Function (NEGF) formalism. A good agreement between the "current vs. voltage" characteristics produced by both methods is demonstrated for three representative device structures, a carbon nanotube field-effect transistor, a GeSe selector for crossbar arrays, and a conductive bridging random-access memory cell. Different approaches to extract the electrostatic contribution from the RT-TDDFT Hamiltonian and to incorporate the result into the NEGF calculations are presented.

Efficient algorithms for large-scale quantum transport calculations.

Massively parallel algorithms are presented in this paper to reduce the computational burden associated with quantum transport simulations from first-principles. The power of modern hybrid computer architectures is harvested in order to determine the open boundary conditions that connect the simulation domain with its environment and to solve the resulting Schrödinger equation. While the former operation takes the form of an eigenvalue problem that is solved by a contour integration technique on the available central processing units (CPUs), the latter can be cast into a linear system of equations that is simultaneously processed by SplitSolve, a two-step algorithm, on general-purpose graphics processing units (GPUs). A significant decrease of the computational time by up to two orders of magnitude is obtained as compared to standard solution methods.

A generalized Poisson solver for first-principles device simulations.

Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated.