Showing 1–42 of 42 results for author: Ruan, X

Sampling-accelerated First-principles Prediction of Phonon Scattering Rates for Converged Thermal Conductivity and Radiative Properties

Authors: Ziqi Guo, Zherui Han, Dudong Feng, Guang Lin, Xiulin Ruan

Abstract: First-principles prediction of thermal conductivity and radiative properties is crucial. However, computing phonon scattering, especially for four-phonon scattering, could be prohibitively expensive, and the thermal conductivity even for silicon was still under-predicted and not converged in the literature. Here we propose a method to estimate scattering rates from a small sample of scattering pro… ▽ More First-principles prediction of thermal conductivity and radiative properties is crucial. However, computing phonon scattering, especially for four-phonon scattering, could be prohibitively expensive, and the thermal conductivity even for silicon was still under-predicted and not converged in the literature. Here we propose a method to estimate scattering rates from a small sample of scattering processes using maximum likelihood estimation. The computational cost of estimating scattering rates and associated thermal conductivity and radiative properties is dramatically reduced by over 99%. This allows us to use an unprecedented q-mesh of 32*32*32 for silicon and achieve a converged thermal conductivity value that agrees much better with experiments. The accuracy and efficiency of our approach make it ideal for the high-throughput screening of materials for thermal and optical applications. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Four phonon-dominated near-field radiation in weakly anharmonic polar materials

Authors: Dudong Feng, Xiaolong Yang, Zherui Han, Xiulin Ruan

Abstract: Inelastic scattering processes typically introduce friction among carriers and reduce the transport properties of photons, phonons, and electrons. However, we predict that in contrast to the role in reducing thermal conductivity, four-phonon scattering dominates near-field radiative heat transfer (NFRHT) in both boron arsenide~(BAs) and boron antimonide. Including four-phonon scattering results in… ▽ More Inelastic scattering processes typically introduce friction among carriers and reduce the transport properties of photons, phonons, and electrons. However, we predict that in contrast to the role in reducing thermal conductivity, four-phonon scattering dominates near-field radiative heat transfer (NFRHT) in both boron arsenide~(BAs) and boron antimonide. Including four-phonon scattering results in a nearly 400-fold increase in the total heat flux between two BAs thin-films compared to three-phonon scattering alone. This non-intuitive enhancement arises from the large number of NFRHT channels activated by four-phonon scattering outcompete the effect of decreased coupling strength of surface phonon polaritons at the resonance frequency. Additionally, we point out that four-phonon scattering to decrease NFRHT in certain other systems. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Predictions and Measurements of Thermal Conductivity of Ceramic Materials at High Temperature

Authors: Zherui Han, Zixin Xiong, William T. Riffe, Hunter B. Schonfeld, Mauricio Segovia, Jiawei Song, Haiyan Wang, Xianfan Xu, Patrick E. Hopkins, Amy Marconnet, Xiulin Ruan

Abstract: The lattice thermal conductivity ($κかっぱ$) of two ceramic materials, cerium dioxide (CeO$_2$) and magnesium oxide (MgO), is computed up to 1500 K using first principles and the phonon Boltzmann Transport Equation (PBTE) and compared to time-domain thermoreflectance (TDTR) measurements up to 800 K. Phonon renormalization and the four-phonon effect, along with high temperature thermal expansion, are int… ▽ More The lattice thermal conductivity ($κかっぱ$) of two ceramic materials, cerium dioxide (CeO$_2$) and magnesium oxide (MgO), is computed up to 1500 K using first principles and the phonon Boltzmann Transport Equation (PBTE) and compared to time-domain thermoreflectance (TDTR) measurements up to 800 K. Phonon renormalization and the four-phonon effect, along with high temperature thermal expansion, are integrated in our \textit{ab initio} molecular dynamics (AIMD) calculations. This is done by first relaxing structures and then fitting to a set of effective force constants employed in a temperature-dependent effective potential (TDEP) method. Both three-phonon and four-phonon scattering rates are computed based on these effective force constants. Our calculated thermal conductivities from the PBTE solver agree well with literature and our TDTR measurements. Other predicted thermal properties including thermal expansion, frequency shift, and phonon linewidth also compare well with available experimental data. Our results show that high temperature softens phonon frequency and reduces four-phonon scattering strength in both ceramics. Compared to MgO, we find that CeO$_2$ has weaker four-phonon effect and renormalization greatly reduces its four-phonon scattering rates. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Journal ref: Phys. Rev. B 108, 184306 (2023)

Is Thermal Conductivity of Graphene Divergent and Higher Than Diamond?

Authors: Zherui Han, Xiulin Ruan

Abstract: The thermal conductivity of monolayer graphene is an outstanding challenge with no consensus reached on its exact value and length convergence so far. We consider four-phonon scattering, phonon renormalization, and an exact solution to phonon Boltzmann transport equation (BTE) from first principles. Using this computational formalism with unprecedented sampling grid, we show that when four-phonon… ▽ More The thermal conductivity of monolayer graphene is an outstanding challenge with no consensus reached on its exact value and length convergence so far. We consider four-phonon scattering, phonon renormalization, and an exact solution to phonon Boltzmann transport equation (BTE) from first principles. Using this computational formalism with unprecedented sampling grid, we show that when four-phonon scattering is included the thermal conductivity is convergent with system size at a room temperature value of 1300 W/(m$\cdot$K), which is lower than that of diamond. On the contrary, considering three-phonon scattering only yields divergence with size due to the momentum-conserving normal processes of flexural phonons. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Temperature-Dependent Full Spectrum Optical Responses of Semiconductors from First Principles

Authors: Zherui Han, Changkyun Lee, Jiawei Song, Haiyan Wang, Peter Bermel, Xiulin Ruan

