Source author record

Peng Yu

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2026arXiv

Quadrature-TreeSHAP: Depth-Independent TreeSHAP and Shapley Interactions

Shapley values are a standard tool for explaining predictions of tree ensembles, with Path-Dependent SHAP being the most widely used variant. Despite substantial progress, existing methods still exhibit trade-offs between depth-dependent runtime, numerical stability, and support for higher-order interactions. To address these challenges, we introduce Quadrature-TreeSHAP, a quadrature-based reformulation of Path-Dependent TreeSHAP that is numerically stable, naturally extends to any-order Shapley interaction values and is practically insensitive to tree depth. Our implementation supports both CPU and GPU and is integrated into XGBoost. Our method is based on a weighted-Banzhaf interaction polynomial, which expresses Banzhaf interaction values as expectations under a feature participation probability $p$. Shapley values and any-order interaction values are then recovered by integrating these polynomials over $p$ from 0 to 1. We evaluate these integrals using Gauss-Legendre quadrature, and show that, in practice, only 8 fixed quadrature points are sufficient to reach machine precision. In fact, Quadrature-TreeSHAP with 8 fixed points achieves greater numerical stability than TreeSHAP. This fixed-point formulation removes depth dependence from the inner computation and enables efficient SIMD execution. We confirm these advantages empirically. On 12 XGBoost benchmarks, Quadrature-TreeSHAP computes Shapley values 1.06x-10.59x faster than TreeSHAP on CPU and 1.84x-6.95x faster than GPUTreeSHAP on GPU. Shapley pairwise interactions are 3.80x-58.11x faster on CPU, with higher-order interactions achieving speedups of up to 1200x compared to TreeSHAP-IQ.

preprint2024arXiv

Hot electrons and electromagnetic effects in the broadband Au, Ag, and Ag-Au nanocrystals: The UV, visible, and NIR plasmons

Energetic and optical properties of plasmonic nanocrystals strongly depend on their sizes, shapes, and composition. Whereas using plasmonic nanoparticles in biotesting has become routine, applications of plasmonics in energy are still early in development. Here, we investigate hot electron (HE) generation and related electromagnetic effects in both mono- and bi-metallic nanorods (NRs) and focus on one promising type of bi-metallic nanocrystals - core-shell Au-Ag nanorods. The spectra of the NRs are broadband, highly tunable with their geometry, and have few plasmon resonances. In this work, we provide a new quantum formalism describing the HE generation in bi-metallic nanostructures. Interestingly, we observe that the HE generation rate at the UV plasmon resonance of Au-Ag NRs appears to be very high. These HEs are highly energetic and suitable for carbon-fuel reactions. Simultaneously, the HE generation at the longitudinal plasmon (L-plasmon) peaks, which can be tuned from the yellow to near-IR, depends on the near-field and electromagnetic Mie effects, limiting the HE efficiencies for the long and large NRs. These properties of the L-plasmon relate to all kinds of NRs (Au, Ag, and Au-Ag). We also consider the generation of the interband d-holes in Au and Ag, since the involvement of the d-band is crucial for the energetic properties of UV plasmons. The proposed formalism is an important development for the description of bi-metallic (or tri-metallic, or more complex) nanostructures, and it paves the way to the efficient application of the plasmonic HEs and hot holes in sensing, nanotechnology, photocatalysis, and electrophotochemistry.

preprint2022arXiv

A machine learning based method to generate random packed isotropic porous media with desired porosity and permeability

