BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Philip B. Allen

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallcond-mat.othercond-mat.dis-nncond-mat.stat-mechcond-mat.str-elmath-phmath.MPphysics.comp-ph

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Non-local thermal transport modeling using the thermal distributor

Thermal transport in a quasi-ballistic regime is determined not only by the local temperature $T(r)$, or its gradient $\nabla T(r)$, but also by temperature distribution at neighboring points. For an accurate description of non-local effects on thermal transport, we employ the thermal distributor, $Θ(r,r')$, which provides the temperature response of the system at point $r$ to the heat input at point $r'$. We determine the thermal distributors from the linearized Peierls-Boltzmann equation (LPBE) and the relaxation time approximation (RTA) of the Peierls-Boltzmann equation and employ them to describe thermal transport in quasi-ballistic graphene devices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlocal Phonon Heat Transport Seen in 1-d Pulses

Phonons are the main heat carriers in semiconductor devices. In small devices, heat is not driven by a local temperature gradient, but by local points of heat input and removal. This complicates theoretical modeling. Study of the propagation of vibrational energy from an initial localized pulse provides insight into nonlocal phonon heat transport. We report simulations of pulse propagation in one dimension. The 1d case has tricky anomalies, but provides the simplest pictures of the evolution from initially ballistic toward longer time diffusive propagation. Our results show surprising details, such as diverse results from different definitions of atomistic local energy, and failure to exhibit pure diffusion at long times. Boltzmann phonon gas theory, including external energy insertion, is applied to this inherentlytime-dependent and nonlocal problem. The solution, using relaxation time approximation for impurity scattering, does not closely agree with the simulated results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accidental Degeneracy in k-space, Geometrical Phase, and the Perturbation of $π$ by Spin-orbit Interactions

Since closed lines of {\it accidental} electronic degeneracies were demonstrated to be possible, even frequent, by Herring in 1937, no further developments arose for eight decades. The earliest report of such a nodal loop in a real material -- aluminum -- is recounted and elaborated on. Nodal loop semimetals have become a focus of recent activity, with emphasis on other issues. Band degeneracies are, after all, the origin of topological phases in crystalline materials. Spin-orbit interaction lifts accidental band degeneracies, with the resulting spectrum being provided here. The geometric phase $γ(C)=\pmπ$ for circuits $C$ surrounding a line of such degeneracy cannot survive completely unchanged. The change depends on how the spin is fixed during adiabatic evolution. For spin fixed along the internal spin-orbit field, $γ(C)$ decreases to zero as the circuit collapses around the line of lifted degeneracy. For spin fixed along a perpendicular axis, the conical intersection persists and $γ(C)=\pmπ$ is unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First-principles study of pyroelectricity in GaN and ZnO

First-principles calculations are made for the primary pyroelectric coefficients of wurtzite GaN and ZnO. The pyroelectricity is attributed to the quasiharmonic thermal shifts of internal strains (internal displacements of cations and anions carrying their Born effective charges). The primary (zero-external-strain) pyroelectricity dominates at low temperatures, while the secondary pyroelectricity (the correction from external thermal strains) becomes comparable with the primary pyroelectricity at high temperatures. Contributions from the acoustic and the optical phonon modes to the primary pyroelectric coefficient are only moderately well described by the corresponding Debye function and Einstein function respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of Fröhlich polaron coupling on renormalized electron bands in polar semiconductors. Results for zincblende GaN

\ni We develop a simple method to study the zero-point and thermally renormalized electron energy $\varepsilon_{\mathbf{k}n}(T)$ for $\mathbf{k}n$ the conduction band minimum or valence maximum in polar semiconductors. We use the adiabatic approximation, including an imaginary broadening parameter $iδ$ to supress noise in the density-functional integrations. Fröhlich polaron methods provide analytic expressions for the contribution of the problematic optical phonon mode. We use this to correct the renormalization obtained from the adiabatic approximation. Test calculations are done for zincblende GaN for an 18x18x18 integration grid. The Fröhlich correction is of order -0.02 eV for the zero-point energy shift of the conduction band minimum, and +0.03 eV for the valence band maximum; the correction to renormalization of the 3.28 eV gap is -0.05 eV, a significant fraction of the total zero point renormalization of -0.15 eV.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phonon thermal conductivity by non-local non-equilibrium molecular dynamics

Non-equilibrium (NE) molecular dynamics (MD), or NEMD, gives a "direct" simulation of thermal conductivity kappa. Heat H(x) is added and subtracted in equal amounts at different places x. After steady state is achieved, the temperature T(x) is found by averaging over finite sections. Usually the aim is to extract a value of dT/dx from a place distant from sources and sinks of heat. This yields an effective kappa(L) for the thermal conductivity, L being the system size. The result is then studied as a function of L, to extract the bulk limit kappa. Here instead, our heat is H(x)~sin(qx), where q=2pi/L. This causes a steady-state temperature T_0 + Delta T sin(2pi x/L). A thermal conductivity kappa(q) is extracted, which is well converged at the chosen q (or L). Bulk conductivity kappa requires taking the q to 0 limit. The method is tested for liquid and crystalline argon. One advantage is reduced computational noise at a given total MD run time. Another advantage is that kappa(q) has a more physical meaning than kappa(L). It can be easily studied using Peierls-Boltzmann transport theory. New formulas for kappa(q) in simplified Debye-type models give new insight about extrapolation to q to 0 or 1/L to 0. In particular, it is shown that kappa(L$ is unlikely to behave as kappa -C/L, and much more likely to behave as kappa-C'/sqrt(L). Convergence problems encountered in computational cells with very large aspect ratios L(parallel)/L(perp) are also analyzed. Some details are contained in the "Supplemental Material" file.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anharmonic Phonon Quasiparticle Theory of Zero-point and Thermal Shifts in Insulators: Heat Capacity, Bulk Modulus, and Thermal Expansion

The Quasi-harmonic (QH) approximation uses harmonic vibrational frequencies omega(H,Q,V), computed at volumes V near the volume where the Born-Oppenheimer (BO) energy is minimum. When this is used in the harmonic free energy, QH approximation gives a good zeroth order theory of thermal expansion, and first order theory of bulk modulus. Here, n-th order means smaller than the leading term by n powers of epsilon, where epsilon is the ratio hbar omega(Q)/E(el) or kT/E(el), and E(el) is an electronic energy scale, typically 2 to 10 eV. Experiment often shows evidence for next order corrections. When such corrections are needed, anharmonic interactions must be included. The most accessible measure of anhamonicity is the quasiparticle (QP) energy, omega(Q,V,T), seen experimentally by vibrational spectroscopy. However, this cannot just be inserted into the harmonic free energy F(H). In this paper, a free energy formula is found which corrects the double-counting of anharmonic interactions that is made when F is approximated by F(H,omega(Q,V,T)). The term "QP thermodynamics" is used for this way of treating anharmonicity. It enables (n+1)-order corrections, if QH theory is accurate to order n. This procedure is used to give corrections to specific heat and volume thermal expansion. The QH formulas for isothermal and adiabatic bulk moduli are clarified, and the route to higher order corrections is indicated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electron Self-Energy and Generalized Drude Formula for Infrared Conductivity of Metals

Goetze and Woelfle (GW) wrote the conductivity in terms of a memory function M as (ine2/m)/(omega+M(omega)), where M=i/tau in the Drude limit. The analytic properties of -M are the same as those of the self-energy of a retarded Green's function. In the approximate treatment of GW, -M closely resembles a self-energy, with differences, e.g., the imaginary part is twice too large. The correct relation between -M and the self-energy is known for the electron-phonon case and is conjectured to be similar for other perturbations. When vertex corrections are ignored there is a known relation. A derivation using Matsubara temperature Green's functions is given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electronic and Nuclear Quantum Effects on the Ice XI/Ice Ih Phase Transition

We study the isotope effect on the temperature of the proton order/disorder phase transition between ice XI and ice Ih, using the quasiharmonic approximation combined with \textit{ab initio} density functional theory calculations. We show that this method is accurate enough to obtain a phase transition temperature difference between light ice (H$_2$O) and heavy ice (D$_2$O) of 6 K as compared to the experimental value of 4 K. More importantly, we are able to explain the origin of the isotope effect on the much debated large temperature difference observed in the phase transition. The source of the difference is directly linked to the physics behind the anomalous isotope effect on the volume of hexagonal ice that was recently explained in [Phys. Rev. Lett. 108, 193003 (2012)]. These results indicate that the same physics might be behind the isotope effects in transition temperatures between other ice phases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Size Effects in Thermal Conduction by Phonons

Heat transport in nanoscale systems is both hard to measure microscopically, and hard to interpret. Ballistic and diffusive heat flow coexist, adding confusion. This paper looks at a very simple case: a nanoscale crystal repeated periodically. This is a popular model for simulation of bulk heat transport using classical molecular dynamics (MD), and is related to transient thermal grating experiments. Nanoscale effects are seen in perhaps their simplest form. The model is solved by an extension of standard quasiparticle gas theory of bulk solids. Both structure and heat flow are constrained by periodic boundary conditions. Diffusive transport is fully included, while ballistic transport by phonons of long mean free path is diminished in a specific way. Heat current J(x) and temperature gradient dT(x')/dx' have a non-local relationship, via k(x-x'), over a distance |x-x'| determined by phonon mean free paths. In MD modeling of bulk conductivity, finite computer resources limit system size. Long mean free paths, comparable to the scale of heating and cooling, cause undesired finite-size effects that have to be removed by extrapolation. The present model allows this extrapolation to be quantified. Calculations based on the Peierls-Boltzmann equation, using a generalized Debye model, show that extrapolation involves fractional powers of 1/L. It is also argued that heating and cooling should be distributed sinusoidally (de/dt~cos(2pi x/L)) to improve convergence of numerics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Special Quasi-ordered Structures: role of short-range order in the semiconductor alloy (GaN)$_{1-x}$(ZnO)$_x$

This paper studies short-range order (SRO) in the semiconductor alloy (GaN)$_{1-x}$(ZnO)$_x$. Monte Carlo simulations performed on a density functional theory (DFT)-based cluster expansion model show that the heterovalent alloys exhibit strong SRO because of the energetic preference for the valence-matched nearest-neighbor Ga-N and Zn-O pairs. To represent the SRO-related structural correlations, we introduce the concept of Special Quasi-ordered Structure (SQoS). Subsequent DFT calculations reveal dramatic influence of SRO on the atomic, electronic and vibrational properties of the (GaN)$_{1-x}$(ZnO)$_x$ alloy. Due to the enhanced statistical presence of the energetically unfavored Zn-N bonds with the strong Zn3$d$-N2$p$ repulsion, the disordered alloys exhibit much larger lattice bowing and band-gap reduction than those of the short-range ordered alloys. Inclusion of lattice vibrations stabilizes the disordered alloy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolutionary Method for Predicting Surface Reconstructions with Variable Stoichiometry

We present a specially designed evolutionary algorithm for the prediction of surface reconstructions. This new technique allows one to automatically explore all the low-energy configurations with variable surface atoms and variable surface unit cells through the whole chemical potential range. The power of evolutionary search is demonstrated by the efficient identification of diamond 2*1 (100) and 2*1 (111) surfaces with a fixed number of surface atom and a fixed cell size. With further variation of surface unit cells, we study the reconstructions of the polar surface MgO (111). Experiment has detected an oxygen trimer (ozone) motif (Plass et al, 1998). We predict a new version of this motif which can be thermodynamically stable at extreme oxygen rich condition. Finally, we perform a variable stoichiometry search for a complex ternary system: semi-polar GaN (101bar1) with and without adsorbed oxygen. The search yields a non-intuitive reconstruction based on N3-trimers. These examples demonstrate that an automated scheme to explore the energy landscape of surfaces will improve our understanding of surface reconstructions. The method presented in this report can be generally applied to binary and multi-component systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Heat Conduction by Phonons across a Film

Quasiparticle theory gives a local relation between heat current and temperature gradient, provided the quasiparticle mean free path is smaller than the scale of variation of temperature. When mean free paths are comparable to sample size, the relation becomes non-local. This non-local relation is formulated for phonon carriers; an explicit form is found in the approximation where current relaxation rates are replaced by quasiparticle relaxation rates. A quasiparticle definition of local temperature is offered. To extract the spatial variation of the temperature gradient requires inverting the non-local relation. A variational principle is constructed. The heat current is bounded above when evaluated for a trial temperature gradient. The true temperature gradient is the one which minimizes heat current. The simplest variational approximation (a constant temperature gradient) gives an approximate formula for the thickness dependence of the heat transport across a film whose thickness is comparable to mean free paths of phonons. The formula interpolates between ballistic and diffusive limits.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improved Callaway model for Lattice Thermal Conductivity

In developing the phonon quasiparticle picture, Peierls discovered that, in a perfect crystal, without Umklapp (U) events, a current-carrying distribution can never relax to a zero-current distribution. Callaway introduced a simplified approximate model version of the Peierls-Boltzmann equation, retaining its the ability to deal separately with Normal (N) and U events. This paper clarifies and improves the Callaway model, and shows that Callaway underestimated the suppression of N-processes in relaxing thermal current. The new result should improve computations of thermal conductivity from relaxation-time studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Temperature and composition dependence of short-range order and entropy, and statistics of bond length: the semiconductor alloy (GaN)$_{1-x}$(ZnO)$_x$

We present total energy and force calculations on the (GaN)$_{1-x}$(ZnO)$_{x}$ alloy. Site-occupancy configurations are generated by Monte Carlo (MC) simulations, based on a cluster expansion (CE) model proposed in a previous study. Surprisingly large local atomic coordinate relaxations are found by density-functional calculations using a 432-atom periodic supercell, for three representative configurations at $x=0.5$. These are used to generate bond length distributions. The configurationally averaged composition- and temperature-dependent short-range order (SRO) parameters of the alloys are discussed. Entropy is approximated in terms of pair distribution statistics and thus related to SRO parameters. This approximate entropy is compared with accurate numerical values from MC. An empirical model for the dependence of bond length on local chemical environments is proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Octahedral Tilt Instability of ReO_3-type Crystals

The octahedron tilt transitions of ABX_3 perovskite-structure materials lead to an anti-polar (or antiferroelectric) arrangement of dipoles, with the low temperature structure having six sublattices polarized along various crystallographic directions. It is shown that an important mechanism driving the transition is long range dipole-dipole forces acting on both displacive and induced parts of the anion dipole. This acts in concert with short range repulsion, allowing a gain of electrostatic (Madelung) energy, both dipole-dipole and charge-charge, because the unit cell shrinks when the hard ionic spheres of the rigid octahedron tilt out of linear alignment.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Numerical study of resistivity of model disordered three-dimensional metals

We calculate the zero-temperature resistivity of model 3-dimensional disordered metals described by tight-binding Hamiltonians. Two different mechanisms of disorder are considered: diagonal and off-diagonal. The non-equilibrium Green function formalism provides a Landauer-type formula for the conductance of arbitrary mesoscopic systems. We use this formula to calculate the resistance of finite-size disordered samples of different lengths. The resistance averaged over disorder configurations is linear in sample length and resistivity is found from the coefficient of proportionality. Two structures are considered: (1) a simple cubic lattice with one s-orbital per site, (2) a simple cubic lattice with two d-orbitals. For small values of the disorder strength, our results agree with those obtained from the Boltzmann equation. Large off-diagonal disorder causes the resistivity to saturate, whereas increasing diagonal disorder causes the resistivity to increase faster than the Boltzmann result. The crossover toward localization starts when the Boltzmann mean free path relative to the lattice constant has a value between 0.5 and 2.0 and is strongly model dependent.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Toward a Theory of Orbiton Dispersion in LaMnO_3

At 750K, LaMnO_3 has a cooperative Jahn-Teller (JT) distortion, with Mn atoms in distorted oxygen octahedra. This lifts the degeneracy of the singly-occupied $e_g$ orbitals of the Mn$^{3+}$ ions, which then become orbitally ordered. We use a minimal model to describe the ordered phase at T=0. The on-site Coulomb repulsion $U$ is set to infinity. There are two electronic orbitals and three oxygen vibrational coordinates per unit cell. In addition to spin excitations and phonons, the model has electronic excitations consisting of mis-orienting orbitals on Mn ions. Neglecting coupling to the oxygen displacements, the gap to such excitations is $2Δ=16g^2/K$ where $g$ is the electron-phonon coupling and $K$ is the oxygen spring constant. When static oxygen displacements are coupled, this excitation becomes a self-trapped exciton with energy $Δ$, half the JT gap. Adding dynamic oxygen displacements in one-phonon approximation introduces dispersion to both the (previously Einstein-like) phonons and the orbital excitons (``orbitons''). One of the phonon branches has zero frequency at $\vec{k}=(0,0,π)$. This is the Goldstone mode of the JT broken symmetry.

0 citations0 reviewsNo code linkedOpen access