Source author record

Babak Sadigh

Babak Sadigh 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

18works
5topics
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

18 published item(s)

preprint2026arXiv

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. The method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

preprint2020arXiv

Self-Supervised Learning and Prediction of Microstructure Evolution with Recurrent Neural Networks

Microstructural evolution is a key aspect of understanding and exploiting the structure-property-performance relation of materials. Modeling microstructure evolution usually relies on coarse-grained simulations with evolution principles described by partial differential equations (PDEs). Here we demonstrate that convolutional recurrent neural networks can learn the underlying physical rules and replace PDE-based simulations in the prediction of microstructure phenomena. Neural nets are trained by self-supervised learning with image sequences from simulations of several common processes, including plane wave propagation, grain growth, spinodal decomposition and dendritic crystal growth. The trained networks can accurately predict both short-term local dynamics and long-term statistical properties of microstructures and is capable of extrapolating beyond the training datasets in spatiotemporal domains and configurational and parametric spaces. Such a data-driven approach offers significant advantages over PDE-based simulations in time stepping efficiency and offers a useful alternative especially when the material parameters or governing PDEs are not well determined.

preprint2016arXiv

Universal scaling of the electronic and the elastic energies of small polarons revealed by high-throughput first-principles calculations

Formation of self-trapped holes (STH) in a comprehensive list of scintillator materials, including halides and chalcogenides, are studied using an accurate and computationally efficient first-principles method, the polaron self-interaction correction (pSIC). The key characteristics of small hole polarons, including their geometries, energies and degree of localization, are found vastly different from halides to oxides to systems with open-shell cations. Nevertheless, we find a universal linear relation between the energy gap separating the bound hole level from the valence band maximum and the elastic energy associated with the lattice displacement field that accompanies the polaron.

preprint2015arXiv

A first-principles study of co-doping in lanthanum bromide

Co-doping of Ce-doped LaBr$_3$ with Ba, Ca, or Sr improves the energy resolution that can be achieved by radiation detectors based on these materials. Here, we present a mechanism that rationalizes of this enhancement that on the basis of first principles electronic structure calculations and point defect thermodynamics. It is shown that incorporation of Sr creates neutral $V_\text{Br}$-Sr$_\text{La}$ complexes that can temporarily trap electrons. As a result, Auger quenching of free carriers is reduced, allowing for a more linear, albeit slower, scintillation light yield response. Experimental Stokes shifts can be related to different Ce$_\text{La}$-Sr$_\text{La}$-$V_\text{Br}$ triple complex configurations. Co-doping with other alkaline as well as alkaline earth metals is considered as well. Alkaline elements are found to have extremely small solubilities on the order of 0.1 ppm and below at 1000 K. Among the alkaline earth metals the lighter dopant atoms prefer interstitial-like positions and create strong scattering centers, which has a detrimental impact on carrier mobilities. Only the heavier alkaline earth elements combine matching ionic radii with sufficiently high solubilities. This provides a rationale for the experimental finding that improved scintillator performance is exclusively achieved using Sr, Ca, or Ba. The present mechanism demonstrates that co-doping of wide gap materials can provide an efficient means for managing charge carrier populations under out-of-equilibrium conditions. In the present case dopants are introduced that manipulate not only the concentrations but the electronic properties of intrinsic defects without introducing additional gap levels. This leads to the availability of shallow electron traps that can temporarily localize charge carriers, effectively deactivating carrier-carrier recombination channels. The principles of this ... [continued]

preprint2015arXiv

A variational polaron self-interaction corrected total-energy functional for charge excitations in wide-band gap insulators

We conduct a detailed investigation of the polaron self-interaction (pSI) error in standard approximations to the exchange-correlation (XC) functional within density-functional theory (DFT). The pSI leads to delocalization error in the polaron wave function and energy, as calculated from the Kohn-Sham (KS) potential in the native charge state of the polaron. This constitutes the origin of the systematic failure of DFT to describe polaron formation in band insulators. It is shown that the delocalization error in these systems is, however, largely absent in the KS potential of the closed-shell neutral charge state. This leads to a modification of the DFT total-energy functional that corrects the pSI in the XC functional. The resulting pSIC-DFT method constitutes an accurate parameter-free {\it ab initio} methodology for calculating polaron properties in insulators at a computational cost that is orders of magnitude smaller than hybrid XC functionals. Unlike approaches that rely on parametrized localized potentials such as DFT+$U$, the pSIC-DFT method properly captures both site and bond-centered polaron configurations. This is demonstrated by studying formation and migration of self-trapped holes in alkali halides (bond-centered) as well as self-trapped electrons in an elpasolite compound (site-centered). The pSIC-DFT approach consistently reproduces the results obtained by hybrid XC functionals parametrized by DFT+$G_0W_0$ calculations. Finally, we generalize the pSIC approach to hybrid functionals, and show that in stark contrast to conventional hybrid calculations of polaron energies, the pSIC-hybrid method is insensitive to the parametrization of the hybrid XC functional. On this basis, we further rationalize the success of the pSIC-DFT approach.