Porous materials are used in many fields, including energy industry, agriculture, medical industry, etc. The generation of digital porous media facilitates the fabrication of real porous media and the analysis of their properties. The past random digital porous media generation methods are unable to generate a porous medium with a specific permeability. A new method is proposed in the present study, which can generate the random packed isotropic porous media with specific porosity and permeability. Firstly, the process of generating the random packed isotropic porous media is detailed. Secondly, the permeability of the generated porous media is calculated with the multi-relaxation time (MRT) lattice Boltzmann method (LBM), which is prepared for the training of convolutional neural network (CNN). Thirdly, 3000 samples on the microstructure of porous media and their permeabilities are used to train the CNN model. The trained model is very effective in predicting the permeability of a porous medium. Finally, our method is elaborated and the choice of target permeability in this method is discussed. With the support of a powerful computer, a porous medium that satisfies the error condition of porosity and permeability can be generated in a short time.

preprint2022arXiv

Comparison of three numerical stabilization techniques of viscoelastic flows: vortex shedding behind a confined cylinder

In this study, the OpenFOAM platform, based on the finite volume method, is applied to investigate the two-dimensional viscoelastic flow past a circular cylinder. The FENE-P model, which considers the bounded elongation of polymer molecules, is chosen to describe the elastic constitutive relationship of the polymer solution. The maximum molecular chain lengths of L = 10, 50, 100, and 200 are considered, which describe the molecular conformation characteristics of the polymer solution. To improve the numerical instability of the viscoelastic flow simulation, three different methods, i.e., the traditional method (Td) with the addition of artificial viscosity, the logarithmic reconstruction method (Log), and the square root tensor method (Sqrt), are evaluated. The results show that the artificial viscosity has a little effect on the accuracy for the simulation with a small molecular chain length (L = 10). However, for long molecular chain lengths such as L = 100 and L = 200, the addition of artificial dissipation tends to overestimate the drag, which indicates that special caution is needed to incorporate the artificial dissipation in the simulation. Moreover, the logarithmic reconstruction method shows a strong grid-dependent characteristics, which may produce unphysical results.

preprint2022arXiv

Direct Light Orbital Angular Momentum Detection in Mid-Infrared based on Type-II Weyl Semimetal TaIrTe4

The capability of direct photocurrent detection of orbital angular momentum (OAM) of light has recently been realized with topological Weyl semimetal, but limited to near infrared wavelength range. The extension of direct OAM detection to midinfrared, a wavelength range that plays important role in a vast range of applications, such as autonomous driving, night vision and motion detection, is challenging and has not yet been realized. This is because most studies of photocurrent responses are not sensitive to the phase information and the photo response is usually very poor in the mid-infrared. In this study, we designed a photodetector based on Type-II Weyl semimetal tantalum iridium tellurides with designed electrode geometries for direct detection of the topological charge of OAM through orbital photogalvanic effect. Our results indicate helical phase gradient of light can be distinguished by a current winding around the optical beam axis with a magnitude proportional to its quantized OAM mode number. The topological enhanced response at mid-infrared of TaIrTe4 further help overcome the low responsivity issues and finally render the direct orbital angular momentum detection capability in mid-infrared. Our study enables on-chip integrated OAM detection, and thus OAM sensitive focal plane arrays in mid-infrared. Such capability triggers new route to explore applications of light carrying OAM, especially that it can crucially promote the performance of many mid-infrared imaging related applications, such as intricate target recognition and night vision.

preprint2022arXiv

Shear-induced droplet mobility within porous surfaces

Droplet mobility under shear flows is important in a wide range of engineering applications, e.g., fog collection, and self-cleaning surfaces. For structured surfaces to achieve superhydrophobicity, the removal of stains adhered within the microscale surface features strongly determines the functional performance and durability. In this study, we numerically investigate the shear-induced mobility of the droplet trapped within porous surfaces. Through simulations covering a wide range of flow conditions and porous geometries, three droplet mobility modes are identified, i.e., the stick-slip, crossover, and slugging modes. To quantitatively characterise the droplet dynamics, we propose a droplet-scale capillary number that considers the driving force and capillary resistance. By comparing against the simulation results, the proposed dimensionless number presents a strong correlation with the leftover volume. The dominating mechanisms revealed in this study provide a basis for further research on enhancing surface cleaning and optimising design of anti-fouling surfaces.

