Researcher profile

Long Wang

Long Wang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

27works

0followers

24topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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.

preprint2026arXiv

A free-fall-based switching criterion for P^3 T N-body methods in collisional stellar systems

The P$^3$T scheme is a hybrid method for simulating gravitational $N$-body systems. It combines a fast particle-tree (PT) algorithm for long-range forces with a high-accuracy particle-particle (PP, direct $N$-body) solver for short-range interactions. Preserving both PT efficiency and PP accuracy requires a robust PT-PP switching criterion. We introduce a simple free-fall-based switching criterion for general stellar systems, alongside the commonly used velocity-dispersion-based ($σ$-based) criterion. Using the \textsc{petar} code with the P$^3$T scheme and slow-down algorithmic regularization for binaries and higher-order multiples, we perform extensive simulations of star clusters to evaluate how each criterion affects energy conservation and binary evolution. For systems in virial equilibrium, we find that the free-fall-based criterion is generally more accurate for low-$σ$ or loose clusters containing binaries, whereas the $σ$-based criterion is better suited for high-$σ$ systems. Under subvirial or fractal initial conditions, both criteria struggle to maintain high energy conservation; however, the free-fall-based criterion improves as the tree timestep is reduced, whereas the $σ$-based degrades due to its low-accuracy treatment of two-body encounters.

preprint2026arXiv

Collapse versus Disruption: The Fate of Compact Stellar Systems in Ultralight Dark Matter Halos

Interference of the ultralight dark matter (ULDM) field generates time-varying gravitational potential fluctuations, which stochastically heat stellar systems embedded in ULDM halos. Small-sized stellar systems are therefore often used to set stringent constraints on ULDM. However, the evolution of systems with sizes well below the ULDM de Broglie wavelength remains poorly explored. Using numerical simulations, we show that the evolution of compact stellar systems in ULDM halos is governed by the interplay between internal stellar relaxation and ULDM-induced heating. We find the following main results. First, in sufficiently compact systems, relaxation-driven core collapse dominates, allowing the system to remain bound and dense, while ULDM-induced stripping of outer stars further accelerates the collapse. Second, in more extended systems, ULDM heating dominates and ultimately disrupts the system. Near the disruption threshold, we identify systems resembling ultra-faint dwarfs like Segue 1. Third, we further introduce a dimensionless parameter to quantify the relative importance of heating and relaxation and finally lead to an evolutionary phase diagram. Our results reveal the rich and nontrivial dynamics of compact stellar systems in ULDM halos, indicating that precise system modeling is essential for robust ULDM constraints.

preprint2022arXiv

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and fractal (f2) for clusters $<100$ Myr, and halo (h) and tidal-tail (t) for clusters $>100$ Myr. The kinematical substructures of f1-type clusters are elongated; these resemble the disrupted cluster Group X. Kinematic tails are distinct in t-type clusters, especially Pleiades. We identify 29 hierarchical groups in four young regions (Alessi 20, IC 348, LP 2373, LP 2442); ten among these are new. The hierarchical groups form filament networks. Two regions (Alessi 20, LP 2373) exhibit global "orthogonal" expansion (stellar motion perpendicular to the filament), which might cause complete dispersal. Infalling-like flows (stellar motion along the filament) are found in UBC 31 and related hierarchical groups in the IC 348 region. Stellar groups in the LP 2442 region (LP 2442 gp 1-5) are spatially well-mixed but kinematically coherent. A merging process might be ongoing in the LP 2442 subgroups. For younger systems ($\lesssim30$ Myr), the mean axis ratio, cluster mass and half-mass radius tend to increase with age values. These correlations between structural parameters may imply two dynamical processes occurring in the hierarchical formation scenario in young stellar groups: (1) filament dissolution and (2) sub-group mergers.

preprint2022arXiv

A tale of caution: the tails of NGC 752 are much longer than claimed

