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

Zheyong Fan

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

ResearcherUnclaimed source record
cond-mat.mes-hallphysics.comp-phcond-mat.mtrl-scicond-mat.stat-mechcond-mat.dis-nncond-mat.str-elquant-ph

Catalog footprint

What is connected

21works
7topics
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

21 published item(s)

preprint2022arXiv

Exactly equivalent thermal conductivity in finite systems from equilibrium and nonequilibrium molecular dynamics simulations

In a previous paper [Physical Review B \textbf{103}, 035417 (2021)], we showed that the equilibrium molecular dynamics (EMD) method can be used to compute the apparent thermal conductivity of finite systems. It has been shown that the apparent thermal conductivity from EMD for a system with domain length $2L$ is equal to that from nonequilibrium molecular dynamics (NEMD) for a system with domain length $L$. Taking monolayer silicence with an accurate machine learning potential as an example, here we show that the thermal conductivity values from EMD and NEMD agree for the same domain length if the NEMD is applied with periodic boundary conditions in the transport direction. Our results thus establish an exact equivalence between EMD and NEMD for thermal conductivity calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improving the accuracy of the neuroevolution machine learning potential for multi-component systems

In a previous paper [Fan Z \textit{et al}. 2021 Phys. Rev. B, \textbf{104}, 104309], we developed the neuroevolution potential (NEP), a framework of training neural network based machine-learning potentials using a natural evolution strategy and performing molecular dynamics (MD) simulations using the trained potentials. The atom-environment descriptor in NEP was constructed based on a set of radial and angular functions. For multi-component systems, all the radial functions between two atoms are multiplied by some fixed factors that depend on the types of the two atoms only. In this paper, we introduce an improved descriptor for multi-component systems, in which different radial functions are multiplied by different factors that are also optimized during the training process, and show that it can significantly improve the regression accuracy without increasing the computational cost in MD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport

We develop a neuroevolution-potential (NEP) framework for generating neural network based machine-learning potentials. They are trained using an evolutionary strategy for performing large-scale molecular dynamics (MD) simulations. A descriptor of the atomic environment is constructed based on Chebyshev and Legendre polynomials. The method is implemented in graphic processing units within the open-source GPUMD package, which can attain a computational speed over $10^7$ atom-step per second using one Nvidia Tesla V100. Furthermore, per-atom heat current is available in NEP, which paves the way for efficient and accurate MD simulations of heat transport in materials with strong phonon anharmonicity or spatial disorder, which usually cannot be accurately treated either with traditional empirical potentials or with perturbative methods.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations

We study the interfacial thermal conductance of grain boundaries (GBs) between monolayer graphene and hexagonal boron nitride (h-BN) sheets using a combined atomistic approach. First, realistic samples containing graphene/h-BN GBs with different tilt angles are generated using the phase-field crystal (PFC) model developed recently [P. Hirvonen \textit{et al.}, Phys. Rev. B \textbf{100}, 165412 (2019)] that captures slow diffusive relaxation inaccessible to molecular dynamics (MD) simulations. Then, large-scale MD simulations using the efficient GPUMD package are performed to assess heat transport and rectification properties across the GBs. We find that lattice mismatch between the graphene and h-BN sheets plays a less important role in determining the interfacial thermal conductance as compared to the tilt angle. In addition, we find no significant thermal rectification effects for these GBs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations

We propose a way to properly interpret the apparent thermal conductivity obtained for finite systems using equilibrium molecular dynamics simulations (EMD) with fixed or open boundary conditions in the transport direction. In such systems the heat current autocorrelation function develops negative values after a correlation time which is proportional to the length of the simulation cell in the transport direction. Accordingly, the running thermal conductivity develops a maximum value at the same correlation time and eventually decays to zero. By comparing EMD with nonequilibrium molecular dynamics (NEMD) simulations, we conclude that the maximum thermal conductivity from EMD in a system with domain length 2L is equal to the thermal conductivity from NEMD in a system with domain length L. This facilitates the use of nonperiodic-boundary EMD for thermal transport in finite samples in close correspondence to NEMD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study