Abstract: From ultraviolet to mid-infrared region, light-matter interaction mechanisms in semiconductors progressively shift from electronic transitions to phononic resonances and are affected by temperature. Here, we present a parallel temperature-dependent treatment of both electrons and phonons entirely from first principles, enabling the prediction of full-spectrum optical responses. At elevated tempera… ▽ More From ultraviolet to mid-infrared region, light-matter interaction mechanisms in semiconductors progressively shift from electronic transitions to phononic resonances and are affected by temperature. Here, we present a parallel temperature-dependent treatment of both electrons and phonons entirely from first principles, enabling the prediction of full-spectrum optical responses. At elevated temperatures, $\textit{ab initio}$ molecular dynamics is employed to find thermal perturbations to electronic structures and construct effective force constants describing potential landscape. Four-phonon scattering and phonon renormalization are included in an integrated manner in this approach. As a prototype ceramic material, cerium dioxide (CeO$_2$) is considered in this work. Our first-principles calculated refractive index of CeO$_2$ agrees well with measured data from literature and our own temperature-dependent ellipsometer experiment. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Effects of hot phonons and thermal stress in micro-Raman spectra of molybdenum disulphide

Authors: Peter Sokalski, Zherui Han, Gabriella Coloyan Fleming, Brandon Smith, Sean E. Sullivan, Rui Huang, Xiulin Ruan, Li Shi

Abstract: Micro-Raman spectroscopy has become an important tool in probing thermophysical behavior in emerging functional materials such as two-dimensional (2D) layered structures. Localized heating by the focused Raman excitation laser beam is expected to produce both stress and nonequilibrium temperature distributions in the material. Here we investigate the effects of hot optical phonons in the Raman spe… ▽ More Micro-Raman spectroscopy has become an important tool in probing thermophysical behavior in emerging functional materials such as two-dimensional (2D) layered structures. Localized heating by the focused Raman excitation laser beam is expected to produce both stress and nonequilibrium temperature distributions in the material. Here we investigate the effects of hot optical phonons in the Raman spectra of molybdenum disulphide (MoS2) and distinguish them from those caused by thermally-induced compressive stress, which causes a Raman frequency blue shift in the focused Raman laser spot. With the stress effect accounted for in micro-Raman measurements, the degree of nonequilibrium between the hot optical phonons and the acoustic phonon bath is found to be much smaller than those observed in prior micro-Raman measurements of suspended graphene. The observation agrees with a first-principles based theoretical prediction of overpopulated zone-center optical phonons compared to other optical phonons in the Brillouin zone and acoustic phonons. The findings provide detailed insight into the energy relaxation processes in this emerging 2D transition metal dichalcogenide (TMD) electronic and optoelectronic material and clarify an important question in micro-Raman measurements of thermal transport in this and other 2D materials. △ Less

Submitted 9 November, 2022; v1 submitted 30 June, 2022; originally announced June 2022.

Comments: 19 pages. 4 figures

Prediction of Hot Zone-center Optical Phonons in Laser Irradiated Molybdenum Disulfide with a Semiconductor Multitemperature Model

Authors: Zherui Han, Peter Sokalski, Li Shi, Xiulin Ruan

Abstract: Previous multitemperature model (MTM) resolving phonon temperatures at the polarization level and measurements have uncovered remarkable nonequilibrium among different phonon polarizations in laser irradiated graphene and metals. Here, we develop a semiconductor-specific MTM (SC-MTM) by including electron-hole pair generation, diffusion, and recombination, and show that a phonon polarization-level… ▽ More Previous multitemperature model (MTM) resolving phonon temperatures at the polarization level and measurements have uncovered remarkable nonequilibrium among different phonon polarizations in laser irradiated graphene and metals. Here, we develop a semiconductor-specific MTM (SC-MTM) by including electron-hole pair generation, diffusion, and recombination, and show that a phonon polarization-level model does not yield observable polarization-based nonequilibrium in laser-irradiated molybdenum disulfide (MoS$_2$). In contrast, appreciable nonequilibrium is predicted between zone-center optical phonons and the other modes. The momentum-based nonequilibrium ratio is found to increase with decreasing laser spot size and interaction with a substrate. This finding is relevant to the understanding of the energy relaxation process in two-dimensional optoelectronic devices and Raman measurements of thermal transport. △ Less

Submitted 10 January, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

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

Authors: Zherui Han, Xiaolong Yang, Sean E. Sullivan, Tianli Feng, Li Shi, Wu Li, Xiulin Ruan

Abstract: 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 co… ▽ More 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). △ Less

Submitted 1 July, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Journal ref: Phys. Rev. Lett. 128, 045901, 2022

FourPhonon: An extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity

Authors: Zherui Han, Xiaolong Yang, Wu Li, Tianli Feng, Xiulin Ruan

Abstract: FourPhonon is a computational package that can calculate four-phonon scattering rates in crystals. It is built within ShengBTE framework, which is a well-recognized lattice thermal conductivity solver based on Boltzmann transport equation. An adaptive energy broadening scheme is implemented for the calculation of four-phonon scattering rates. In analogy with $thirdorder.py$ in ShengBTE, we also pr… ▽ More FourPhonon is a computational package that can calculate four-phonon scattering rates in crystals. It is built within ShengBTE framework, which is a well-recognized lattice thermal conductivity solver based on Boltzmann transport equation. An adaptive energy broadening scheme is implemented for the calculation of four-phonon scattering rates. In analogy with $thirdorder.py$ in ShengBTE, we also provide a separate python script, $Fourthorder.py$, to calculate fourth-order interatomic force-constants. The extension module preserves all the nice features of the well-recognized lattice thermal conductivity solver ShengBTE, including good parallelism and straightforward workflow. In this paper, we discuss the general theory, program design, and example calculations on Si, BAs and $\mathrm{LiCoO_2}$. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Journal ref: Comput Phys Commun 270 (2022) 108179

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

Authors: Zhen Tong, Joseph Peoples, Xiangyu Li, Xiaolong Yang, Hua Bao, Xiulin Ruan

Abstract: 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… ▽ More 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. △ Less

Submitted 21 January, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 5 pages, 4 figures

An iterative machine learning approach for discovering unexpected thermal conductivity enhancement in aperiodic superlattices

Authors: Prabudhya Roy Chowdhury, Xiulin Ruan

