The Most Trustworthy Path Selection in Online Social Networks

In online social networks, it is crucial for a service consumer to find the most trustworthy path to a target service provider from numerous social trust paths between them. The selection of the most trustworthy path (namely, optimal social trust path (OSTP)) with multiple end-to-end quality of trust (QoT) constraints has been proved to be NP-Complete. Heuristic algorithms with polynomial and pseudo-polynomial-time complexities are often used to deal with this challenging problem. However, existing solutions cannot guarantee the efficiency of searching; that is, they can hardly avoid obtaining partial optimal solutions during searching process. Quantum annealing uses delocalization and tunneling to avoid falling into local minima without sacrificing execution time. It has been proved to be a promising way to many optimization problems in recently published literature. In this paper, for the first time, QA based OSTP algorithm (QA OSTP) is applied to the selection of the most trustworthy path. The experiment results show that QA based algorithm has better performance than its heuristic opponents.


Introduction
Online social networks (OSNs) [1,2] have gained much attention from the world recently.According to Nielson's report in 2009, more than two-thirds of the global online population visit and participate in online social networks [3].In such kind of networks, the network structure is made up of participants such as individuals or organizations and links between participants such as interactions, relationships, and connections.Since there may be over thousands of paths between a pair of interactive participants in large scale OSNs [4], a primary concern arise spontaneously, that is, for a source participant (e.g., a service consumer), how to exactly find the most trustworthy path to a target participant (e.g., a service provider) within a tolerable period of time?Unfortunately, the selection of the most trustworthy path, in other words, optimal social trust path (OSTP) selection with multiple end-toend quality of trust (QoT) constraints, has been proved to be NP-Complete [5].Evaluating the trustworthiness of the target participant along all these social trust paths requires large amount of computation time [6].Polynomial and pseudopolynomial-time heuristic algorithms are often used to deal with this problem.However, existing solutions cannot guarantee the quality of final configuration or the speed of searching.Thus, this field is attracting more and more attention recently with the rapid development of OSNs.

Related Work
2.1.Multiple QoS Constrained Path Selection.Korkmaz and Krunz [7] firstly proposed an approximation algorithm H MCOP to determine a feasible path that satisfies a set of constraints while maintaining high utilization of network resources.In their algorithm, both multiconstrained values and values of QoS parameters values are aggregated based on where  ≥ 1;   () is the aggregated value of the th QoS parameter of path ;   is the th QoS constraint value of the selected path between source node  and destination node .
H MCOP firstly uses Dijkstra's algorithm to find the path with the minimum   from  and , which intends to find whether there exists a feasible solution satisfying all end-toend QoS parameters in a subnetwork.During this process, the aggregated value of each QoS parameter for the identified path from  to  is computed and recorded at each intermediated node along the path.If there exists at least one feasible solution, then these aggregated values are used in another search from  to , which intends to identify a feasible path from  to  with the minimal cost.
Yu et al. [8] studied the problem of service selection with multiple QoS constraints and proposed an approximation algorithm, MCSP-K based on H MCOP, which checks the number of paths it currently keeps and removes the path with the maximum   if the path number reaches .In their algorithm, the paths with  minimum   will be kept at each intermediate node, which will ensure that MCSP-K will never prune out a feasible path if one exists.In the service candidate graph therein, there is a link between any two nodes in adjacent service sets.If these requirements cannot be satisfied in a network, MCSP K will search all the paths from source node to each intermediate that will lead the time complexity to be exponential.Thus, this algorithm does not fit for large scale online social networks.