preprint2022arXiv

Subcritical insability of viscoelastic flow over a circular cylinder: A numerical study

In this paper, we discuss whether the instability of viscoelastic flow around a circular cylinder is subcritical or supercritical by numerical simulation. The Oldroyd-B model is selected to describe the viscoelastic constitutive relationship. The Log-conformation reformulation is employed to stabilize numerical simulation. The parameter ranges investigated are the Reynolds numbers ($Re$) spanning from 5 to 100 and the Weissenberg number ($Wi$) spanning from 0 to 10, with a fixed viscosity ratio of $β= 0.9$. Simulations are performed under two paths, i.e., 1) increasing $Re$ (or $Wi$) and 2) decreasing $Re$ (or $Wi$) slowly and gradually from one state to the next. The results show that the statistical solutions such as the time-averaged velocity obtained along the two paths are not identical over certain parameter range, which is around the transition point from the steady to unsteady flow. This loading path dependence behaviour indicates that the flow instability is subcritical.

preprint2022arXiv

The reconstructed thermal lattice Boltzmann flux solver and its applications for simulations of thermal flows

The thermal lattice Boltzmann flux solver (TLBFS) has been proposed to overcome the drawbacks of the thermal lattice Boltzmann models. However, as a weakly compressible model, its mechanism of good numerical stability for high Rayleigh number thermal flows is still unclear. To reveal the mechanism, the present paper firstly derives the macroscopic equations of TLBFS (MEs-TLBFS) with actual numerical dissipation terms by approximating its computational process. By solving MEs-TLBFS with the finite volume method, the reconstructed TLBFS (RTLBFS) is proposed. Detailed analyses prove that these actual numerical dissipation terms are the mechanism of the good numerical stability of TLBFS for high Rayleigh number thermal flows. More detailed numerical tests indicate RTLBFS has similar performances as TLBFS for stability, accuracy, and efficiency. Moreover, the present RTLBFS shows a clear mechanism to achieve good numerical stability for high Rayleigh number flows.

preprint2022arXiv

Wake asymmetry weakening in viscoelastic fluids: Numerical discovery and mechanism exploration

Viscoelasticity weakens the asymmetry of laminar shedding flow behind a blunt body in a free domain. In the present study, this finding is confirmed by four unsteady viscoelastic flows with asymmetric flow configuration, i.e., flow over an inclined flat plate with various angles of incidence, flow over a rotating circular cylinder, flow over a circular cylinder with asymmetric slip boundary distribution, and flow over an inclined row of eight equally closely spaced circular cylinders (which can be considered as a single large blunt body) through direct numerical simulation combined with the Peterlin approximation of the finitely extensible nonlinear elastic (FENE-P) model. At high Weissenberg number, an arc shape region with high elastic stress, which is similar to shock wave, forms in the frontal area of the blunt body. This region acts as a stationary shield to separate the flow into different regions. Thus, the free stream resembles to pass this shield instead of the original blunt body. As this shield has symmetric feature, the wake flow restores symmetry.

preprint2022arXiv

Wall-induced translation of a rotating particle in viscoelastic fluid

Shear-thinning and viscoelasticity are two non-Newtonian fluid properties widely existing in biological fluids. In this study, we found that the translation motion of a rotating particle near a wall speed up firstly, and then slows down with enhancement of fluid viscoelasticity, which is different from the behavior reported in shear thinning fluid (Chen et al. J. Fluid Mech. 2021, 927). Our research is carried out by numerical simulation of Navier-Stokes equations combined with Oldroyd-B constitutive model. This work is expected to be helpful to understand the movement of a rotating sphere near a wall in complex fluids comprehensively.

preprint2021arXiv

Marangoni Convection-Driven Laser Fountains and Waves on Free Surfaces of Liquids