Single-walled carbon nanotubes (SWCNTs) in their pristine form have high thermal conductivity whose further improvement has attracted a lot of interest. Some theoretical studies have suggested that the thermal conductivity of a $(10,10)$ SWCNT is dramatically enhanced by C$_{60}$ fullerene encapsulation. However, recent experiments on SWCNT bundles show that fullerene encapsulation leads to a reduction rather than an increase in thermal conductivity. Here, we employ three different molecular dynamics methods to study the influence of C$_{60}$ encapsulation on heat transport in a $(10,10)$ SWCNT. All the three methods consistently predict a reduction of the thermal conductivity of $(10,10)$ SWCNT upon C$_{60}$ encapsulation by $20\%-30\%$, in agreement with experimental results on bundles of SWCNTs. We demonstrate that there is a simulation artifact in the Green-Kubo method which gives anomalously large thermal conductivity from artificial convection. Our results show that the C$_{60}$ molecules conduct little heat compared to the outer SWCNT and reduce the phonon mean free paths of the SWCNT by inducing extra phonon scattering. We also find that the thermal conductivity of a $(10,10)$ SWCNT monotonically decreases with increasing filling ratio of C$_{60}$ molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties

Silicon is an important material and many empirical interatomic potentials have been developed for atomistic simulations of it. Among them, the Tersoff potential and its variants are the most popular ones. However, all the existing Tersoff-like potentials fail to reproduce the experimentally measured thermal conductivity of diamond silicon. Here we propose a modified Tersoff potential and develop an efficient open source code called GPUGA (graphics processing units genetic algorithm) based on the genetic algorithm and use it to fit the potential parameters against energy, virial and force data from quantum density functional theory calculations. This potential, which is implemented in the efficient open source GPUMD (graphics processing units molecular dynamics) code, gives significantly improved descriptions of the thermal conductivity and phonon dispersion of diamond silicon as compared to previous Tersoff potentials and at the same time well reproduces the elastic constants. Furthermore, we find that quantum effects on the thermal conductivity of diamond silicon at room temperature are non-negligible but small: using classical statistics underestimates the thermal conductivity by about 10\% as compared to using quantum statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime, and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it to the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification. Our findings are important to obtain correct results from molecular dynamics simulations of nanoscale heat transport as the accuracy of the interatomic potentials is rapidly improving.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dominant source of disorder in graphene: Charged impurities or ripples?

Experimentally produced graphene sheets exhibit a wide range of mobility values. Both extrinsic charged impurities and intrinsic ripples (corrugations) have been suggested to induce long-range disorder in graphene and could be a candidate for the dominant source of disorder. Here, using large-scale molecular dynamics and quantum transport simulations, we find that the hopping disorder and the gauge and scalar potentials induced by the ripples are short-ranged, in strong contrast with predictions by continuous models, and the transport fingerprints of the ripple disorder are very different from those of charged impurities. We conclude that charged impurities are the dominant source of disorder in most graphene samples, whereas scattering by ripples is mainly relevant in the high carrier density limit of ultraclean graphene samples (with a charged impurity concentration < 10 ppm) at room and higher temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale modeling of polycrystalline graphene: A comparison of structure and defect energies of realistic samples from phase field crystal models

We extend the phase field crystal (PFC) framework to quantitative modeling of polycrystalline graphene. PFC modeling is a powerful multiscale method for finding the ground state configurations of large realistic samples that can be further used to study their mechanical, thermal or electronic properties. By fitting to quantum-mechanical density functional theory (DFT) calculations, we show that the PFC approach is able to predict realistic formation energies and defect structures of grain boundaries. We provide an in-depth comparison of the formation energies between PFC, DFT and molecular dynamics (MD) calculations. The DFT and MD calculations are initialized using atomic configurations extracted from PFC ground states. Finally, we use the PFC approach to explicitly construct large realistic polycrystalline samples and characterize their properties using MD relaxation to demonstrate their quality.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations

We derive expressions of interatomic force and heat current for many-body potentials such as the Tersoff, the Brenner, and the Stillinger-Weber potential used extensively in molecular dynamics simulations of covalently bonded materials. Although these potentials have a many-body nature, a pairwise force expression that follows Newton's third law can be found without referring to any partition of the potential. Based on this force formula, a stress applicable for periodic systems can be unambiguously defined. The force formula can then be used to derive the heat current formulas using a natural potential partitioning. Our heat current formulation is found to be equivalent to most of the seemingly different heat current formulas used in the literature, but to deviate from the stress-based formula derived from two-body potential. We validate our formulation numerically on various systems descried by the Tersoff potential, namely three-dimensional silicon and diamond, two-dimensional graphene, and quasi-one-dimensional carbon nanotube. The effects of cell size and time used in the simulation are examined.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impurity-induced quantum phase transition in finite Heisenberg spin chains: Criteria for existence and stability