Understanding the exact extent and content of tidal tails of open clusters provide useful clues on how field stars populate the Milky Way. We reanalyse, using Gaia EDR3 data, the tails around the open cluster NGC 752. Compared to previous analyses, we look at a much wider region around the cluster and use first the convergent point method, coupled with a clustering analysis using DBSCAN. We find that the cluster, located 433 pc away and well described by a Plummer profile, has very long and asymmetric tails, extending more than 260 pc on the sky (from tip to tip) - four times larger than previously thought - and contains twice as many stars. Numerical models computed with PETAR serve as guide and confirm our analysis. The tails follows the predictions from models, but the trailing tail appears slightly distorted, possibly indicating that the cluster had a complicated history of galactic encounters. Applying an alternative method to the newly developed compact convergent method, we potentially trace the cluster's tidal tails to their full extent, covering several thousands of parsecs and more than 1,000 stars. Our analysis therefore opens a new window on the study of open clusters, whose potential will be fully unleashed with future Gaia data releases.

preprint2022arXiv

Enhanced phase noise reduction in localized two-way optical frequency comparison

High-stability optical frequency comparison over fiber link enables the establishment of ultrastable optical clock networks, having the potential to promote a series of applications, including metrology, geodesy, and astronomy. In this article, we theoretically analyze and experimentally demonstrate a timedelayed local two-way (TD-LTW) optical frequency comparison scheme with improved comparison stability, showing that the fractional instability of optical frequency comparison over a 50- km transfer link can be reduced from $1.30\times10^{-15}$ to $5.25\times10^{-16}$ at the 1 s integration time with an improvement factor of 2.48. Additionally, we also for the first time model and experimentally verify the effect of the inhomogeneous phase noise along the fiber link on the system performance. We believe that the theory and technique proposed here will be helpful in developing the high-stability optical clock networks over large-area fiber links.

preprint2022arXiv

Generative Wind Power Curve Modeling Via Machine Vision: A Self-learning Deep Convolutional Network Based Method

This paper develops a novel self-training U-net (STU-net) based method for the automated WPC model generation without requiring data pre-processing. The self-training (ST) process of STU-net has two steps. First, different from traditional studies regarding the WPC modeling as a curve fitting problem, in this paper, we renovate the WPC modeling formulation from a machine vision aspect. To develop sufficiently diversified training samples, we synthesize supervisory control and data acquisition (SCADA) data based on a set of S-shape functions depicting WPCs. These synthesized SCADA data and WPC functions are visualized as images and paired as training samples(I_x, I_wpc). A U-net is then developed to approximate the model recovering I_wpc from I_x. The developed U-net is applied into observed SCADA data and can successfully generate the I_wpc. Moreover, we develop a pixel mapping and correction process to derive a mathematical form f_wpc representing I_wpcgenerated previously. The proposed STU-net only needs to train once and does not require any data preprocessing in applications. Numerical experiments based on 76 WTs are conducted to validate the superiority of the proposed method by benchmarking against classical WPC modeling methods. To demonstrate the repeatability of the presented research, we release our code at https://github.com/IkeYang/STU-net.

preprint2022arXiv

Gravitational wave of intermediate-mass black holes in Population III star clusters

Previous theoretical studies suggest that the Population III (Pop3) stars tend to form in extremely metal poor gas clouds with approximately $10^5 M_\odot$ embedded in mini dark matter halos. Very massive stars can form via multiple collisions in Pop3 star clusters and eventually evolve to intermediate-mass black holes (IMBHs). In this work, we conduct star-by-star $N$-body simulations for modelling the long-term evolution of Pop3 star clusters. We find that if the mini dark matter halos can survive today, these star clusters can avoid tidal disruption by the galactic environment and can efficiently produce IMBH-BH mergers among a wide range of redshift from 0 to 20. The average gravitational wave event rate is estimated to be $0.1-0.8~\mathrm{yr}^{-1} \mathrm{Gpc}^{-3}$, and approximately $40-80$ percent of the mergers occur at high redshift ($z>6$). The characteristic strain shows that a part of low-redshift mergers can be detected by LISA, TianQin, and Taiji, whereas most mergers can be covered by DECIGO and advanced LIGO/VIRGO/Kagra. Mergers with pair-instability BHs have a rate of approximately $0.01-0.15$~yr$^{-1}$~Gpc$^{-3}$, which can explain the GW190521-like events.