It is well accepted that an outward Marangoni convection from a low surface tension region will make the surface depressed. Here, we report that this established perception is only valid for thin liquid films. Using surface laser heating, we show that in deep liquids a laser beam actually pulls up the fluid above the free surface generating fountains with different shapes. Whereas with decreasing liquid depth a transition from fountain to indentation with fountain in-indentation is observed. Further, high-speed imaging reveals a transient surface process before steady elevation is formed, and this dynamic deformation is subsequently utilized to resonantly excite giant surface waves by a modulated laser beam. Computational fluid dynamics models reveal the underlying flow patterns and quantify the depth-dependent and time-resolved surface deformations. Our discoveries and techniques have upended the century-old perception and opened up a new regime of interdisciplinary research and applications of Marangoni-induced interface phenomena and optocapillary fluidic surfaces-the control of fluids with light.

preprint2019arXiv

Broadband mid-infrared perfect absorber using fractal Gosper curve

Designing broadband metamaterial perfect absorbers is challenging due to the intrinsically narrow bandwidth of surface plasmon resonances. Here, the paper reports an ultra-broadband metamaterial absorber by using space filling Gosper curve. The optimized result shows an average absorptivity of 95.78% from 2.64 to 9.79 μm across the entire mid-infrared region. Meanwhile, the absorber shows insensitivity to the polarization angle and the incident angle of the incident light. The underlying physical principles, used in our broadband absorber, involve a fractal geometry with multiple scales and a dissipative plasmonic crystal. The broadband perfect absorption can be attributed to multiple electric resonances at different wavelengths supported by a few segments in the defined Gosper curve.

preprint2016arXiv

A minimal, "hydrogen atom" version of an inversion-breaking Weyl semimetal

The recent explosion of research interest in Weyl semimetals has led to many proposed Weyl semimetal candidates and a few experimental observations of a Weyl semimetal in real materials. Through this experience, we have come to appreciate that typical Weyl semimetals host many Weyl points. For instance, the first Weyl semimetal observed in experiment, TaAs, hosts 24 Weyl points. Similarly, the Mo$_x$W$_{1-x}$Te$_2$ series, recently under study as the first Type II Weyl semimetal, has eight Weyl points. However, it is well-understood that for a Weyl semimetal without inversion symmetry but with time-reversal symmetry, the minimum number of Weyl points is four. Realizing such a minimal Weyl semimetal is fundamentally relevant because it would offer the simplest "hydrogen atom" example of an inversion-breaking Weyl semimetal. At the same time, transport experiments and device applications may be simpler in a system with as few Weyl points as possible. Recently, TaIrTe$_4$ has been predicted to be a minimal, inversion-breaking Weyl semimetal. However, crucially, the Weyl points and Fermi arcs live entirely above the Fermi level, making them inaccessible to conventional angle-resolved photoemission spectroscopy (ARPES). Here we use pump-probe ARPES to directly access the band structure above the Fermi level in TaIrTe$_4$. We directly observe Weyl points and topological Fermi arcs, showing that TaIrTe$_4$ is a Weyl semimetal. We find that, in total, TaIrTe$_4$ has four Weyl points, providing the first example of a minimal inversion-breaking Weyl semimetal. Our results hold promise for accessing exotic transport phenomena arising in Weyl semimetals in a real material.

preprint2016arXiv

Discovery of a new type of topological Weyl fermion semimetal state in Mo$_x$W$_{1-x}$Te$_2$

The recent discovery of a Weyl semimetal in TaAs offers the first Weyl fermion observed in nature and dramatically broadens the classification of topological phases. However, in TaAs it has proven challenging to study the rich transport phenomena arising from emergent Weyl fermions. The series Mo$_x$W$_{1-x}$Te$_2$ are inversion-breaking, layered, tunable semimetals already under study as a promising platform for new electronics and recently proposed to host Type II, or strongly Lorentz-violating, Weyl fermions. Here we report the discovery of a Weyl semimetal in Mo$_x$W$_{1-x}$Te$_2$ at $x = 25\%$. We use pump-probe angle-resolved photoemission spectroscopy (pump-probe ARPES) to directly observe a topological Fermi arc above the Fermi level, demonstrating a Weyl semimetal. The excellent agreement with calculation suggests that Mo$_x$W$_{1-x}$Te$_2$ is the first Type II Weyl semimetal. We also find that certain Weyl points are at the Fermi level, making Mo$_x$W$_{1-x}$Te$_2$ a promising platform for transport and optics experiments on Weyl semimetals.

