Source author record

Boris Kozinsky

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

16 published item(s)

preprint2026arXiv

Mixture of Experts Framework in Machine Learning Interatomic Potentials for Atomistic Simulations

First-principles atomistic simulations are essential for understanding complex material phenomena but are fundamentally limited by their computational cost. While Machine Learning Interatomic Potentials (MLIPs) have drastically improved cost for a given accuracy, their inference cost remains a bottleneck for massive systems or long timescales. To address this, we introduce a multifidelity "Mixture-of-Experts" framework based on the E(3)-equivariant Allegro architecture. Our method spatially partitions the simulation domain into a chemically complex region (e.g., reactive interfaces) and a simple region (e.g., bulk lattice), assigning models of varying capacity to each. Among the challenges in such static domain decomposition, the mechanical mismatch between models at the interface is particularly critical, as it can generate artificial stress fields and instability. We address this challenge with a co-training strategy in which the loss function includes agreement constraints -- penalties on per-atom energy and force discrepancies between models evaluated on shared bulk environments -- forcing the independent models to learn a consistent physical description of the bulk material. We validate this approach on a realistic Pt+CO catalytic system, demonstrating that the co-trained models maintain exact energy conservation, align their bulk mechanical response (e.g., equation of state and bulk modulus), and achieve predictive accuracy comparable to a full high-fidelity simulation at more than twice the computational speed.

preprint2022arXiv

CIDER: An Expressive, Nonlocal Feature Set for Machine Learning Density Functionals with Exact Constraints

Machine learning (ML) has recently gained attention as a means to develop more accurate exchange-correlation (XC) functionals for density functional theory, but functionals developed thus far need to be improved on several metrics, including accuracy, numerical stability, and transferability across chemical space. In this work, we introduce a set of nonlocal features of the density called the CIDER formalism, which we use to train a Gaussian process model for the exchange energy that obeys the critical uniform scaling rule for exchange. The resulting CIDER exchange functional is significantly more accurate than any semi-local functional tested here, and it has good transferability across main-group molecules. This work therefore serves as an initial step toward more accurate exchange functionals, and it also introduces useful techniques for developing robust, physics-informed XC models via ML.

preprint2022arXiv

Engineering ideal helical topological networks in stanene via Zn decoration

The xene family of topological insulators plays a key role in many proposals for topological electronic, spintronic, and valleytronic devices. These proposals rely on applying local perturbations, including electric fields and proximity magnetism, to induce topological phase transitions in xenes. However, these techniques lack control over the geometry of interfaces between topological regions, a critical aspect of engineering topological devices. We propose adatom decoration as a method for engineering atomically precise topological edge modes in xenes. Our first-principles calculations show that decorating stanene with Zn adatoms exclusively on one of two sublattices induces a topological phase transition from the quantum spin Hall (QSH) to quantum valley Hall (QVH) phase and confirm the existence of spin-valley polarized edge modes propagating at QSH/QVH interfaces. We conclude by discussing technological applications of these edge modes that are enabled by the atomic precision afforded by recent advances in adatom manipulation technology.

preprint2022arXiv

Learning Local Equivariant Representations for Large-Scale Atomistic Dynamics

A simultaneously accurate and computationally efficient parametrization of the energy and atomic forces of molecules and materials is a long-standing goal in the natural sciences. In pursuit of this goal, neural message passing has lead to a paradigm shift by describing many-body correlations of atoms through iteratively passing messages along an atomistic graph. This propagation of information, however, makes parallel computation difficult and limits the length scales that can be studied. Strictly local descriptor-based methods, on the other hand, can scale to large systems but do not currently match the high accuracy observed with message passing approaches. This work introduces Allegro, a strictly local equivariant deep learning interatomic potential that simultaneously exhibits excellent accuracy and scalability of parallel computation. Allegro learns many-body functions of atomic coordinates using a series of tensor products of learned equivariant representations, but without relying on message passing. Allegro obtains improvements over state-of-the-art methods on the QM9 and revised MD-17 data sets. A single tensor product layer is shown to outperform existing deep message passing neural networks and transformers on the QM9 benchmark. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular dynamics simulations based on Allegro recover structural and kinetic properties of an amorphous phosphate electrolyte in excellent agreement with first principles calculations. Finally, we demonstrate the parallel scaling of Allegro with a dynamics simulation of 100 million atoms.

preprint2022arXiv

Micron-scale heterogeneous catalysis with Bayesian force fields from first principles and active learning

