Source author record

Xiulin Ruan

Xiulin Ruan 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

14works
3topics
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

14 published item(s)

preprint2021arXiv

Atomistic metrics of BaSO$_4$ as an ultra-efficient radiative cooling material: a first-principles prediction

Radiative cooling has recently revived due to its significant potential as an environmentally friendly cooling technology. However, the design of particle-matrix cooling nanocomposites was generally carried out via tedious trial-and-error approaches, and the atomistic physics for efficient radiative cooling was not well understood. In this work, we identify the atomistic metrics of Barium Sulfate (BaSO$_4$) nanocomposite, which is an ultra-efficient radiative cooling material, using a predictive first-principles approach coupled with Monte Carlo simulations. Our results show that BaSO$_4$-acrylic nanocomposites not only attain high total solar reflectance of 92.5% (0.28 - 4.0 um), but also simultaneously demonstrate high normal emittance of 96.0% in the sky window region (8 - 13 um), outperforming the commonly used $α$-quartz ($α$-SiO$_2$). We identify two pertinent characters of ultra-efficient radiative cooling paints: i) a balanced band gap and refractive index, which enables strong scattering while negating absorption in the solar spectrum, and ii) a sufficient number of infrared-active optical resonance phonon modes resulting in abundant Reststrahlen bands and high emissivity in the sky window. The first principles approach and the resulted physical insights in this work pave the way for further search of ultra-efficient radiative cooling materials.

preprint2021arXiv

Raman Linewidth Contributions from Four-Phonon and Electron-Phonon Interactions in Graphene

The Raman peak position and linewidth provide insight into phonon anharmonicity and electron-phonon interactions (EPI) in materials. For monolayer graphene, prior first-principles calculations have yielded decreasing linewidth with increasing temperature, which is opposite to measurement results. Here, we explicitly consider four-phonon anharmonicity, phonon renormalization, and electron-phonon coupling, and find all to be important to successfully explain both the $G$ peak frequency shift and linewidths in our suspended graphene sample at a wide temperature range. Four-phonon scattering contributes a prominent linewidth that increases with temperature, while temperature dependence from EPI is found to be reversed above a doping threshold ($\hbarω_G/2$, with $ω_G$ being the frequency of the $G$ phonon).

preprint2019arXiv

Role of higher-order phonon scattering in the zone-center optical phonon linewidth and the Lorenz oscillator model

Zone-center optical phonon linewidth is a key parameter for infrared and Raman spectra as well as the Lorenz oscillator model. While three-phonon scattering was often assumed to be the leading contribution, in this work we find, surprisingly, that higher-order phonon scattering universally plays a significant or even dominant role over three-phonon scattering at room temperature, and more so at elevated temperatures, for a wide range of materials including diamond, Si, Ge, boron arsenide (BAs), cubic silicon carbide (3C-SiC), and $α$-quartz. This is enabled by the large fourphonon scattering phase space of zone-center optical phonons, and distinct from heat conduction where at room temperature four-phonon scattering is still secondary to three-phonon scattering. Moreover, our results imply that five-phonon and even higher-order scattering may be significant for some large band-gap materials, e.g., BAs. Our predicted infrared optical properties through the Lorenz oscillator model, after including four-phonon scattering, show much better agreement with experimental measurements than those three-phonon based predictions. This work will raise broad interest of studying high-order scattering in various areas beyond heat conduction.

preprint2019arXiv

The use of strain and grain boundaries to tailor phonon transport properties: A first principles study of 2H-phase $CuAlO_{2}$ (Part II)

Transparent oxide materials, such as $CuAlO_{2}$, a p-type transparent conducting oxide (TCO), have recently been studied for high temperature thermoelectric power generators and coolers for waste heat. TCO materials are generally low cost and non-toxic. The potential to engineer them through strain and nano-structuring are two promising avenues toward continuously tuning the electronic and thermal properties to achieve high zT values and low cost/kW-hr devices. In this work, the strain-dependent lattice thermal conductivity of 2H $CuAlO_{2}$ is computed by solving the phonon Boltzmann transport equation with interatomic force constants extracted from first-principles calculations. While the average bulk thermal conductivity is around 32 W/(K-m) at room temperature, it drops to between 5-15 W/(K-m) for typical experimental grain sizes from 3nm to 30nm at room temperature. We find that strain can offer both an increase as well as a decrease in the thermal conductivity as expected, however the overall inclusion of small grain sizes dictates the potential for low thermal conductivity in this material.

preprint2019arXiv

Thermal boundary resistance predictions with non-equilibrium Green's function and molecular dynamics simulations