Abstract: While machine learning (ML) has shown increasing effectiveness in optimizing materials properties under known physics, its application in challenging conventional wisdom and discovering new physics still remains challenging due to its interpolative nature. In this work, we demonstrate the potential of using ML for such applications by implementing an adaptive ML-accelerated search process that can… ▽ More While machine learning (ML) has shown increasing effectiveness in optimizing materials properties under known physics, its application in challenging conventional wisdom and discovering new physics still remains challenging due to its interpolative nature. In this work, we demonstrate the potential of using ML for such applications by implementing an adaptive ML-accelerated search process that can discover unexpected lattice thermal conductivity ($κかっぱ_l$) enhancement instead of reduction in aperiodic superlattices (SLs) as compared to periodic superlattices. We use non-equilibrium molecular dynamics (NEMD) simulations for high-fidelity calculations of $κかっぱ_l$ for a small fraction of SLs in the search space, along with a convolutional neural network (CNN) which can rapidly predict $κかっぱ_l$ for a large number of structures. To ensure accurate prediction by the CNN for the target unknown structures, we iteratively identify aperiodic SLs containing structural features which lead to locally enhanced thermal transport, and include them as additional training data for the CNN in each iteration. As a result, our CNN can accurately predict the high $κかっぱ_l$ of aperiodic SLs that are absent from the initial training dataset, which allows us to identify the previously unseen exceptional structures. The identified RML structures exhibit increased coherent phonon contribution to thermal conductivity owing to the presence of closely spaced interfaces. Our work describes a general purpose machine learning approach for identifying low-probability-of-occurrence exceptional solutions within an extremely large subspace and discovering the underlying physics. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Quantum versus Classical Regime in Circuit Quantum Acoustodynamics

Authors: Gang-hui Zeng, Yang Zhang, Aleksey N. Bolgar, Dong He, Bin Li, Xin-hui Ruan, Lan Zhou, Le-Mang Kuang, Oleg V. Astafiev, Yu-xi Liu, Z. H. Peng

Abstract: We experimentally study a circuit quantum acoustodynamics system, which consists of a superconducting artificial atom, coupled to both a two-dimensional surface acoustic wave resonator and a one-dimensional microwave transmission line. The strong coupling between the artificial atom and the acoustic wave resonator is confirmed by the observation of the vacuum Rabi splitting at the base temperature… ▽ More We experimentally study a circuit quantum acoustodynamics system, which consists of a superconducting artificial atom, coupled to both a two-dimensional surface acoustic wave resonator and a one-dimensional microwave transmission line. The strong coupling between the artificial atom and the acoustic wave resonator is confirmed by the observation of the vacuum Rabi splitting at the base temperature of dilution refrigerator. We show that the propagation of microwave photons in the microwave transmission line can be controlled by a few phonons in the acoustic wave resonator. Furthermore, we demonstrate the temperature effect on the measurements of the Rabi splitting and temperature induced transitions from high excited dressed states. We find that the spectrum structure of two-peak for the Rabi splitting becomes into those of several peaks, and gradually disappears with the increase of the environmental temperature $T$. The quantum-to-classical transition is observed around the crossover temperature $T_{c}$, which is determined via the thermal fluctuation energy $k_{B}T$ and the characteristic energy level spacing of the coupled system. Experimental results agree well with the theoretical simulations via the master equation of the coupled system at different effective temperatures. △ Less

Submitted 18 January, 2021; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: 14 pages, 6 figures

Anisotropic ultrafast spin/valley dynamics in WTe2 films

Authors: Yequan Chen, Zhendong Chen, Yongda Chen, Liming Chen, Jiai Ning, Bo Liu, Chunchen Zhang, Xianyang Lu, Xuezhong Ruan, Wenqing Liu, Peng Wang, Fengqi Song, Chunfeng Zhang, Fengqiu Wang, Jing Wu, Liang He, Xuefeng Wang, Rong Zhang, Yongbing Xu

Abstract: WTe2 Weyl semimetal hosts the natural broken inversion symmetry and strong spin orbit coupling, making it promising for exotic spin/valley dynamics within a picosecond timescale. Here, we unveil an anisotropic ultrafast spin/valley dynamics in centimeter-scale, single-crystalline Td-WTe2 films using a femtosecond pump-probe technique at room temperature. We observe a transient (~0.8 ps) intra-vall… ▽ More WTe2 Weyl semimetal hosts the natural broken inversion symmetry and strong spin orbit coupling, making it promising for exotic spin/valley dynamics within a picosecond timescale. Here, we unveil an anisotropic ultrafast spin/valley dynamics in centimeter-scale, single-crystalline Td-WTe2 films using a femtosecond pump-probe technique at room temperature. We observe a transient (~0.8 ps) intra-valley transition and a subsequent polarization duration (~5 ps) during the whole spin/valley relaxation process. Furthermore, the relaxation exhibits the remarkable anisotropy of approximately six-fold and two-fold symmetries due to the intrinsic anisotropy along the crystalline orientation and the extrinsic matrix element effect, respectively. Our results offer a prospect for the ultrafast manipulation of spin/valleytronics in topological quantum materials for dissipationless high-speed spin/valleytronic devices. △ Less

Submitted 3 September, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Comments: 21 pages, 4 figures

Journal ref: J. Appl. Phys. 131, 163903 (2022)

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

Authors: Evan Witkoske, Zhen Tong, Yining Feng, Xiulin Ruan, Mark Lundstrom, Na Lu

Abstract: 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 therma… ▽ More 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. △ Less

Submitted 13 December, 2019; originally announced December 2019.

Computing ground states of Bose-Einstein Condensates with higher order interaction via a regularized density function formulation

Authors: Weizhu Bao, Xinran Ruan

Abstract: We propose and analyze a new numerical method for computing the ground state of the modified Gross-Pitaevskii equation for modeling the Bose-Einstein condensate with a higher order interaction by adapting the density function formulation and the accelerated projected gradient method. By reformulating the energy functional $E(φふぁい)$ with $φふぁい$, the wave function, in terms of the density $ρろー=|φふぁい|^2$, the o… ▽ More We propose and analyze a new numerical method for computing the ground state of the modified Gross-Pitaevskii equation for modeling the Bose-Einstein condensate with a higher order interaction by adapting the density function formulation and the accelerated projected gradient method. By reformulating the energy functional $E(φふぁい)$ with $φふぁい$, the wave function, in terms of the density $ρろー=|φふぁい|^2$, the original non-convex minimization problem for defining the ground state is then reformulated to a convex minimization problem. In order to overcome the semi-smoothness of the function $\sqrtρろー$ in the kinetic energy part, a regularization is introduced with a small parameter $0<\varepsilon\ll1$. Convergence of the regularization is established when $\varepsilon\to0$. The regularized convex optimization problem is discretized by the second order finite difference method. The convergence rates in terms of the density and energy of the discretization are established. The accelerated projected gradient method is adapted for solving the discretized optimization problem. Numerical results are reported to demonstrate the efficiency and accuracy of the proposed numerical method. Our results show that the proposed method is much more efficient than the existing methods in the literature, especially in the strong interaction regime. △ Less

