BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Haifeng Song

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.str-elmath-phmath.MPquant-phphysics.app-phphysics.chem-phphysics.comp-ph

Catalog footprint

What is connected

12works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

12 published item(s)

preprint2022arXiv

Melting curve of magnesium up to 460 GPa from ab initio molecular dynamics simulations

Based on ab initio molecular dynamics simulations, we determined the melting curve of magnesium (Mg) up to ~460 GPa using the solid-liquid coexistence method. Between ~30 and 100 GPa, our melting curve is noticeably lower than those from static experiments, but is in good agreement with recent shock experiments. Up to ~450 GPa, our melting curve is generally consistent with the melting points from first-principles calculations using the small-cell coexistence method. We found that, at high pressures of a few hundred GPa, due to the strong softening of interatomic interactions in the liquid phase, solid-liquid coexistence simulations of Mg show some characteristics distinctively different from other metal systems, such as aluminum. For example, at a given volume, the pressure and temperature range for maintaining a stable solid-liquid coexistence state can be very small. The strong softening in the liquid phase also causes the unusual behavior of reentrant melting to occur at very high pressures. The onset of reentrant melting is predicted at ~305 GPa, close to that at ~300 GPa from the small-cell coexistence method. We show that the calculated melting points, considering reentrant melting, can be excellently fitted to a low-order Kechin equation, thereby making it possible for us to obtain a first-principles melting curve of Mg at pressures above 50 GPa for the first time. Similar characteristics in solid-liquid coexistence simulations, as well as reentrant melting, are also expected for other systems with strong softening in the liquid phase at high pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

First-Principles Study of Hydrogen Behaviors in $α$-Pu$_{2}$O$_{3}$

The in-depth understanding of hydrogen permeation through plutonium-oxide overlayers is the prerequisite to evaluate the complex hydriding induction period of Pu. In this work, the incorporation, diffusion and dissolution of hydrogen in $α$-Pu$_{2}$O$_{3}$ are investigated by the first-principles calculations and $\textit{ab initio}$ thermodynamic method based on DFT+U and DFT-D3 schemes. Our study reveals that the hydrogen incorporation is endothermic and the separated H atoms prefer to recombine as H$_{2}$ molecules rather than reacting with $α$-Pu$_{2}$O$_{3}$. The H and H$_{2}$ diffusion are both feasible, generally, H will recombine first as H$_{2}$ and then migrate. Both pressure P$_{H2}$ and temperature can promote the hydrogen dissolution in $α$-Pu$_{2}$O$_{3}$. The single H$_{2}$ molecule incorporation and (H+H$_{2}$) mixed dissolution will successively appear when increasing P$_{H2}$. Compared to PuO$_{2}$, this work indicates that Pu sesquioxide is hardly reduced by hydrogen, but the porous $α$-Pu$_{2}$O$_{3}$ facilitates hydrogen transport in Pu oxide layers. We presents the microscopic picture of hydrogen behaviors in the defect-free $α$-Pu$_{2}$O$_{3}$, which could shed some light on the study of the hydriding induction period of Pu.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spin-Triplet Excitonic Insulator: The Case of Graphone

While various excitonic insulators have been studied in the literature, due to the perceived too-small spin splitting, spin-triplet excitonic insulator is rare. In two-dimensional systems such as a graphone, however, it is possible, as revealed by first-principles calculations coupled with Bethe-Salpeter equation. The critical temperature, given by an effective Hamiltonian, is 11.5 K. While detecting excitonic insulators is still a daunting challenge, the condensation of triplet excitons will result in spin superfluidity, which can be directly measured by a transport experiment. Nonlocal dielectric screening also leads to an unexpected phenomenon, namely, an indirect-to-direct transition crossover between single-particle band and exciton dispersion in graphone, which offers yet another test by experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

UO2/BeO interfacial thermal resistance and its effect on fuel thermal conductivity

