Abstract

We investigate the multiple attribute group material selection problems in which the attribute values take the form of interval 2-tuple linguistic information. Firstly, some operational laws and possibility degree of interval 2-tuple linguistic variables are introduced. Then, we develop some interval 2-tuple linguistic aggregation operators called interval 2-tuple hybrid harmonic mean (ITHHM) operator, induced interval 2-tuple ordered weighted harmonic mean (I-ITOWHM) operator, and induced interval 2-tuple hybrid harmonic mean (I-ITHHM) operator and study some desirable properties of the I-ITOWHM operator. In particular, all these operators can be reduced to aggregate 2-tuple linguistic variables. Based on the I-ITHHM and the ITWHM (interval 2-tuple weighted harmonic mean) operators, an approach to multiple attribute group decision-making with interval 2-tuple linguistic information is proposed. Finally, a practical application to material selection problem is given to verify the developed approach and to demonstrate its practicality and effectiveness.

1. Introduction

The 2-tuple linguistic representation model, characterized by a linguistic term and a numeric value, was developed by Herrera and Martínez [1] based on the concept of symbolic translation. It has exact characteristic in linguistic information processing and can effectively avoid information distortion and loss which occur formerly in the linguistic computing process. In many real decision-making cases, experts tend to give opinions according to their experience and knowledge, and it is more suitable to provide assessments by means of linguistic terms rather than numerical ones [2, 3], due to the complexity of the objects and the vagueness of human thinking. Therefore, many researchers have investigated linguistic multiple attribute group decision-making (MAGDM) problems and proposed lots of methods to deal with linguistic evaluation information. These linguistic computational models can be mainly classified into three types [3, 4]: the method based on membership functions, the method based on linguistic symbols, and the method based on linguistic 2-tuples. However, one important weakness of the former two linguistic computational models is that they performed the retranslation step as an approximation process to express the results in the original expression domain (initial term set) thus provoking a lack of accuracy [5]. In contrast, the 2-tuple linguistic method is a symbolic model that extends the use of indexes modifying the fuzzy linguistic approach representation by adding a parameter to the basic linguistic representation in order to improve the accuracy of the linguistic computations and the interpretability of the results. Due to its distinguished power and efficiency in dealing with linguistic assessments of decision-makers, the 2-tuple linguistic model has been extensively applied to decision-making processes since its appearance [610].

In recent years, much progress has been made in research relating to 2-tuple aggregation operators since information aggregation plays a significant part in the MAGDM process. For example, Herrera and Martínez [1] defined the 2-tuple arithmetic mean operator, the 2-tuple weighted averaging operator, and the 2-tuple ordered weighted averaging operator. Wei [11] proposed the extended 2-tuple weighted geometric (ET-WG) and the extended 2-tuple ordered weighted geometric (ET-OWG) operators and Wei [12] developed some generalized aggregating operators with 2-tuple linguistic information, such as the generalized 2-tuple weighted average (G-2TWA), the generalized 2-tuple ordered weighted average (G-2TOWA), and the induced generalized 2-tuple ordered weighted average (IG-2TOWA) operators. Yang and Chen [13] considered the linguistic group decision-making problem with interdependent attributes and introduced some 2-tuple linguistic aggregation operators by using the Choquet integral: the 2-tuple correlated averaging (TCA) operator, the 2-tuple correlated geometric (TCG) operator, and the generalized 2-tuple correlated averaging (GTCA) operator. Xu and Wang [14] developed the 2-tuple linguistic power average (2TLPA), the 2-tuple linguistic weighted power average (2TLWPA), and the 2-tuple linguistic power ordered weighted average (2TLPOWA) operators and studied their desired properties, such as commutativity, idempotency, and boundary. Wan [15] developed some hybrid arithmetic aggregation operators with 2-tuple linguistic information, involving the 2-tuple hybrid weighted arithmetic average (THWA) operator, the 2-tuple hybrid linguistic weighted arithmetic average (T-HLWA) operator, and the extended 2-tuple hybrid linguistic weighted arithmetic average (ET-HLWA) operator. Recently, Sonia et al. [16] introduced some linguistic aggregation operators with conservation of interaction between criteria, which include the 2-tuple Choquet integral averaging (TCIA) operator, the 2-tuple ordered Choquet integral averaging (TOCIA) operator, and the combined 2-tuple Choquet integral averaging operator. Considering the interactive phenomenon among experts (or attributes) in MAGDM, Lin et al. [17] proposed a generalized interval 2-tuple linguistic Shapley chi-square averaging operator for facility location selection. Li and Liu [18] proposed some new aggregation operators of 2-tuple linguistic information based on Heronian mean, and Ju et al. [19] gave some new Shapley 2-tuple linguistic Choquet aggregation operators for MAGDM. The linguistic proportional 2-tuple power average operator was introduced by Jiang et al. [20] to aggregate linguistic values of unbalanced linguistic term sets considering the relationship among the aggregated values. Besides, the 2-tuple linguistic extended Bonferroni mean (EBM) aggregation operators and the 2-tuple linguistic partition Bonferroni mean (PBM) aggregation operators were developed in [21] and [22], respectively, for dealing with MAGDM problems.