Submitted 24 August, 2019; originally announced August 2019.

Comments: 26 pages, 5 figures

MSC Class: 35Q55; 65N06; 65N25; 90C30

Four-phonon scattering diminishes the optical phonon contribution and isotope effect to thermal conductivity of III-V semiconductors

Authors: Xiaolong Yang, Tianli Feng, Ju Li, Xiulin Ruan

Abstract: Recent studies reveal that four-phonon scattering is generally important in determining thermal conductivities of solids. However, these studies have been focused on materials where thermal conductivity $κかっぱ$ is dominated by acoustic phonons, and the impact of four phonon scattering, although significant, is still generally smaller than three-phonon scattering. In this work, taking AlSb as example,… ▽ More Recent studies reveal that four-phonon scattering is generally important in determining thermal conductivities of solids. However, these studies have been focused on materials where thermal conductivity $κかっぱ$ is dominated by acoustic phonons, and the impact of four phonon scattering, although significant, is still generally smaller than three-phonon scattering. In this work, taking AlSb as example, we demonstrated that four-phonon scattering is even more critical to three-phonon scattering as it diminishes optical phonon thermal transport, and therefore significantly reduces the thermal conductivities of materials in which optical branches have long three-phonon lifetimes. Also, our calculations show that four-phonon scattering can play an extremely important role in weakening the isotope effect on $κかっぱ$. Specifically, four-phonon scattering reduces the room-temperature $κかっぱ$ of the isotopically pure and natural-occurring AlSb by 70$\%$ and 50$\%$, respectively. The reduction for isotopically pure and natural-occurring c-GaN is about 34$\%$ and 27$\%$, respectively. For isotopically-pure w-GaN, the reduction is about 13$\%$ at room temperature and 25$\%$ at 400 K. These results provided important guidance for experimentalists for achieving high thermal conductivities in III-V compounds for applications in semiconductor industry. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. B 100, 245203 (2019)

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

Authors: Xiaolong Yang, Tianli Feng, Joon Sang Kang, Yongjie Hu, Ju Li, Xiulin Ruan

Abstract: 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 el… ▽ More 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. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. B 101, 161202 (2020)

A comprehensive first-principles analysis of phonon thermal conductivity and electron-phonon coupling in different metals

Authors: Zhen Tong, Shouhang Li, Xiulin Ruan, Hua Bao

Abstract: Separating electron and phonon thermal conductivity components is imperative for understanding the principle thermal transport mechanisms in metals and highly desirable in many applications. In this work, we predict the mode-dependent electron and phonon thermal conductivities of 18 different metals at room-temperature from first-principles. Our first-principles predictions, in general, agree well… ▽ More Separating electron and phonon thermal conductivity components is imperative for understanding the principle thermal transport mechanisms in metals and highly desirable in many applications. In this work, we predict the mode-dependent electron and phonon thermal conductivities of 18 different metals at room-temperature from first-principles. Our first-principles predictions, in general, agree well with experimental data. We find that the phonon thermal conductivity is in the range of 2 - 18 $W/mK$, which accounts for 1% - 40% of the total thermal conductivity. It is also found that the phonon thermal conductivities in transition metals and transition-intermetallic-compounds (TICs) are non-negligible compared to noble metals due to their high phonon group velocities. Besides, the electron-phonon coupling effect on phonon thermal conductivity in transition metals and intermetallic compounds is stronger than that of nobles, which is attributed to the larger electron-phonon coupling constant with a high electron density of state within Fermi window and high phonon frequency. The noble metals have higher electron thermal conductivities compared to transition metals and TICs, which is mainly due to the weak electron-phonon coupling in noble metals. It is also shown that the Lorenz ratios of transition metals and transition-intermetallic-compounds hold larger deviations from the Sommerfeld value $L_0=2.44 \times 10^{-8} W ΩおめがK^{-2}$. We also find the mean free paths (MFPs) for phonon (within 10 nm) are smaller than those of electron (5 - 25 nm). The electrical conductivity and electron thermal conductivity are strongly related to the MFPs of the electron. △ Less

Submitted 2 July, 2019; originally announced July 2019.

Comments: 23pages,10 figures, 2 tables

Journal ref: Phys. Rev. B 100, 144306 (2019)

Dramatic increase in the thermal boundary conductance and radiation limit from a Nonequilibrium Landauer Approach

Authors: Jingjing Shi, Xiaolong Yang, Timothy S. Fisher, Xiulin Ruan

Abstract: Thermal boundary conductance (TBC) is critical in many thermal and energy applications. A decades-old puzzle has been that many of the measured TBCs, such as those well characterized across Al/Si and ZnO/GaN interfaces, significantly exceed theoretical results or even the absolute upper limit called the ``radiation limit", suggesting the failure of the theory. Here, we identify that for high-trans… ▽ More Thermal boundary conductance (TBC) is critical in many thermal and energy applications. A decades-old puzzle has been that many of the measured TBCs, such as those well characterized across Al/Si and ZnO/GaN interfaces, significantly exceed theoretical results or even the absolute upper limit called the ``radiation limit", suggesting the failure of the theory. Here, we identify that for high-transmission interfaces, the commonly assumed phonon local thermal equilibrium adjacent to the interface fails, and the measurable phonon temperatures are not their emission temperature. We hence develop a ``nonequilibrium Landauer approach" and define the unique ``dressed" and ``intrinsic" TBCs. Combining our approach even with a simple diffuse mismatch model (DMM) nearly doubles the theoretical TBCs across the Al/Si and ZnO/GaN interfaces, and the theoretical results agree with experiments for the first time. The radiation limit is also redefined and found to increase over 100\% over the original radiation limit, and it can now well bound all the experimental data. △ Less

