Source author record

Axel Van de Walle

Axel Van de Walle 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

13works
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

13 published item(s)

preprint2022arXiv

Rapid screening of high-throughput ground state predictions

High-through computational thermodynamic approaches are becoming an increasingly popular tool to uncover novel compounds. However, traditional methods tend to be limited to stability predictions of stoichiometric phases at absolute zero. Such methods thus carry the risk of identifying an excess of possible phases that do not survive to temperatures of practical relevance. We demonstrate how the Calphad formalism, informed by simple first-principles input can be simply used to overcome this problem at a low computational cost and deliver quantitatively useful phase diagram predictions at all temperatures. We illustrate the method by re-assessing prior compound formation predictions and reconcile these findings with long-standing experimental evidence to the contrary.

preprint2020arXiv

Re-entrant melting of sodium, magnesium, and aluminum: General trend

Re-entrant melting (in which a substance's melting point starts to decrease beyond a certain pressure) is believed to be an unusual phenomenon. Among the elements, it has so far only been observed in a very limited number of species, e.g., the alkali metals. Our density functional theory calculations reveal that this behavior actually extends beyond alkali metals to include magnesium, which also undergoes re-entrant melting, though at the much higher pressure of ~300 GPa. We find that the origin of re-entrant melting is the faster softening of interatomic interactions in the liquid phase than in the solid, as pressure rises. We propose a simple approach to estimate pressure-volume relations and show that this characteristic softening pattern is widely observed in metallic elements. We verify this prediction in the case of aluminum by finding re-entrant melting at ~4000 GPa. These results suggest that re-entrant melting may be a more universal feature than previously thought.

preprint2014arXiv

Ab initio calculation of anisotropic interfacial excess free energies

We describe a simple method to determine, from ab initio calculations, the complete orientation-dependence of interfacial free energies in solid-state crystalline systems. We illustrate the method with an application to precipitates in the Al-Ti alloy system. The method combines the cluster expansion formalism in its most general form (to model the system's energetics) with the inversion of the well-known Wulff construction (to recover interfacial energies from equilibrium precipitate shapes). Although the inverse Wulff construction only provides the relative magnitude of the various interfacial free energies, absolute free energies can be recovered from a calculation of a single, conveniently chosen, planar interface. The method is able to account for essentially all sources of entropy (arising from phonons, bulk point defects, as well as interface roughness) and is thus able to transparently handle both atomically smooth and rough interfaces. The approach expresses the resulting orientation-dependence of the interfacial properties using symmetry-adapted bases for general orientation-dependent quantities. As a by-product, this paper thus provides a simple and general method to generate such basis functions, which prove useful in a variety of other applications, for instance to represent the anisotropy of the so-called constituent strain elastic energy.

preprint2013arXiv

Accelerated Molecular Dynamics through stochastic iterations to strengthen yield of path hopping over upper states (SISYPHUS)

We present a new method, called SISYPHUS (Stochastic Iterations to Strengthen Yield of Path Hopping over Upper States), for extending accessible time-scales in atomistic simulations. The method proceeds by separating phase space into basins, and transition regions between the basins based on a general collective variable (CV) criterion. The transition regions are treated via traditional molecular dynamics (MD) while Monte Carlo (MC) methods are used to (i) estimate the expected time spent in each basin and (ii) thermalize the system between two MD episodes. In particular, an efficient adiabatic switching based scheme is used to estimate the time spent inside the basins. The method offers various advantages over existing approaches in terms of (i) providing an accurate real time scale, (ii) avoiding reliance on harmonic transition state theory and (iii) avoiding the need to enumerate all possible transition events. Applications of SISYPHUS to low temperature vacancy diffusion in BCC Ta and adatom island ripening in FCC Al are presented. A new CV appropriate for such condensed phases, especially for transitions involving collective motions of several atoms, is also introduced.

preprint2013arXiv

Realistic time-scale fully atomistic simulations of surface nucleation of dislocations in pristine nanopillars

We use our recently proposed accelerated dynamics algorithm (Tiwary & van de Walle, 2011) to calculate temperature and stress dependence of activation free energy for surface nucleation of dislocations in pristine Gold nanopillars under realistic loads. While maintaining fully atomistic resolution, we achieve the fraction of a second time-scale regime. We find that the activation free energy depends significantly on the driving force (stress or strain) and temperature, leading to very high activation entropies. We also perform compression tests on Gold nanopillars for strain rates varying between 7 orders of magnitudes, reaching as low as 10^3/s. Our calculations show the quantitative effects on the yield point of unrealistic strain-rate Molecular Dynamics calculations: we find that while the failure mechanism for <001> compression of Gold nanopillars remains the same across the entire strain-rate range, the elastic limit (defined as stress for nucleation of the first dislocation) depends significantly on the strain-rate. We also propose a new methodology that overcomes some of the limits in our original accelerated dynamics scheme (and accelerated dynamics methods in general). We lay out our methods in sufficient details so as to be used for understanding and predicting deformation mechanism under realistic driving forces for various problems.

preprint2012arXiv

Ab initio thermodynamics of intrinsic oxygen vacancies in ceria

Nonstoichiometric ceria(CeO$_{2-δ}$) is a candidate reaction medium to facilitate two step water splitting cycles and generate hydrogen. Improving upon its thermodynamic suitability through doping requires an understanding of its vacancy thermodynamics. Using density functional theory(DFT) calculations and a cluster expansion based Monte Carlo simulations, we have studied the high temperature thermodynamics of intrinsic oxygen vacancies in ceria. The DFT+$U$ approach was used to get the ground state energies of various vacancy configurations in ceria, which were subsequently fit to a cluster expansion Hamiltonian to efficiently model the configurational dependence of energy. The effect of lattice vibrations was incorporated through a temperature dependent cluster expansion. Lattice Monte Carlo simulations using the cluster expansion Hamiltonian were able to detect the miscibility gap in the phase diagram of ceria. The inclusion of vibrational and electronic entropy effects made the agreement with experiments quantitative. The deviation from an ideal solution model was quantified by calculating as a function of nonstoichiometry, a) the solid state entropy from Monte Carlo simulations and b) Warren-Cowley short range order parameters of various pair clusters.