In some situations, however, the input arguments take the form of interval 2-tuple linguistic values because of time pressure, lack of knowledge or data, and decision-makers’ limited attention and information processing capabilities [2325]. Furthermore, all decision-makers cannot easily express their evaluations on the established criteria with one given linguistic term set. That is, some decision-makers feel that the cardinality of the linguistic term set is too small to fully express their judgments on the alternatives, while others think it is so big that the evaluations on the alternatives are out of their ability. Therefore, Zhang [26] introduced the interval 2-tuple linguistic representation model and developed some interval 2-tuple linguistic aggregation operators, such as the interval 2-tuple weighted average (ITWA) and the interval 2-tuple ordered weighted average (ITOWA) operators. Zhang [27] further developed some new interval 2-tuple linguistic aggregation operators, including the interval 2-tuple weighted geometric (ITWG), the interval 2-tuple ordered weighted geometric (ITOWG), the generalized interval 2-tuple weighted average (GITWA), and the generalized interval 2-tuple ordered weighted average (GITOWA) operators. Liu et al. [28] put forward some Bonferroni mean operators under the interval-valued 2-tuple linguistic context and investigated their properties and special cases. Wu et al. [29] presented the 2-tuple linguistic generalized power average (2TLGPA), the 2-tuple linguistic generalized power ordered weighted average (2TLGPOWA), and the interval 2-tuple linguistic generalized power average (I2TLGPA) operators for MAGDM. Motivated by the idea of harmonic mean operators [30, 31], in this paper, we develop some interval 2-tuple linguistic harmonic mean operators, such as the interval 2-tuple hybrid harmonic mean (ITHHM) operator, the induced interval 2-tuple ordered weighted harmonic mean (I-ITOWHM) operator, and the induced interval 2-tuple hybrid harmonic mean (I-ITHHM) operator. And a new approach based on the I-ITHHM and the interval 2-tuple weighted harmonic mean (ITWHM) operators for MAGDM with interval 2-tuple linguistic information is presented.

The remainder of this paper is set out as follows. In Section 2, we introduce some basic concepts and operation laws of interval 2-tuple linguistic variables. In Section 3, we develop some new aggregation operators called the ITHHM operator, the I-ITOWHM operator, and the I-ITHHM operator and investigate some desirable properties of the I-ITOWHM operator. In Section 4, we develop an approach based on the I-ITHHM and the ITWHM operators to interval 2-tuple linguistic MAGDM. A material selection example is given in Section 5 to verify the developed approach and to demonstrate its feasibility and practicality. Finally, conclusions and future directions are provided in Section 6.

2. Preliminaries

2.1. 2-Tuple Linguistic Variables

A linguistic variable is a variable whose values are expressed in linguistic terms. In other words, it is a variable whose values are not numbers but words or sentences in a natural or artificial language. The concept of linguistic variables is very useful in dealing with situations which are too complex or too ill defined to be reasonably described by traditional quantitative expressions [32]. For identifying the diversity of each evaluation item and facilitating the computation, linguistic terms often possess some characteristics such as finite set, odd cardinality, semantic symmetric, ordinal level, and compensative operation [33]. Let be a linguistic term set with odd cardinality, where represents a possible value for a linguistic variable. It is required that the linguistic term set should satisfy the following characteristics [1]:(1)Negation operator: such that .(2)The set which is ordered: , if .(3)Max operator: , if .(4)Min operator: , if .