Submitted 14 April, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Phonon anharmonic frequency shift induced by four-phonon scattering calculated from first principles

Authors: Tianli Feng, Xiaolong Yang, Xiulin Ruan

Abstract: Phonon energies at finite temperatures shift away from their harmonic values due to anharmonicity. In this paper, we have realized the rigorous calculation of phonon energy shifts of silicon by three and four-phonon scattering from first principles. The anharmonic fourth-order force constants are calculated by considering up to the fifth nearest neighbors. The results agree reasonably well with av… ▽ More Phonon energies at finite temperatures shift away from their harmonic values due to anharmonicity. In this paper, we have realized the rigorous calculation of phonon energy shifts of silicon by three and four-phonon scattering from first principles. The anharmonic fourth-order force constants are calculated by considering up to the fifth nearest neighbors. The results agree reasonably well with available data from inelastic neutron scattering throughout the Brillouin zone. Surprisingly, the frequency shifts of optical phonon modes near the $Γがんま$ point are sensitive to the cutoff radius of the fourth-order force constants, in contrast to the four-phonon scattering rates, which nearly saturate when considering the second nearest neighbors. We have also compared the results with ab initio molecular dynamics simulations and found that the higher order of anharmonicity is important for optical phonons. Our work provides critical insight into the anharmonic phonon frequency shift and will have significant impact on the thermal and optical applications. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Enhancement of intrinsic magnetic damping in defect-free epitaxial Fe3O4 thin films

Authors: Xianyang Lu, Lewis J. Atkinson, Balati Kuerbanjiang, Bo Liu, Guanqi Li, Yichuan Wang, Junlin Wang, Xuezhong Ruan, Jing Wu, Richard F. L. Evans, Vlado K. Lazarov, Roy W. Chantrell, Yongbing Xu

Abstract: We have investigated the magnetic damping of precessional spin dynamics in defect-controlled epitaxial grown Fe$_3$O$_4$(111)/Yttria-stabilized Zirconia (YSZ) nanoscale films by all-optical pump-probe measurements. The intrinsic damping constant of the defect-free Fe$_3$O$_4$ film is found to be strikingly larger than that of the as-grown Fe$_3$O$_4$ film with structural defects. We demonstrate th… ▽ More We have investigated the magnetic damping of precessional spin dynamics in defect-controlled epitaxial grown Fe$_3$O$_4$(111)/Yttria-stabilized Zirconia (YSZ) nanoscale films by all-optical pump-probe measurements. The intrinsic damping constant of the defect-free Fe$_3$O$_4$ film is found to be strikingly larger than that of the as-grown Fe$_3$O$_4$ film with structural defects. We demonstrate that the population of the first-order perpendicular standing spin wave (PSSW) mode, which is exclusively observed in the defect-free film under sufficiently high external magnetic fields, leads to the enhancement of the magnetic damping of the uniform precession (Kittel) mode. We propose a physical picture in which the PSSW mode acts as an additional channel for the extra energy dissipation of the Kittel mode. The energy transfer from Kittel mode to PSSW mode increases as in-plane magnetization precession becomes more uniform, resulting in the unique intrinsic magnetic damping enhancement in the defect-free Fe$_3$O$_4$ film. △ Less

Submitted 16 May, 2019; v1 submitted 30 October, 2017; originally announced October 2017.

Comments: 5 pages, 4 figures

Journal ref: Appl. Phys. Lett. 114, 192406, 2019

Spectral analysis of non-equilibrium molecular dynamics: spectral phonon temperature and phonon local non-equilibrium in thin films and across interfaces

Authors: Tianli Feng, Wenjun Yao, Zuyuan Wang, Jingjing Shi, Chuang Li, Bingyang Cao, Xiulin Ruan

Abstract: Although extensive experimental and theoretical works have been conducted to understand the ballistic and diffusive phonon transport in nanomaterials recently, direct observation of temperature and thermal nonequilibrium of different phonon modes has not been realized. Herein, we have developed a method within the framework of molecular dynamics to calculate the temperatures of phonon in both real… ▽ More Although extensive experimental and theoretical works have been conducted to understand the ballistic and diffusive phonon transport in nanomaterials recently, direct observation of temperature and thermal nonequilibrium of different phonon modes has not been realized. Herein, we have developed a method within the framework of molecular dynamics to calculate the temperatures of phonon in both real and phase spaces. Taking silicon thin film and graphene as examples, we directly obtained the spectral phonon temperature (SPT) and observed the local thermal nonequilibrium between the ballistic and diffusive phonons. Such nonequilibrium also generally exists across interfaces and is surprisingly large, and it provides an additional thermal interfacial resistance mechanism. Our SPT results directly show that the vertical thermal transport across the dimensionally mismatched graphene/substrate interface is through the coupling between flexural acoustic phonons of graphene and the longitudinal phonons in the substrate with mode conversion. In the dimensionally matched interfaces, e.g. graphene/graphene junction and graphene/boron nitride planar interfaces, strong coupling occurs between the acoustic phonon modes on both sides, and the coupling decreases with interfacial mixing. The SPT method together with the spectral heat flux can eliminate the size effect of the thermal conductivity prediction induced from ballistic transport. Our work shows that in thin films and across interfaces, phonons are in local thermal nonequilibrium. △ Less

Submitted 31 March, 2017; originally announced March 2017.

Comments: 14 pages, 10 figures

Journal ref: Phys. Rev. B 95, 195202 (2017)

Compressive mechanical response of graphene foams and their thermal resistance with copper interfaces

Authors: Wonjun Park, Xiangyu Li, Nirajan Mandal, Xiulin Ruan, Yong P. Chen