UO2/BeO interfacial thermal resistance (ITR) is calculated by diffuse mismatch model (DMM) and the effects of ITR on UO2-BeO thermal conductivity are investigated. ITR predicted by DMM is on the order of 10-9 m2K/W. Using this ITR, UO2-BeO thermal conductivities are calculated by theoretical models and compared with experimental data. The results indicate that DMM prediction is applicable to the interface between UO2 and dispersed BeO, while not applicable to the interface between UO2 and continuous BeO. If the thermal conductivity of UO2 containing continuous BeO was to be in agreement with experimental data, its ITR should be on the order of 10-6 - 10-5 m2K/W. Therefore, the vibrational mismatch between UO2 and BeO considered by DMM is the major mechanism for attenuating the heat flux through UO2/dispersed-BeO interface, but not for UO2/continuous-BeO interface. Furthermore, it is found that the presence of ITR leads to the dependence of the thermal conductivity of UO2 containing dispersed BeO on BeO size. With the decrease in BeO size, UO2-BeO thermal conductivity decreases. When BeO size is smaller than a critical value, UO2-BeO thermal conductivity becomes even smaller than UO2 thermal conductivity. For UO2 containing continuous BeO, the thermal conductivity decreases with the decrease in the size of UO2 granule surrounded by BeO, but not necessarily smaller than UO2 thermal conductivity. Under a critical temperature, UO2-BeO thermal conductivity is always larger than UO2 thermal conductivity. Above the critical temperature, UO2-BeO thermal conductivity is larger than UO2 thermal conductivity only when UO2 granule size is large enough. The conditions for achieving the targeted enhancement of UO2 thermal conductivity by doping with BeO are derived. These conditions can be used to design and optimize the distribution, content, size of BeO, and the size of UO2 granule.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics

Wavefunction extrapolation greatly reduces the number of self-consistent field (SCF) iterations and thus the overall computational cost of Born-Oppenheimer molecular dynamics (BOMD) that is based on the Kohn-Sham density functional theory. Going against the intuition that the higher order of extrapolation possesses a better accuracy, we demonstrate, from both theoretical and numerical perspectives, that the extrapolation accuracy firstly increases and then decreases with respect to the order, and an optimal extrapolation order in terms of minimal number of SCF iterations always exists. We also prove that the optimal order tends to be larger when using larger MD time steps or more strict SCF convergence criteria. By example BOMD simulations of a solid copper system, we show that the optimal extrapolation order covers a broad range when varying the MD time step or the SCF convergence criterion. Therefore, we suggest the necessity for BOMD simulation packages to open the user interface and to provide more choices on the extrapolation order. Another factor that may influence the extrapolation accuracy is the alignment scheme that eliminates the discontinuity in the wavefunctions with respect to the atomic or cell variables. We prove the equivalence between the two existing schemes, thus the implementation of either of them does not lead to essential difference in the extrapolation accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Three-dimensional bulk electronic structure of the Kondo lattice CeIn3 revealed by photoemission

We show the three-dimensional electronic structure of the Kondo lattice CeIn3 using soft x-ray angle resolved photoemission spectroscopy in the paramagnetic state. For the first time, we have directly observed the three-dimensional topology of the Fermi surface of CeIn3 by photoemission. The Fermi surface has a complicated hole pocket centred at the Γ-Z line and an elliptical electron pocket centred at the R point of the Brillouin zone. Polarization and photon-energy dependent photoemission results both indicate the nearly localized nature of the 4f electrons in CeIn3, consistent with the theoretical prediction by means of the combination of density functional theory and single-site dynamical meanfield theory. Those results illustrate that the f electrons of CeIn3, which is the parent material of CeMIn5 compounds, are closer to the localized description than the layered CeMIn5 compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale dynamics of open three-level quantum systems with two quasi-degenerate levels

We consider a three-level quantum system interacting with a bosonic thermal reservoir. Two energy levels of the system are nearly degenerate but well separated from the third one. The system-reservoir interaction constant is larger than the energy difference of the degenerate levels, but it is smaller than the separation between the latter and the remaining level. We show that the quasi-degeneracy of energy levels leads to the existence of a manifold of quasi-stationary states, and the dynamics exhibits two characteristic time scales. On the first, shorter one, initial states approach the quasi-stationary manifold. Then, on the much longer second time scale, the final unique equilibrium is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The different roles of Pu-oxide overlayers in the hydrogenation of Pu-metal: An ab initio molecular dynamics study based on vdW-DFT+U