preprint2016arXiv

Large-area and high-quality 2D transition metal telluride

Atomically thin transitional metal ditellurides like WTe2 and MoTe2 have triggered tremendous research interests because of their intrinsic nontrivial band structure. They are also predicted to be 2D topological insulators and type-II Weyl semimetals. However, most of the studies on ditelluride atomic layers so far rely on the low-yield and time-consuming mechanical exfoliation method. Direct synthesis of large-scale monolayer ditellurides has not yet been achieved. Here, using the chemical vapor deposition (CVD) method, we demonstrate controlled synthesis of high-quality and atom-thin tellurides with lateral size over 300 μm. We found that the as-grown WTe2 maintains two different stacking sequences in the bilayer, where the atomic structure of the stacking boundary is revealed by scanning transmission electron microscope (STEM). The low-temperature transport measurements revealed a novel semimetal-to-insulator transition in WTe2 layers and an enhanced superconductivity in few-layer MoTe2. This work paves the way to the synthesis of atom-thin tellurides and also quantum spin Hall devices.

preprint2016arXiv

Measuring Chern numbers above the Fermi level in the Type II Weyl semimetal Mo$_x$W$_{1-x}$Te$_2$

It has recently been proposed that electronic band structures in crystals give rise to a previously overlooked type of Weyl fermion, which violates Lorentz invariance and, consequently, is forbidden in particle physics. It was further predicted that Mo$_x$W$_{1-x}$Te$_2$ may realize such a Type II Weyl fermion. One crucial challenge is that the Weyl points in Mo$_x$W$_{1-x}$Te$_2$ are predicted to lie above the Fermi level. Here, by studying a simple model for a Type II Weyl cone, we clarify the importance of accessing the unoccupied band structure to demonstrate that Mo$_x$W$_{1-x}$Te$_2$ is a Weyl semimetal. Then, we use pump-probe angle-resolved photoemission spectroscopy (pump-probe ARPES) to directly observe the unoccupied band structure of Mo$_x$W$_{1-x}$Te$_2$. For the first time, we directly access states $> 0.2$ eV above the Fermi level. By comparing our results with $\textit{ab initio}$ calculations, we conclude that we directly observe the surface state containing the topological Fermi arc. Our work opens the way to studying the unoccupied band structure as well as the time-domain relaxation dynamics of Mo$_x$W$_{1-x}$Te$_2$ and related transition metal dichalcogenides.

preprint2016arXiv

Time Resource Networks

The problem of scheduling under resource constraints is widely applicable. One prominent example is power management, in which we have a limited continuous supply of power but must schedule a number of power-consuming tasks. Such problems feature tightly coupled continuous resource constraints and continuous temporal constraints. We address such problems by introducing the Time Resource Network (TRN), an encoding for resource-constrained scheduling problems. The definition allows temporal specifications using a general family of representations derived from the Simple Temporal network, including the Simple Temporal Network with Uncertainty, and the probabilistic Simple Temporal Network (Fang et al. (2014)). We propose two algorithms for determining the consistency of a TRN: one based on Mixed Integer Programing and the other one based on Constraint Programming, which we evaluate on scheduling problems with Simple Temporal Constraints and Probabilistic Temporal Constraints.

preprint2016arXiv

Topological Weyl phase transition in Mo$_x$W$_{1-x}$Te$_2$

