Source author record

J. Bernholc

J. Bernholc appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci cond-mat.mes-hall physics.comp-ph quant-ph

Catalog footprint

What is connected

6works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2026arXiv

pyRMG: A Framework for High-Throughput, Large-Cell DFT Calculations on Supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code's grid-decomposition approach scales nearly linearly with the number of GPUs, even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global Fast Fourier Transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE computational materials science Python packages, pyRMG automates input generation and convergence checking, and integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. We demonstrate pyRMG for a high-throughput study of interfacial strain and twist-angle effects in lattice-matched, 2D Bi$_2$Se$_3$/NbSe$_2$ heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

preprint2020arXiv

Study of anharmonicity in Zirconium Hydrides using inelastic neutron scattering and ab-initio computer modeling

The anharmonic phenomena in Zirconium Hydrides and Deuterides, including ε-ZrH2, γ-ZrH, and γ-ZrD, have been investigated from aspects of inelastic neutron scattering (INS) and lattice dynamics calculations within the framework of density functional theory (DFT). The observed multiple sharp peaks below harmonic multi-phonon bands in the experimental spectra of all three materials did not show up in the simulated INS spectra based on the harmonic approximation, indicating the existence of strong anharmonicity in those materials and the necessity of further explanations. We present a detailed study on the anharmonicity of zirconium hydrides/deuterides by exploring the 2D potential energy surface of hydrogen/deuterium atoms, and solving the corresponding 2D single-particle Schrodinger equation to get the eigenfrequencies. The obtained results well describe the experimental INS spectra and show harmonic behavior in the fundamental modes and strong anharmonicity at higher energies.

preprint2015arXiv

Parallel implementation of electronic structure eigensolver using a partitioned folded spectrum method

A parallel implementation of an eigensolver designed for electronic structure calculations is presented. The method is applicable to computational tasks that solve a sequence of eigenvalue problems where the solution for a particular iteration is similar but not identical to the solution from the previous iteration. Such problems occur frequently when performing electronic structure calculations in which the eigenvectors are solutions to the Kohn-Sham equations. The eigenvectors are represented in some type of basis but the problem sizes are normally too large for direct diagonalization in that basis. Instead a subspace diagonalization procedure is employed in which matrix elements of the Hamiltonian operator are generated and the eigenvalues and eigenvectors of the resulting reduced matrix are obtained using a standard eigensolver from a package such as LAPACK or SCALAPACK. While this method works well and is widely used, the standard eigensolvers scale poorly on massively parallel computer systems for the matrix sizes typical of electronic structure calculations. We present a new method that utilizes a partitioned folded spectrum scheme (PFSM) that takes into account the iterative nature of the problem and performs well on massively parallel systems. Test results for a range of problems are presented that demonstrate an equivalent level of accuracy when compared to the standard eigensolvers, while also executing up to an order of magnitude faster. Unlike O(N) methods, the technique works equally well for metals and systems with unoccupied orbitals as for insulators and semiconductors. Timing and accuracy results are presented for a range of systems, including a 512 atom diamond cell, a cluster of 13 C60 molecules, bulk copper, a 216 atom silicon cell with a vacancy, using 40 unoccupied states/atom, and a 4000 atom aluminum supercell.

preprint2010arXiv

Band gap control of small bundles of carbon nanotubes using applied electric fields: A density functional theory study

Electrostatic screening between carbon nanotubes (CNTs) in a small CNT bundle leads to a switching behavior induced by electric field perpendicular to the bundle axis. Using a first-principles method, we investigate the electronic structures of bundles consisting of two or three CNTs and the effects of the electric field applied perpendicular to the bundle axis. The applied field causes band gap closure in semiconducting bundles, while a gap opening occurs in metallic ones, which enables considerable modulation of bundle conductivity. The modulation effect originates from symmetry breaking due to electrostatic screening between the adjacent tube walls.

preprint2010arXiv

Quantum-interference-controlled three-terminal molecular transistors based on a single ring-shaped-molecule connected to graphene nanoribbon electrodes

We study all-carbon-hydrogen molecular transistors where zigzag graphene nanoribbons play the role of three metallic electrodes connected to a ring-shaped 18-annulene molecule. Using the nonequilibrium Green function formalism combined with density functional theory, recently extended to multiterminal devices, we show that the proposed nanostructures exhibit exponentially small transmission when the source and drain electrodes are attached in a configuration that ensures destructive interference of electron paths around the ring. The third electrode, functioning either as an attached infinite-impedance voltage probe or as an "air-bridge" top gate covering half of molecular ring, introduces dephasing that brings the transistor into the "on" state with its transmission in the latter case approaching the maximum limit for a single conducting channel device. The current through the latter device can also be controlled in the far-from-equilibrium regime by applying a gate voltage.

preprint2009arXiv

First-principles methodology for quantum transport in multiterminal junctions

We present a generalized approach for computing electron conductance and I-V characteristics in multiterminal junctions from first-principles. Within the framework of Keldysh theory, electron transmission is evaluated employing an O(N) method for electronic-structure calculations. The nonequilibrium Green function for the nonequilibrium electron density of the multiterminal junction is computed self-consistently by solving Poisson equation after applying a realistic bias. We illustrate the suitability of the method on two examples of four-terminal systems, a radialene molecule connected to carbon chains and two crossed carbon chains brought together closer and closer. We describe charge density, potential profile, and transmission of electrons between any two terminals. Finally, we discuss the applicability of this technique to study complex electronic devices.