• The effect of Sm doping on the fcc-Al nucleation was investigated in Al-Sm liquids with low Sm concentrations (xSm) with molecular dynamics simulations. The nucleation in the moderately undercooled liquid is achieved by the recently developed persistent-embryo method. Systematically computing the nucleation rate with different xSm (xSm=0%, 1%, 2%, 3%, 5%) at 700 K, we found Sm dopant reduces the nucleation rate by up to 25 orders of magnitudes with only 5% doping concentration. This effect is mostly associated with the increase in the free energy barrier with a minor contribution from suppression of the attachment to the nucleus caused by Sm doping.
  • A metastable {\epsilon}-Al60Sm11 phase appears during the initial devitrification of as-quenched Al-10.2 at.% Sm glasses. The {\epsilon} phase is nonstoichiometric in nature since Al occupation is observed on the 16f Sm lattice sites. Scanning transmission electron microscopic images reveal profound spatial correlation of Sm content on these sites, which cannot be explained by the "average crystal" description from Rietveld analysis of diffraction data. Thermodynamically favorable configurations, established by Monte Carlo (MC) simulations based on a cluster-expansion model, also give qualitatively different correlation functions from experimental observations. On the other hand, molecular dynamics simulations of the growth of {\epsilon}-Al60Sm11 in undercooled liquid show that when the diffusion range of Sm is limited to ~ 4 {\AA}, the correlation function of the as-grown crystal structure agrees well with that of the STEM images. Our results show that kinetic effects, especially the limited diffusivity of Sm atoms plays the fundamental role in determining the nonstoichiometric site occupancies of the {\epsilon}-Al60Sm11 phase during the crystallization process.
  • The crystal nucleation from liquid in most cases is too rare to be accessed within the limited timescales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. Thus, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.
  • In this paper, we present a coupled experimental/theoretical investigation of pressure effect on the ferromagnetism of LaCrGe3 and LaCrSb3 compounds. The magnetic, electronic, elastic and mechanical properties of LaCrGe3 and LaCrSb3 at ambient condition are studied by first-principles density functional theory calculations. The pressure dependences of the magnetic properties of LaCrGe3 and LaCrSb3 are also investigated. The ferromagnetism in LaCrGe3 is rather fragile with a ferro- to paramagnetic transition at a relatively small pressure (around 7 GPa from our calculations, and 2 GPa in experiments). The key parameter controlling the magnetic properties of LaCrGe3 is found to be the proximity of the Cr DOS to the Fermi surface, a proximity that is strongly correlated to the distance between Cr atoms along the c-axis, suggesting that there would be a simple way to suppress magnetism in systems with one dimensional arrangement of magnetic atoms. By contrast, the ferromagnetism in LaCrSb3 is not fragile. Our calculation results are consistent with our experimental results and demonstrate the feasibility of using first-principles calculations to aid experimental explorations in screening for materials with fragile magnetism.
  • We performed a systematic search for low-energy structures of binary iron silicide over a wide range of compositions using the crystal structure prediction method based on adaptive genetic algorithm. 36 structures with formation energies within 50 meV/atom (11 of them are within 20 meV) above the convex hull formed by experimentally known stable structures are predicted. Magnetic properties of these low-energy structures are investigated. Some of these structures can be promising candidates for rare-earth-free permanent magnet.
  • The atomic structure of Ni64.5Zr35.5 metallic glass has been investigated by molecular dynamics (MD) simulations. The calculated structure factors from the MD glassy sample at room temperature agree well with the X-ray diffraction (XRD) and neutron diffraction (ND) experimental data. Using the pairwise cluster alignment and clique analysis methods, we show that there are three types dominant short-range order (SRO) motifs around Ni atoms in the glass sample of Ni64.5Zr35.5, i.e., Mixed-Icosahedron(ICO)-Cube, Twined-Cube and icosahedron-like clusters. Furthermore, chemical order and medium-range order (MRO) analysis show that the Mixed-ICO-Cube and Twined-Cube clusters exhibit the characteristics of the crystalline B2 phase. Our simulation results suggest that the weak glass-forming ability (GFA) of Ni64.5Zr35.5 can be attributed to the competition between the glass forming ICO SRO and the crystalline Mixed-ICO-Cube and Twined-Cube motifs.
  • A combination of theoretical calculation and the experimental synthesis to explore the new ternary compound is demonstrated in the Sr-Pt-Bi system. Since Pt-Bi is considered as a new critical charge-transfer pair for superconductivity, it inspired us to investigate the Sr-Pt-Bi system. With a thorough calculation of all the known stable/metastable compounds in the Sr-Pt-Bi system and crystal structure predictions, the thermodynamic stability of hypothetical stoichiometry, SrPtBi2, is determined. Followed by the high-temperature synthesis and crystallographic analysis, the first ternary bismuthide in Sr-Pt-Bi, SrPtBi2 was prepared and the stoichiometry was confirmed experimentally. SrPtBi2 crystallizes in the space group Pnma (S.G. 62, Pearson Symbol oP48), which matches well with theoretical prediction using an adaptive genetic algorithm (AGA). Using first-principles calculations, we demonstrate that the orthorhombic structure has lower formation energies than other 112 structure types, such as tetragonal BaMnBi2 (CuSmP2) and LaAuBi2 (CuHfSi2) structure types. The bonding analysis indicates the Pt-Bi interactions play a critical role in structural stability. The physical properties measurements show the metallic properties with low electron-phonon interactions at the low temperature, which agrees with the electronic structure assessment.
  • The highest pressure form of the major Earth-forming mantle silicate is MgSiO3 post-perovskite (PPv). Understanding the fate of PPv at TPa pressures is the first step for understanding the mineralogy of super-Earths-type exoplanets, arguably the most interesting for their similarities with Earth. Modeling their internal structure requires knowledge of stable mineral phases, their properties under compression, and major element abundances. Several studies of PPv under extreme pressures support the notion that a sequence of pressure induced dissociation transitions produce the elementary oxides SiO2 and MgO as the ultimate aggregation form at ~3 TPa. However, none of these studies have addressed the problem of mantle composition, particularly major element abundances usually expressed in terms of three main variables, the Mg/Si and Fe/Si ratios and the Mg#, as in the Earth. Here we show that the critical compositional parameter, the Mg/Si ratio, whose value in the Earth's mantle is still debated, is a vital ingredient for modeling phase transitions and internal structure of super-Earth mantles. Specifically, we have identified new sequences of phase transformations, including new recombination reactions that depend decisively on this ratio. This is a new level of complexity that has not been previously addressed, but proves essential for modeling the nature and number of internal layers in these rocky mantles.
  • We report the growth, structure and characterization of BaZnGa, identifying it as the sole known ternary compound in the Ba-Zn-Ga system. Single crystals of BaZnGa can be grown out of excess Ba-Zn and adopt a tI36 structure type. There are three unique Ba sites and three M\,=\,Zn/Ga sites. Using DFT calculations we can argue that whereas one of these three M sites is probably solely occupied by Ga, the other two M sites, most likely, have mixed Zn/Ga occupancy. Temperature dependent resistivity and magnetization measurements suggest that BaZnGa is a poor metal with no electronic or magnetic phase transitions between 2\,K and 300\,K.
  • Co3Si was recently reported to exhibit remarkable magnetic properties in the nanoparticle form [Appl. Phys. Lett. 108, 152406 (2016)], yet better understanding of this material is to be promoted. Here we report a study on the crystal structures of Co3Si using adaptive genetic algorithm, and discuss its electronic and magnetic properties from first-principles calculations. Several competing phases of Co3Si have been revealed from our calculations. We show that the hexagonal Co3Si structure reported in experiments has lower energy in non-magnetic state than ferromagnetic state at zero temperature. The ferromagnetic state of the hexagonal structure is dynamically unstable with imaginary phonon modes and transforms to a new orthorhombic structure, which is confirmed by our structure searches to have the lowest energy for both Co3Si and Co3Ge. Magnetic properties of the experimental hexagonal structure and the lowest-energy structures obtained from our structure searches are investigated in detail.
  • The dependence of the magnetocrystalline anisotropy energy (MAE) of MCo5 (M = Y, La) on the Coulomb correlations and strength of spin orbit (SO) interaction within the GGA + U scheme is investigated. A range of parameters suitable for the satisfactory description of key magnetic properties is determined. The origin of MAE in these materials is mostly related to the large orbital moment anisotropy of Co atoms on the 2c crystallographic site. Dependence of relativistic effects on Coulomb correlations, applicability of the second order perturbation theory for the description of MAE and effective screening of the SO interaction in these systems are discussed using a generalized virial theorem.
  • Transport and magnetic studies of PbTaSe$_2$ under pressure suggest existence of two superconducting phases with the low temperature phase boundary at $\sim 0.25$ GPa that is defined by a very sharp, first order, phase transition. The first order phase transition line can be followed via pressure dependent resistivity measurements, and is found to be near 0.12 GPa near room temperature. Transmission electron microscopy and x-ray diffraction at elevated temperatures confirm that this first order phase transition is structural and occurs at ambient pressure near $\sim 425$ K. The new, high temperature / high pressure phase has a similar crystal structure and slightly lower unit cell volume relative to the ambient pressure, room temperature structure. Based on first-principles calculations this structure is suggested to be obtained by shifting the Pb atoms from the $1a$ to $1e$ Wyckoff position without changing the positions of Ta and Se atoms. PbTaSe$_2$ has an exceptionally pressure sensitive, structural phase transition with $\Delta T_s/\Delta P \approx - 1700$ K/GPa near 4 K, this first order transition causes an $\sim 1$ K ($\sim 25 \%$) step - like decrease in $T_c$ as pressure is increased through 0.25 GPa.
  • Structures and magnetic properties of Fe16-xCoxN2 are studied using adaptive genetic algorithm and first-principles calculations. We show that substituting Fe by Co in Fe16N2 with Co/Fe ratio smaller than 1 can greatly improve the magnetic anisotropy of the material. The magnetocrystalline anisotropy energy from first-principles calculations reaches 3.18 MJ/m3 (245.6 {\mu}eV per metal atom) for Fe12Co4N2, much larger than that of Fe16N2 and is one of the largest among the reported rare-earth free magnets. From our systematic crystal structure searches, we show that there is a structure transition from tetragonal Fe16N2 to cubic Co16N2 in Fe16-xCoxN2 as the Co concentration increases, which can be well explained by electron counting analysis. Different magnetic properties between the Fe-rich (x < 8) and Co-rich (x > 8) Fe16-xCoxN2 is closely related to the structural transition.
  • The temperature-pressure phase diagram of the ferromagnet LaCrGe$_3$ is determined for the first time from a combination of magnetization, muon-spin-rotation and electrical resistivity measurements. The ferromagnetic phase is suppressed near $2.1$~GPa, but quantum criticality is avoided by the appearance of a magnetic phase, likely modulated, AFM$_Q$. Our density functional theory total energy calculations suggest a near degeneracy of antiferromagnetic states with small magnetic wave vectors $Q$ allowing for the potential of an ordering wave vector evolving from $Q=0$ to finite $Q$, as expected from the most recent theories on ferromagnetic quantum criticality. Our findings show that LaCrGe$_3$ is a very simple example to study this scenario of avoided ferromagnetic quantum criticality and will inspire further study on this material and other itinerant ferromagnets.
  • Using a motif-network search scheme, we studied the tetrahedral structures of the dilithium/disodium transition metal orthosilicates A2MSiO4 with A = Li or Na and M = Mn, Fe or Co. In addition to finding all previously reported structures, we discovered many other different tetrahedral-network-based crystal structures which are highly degenerate in energy. These structures can be classified into structures with 1D, 2D and 3D M-Si-O frameworks. A clear trend of the structural preference in different systems was revealed and possible indicators that affect the structure stabilities were introduced. For the case of Na systems which have been much less investigated in the literature relative to the Li systems, we predicted their ground state structures and found evidence for the existence of new structural motifs.
  • We develop a variational scheme called "Gutzwiller renormalization group" (GRG), which enables us to calculate the ground state of Anderson impurity models (AIM) with arbitrary numerical precision. Our method can exploit the low-entanglement property of the ground state in combination with the framework of the Gutzwiller wavefunction, and suggests that the ground state of the AIM has a very simple structure, which can be represented very accurately in terms of a surprisingly small number of variational parameters. We perform benchmark calculations of the single-band AIM that validate our theory and indicate that the GRG might enable us to study complex systems beyond the reach of the other methods presently available and pave the way to interesting generalizations, e.g., to nonequilibrium transport in nanostructures.
  • The insulating ground state of the 5d transition metal oxide CaIrO3 has been classified as a Mott-type insulator. Based on a systematic density functional theory (DFT) study with local, semilocal, and hybrid exchange-correlation functionals, we reveal that the Ir t2g states exhibit large splittings and one-dimensional electronic states along the c axis due to a tetragonal crystal field. Our hybrid DFT calculation adequately describes the antiferromagnetic (AFM) order along the c direction via a superexchange interaction between Ir4+ spins. Furthermore, the spin-orbit coupling (SOC) hybridizes the t2g states to open an insulating gap. These results indicate that CaIrO3 can be represented as a spin-orbit Slater insulator, driven by the interplay between a long-range AFM order and the SOC. Such a Slater mechanism for the gap formation is also demonstrated by the DFT + dynamical mean field theory calculation, where the metal-insulator transition and the paramagnetic to AFM phase transition are concomitant with each other.
  • We develop a new implementation of the Gutzwiller approximation (GA) and interface it with the local density approximation (LDA). This formulation enables us to study complex $4f$ and $5f$ systems. We perform calculations of praseodymium and $\alpha$-plutonium under pressure, which compare very well with the experiments. Our study of praseodymium indicates that both structure change and $f$-delocalization are important to obtain the correct phase diagram and, in particular, the pressure-induced volume-collapse transition. Our calculations of $\alpha$-plutonium indicate that, even though the $f$ electrons are delocalized in this phase, the electron-correlations affect substantially its electronic structure and thermodynamical properties.
  • The structures and magnetic properties of the Co-Zr-B alloys near the Co5Zr composition were studied using adaptive genetic algorithm and first-principles calculations to guide further experimental effort on optimizing their magnetic performances. Through extensive structural searches, we constructed the contour maps of the energetics and magnetic moments of the Co-Zr-B magnet alloys as a function of composition. We found that the Co-Zr-B system exhibits the same structural motif as the "Co11Zr2" polymorphs, which plays a key role in achieving high coercivity. Boron atoms can either substitute selective cobalt atoms or occupy the interstitial sites. First-principles calculation shows that the magnetocrystalline anisotropy energies can be significantly improved through proper boron doping.
  • We use angle-resolved photoemission spectroscopy to study heavy fermion superconductor Ce2RhIn8. The Fermi surface is rather complicated and consists of several hole and electron pock- ets. We do not observe kz dispersion of Fermi sheets, which is consistent with 2D character of the electronic structure. Comparison of the ARPES data and band structure calculations points to a localized picture of f electrons. Our findings pave the way for understanding the transport and thermodynamical properties of this material.
  • Gutzwiller wavefunction is a physically well motivated trial wavefunction for describing correlated electron systems. In this work, a new approximation is introduced to facilitate evaluation of the expectation value of any operator within the Gutzwiller wavefunction formalism. The basic idea is to make use of a specially designed average over Gutzwiller wavefunction coefficients expanded in the many-body Fock space to approximate the ratio of expectation values between a Gutzwiller wavefunction and its underlying noninteracting wavefunction. To check with the standard Gutzwiller approximation (GA), we test its performance on single band systems and find quite interesting properties. On finite systems, we noticed that it gives superior performance than GA, while on infinite systems it asymptotically approaches GA. Analytic analysis together with numerical tests are provided to support this claimed asymptotic behavior. At the end, possible improvements on the approximation and its generalization towards multiband systems are illustrated and discussed.
  • Phase selection in deeply undercooled liquids and devitrified glasses during heating involves complex interplay between the barriers to nucleation and the ability for these nuclei to grow. During the devitrification of glassy alloys, complicated metastable structures often precipitate instead of simpler, more stable compounds. Here, we access this unconventional type of phase selections by investigating an Al-10%Sm system, where a complicated cubic structure first precipitates with a large lattice parameter of 1.4 nm. We not only solve the structure of this "big cubic" phase containing ~140 atoms but establish an explicit interconnection between the structural orderings of the amorphous alloy and the cubic phase, which provides a low-barrier nucleation pathway at low temperatures. The surprising rapid growth of the crystal is attributed to its high tolerance to point defects, which minimize the short-scale atomic rearrangements to form the crystal. Our study suggests a new scenario of devitrification, where phase transformation proceeds initially without partitioning through a complex intermediate crystal phase.
  • We perform first-principle calculations of elemental cerium and compute its pressure-temperature phase diagram, finding good quantitative agreement with the experiments. Our calculations indicate that, while a signature of the volume-collapse transition appears in the free energy already at low temperatures, at larger temperatures this signature is enhanced because of the entropic effects, and originates an actual thermodynamical instability. Furthermore, we find that the catalyst determining this feature is --- in all temperature regimes --- a pressure-induced effective reduction of the $f$-level degeneracy due to the spin-orbit coupling. Our analysis suggests also that the lattice vibrations might be crucial in order to capture the behavior of the pressure-temperature transition line at large temperatures.
  • The underlying structural order that transcends the liquid, glass and crystalline states is identified using an efficient genetic algorithm (GA). GA identifies the most common energetically favorable packing motif in crystalline structures close to the alloy's Al-10 at.% Sm composition. These motifs are in turn compared to the observed packing motifs in the actual liquid structures using a cluster-alignment method which reveals the average topology. Conventional descriptions of the short-range order, such as Voronoi tessellation, are too rigid in their analysis of the configurational poly-types when describing the chemical and topological ordering during transition from undercooled metallic liquids to crystalline phases or glass. Our approach here brings new insight into describing mesoscopic order-disorder transitions in condensed matter physics.
  • Plutonium is the most exotic and mysterious element in the periodic table. It has 6 metallic phases and peculiar physical properties not yet understood. One of the most intriguing properties of Pu is that relatively small changes of temperature can induce transitions between different structures, that are accompanied by very large changes of equilibrium volumes. This fact has stimulated extensive theoretical and experimental studies. In spite of this, a convincing explanation of the metallurgic properties of Pu based on fundamental principles is still lacking, and none of the previous theories has been able to describe simultaneously the energetics and the $f$ electronic structure of all of the phases of Pu on the same footing. Here we provide a bird's eye view of Pu by studying the zero-temperature pressure-volume phase diagram and the $f$ electronic structure of all of its crystalline phases from first principles. In particular, we clarify the way in which the $f$-electron correlations determine its unusual energetics. Our theoretical energetics and ground-state $f$ electronic structure are both in good quantitative agreement with the experiments.