The 2-tuple linguistic representation model was firstly presented by Herrera and Martínez [1] based on the concept of symbolic translation. It is used to represent the linguistic information by means of a linguistic 2-tuple, , where is a linguistic term from the predefined linguistic term set and is a numerical value representing the symbolic translation. In the classical 2-tuple linguistic approach, the range of value is relevant to the granularity of the linguistic term sets. Here, is the result of an aggregation of the indices of a set of labels assessed in the linguistic term set . To overcome this restriction, Tai and Chen [33] further proposed a generalized 2-tuple linguistic model and translation functions.

Definition 1. Let be a linguistic term set and let be a value representing the result of a symbolic aggregation operation. Then, the generalized translation function used to obtain the 2-tuple equivalent to can be defined as follows [33]: where is the usual rounding operation. The interval of is determined by the number of linguistic terms in . For example, if contains seven linguistic terms, then and .

Definition 2. Let be a linguistic term set and let be a 2-tuple. There exists a function , which is able to convert the 2-tuple into its equivalent numerical value . The reverse function is defined as follows [33]:It is obvious that the conversion of a linguistic term into a linguistic 2-tuple consists of adding a value 0 as symbolic translation [1]:

2.2. Interval 2-Tuple Linguistic Variables

The interval 2-tuple linguistic representation model was put forward by Zhang [26] as a generalization of the 2-tuple linguistic variable.

Definition 3. Let be a linguistic term set. An interval 2-tuple linguistic variable is composed of two 2-tuples, denoted by , where , and and represent the linguistic label of the linguistic term set and symbolic translation, respectively. The interval 2-tuple that expresses the equivalent information to an interval value is derived by the following function [26, 27]:On the contrary, there is always a function such that an interval 2-tuple can be converted into an interval value as follows: In particular, if and , then the interval 2-tuple linguistic variable reduces to a 2-tuple linguistic variable.

Definition 4. Consider any three interval 2-tuples , , and and let ; then, their operations can be defined as follows [34]: Motivated by the formulas proposed by Xu [35, 36], the comparison of linguistic information represented by interval 2-tuples is implemented on the basis of the possibility degree of the interval 2-tuple linguistic variables. It can be defined as follows.

Definition 5. Let and be two interval 2-tuples, and let and let ; then, the possibility degree of is defined asSimilarly, the possibility degree of is defined as

From Definition  5, we can easily get the following results:(1);(2) . In particular,

3. Interval 2-Tuple Linguistic Harmonic Mean Operators

3.1. Interval 2-Tuple Hybrid Harmonic Mean Operator

Definition 6. Let WHM: , if where is a collection of positive real numbers, is the weight vector of , with and , and is the set of all positive real numbers; then, WHM is called the weighted harmonic mean operator [37].

Zhang [27] extended the WHM operator to accommodate the situations where the input arguments are interval 2-tuple linguistic variables. For convenience, we let be the set of all interval 2-tuple linguistic variables.

Definition 7. Let be a set of interval 2-tuples, and let ITWHM: , ifwhere is the weight vector of , with and ; then, ITWHM is called the interval 2-tuple weighted harmonic mean (ITWHM) operator.

Based on the OWA and the ITWHM operators, Zhang [27] developed the interval 2-tuple ordered weighted harmonic mean (ITOWHM) operator as follows.

Definition 8. Let be a set of interval 2-tuples. An ITOWHM operator of dimension is a mapping ITOWHM: , which has an associated weight vector with and , such that where is the th largest of .

Obviously, the fundamental characteristic of the ITWHM operator is that it considers the importance of each given interval 2-tuple linguistic variable, whereas the fundamental characteristic of the ITOWHM operator is the reordering step, and it weights all the ordered positions of interval 2-tuple linguistic variables instead of weighting the given interval 2-tuples themselves. In the following, by combining the advantages of the ITWHM and ITOWHM operators, we develop an interval 2-tuple hybrid harmonic mean (ITHHM) operator that weights both the given interval 2-tuple linguistic variables and their ordered positions.