preprint2014arXiv

Efficacy of the DFT+U formalism for modeling hole polarons in perovskite oxides

We investigate the formation of self-trapped holes (STH) in three prototypical perovskites (SrTiO3, BaTiO3, PbTiO3) using a combination of density functional theory (DFT) calculations with local potentials and hybrid functionals. First we construct a local correction potential for polaronic configurations in SrTiO3 that is applied via the DFT+U method and matches the forces from hybrid calculations. We then use the DFT+U potential to search the configuration space and locate the lowest energy STH configuration. It is demonstrated that both the DFT+U potential and the hybrid functional yield a piece-wise linear dependence of the total energy on the occupation of the STH level suggesting that self-interaction effects have been properly removed. The DFT+U model is found to be transferable to BaTiO3 and PbTiO3, and formation energies from DFT+U and hybrid calculations are in close agreement for all three materials. STH formation is found to be energetically favorable in SrTiO3 and BaTiO3 but not in PbTiO3, which can be rationalized by considering the alignment of the valence band edges on an absolute energy scale. In the case of PbTiO3 the strong coupling between Pb 6s and O 2p states lifts the valence band minimum (VBM) compared to SrTiO3 and BaTiO3. This reduces the separation between VBM and STH level and renders the STH configuration metastable with respect to delocalization (band hole state). We expect that the present approach can be adapted to study STH formation also oxides with different crystal structures and chemical composition.

preprint2014arXiv

Origin of resolution enhancement by co-doping of scintillators: Insight from electronic structure calculations

It was recently shown that the energy resolution of Ce-doped LaBr$_3$ scintillator radiation detectors can be crucially improved by co-doping with Sr, Ca, or Ba. Here we outline a mechanism for this enhancement on the basis of electronic structure calculations. We show that (i) Br vacancies are the primary electron traps during the initial stage of thermalization of hot carriers, prior to hole capture by Ce dopants; (ii) isolated Br vacancies are associated with deep levels; (iii) Sr doping increases the Br vacancy concentration by several orders of magnitude; (iv) $\text{Sr}_\text{La}$ binds to $V_\text{Br}$ resulting in a stable neutral complex; and (v) association with Sr causes the deep vacancy level to move toward the conduction band edge. The latter is essential for reducing the effective carrier density available for Auger quenching during thermalization of hot carriers. Subsequent de-trapping of electrons from $\text{Sr}_\text{La}-V_\text{La}$ complexes then can activate Ce dopants that have previously captured a hole leading to luminescence. This mechanism implies an overall reduction of Auger quenching of free carriers, which is expected to improve the linearity of the photon light yield with respect to the energy of incident electron or photon.

preprint2013arXiv

Free energy generalization of the Peierls potential in iron

In body-centered cubic (bcc) crystals, ${1}{2}111$ screw dislocations exhibit high intrinsic lattice friction as a consequence of their non-planar core structure, which results in a periodic energy landscape known as the Peierls potential, $U_P$. The main features determining plastic flow, including its stress and temperature dependences, can be derived directly from this potential, hence its importance. In this Letter, we use thermodynamic integration to provide a full thermodynamic extension of $U_P$ for bcc Fe. We compute the Peierls free energy path as a function of stress and temperature and show that the critical stress vanishes at 700K, supplying the qualitative elements that explain plastic behavior in the athermal limit.

preprint2013arXiv

Low-temperature criticality of martensitic transformations of Cu nanoprecipitates in α-Fe

Nanoprecipitates form during nucleation of multiphase equilibria in phase segregating multicomponent systems. In spite of their ubiquity, their size-dependent physical chemistry, in particular at the boundary between phases with incompatible topologies, is still rather arcane. Here we use extensive atomistic simulations to map out the size-temperature phase diagram of Cu nanoprecipitates in α-Fe. The growing precipitates undergo martensitic transformations from the body-centered cubic (BCC) phase to multiply-twinned 9R structures. At high temperatures, the transitions exhibit strong first-order character and prominent hysteresis. Upon cooling the discontinuities become less pronounced and the transitions occur at ever smaller cluster sizes. Below 300 K the hysteresis vanishes while the transition remains discontinuous with a finite but diminishing latent heat. This unusual size-temperature phase diagram results from the entropy generated by the soft modes of the BCC-Cu phase, which are stabilized through confinement by the α-Fe lattice.