preprint2022arXiv

Investigating the UV-excess in star clusters with $N$-body simulations: predictions for future CSST observations

We study the origin of the UV-excess in star clusters by performing N-body simulations of six clusters with N=10k and N=100k (single stars & binary systems) and metallicities of Z=0.01, 0.001, and 0.0001, using PETAR. All models initially have a 50 percent primordial binary fraction. Using GalevNB we convert the simulated data into synthetic spectra and photometry for the China Space Station Telescope (CSST) and Hubble Space Telescope (HST). From the spectral energy distributions we identify three stellar populations that contribute to the UV-excess: (1) second asymptotic giant branch stars, which contribute to the UV flux at early times; (2) naked helium stars, and (3) white dwarfs, which are long-term contributors to the FUV spectra. Binary stars consisting of a white dwarf and a main-sequence star are cataclysmic variable (CV) candidates. The magnitude distribution of CV candidates is bimodal up to 2 Gyr. The bright CV population is particularly bright in FUV-NUV. The FUV-NUV color of our model clusters is 1-2 mag redder than the UV-excess globular clusters in M 87 and in the Milky Way. This discrepancy may be induced by helium enrichment in observed clusters. Our simulations are based on simple stellar evolution; we do not include the effects of variations in helium and light elements or multiple stellar populations. A positive radial color gradient is present in CSST NUV-y for main-sequence stars of all models with a color difference of 0.2-0.5 mag, up to 4 half-mass radii. The CSST NUV-g color correlates strongly with HST FUV-NUV for NUV-g>1 mag, with the linear relation $FUV-NUV=(1.09\pm0.12)\times(NUV-g)+(-1.01\pm0.22)$. This allows for conversion of future CSST NUV-g colors into HST FUV-NUV colors, which are sensitive to UV-excess features. We find that CSST will be able to detect UV-excess in galactic/extra-galactic star clusters with ages >200 Myr.

preprint2022arXiv

Remarks on nef and movable cones of hypersurfaces in Mori dream spaces

We investigate nef and movable cones of hypersurfaces in Mori dream spaces. The first result is: Let $Z$ be a smooth Mori dream space of dimension at least four whose extremal contractions are of fiber type of relative dimension at least two and let $X$ be a smooth ample divisor in $Z$, then $X$ is a Mori dream space as well. The second result is: Let $Z$ be a Fano manifold of dimension at least four whose extremal contractions are of fiber type and let $X$ be a smooth anti-canonical hypersurface in $Z$, which is a smooth Calabi--Yau variety, then the unique minimal model of $X$ up to isomorphism is $X$ itself, and moreover, the movable cone conjecture holds for $X$, namely, there exists a rational polyhedral cone which is a fundamental domain for the action of birational automorphisms on the effective movable cone of $X$. The third result is: Let $P:= \mathbb{P}^n \times \cdots \times \mathbb{P}^n$ be the $N$-fold self-product of the $n$-dimensional projective space. Let $X$ be a general complete intersection of $n+1$ hypersurfaces of multidegree $(1, \dots, 1)$ in $P$ with $\dim X \geq 3$. Then $X$ has only finitely many minimal models up to isomorphism, and moreover, the movable cone conjecture holds for $X$.

preprint2022arXiv

Robust optimal policies for team Markov games