Definition 9. Let be a set of interval 2-tuples. An ITHHM operator of dimension is a mapping ITHHM: , which has an associated weight vector with and , such thatwhere is the th largest of the weighted interval 2-tuples , is the weight vector of , with and , and is the balancing coefficient.
In particular, if , then the ITHHM operator is reduced to the ITOWHM operator; if , then the ITHHM operator is reduced to the ITWHM operator. Thus, the ITHHM operator generalizes both the ITWHM and ITOWHM operators and reflects the importance degrees of both the given argument and its ordered position. Moreover, if the interval 2-tuples are degenerated to the 2-tuples , then the ITHHM operator is reduced to the 2-tuple linguistic hybrid harmonic (2TLHH) operator [38].
To rank these interval 2-tuple linguistic arguments , we first compare each argument with all arguments by using (10) and let . Then, we can construct a complementary matrix , where , , and , Summing all elements in each line of the matrix , we have . Then, we can rank the arguments in descending order in accordance with the values of .

Example 10. Assuming is a linguistic term set, , , , and and let be the weight vector of . Then, we get the weighted interval 2-tuples:To rank these arguments, we first compare each argument with all arguments by using (10) and then construct a complementary matrixSumming all elements in each line of the matrix , we haveThen, we rank the arguments in descending order in accordance with the values of :Suppose that the weight vector of the ITHHM operator is ; then by (15), we get

3.2. Induced Interval 2-Tuple Ordered Weighted Harmonic Mean Operator

Definition 11. An induced ordered weighted harmonic mean (IOWHM) operator [39] is defined as follows: where is a weight vector such that and , is the value of the IOWA pair having the th largest , and in is referred to as the order inducing variable.
In the following, we will develop an induced interval 2-tuple ordered weighted harmonic mean (I-ITOWHM) operator.

Definition 12. Let be a set of interval 2-tuples and let be an associated weight vector, with . The induced interval 2-tuple ordered weighted harmonic mean (I-ITOWHM) operator is defined as where is the value of the IOWA pair having the th largest and in is referred to as the order inducing variable.
If there is a tie between and with respect to the order inducing variables such that , in this case, we replace the argument component of each of and by their average in the process of aggregation. If items are tied, then we replace these by replicas of their average. The weight vector can be determined by using some weight determining methods like the normal distribution based method [40].
In particular, if , then the I-ITOWHM operator is reduced to the interval 2-tuple harmonic mean (ITHM) operator:If , for all , where is the ordered position of , the I-ITOWHM operator is reduced to the ITWHM operator; if , for all , then the I-ITOWHM operator is reduced to the ITOWHM operator. Moreover, if the interval 2-tuples are degenerated to the 2-tuples , then the I-ITOWHM operator is reduced to the induced 2-tuple ordered weighted harmonic mean (I-TOWHM) operator:where is the value of the IOWA pair having the th largest , and in is referred to as the order inducing variable. If there is a tie between and with respect to the order inducing variables such that , then we replace the argument component of each of and by their average in the process of aggregation. If items are tied, then we replace these by replicas of their average.
The I-ITOWHM operator has the following properties similar to those of the IOWHM [39] and the IOWA operators [41].

Theorem 13 (commutativity). Considerwhere is any permutation of .

Proof. LetSince is any permutation of , we have ; then,

Theorem 14 (idempotency). If for all , then

Proof. Since for all , we have

Theorem 15 (monotonicity). If for all , then

Proof. LetSince for all , it follows that ; then,

Example 16. Assuming is a linguistic term set, , , , and are four IOWA pairs. Then, we desire to aggregate using the weight vector . Performing the ordering of the IOWA pairs with respect to the first component, we getThis ordering includes the ordered interval 2-tuple arguments:Thus, we get an aggregated value

3.3. Induced Interval 2-Tuple Hybrid Harmonic Mean Operator

Inspired by the induced hybrid averaging (IHA) operator [42], in the following, we further generalize the ITHHM operator by using order inducing variables, obtaining the induced interval 2-tuple hybrid harmonic mean (I-ITHHM) operator.

Definition 17. Let be a set of interval 2-tuples. An I-ITHHM operator of dimension is a mapping I-ITHHM: , which has an associated weight vector with and , such thatwhere is the weighted interval 2-tuple of the IOWA pair having the th largest is the order inducing variable, is the weight vector of , with and , and is the balancing coefficient.
In particular, if , then the I-ITHHM operator is reduced to the I-ITOWHM operator; if , then the I-ITHHM operator is reduced to the ITWHM operator. If , for all , then the I-ITHHM operator is reduced to the ITHHM operator. Moreover, if the interval 2-tuples are degenerated to the 2-tuples , then the I-ITHHM operator is reduced to the induced 2-tuple hybrid harmonic mean (I-THHM) operator:where is the weighted 2-tuple of the IOWA pair having the th largest is the order inducing variable, is the weight vector of , with and , and is the balancing coefficient.

