
The effect of Sm doping on the fccAl nucleation was investigated in AlSm
liquids with low Sm concentrations (xSm) with molecular dynamics simulations.
The nucleation in the moderately undercooled liquid is achieved by the recently
developed persistentembryo 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 asquenched Al10.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 clusterexpansion 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 asgrown 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 glassforming CuZr
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 firstprinciples 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 caxis, 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 firstprinciples
calculations to aid experimental explorations in screening for materials with
fragile magnetism.

We performed a systematic search for lowenergy 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 lowenergy structures are investigated. Some of
these structures can be promising candidates for rareearthfree 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 Xray 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 shortrange order (SRO) motifs around Ni atoms in the glass
sample of Ni64.5Zr35.5, i.e., MixedIcosahedron(ICO)Cube, TwinedCube and
icosahedronlike clusters. Furthermore, chemical order and mediumrange order
(MRO) analysis show that the MixedICOCube and TwinedCube clusters exhibit
the characteristics of the crystalline B2 phase. Our simulation results suggest
that the weak glassforming ability (GFA) of Ni64.5Zr35.5 can be attributed to
the competition between the glass forming ICO SRO and the crystalline
MixedICOCube and TwinedCube motifs.

A combination of theoretical calculation and the experimental synthesis to
explore the new ternary compound is demonstrated in the SrPtBi system. Since
PtBi is considered as a new critical chargetransfer pair for
superconductivity, it inspired us to investigate the SrPtBi system. With a
thorough calculation of all the known stable/metastable compounds in the
SrPtBi system and crystal structure predictions, the thermodynamic stability
of hypothetical stoichiometry, SrPtBi2, is determined. Followed by the
hightemperature synthesis and crystallographic analysis, the first ternary
bismuthide in SrPtBi, 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 firstprinciples 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 PtBi
interactions play a critical role in structural stability. The physical
properties measurements show the metallic properties with low electronphonon
interactions at the low temperature, which agrees with the electronic structure
assessment.

The highest pressure form of the major Earthforming mantle silicate is
MgSiO3 postperovskite (PPv). Understanding the fate of PPv at TPa pressures is
the first step for understanding the mineralogy of superEarthstype
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 superEarth 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 BaZnGa system. Single crystals
of BaZnGa can be grown out of excess BaZn 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 firstprinciples 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
nonmagnetic 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
lowestenergy 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 xray 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
firstprinciples 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 Fe16xCoxN2 are studied using adaptive
genetic algorithm and firstprinciples 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 firstprinciples 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 rareearth free magnets. From our systematic crystal
structure searches, we show that there is a structure transition from
tetragonal Fe16N2 to cubic Co16N2 in Fe16xCoxN2 as the Co concentration
increases, which can be well explained by electron counting analysis. Different
magnetic properties between the Ferich (x < 8) and Corich (x > 8) Fe16xCoxN2
is closely related to the structural transition.

The temperaturepressure phase diagram of the ferromagnet LaCrGe$_3$ is
determined for the first time from a combination of magnetization,
muonspinrotation 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 motifnetwork 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 tetrahedralnetworkbased
crystal structures which are highly degenerate in energy. These structures can
be classified into structures with 1D, 2D and 3D MSiO 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
lowentanglement 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 singleband 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 Motttype insulator. Based on a systematic density functional
theory (DFT) study with local, semilocal, and hybrid exchangecorrelation
functionals, we reveal that the Ir t2g states exhibit large splittings and
onedimensional 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 spinorbit coupling (SOC) hybridizes the t2g states to
open an insulating gap. These results indicate that CaIrO3 can be represented
as a spinorbit Slater insulator, driven by the interplay between a longrange
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
metalinsulator 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 pressureinduced volumecollapse transition. Our
calculations of $\alpha$plutonium indicate that, even though the $f$ electrons
are delocalized in this phase, the electroncorrelations affect substantially
its electronic structure and thermodynamical properties.

The structures and magnetic properties of the CoZrB alloys near the Co5Zr
composition were studied using adaptive genetic algorithm and firstprinciples
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 CoZrB magnet alloys as a
function of composition. We found that the CoZrB 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. Firstprinciples calculation shows that
the magnetocrystalline anisotropy energies can be significantly improved
through proper boron doping.

We use angleresolved 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 manybody 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 Al10%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 lowbarrier nucleation pathway at
low temperatures. The surprising rapid growth of the crystal is attributed to
its high tolerance to point defects, which minimize the shortscale 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 firstprinciple calculations of elemental cerium and compute its
pressuretemperature phase diagram, finding good quantitative agreement with
the experiments. Our calculations indicate that, while a signature of the
volumecollapse 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 pressureinduced effective reduction of the $f$level
degeneracy due to the spinorbit coupling. Our analysis suggests also that the
lattice vibrations might be crucial in order to capture the behavior of the
pressuretemperature 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 Al10 at.% Sm composition. These motifs are in
turn compared to the observed packing motifs in the actual liquid structures
using a clusteralignment method which reveals the average topology.
Conventional descriptions of the shortrange order, such as Voronoi
tessellation, are too rigid in their analysis of the configurational polytypes
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 orderdisorder 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 zerotemperature pressurevolume 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 groundstate $f$
electronic structure are both in good quantitative agreement with the
experiments.