Harmonic Centrality is a variant of Closeness Centrality. The average shortest distance measurement proposed by harmonic centrality is compatible with infinite values which would occur in a disconnected graph. Harmonic centrality was first proposed by M. Marchiori and V. Latora in 2000, and then by A. Dekker and Y. Rochat in 2005 and 2009:
Harmonic centrality ranges from 0 to 1; higher scores indicate that a node is closer to other nodes in the graph.
The shortest distance between two nodes is defined as the number of edges in the shortest path connecting them. Please refer to Closeness Centrality for more details.
The harmonic mean is the reciprocal of the arithmetic mean of the reciprocals of the variables. The formula for calculating the arithmetic mean A and the harmonic mean H is as follows:

A classic application of harmonic mean is to calculate the average speed when traveling back and forth at different speeds. Suppose there is a round trip, the forward and backward speeds are 30 km/h and 10 km/h respectively. What is the average speed for the entire trip?
The arithmetic mean A = (30+10)/2 = 20 km/h is not appropriate in this case. Since the backward journey takes three times as long as the forward, during most time of the entire trip the speed stays at 10 km/h, so we expect the average speed to be closer to 10 km/h.
Assuming the one-way distance is 1, the average speed that takes travel time into consideration is 2/(1/30+1/10) = 15 km/h. This value, the harmonic mean, is adjusted by the time spent during each journey.
Harmonic centrality score of a node defined by this algorithm is the inverse of the harmonic mean of the shortest distances from the node to all other nodes. The formula is:

where x is the target node, y is any node in the graph other than x, k-1 is the number of y, d(x,y) is the shortest distance between x and y, d(x,y) = +∞ when x and y are not reachable to each other, in this case 1/d(x,y) = 0.

The harmonic centrality of node a in the above graph is (1 + 1/2 + 1/+∞ + 1/+∞) / 4 = 0.375, and the harmonic centrality of node d is (1/+∞ + 1/+∞ + 1/+∞ + 1) / 4 = 0.25.
NOTEHarmonic Centrality algorithm consumes considerable computing resources. For a graph with V nodes, it is recommended to perform (uniform) sampling when V > 10,000, and the suggested number of samples is the base-10 logarithm of the number of nodes (
log(V)).
For each execution of the algorithm, sampling is performed only once, centrality score of each node is computed based on the shortest distance between the node and all sample nodes.

Run the following statements on an empty graph to define its structure and insert data:
ALTER GRAPH CURRENT_GRAPH ADD NODE { user () }; ALTER GRAPH CURRENT_GRAPH ADD EDGE { vote ()-[{score uint32}]->() }; INSERT (A:user {_id: "A"}), (B:user {_id: "B"}), (C:user {_id: "C"}), (D:user {_id: "D"}), (E:user {_id: "E"}), (F:user {_id: "F"}), (G:user {_id: "G"}), (H:user {_id: "H"}), (A)-[:vote {score: 2}]->(B), (A)-[:vote {score: 3}]->(E), (B)-[:vote {score: 4}]->(B), (B)-[:vote {score: 2}]->(C), (C)-[:vote {score: 3}]->(A), (D)-[:vote {score: 1}]->(A), (F)-[:vote {score: 1}]->(G);
To load the entire graph to the HDC server hdc-server-1 as my_hdc_graph:
CREATE HDC GRAPH my_hdc_graph ON "hdc-server-1" OPTIONS { nodes: {"*": ["*"]}, edges: {"*": ["*"]}, direction: "undirected", load_id: true, update: "static" }
Algorithm name: harmonic_centrality
Name | Type | Spec | Default | Optional | Description |
|---|---|---|---|---|---|
ids | []_id | / | / | Yes | Specifies nodes for computation by their _id. If unset, computation includes all nodes. |
uuids | []_uuid | / | / | Yes | Specifies nodes for computation by their _uuid. If unset, computation includes all nodes. |
direction | String | in, out | / | Yes | Specifies that all edges in the shortest paths must be either incoming (in) or outgoing (out). |
edge_schema_property | []"<@schema.?><property>" | / | / | Yes | Specifies numeric edge properties used as weights by summing their values. Only properties of numeric type are considered, and edges without these properties are ignored. |
impl_type | String | dijkstra, delta_stepping, spfa, beta | beta | Yes | Specifies the algorithm used to compute weighted shortest paths: Dijkstra, Delta-Stepping, SPFA or the default (beta) Ultipa algorithm. Valid only when edge_schema_property is specified. |
sample_size | Integer | -1, -2, [1, |V|] | -2 | Yes | Specifies the sampling strategy for computation:
|
return_id_uuid | String | uuid, id, both | uuid | Yes | Includes _uuid, _id, or both in the results to represent nodes. |
limit | Integer | ≥-1 | -1 | Yes | Limits the number of results returned. Set to -1 to include all results. |
order | String | asc, desc | / | Yes | Sorts the results by harmonic_centrality. |
CALL algo.harmonic_centrality.write("my_hdc_graph", { return_id_uuid: "id", order: "desc" }, { file: { filename: "harmonic" } })
Result:
File: harmonic_id,harmonic_centrality A,0.571429 B,0.428571 C,0.428571 D,0.357143 E,0.357143 F,0.142857 G,0.142857 H,0
Writes the harmonic_centrality values from the results to the specified node property. The property type is float.
CALL algo.harmonic_centrality.write("my_hdc_graph", {}, { db: { property: "hc" } })
CALL algo.harmonic_centrality.run("my_hdc_graph", { return_id_uuid: "id", ids: ["A", "B"], edge_schema_property: "score" }) YIELD hc RETURN hc
Result:
| _id | harmonic_centrality |
|---|---|
| A | 0.309523 |
| B | 0.219048 |
CALL algo.harmonic_centrality.stream("my_hdc_graph", { direction: "in", return_id_uuid: "id" }) YIELD hc FILTER hc.harmonic_centrality = 0 RETURN hc
Result:
| _id | harmonic_centrality |
|---|---|
| D | 0 |
| F | 0 |
| H | 0 |