Abstract: We report compressive mechanical response of graphene foams (GFs) and the thermal resistance ($R_{TIM}$) between copper (Cu) and GFs, where GFs were prepared by the chemical vapor deposition (CVD) method. We observe that Young's modulus ($E_{GF}$) and compressive strength ($σしぐま_{GF}$) of GFs have a power law dependence on increasing density ($ρろー_{GF}$) of GFs. The maximum efficiency of absorbed energ… ▽ More We report compressive mechanical response of graphene foams (GFs) and the thermal resistance ($R_{TIM}$) between copper (Cu) and GFs, where GFs were prepared by the chemical vapor deposition (CVD) method. We observe that Young's modulus ($E_{GF}$) and compressive strength ($σしぐま_{GF}$) of GFs have a power law dependence on increasing density ($ρろー_{GF}$) of GFs. The maximum efficiency of absorbed energy ($ηいーた_{max}$) for all GFs during the compression is larger than ~0.39. We also find that a GF with a higher $ρろー_{GF}$ shows a larger $ηいーた_{max}$. In addition, we observe that the measured $R_{TIM}$ of Cu/GFs at room temperature with a contact pressure of 0.25 MP applied increases from ~50 to ~90 $mm^2K/W$ when $ρろー_{GF}$ increases from 4.7 to 31.9 $mg/cm^3$. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Journal ref: APL Materials 5, 036102 (2017)

Optical generation and detection of local non-equilibrium phonons in suspended graphene

Authors: Sean Sullivan, Ajit Vallabhaneni, Iskandar Kholmanov, Xiulin Ruan, Jayathi Murthy, Li Shi

Abstract: The measured frequencies and intensities of different first- and second- order Raman peaks of suspended graphene are used to show that optical phonons and different acoustic phonon polarizations are driven out of local equilibrium inside a sub-micron laser spot. The experimental results are correlated with a first principles-based multiple temperature model to suggest a considerably lower equivale… ▽ More The measured frequencies and intensities of different first- and second- order Raman peaks of suspended graphene are used to show that optical phonons and different acoustic phonon polarizations are driven out of local equilibrium inside a sub-micron laser spot. The experimental results are correlated with a first principles-based multiple temperature model to suggest a considerably lower equivalent local temperature of the flexural phonons than those of other phonon polarizations. The finding reveals weak coupling between the flexural modes with hot electrons and optical phonons. Since the ultrahigh intrinsic thermal conductivity of graphene has been largely attributed to contributions from the flexural phonons, the observed local non-equilibrium phenomena have important implications for understanding energy dissipation processes in graphene-based electronic and optoelectronic devices, as well as in Raman measurements of thermal transport in graphene and other two-dimensional materials. △ Less

Submitted 11 January, 2017; originally announced January 2017.

Comments: 10 pages, 5 figures

Uncertainty quantification of thermal conductivities from equilibrium molecular dynamics simulations

Authors: Zuyuan Wang, Salar Safarkhani, Guang Lin, Xiulin Ruan

Abstract: 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,… ▽ More 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. △ Less

Submitted 19 December, 2016; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: Main article: 14 pages, 11 figures; Supplemental information: 4 pages, 3 figures

Thermal Transport at the Nanoscale - A Fourier's Law vs. Phonon Boltzmann Equation Study

Authors: Jan Kaiser, Tianli Feng, Jesse Maassen, Xufeng Wang, Xiulin Ruan, Mark Lundstrom

Abstract: Steady-state thermal transport in nanostructures with dimensions comparable to the phonon mean-free-path is examined. Both the case of contacts at different temperatures with no internal heat generation and contacts at the same temperature with internal heat generation are considered. Fourier's Law results are compared to finite volume method solutions of the phonon Boltzmann equation in the gray… ▽ More Steady-state thermal transport in nanostructures with dimensions comparable to the phonon mean-free-path is examined. Both the case of contacts at different temperatures with no internal heat generation and contacts at the same temperature with internal heat generation are considered. Fourier's Law results are compared to finite volume method solutions of the phonon Boltzmann equation in the gray approximation. When the boundary conditions are properly specified, results obtained using Fourier's Law without modifying the bulk thermal conductivity are in essentially exact quantitative agreement with the phonon Boltzmann equation in the ballistic and diffusive limits. The errors between these two limits are examined in this paper. For the four cases examined, the error in the apparent thermal conductivity as deduced from a correct application of Fourier's Law is less than 6%. We also find that the Fourier's Law results presented here are nearly identical to those obtained from a widely-used ballistic-diffusive approach, but analytically much simpler. Although limited to steady-state conditions with spatial variations in one dimension and to a gray model of phonon transport, the results show that Fourier's Law can be used for linear transport from the diffusive to the ballistic limit. The results also contribute to an understanding of how heat transport at the nanoscale can be understood in terms of the conceptual framework that has been established for electron transport at the nanoscale. △ Less

Submitted 9 December, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: 23 pages, 11 figures

Effects of the non-parabolic kinetic energy on non-equilibrium polariton condensates

Authors: Florian Pinsker, Xinran Ruan, Tristram Alexander

Abstract: In the study of non-equilibrium polariton condensates it is usually assumed that the dispersion relation of polaritons is parabolic in nature. We show that considering the true non-parabolic kinetic energy of polaritons leads to significant changes in the behaviour of the condensate due to the curvature of the dispersion relation and the possibility of transfer of energy to high wavenumber compone… ▽ More In the study of non-equilibrium polariton condensates it is usually assumed that the dispersion relation of polaritons is parabolic in nature. We show that considering the true non-parabolic kinetic energy of polaritons leads to significant changes in the behaviour of the condensate due to the curvature of the dispersion relation and the possibility of transfer of energy to high wavenumber components in the condensate spatial profile. We present explicit solutions for plane waves and linear excitations, and identify the differences in the theoretical predictions between the parabolic and non-parabolic mean-field models, showing the possibility of symmetry breaking in the latter. We then consider the evolution of wavepackets and show that self-localisation effects may be observed due to the curvature of the dispersion relation. Finally, we revisit the dynamics of dark soliton trains and show that additional localized density excitations may emerge in the dynamics due to the excitation of high frequency components, mimicking the appearance of near-bright solitary waves over short timescales. △ Less

Submitted 4 April, 2017; v1 submitted 7 June, 2016; originally announced June 2016.

Comments: to appear in Scientific reports

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

Authors: Tianli Feng, Bo Qiu, Xiulin Ruan

Abstract: 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-p… ▽ More 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}$. △ Less

Submitted 24 December, 2015; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: Published: Physical Review B 92, 235206 (2015)

Mean-field regime and Thomas-Fermi approximations of trapped Bose-Einstein condensates with higher order interactions in one and two dimensions

Authors: Xinran Ruan, Yongyong Cai, Weizhu Bao