The non-equilibrium Green's function (NEGF) method with Büttiker probe scattering self-energies is assessed by comparing its predictions for the thermal boundary resistance with molecular dynamics (MD) simulations. For simplicity, the interface of Si/heavy-Si is considered, where heavy-Si differs from Si only in the mass value. With Büttiker probe scattering parameters tuned against MD in homogeneous Si, the NEGF-predicted thermal boundary resistance quantitatively agrees with MD for wide mass ratios. Artificial resistances that the unaltered Landauer approach yield at virtual interfaces in homogeneous systems are absent in the present NEGF approach. Spectral information result from NEGF in its natural representation without further transformations. The spectral results show that the scattering between different phonon modes plays a crucial role in thermal transport across interfaces. Büttiker probes provide an efficient and reliable way to include anharmonicity in phonon related NEGF. NEGF including the Büttiker probes can reliably predict phonon transport across interfaces and at finite temperatures.

preprint2016arXiv

Quantum mechanical prediction of four-phonon scattering rates and reduced thermal conductivity of solids

Recently, first principle-based prediction of lattice thermal conductivity $κ$ from the perturbation theory has achieved significant success. However, it only includes three-phonon scattering due to the assumption that four-phonon and higher-order processes are generally unimportant. Also, directly evaluating the scattering rates of four-phonon and higher-order processes has been a long-standing challenge. In this work, however, we have developed a formalism to explicitly determine quantum mechanical scattering probability matrices for four-phonon scattering in the full Brillouin Zone, and by mitigating the computational challenge we have directly calculated four-phonon scattering rates. We find that four-phonon scattering rates are comparable to three-phonon scattering rates at medium and high temperatures, and they increase quadratically with temperature. As a consequence, $κ$ of Lennard-Jones argon is reduced by more than 60% at 80 K when four-phonon scattering is included. Also, in less anharmonic materials -- diamond, silicon, and germanium, $κ$ is still reduced considerably at high temperature by four-phonon scattering. Also, the thermal conductivity of optical phonons is dominated by the fourth and higher orders phonon scattering even at low temperature.

preprint2016arXiv

Uncertainty quantification of thermal conductivities from equilibrium molecular dynamics simulations

Equilibrium molecular dynamics (EMD) simulations along with the Green-Kubo formula have been widely used to calculate lattice thermal conductivities. Previous studies, however, focused primarily on the calculated thermal conductivities, with the uncertainty of the thermal conductivities remaining poorly understood. In this paper, we study the quantification of the uncertainty by using solid argon, silicon, and germanium as model material systems, and examine the origin of the observed uncertainty. We find that the uncertainty increases with the upper limit of the correlation time, $t_{\text{corre, UL}}$, and decreases with the total simulation time, $t_{\text{total}}$, whereas the velocity initialization seed, simulation domain size, temperature, and type of material have minimal effects. The relative uncertainties of the thermal conductivities for solid argon, silicon, and germanium under different simulation conditions all follow a similar trend, which can be fit with a "universal" square-root relation, as $σ_{k_x}/k_{x, \text{ave}} = 2(t_{\text{total}}/t_{\text{corre, UL}})^{-0.5}$. We have also conducted statistical analysis of the EMD-predicted thermal conductivities and derived a formula that correlates the relative error bound ($Q$), confidence level ($P$), $t_{\text{corre, UL}}$, $t_{\text{total}}$, and number of independent simulations ($N$). We recommend choosing $t_{\text{corre, UL}}$ to be 5-10 times the effective phonon relaxation time, $τ_{\text{eff}}$, and choosing $t_{\text{total}}$ and $N$ based on the desired relative error bound and confidence level. This study provides new insights into understanding the uncertainty of EMD-predicted thermal conductivities. It also provides a guideline for running EMD simulations to achieve a desired relative error bound with a desired confidence level and for reporting EMD-predicted thermal conductivities.

preprint2015arXiv

Coupling between Phonon-Phonon and Phonon-Impurity Scattering: A Critical Revisit of the Spectral Matthiessen's Rule

The spectral Matthiessen's rule is commonly used to calculate the total phonon scattering rate when multiple scattering mechanisms exist. Here we predict the spectral phonon relaxation time $τ$ of defective bulk silicon using normal mode analysis based on molecular dynamics, and show that the spectral Matthiessen's rule is not accurate due to the neglect of the coupling between anharmonic phonon-phonon scattering $τ_a^{-1}$ and phonon-impurity scattering $τ_i^{-1}$. As a result, the spectral Matthiessen's rule underestimates the total phonon scattering rate, and hence overestimates the thermal conductivity $κ$ of mass-doped and Ge-doped silicon by about 20-40%. We have also directly estimated this coupling scattering rate, so called coupled five-phonon scattering $τ_{\rm couple}^{-1}$, and achieved good agreement between $τ_a^{-1}+τ_i^{-1}+τ_{\rm couple}^{-1}$ and the total scattering rate $τ_{tot}^{-1}$.

preprint2015arXiv

First-Principles Simulation of Electron Mean-Free-Path Spectra and Thermoelectric Properties in Silicon