Heuristic Algorithm for OSTP.
Until now, there are only a few works that are proposed to address the problem in online social networks where some significant influence factors including trust, recommendation roles, and social relationships are taken into account.
Liu et al. [5] developed a novel efficient heuristic algorithm for OSTP selection named MFPB HOSTP in online social networks based on Dijkstra's algorithm [9].In MFPB HOSTP, they first proposed the objective function given in where   ,   ,   are QoT parameters and represent trust information between participants, social relationships between adjacent participants, and recommendation roles of a participant, respectively.They adopted the Backward Search procedure to identify the path with the minimal  from  to  to investigate whether there exists a feasible solution where all end-to-end QoT constraints can be satisfied in the subnetwork and to record the aggregated QoT parameters of the path identified from  to each intermediate node.If a feasible solution exists, MFPB HOSTP then adopts the Forward Search procedure to find a near-optimal path from  to .This process adopts the information provided by Backward Search to identify whether there is another path that satisfies QoT constraints.In this process, MFPB HOSTP first searches the path with maximal F value from .Assume one of the adjacent nodes of ,  is selected based on Dijkstra's shortest path algorithm as the utility of the path from  to it is maximal.Let  ()  →  denote the backward local path from  to  identified in the Backward Search procedure.Then a foreseen path from  to  via  is formed.If this foreseen path is feasible, then MFPB HOSTP chooses the next node from  with the maximal F value that is calculated based on Dijkstra's shortest path algorithm.Otherwise, MFPB HOSTP does not search the path from  and the link  →  will be deleted from the subnetwork.Subsequently, MFPB HOSTP performs the Forward Search procedure to search the path from  in the subnetwork without the link  → .MFPB HOSTP is one of the most promising algorithms in solving the OSTP selection problem as it outperforms prior exiting algorithms in both efficiency and the quality of delivered solutions [5].

Quantum Annealing (QA).
In statistical mechanics, a physical process called simulated annealing (SA) is often performed in order to relax the system to the state with the minimum energy.In the basic form of SA, it first generates an initial solution as the current feasible solution using Metropolis algorithm [10,11].Then another solution is selected in the neighborhood of the current solution and replaces the current solution with the new one with the following transition probability given by Metropolis criterion: where (,   ) min {1,  (()−())/ } otherwise, where (, ) represent the probability with which configuration  is derived from  and () ⊆  is the neighborhood set of .With the temperature decreasing, only the better deterioration configuration can be accepted.Simulated annealing is a particularly promising minimization technique.For example, it has proved effective in finding the global minimum of multidimensional functions having large numbers of local minima [12].
It is worth mentioning that introducing thermal fluctuation is not the only way to perform annealing; QA depends on quantum fluctuation instead [13].A prominent advantage of quantum fluctuation over thermal fluctuation originates from the fundamental property of the quantum theory, namely, the possibility of tunneling through classically impenetrable potential barriers between energy valleys.Consequently, methods of quantum search, in principle, could be more efficient than the classical search methods [14][15][16][17].Furthermore, the effect of quantum tunneling is shown to be crucial for solving many computationally difficult problems, including the class of nondeterministic polynomial time problems.A practical implementation of QA will need to solve the timedependent Schrödinger equation in a very large and exponentially growing Hilbert space, which can only rely on a robust quantum computer.Recent research in this area has been carried out by Green's Function Monte Carlo (GFMC) simulations using quantum-classical mapping with the aid of a Suzuki-Trotter transformation [18,19], inspired from which we apply QA to the OSTP in OSNs.To the best of our knowledge, this is the first application of QA to the OSTP, and experiment results show that QA based searching algorithm has better performance than its heuristic opponent.
The structure of the present paper is as follows.In Section 3, we describe our proposed algorithm in detail.Specifically, in Section 3.1 Quantum Hamiltonian is given, which is then approximated by a classical one with the aid of a Suzuki-Trotter transformation.The QA based OSTP algorithm (QA OSTP) is presented and efficiency concerns are addressed in Sections 3.2 and 3.3, respectively.Section 4 contains experimental results and analysis, which show that QA OSTP outperforms MFPB HOSTP.Finally, we conclude this paper in Section 5 with a summary and research prospects.