Topological phases of matter exhibit phase transitions between distinct topological classes. These phase transitions are exotic in that they do not fall within the traditional Ginzburg-Landau paradigm but are instead associated with changes in bulk topological invariants and associated topological surface states. In the case of a Weyl semimetal this phase transition is particularly unusual because it involves the creation of bulk chiral charges and the nucleation of topological Fermi arcs. Here we image a topological phase transition to a Weyl semimetal in Mo$_x$W$_{1-x}$Te$_2$ with changing composition $x$. Using pump-probe ultrafast angle-resolved photoemission spectroscopy (pump-probe ARPES), we directly observe the nucleation of a topological Fermi arc at $x_c \sim 7\%$, showing the critical point of a topological Weyl phase transition. For Mo dopings $x < x_c$, we observe no Fermi arc, while for $x > x_c$, the Fermi arc gradually extends as the bulk Weyl points separate. Our results demonstrate for the first time the creation of magnetic monopoles in momentum space. Our work opens the way to manipulating chiral charge and topological Fermi arcs in Weyl semimetals for transport experiments and device applications.

preprint2015arXiv

Controlled Synthesis of Organic/Inorganic van der Waals Solid for Tunable Light-matter Interactions

Van der Waals (vdW) solids, as a new type of artificial materials that consist of alternating layers bonded by weak interactions, have shed light on fascinating optoelectronic device concepts. As a result, a large variety of vdW devices have been engineered via layer-by-layer stacking of two-dimensional materials, although shadowed by the difficulties of fabrication. Alternatively, direct growth of vdW solids has proven as a scalable and swift way, highlighted by the successful synthesis of graphene/h-BN and transition metal dichalcogenides (TMDs) vertical heterostructures from controlled vapor deposition. Here, we realize high-quality organic and inorganic vdW solids, using methylammonium lead halide (CH3NH3PbI3) as the organic part (organic perovskite) and 2D inorganic monolayers as counterparts. By stacking on various 2D monolayers, the vdW solids behave dramatically different in light emission. Our studies demonstrate that h-BN monolayer is a great complement to organic perovskite for preserving its original optical properties. As a result, organic/h-BN vdW solid arrays are patterned for red light emitting. This work paves the way for designing unprecedented vdW solids with great potential for a wide spectrum of applications in optoelectronics.

preprint2015arXiv

Unoccupied electronic structure and signatures of topological Fermi arcs in the Weyl semimetal candidate Mo$_x$W$_{1-x}$Te$_2$

Weyl semimetals have sparked intense research interest, but experimental work has been limited to the TaAs family of compounds. Recently, a number of theoretical works have predicted that compounds in the Mo$_x$W$_{1-x}$Te$_2$ series are Weyl semimetals. Such proposals are particularly exciting because Mo$_x$W$_{1-x}$Te$_2$ has a quasi two-dimensional crystal structure well-suited to many transport experiments, while WTe$_2$ and MoTe$_2$ have already been the subject of numerous proposals for device applications. However, with available ARPES techniques it is challenging to demonstrate a Weyl semimetal in Mo$_x$W$_{1-x}$Te$_2$. According to the predictions, the Weyl points are above the Fermi level, the system approaches two critical points as a function of doping, there are many irrelevant bulk bands, the Fermi arcs are nearly degenerate with bulk bands and the bulk band gap is small. Here, we study Mo$_x$W$_{1-x}$Te$_2$ for $x = 0.07$ and 0.45 using pump-probe ARPES. The system exhibits a dramatic response to the pump laser and we successfully access states $> 0.2$eV above the Fermi level. For the first time, we observe direct, experimental signatures of Fermi arcs in Mo$_x$W$_{1-x}$Te$_2$, which agree well with theoretical calculations of the surface states. However, we caution that the interpretation of these features depends sensitively on free parameters in the surface state calculation. We comment on the prospect of conclusively demonstrating a Weyl semimetal in Mo$_x$W$_{1-x}$Te$_2$.