Abstract: We derive rigorously one- and two-dimensional mean-field equations for cigar- and pancake-shaped Bose-Einstein condensates (BEC) with higher order interactions (HOI). We show how the higher order interaction modifies the contact interaction of the strongly confined particles. Surprisingly, we find that the usual Gaussian profile assumption for the strongly confining direction is inappropriate for… ▽ More We derive rigorously one- and two-dimensional mean-field equations for cigar- and pancake-shaped Bose-Einstein condensates (BEC) with higher order interactions (HOI). We show how the higher order interaction modifies the contact interaction of the strongly confined particles. Surprisingly, we find that the usual Gaussian profile assumption for the strongly confining direction is inappropriate for the cigar-shaped BEC case, and a Thomas-Fermi type profile should be adopted instead. Based on the derived mean field equations, the Thomas-Fermi densities are analyzed in presence of the contact interaction and HOI. For both box and harmonic traps in one, two and three dimensions, we identify the analytical Thomas-Fermi densities, which depend on the competition between the contact interaction and the HOI. △ Less

Submitted 31 October, 2015; originally announced November 2015.

Comments: 7 figures

Journal ref: J. Phys. B: At. Mol. Opt. Phys., 49 (2016), 125304

High-performance Thermal Interface Material Based on Few-layer Graphene Composite

Authors: Wonjun Park, Yufen Guo, Xiangyu Li, Jiuning Hu, Liwei Liu, Xiulin Ruan, Yong P. Chen

Abstract: We developed high-performance thermal interface materials (TIMs) based on few-layer graphene (FLG) composite, where FLG was prepared by the interlayer catalytic exfoliation (ICE) method. We experimentally demonstrated feasibility of FLG composites as TIMs by investigating their thermal and mechanical properties, and reliability. We measured the thermal interface resistance ($R_{int}$) between FLG… ▽ More We developed high-performance thermal interface materials (TIMs) based on few-layer graphene (FLG) composite, where FLG was prepared by the interlayer catalytic exfoliation (ICE) method. We experimentally demonstrated feasibility of FLG composites as TIMs by investigating their thermal and mechanical properties, and reliability. We measured the thermal interface resistance ($R_{int}$) between FLG composite TIMs (FLGTs) and copper and to be 3.2$\pm$1.7 and 4.3$\pm$1.4 $mm^2$K/W for 5 vol.% and 10 vol.% FLGTs at 330 K, respectively, comparable to or even lower than that of many commercial TIMs. In addition, the thermal conductivity ($κかっぱ_{TIM}$) of FLGTs is increased by an enhancement factor ($βべーた$) of ~17 as the FLG concentration increases from 0 to 10 vol.%. We also characterized Vickers hardness and glass transition temperature ($T_g$) of our FLGTs. We find that our FLGTs are thermally and mechanically reliable within practical operating temperature and pressure ranges. △ Less

Submitted 24 February, 2017; v1 submitted 31 October, 2015; originally announced November 2015.

Journal ref: J. Phys. Chem. C 119, 26753 (2015)

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

Authors: Tianli Feng, Xiulin Ruan

Abstract: 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 c… ▽ More 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. △ Less

Submitted 6 January, 2016; v1 submitted 2 October, 2015; originally announced October 2015.

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

Authors: Bo Qiu, Zhiting Tian, Ajit Vallabhaneni, Bolin Liao, Jonathan M. Mendoza, Oscar D. Restrepo, Xiulin Ruan, Gang Chen

Abstract: 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 MF… ▽ More 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. △ Less

Submitted 4 February, 2015; v1 submitted 16 September, 2014; originally announced September 2014.

Electrical and thermal conductivities of reduced graphene oxide/polystyrene composites

Authors: Wonjun Park, Jiuning Hu, Luis A. Jauregui, Xiulin Ruan, Yong P. Chen

Abstract: The author reports an experimental study of electrical and thermal transport in reduced graphene oxide (RGO)/polystyrene (PS) composites. The electrical conductivity ($σしぐま$) of RGO/PS composites with different RGO concentrations at room temperature shows a percolation behavior with the percolation threshold of ~ 0.25 vol.%. Their temperature-dependent electrical conductivity follows Efros-Shklovskii… ▽ More The author reports an experimental study of electrical and thermal transport in reduced graphene oxide (RGO)/polystyrene (PS) composites. The electrical conductivity ($σしぐま$) of RGO/PS composites with different RGO concentrations at room temperature shows a percolation behavior with the percolation threshold of ~ 0.25 vol.%. Their temperature-dependent electrical conductivity follows Efros-Shklovskii (ES) variable range hopping (VRH) conduction in the temperature range of 30 to 300 K. The thermal conductivity ($κかっぱ$) of composites is enhanced by ~ 90 % as the concentration is increased from 0 to 10 vol.%. The thermal conductivity of composites approximately linearly increases with increasing temperature from 150 to 300 K. Composites with a higher concentration show a stronger temperature dependence in the thermal conductivity. △ Less

Submitted 24 February, 2017; v1 submitted 10 March, 2014; originally announced March 2014.

Journal ref: Appl. Phys. Lett. 104, 113101 (2014)

Graphene/liquid crystal based terahertz phase shifters

Authors: Yang Wu, Xuezhong Ruan, Chih-Hsin Chen, Young Jun Shin, Youngbin Lee, Jing Niu, Jingbo Liu, Yuanfu Chen, Kun-Lin Yang, Xinhai Zhang, Jong-Hyun Ahn, Hyunsoo Yang

Abstract: Due to its high electrical conductivity and excellent transmittance at terahertz frequencies, graphene is a promising candidate as transparent electrodes for terahertz devices. We demonstrate a liquid crystal based terahertz phase shifter with the graphene films as transparent electrodes. The maximum phase shift is 10.8 degree and the saturation voltage is 5 V with a 50 um liquid crystal cell. The… ▽ More Due to its high electrical conductivity and excellent transmittance at terahertz frequencies, graphene is a promising candidate as transparent electrodes for terahertz devices. We demonstrate a liquid crystal based terahertz phase shifter with the graphene films as transparent electrodes. The maximum phase shift is 10.8 degree and the saturation voltage is 5 V with a 50 um liquid crystal cell. The transmittance at terahertz frequencies and electrical conductivity depending on the number of graphene layer are also investigated. The proposed phase shifter provides a continuous tunability, fully electrical controllability, and low DC voltage operation. △ Less