Quantum-mechanically accurate reactive molecular dynamics (MD) at the scale of billions of atoms has been achieved for the heterogeneous catalytic system of H$_2$/Pt(111) using the FLARE Bayesian force field. This achievement provides accelerated time-to-solution from first principles, with Bayesian active learning enabling efficient and autonomous training of the machine learning model. The resulting model is then deployed in LAMMPS on GPUs using the Kokkos performance portability library. The Bayesian force field provides quantitative uncertainty of predictions on every atomic environment, critical for detecting configurations in large reactive simulations that are outside of the training set. Scaling benchmarks were performed using real-application MD of the H$_2$/Pt(111) heterogeneous catalysis on the Summit supercomputer, with simulations reaching 0.5 trillion atoms on 4556 GPU nodes.

preprint2022arXiv

Microscopic picture of paraelectric perovskites from structural prototypes

We show with first-principles molecular dynamics the persistence of intrinsic $\langle111\rangle$ Ti off-centerings for BaTiO$_3$ in its cubic paraelectric phase. Intriguingly, these are inconsistent with the Pm$\bar 3$m space group often used to atomistically model this phase using density functional theory or similar methods. Therefore we deploy a systematic symmetry analysis to construct representative structural models in the form of supercells that satisfy a desired point symmetry but are built from the combination of lower-symmetry primitive cells. We define as structural prototypes the smallest of these that are both energetically and dynamically stable. Remarkably, two 40-atom prototypes can be identified for paraelectric BaTiO$_3$; these are also common to many other $AB$O$_3$ perovskites. These prototypes can offer structural models of paraelectric phases that can be used for the computational engineering of functional materials. Last, we show that the emergence of B-cation off-centerings and the primitive-cell phonon instabilities are controlled by the equilibrium volume, in turn dictated by the filler A cation.

preprint2022arXiv

Phoebe: a High-Performance Framework for Solving Phonon and Electron Boltzmann Transport Equations

Understanding the electrical and thermal transport properties of materials is critical to the design of electronics, sensors and energy conversion devices. Computational modeling can accurately predict materials properties but, in order to be reliable, require accurate descriptions of electron and phonon states and their interactions. While first-principles methods are capable of describing the energy spectrum of each carrier, using them to compute transport properties is still a formidable task, both computationally demanding and memory intensive, requiring integration of fine microscopic scattering details for estimation of macroscopic transport properties. To address this challenge, we present Phoebe - a newly developed software package that includes the effects of electron-phonon, phonon-phonon, boundary, and isotope scattering in computations of electrical and thermal transport properties of materials with a variety of available methods and approximations. This open source C++ code combines MPI-OpenMP hybrid parallelization with GPU acceleration and distributed memory structures to manage computational cost, allowing Phoebe to effectively take advantage of contemporary computing infrastructures. We demonstrate that Phoebe accurately and efficiently predicts a wide range of transport properties, opening avenues for accelerated computational analysis of complex crystals.

preprint2021arXiv

E(3)-Equivariant Graph Neural Networks for Data-Efficient and Accurate Interatomic Potentials

This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

preprint2020arXiv

Accurate and scalable multi-element graph neural network force field and molecular dynamics with direct force architecture

Recently, machine learning (ML) has been used to address the computational cost that has been limiting ab initio molecular dynamics (AIMD). Here, we present GNNFF, a graph neural network framework to directly predict atomic forces from automatically extracted features of the local atomic environment that are translationally-invariant, but rotationally-covariant to the coordinate of the atoms. We demonstrate that GNNFF not only achieves high performance in terms of force prediction accuracy and computational speed on various materials systems, but also accurately predicts the forces of a large MD system after being trained on forces obtained from a smaller system. Finally, we use our framework to perform an MD simulation of Li7P3S11, a superionic conductor, and show that resulting Li diffusion coefficient is within 14% of that obtained directly from AIMD. The high performance exhibited by GNNFF can be easily generalized to study atomistic level dynamics of other material systems.

preprint2015arXiv

AiiDA: Automated Interactive Infrastructure and Database for Computational Science