In stochastic dynamic environments, team Markov games have emerged as a versatile paradigm for studying sequential decision-making problems of fully cooperative multi-agent systems. However, the optimality of the derived policies is usually sensitive to model parameters, which are typically unknown and required to be estimated from noisy data in practice. To mitigate the sensitivity of optimal policies to these uncertain parameters, we propose a robust model of team Markov games in this paper, where agents utilize robust optimization approaches to update strategies. This model extends team Markov games to the scenario of incomplete information and meanwhile provides an alternative solution concept of robust team optimality. To seek such a solution, we develop a robust iterative learning algorithm of team policies and prove its convergence. This algorithm, compared with robust dynamic programming, not only possesses a faster convergence rate, but also allows for using approximation calculations to alleviate the curse of dimensionality. Moreover, some numerical simulations are presented to demonstrate the effectiveness of the algorithm by generalizing the game model of sequential social dilemmas to uncertain scenarios.

preprint2022arXiv

SIRIUS Project. IV. The formation history of the Orion Nebula Cluster driven by clump mergers

The Orion Nebula Cluster (ONC) is an excellent example for understanding the formation of star clusters. Recent studies have shown that ONC has three distinct age populations and anisotropy in velocity dispersions, which are key characteristics for understanding the formation history of the ONC. In this study, we perform a smoothed-particle hydrodynamics/$N$-body simulation of star cluster formation from a turbulent molecular cloud. In this simulation, stellar orbits are integrated using a high-order integrator without gravitational softening; therefore, we can follow the collisional evolution of star clusters. We find that hierarchical formation causes episodic star formation that is observed in the ONC. In our simulation, star clusters evolve due to mergers of subclumps. The mergers bring cold gas with the clumps into the forming cluster. This enhances the star formation in the cluster centre. The dense cold gas in the cluster centre continues to form stars until the latest time. This explains the compact distribution of the youngest stars observed in the ONC. Subclump mergers also contribute to the anisotropy in the velocity dispersions and the formation of runaway stars. However, the anisotropy disappears within 0.5 Myr. The virial ratio of the cluster also increases after a merger due to the runaways. These results suggest that the ONC recently experienced a clump merger. We predict that most runaways originated from the ONC have already been found, but walkaways have not.

preprint2022arXiv

SIRIUS Project. V. Formation of off-center ionized bubbles associated with Orion Nebula Cluster

Massive stars born in star clusters terminate star cluster formation by ionizing the surrounding gas. This process is considered to be prevalent in young star clusters containing massive stars. The Orion Nebula is an excellent example associated with a forming star cluster including several massive stars (the Orion Nebula Cluster; ONC) and a 2-pc size H{\sc ii} region (ionized bubble) opening toward the observer; however, the other side is still covered with dense molecular gas. Recent astrometric data acquired by the Gaia satellite revealed the stellar kinematics in this region. By comparing this data with star cluster formation simulation results, we demonstrate that massive stars born in the ONC center were ejected via three-body encounters. Further, orbit analysis indicates that $θ^2$ Ori A, the second massive star in this region, was ejected from the ONC center toward the observer and is now returning to the cluster center. Such ejected massive stars can form a hole in the dense molecular cloud and contribute to the formation of the 2-pc bubble. Our results demonstrate that the dynamics of massive stars are essential for the formation of star clusters and H{\sc ii} regions that are not always centered by massive stars.

preprint2022arXiv

The Geometry of Zero-Determinant Strategies

The advent of Zero-Determinant (ZD) strategies has reshaped the study of reciprocity and cooperation in the iterated Prisoner's Dilemma games. The ramification of ZD strategies has been demonstrated through their ability to unilaterally enforce a linear relationship between their own average payoff and that of their co-player. Common practice conveniently represents this relationship by a straight line in the parametric plot of pairwise payoffs. Yet little attention has been paid to studying the actual geometry of the strategy space of all admissible ZD strategies. Here, our work offers intuitive geometric relationships between different classes of ZD strategies as well as nontrivial geometric interpretations of their specific parameterizations. Adaptive dynamics of ZD strategies further reveals the unforeseen connection between general ZD strategies and the so-called equalizers that can set any co-player's payoff to a fixed value. We show that the class of equalizers forming a hyperplane is the critical equilibrium manifold, only part of which is stable. The same hyperplane is also a separatrix of the cooperation-enhancing region where the optimum response is to increase cooperation for each of the four payoff outcomes. Our results shed light on the simple but elegant geometry of ZD strategies that is previously overlooked.