Based on the van der Waals density functional theory (vdW-DFT)+U scheme, we carry out the ab initio molecular dynamics (AIMD) study of the interaction dynamics for H$_{2}$ impingement against the stoichiometric PuO$_{2}$(111), the reduced PuO$_{2}$(111), and the stoichiometric $α$-Pu$_{2}$O$_{3}$(111) surfaces. The hydrogen molecular physisorption states, which can not be captured by pure DFT+\textit{U} method, are obtained by employing the vdW-DFT+\textit{U} scheme. We show that except for the weak physisorption, PuO$_{2}$(111) surfaces are so difficult of access that almost all of the H$_{2}$ molecules will bounce back to the vacuum when their initial kinetic energies are not sufficient. Although the dissociative adsorption of H$_{2}$ on PuO$_{2}$(111) surfaces is found to be very exothermic, the collision-induced dissociation barriers of H$_{2}$ are calculated to be as high as $3.2$ eV and $2.0$ eV for stoichiometric and reduced PuO$_{2}$ surfaces, respectively. Unlike PuO$_{2}$, our AIMD study directly reveals that the hydrogen molecules can penetrate into $α$-Pu$_{2}$O$_{3}$(111) surface and diffuse easily due to the $25$\ native O vacancies located along the $\langle $111$\rangle $ diagonals of $α$-Pu$_{2}$O$_{3}$ matrix. By examining the temperature effect and the internal vibrational excitations of H$_{2}$, we provide a detailed insight into the interaction dynamics of H$_{2}$ in $α$-Pu$_{2}$O$_{3}$. The optimum pathways for hydrogen penetration and diffusion, the corresponding energy barriers ($1.0$ eV and $0.53$ eV, respectively) and rate constants are systematically calculated. Overall, our study fairly reveals the different interaction mechanisms between H$_{2}$ and Pu-oxide surfaces, which have strong implications to the interpretation of experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A rigorous derivation of the transition rate for the dissipative two-state system

The transition rate for a two-state system interacting with a bosonic heat bath, from the initial state `up' at time t=0 to `down' at time t>0, was derived formally in the seminal paper [12] by Leggett, Chakravarty, Dorsey, Fisher, Garg and Zwerger. The expression was obtained for small tunneling parameter of the two-level system using time-dependent perturbation theory to second order. It contains all orders in the parameter of coupling between the system and the bath. We give a rigorous derivation of the transition rate, using a quantum resonance approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First-principles study of surface properties of PuO2: Effects of thickness and O-vacancy on surface stability and chemical activity

The (111), (110), and (001) surfaces properties of PuO2 are studied by using density-functional theory+U method. The total-energy static calculations determine the relative order of stability for low-index PuO2 surfaces, namely, O-terminated (111) > (110) > defective (001) > polar (001). The effect of thickness is shown to modestly modulate the surface stability and chemical activity of the (110) surface. The high work function of 6.19 eV indicates the chemical inertia of the most stable (111) surface, and the surface O-vacancy with concentration C_V=25% can efficiently lower the work function to 4.35 eV, which is a crucial indicator of the difference in the surface chemical activities between PuO2 and α-Pu2O3. For the polar (001) surface, 50% on-surface O-vacancy can effectively quench the dipole moment and stabilize the surface structure, where the residual surface oxygen atoms are arranged in a zigzag manner along the <100> direction. We also investigate the relative stability of PuO2 surfaces in an oxygen environment. Under oxygen-rich conditions, the stoichiometric O-terminated (111) is found to be the most stable surface. Whereas under O-reducing conditions, the on-surface O-vacancy of C_V = 1/9 is stable, and for high reducing conditions, the (111) surface with nearly one monolayer subsurface oxygen removed (C_V = 8/9) becomes most stable.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The environmental dependence of redox energetics of PuO2 and α-Pu2O3: A quantitative solution from DFT+U calculations

We report a comprehensive density functional theory (DFT) + $U$ study of the energetics of charged and neutral oxygen defects in both PuO$_{2}$ and $α$-Pu$_{2}$O$_{3}$, and present a quantitative determination of the equilibrium compositions of reduced PuO$_{2}$ (PuO$_{2-x}$) as functions of environmental temperature and partial pressure of oxygen, which shows fairly agreement with corresponding high-temperature experiments. Under ambient conditions, the endothermic reduction of PuO$_{2}$ to $α$-Pu$_{2}$O$_{3} $ is found to be facilitated by accompanying volume expansion of PuO$_{2-x}$ and the possible migration of O-vacancy, whereas further reduction of $α$-Pu$_{2}$O$_{3}$ is predicted to be much more difficult. In contrast to the endothermic oxidation of PuO$_{2}$,\ the oxidation of $α$-Pu$_{2} $O$_{3}$ is a stable exothermic process.

0 citations0 reviewsNo code linkedOpen access