Computational science has seen in the last decades a spectacular rise in the scope, breadth, and depth of its efforts. Notwithstanding this prevalence and impact, it is often still performed using the renaissance model of individual artisans gathered in a workshop, under the guidance of an established practitioner. Great benefits could follow instead from adopting concepts and tools coming from computer science to manage, preserve, and share these computational efforts. We illustrate here our paradigm sustaining such vision, based around the four pillars of Automation, Data, Environment, and Sharing. We then discuss its implementation in the open-source AiiDA platform (http://www.aiida.net), that has been tuned first to the demands of computational materials science. AiiDA's design is based on directed acyclic graphs to track the provenance of data and calculations, and ensure preservation and searchability. Remote computational resources are managed transparently, and automation is coupled with data storage to ensure reproducibility. Last, complex sequences of calculations can be encoded into scientific workflows. We believe that AiiDA's design and its sharing capabilities will encourage the creation of social ecosystems to disseminate codes, data, and scientific workflows.

preprint2015arXiv

Effects of sublattice symmetry and frustration on ionic transport in garnet solid electrolytes

We use rigorous group-theoretic techniques and molecular dynamics to investigate the connection between structural symmetry and ionic conductivity in the garnet family of solid Li-ion electrolytes. We identify new ordered phases and order-disorder phase transitions that are relevant for conductivity optimization. Ionic transport in this materials family is controlled by the frustration of the Li sublattice caused by incommensurability with the host structure at non-integer Li concentrations, while ordered phases explain regions of sharply lower conductivity. Disorder is therefore predicted to be optimal for ionic transport in this and other conductor families with strong Li interaction.

preprint2015arXiv

Electron-phonon-averaged approximation for first-principles computations of electron relaxation times and transport properties in semiconductor materials

We present a simple and efficient approximation to the electron-phonon scattering rate suitable for high-throughput screening of candidate materials for thermoelectric devices, based on electronic transport. The method is applied to calculate the electronic transport coefficients of half-Heusler compounds, showing agreement with experimental data. By directly computing electrical and the electronic part of the thermal conductivities, we find deviations from the Wiedemann-Franz law in these compounds at high temperatures and low carrier concentrations.

preprint2014arXiv

Electron-Phonon Interactions and the Intrinsic Electrical Resistivity of Graphene

We present a first-principles study of the temperature- and density-dependent intrinsic electrical resistivity of graphene. We use density-functional theory and density-functional perturbation theory together with very accurate Wannier interpolations to compute all electronic and vibrational properties and electron-phonon coupling matrix elements; the phonon-limited resistivity is then calculated within a Boltzmann-transport approach. An effective tight-binding model, validated against first-principles results, is also used to study the role of electron-electron interactions at the level of many-body perturbation theory. The results found are in excellent agreement with recent experimental data on graphene samples at high carrier densities and elucidate the role of the different phonon modes in limiting electron mobility. Moreover, we find that the resistivity arising from scattering with transverse acoustic phonons is 2.5 times higher than that from longitudinal acoustic phonons. Last, high-energy, optical, and zone-boundary phonons contribute as much as acoustic phonons to the intrinsic electrical resistivity even at room temperature and become dominant at higher temperatures.

preprint2013arXiv

BoltzWann: A code for the evaluation of thermoelectric and electronic transport properties with a maximally-localized Wannier functions basis

We present a new code to evaluate thermoelectric and electronic transport properties of extended systems with a maximally-localized Wannier function basis set. The semiclassical Boltzmann transport equations for the homogeneous infinite system are solved in the constant relaxation-time approximation and band energies and band derivatives are obtained via Wannier interpolations. Thanks to the exponential localization of the Wannier functions obtained, very high accuracy in the Brillouin zone integrals can be achieved with very moderate computational costs. Moreover, the analytical expression for the band derivatives in the Wannier basis resolves any issues that may occur when evaluating derivatives near band crossings. The code is tested on binary and ternary skutterudites CoSb_3 and CoGe_{3/2}S_{3/2}.

preprint2013arXiv

High-Throughput Screening of Perovskite Alloys for Piezoelectric Performance and Formability

We screen a large chemical space of perovskite alloys for systems with the right properties to accommodate a morphotropic phase boundary (MPB) in their composition-temperature phase diagram, a crucial feature for high piezoelectric performance. We start from alloy end-points previously identified in a high-throughput computational search. An interpolation scheme is used to estimate the relative energies between different perovskite distortions for alloy compositions with a minimum of computational effort. Suggested alloys are further screened for thermodynamic stability. The screening identifies alloy systems already known to host a MPB, and suggests a few new ones that may be promising candidates for future experiments. Our method of investigation may be extended to other perovskite systems, e.g., (oxy-)nitrides, and provides a useful methodology for any application of high-throughput screening of isovalent alloy systems.

preprint2011arXiv

Electronic, vibrational and transport properties of pnictogen substituted ternary skutterudites

First principles calculations are used to investigate electronic band structure and vibrational spectra of pnictogen substituted ternary skutterudites. We compare the results with the prototypical binary composition CoSb$_3$ to identify the effects of substitutions on the Sb site, and evaluate the potential of ternary skutterudites for thermoelectric applications. Electronic transport coefficients are computed within the Boltzmann transport formalism assuming a constant relaxation time, using a new methodology based on maximally localized Wannier function interpolation. Our results point to a large sensitivity of the electronic transport coefficients to carrier concentration and to scattering mechanisms associated with the enhanced polarity. The ionic character of the bonds is used to explain the detrimental effect on the thermoelectric properties.