preprint2013arXiv

Thermodynamic and mechanical properties of copper precipitates in alpha-iron from atomistic simulations

The Fe-Cu system has attracted much attention over the last several decades due to its technological importance as a model alloy for Cu steels. In spite of these efforts several aspects of its phase diagram remain unexplained. Here we use atomistic simulations to characterize the polymorphic phase diagram of Cu precipitates in body-centered cubic (BCC) Fe and establish a consistent link between their thermodynamic and mechanical properties in terms of thermal stability, shape, and strength. The size at which Cu precipitates transform from BCC to a close-packed 9R structure is found to be strongly temperature dependent, ranging from approximately 4 nm in diameter (~2,700 atoms) at 200 K to about 8 nm (~22,800 atoms) at 700 K. These numbers are in very good agreement with the interpretation of experimental data given Monzen et al. [Phil. Mag. A 80, 711 (2000)]. The strong temperature dependence originates from the entropic stabilization of BCC Cu, which is mechanically unstable as a bulk phase. While at high temperatures the transition exhibits first-order characteristics, the hysteresis, and thus the nucleation barrier, vanish at temperatures below approximately 300\,K. This behavior is explained in terms of the mutual cancellation of the energy differences between core and shell (wetting layer) regions of BCC and 9R nanoprecipitates, respectively. The proposed mechanism is not specific for the Fe--Cu system but could generally be observed in immiscible systems, whenever the minority component is unstable in the lattice structure of the host matrix. Finally, we also study the interaction of precipitates with screw dislocations as a function of both structure and orientation. The results provide a coherent picture of precipitate strength that unifies previous calculations and experimental observations.

preprint2012arXiv

A scalable parallel Monte Carlo algorithm for atomistic simulations of precipitation in alloys

We present an extension of the semi-grandcanonical (SGC) ensemble that we refer to as the variance-constrained semi-grandcanonical (VC-SGC) ensemble. It allows for transmutation Monte Carlo simulations of multicomponent systems in multiphase regions of the phase diagram and lends itself to scalable simulations on massively parallel platforms. By combining transmutation moves with molecular dynamics steps structural relaxations and thermal vibrations in realistic alloys can be taken into account. In this way, we construct a robust and efficient simulation technique that is ideally suited for large-scale simulations of precipitation in multicomponent systems in the presence of structural disorder. To illustrate the algorithm introduced in this work, we study the precipitation of Cu in nanocrystalline Fe.

preprint2012arXiv

Calculations of excess free energies of precipitates via direct thermodynamic integration across phase boundaries

We describe a technique for constraining macroscopic fluctuations in thermodynamic variables well-suited for Monte Carlo (MC) simulations of multiphase equilibria. In particular for multicomponent systems this amounts to a statistical ensemble that implements constraints on both the average composition as well as its fluctuations. The variance-constrained semi-grandcanonical (VC-SGC) ensemble allows for MC simulations, in which single-phase systems can be reversibly switched into multiphase equilibria allowing the calculation of excess free energies of precipitates of complex shapes by thermodynamic integration. The basic features as well as the scaling and convergence properties of this technique are demonstrated by application to an Ising model. Finally, the VC-SGC MC simulation technique is used to calculate alpha/alpha' interface free energies in Fe-Cr alloys as a function of orientation and temperature taking into account configurational, vibrational, and structural degrees of freedom.

preprint2012arXiv

Electronic structure of LaBr3 from quasi-particle self-consistent GW calculations

Rare-earth based scintillators in general and lanthanum bromide (LaBr_3) in particular represent a challenging class of materials due to pronounced spin-orbit coupling and subtle interactions between d and f states that cannot be reproduced by standard density functional theory (DFT). Here a detailed investigation of the electronic band structure of LaBr_3 using the quasi-particle self-consistent GW (QPscGW) method is presented. This parameter-free approach is shown to yield an excellent description of the electronic structure of LaBr_3. Specifically it is able to reproduce the band gap, the correct level ordering and spacing of the 4f and 5d states, as well as the spin-orbit splitting of La-derived states. The QPscGW results are subsequently used to benchmark several computationally less demanding techniques including DFT+U, hybrid exchange-correlation functionals, and the G_0W_0 method. Spin-orbit coupling is included self-consistently at each QPscGW iteration and maximally localized Wannier functions are used to interpolate quasi-particle energies. The QPscGW results provide an excellent starting point for investigating the electronic structure of excited states, charge self-trapping, and activator ions in LaBr_3 and related materials.

preprint2010arXiv

Are there stable long-range ordered Fe(1-x)Cr(x) compounds?