Algorithm Description
In order to give prominence to the main problem of OSTP in online social networks, the Enron email dataset [20], a widely used system in the investigation of social networks [5,[21][22][23], is selected for our research, and the related parameter sets such as QoT set and set of weight of QoT parameters are the same as those in [5] for comparison.
In order to solve OSTP in large scale online social networks using quantum annealing, we first need to map the problem onto a highly constrained Ising model [24,25].Then searching the optimal social trust path is corresponding to finding the ground state with the lowest energy in an Ising model.Moreover, some assumptions should be given first as follows to eliminate some secondary factors that will increase unnecessary complexity and may influence the performance of the online social networks.
Assumption 1.Before searching the OSTP, three QoT parameters, that is, the values of trust, the social intimacy degree between participants, and the role impact factor of participants mentioned in [5], have been already obtained through mining techniques [22].Assumption 2. Since the transverse Ising spin glass (TISG) model is the simplest model in which quantum effects in a random system can and have been studied extensively and systematically [26], here we focus only on the TISG.Assumption 3. The OSNs studied here are symmetric.Symmetric means that each pair of adjacent participants, truster and trustee, can reverse their roles without changing their trust values.For example, in Figure 1, we set  AC =  CA ,  AC =  CA .

Problem Representation.
In a TISG model, phase space is spanned by all the set of spin variables {  }, each of which corresponds to a possible configuration.
(1) Quantum Hamiltonian.In QA, the total Hamiltonian for the OSTP can be expressed as a time-dependent Hamiltonian: where  pot represents the classical potential energy of a given configuration and  kin is a suitable kinetic energy operator providing the necessary quantum fluctuations to escape local minima.In QA, we seek to minimize  pot as side effect of minimizing  OSTP .A suitable configuration is reached if and only if  pot is zero.
In TISG model, the total Hamiltonian in ( 5) can be rewritten as where  TISG denotes potential energy of TISG model and  TF () is a fictitious kinetic energy introduced typically by the time-dependent transverse field.
A TISG model consists of a set of spins, each of which can only be in one of two states.Each of these spin variable usually takes on the value of either +1 or −1, also known as an up-spin and a down-spin, respectively.Formally, for online social networks with  participants each configuration of the system (a feasible social trust path) is associated to a  ×  matrix  with 0/1 entries in the following way: for each pair of the participants  and , if the directed social trust path (an ordered sequence of participants in online social networks) goes through the link between  and , then   = 1, or else   = 0.Here we renumber the independent variables as   ( = 1, . . ., ) and the other dependent variables can be expressed by   .Then, in quantum mechanics, the quantum Hamiltonian of the OSTP can be expressed as a -spins TISG model through the transformation   → (1 +    )/2, where    is the Pauli matrix of qubit .Moreover, the possible social trust paths can be represented by different quantum states of the  spins where  and  represent two randomly selected participants, that is, two qubits in TISG model, V  is the off-resonance term related to ,    and    describe the Pauli matrix of qubit  and , and   denotes the spin-spin coupling between qubit  and .
We use Figure 2 and Table 1 as a demonstration to describe this kind of correspondence.
In Figure 2, OSNs are composed by four participants, where all the parameters are the same as those defined in [5].There are two possible social trust paths between participant 1 and participant 4, that is, 1-3-4 and 1-2-4.The utility values of these two paths determined by the utility function in  [5] are 0.384 and 0.744, respectively.Thus even path 1-3-4 has larger trust value and , participant 1 will choose 1-2-4 as its preferred social trust path according to the weights specified by itself.Note that the realization of quantum annealing requires introducing an artificial and adjustable quantum kinetic operator that can provide the quantum fluctuations to escape the local minima.Moreover, the quantum annealing process is required to be slow enough to approximate the adiabatic evolution.Reasonably, the choice of  kin should encompass one important question that is determining which configurations are to become direct neighbors of a given configuration.Here we use an effective tactics named "plus-minus" [27] as disturbance mechanism to realize a neighborhood of a given configuration.Then in TISG model, to implement QA, a fictitious kinetic energy is introduced typically by the timedependent transverse field: where Γ() is the time-dependent power of the transverse field, each  ± ⟨,⟩ flips an Ising spin variable at position (, ) and at the symmetric position (, ), that is,  ± ⟨,⟩ =  ± ,  ± , .Then we have the following form of the total quantum Hamiltonian of OSTP: Initially the strength of the transverse field Γ() is chosen to be very large, and  OSTP is dominated by the third term of (9).Then Γ() is gradually and monotonically decreased toward zero, leaving eventually only the first two terms.Note that every state of TISG model can be described as state vector |()⟩, and it will evolve with time and should follow the RT Schrödinger equation, With Γ() decreasing, accordingly the state vector |()⟩ is expected to evolve from the trivial initial ground state of transverse-field term (6) to the nontrivial ground state of (7), which is the solution of the OSTP.
Then an important issue arises, that is, how slowly we should decrease Γ() to keep the state vector arbitrarily close to the instantaneous ground state of  OSTP .As mentioned in Section 2.3, we will not attempt an actual Schrödinger annealing evolution of the quantum Hamiltonian due to the large Hilbert space.On the contrary, we address the quantum problem by GFMC-QA [25], where annealing will take place in the fictitious time represented by the number of Monte Carlo steps.However, in order to figure out OSTP by GFMC-QA, Suzuki-Trotter transformation should be performed in advance, which requires calculation of the matrix elements of an exponential operator between arbitrary configurations |⟩ and ⟨  | of the system, a complicated issue for  Kin .Moreover, since the energy gap between the ground state and the first excited state is large at the beginning and decreases with the annealing time, hyperbolic interpolation makes the annealing process more efficient and smoother than the linear one does [18]; we make the strength of the transverse field as following hyperbolic interpolation: where  is the total annealing time and ,  are two control operators.This form is trivially Trotter-discretized [18,19], since the spin-flip term acts independently on single spins at each time slot.
(2) Suzuki-Trotter Transformation.In GFMC-QA, a quantum Hamiltonian is approximated by a classical one with the aid of a Suzuki-Trotter transformation.This is possible because of an analogy with a standard TISG model in a transverse field [13].The transformation maps the quantum Hamiltonian to an effective classical Hamiltonian  similar to the one mentioned in [24], and then (9) can be rewritten as where  can be viewed as consisting of  replicas {Υ , ,  = 1, . . ., } of the classical potential energy of a given configuration  pot ({Υ  }), with an interaction of a combined kinetic energy between them, and Υ , denotes the th spin of the th replica.The term  Γ is the coupling among the replicas, which can be written as where  is the temperature, at which each replica is simulated.

Proposed QA OSTP.
In this part, we will describe the components of our algorithm in detail.Notations that are used in QA OSTP are shown in Table 2. Flow chart of QA OSTP is given in Figure 3.The first diamond (outermost loop) is controlled by .Here a linear annealing schedule consisting of the initial temperature  0 and Max steps is selected.Each Monte Carlo step for QA OSTP consists of a loop starting from the second diamond where  is a tunable multiplier;  *  moves are conducted at each step after which the control parameter is decreased.QA OSTP keeps making the next Monte Carlo Step each time until the termination condition is satisfied.The replicas are always connected to each other in numerical order in the same way throughout the search for the purpose of spin products.The random disturbing means only changing the order in which replicas are selected for search.We find this as an efficient scheme that can promote diversity in the population of configurations.

Theoretical Analysis.
The convergence conditions for the implementation of QA OSTP with the quantum Monte Carlo evolution are investigated in this part.In our method, the -dimensional TISG is mapped to a ( + 1)-dimensional classical Ising system.By using the Suzuki-Trotter transformation, at temperature , when the length of the extra dimension  → ∞, the system partition function can be defined as where  is Boltzmann's constant,  ()  represents a classical Ising spin at site  on the kth Trotter slice, and  is the space of discrete configuration.
A Monte Carlo step is characterized by the transition probability from configuration   to configuration    at time step : where   is the present configuration,    is the next candidate configuration which is generated with defined probability, both of which belong to the space of discrete configuration , and (   ,   ) is the generation probability with which    can be generated from   .The acceptance probability Λ(   ,   ; ) of QA OSTP can be written as where  0 (  ) means the cost function whose global minimum is the optimal social trust path of OSTP problem in OSNs,  1 (  ) represents the kinetic energy, Θ(  ; ) is the stationary distribution of the homogeneous Markov chain defined by matrix I at a given  1 (), and Φ(⋅) is the monotone increasing acceptance function satisfying 0 ≤ Φ(⋅) ≤ 1.
To derive the convergence conditions for the implementation of QA OSTP with the quantum Monte Carlo evolution, we should prove that inhomogeneous Markov chain associated with QA OSTP is strongly ergodic under appropriate conditions [28].Our main results are summarized in the following theorems.
Theorem 4. For a causal system, the transition matrix I has the following lower bound: where  is the minimum nonvanishing value of (   ,   ), and  0 and  1 are the maximum changes in a single step in  0 (  ) and  1 (  ), respectively.
Proof.Directly following the definition of the transition probability and the property of the acceptance function in (15), (16), and ( 17), for both positive  0 and  1 , if Θ(   ; )/Θ (  ; ) < 1, we get else, we have Theorem 4 is thereby proved.
Theorem 5.The inhomogeneous Markov chain generated by (16) and ( 17) is strongly ergodic and converges to the equilibrium state, which is corresponding to the term exp(− 0 (  )/  0 ) in (17).
Proof.In order to prove strong ergodicity, we refer to the conditions for strong ergodicity [28].If there exists the transition matrix N on  such that N(   ,   ) = N(   ,    ) for any   ,    ,    ∈ , then for a causal system, we have lim Since the Markov chain is proved to be weakly ergodic in [28], we consider for any  > 0, (, ) = I ,  0 , where  0 belongs to the set of probability distributions on   ∈ .For Therefore, the inhomogeneous Markov chain generated by ( 16) and ( 17) is strongly ergodic and thereby is proved.

Experiment Settings.
In our experiments, if no otherwise specified, all the related parameters are set following reference [5].In order to evaluate our proposed algorithm, we compare QA OSTP with MFPB HOSTP in terms of two key factors, that is, execution time and the utility of the selected social trust path.
For QA OSTP, we implement a similar GFMC that was used in [25] at a fixed low temperature  (we used  = 10/3).The quantum model is mapped onto a classical model with an extra imaginary-time dimension, consisting of  ferromagnetically coupled replicas of the original spin problem, at temperature  [25] (we used  = 30).Since QA requires initial configurations equilibrated at temperature , an obvious choice is to take  = 100 [25].Finally, the transverse field Γ is annealed hyperbolically in a MC time  from an initial value Γ 0 = 300 to a final value of zero.In QA, we used exclusively "plus-minus" tactics, with a static neighborhood pruning [29], which restricts the attempted neighborhood realization by allowing only a fixed number  (we used  = 20) of nearest neighbors of participant  to be candidates for   .Our MC step consisted of  *  attempted operations of "plusminus" tactics (for QA, in each of the  replicas).In QA, we averaged the best social trust path utility found over up to 100 independent searches.In our experiments, the three QoT parameters are randomly generated.The end-to-end QoT constraints specified by a source participant are set as   ≥ 0.05,   ≥ 0.001, and   ≥ 0.3, respectively.We first randomly select 80 pairs of source and target participants from the Enron email dataset with 87,474 nodes and 30,0511 links.Moreover, following the small world characteristic we set the maximal length of a social trust path 6 hops.Then we number the different network scales from 1 to 25, with number of nodes varying from 50 to 400 and number of links varying from 63 to 2356, respectively.
Both algorithms were implemented in Matlab 7.0 and run on a PC with a 3 GHz Intel processor and 3 GB of RAM with Windows 7. From the figure, we can observe that our QA OSTP can always find utilities that are not worse than those of MFPB HOSTP (e.g., case S1 and S2).This is because in QA OSTP, quantum fluctuation is adopted to avoid the local minima and quantum mechanics works with wave functions that can sample different regions of phase space equally well.While for MFPB HOSTP, although subnetworks are extracted through exhaustive searching before the algorithm execution and Backward Search scheme is used to estimate whether there exists a feasible solution in a subnetwork, due to the intrinsic characteristic of heuristic algorithms, it does not guarantee that the best social trust path will be found.Moreover with the network scale growing larger, when the social trust path with the maximal utility is not a feasible solution, the heuristic search can hardly find a near optimal solution and usually returns an infeasible one even when a feasible solution exists (e.g., case S3).Thus, in any case, QA OSTP shows better performance.In particular, we find that the mean value of utilities of QA OSTP is 42.85% more than that of MFPB HOSTP in Figure 4(a), 58.50% more in Figure 4(b), 55.84% more in Figure 4(c), and 63.12% more in Figure 4(d).

Performance Analysis
Figure 5 shows the average execution time of algorithms with different weights of QoT parameters and different network scales.Note that, for MFPB HOSTP, execution time should include the exhaustive searching time for extracting subnetworks, which is not taken into consideration in reference [5].From Figure 5, we can see when the network scale is not large, both MFPB HOSTP and QA OSTP perform well and the difference is trivial because the searching space is relatively small.However, with the network scale expanding, we can observe that QA OSTP can outperform MFPB HOSTP in execution time.This is an interesting result since from conventional point, annealing process may be a little time-consuming.The reason we think is that the quantum Hamiltonian in QA OSTP is approximated with the aid of a Suzuki-Trotter transformation in GFMC-QA and quantum tunneling can avoid some unnecessary searching in MFPB HOSTP, which can accelerate the annealing process.From Figure 5, we can see that the average execution time of our proposed algorithm is only 33.47% of that of MFPB HOSTP in Figure 5(a), 28.36% in Figure 5(b), 31.99% in Figure 5(c), and 26.78% in Figure 5(d).
Based on the above experiments conducted with different scales and parameters, we can observe that QA OSTP is a promising algorithm and it performs better than its heuristic opponent in both the quality of the selected social trust path and the execution time.

Conclusion
A novel quantum annealing based OSTP algorithm, that is, QA OSTP for selection of the most trustworthy path to service provider in online social networks, was proposed.To the best our knowledge, this is the first application of quantum annealing to the challenging NP-Complete OSTP problem in online social networks.Since quantum mechanics work with wave functions that can sample different regions of phase space equally well and quantum systems can tunnel through classically impenetrable potential barriers between energy valleys, a process that might prove more effective than waiting for them to be overcome thermally as in SA, QA OSTP is able to outperform its heuristic opponents and even find configuration of comparable quality to the best algorithm.Thus, QA OSTP is proven to be a very promising tool for solving the OSTP in online social networks.
As for the future work, understanding how quantum mechanics can quantitatively improve the quality of solution of OSTP is still an important open issue.Moreover, since GFMC simulations using quantum-classical mapping with the aid of a Suzuki-Trotter transformation only simulate the equilibrium behavior at finite temperature, we plan to devise another effective and alternative scheme to solve the infinite time Schrödinger equation with stochastic processes.

4 Figure 4 :
Figure 4: Comparison of path utilities of networks.

4 Figure 5 :
Figure 5: Comparison of execution time of algorithms.

Table 1 :
Possible social trust path and corresponding state vector in Figure2.

Table 2 :
Notations used in QA OSTP.