The mean-free-paths (MFPs) of energy carriers are of critical importance to the nano-engineering of better thermoelectric materials. Despite significant progress in the first-principles-based understanding of the spectral distribution of phonon MFPs in recent years, the spectral distribution of electron MFPs remains unclear. In this work, we compute the energy dependent electron scatterings and MFPs in silicon from first-principles. The electrical conductivity accumulation with respect to electron MFPs is compared to that of the phonon thermal conductivity accumulation to illustrate the quantitative impact of nanostructuring on electron and phonon transport. By combining all electron and phonon transport properties from first-principles, we predict the thermoelectric properties of the bulk and nanostructured silicon, and find that silicon with 20 nm nanograins can result in more than five times enhancement in their thermoelectric figure of merit as the grain boundaries scatter phonons more significantly than that of electrons due to their disparate MFP distributions.

preprint2012arXiv

Tunable Thermal Transport and Thermal Rectification in Strained Graphene Nanoribbons

Using molecular dynamics(MD) simulations, we study thermal transport in graphene nanoribbons (GNR) subjected to uniform uniaxial and nonuniform strain fields. We predict significant thermal rectification (over 70%) in a rectangular armchair GNR by applying a transverse force asymmetrically. The heat flux is larger from the less stressed region to the more stressed region. Furthermore, we develop a theoretical framework based on the non-equilibrium thermodynamics to discuss when thermal rectification under a stress gradient can occur. We conclude with a discussion of details relevant to experiment.

preprint2011arXiv

Molecular dynamics simulations of thermal conductivity and spectral phonon relaxation time in suspended and supported graphene

We perform molecular dynamics (MD) simulations with phonon spectral analysis aiming at understanding the two dimensional (2D) thermal transport in suspended and supported graphene. Within the framework of equilibrium MD simulations, we perform spectral energy density (SED) analysis to obtain the lifetime of individual phonon modes. The per-mode contribution to thermal conductivity is then calculated to obtain the lattice thermal conductivity in the temperature range 300-650 K. In contrast to prior studies, our results suggest that the contribution from out-of-plane acoustic (or ZA) branch to thermal conductivity is around 25-30% in suspended single-layer graphene (SLG) at room temperature. The thermal conductivity is found to reduce when SLG is put on amorphous SiO2 substrate. Such reduction is attributed to the strengthened scattering in all phonon modes in the presence of the substrate. Among them, ZA modes are mostly affected with their contribution to thermal conductivity reduced to around 15%. As a result, thermal transport is dominated by in-plane acoustic phonon modes in supported SLG.

preprint2011arXiv

Nonlinear thermal transport and negative differential thermal conductance in graphene nanoribbons

We employ classical molecular dynamics to study the nonlinear thermal transport in graphene nanoribbons (GNRs). For GNRs under large temperature biases beyond linear response regime, we have observed the onset of negative differential thermal conductance (NDTC). NDTC is tunable by varying the manner of applying the temperature biases. NDTC is reduced and eventually disappears when the length of the GNR increases. We have also observed NDTC in triangular GNRs, where NDTC exists only when the heat current is from the narrower to the wider end. These effects may be useful in nanoscale thermal managements and thermal signal processing utilizing GNRs.

preprint2010arXiv

Thermal Conductivity and Thermal Rectification in Graphene Nanoribbons: a Molecular Dynamics Study

We have used molecular dynamics to calculate the thermal conductivity of symmetric and asymmetric graphene nanoribbons (GNRs) of several nanometers in size (up to ~4 nm wide and ~10 nm long). For symmetric nanoribbons, the calculated thermal conductivity (e.g. ~2000 W/m-K @400K for a 1.5 nm {\times} 5.7 nm zigzag GNR) is on the similar order of magnitude of the experimentally measured value for graphene. We have investigated the effects of edge chirality and found that nanoribbons with zigzag edges have appreciably larger thermal conductivity than nanoribbons with armchair edges. For asymmetric nanoribbons, we have found significant thermal rectification. Among various triangularly-shaped GNRs we investigated, the GNR with armchair bottom edge and a vertex angle of 30° gives the maximal thermal rectification. We also studied the effect of defects and found that vacancies and edge roughness in the nanoribbons can significantly decrease the thermal conductivity. However, substantial thermal rectification is observed even in the presence of edge roughness.

preprint2010arXiv

Tuning the thermal conductivity of graphene nanoribbons by edge passivation and isotope engineering: a molecular dynamics study

Using classical molecular dynamics simulation, we have studied the effect of edge-passivation by hydrogen (H-passivation) and isotope mixture (with random or supperlattice distributions) on the thermal conductivity of rectangular graphene nanoribbons (GNRs) (of several nanometers in size). We found that the thermal conductivity is considerably reduced by the edge H-passivation. We also find that the isotope mixing can reduce the thermal conductivities, with the supperlattice distribution giving rise to more reduction than the random distribution. These results can be useful in nanoscale engineering of thermal transport and heat management using GNRs.