preprint2012arXiv

Intrinsic Defects and Dopability of Zinc Phosphide

Zinc Phosphide ($Zn_3P_2$) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type $Zn_3P_2$ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. This is consistent with experimental observations of both the tendency of conductivity to rise with phosphorus partial pressure, and with current partial successes in n-type doping in very zinc-rich growth conditions.

preprint2011arXiv

First Principles Phase Diagram Calculations for the Octahedral-Interstitial System $HfO_{X}$, $0 \leq X \leq 1/2$

First principles based phase diagram calculations were performed for the hexagonal closest packed octahedral-interstitial solid solution system $αHfO_{X} ($αHf[ ]_{1-X}O_{X}$; [ ]=Vacancy; $0 \leq X \leq 1/2$). The cluster expansion method was used to do a ground state analysis, and to calculate the phase diagram. The predicted diagram has four ordered ground-states in the range $0 \leq X \leq 1/2$, but one of these, at X=5/12, is predicted to disproportionate at T$\approx 220K$. At X$\approx 1/3 ($Hf_{3}O$) and X$\approx 1/2 ($Hf_{2}O$), order-parameter vs temperature plots predict Devil's Staircases of closely related ordered structures.

preprint2011arXiv

First Principles Phase Diagram Calculations for the Octahedral-Interstitial System ZrO$_{X}$, $0 \leq X \leq 1/2$

First principles based phase diagram calculations were performed for the octahedral-interstitial solid solution system αZrOX (αZr[ ]_(1-X)OX; [ ]=Vacancy; 0 \leq X \leq 1/2). The cluster expansion method was used to do a ground state analysis, and to calculate the phase diagram. The predicted diagram has four ordered ground-states in the range 0 \leq X \leq 1/2, but one of these, at X=5/12, is predicted to disproportionate at T \approx 20K, well below the experimentally investigated range T \approx 420K. Thus, at T \succeq 420K, the first-principles based calculation predicts three ordered phases rather than the four that have been reported by experimentalists.

preprint2011arXiv

Hybrid deterministic and stochastic approach for efficient atomistic simulations at long time scales

We propose a hybrid deterministic and stochastic approach to achieve extended time scales in atomistic simulations that combines the strengths of molecular dynamics (MD) and Monte Carlo (MC) simulations in an easy-to-implement way. The method exploits the rare event nature of the dynamics similar to most current accelerated MD approaches but goes beyond them by providing, without any further computational overhead, (a) rapid thermalization between infrequent events, thereby minimizing spurious correlations, and (b) control over accuracy of time-scale correction, while still providing similar or higher boosts in computational efficiency. We present two applications of the method: (a) Vacancy-mediated diffusion in Fe yields correct diffusivities over a wide range of temperatures and (b) source-controlled plasticity and deformation behavior in Au nanopillars at realistic strain rates (10^4/s and lower), with excellent agreement with previous theoretical predictions and in situ high-resolution transmission electron microscopy observations. The method gives several orders-of-magnitude improvements in computational efficiency relative to standard MD and good scalability with the size of the system.

preprint2010arXiv

Entropic Stabilization and Retrograde Solubility in Zn4Sb3

Zn4Sb3 is shown to be entropically stabilized versus decomposition to Zn and ZnSb though the effects of configurational disorder and phonon free energy. Single phase stability is predicted for a range of compositions and temperatures. Retrograde solubility of Zn is predicted on the two-phase boundary region between Zn4Sb3 and Zn. The complex temperature dependent solubility can be used to explain the variety of nanoparticle formation observed in the system: formation of ZnSb on the Sb rich side, Zn on the far Zn rich side and nano-void formation due to Zn precipitates being reabsorbed at lower temperatures.

preprint2010arXiv

Interatomic potentials for mixed oxide (MOX) nuclear fuels

We extend our recently developed interatomic potentials for UO_{2} to the mixed oxide fuel system (U,Pu,Np)O_{2}. We do so by fitting against an extensive database of ab initio results as well as to experimental measurements. The applicability of these interactions to a variety of mixed environments beyond the fitting domain is also assessed. The employed formalism makes these potentials applicable across all interatomic distances without the need for any ambiguous splining to the well-established short-range Ziegler-Biersack-Littmark universal pair potential. We therefore expect these to be reliable potentials for carrying out damage simulations (and Molecular Dynamics simulations in general) in nuclear fuels of varying compositions for all relevant atomic collision energies.

preprint2009arXiv

Building effective models from sparse but precise data

A common approach in computational science is to use a set of of highly precise but expensive calculations to parameterize a model that allows less precise, but more rapid calculations on larger scale systems. Least-squares fitting on a model that underfits the data is generally used for this purpose. For arbitrarily precise data free from statistic noise, e.g. ab initio calculations, we argue that it is more appropriate to begin with a ensemble of models that overfit the data. Within a Bayesian framework, a most likely model can be defined that incorporates physical knowledge, provides error estimates for systems not included in the fit, and reproduces the original data exactly. We apply this approach to obtain a cluster expansion model for the Ca[Zr,Ti]O3 solid solution.