Example 18. Assuming is a linguistic term set, , , , and are four IOWA pairs. Let be the weight vector of . Performing the ordering of the IOWA pairs with respect to the first component, we getThen, we get the weighted interval 2-tuplesSuppose that the weight vector of the I-ITHHM operator is ; then, by (36), we get

4. An Approach to MAGDM with Interval 2-Tuple Linguistic Information

In this section, we will develop an approach based on the proposed interval 2-tuple linguistic harmonic mean operators for the MAGDM problems in which both the attribute weights and the expert weights take the form of real numbers and the attribute values take the form of interval 2-tuple linguistic variables.

Suppose that a MAGDM problem has decision-makers , alternatives , and decision attributes . Each decision-maker is given a weight satisfying to reflect his/her relative importance in the group decision-making process. Let be the linguistic decision matrix of the th decision-maker, where is the linguistic information provided by on the assessment of with respect to. Let be the weight vector of attributes, where , and . In addition, decision-makers may use different linguistic term sets to express their preference values.

In the following, we apply the I-ITHHM and the ITWHM operators to multiple attribute group decision-making with interval 2-tuple linguistic information.

Step 1. Convert the linguistic decision matrix into interval 2-tuple linguistic decision matrix , where and .
Suppose that provides its assessments in a set of five linguistic terms and the linguistic term set is denoted as The linguistic information provided in the decision matrix can be converted into its corresponding interval 2-tuple linguistic assessments according to the following ways:(i)A certain grade such as Poor, which can be written as ;(ii)an interval such as Poor-Medium, which means that the assessment of an alternative with respect to the attribute under consideration is between Poor and Medium which can be written as .

Step 2. Utilize the I-ITHHM operator which has an associated weight vector with and to aggregate all the decision matrices into a collective decision matrix , where is the weighted interval 2-tuple of the IOWA pair having the th largest .

Step 3. Utilize the decision information given in matrix and the ITWHM operatorto derive the collective overall preference value of the alternative .

Step 4. To rank these collective overall preference values , we first compare each with all by using (10). Let , and then we develop a complementary matrix as , where , , and ,
Summing all elements in each line of the matrix , we have , Then, we rank in descending order in accordance with the values of .

Step 5. Rank all the alternatives and select the best one(s) in accordance with the collective overall preference values .

Step 6. End.

5. An Illustrative Example

In what follows, an illustrative example adapted from [43] is presented to illustrate the application of the developed method for material selection problem. This example problem is related with selection of a suitable work material for a product which needs to be designed for operating in a high-temperature environment. After preliminary screening, five alternative materials, , have been designated for further evaluation. An expert committee of four decision-makers, that is, DM1, DM2, DM3, and DM4, has been created in order to evaluate and select the most appropriate material for the application. The attributes which have been considered for the analysis are: tensile strength;: Young’s modulus;: density;: corrosion resistance.The weight vector of these attributes is . The four decision-makers, whose weight vector , employ different linguistic term sets to assess the suitability of the material with respect to the above selection attributes. Specifically, the decision-maker DM1 provides his assessments by using the linguistic term set ; DM2 provides his assessments using ; DM3 provides his assessments using ; DM4 provides his assessments using . These linguistic term sets are denoted as follows:

The linguistic assessments of the five alternatives on each attribute provided by the four decision-makers are presented in Table 1.

Then, we utilize the method being proposed to get the most desirable alternative(s).

Step 1. Convert the linguistic decision matrix shown in Table 1 into interval 2-tuple linguistic decision matrix , which is depicted in Table 2.

Step 2. Utilize the decision information given in matrixes and the I-ITHHM operator to get a collective decision matrix, as presented in Table 3.
In this example, the weight vector of the I-ITHHM operator is according to the normal distribution based method [40].

Step 3. Utilize the decision information given in matrix and the ITWHM operator to obtain the collective overall preference value of the alternative :

Step 4. Calculate the values of the collective overall preference values by using (10), and the results obtained are given as follows:

Step 5. Rank all the alternative in accordance with the values of in descending order: , and thus the most desirable alternative is .

Next, we use the ITHHM and the I-ITOWHM operators in Step to rank the alternatives. The ranking results are shown in Table 4. It is easy to see from Table 4 that the best selection is the alternative for all the aggregation operators used in this example. It should be noted that, depending on the aggregation operators used, the ranking orders of the alternative materials may be different. Therefore, according to the particular type of aggregation operator used, the results may lead to different decisions. But, in this example, it seems clear that is the optimal choice. Hence, the committee can recommend that material is the most suitable selection for the product.

In addition, to further evaluate the proposed method, we use the above material selection problem to analyze some comparable methods, which are based on the aggregation operators proposed by Xu [30], Park et al. [31], and Wan [15], respectively. In the first method, the fuzzy hybrid harmonic mean (FHHM) operator is utilized in the aggregation stage and the fuzzy weighted harmonic mean (FWHM) operator is used in the exploitation stage. Similarly, the second method employs the uncertain linguistic hybrid harmonic mean (ULHHM) operator in the aggregation stage and the uncertain linguistic weighted harmonic mean (ULWHM) operator in the exploitation stage; the third method employs the 2-tuple hybrid weighted arithmetic average (THWA) operator in the aggregation stage and the 2-tuple weighted averaging (TWA) operator in the exploitation stage. The ranking results of the four alternatives derived by using these methods are presented in Table 4.

According to Table 4, the ranking orders of the alternatives obtained by the method proposed in this paper are exactly the same as those determined by the methods based on previous operators. Thus, the proposed method is validated. However, in comparison with the listed methods, the proposed approach using interval 2-tuple linguistic harmonic mean operators is more reasonable and flexible for solving MAGDM problems because of the following:(i)It has exact characteristic in linguistic information processing and can effectively avoid the loss and distortion of information which occur formerly in other types of linguistic computational models.(ii)The uncertainty and diversity of decision-makers’ assessment information can be well reflected and modelled using interval 2-tuple linguistic variables. Moreover, the linguistic term sets with different granularity of uncertainty can be used by decision-makers for assessing alternatives.(iii)We can represent more complex group decision-making processes that include psychological factors such as time pressure and personal affects to each alternative, by using order inducing variables in the aggregation stage.

6. Conclusions

In this paper, we have developed some new harmonic aggregation operators including the interval 2-tuple hybrid harmonic mean (ITHHM) operator, the induced interval 2-tuple ordered weighted harmonic mean (I-ITOWHM) operator, and the induced interval 2-tuple hybrid harmonic mean (I-ITHHM) operator. It has been shown that both the ITWHM and ITOWHM operators are the special cases of the ITHHM operator, and if all the input interval 2-tuple data are reduced to the 2-tuple data, then the developed operators are reduced to the 2TLHH operator, the I-TOWHM operator, and the I-THHM operator, respectively. We have studied some desired properties of the I-ITOWHM operator, such as commutativity, idempotency, and monotonicity, and applied the I-ITHHM and the ITWHM operators to multiple attribute group decision-making with interval 2-tuple linguistic information. Finally, a material selection example has been given to verify the developed method and to demonstrate its practicality and effectiveness.

In the future, we expect to present further extensions to the proposed approach by adding new characteristics in the decision process and consider the potential applications of the developed interval 2-tuple linguistic harmonic mean operators to other fields. First, in many real-world situations, decision-makers may hesitate among several possible linguistic values or think of richer expressions for assessing an alternative because of uncertainty. Thus, extending the proposed decision model by using the hesitant fuzzy linguistic term sets [4446] is recommended in future research to solve the material selection problem more efficiently. Second, the proposed approach for group decision-making based on interval 2-tuple linguistic harmonic mean operators is a general method, which can be easily applied to deal with other decision-making problems such as robot evaluation and selection, green supply chain management, and electric vehicle charging station planning.

Competing Interests

The authors declare no competing interests.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (no. 71402090), the National Social Science Foundation (no. 15CGL003), the Program for Professor of Special Appointment (Young Eastern Scholar) at Shanghai Institutions of Higher Learning (no. QD2015019), and the Shanghai Science and Technology Innovation Action Plan Soft Science Foundation (no. 16692103800).