preprint2021arXiv

Passive Optical Phase Stabilization on a Ring Fiber Network

Optical frequency transfer provides the means for high-fidelity frequency transfer across thousands of kilometers. A critical step in the further developing optical frequency transfer is its capability to transfer a high spectral purity feature from ultrastable lasers or optical clocks to any remote locations and, at the same time, its adaptability to incorporate the optical frequency transfer technique into any existing communication networks with different topologies. Here we for the first time report a technique that delivers optical-frequency signals to multiple independent remote hubs along a ring optical-fiber network with passive phase stabilization. The technique automatically corrects optical-fiber length fluctuations of arbitrary hubs along the loop by mixing and shifting optical signals. Without the help of an active phase tracker and a compensator, it could significantly mitigate some technical problems such as the limited compensation speed and phase recovery time, the phase jitter contamination caused by the servo bump in conventional phase noise cancellation. Moreover, by transmitting optical signals along both directions using the same optical source, it can improve the signal-to-noise ratio at each hub. This technique maintains the same delay-limited phase noise correction capability as in conventional techniques and, furthermore, improves the phase jitter by a factor of 3, opening a way to a broad distribution of an ultrastable frequency reference with high spectral purity and enabling a wide range of applications beyond metrology over a ring fiber network with the naturally impressive reliability and scalability.

preprint2020arXiv

A slow-down time-transformed symplectic integrator for solving the few-body problem