A quantum phase transition may occur in a system at zero temperature when a controlling parameter is tuned towards a critical point. An important question is whether such a critical point exists in a particular system and how stable it is. Here, we identify the critical point of a quantum phase transition as a singular point in the affine algebraic variety of the characteristic equation for the Hamiltonian describing the system, with an unstable critical point being associated with an isolated singular point which has a finite Tjurina number. The theory is illustrated by studying a model system of zero-dimensional (finite) Heisenberg spin chain with an impurity, which exhibits a nontrivial first-order quantum phase transition. Both analytical and numerical calculations show that the quantum phase transition always exists when the impurity has a $Z_2$ symmetry but only remains in systems with an even number of spin sites when the $Z_2$ symmetry is broken.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Silicon and silicon-nitrogen impurities in graphene: structure, energetics and effects on electronic transport

We theoretically study the atomic structure and energetics of silicon and silicon-nitrogen impurities in graphene. Using density-functional theory, we get insight into the atomic structures of the impurities, evaluate their formation energies and assess their abundance in realistic samples. We find that nitrogen, as well as oxygen and hydrogen, are trapped at silicon impurities, considerably altering the electronic properties of the system. Furthermore, we show that nitrogen doping can induce local magnetic moments resulting in spin-dependent transport properties, even though neither silicon nor nitrogen impurities are magnetic by themselves. To simulate large systems with many randomly distributed impurities, we derive tight-binding models that describe the effects of the impurities on graphene π electron structure. Then by using the linear-scaling real-space Kubo-Greenwood method, we evaluate the transport properties of large-scale systems with random distribution of impurities, and find the fingerprint-like scattering cross sections for each impurity type. The transport properties vary widely, and our results indicate that some of the impurities can even induce strong localization in realistic graphene samples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anderson localization in two-dimensional graphene with short-range disorder: One-parameter scaling and finite-size effects

We study Anderson localization in graphene with short-range disorder using the real-space Kubo-Greenwood method implemented on graphics processing units. Two models of short-range disorder, namely, the Anderson on-site disorder model and the vacancy defect model, are considered. For graphene with Anderson disorder, localization lengths of quasi-one-dimensional systems with various disorder strengths, edge symmetries, and boundary conditions are calculated using the real-space Kubo-Greenwood formalism, showing excellent agreement with independent transfer matrix calculations and superior computational efficiency. Using these data, we demonstrate the applicability of the one-parameter scaling theory of localization length and propose an analytical expression for the scaling function, which provides a reliable method of computing the two-dimensional localization length. This method is found to be consistent with another widely used method which relates the two-dimensional localization length to the elastic mean free path and the semiclassical conductivity. Abnormal behavior at the charge neutrality point is identified and interpreted to be caused by finite-size effects when the system width is comparable to or smaller than the elastic mean free path. We also demonstrate the finite-size effect when calculating the two-dimensional conductivity in the localized regime and show that a renormalization group beta function consistent with the one-parameter scaling theory can be extracted numerically. For graphene with vacancy disorder, we show that the proposed scaling function of localization length also applies. Lastly, we discuss some ambiguities in calculating the semiclassical conductivity around the charge neutrality point due to the presence of resonant states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electronic and transport properties in geometrically disordered graphene antidot lattices

A graphene antidot lattice, created by a regular perforation of a graphene sheet, can exhibit a considerable band gap required by many electronics devices. However, deviations from perfect periodicity are always present in real experimental setups and can destroy the band gap. Our numerical simulations, using an efficient linear-scaling quantum transport simulation method implemented on graphics processing units, show that disorder that destroys the band gap can give rise to a transport gap caused by Anderson localization. The size of the defect induced transport gap is found to be proportional to the radius of the antidots and inversely proportional to the square of the lattice periodicity. Furthermore, randomness in the positions of the antidots is found to be more detrimental than randomness in the antidot radius. The charge carrier mobilities are found to be very small compared to values found in pristine graphene, in accordance with recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Obtaining localization properties efficiently using the Kubo-Greenwood formalism

We establish, through numerical calculations and comparisons with a recursive Green's function based implementation of the Landauer-Büttiker formalism, an efficient method for studying Anderson localization in quasi-one-dimensional and two-dimensional systems using the Kubo-Greenwood formalism. Although the recursive Green's function method can be used to obtain the localization length of a mesoscopic conductor, it is numerically very expensive for systems that contain a large number of atoms transverse to the transport direction. On the other hand, linear-scaling has been achieved with the Kubo-Greenwood method, enabling the study of effectively two-dimensional systems. While the propagating length of the charge carriers will eventually saturate to a finite value in the localized regime, the conductances given by the Kubo-Greenwood method and the recursive Green's function method agree before the saturation. The converged value of the propagating length is found to be directly proportional to the localization length obtained from the exponential decay of the conductance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations

In this paper, we develop a highly efficient molecular dynamics code fully implemented on graphics processing units for thermal conductivity calculations using the Green-Kubo formula. We compare two different schemes for force evaluation, a previously used thread-scheme where a single thread is used for one particle and each thread calculates the total force for the corresponding particle, and a new block-scheme where a whole block is used for one particle and each thread in the block calculates one or several pair forces between the particle associated with the given block and its neighbor particle(s) associated with the given thread. For both schemes, two different classical potentials, namely, the Lennard-Jones potential and the rigid-ion potential are implemented. While the thread-scheme performs a little better for relatively large systems, the block-scheme performs much better for relatively small systems. The relative performance of the block-scheme over the thread-scheme also increases with the increasing of the cutoff radius. We validate the implementation by calculating lattice thermal conductivities of solid argon and lead telluride.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient linear-scaling quantum transport calculations on graphics processing units and applications on electron transport in graphene

We implement, optimize, and validate the linear-scaling Kubo-Greenwood quantum transport simulation on graphics processing units by examining resonant scattering in graphene. We consider two practical representations of the Kubo-Greenwood formula: a Green-Kubo formula based on the velocity auto-correlation and an Einstein formula based on the mean square displacement. The code is fully implemented on graphics processing units with a speedup factor of up to 16 (using double-precision) relative to our CPU implementation. We compare the kernel polynomial method and the Fourier transform method for the approximation of the Dirac delta function and conclude that the former is more efficient. In the ballistic regime, the Einstein formula can produce the correct quantized conductance of one-dimensional graphene nanoribbons except for an overshoot near the band edges. In the diffusive regime, the Green-Kubo and the Einstein formalisms are demonstrated to be equivalent. A comparison of the length-dependence of the conductance in the localization regime obtained by the Einstein formula with that obtained by the non-equilibrium Green's function method reveals the challenges in defining the length in the Kubo-Greenwood formalism at the strongly localized regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A molecular dynamics investigation of the mechanical properties of graphene nanochain

In this paper, we investigate, by molecular dynamics simulations, the mechanical properties of a new carbon nanostructure, termed graphene nanochain, constructed by sewing up pristine or twisted graphene nanoribbons (GNRs) and interlocking the obtained nanorings. The obtained tensile strength of defect-free nanochain is a little lower than that of pristine GNRs and the fracture point is earlier than that of the GNRs. The effects of length, width and twist angle of the constituent GNRs on the mechanical performance are analyzed. Furthermore, defect effect is investigated and in some high defect coverage cases, an interesting mechanical strengthening-like behavior is observed. This structure supports the concept of long-cable manufacturing and advanced material design can be achieved by integration of nanochain with other nanocomposites. The technology used to construct the nanochain is experimentally feasible, inspired by the recent demonstrations of atomically precise fabrications of GNRs with complex structures [Phys. Rev. Lett,2009,\textbf{102},205501; Nano Lett., 2010, \textbf{10},4328; Nature,2010,\textbf{466},470]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Braid Matrices and Quantum Gates for Ising Anyons Topological Quantum Computation

We study various aspects of the topological quantum computation scheme based on the non-Abelian anyons corresponding to fractional quantum hall effect states at filling fraction 5/2 using the Temperley-Lieb recoupling theory. Unitary braiding matrices are obtained by a normalization of the degenerate ground states of a system of anyons, which is equivalent to a modification of the definition of the 3-vertices in the Temperley-Lieb recoupling theory as proposed by Kauffman and Lomonaco. With the braid matrices available, we discuss the problems of encoding of qubit states and construction of quantum gates from the elementary braiding operation matrices for the Ising anyons model. In the encoding scheme where 2 qubits are represented by 8 Ising anyons, we give an alternative proof of the no-entanglement theorem given by Bravyi and compare it to the case of Fibonacci anyons model. In the encoding scheme where 2 qubits are represented by 6 Ising anyons, we construct a set of quantum gates which is equivalent to the construction of Georgiev.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constructing quantum circuits for maximally entangled multi-qubit states using the genetic algorithm

Numerical optimization methods such as hillclimbing and simulated annealing have been applied to search for highly entangled multi-qubit states. Here the genetic algorithm is applied to this optimization problem -- to search not only for highly entangled states, but also for the corresponding quantum circuits creating these states. Simple quantum circuits for maximally (highly) entangled states are discovered for 3, 4, 5, and 6-qubit systems; and extension of the method to systems with more qubits is discussed. Among other results we have found explicit quantum circuits for maximally entangled 5 and 6-qubit circuits, with only 8 and 13 quantum gates respectively. One significant advantage of our method over previous ones is that it allows very simple construction of quantum circuits based on the quantum states found.

0 citations0 reviewsNo code linkedOpen access