Researcher profile

Zheyong Fan

Zheyong Fan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 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.

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.

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.

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.

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.

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.

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.

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.

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.