Submitted 6 September, 2013; originally announced September 2013.

Journal ref: Optics Express 21, 21395 (2013)

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

Authors: Bo Qiu, Xiulin Ruan

Abstract: 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 calculat… ▽ More 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. △ Less

Submitted 20 November, 2011; originally announced November 2011.

Journal ref: Appl. Phys. Lett. 100, 193101 (2012)

Nonlinear thermal transport and negative differential thermal conductance in graphene nanoribbons

Authors: Jiuning Hu, Yan Wang, Ajit Vallabhaneni, Xiulin Ruan, Yong P. Chen

Abstract: 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 th… ▽ More 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. △ Less

Submitted 28 September, 2011; v1 submitted 5 April, 2011; originally announced April 2011.

Comments: 5 pages, 4 figures

Journal ref: Applied Physics Letters 99 ,113101 (2011)

Tunable Thermal Transport and Thermal Rectification in Strained Graphene Nanoribbons

Authors: K. G. S. H. Gunawardana, Kieran Mullen, Jiuning Hu, Yong P. Chen, Xiulin Ruan

Abstract: 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… ▽ More 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. △ Less

Submitted 22 May, 2012; v1 submitted 12 November, 2010; originally announced November 2010.

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

Authors: Jiuning Hu, Xiulin Ruan, Yong P. Chen

Abstract: 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 gr… ▽ More 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. △ Less

Submitted 6 August, 2010; originally announced August 2010.

Comments: 13 pages, 5 figures, slightly expanded from the published version on Nano Lett. with some additional notes

Journal ref: Jiuning Hu, Xiulin Ruan, Yong P. Chen, Nano Lett. 9 (2009) 2730-2735

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

Authors: Jiuning Hu, Stephen Schiffli, Ajit Vallabhaneni, Xiulin Ruan, Yong P. Chen

Abstract: 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 t… ▽ More 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. △ Less

Submitted 6 August, 2010; originally announced August 2010.

Comments: 4 pages, 4 figures

Spin dynamics of electrons in the first excited subband of a high-mobility low-density 2D electron system

Authors: Haihui Luo, Xuan Qian, Xuezhong Ruan, Yang Ji

Abstract: We report on time-resolved Kerr rotation measurements of spin coherence of electrons in the first excited subband of a high-mobility low-density two-dimensional electron system in a GaAs/Al0.35Ga0.65As heterostructure. While the transverse spin lifetime (T2*) of electrons decreases monotonically with increasing magnetic field, it has a non-monotonic dependence on the temperature, with a peak val… ▽ More We report on time-resolved Kerr rotation measurements of spin coherence of electrons in the first excited subband of a high-mobility low-density two-dimensional electron system in a GaAs/Al0.35Ga0.65As heterostructure. While the transverse spin lifetime (T2*) of electrons decreases monotonically with increasing magnetic field, it has a non-monotonic dependence on the temperature, with a peak value of 596 ps at 36 K, indicating the effect of inter-subband electron-electron scattering on the electron spin relaxation. The spin lifetime may be long enough for potential device application with electrons in excited subbands. △ Less

Submitted 23 February, 2009; originally announced February 2009.

Effect of electron-electron scattering on spin dephasing in a high-mobility low-density twodimensional electron gas

Authors: X. Z. Ruan, H. H. Luo, Yang Ji, Z. Y. Xu, V. Umansky

Abstract: Utilizing time-resolved Kerr rotation techniques, we have investigated the spin dynamics of a high mobility, low density two dimensional electron gas in a GaAs/Al0:35Ga0:65As heterostructure in dependence on temperature from 1.5 K to 30 K. It is found that the spin relaxation/dephasing time under a magnetic field of 0.5 T exhibits a maximum of 3.12 ns around 14 K, superimposed on an increasing b… ▽ More Utilizing time-resolved Kerr rotation techniques, we have investigated the spin dynamics of a high mobility, low density two dimensional electron gas in a GaAs/Al0:35Ga0:65As heterostructure in dependence on temperature from 1.5 K to 30 K. It is found that the spin relaxation/dephasing time under a magnetic field of 0.5 T exhibits a maximum of 3.12 ns around 14 K, superimposed on an increasing background with rising temperature. The appearance of the maximum is ascribed to that at the temperature where the crossover from the degenerate to the nondegenerate regime takes place, electron-electron Coulomb scattering becomes strongest, and thus inhomogeneous precession broadening due to D'yakonov-Perel'(DP) mechanism becomes weakest. These results agree with the recent theoretical predictions [Zhou et al., PRB 75, 045305 (2007)], verifying the importance of electron-electron Coulomb scattering to electron spin relaxation/dephasing. △ Less

Submitted 17 April, 2008; originally announced April 2008.

Comments: 4 pages, 2 figures

Evidence for the intrinsic ferromagnetism of semiconductor (Ga,Cr)As revealed by magnetic circular dichroism

Authors: H. D. Gan, H. Z. Zheng, J. F. Bi, Y. Ji, B. Q. Sun, G. R. Li, X. Z. Ruan, J. Lu, J. H. Zhao

Abstract: To clarify whether or not (Ga,Cr)As is an intrinsic diluted magnetic semiconductor, a systematic study of the magnetic circular dichroism (MCD) was carried out for a series of (Ga,Cr)As epilayers grown by the low-temperature molecular-beam epitaxy technique. The present work provides unambiguous evidence for the intrinsic ferromagnetism of the (Ga,Cr)As epilayers with all the necessary propertie… ▽ More To clarify whether or not (Ga,Cr)As is an intrinsic diluted magnetic semiconductor, a systematic study of the magnetic circular dichroism (MCD) was carried out for a series of (Ga,Cr)As epilayers grown by the low-temperature molecular-beam epitaxy technique. The present work provides unambiguous evidence for the intrinsic ferromagnetism of the (Ga,Cr)As epilayers with all the necessary properties in accordance with that of a diluted magnetic semiconductor, especially the hysteresis characteristics, which is opened up in the magnetic field dependence of MCD. △ Less

Submitted 1 November, 2007; originally announced November 2007.

Comments: 13 pages