The heat of formation of Fe-Cr alloys undergoes an anomalous change of sign at small Cr concentrations. This observation raises the question whether there are intermetallic phases present in this composition range. Here we report the discovery of several long-range ordered structures that represent ground state phases at zero Kelvin. In particular we have identified a structure at 3.7% Cr with an embedding energy which is 49 meV/Cr atom below the solid solution. This implies there is an effective long-range attractive interaction between Cr atoms. We propose that the structures found in this study complete the low temperature-low Cr region of the phase diagram.

preprint2010arXiv

Pressure-induced phase transition in the electronic structure of palladium nitride

We present a combined theoretical and experimental study of the electronic structure and equation of state (EOS) of crystalline PdN2. The compound forms above 58 GPa in the pyrite structure and is metastable down to 11 GPa. We show that the EOS cannot be accurately described within either the local density or generalized gradient approximations. The Heyd-Scuseria-Ernzerhof exchange-correlation functional (HSE06), however, provides very good agreement with experimental data. We explain the strong pressure dependence of the Raman intensities in terms of a similar dependence of the calculated band gap, which closes just below 11 GPa. At this pressure, the HSE06 functional predicts a first-order isostructural transition accompanied by a pronounced elastic instability of the longitudinal-acoustic branches that provides the mechanism for the experimentally observed decomposition. Using an extensive Wannier function analysis, we show that the structural transformation is driven by a phase transition of the electronic structure, which is manifested by a discontinuous change in the hybridization between Pd-d and N-p electrons as well as a conversion from single to triple bonded nitrogen dimers. We argue for the possible existence of a critical point for the isostructural transition, at which massive fluctuations in both the electronic as well as the structural degrees of freedom are expected.

preprint2010arXiv

Short-range order and precipitation in Fe-rich Fe-Cr alloys: Atomistic off-lattice Monte Carlo simulations

Short-range order (SRO) in Fe-rich Fe-Cr alloys is investigated by means of atomistic off-lattice Monte Carlo simulations in the semi-grand canonical ensemble using classical interatomic potentials. The SRO parameter defined by Cowley [Phys. Rev. B 77, 669 (1950)] is used to quantify the degree of ordering. In agreement with experiments a strong ordering tendency in the Cr distribution at low Cr concentrations (~< 5%) is observed, as manifested in negative values of the SRO parameters. For intermediate Cr concentrations (5% ~< c_Cr ~< 15%) the SRO parameter for the alpha-phase goes through a minimum, but at the solubility limit the alpha-phase still displays a rather strong SRO. In thermodynamic equilibrium for concentrations within the two-phase region the SRO parameter measured over the entire sample therefore comprises the contributions from both the alpha and alpha-prime phases. If both of these contributions are taken into account, it is possible to quantitatively reproduce the experimental results and interpret their physical implications. It is thereby shown that the inversion of the SRO observed experimentally is due to the formation of stable (supercritical) alpha-prime precipitates. It is not related to the loss of SRO in the alpha-phase or to the presence of unstable (subcritical) Cr precipitates in the alpha-phase.

preprint2010arXiv

The Urbach tail in silica glass from first principles

We present density-functional theory calculations of the optical absorption spectra of silica glass for temperatures up to 2400 K. The calculated spectra exhibit exponential tails near the fundamental absorption edge that follow the Urbach rule, in good agreement with experiments. We also discuss the accuracy of our results by comparing to hybrid exchange correlation functionals. By deriving a simple relationship between the exponential tails of the absorption coefficient and the electronic density-of-states, we establish a direct link between the photoemission and the absorption spectra near the absorption edge. This relationship is subsequently employed to determine the lower bound to the Urbach frequency regime. Most interestingly, in this frequency interval, the optical absorption is Poisson distributed with very large statistical fluctuations. Finally, We determine the upper bound to the Urbach frequency regime by identifying the frequency at which transition to Poisson distribution takes place.

preprint2010arXiv

Thermodynamic ground states of platinum metal nitrides

The thermodynamic stabilities of various phases of the nitrides of the platinum metal elements are systematically studied using density functional theory. It is shown that for the nitrides of Rh, Pd, Ir and Pt two new crystal structures, in which the metal ions occupy simple tetragonal lattice sites, have lower formation enthalpies at ambient conditions than any previously proposed structures. The region of stability with respect to those structures extends to 17 GPa for PtN2. Calculations show that the PtN2 simple tetragonal structures at this pressure are thermodynamically stable also with respect to phase separation. The fact that the local density and generalized gradient approximations predict different values of the absolute formation enthalpies as well different relative stabilities between simple tetragonal and the pyrite or marcasite structures are further discussed.