An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which combines the time-transformed explicit symplectic integrator (Preto & Tremaine 1999; Mikkola & Tanikawa 1999) and the slow-down method (Mikkola & Aarseth 1996). The former conserves the Hamiltonian and the angular momentum for a long-term evolution, while the latter significantly reduces the computational cost for a weakly perturbed binary. In this work, the Hamilton equations of this algorithm are analyzed in detail. We mathematically and numerically show that it can correctly reproduce the secular evolution like the orbit averaged method and also well conserve the angular momentum. For a weakly perturbed binary, the method is possible to provide a few order of magnitude faster performance than the classical algorithm. A publicly available code written in the c++ language, SDAR, is available on GitHub (https://github.com/lwang-astro/SDAR). It can be used either as a stand alone tool or a library to be plugged in other $N$-body codes. The high precision of the floating point to 62 digits is also supported.

preprint2020arXiv

A theoretical framework of the scaled Gaussian stochastic process in prediction and calibration

Model calibration or data inversion is one of fundamental tasks in uncertainty quantification. In this work, we study the theoretical properties of the scaled Gaussian stochastic process (S-GaSP), to model the discrepancy between reality and imperfect mathematical models. We establish the explicit connection between Gaussian stochastic process (GaSP) and S-GaSP through the orthogonal series representation. The predictive mean estimator in the S-GaSP calibration model converges to the reality at the same rate as the GaSP with a suitable choice of the regularization and scaling parameters. We also show the calibrated mathematical model in the S-GaSP calibration converges to the one that minimizes the $L_2$ loss between the reality and mathematical model, whereas the GaSP model with other widely used covariance functions does not have this property. Numerical examples confirm the excellent finite sample performance of our approaches compared to a few recent approaches.

preprint2020arXiv

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

preprint2020arXiv

Dual Temporal Memory Network for Efficient Video Object Segmentation

Video Object Segmentation (VOS) is typically formulated in a semi-supervised setting. Given the ground-truth segmentation mask on the first frame, the task of VOS is to track and segment the single or multiple objects of interests in the rest frames of the video at the pixel level. One of the fundamental challenges in VOS is how to make the most use of the temporal information to boost the performance. We present an end-to-end network which stores short- and long-term video sequence information preceding the current frame as the temporal memories to address the temporal modeling in VOS. Our network consists of two temporal sub-networks including a short-term memory sub-network and a long-term memory sub-network. The short-term memory sub-network models the fine-grained spatial-temporal interactions between local regions across neighboring frames in video via a graph-based learning framework, which can well preserve the visual consistency of local regions over time. The long-term memory sub-network models the long-range evolution of object via a Simplified-Gated Recurrent Unit (S-GRU), making the segmentation be robust against occlusions and drift errors. In our experiments, we show that our proposed method achieves a favorable and competitive performance on three frequently-used VOS datasets, including DAVIS 2016, DAVIS 2017 and Youtube-VOS in terms of both speed and accuracy.

preprint2020arXiv

New results on EP elements in rings with involution

In this paper, we mainly give characterizations of EP elements in terms of equations. In addition, a related notion named a central EP element is defined and investigated. Finally, we focus on characterizations of a generalized EP element, i.e. *-DMP element.

preprint2020arXiv

Online distributed algorithms for seeking generalized Nash equilibria in dynamic environments

In this paper, we study the distributed generalized Nash equilibrium seeking problem of non-cooperative games in dynamic environments. Each player in the game aims to minimize its own time-varying cost function subject to a local action set. The action sets of all players are coupled through a shared convex inequality constraint. Each player can only have access to its own cost function, its own set constraint and a local block of the inequality constraint, and can only communicate with its neighbours via a connected graph. Moreover, players do not have prior knowledge of their future cost functions. To address this problem, an online distributed algorithm is proposed based on consensus algorithms and a primal-dual strategy. Performance of the algorithm is measured by using dynamic regrets. Under mild assumptions on graphs and cost functions, we prove that if the deviation of variational generalized Nash equilibrium sequence increases within a certain rate, then the regrets, as well as the violation of inequality constraint, grow sublinearly. A simulation is presented to demonstrate the effectiveness of our theoretical results.

preprint2020arXiv

PeTar: a high-performance N-body code for modeling massive collisional stellar systems

The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been performed by using the NBODY6++GPU code. Such models took half a year computational time on a GPU based super-computer. In this work, we develop a new N-body code, PeTar, by combining the methods of Barnes-Hut tree, Hermite integrator and slow-down algorithmic regularization (SDAR). The code can accurately handle an arbitrary fraction of multiple systems (e.g. binaries, triples) while keeping a high performance by using the hybrid parallelization methods with MPI, OpenMP, SIMD instructions and GPU. A few benchmarks indicate that PeTar and NBODY6++GPU have a very good agreement on the long-term evolution of the global structure, binary orbits and escapers. On a highly configured GPU desktop computer, the performance of a million-body simulation with all stars in binaries by using PeTar is 11 times faster than that of NBODY6++GPU. Moreover, on the Cray XC50 supercomputer, PeTar well scales when number of cores increase. The ten million-body problem, which covers the region of ultra compact dwarfs and nuclearstar clusters, becomes possible to be solved.

preprint2020arXiv

Step-size effect in the time-transformed leapfrog integrator on elliptic and hyperbolic orbits

A drift-kick-drift (DKD) type leapfrog symplectic integrator applied for a time-transformed separable Hamiltonian (or time-transformed symplectic integrator; TSI) has been known to conserve the Kepler orbit exactly. We find that for an elliptic orbit, such feature appears for an arbitrary step size. But it is not the case for a hyperbolic orbit: when the half step size is larger than the conjugate momenta of the mean anomaly, a phase transition happens and the new position jumps to the nonphysical counterpart of the hyperbolic trajectory. Once it happens, the energy conservation is broken. Instead, the kinetic energy minus the potential energy becomes a new conserved quantity. We provide a mathematical explanation for such phenomenon. Our result provides a deeper understanding of the TSI method, and a useful constraint of the step size when the TSI method is used to solve the hyperbolic encounters. This is particular important when an (Bulirsch-Stoer) extrapolation integrator is used together, which requires the convergence of integration errors.

preprint2020arXiv

Underactuation Design for Tendon-driven Hands via Optimization of Mechanically Realizable Manifolds in Posture and Torque Spaces

Grasp synergies represent a useful idea to reduce grasping complexity without compromising versatility. Synergies describe coordination patterns between joints, either in terms of position (joint angles) or effort (joint torques). In both of these cases, a grasp synergy can be represented as a low-dimensional manifold lying in the high-dimensional joint posture or torque space. In this paper, we use the term \textit{Mechanically Realizable Manifolds} to refer to the subset of such manifolds (in either posture or torque space) that can be achieved via mechanical coupling of the joints in underactuated hands. We present a method to optimize the design parameters of an underactuated hand in order to shape the Mechanically Realizable Manifolds to fit a pre-defined set of desired grasps. Our method guarantees that the resulting synergies can be physically implemented in an underactuated hand, and will enable the resulting hand to both reach the desired grasp postures and achieve quasistatic equilibrium while loading the grasps. We demonstrate this method on three concrete design examples motivated by a real use case, and evaluate and compare their performance in practice.

preprint2020arXiv

Very high redshift quasars and the rapid emergence of super-massive black holes

The observation of quasars at very high redshift such as Poniuaena is a challenge for models of super-massive black hole (SMBH) formation. This work presents a study of SMBH formation via known physical processes in star-burst clusters formed at the onset of the formation of their hosting galaxy. While at the early stages hyper-massive star-burst clusters reach the luminosities of quasars, once their massive stars die, the ensuing gas accretion from the still forming host galaxy compresses its stellar black hole (BH) component to a compact state overcoming heating from the BH--BH binaries such that the cluster collapses, forming a massive SMBH-seed within about a hundred Myr. Within this scenario the SMBH--spheroid correlation emerges near-to-exactly. The highest-redshift quasars may thus be hyper-massive star-burst clusters or young ultra-compact dwarf galaxies (UCDs), being the precursors of the SMBHs that form therein within about 200 Myr of the first stars. For spheroid masses <10^9.6 Msun a SMBH cannot form and instead only the accumulated nuclear cluster remains. The number evolution of the quasar phases with redshift is calculated and the possible problem of missing quasars at very high redshift is raised. SMBH-bearing UCDs and the formation of spheroids are discussed critically in view of the high redshift observations. A possible tension is found between the high star-formation rates (SFRs) implied by downsizing and the observed SFRs, which may be alleviated within the IGIMF theory and if the downsizing times are somewhat longer.

preprint2019arXiv

Accelerated FDPS --- Algorithms to Use Accelerators with FDPS

In this paper, we describe the algorithms we implemented in FDPS to make efficient use of accelerator hardware such as GPGPUs. We have developed FDPS to make it possible for many researchers to develop their own high-performance parallel particle-based simulation programs without spending large amount of time for parallelization and performance tuning. The basic idea of FDPS is to provide a high-performance implementation of parallel algorithms for particle-based simulations in a "generic" form, so that researchers can define their own particle data structure and interparticle interaction functions and supply them to FDPS. FDPS compiled with user-supplied data type and interaction function provides all necessary functions for parallelization, and using those functions researchers can write their programs as though they are writing simple non-parallel program. It has been possible to use accelerators with FDPS, by writing the interaction function that uses the accelerator. However, the efficiency was limited by the latency and bandwidth of communication between the CPU and the accelerator and also by the mismatch between the available degree of parallelism of the interaction function and that of the hardware parallelism. We have modified the interface of user-provided interaction function so that accelerators are more efficiently used. We also implemented new techniques which reduce the amount of work on the side of CPU and amount of communication between CPU and accelerators. We have measured the performance of N-body simulations on a systems with NVIDIA Volta GPGPU using FDPS and the achieved performance is around 27 \% of the theoretical peak limit. We have constructed a detailed performance model, and found that the current implementation can achieve good performance on systems with much smaller memory and communication bandwidth.

preprint2019arXiv

Evolutionary dynamics with game transitions

The environment has a strong influence on a population's evolutionary dynamics. Driven by both intrinsic and external factors, the environment is subject to continual change in nature. To capture an ever-changing environment, we consider a model of evolutionary dynamics with game transitions, where individuals' behaviors together with the games they play in one time step influence the games to be played next time step. Within this model, we study the evolution of cooperation in structured populations and find a simple rule: weak selection favors cooperation over defection if the ratio of the benefit provided by an altruistic behavior, $b$, to the corresponding cost, $c$, exceeds $k-k'$, where $k$ is the average number of neighbors of an individual and $k'$ captures the effects of the game transitions. Even if cooperation cannot be favored in each individual game, allowing for a transition to a relatively valuable game after mutual cooperation and to a less valuable game after defection can result in a favorable outcome for cooperation. In particular, small variations in different games being played can promote cooperation markedly. Our results suggest that simple game transitions can serve as a mechanism for supporting prosocial behaviors in highly-connected populations.

preprint2019arXiv

MOCCA-SURVEY Database I: Dissolution of tidally filling star clusters harboring black hole subsystem

We investigate the dissolution process of star clusters embedded in an external tidal field and harboring a subsystem of stellar-mass black hole. For this purpose we analyzed the MOCCA models of real star clusters contained in the Mocca Survey Database I. We showed that the presence of a stellar-mass black hole subsystem in tidally filling star cluster can lead to abrupt cluster dissolution connected with the loss of cluster dynamical equilibrium. Such cluster dissolution can be regarded as a third type of cluster dissolution mechanism. We additionally argue that such a mechanism should also work for tidally under-filling clusters with a top-heavy initial mass function.

Long Wang

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

27 published item(s)

A free-fall-based switching criterion for P^3 T N-body methods in collisional stellar systems

Collapse versus Disruption: The Fate of Compact Stellar Systems in Ultralight Dark Matter Halos

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

A tale of caution: the tails of NGC 752 are much longer than claimed

Enhanced phase noise reduction in localized two-way optical frequency comparison

Generative Wind Power Curve Modeling Via Machine Vision: A Self-learning Deep Convolutional Network Based Method

Gravitational wave of intermediate-mass black holes in Population III star clusters

Investigating the UV-excess in star clusters with $N$-body simulations: predictions for future CSST observations

Remarks on nef and movable cones of hypersurfaces in Mori dream spaces

Robust optimal policies for team Markov games

SIRIUS Project. IV. The formation history of the Orion Nebula Cluster driven by clump mergers

SIRIUS Project. V. Formation of off-center ionized bubbles associated with Orion Nebula Cluster

The Geometry of Zero-Determinant Strategies

Passive Optical Phase Stabilization on a Ring Fiber Network

A slow-down time-transformed symplectic integrator for solving the few-body problem

A theoretical framework of the scaled Gaussian stochastic process in prediction and calibration

BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

Dual Temporal Memory Network for Efficient Video Object Segmentation

New results on EP elements in rings with involution

Online distributed algorithms for seeking generalized Nash equilibria in dynamic environments

PeTar: a high-performance N-body code for modeling massive collisional stellar systems

Step-size effect in the time-transformed leapfrog integrator on elliptic and hyperbolic orbits

Underactuation Design for Tendon-driven Hands via Optimization of Mechanically Realizable Manifolds in Posture and Torque Spaces

Very high redshift quasars and the rapid emergence of super-massive black holes

Accelerated FDPS --- Algorithms to Use Accelerators with FDPS

Evolutionary dynamics with game transitions

MOCCA-SURVEY Database I: Dissolution of tidally filling star clusters harboring black hole subsystem