An Asymmetric Popularity-Similarity Optimization Method for Embedding Directed Networks into Hyperbolic Space

Network embedding is a frontier topic in current network science. The scale-free property of complex networks can emerge as a consequence of the exponential expansion of hyperbolic space. Some embedding models have recently been developed to explore hyperbolic geometric properties of complex networks—in particular, symmetric networks. Here, we propose a model for embedding directed networks into hyperbolic space. In accordance with the bipartite structure of directed networks and multiplex node information, the method replays the generation law of asymmetric networks in hyperbolic space, estimating the hyperbolic coordinates of each node in a directed network by the asymmetric popularity-similarity optimization method in the model. Additionally, the experiments in several real networks show that our embedding algorithm has stability and that the model enlarges the application scope of existing methods.


Introduction
Complex networks can largely simplify real systems and preserve the essential information of the interaction structure. It thus is an ideal tool for investigating complex systems. However, complex networks are models with nongeometric properties, which include a large set of tools and methodologies developed in geometry that cannot be applied to complex networks. In this context, there is a wave of studies exploring geometric properties of complex networks, aiming at mapping complex networks into latent variables (that is, hidden variables and Gaussian latent variables) or a low-dimensional metric space (that is, Euclidean space or hyperbolic space) [1][2][3].
Advances in network geometry have shown that structural properties observed in scale-free networks derived from real complex systems can emerge as the geometrical properties [4,5]. Hyperbolic geometry is a branch of non-Euclidean geometry, and it has many applications in practical engineering techniques. More importantly, the random geometric model and hyperbolic embedding model of growth [3,6] proposed later can easily explain the heterogeneity and high clustering of scale-free networks and even give the clear meaning of each coordinate. ese models can not only simulate the growth of networks but also explain the dynamic process of the classic BA network model in complex networks. Related studies are thus becoming increasingly popular, and network models with geometric properties have been used successfully in many fields in network science and other disciplines, including brain science [7,8], international trade [9], route transfer [10][11][12], and protein formation mechanisms [13,14].
When this research framework was proposed, it attracted wide attention from scholars, and then models and hyperbolic spatial embedding methods were developed for it [3,6,[15][16][17].
ese models have a great performance in studying the potential structure of networks. In particular, in the popularity-similarity optimization method, the hyperbolic property of complex networks is combined with the hidden space property, and the spatial position of nodes is estimated by statistical inference. However, these models cannot describe real systems completely. One drawback is that those models ignore the directionality of links. e relationships between nodes may be unequal in most real networks, which is the so-called asymmetry property of links. Although the asymmetry property may bring many challenges for network embedding, ignoring the asymmetry of directed networks will lose much important information and cannot lead to a full representation of the structure and function of real systems.
Another branch of study is network representation learning or graph embedding through machine learning methods and matrix analysis [1,[18][19][20][21][22], including directed networks [23,24]. ose models have provided many useful modeling inspirations for us, such as motifs, random walks, centrality information with graph convolutions [25], the high-order proximity [26], and the spatial and temporal proximity [27]. Unfortunately, an obvious limitation is that these models cannot identify the real meaning of the spatial coordinates when obtaining a vector representation of the nodes. More importantly, since it has been demonstrated that a scale-free degree distribution is a basic condition for the embedding of hyperbolic space [3], it is urgent to develop a directed network embedding model in hyperbolic space.
To further address this problem, we explore the intrinsic relationship between directed links and the network topology. Interestingly, in a nontrivial way, directed networks have a hidden bipartite structure [28]. We checked that all kinds of directed networks have such structures, and this phenomenon is the universal law of directed networks. e contribution of our work is that we offer a directed network embedding scheme based on node information multiplexing and identifying a potential topological structure (bipartite structure) as well as a new idea for the dimensionality reduction of directed network data. In addition, using visualization technology, we provide a new snapshot of directed networks in the hyperbolic space, which enables us to show the nodes' status and macrolevel structure and features.
In this paper, we propose an asymmetric popularitysimilarity optimization method for embedding directed networks into hyperbolic space through the bipartite structure of directed networks, and our methodology is grounded in the network topology information and its characteristics. From this perspective, we first introduce the concept of the bipartite structure of directed networks firstly. Based on this, the mapping model for mapping nodes from directed networks to hyperbolic space is discussed in Section 2. Four real-world directed networks are used to test the applicability of our mapping method in Section 3. Section 4 closes the article with concluding remarks.

Data.
In the simulation, we use four empirical network datasets and give the full description of the data source, including the field of Caenorhabditis elegans (C. elegans) neural systems, the international trade, email relationships, and the international migration. (1) Email dataset is the internal email communication network between employees of a mid-sized manufacturing company [29,30]. e network is directed and nodes represent employees and edges between two nodes are individual emails. (2) International migration dataset is a weighted directed network with 153 nodes and can be required from the world bank (http://www.worldbank.org/). Adjacency matrix a ij describes the immigration from the country i to j with the migration flow w ij . (3) C. elegans neural dataset, a directed network, describes the neural interconnection via chemical synapses and gap junctions, which can be obtained from the Wormatlas database [31]. (4) e international trade dataset describes the trade relationships among countries (or regions), which are obtained from the United Nations commodity trade statistics database (https://comtrade.un.org/). If country i imports to country j with weights w, adjacency matrix is a ij � w.

Complex Network and Hyperbolic Space.
Hyperbolic space is an isotropic space with negative curvature that cannot be embedded in any Euclidean space. e topological geometry and hyperbolic geometry of complex networks are intimately related, which has been well explored in mathematics [32], because the shortest paths in networks, those defining chemical distances, closely follow their hyperbolic geodesics in the latent space. Poincaré disks with radius R � 2ln (N/πμκ 2 0 ) can be used to represent hyperbolic spaces according to Ref. [3]. e main property of hyperbolic geometry is the exponential expansion of space: the area A(R) of a two-dimensional hyperbolic disc grows with R as A(R) ∼ e R . From the perspective of hyperbolic geometry, the emergence of scale-free networks produces two exponentials: node density exponentially increases with the distance r from the center of the disk, and the average degree exponentially declines with r. In early studies, the hyperbolic geometric model and hidden geometric models perform equivalently for embedding scale-free network topology into a metric space. A hidden geometric model assigns to nodes their expected degrees κ and θ, where κ is derived from the power- is the minimum expected degree of node, k is the average degree of the network, and c is the power exponent. θ is chosen uniformly at random from [0, 2π]. en, each pair of nodes, with probability f(χ) � (1 + χ β ) − 1 , is connected, where the effective distance χ ij � d ij /μκ i κ j � ΔθN/2μπκ i κ j , where the angular distance d ij � ΔθN/(2π), and the parameters β and μ are constrained by μ � β sin(π/β)/2πk.
To build a hyperbolic geometry model for complex networks, −1 is set as the curvature of the hyperbolic space. Poincaré model assigns to nodes their radial r coordinates with density ρ(r) � α sin hαr/cos hαR − 1 and angular coordinates θ with density ρ(θ) � 1/2π. e geodesic distance x between nodes (r i , θ i ) and (r j , θ j ) satisfies cos hx � cos hr i cos hr j − sin hr i sin hr j cos Δθ. Nodes are connected by probability f(x) � (1 + e (x− R)β/2 ) − 1 . e process of constructing complex networks reflects the competition of two forces in this model: popularity (r) and similarity (θ).
In addition, note that the hidden metric space model and hyperbolic geometric model of a network perform equivalently through the κ − r transformation: r � R − 2 ln(κ/κ 0 ).

Interaction between Links with
Direction and a bipartite Structure. Directed links, as important linking features of the network mode, are increasingly used to incorporate the dynamics of the evolution and node status of real systems. e asymmetry property increases the difficulty of network embedding. e asymmetry of links can be identified in the topological information of complex systems, that is, a bipartite structure [33]. Different from directed networks, each node is split into two nonoverlapping groups, according to its unique feature, in bipartite structures. Modifying modeling methods of directed networks, therefore, is critical in overcoming difficulties.
Nodes in directed networks can be split into two parts, and each side of a link comes from different sets. Such a way of modeling is called the directed network with a bipartite structure, as shown in Figures 1(a) and 1(b). Specifically, each node v i in a directed network contains a i in set A and b i in set B in a bipartite structure. A directed link goes from v i to v j (v i ⟶ v j ) and can be mapped by (a i ⟶ b j ). By doing so, the number of nodes becomes twice as large especially, but this does not affect subsequent work, and a directed network will be converted to a bipartite network. Additionally, the method reconstructs geometric directed networks through dividing nodes into two categories, but these two categories are not independent completely and its intrinsic relationships will be considered by the parameter β during the modeling.
Since the bipartite structure has been discovered to be a hidden geometric property in the latent metrics space [16], we also expect the origin node and the end node to be able to form directed links under the multiscale node multiplexing perspective. Relationships between direction and topology provide a powerful idea for spatial mapping and directed link prediction. Unlike in the bipartite network, node i in two sets (a i , b i ) is one-to-one and they represent a certain kind of attribute of node i. For directed networks, for example, in the case of international trade networks, countries with stronger export capabilities are more likely to become trading partners with ones that have a high level of import. To show this, we develop a new method for displaying the asymmetry of the adjacency matrix of networks in hyperbolic space with a hybrid of direction and topology.

A Geometric Model of a Directed Network in Hyperbolic
Space.
e directed network embedding model describes how generative geometric directed networks are embedded in hyperbolic space. e model multiplexes node information as an embedding foundation by identifying the bipartite structure of directed networks and considers the trade-off between node popularity, represented by the radial coordinate, and similarity, represented by the angular coordinate distance, to be the definition of the connection probability.
In the binary network embedding model, the popularity is higher and the probability is greater between nodes with greater similarity. Unlike binary networks, directed networks should be linked by trading off four types of forces, namely, the similarity and popularity in each set (set A or set B). Radial coordinates (r ai and r bi ) and angular coordinates (θ ai and θ bi ) represent the popularity and similarity of nodes, respectively. e connection of a pair of nodes in directed networks would balance four forces: out-popularity r ai , outsimilarity θ bi , in-popularity r bj , and in-similarity θ aj . e flow diagram of the directed network embedding algorithm is shown in Figure 1.
Computational implementations of the directed embedding model are as follows.
Step 1. Obtaining model parameters from directed networks with bipartite structures. e scale-free property is ubiquitously observed in realworld networks, and hyperbolic geometry captures such features of complex networks. According to the principle that nodes multiplex information, basic properties can be computed: the average degree k * , the degree distribution (P(k * ) ∼ k −c * * ), and exponent of the power-law distribution c * in set * ∈ A, B { }. However, not all of the real-world networks have the directed scale-free property. Besides, most networks are heterogeneous in node degrees and weights, and those heterogeneities enable us to address the problem: finding the backbone structure of the complex network by filtering links [17,34] to capture the power-law property in directed networks. e method of constructing backbone networks in the paper is that the link with the smallest weight will be deleted in turns and stopping by the threshold. To find the appropriate value of the threshold, we plot the fraction of remaining nodes in the backbone N B /N 0 vs. the fraction of remaining links L B /L 0 for each deleting step. According to Ref. [9], the best choice of threshold is the point in the plane that maximizes the vertical distance to the diagonal, so that most nodes remain in the system and most links with small weight will be removed.
Another significant parameter of our model is β ∈ (1, +∞). β controls the clustering property of complex networks and is intractable due to its nontrivial dependence on the topological structure [16]. Estimating parameter in the symmetric network embedding model is performed by comparing the real networks to synthetic networks generated with the model using different values for several topological properties [9,11,35]. However, this method fails to compute the values of directed networks in experiments. Note that the embedding model splits each node into two sets, but these two sets should intuitively be dependent to a certain extent. erefore, the β in the directed network embedding model needs to take into account the clustering property and intrinsic relationships between two set nodes.
Common neighbors meet the conditions discussed above. For one thing, the clustering definition in bipartite structures derived from common neighbors to enable representing an intrinsic relationship between two sets, such as the quadrilateral-based clustering coefficient and the 4-loop density [5]. For another, empirical results indicate that the higher values of β favor connections at smaller angular distances, and the number of common neighbors m grows asymptotically as a power-law function of angular similarity: m ∼ Δθ − β [16]. Here, we use the number of common neighbors as the coarse-grained representation of the angle similarity to estimate β (that is, We, finally, take the area under the curve (AUC) as a function of β * (namely, AUC(β * )) to determine the parameter by the performance of mapping effect in our model. Specifically, for international trade networks, international migration networks, and email networks, β � 3.5; the C. elegans networks use β � 3.42.
Geometric models-hyperbolic geometric models and hidden geometric models-have successfully captured the natural geometry underlying real complex networks. Similarly, a directed network embedding into hyperbolic space also has two equivalent models. Next, two geometric models of directed networks will be introduced.
In the hyperbolic model of directed networks, scale-free directed networks are generated by scattering N nodes randomly into a hyperbolic disk of radius R. Each node v i is assigned radial coordinates r and angular coordinates θ: v ai (r ai , θ ai ) are in set A, and v bi (r bi , θ bi ) are in set B. e connection probability (p ij ) of the directed link v ij from a i to b j in the hyperbolic space represented by the Poincaré disk is defined by the Fermi-Dirac distribution: f(x ai,bj ) � (1 + e β(x ai,bj − R)/2 ) − 1 . e hyperbolic distance x ai,bj between a i and b j is given by the hyperbolic cosine function: cos hx ai,bj � cos hr ai · cos hr bj − sin hr ai · sin hr bj cos Δθ. Note that hyperbolic distance can be well approximated by r ai + r bj + 2ln(d ai,bj /2).
In the hidden geometric model of directed networks, connecting probability between nodes is any integral function: where the effective distance between the nodes χ � d ai,bj /μκ ai κ bj , and μ � β sin(π/β)/2πk * is parameter of model. κ ai and κ bi represent the expect degree of nodes in set A and set B, respectively. Expect degree κ * i for every node is drawn from the power-law distribution: Obtaining the expected degree (κ) by maximum likelihood estimation (MLE) Obtaining similarity (θ) by localized metropolishastings (LMH)

Row data Basic property of networks
Average degree: 〈k A 〉 Degree distribution: Set A: Power exponent: γ B Set B: Set A: origin node Set B: origin node Bipartite structure

Multiplex node information
Directed network e asymmetric popularity (r)-similarity (θ) method Assigns to node i expected degree and angular coordinate (κ ai, κ bj , θ ai , θ bj )
Step 3. Spatial positions and parameter estimation.
Given a snapshot of a real directed network consisting of N nodes, we use the hidden geometric model to represent a directed network. An asymmetric popularity-similarity optimization method is proposed by us to compute the radial (popularity) r * i and angular (similarity) θ * i coordinates for each node i. In particular, the embedding process of asymmetric links aims to find the coordinates of every node such that the likelihood that the given the resulting directed network topology is generated by the model described above is maximal.
Form a statistical perspective, inferring coordinates is to find the best match the hyperbolic model by a given adjacency matrix. e maximum likelihood estimation method has been widely applied to infer coordinates in hyperbolic embedding [9,11,36]. e hidden variables take particular values {κ * i , θ * i } in the network and can be found from the observed adjacency matrix a ij using Bayes' rule as follows: where . e prior probability of the hidden variables is given by and the likelihood is where a ij is adjacency matrix; if there is a link between i and j, a ij � 1, otherwise a ij � 0. e hidden variables are {κ * i , θ * i } and the effective distance χ � NkΔθ/β sin(π/β)κ ai κ bi . e obtained hidden coordinates are optimized by maximizing the likelihood in equation (3) and its logarithm: where the constant C is independent of κ * i and θ * i .
Next, node positions can be obtained by the likelihood function L. Firstly, inferring the radial coordinates is relatively easy. We derive the analytical solutions by partial derivation of the equation (6) with respect to the expected degree κ * i : By doing this, we can obtain the parameter κ * i (the expected degree of node i) from the expected degree κ ∼ * i as follows: κ ∼ * i � k * i − c * i /β and κ * i � max k * i − c * i /β, κ * 0 }. Similar to r, the angular coordinates are obtained by maximizing the likelihood function equation (6). However, it is difficult to obtain an analytical solution, and numerical methods have therefore been used to estimate angular coordinates, including the standard Metropolis-Hastings algorithm (SMH) [37] and the localized Metropolis-Hastings algorithm (LMH) [11]. Since the LMH method performs well in estimating and computes angular coordinates exactly it in a distributed manner without knowing the global network topology [11], we applied the LMH method to infer angular coordinates. e global log-likelihood equation (5) is represented by the following equation: In the LMH method, the local contribution lnL i of every node i to the global log-likelihood lnL will be defined firstly as ln L i � i≠j a ij ln p(χ) + i≠j (1 − a ij )ln[1 − p(χ)], where ln L � 1/2lnL i . Next, nodes are visited one by one. When the particular node i is visited, all other nodes have fixed the positions, and the angular position of i is moved according to the fitness, which maximizes the local log-likelihood at each node visit. e angular position is sampled at intervals with Δθ � 1/N.

Validating the Asymmetric Popularity-Similarity Optimization Method.
To assess how well the asymmetric popularity-similarity optimization (A-PSO) method performs, we examine the embedding accuracy by comparing the topology inferred by the A-PSO method to real-world directed networks. Several experiments are conducted to analyze the effectiveness of the following methods: the C. elegans neural network, the international trade network, the email network, and the international migration network. Additionally, we analyze the robustness and performance of the embedding model by repeated experiments, and the results are stable. We take two measures: (1) comparing the original network to the topological structure, i.e., in terms of the degree distribution, clustering coefficient distribution, and betweenness centrality distribution; and (2) performing a global test: linking the empirical connection probability Complexity with data to the theoretical prediction with hyperbolic distance equation (1). e first test to evaluate embedding accuracy is that we experiment and compare the first-order neighbors of nodes and the higher-order cases, including degree, clustering coefficient, and betweenness centrality, according to standard practice. We represent degree cumulative distribution, clustering cumulative distribution, and betweenness cumulative distribution in Figure 2 and observe a good match between the properties of the synthetic networks constructed by our models and real directed networks. e second check is that we compute the connection probability from empirical network data and compare it with the theoretical prediction given by equation (1). e observed closeness of the empirical and theoretical connection probabilities in Figure 3 suggests that hyperbolic metric spaces are reasonable representations of the directed network. In addition, we also find the area under the curve (AUC), which is widely used in the field of link prediction as an accuracy evaluation index. AUC statistics measure how well our model can reproduce networks in hyperbolic space, and the best possible result corresponds to AUC � 1. From Table 1, the embedding results with high AUC suggest that directed networks can be reproduced in hyperbolic space, especially socioeconomic networks, whose AUC values are well above 0.88: the international trade network and the email network.

Linking the Hyperbolic Property and the Empirical
Network. Using the A-PSO method for embedding directed networks into hyperbolic space in the above sections, we showed the association between the theoretical model and empirical networks. In this section, we first provide comparisons between hyperbolic properties (popularity, distance, and the core-periphery structure) and economic measures to show how they correlate monetary macroeconomic indicators and represent economic systems and their time evolution in hyperbolic space. Next, we analyze the relationship between disk partitions and real-function areas by comparing different measures of distance (similarity, real distance, and hyperbolic distance) in the C. elegans neural system.

Economic Systems in Hyperbolic Space.
e international trade network is a complex system of states and their trade relations, and the aim of examining it is to understand international trade in terms of quantities and mechanisms [38][39][40]. e three key characteristics of the international trading system-globalization, stratification, and localization-lay the foundation for the geometric interpretation of the international trading system [9]. From the perspective of the directed network generation mechanism, closer countries have a greater chance of becoming connected by national status and similarity in the underlying trade space. rough statistical inference techniques and the structural information of networks, a directed international trade network can be embedded in hyperbolic space. Such an embedding model uses two dimensions of national information for each country, import popularity and export popularity. Extracting a higher-dimensional national status measure can help to analyze how the multitrade relations change the world trade and even reshape the world trading system. But the international trade networks often have not the directed scale-free property; we use the filtering links method mentioned in Step 1 to address a problem, as shown in Figure 4(d).
To study this issue further, we apply our model to visualize international trade systems and infer the popularity and similarity properties (r ai , r bj , θ ai , θ bj ) of nodes, according to the data from 20 years (1996-2016) of the international trade. Recently, studies have discussed the correlation of radial coordinates as equivalent to correlations among node degrees and the national economic scale [38,39]. Similarly, radial coordinates and GDP have a significant negative correlation (approximately −0.5), which indicates that the spatial position of countries can be used to identify the national economic size, especially the country's export capacity (the correlation coefficient of imports and GDP is higher than that of exports). at is, export status can be used as a coarse-grained measure of the national economy size. Moreover, the correlation coefficient between angular distance and geographical distance is not significant (approximately 0.1), which indicates that the international trade network does not have geographical clustering in the hyperbolic space. e core-periphery structure of hyperbolic network shows that only approximately 35% of its nodes are in the center position of the hyperbolic space. e long-term evolution of the international trade system based on the network hyperbolic embedding method is shown in Figure 3(c). From the perspective of exporting ability, it presents world energy commodity trade as an imbalanced, diversified, and multipolar development. e United States and Russia have always occupied a central position, reflecting the fact that energy resources are the decisive factor in exporting capacity. In particular, with the depletion of the North Sea oil field, the UK's export position in energy trade has been gradually marginalized. Interestingly, Saudi Arabia's energy export status has gradually been marginalized due to the change in the direction of Saudi energy policy. Asia, Africa, and the European continent have become active areas and injected new vitality into the energy trade market. e European Community, China, and India have moved to central positions following the increase of their importing dominance. While India has moved towards a more central position of superpower status during the last few years, the United States, which is a leader among trade superpowers, has been at the core of international trade.

Complex Ecosystems in Hyperbolic Space.
Caenorhabditis elegans is a soil-dwelling nematode that is evolutionarily rudimentary. It contains approximately 300 neurons, and neural interconnections are made via chemical synapses and gap junctions. Despite a century of investigation, knowledge of nematode neuronal networks is incomplete [41,42]. Here, we use a directed hyperbolic network embedded framework and the potential geometry    relative position between the neurons, and the geometric distances are hyperbolic distance and angular distance. However, hyperbolic distance in our model is asymmetric, and it is difficult for some real tasks, such as node clustering analysis and node centrality computation. To further compare distances in different metric spaces, we define the symmetric distance as follows: x ij � x ai,bj + x bi,aj . e results are shown in Figures 5(a) and 5(b). From the results, we find that the angular distance is similar to the position distance. e difference between hyperbolic distances and position distance, on the contrary, is larger. is shows not only that the effective distance between neurons includes the distances of other dimensions than position distance but also that the effective distance of the nervous system is the result of a nonlinear hybrid of topological information and spatial information.

Complexity
Realistic networks can reproduce the properties of clustering, small-worldness, scale-freeness, and rich-club in the hyperbolic space. Another important feature of complex networks that is commonly observed is the community structure, in which the links of the inner community are dense, and the links between communities are sparse [43]. Since the connection probability is a decreasing function of the hyperbolic distance, there are no angular regions containing a cluster of spatially close nodes that are more densely connected to each other than they are to the rest of the network. Nevertheless, the angular distance of the hyperbolic disc indicates the similarity of nodes, and the partition of the disc is the potential module defined in geometric space [13]. To obtain the community structure of neural networks in the hyperbolic space after the nodes have acquired settled coordinates, the best partition of the Poincaré disc zoning is defined as the proportion that minimizes the average distance from one partition to the rest.
Since there are 10 neural functional areas, we divide the hyperbolic disc into 10 parts of 36 degrees each and compare them with real neural functional areas. Note that the angular coordinates of nodes are distributed in [0, 360 ∘ ] and that a key factor is the difference in initial position selection, which  Figure 5: e result of mapping the C elegans neural network in hyperbolic space. (a, b) e relationships among the hyperbolic distance, the angular distance, and the spatial distance. e abscissa indicates the distance between the nodes, and the ordinate is the statistic for the number of distances. (c) e value of η according to the starting point of different partitions, which is used to find the optimal disc partitioning result of η. e blue and pink lines represent the results of η, which are calculated from the angular distance and the hyperbolic distance, respectively. may change the effect of the partition. To solve this problem, we define the index η, which represents the ratio of the distances inside the community (Poincaré disc zoning) to the external distances (the rest of the zoning), to describe the optimal partition, η, as follows: where n is the number of partition and C i represents the average distance between nodes inside and outside the community. e range of η is [0, 1]. η is closer to 1, the better the community division.
In doing so, we obtain the optimal partition of hyperbolic space by selecting different starting partition positions, taking the minimum distance between nodes as the reference for the initial partition span and maximizing η. We further considered hyperbolic distance and angular distance as the basis for calculating η results, as shown in Figure 5(c). Due to the strong spatial position dependence of neuron function, the common information entropy of the community structure and the functional partition is very low. An interesting finding is that some neurons with higher rankings of out-directed popularity belong to the same neurological functional area, the lateral ganglion, such as RIAL, RIAR, SAAVR, RMDR, and SMDVR, but the partitions in the hyperbolic space are more scattered. is indicates that there are a large number of long-range links among the neurons of C. elegans; for example, the neuron RIAL, in the front, and VD12, in the trailing regions, have similar angular distances in the hyperbolic disc, which is also the reason why the topological partition is not consistent with the functional partition. Additionally, the functional classification of the hyperbolic embedding is a complement to the functional partitioning.

Discussion
In this work, we developed a mapping model of directed networks in hyperbolic space and highlighted the bipartite structure of directed networks. We especially focused on two main issues: (1) how to identify asymmetrical links from topological information, and (2) how to embed directed networks into the hyperbolic space and whether it is feasible to use empirical data to test the model. e results show that the directed links are hidden in the topological information in a nontrivial way. Based on empirical data, we mapped some real directed networks into hyperbolic space, including economic systems and biological ecosystems. We found that our method can reveal the topological features of directed networks, such as the degree of nodes, degree distribution, clustering coefficient, and core-periphery structure.
Furthermore, we analyzed the importance of nodes and evolutionary rules through visualization technology. e results show that the spatial positions of nodes in the hyperbolic space can be used to quantify the importance of nodes. e position change of nodes in the space is consistent with the evolution of state status in the economic system. More importantly, similar phenomena are observed in both socioeconomic systems and nervous systems; namely, there are significant discrepancies between spatial aggregation and the community structure of the topology. is is reflected in two ways: one is that the common information entropy between the community structure of hyperbolic space and the spatial distribution of nodes is particularly low, and the other is that the Pearson correlation coefficient between the node hyperbolic distance and the true distance is not significant. Long-range links weaken the spatial agglomeration of system functions, which means that the hyperbolic embedding reflects the clustering system functions. e effective distance of the network is the result of the trade-off between topological information and geometric features.
It needs to be reiterated that, unlike the binary network, the basis for the construction of the directed network is the trade-off between the four coordinates r ai , r bi , θ ai , and θ bi . Undirected network embedding is only based on the tradeoff between popularity and similarity, which may bias the prediction of the network links at some point. For example, consider the extreme case: if two nodes have large out-degrees and the in-degree is 0, the two nodes must be without links.
e chance of their being linked when undirected networks are embedded is large. Last but not least, asymmetric links increase the difficulty in embedding, which is reflected in the embedding accuracy to some extent. e frontiers of network geometry point that the community structure can be used as a coarse version of its embedding in a hidden space with hyperbolic geometry [44,45]. From the perspective of the mesoscale structure, the contribution of community structure should be incorporated into the embedded model in future studies, which may improve the embedding accuracy and reduce the algorithm complexity.
the scatter points are gathered near two lines, which can be spliced into a line after translation along the X-axis. However, the scattered points in Figure 6(b) are disordered, which shows that the estimation deviation of the node angle in set A is large. e reason for this result is more likely related to the process of node angle estimation and the challenge of asymmetry.
Data Availability e data used in this article are from public databases and can be downloaded from the relevant website. Code is available from the authors at bnuzonywu@163.com.

Conflicts of Interest
e authors declare that they have no conflicts of interest.