text:
Although historically materials discovery has been driven by a laborious trial-and-error process, knowledge-driven materials design can now be enabled by the rational combination of Machine Learning methods and materials databases. Here, data from the AFLOW repository for ab initio calculations is combined with Quantitative Materials Structure-Property Relationship models to predict important properties: metal/insulator classification, band gap energy, bulk/shear moduli, Debye temperature and heat capacities. The prediction’s accuracy compares well with the quality of the training data for virtually any stoichiometric inorganic crystalline material, reciprocating the available thermomechanical experimental data. The universality of the approach is attributed to the construction of the descriptors: Property-Labelled Materials Fragments. The representations require only minimal structural input allowing straightforward implementations of simple heuristic design rules.
Advances in materials science are often slow and fortuitous1. Coupling the field’s combinatorial challenges with the demanding efforts required for materials characterization makes progress uniquely difficult. The number of materials currently characterized, either experimentally2,3 or computationally4, pales in comparison with the anticipated potential diversity. Only considering naturally occurring elements, 9,000 crystal structure prototypes2,3, and stoichiometric compositions, there are roughly 3 × 1011 potential quaternary compounds and 1013 quinary combinations. Indeed, it has been estimated that the total number of theoretical materials can be as large as 10100 (ref. 5). Exacerbating the issue, standard materials characterization practices, such as calculating the band structure, can become expensive when considering finite-size scaling, charge corrections6, and going beyond standard density functional theory (DFT) with Green’s function methods such as the fully self-consistent GW approximation7,8. Ultimately, brute force exploration of this search space, even in high-throughput fashion1,9, is entirely impractical.
To circumvent the issue, many knowledge-based structure–property relationships have been conjectured over the years to aid in the search for novel functional materials–ranging from the simplest empirical relationships10 to complex advanced models11,12,13,14,15,16,17. For instance, many (semi-)empirical rules have been developed to predict band gap energies, such as those incorporating (optical18) electronegativity19. More sophisticated Machine Learning (ML) models were also developed for chalcopyrite semiconductors20, perovskites21, and binary compounds22. Unfortunately, many of these models are limited to a single family of materials, with narrow applicability outside of their training scope.
The development of such structure–property relationships has become an integral practice in the drug industry, which faces a similar combinatorial challenge. The number of potential organic molecules is estimated to be anywhere between 1013 and 10180 (ref. 23). In computational medicinal chemistry, Quantitative Structure-Activity Relationship modelling coupled with virtual screening of chemical libraries have been largely successful in the discovery of novel bioactive compounds24.
Here, we introduce fragment descriptors of materials structure. The combination of these descriptors with ML approaches affords the development of universal models capable of accurate prediction of properties for virtually any stoichiometric inorganic crystalline material. First, the algorithm for descriptor generation is described, along with implementation of ML methods for Quantitative Materials Structure-Property Relationship (QMSPR) modelling. Next, the effectiveness of this approach is assessed through prediction of eight critical electronic and thermomechanical properties of materials, including the metal/insulator classification, band gap energy, bulk and shear moduli, Debye temperature, heat capacities (at constant pressure and volume) and thermal expansion coefficient. The impact and interaction among the most significant descriptors as determined by the ML algorithms are highlighted. As a proof-of-concept, the QMSPR models are then employed to predict thermomechanical properties for compounds previously uncharacterized, and the predictions are validated via the AEL–AGL integrated framework (Automatic Elasticity Library-Automatic GIBBS Library)25,26. Such predictions are of particular value as proper calculation pathways for thermomechanical properties in the most efficient scenarios still require analysis of multiple DFT-runs, elevating the cost of already expensive calculations. Finally, ML-predictions and calculations are both compared to experimental values which ultimately corroborate the validity of the approach.
Other investigations have predicted a subset of the target properties discussed here by building ML approaches where computationally obtained quantities, such as the cohesive energy, formation energy and energy above the convex hull, are part of the input data27. The approach presented here is orthogonal. Once trained, our proposed models achieve comparable accuracies without the need of further ab initio data. All necessary input properties are either tabulated or derived directly from the geometrical structures. There are advantages: (i) a priori, after the training, no further calculations need to be performed, (ii) a posteriori, the modelling framework becomes independent of the source or nature of the training data, for example, calculated versus experimental. The latter allows for rapid extension of predictions to online applications—given the geometry of a cell and the species involved, eight ML predicted properties are returned (aflow.org/aflow-ml).
Many cheminformatics investigations have demonstrated the critical importance of molecular descriptors, which are known to influence model accuracy more than the choice of the ML algorithm28,29. For the purposes of this investigation, fragment descriptors typically used for organic molecules were adapted to serve for materials characterization30. Molecular systems can be described as graphs whose vertices correspond to atoms and edges to chemical bonds. In this representation, fragment descriptors characterize subgraphs of the full 3D molecular network. Any molecular graph invariant can be uniquely represented as a linear combination of fragment descriptors. They offer several advantages over other types of chemical descriptors31, including simplicity of calculation, storage and interpretation32. However, they also come with a few disadvantages. Models built with fragment descriptors perform poorly when presented with new fragments for which they were not trained. In addition, typical fragments are constructed solely with information of the individual atomic symbols (for example, C, N, Na). Such a limited context would be insufficient for modelling the complex chemical interactions within materials.
Mindful of these constraints, fragment descriptors for materials were conceptualized by differentiating atoms not by their symbols but by a plethora of well-tabulated chemical and physical properties33. Descriptor features comprise of various combinations of these atomic properties. From this perspective, materials can be thought of as ‘coloured’ graphs, with vertices decorated according to the nature of the atoms they represent34. Partitions of these graphs form Property-Labelled Materials Fragments (PLMF).
Figure 1 shows the scheme for constructing PLMFs. Given a crystal structure, the first step is to determine the atomic connectivity within it. In general, atomic connectivity is not a trivial property to determine within materials. Not only must the potential bonding distances among atoms be considered, but also whether the topology of nearby atoms allows for bonding. Therefore, a computational geometry approach is employed to partition the crystal structure (Fig. 1a) into atom-centred Voronoi-Dirichlet polyhedra35,36 (Fig. 1b). This partitioning scheme was found to be invaluable in the topological analysis of metal organic frameworks, molecules, and inorganic crystals37. Connectivity between atoms is established by satisfying two criteria: (i) the atoms must share a Voronoi face (perpendicular bisector between neighbouring atoms), and (ii) the interatomic distance must be shorter than the sum of the Cordero covalent radii38 to within a 0.25 Å tolerance. Here, only strong interatomic interactions are modelled, such as covalent, ionic, and metallic bonding, ignoring van der Waals interactions. Owing to the ambiguity within materials, the bond order (single/double/triple bond classification) is not considered. Taken together, the Voronoi centres that share a Voronoi face and are within the sum of their covalent radii form a three-dimensional graph defining the connectivity within the material.
Figure 1: Schematic representing the construction of the Property-Labelled Materials Fragments (PLMF).
The crystal structure (a) is analysed for atomic neighbours via Voronoi tessellation (b). After property labelling, the resulting periodic graph (c) is decomposed into simple subgraphs (d).
In the final steps of the PLMF construction, the full graph and corresponding adjacency matrix (Fig. 1c) are constructed from the total list of connections. The adjacency matrix A of a simple graph (material) with n vertices (atoms) is a square matrix (n × n) with entries aij=1 if atom i is connected to atom j, and aij=0 otherwise. This adjacency matrix reflects the global topology for a given system, including interatomic bonds and contacts within the crystal. The full graph is partitioned into smaller subgraphs, corresponding to individual fragments (Fig. 1d). Although there are several subgraphs to consider in general, the length l is restricted to a maximum of three, where l is the largest number of consecutive, non-repetitive edges in the subgraph. This restriction serves to curb the complexity of the final descriptor vector. In particular, there are two types of fragments. Path fragments are subgraphs of at most l=3 that encode any linear strand of up to four atoms. Only the shortest paths between atoms are considered. Circular fragments are subgraphs of l=2 that encode the first shell of nearest neighbour atoms. In this context, circular fragments represent coordination polyhedra, or clusters of atoms with anion/cation centres each surrounded by a set of its respective counter ion. Coordination polyhedra are used extensively in crystallography and mineralogy39.
The PLMFs are differentiated by local (standard atomic/elemental) reference properties33, which include: (i) general properties: the Mendeleev group and period numbers (gP, pP), number of valence electrons (NV); (ii) measured properties33: atomic mass (matom), electron affinity (EA), thermal conductivity (λ), heat capacity (C), enthalpies of atomization (ΔHat), fusion (ΔHfusion) and vaporization (ΔHvapor), first three ionization potentials (IP1,2,3); and (iii) derived properties: effective atomic charge (Zeff), molar volume (Vmolar), chemical hardness (η)33,40, covalent (rcov)38, absolute41, and van der Waals radii33, electronegativity (χ) and polarizability (αP). Pairs of properties are included in the form of their multiplication and ratio, as well as the property value divided by the atomic connectivity (number of neighbours in the adjacency matrix). For every property scheme q, the following quantities are also considered: minimum (min(q)), maximum (max(q)), total sum (∑q), average (avg(q)) and standard deviation (std(q)) of q among the atoms in the material.
To incorporate information about shape, size and symmetry of the crystal unit cell, the following crystal-wide properties are incorporated: lattice parameters (a, b, c), their ratios (a/b, b/c, a/c), angles (α, β, γ), density, volume, volume per atom, number of atoms, number of species (atom types), lattice type, point group and space group.
All aforementioned descriptors (fragment-based and crystal-wide) can be concatenated together to represent each material uniquely. After filtering out low variance (<0.001) and highly correlated (r2>0.95) features, the final feature vector captures 2,494 total descriptors.
where the first set of indices count over all pairs of atoms and the second is restricted to all pairs i, j of bonded atoms.
In training the models, the same ML method and descriptors are employed without any hand tuning or variable selection. Specifically, models are constructed using the gradient boosting decision tree (GBDT) technique44. All models were validated through y-randomization (label scrambling). Five-fold cross validation is used to assess how well each model will generalize to an independent data set. Hyperparameters are determined with grid searches on the training set and 10-fold cross validation.
The GBDT method44 evolved from the application of boosting methods45 to regression trees46. The boosting method is based on the observation that finding many weakly accurate prediction rules can be a lot easier than finding a single, highly accurate rule47. The boosting algorithm calls this ‘weak’ learner repeatedly, at each stage feeding it a different subset of the training examples. Each time it is called, the weak learner generates a new weak prediction rule. After many iterations, the boosting algorithm combines these weak rules into a single prediction rule aiming to be much more accurate than any single weak rule.
where h(x; am) are the weak learners (decision trees in this case) characterized by parameters am, and M is the total count of decision trees obtained through boosting.
Gradient boosting attempts to solve this minimization problem numerically via steepest descent. The steepest descent direction is the negative gradient of the loss function evaluated at the current model Fm−1, where the step length is chosen using line search.
The influences are normalized so that they add to one. Influences capture the importance of the variable, but not the direction of the response (positive or negative).
Eight predictive models are developed in this work, including: a binary classification model that predicts if a material is a metal or an insulator and seven regression models that predict: the band gap energy (EBG) for insulators, bulk modulus (BVRH), shear modulus (GVRH), Debye temperature (θD), heat capacity at constant pressure (Cp), heat capacity at constant volume (CV), and thermal expansion coefficient (αV).
Figure 2 shows the overall application work-flow. A novel candidate material is first classified as a metal or an insulator. If the material is classified as an insulator, EBG is predicted, whereas classification as a metal implies that the material has no EBG. The six thermomechanical properties are then predicted independent of the material’s metal/insulator classification. The integrated modelling work-flow has been implemented as a web application at aflow.org/aflow-ml, requiring only the atomic species and positions as input for predictions.
Figure 2: Outline of the modelling work-flow.
ML models are represented by orange diamonds. Target properties predicted by these models are highlighted in green.
Although all three models were trained independently, the accuracy of the EBG regression model is inherently dependent on the accuracy of the metal/insulator classification model in this work-flow. However, the high accuracy of the metal/insulator classification model suggests this not to be a practical concern.
One technique for assessing model quality is fivefold cross validation, which gauges how well the model is expected to generalize to an independent data set. For each model, the scheme involves randomly partitioning the set into five groups and predicting the value of each material in one subset while training the model on the other four subsets. Hence, each subset has the opportunity to play the role of the ‘test set’. Furthermore, any observed deviations in the predictions are addressed. For further analysis, all predicted and calculated results are available in Supplementary Note 2.
The accuracy of the metal/insulator classifier is reported as the area under the curve (AUC) of the receiver operating characteristic (ROC) plot (Fig. 3a). The ROC curve illustrates the model’s ability to differentiate between metallic and insulating input materials. It plots the prediction rate for insulators (correctly versus incorrectly predicted) throughout the full spectrum of possible prediction thresholds. An area of 1.0 represents a perfect test, whereas an area of 0.5 characterizes a random guess (the dashed line). The model shows excellent external predictive power with the area under the curve at 0.98, an insulator-prediction success rate (sensitivity) of 0.95, a metal-prediction success rate (specificity) of 0.92, and an overall classification rate of 0.93. For the complete set of 26,674 materials, this corresponds to 2,103 misclassified materials, including 1,359 misclassified metals and 744 misclassified insulators. Evidently, the model exhibits positive bias toward predicting insulators, where bias refers to whether a ML model tends to over- or under-estimate the predicted property. This low false-metal rate is fortunate as the model is unlikely to misclassify a novel, potentially interesting semiconductor as a metal. Overall, the metal classification model is robust enough to handle the full complexity of the periodic table.
Figure 3: Five-fold cross validation plots for the eight ML models predicting electronic and thermomechanical properties.
(a) Receiver operating characteristic (ROC) curve for the classification ML model. (b–h) Predicted versus calculated values for the regression ML models: (b) band gap energy (EBG), (c) bulk modulus (BVRH), (d) shear modulus (GVRH), (e) Debye temperature (θD), (f) heat capacity at constant pressure (CP), (g) heat capacity at constant volume (CV) and (h) thermal expansion coefficient (αV).
The results of the fivefold cross validation analysis for the band gap energy (EBG) regression model are plotted in Fig. 3b. In addition, a statistical profile of these predictions, along with that of the six thermomechanical regression models, is provided in Table 1, which includes metrics such as the root-mean-square error (RMSE), mean absolute error (MAE), and coefficient of determination (r2). Similar to the classification model, the EBG model exhibits a positive predictive bias. The biggest errors come from materials with narrow band gaps, that is, the scatter in the lower left corner in Fig. 3b. These materials predominantly include complex fluorides and nitrides. N2H6Cl2 (ICSD #23145) exhibits the worst prediction accuracy with signed error SE=3.78 eV48. The most underestimated materials are HCN (ICSD #76419) and N2H6Cl2 (ICSD #240903) with SE=−2.67 and −3.19 eV49,50, respectively. This is not surprising considering that all three are molecular crystals. Such systems are anomalies in the ICSD, and fit better in other databases, such as the Cambridge Structural Database51. Overall, 10,762 materials are predicted within 25% accuracy of calculated values, whereas 824 systems have errors over 1 eV.
Table 1: Statistical summary of the fivefold cross validated predictions for the seven regression models.
Figure 3c–h and Table 1 showcase the results of the fivefold cross validation analysis for the six thermomechanical regression models. For both bulk (BVRH) and shear (GVRH) moduli, over 85% of materials are predicted within 20 GPa of their calculated values. The remaining models also demonstrate high accuracy, with at least 90% of the full training set (>2,546 systems) predicted to within 25% of the calculated values. Significant outliers in predictions of the bulk modulus include graphite (ICSD #187640, SE=100 GPa, likely due to extreme anisotropy) and two theoretical high-pressure boron nitrides (ICSD #162873 and #162874, under-predicted by over 110 GPa)52,53. Other theoretical systems are ill-predicted throughout the six properties, including ZN (ICSD #161885), CN2 (ICSD #247676), C3N4 (ICSD #151782) and CH (ICSD #187642)52,54,55,56. Predictions for the GVRH, Debye temperature (θD), and thermal expansion coefficient (αV) tend to be slightly underestimated, particularly for higher calculated values. In addition, mild scattering can be seen for θD and θV, but not enough to have a significant impact on the error or correlation metrics.
Despite minimal deviations, both RMSE and mean absolute error are within 4% of the ranges covered for each property, and the predictions demonstrate excellent correlation with the calculated properties. Note the tight clustering of points just below 3 kB per atom for the heat capacity at constant volume (CV). This is due to CV saturation in accordance with the Dulong-Petit law occurring at or below 300 K for many compounds.
Model interpretation is of paramount importance in any ML study. The significance of each descriptor is determined in order to gain insight into structural features that impact molecular properties of interest. Interpretability is a strong advantage of decision tree methods, particularly with the GBDT approach. One can quantify the predictive power of a specific descriptor by analysing the reduction of the RMSE at each node of the tree.
Partial dependence plots offer yet another opportunity for GBDT model interpretation. Similar to the descriptor significance analysis, partial dependence resolves the effect of a variable (descriptor) on a property, but only after marginalising over all other explanatory variables57. The effect is quantified by the change of that property as relevant descriptors are varied. The plots themselves highlight the most important interactions among relevant descriptors as well as between properties and their corresponding descriptors. Although only the most important descriptors are highlighted and discussed, an exhaustive list of relevant descriptors and their relative contributions can be found in Supplementary Note 1.
For the metal/insulator classification model, the descriptor significance analysis shows that two descriptors have the highest importance (equally), namely and . is the ratio between the fusion enthalpy (ΔHfusion) and the thermal conductivity (λ) averaged over all atoms in the material, and is the ratio between the molar volume (Vmolar) and the covalent radius (rcov) averaged over all atoms in the material. Both descriptors are simple node-specific features. The presence of these two prominent descriptors accounts for the high accuracy of the classification model.
Figure 4 shows the projection of the full dataset onto the dual-descriptor space of and . In this 2D space, metals and insulators are substantially partitioned. To further resolve this separation, the plot is split into four quadrants (see dashed lines) with an origin approximately at , . Insulators are predominately located in quadrant I. There are several clusters (one large and several small) parallel to the x axis. Metals occupy a compact square block in quadrant III within intervals and . Quadrant II is mostly empty with a few materials scattered about the origin. In the remaining quadrant (IV), materials have mixed character.
Figure 4: Semi-log scatter plot of the full data set (26,674 unique materials) in a dual-descriptor space.
avg(ΔHfusionλ−1) versus . Insulators and metals are coloured in red and blue, respectively.
Analysis of the projection shown in Fig. 4 suggests a simple heuristic rule: all materials within quadrant I are classified as insulators (EBG>0), and all materials outside of this quadrant are metals. Remarkably, this unsupervised projection approach achieves a very high classification accuracy of 86% for the entire dataset of 26,674 materials. The model misclassifies only 3,621 materials: 2,414 are incorrectly predicted as insulators and 1,207 are incorrectly predicted as metals. This example illustrates how careful model analysis of the most significant descriptors can yield simple heuristic rules for materials design.
The regression model for the band gap energy (EBG) is more complex. There are a number of descriptors in the model with comparable contributions, and thus, all individual contributions are small. This is expected as a number of conditions can affect EBG. The most important are and avg(Cλ−1) with significance scores of 0.075 and 0.071, respectively, where χ is the electronegativity, Zeff is the effective nuclear charge, C is the specific heat capacity and λ is the thermal conductivity of each atom.
Figure 5 shows partial dependence plots focusing on (ΔIPbond) as an example. It is derived from edge fragments of bonded atoms (l=1) and defined as an absolute difference in ionization potentials averaged over the material. In other words, it is a measure of bond polarity, similar to electronegativity. Figure 5a shows a steady monotonic increase in EBG for larger values of (ΔIPbond). The effect is small, but captures an expected physical principle: polar inorganic materials (for example, oxides, fluorides) tend to have larger EBG.
Figure 5: Partial dependence plots of the EBG, BVRH and θD models.
(a) Partial dependence of EBG on the avg(ΔIPbond) descriptor. For EBG, the 2D interaction between std(ΔIPbond) and avg(ΔIPbond) and between ρ (density) and avg(ΔIPbond) are illustrated in panels (b,c), respectively. (d) Partial dependence of the BVRH on the crystal volume per atom descriptor. For θD, the 2D interaction between avg(ΔEAbond) and and between crystal lattice parameters b and c are illustrated in panels (e,f), respectively.
Given the number of significant interactions involved with this phenomenon, tailoring EBG involves the optimization of a highly non-convex, multidimensional object. Figure 5b illustrates a 2D slice of this object as std(ΔIPbond) and avg(ΔIPbond) vary simultaneously. Like avg(ΔIPbond), std(ΔIPbond) is the s.d. of the set of absolute differences in IP among all bonded atoms. In the context of these two variables, EBG responds to deviations in ΔIPbond among the set of bonded atoms, but remains constant across shifts in avg(ΔIPbond). This suggests an opportunity to tune EBG by considering another composition that varies the deviations among bond polarities. Alternatively, a desired EBG can be maintained by considering another composition that preserves the deviations among bond polarities, even as the overall average shifts. Similarly, Fig. 5c shows the partial dependence on both the density (ρ) and avg(ΔIPbond). Contrary to the previous trend, larger avg(ΔIPbond) values correlate with smaller EBG, particularly for low density structures. Materials with higher density and lower avg(ΔIPbond) tend to have higher EBG. Considering the elevated response (compared with Fig. 5b), the inverse correlation of EBG with the average bond polarity in the context of density suggests an even more effective means of tuning EBG.
A descriptor analysis of the thermomechanical property models reveals the importance of one descriptor in particular, the volume per atom of the crystal. This conclusion certainly resonates with the nature of these properties, as they generally correlate with bond strength26. Figure 4d exemplifies such a relationship, which shows the partial dependence plot of the bulk modulus (BVRH) on the volume per atom. Tightly bound atoms are generally indicative of stronger bonds. As the interatomic distance increases, properties like BVRH generally reduce.
Two of the more interesting dependence plots are also shown in Fig. 5e,f, both of which offer opportunities for tuning the Debye temperature (θD). Figure 5e illustrates the interactions among two descriptors, the absolute difference in electron affinities among bonded atoms averaged over the material (avg(ΔEAbond)), and the s.d. of the set of ratios of the enthalpies of vaporization (ΔHvapor) and atomization ΔHatom) for all atoms in the material . Within these dimensions, two distinct regions emerge of increasing/decreasing θD separated by a sharp division at about . Within these partitions, there are clusters of maximum gradient in θD—peaks within the left partition and troughs within the right. The peaks and troughs alternate with varying . Although is not an immediately intuitive descriptor, the alternating clusters may be a manifestation of the periodic nature of ΔHvapor and ΔHatom (ref. 58). As for the partitions themselves, the extremes of avg(ΔEAatom) characterize covalent and ionic materials, as bonded atoms with similar EA are likely to share electrons, whereas those with varying EA prefer to donate/accept electrons. Considering that EA is also periodic, various opportunities for carefully tuning θD should be available.
Finally, Fig. 5f shows the partial dependence of θD on the lattice parameters b and c. It resolves two notable correlations: (i) uniformly increasing the cell size of the system decreases θD, but (ii) elongating the cell increases it. Again, (i) can be attributed to the inverse relationship between volume per atom and bond strength, but does little to address (ii). Nevertheless, the connection between elongated, or layered, systems and the Debye temperature is certainly not surprising—anisotropy can be leveraged to enhance phonon-related interactions associated with thermal conductivity59 and superconductivity60,61,62. Although the domain of interest is quite narrow, the impact is substantial, particularly in comparison with that shown in Fig. 5e.
Although the expected performances of the ML models can be projected through fivefold cross validation, there is no substitute for validation against an independent dataset. The ML models for the thermomechanical properties were leveraged to make predictions for materials previously uncharacterized, and these predictions were subsequently validated via the AEL–AGL integrated framework25,26. Figure 6 illustrates the models’ performance on the set of 770 additional materials, with relevant statistics displayed in Table 2. For further analysis, all predicted and calculated results are available in Supplementary Note 3.
Figure 6: Model performance evaluation for the six ML models predicting thermomechanical properties of 770 characterized materials.
Predicted versus calculated values for the regression ML models: (a) bulk modulus (BVRH), (b) shear modulus (GVRH), (c) Debye temperature (θD), (d) heat capacity at constant pressure (CP), (e) heat capacity at constant volume (CV), and (f) thermal expansion coefficient (αV).
Table 2: Statistical summary of the predictions for the six thermomechanical regression models.
Comparing with the results of the generalizability analysis shown in Fig. 3 and Table 1, the overall errors are consistent with fivefold cross validation. Five out of six models have r2 of 0.9 or higher. However, the r2 value for the thermal expansion coefficient (αV) is lower than forecasted. The presence of scattering suggests the need for a larger training set—as new, much more diverse materials were likely introduced in the test set. This is not surprising considering the number of variables that can affect thermal expansion63. Otherwise, the accuracy of these predictions confirm the effectiveness of the PLMF representation, which is particularly compelling considering: (i) the limited diversity training dataset (only ∼11% as large as that available for predicting the electronic properties) and (ii) the relative size of the test set (over a quarter the size of the training set).
In the case of the bulk modulus (BVRH), 665 systems (86% of test set) are predicted within 25% of calculated values. Only the predictions of four materials, Bi (ICSD #51674), PrN (ICSD #168643), Mg3Sm (ICSD #104868), and ZrN (ICSD #161885), deviate beyond 100 GPa from calculated values. Bi is a high-pressure phase (Bi-III) with a caged, zeolite-like structure64. The structures of zirconium nitride (wurtzite phase) and praseodymium nitride (B3 phase) were hypothesized and investigated via DFT calculations54,65 and have yet to be observed experimentally.
For the shear modulus (GVRH) 482 materials (63% of the test set) are predicted within 25% of calculated values. Just one system, C3N4 (ICSD #151781), deviates beyond 100 GPa from its calculated value. The Debye temperature (θD) is predicted to within 50 K accuracy for 540 systems (70% of the test set). BeF2 (ICSD #173557), yet another cage (sodalite) structure66, has among the largest errors in three models including θD (SE=−423 K) and both heat capacities (Cp: SE=0.65 kB per atom; CV: SE=0.61 kB per atom). Similar to other ill-predicted structures, this polymorph is theoretical, and has yet to be synthesized.
A comparison between calculated, predicted and experimental results is presented in Fig. 7, with relevant statistics summarized in Table 3. Data are considered for the bulk modulus B, shear modulus G, and (acoustic) Debye temperature θa for 45 well-characterized materials with diamond (SG# 227, AFLOW prototype A_cF8_227_a), zincblende (SG# 216, AB_cF8_216_c_a), rocksalt (SG# 225, AB_cF8_225_a_b), and wurtzite (SG# 186, AB_hP4_186_b_b) structures67,68. Experimental B and G are compared with the BVRH and GVRH values predicted here, and θa is converted to the traditional Debye temperature , where n is the number of atoms in the unit cell. All relevant values are listed in Supplementary Note 4.
Figure 7: Comparison of the AEL—AGL calculations and ML predictions with experimental values for three thermomechanical properties.
(a) bulk modulus (B), (b) shear modulus (G), and (c) Debye temperature (θD).
Table 3: Statistical summary of the AEL–AGL calculations and ML predictions versus experimental values for three thermomechanical properties.
Excellent agreement is found between experimental and calculated values, but more importantly, between experimental and predicted results. With error metrics close to or under expected tolerances from the generalizability analysis, the comparison highlights effective experimental confidence in the approach. The experiments/prediction validation is clearly the ultimate objective of the research presented here.
Traditional trial-and-error approaches have proven ineffective in discovering practical materials. Computational models developed with ML techniques may provide a truly rational approach to materials design. Typical high-throughput DFT screenings involve exhaustive calculations of all materials in the database, often without consideration of previously calculated results. Even at high-throughput rates, an average DFT calculation of a medium size structure (∼50 atoms per unit cell) takes ∼1,170 CPU-hours of calculations or about 37 h on a 32-CPU cores node. However, in many cases, the desired range of values for the target property is known. For instance, the optimal band gap energy and thermal conductivity for optoelectronic applications will depend on the power and voltage conditions of the device63,69. Such cases offer an opportunity to leverage previous results and savvy ML models, such as those developed in this work, for rapid pre-screening of potential materials. Researchers can quickly narrow the list of candidate materials and avoid many extraneous DFT calculations—saving money, time and computational resources. This approach takes full advantage of previously calculated results, continuously accelerating materials discovery. With prediction rates of about 0.1 s per material, the same 32-CPU cores node can screen over 28 million material candidates per day with this framework.
Furthermore, interaction diagrams as depicted in Fig. 5 offer a pathway to design materials that meet certain constraints and requirements. For example, substantial differences in thermal expansion coefficients among the materials used in high-power, high-frequency optoelectronic applications leads to bending and cracking of the structure during the growth process63,69. Not only would this work-flow facilitate the search for semiconductors with large band gap energies, high Debye temperatures (thermal conductivity), but also materials with similar thermal expansion coefficients.
Although the models themselves demonstrate excellent predictive power with minor deviations, outlier analysis reveals theoretical structures to be among the worst offenders. This is not surprising, as the true stability conditions (for example, high-pressure/high-temperature) have yet to be determined, if they exist at all. The ICSD estimates that structures for over 7,000 materials (or roughly 4%) come from calculations rather than actual experiment. Such discoveries exemplify yet another application for ML modelling, rapid/robust curation of large data sets.
To improve large-scale high-throughput computational screening for the identification of materials with desired properties, fast and accurate data mining approaches should be incorporated into the standard work-flow. In this work, we developed a universal QMSPR framework for predicting electronic properties of inorganic materials. Its effectiveness is validated through the prediction of eight key materials properties for stoichiometric inorganic crystalline materials, including the metal/insulator classification, band gap energy, bulk and shear moduli, Debye temperature, heat capacity (at constant pressure and volume) and thermal expansion coefficient. Its applicability extends to all 230 space groups and the vast majority of elements in the periodic table. All models are freely available at aflow.org/aflow-ml.
Two independent data sets were prepared for the creation and validation of the ML models. The training set includes electronic4,70,71,72,73,74 and thermomechanical properties25,26 for a broad diversity of compounds already characterized in the AFLOW database. This set is used to build and analyse the ML models, one model per property. The constructed thermomechanical models are then employed to make predictions of previously uncharacterized compounds in the AFLOW database. Based on these predictions and consideration of computational cost, several compounds are selected to validate the models’ predictive power. These compounds and their computed properties define the test set. The compounds used in both data sets are specified in Supplementary Notes 2 and 3, respectively.
Band gap energy data for 49,934 materials were extracted from the AFLOW repository4,70,71,72,73,74, representing ∼60% of the known stoichiometric inorganic crystalline materials listed in the Inorganic Crystal Structure Database (ICSD)2,3. Although these band gap energies are generally underestimated with respect to experimental values75, DFT+U is robust enough to differentiate between metallic (no EBG) and insulating (EBG>0) systems76. In addition, errors in band gap energy prediction are typically systematic. Therefore, the band gap energy values can be corrected ad hoc with fitting schemes77,78. Prior to model development, both ICSD and AFLOW data were curated: duplicate entries, erroneous structures, and ill-converged calculations were corrected or removed. Noble gases crystals are not considered. The final data set consists of 26,674 unique materials (12,862 with no EBG and 13,812 with EBG>0), covering the seven lattice systems, 230 space groups and 83 elements (H-Pu, excluding noble gases, Fr, Ra, Np, At and Po). All referenced DFT calculations were performed with the Generalized Gradient Approximation PBE79 exchange-correlation functional and projector-augmented wavefunction potentials80,81 according to the AFLOW Standard for High-Throughput Computing76. The Standard ensures reproducibility of the data, and provides visibility/reasoning for any parameters set in the calculation, such as accuracy thresholds, calculation pathways, and mesh dimensions. II. Thermomechanical properties data for just over 3,000 materials were extracted from the AFLOW repository26. These properties include the bulk modulus, shear modulus, Debye temperature, heat capacity at constant pressure, heat capacity at constant volume, and thermal expansion coefficient, and were calculated using the AEL–AGL integrated framework25,26. The AEL (AFLOW Elasticity Library) method26 applies a set of independent normal and shear strains to the structure, and then fits the calculated stress tensors to obtain the elastic constants82. These can then be used to calculate the elastic moduli in the Voigt and Reuss approximations, as well as the Voigt-Reuss-Hill (VRH) averages that are the values of the bulk and shear moduli modelled in this work. The AGL (AFLOW GIBBS Library) method25 fits the energies from a set of isotropically compressed and expanded volumes of a structure to a quasiharmonic Debye-Grüneisen model83 to obtain thermomechanical properties, including the bulk modulus, Debye temperature, heat capacity and thermal expansion coefficient. AGL has been combined with AEL in a single work-flow, so that it can utilize the Poisson ratios obtained from AEL to improve the accuracy of the thermal properties predictions26. After a similar curation of ill-converged calculations, the final data set consists of 2,829 materials. It covers the seven lattice systems, includes unary, binary and ternary compounds, and spans broad ranges of each thermomechanical property, including high thermal conductivity systems such as C (ICSD #182729), BN (ICSD #162874), BC5 (ICSD #166554), CN2 (ICSD #247678), MnB2 (ICSD #187733) and SiC (ICSD #164973), as well as low thermal conductivity systems such as Hg33(Rb,K)3 (ICSD #410567 and #410566), Cs6Hg40 (ICSD #240038), Ca16Hg36 (ICSD #107690), CrTe (ICSD #181056) and Cs (ICSD #426937). Many of these systems additionally exhibit extreme values of the bulk and shear moduli, such as C (high bulk and shear moduli) and Cs (low bulk and shear moduli). Interesting systems such as RuC (ICSD #183169) and NbC (ICSD #189090) with a high bulk modulus (BVRH=317.92 GPa, 263.75 GPa) but low shear modulus (GVRH=16.11 GPa, 31.86 GPa) also populate the set.
Although nearly all ICSD compounds are characterized electronically within the AFLOW database, most have not been characterized thermomechanically owing to the added computational cost. This presented an opportunity to validate the ML models. Of the remaining compounds, several were prioritized for immediate characterization via the AEL–AGL integrated framework25,26. In particular, focus was placed on systems predicted to have a large bulk modulus, as this property is expected to scale well with the other aforementioned thermomechanical properties25,26. The set also includes various other small cell, high symmetry systems expected to span the full applicability domains of the models. This effort resulted in the characterization of 770 additional compounds.
All the ab initio data are freely available to the public as part of the AFLOW online repository and can be accessed through aflow.org following the REST-API interface70.
How to cite this article: Isayev, O. et al. Universal fragment descriptors for predicting properties of inorganic crystals. Nat. Commun. 8, 15679 doi: 10.1038/ncomms15679 (2017).
Curtarolo, S. et al. The high-throughput highway to computational materials design . Nat. Mater. 12, 191–201 (2013).
Bergerhoff, G., Hundt, R., Sievers, R. & Brown, I. D. The inorganic crystal structure data base . J. Chem. Inf. Comput. Sci. 23, 66–69 (1983).
Belsky, A., Hellenbrandt, M., Karen, V. L. & Luksch, P. New developments in the Inorganic Crystal Structure Database (ICSD): accessibility in support of materials research and design . Acta Crystallogr. Sect. B 58, 364–369 (2002).
Curtarolo, S. et al. AFLOWLIB.ORG: a distributed materials properties repository from high-throughput ab initio calculations . Comput. Mater. Sci. 58, 227–235 (2012).
Walsh, A. Inorganic materials: the quest for new functionality . Nat. Chem. 7, 274–275 (2015).
Castleton, C. W. M., Höglund, A. & Mirbt, S. Managing the supercell approximation for charged defects in semiconductors: Finite-size scaling, charge correction factors, the band-gap problem, and the ab initio dielectric constant . Phys. Rev. B 73, 035215 (2006).
Lindgren, I. in Relativistic Many-Body Theory: A New Field-Theoretical Approach Vol. 63, 199–210Springer (2011).
van Schilfgaarde, M., Kotani, T. & Faleev, S. Quasiparticle Self-Consistent GW Theory . Phys. Rev. Lett. 96, 226402 (2006).
Koinuma, H. & Takeuchi, I. Combinatorial solid-state chemistry of inorganic materials . Nat. Mater. 3, 429–438 (2004).
Mizutani, U. Hume-Rothery Rules for Structurally Complex Alloy Phases CRC Press (2011).
Isayev, O. et al. Materials cartography: representing and mining materials space using structural and electronic fingerprints . Chem. Mater. 27, 735–743 (2015).
Ghiringhelli, L. M., Vybiral, J., Levchenko, S. V., Draxl, C. & Scheffler, M. Big data of materials science: critical role of the descriptor . Phys. Rev. Lett. 114, 105503 (2015).
Pyzer-Knapp, E. O., Li, K. & Aspuru-Guzik, A. Learning from the Harvard Clean Energy Project: the use of neural networks to accelerate materials discovery . Adv. Func. Mater 25, 6495–6502 (2015).
Rajan, K. Materials informatics: the materials ‘gene’ and big data . Annu. Rev. Mater. Res. 45, 153–169 (2015).
Carrete, J., Mingo, N., Wang, S. & Curtarolo, S. Nanograined half-heusler semiconductors as advanced thermoelectrics: an ab initio high-throughput statistical study . Adv. Func. Mater 24, 7427–7432 (2014).
van Roekeghem, A., Carrete, J., Oses, C., Curtarolo, S. & Mingo, N. High-throughput computation of thermal conductivity of high-temperature solid phases: the case of oxide and fluoride perovskites . Phys. Rev. X 6, 041061 (2016).
Furmanchuk, A., Agrawal, A. & Choudhary, A. Predictive analytics for crystalline materials: bulk modulus . RSC Adv 6, 95246–95251 (2016).
Duffy, J. A. Variable electronegativity of oxygen in binary oxides: possible relevance to molten fluorides . J. Chem. Phys. 67, 2930–2931 (1977).
Di Quarto, F., Sunseri, C., Piazza, S. & Romano, M. C. Semiempirical correlation between optical band gap values of oxides and the difference of electronegativity of the elements. its importance for a quantitative use of photocurrent spectroscopy in corrosion studies . J. Phys. Chem. B 101, 2519–2525 (1997).
Zeng, Y., Chua, S. J. & Wu, P. On the prediction of ternary semiconductor properties by artificial intelligence methods . Chem. Mater. 14, 2989–2998 (2002).
Pilania, G. et al. Machine learning bandgaps of double perovskites . Sci. Rep 6, 19375 (2016).
Gu, T., Lu, W., Bao, X. & Chen, N. Using support vector regression for the prediction of the band gap and melting point of binary and ternary compound semiconductors . Solid State Sci. 8, 129–136 (2006).
Gorse, A.-D. Diversity in medicinal chemistry space . Curr. Top. Med. Chem. 6, 3–18 (2006).
Varnek A., Tropsha A. (eds) Chemoinformatics Approaches to Virtual Screening RSC (2008).
Toher, C. et al. High-throughput computational screening of thermal conductivity, Debye temperature, and Grüneisen parameter using a quasiharmonic Debye Model . Phys. Rev. B 90, 174107 (2014).
Toher, C. et al. Combining the AFLOW GIBBS and Elastic Libraries to efficiently and robustly screen thermomechanical properties of solids . Phys. Rev. Mater. (in the press).
de Jong, M. et al. A statistical learning framework for materials science: application to elastic moduli of k-nary inorganic polycrystalline compounds . Sci. Rep. 6, 34256 (2016).
Young, S. S., Yuan, F. & Zhu, M. Chemical descriptors are more important than learning algorithms for modelling . Mol. Informatics 31, 707–710 (2012).
Polishchuk, P. G., Kuz’min, V. E., Artemenko, A. G. & Muratov, E. N. Universal approach for structural interpretation of QSAR/QSPR models . Mol. Informatics 32, 843–853 (2013).
Ruggiu, F., Marcou, G., Varnek, A. & Horvath, D. ISIDA property-labelled fragment descriptors . Mol. Informatics 29, 855–868 (2010).
Todeschini, R. & Consonni, V. Handbook of molecular descriptors . Methods and Principles in Medicinal Chemistry Wiley-VCH Verlag GmbH (2000).
Varnek, A., Fourches, D., Hoonakker, F. & Solov’ev, V. P. Substructural fragments: an universal language to encode reactions, molecular and supramolecular structures . J. Comput. Aided Mol. Des. 19, 693–703 (2005).
Lide, D. R. CRC Handbook of Chemistry and Physics 85th edn CRC Press (2004).
Varnek, A. et al. ISIDA—platform for virtual screening based on fragment and pharmacophoric descriptors . Curr. Comput. Aided Drug Des. 4, 191–198 (2008).
Blatov, V. A. Voronoi-Dirichlet polyhedra in crystal chemistry: theory and applications . Crystallogr. Rev. 10, 249–318 (2004).
Carlucci, L., Ciani, G., Proserpio, D. M., Mitina, T. G. & Blatov, V. A. Entangled two-dimensional coordination networks: a general survey . Chem. Rev. 114, 7557–7580 (2014).
Zolotarev, P. N., Arshad, M. N., Asiri, A. M., Al-amshany, Z. M. & Blatov, V. A. A possible route toward expert systems in supramolecular chemistry: 2-periodic h-bond patterns in molecular crystals . Cryst. Growth Des. 14, 1938–1949 (2014).
Cordero, B. et al. Covalent radii revisited . Dalton Trans. 2008, 2832–2838 (2008).
Pauling, L. The Nature of the Chemical Bond and the Structure of Molecules and Crystals: An Introduction to Modern Structural Chemistry Cornell University Press (1960).
Parr, R. G. & Pearson, R. G. Absolute hardness: companion parameter to absolute electronegativity . J. Am. Chem. Soc. 105, 7512–7516 (1983).
Ghosh, D. C. & Biswas, R. Theoretical calculation of absolute radii of atoms and ions. Part 1. The atomic radii . Int. J. Mol. Sci. 3, 87–113 (2002).
Galvez, J., Garcia-Domenech, R., de Julian-Ortiz, J. V. & Soler, R. . Topological approach to drug design . J. Chem. Inf. Comput. Sci. 35, 272–284 (1995).
Kier, L. B. & Hall, L. H. Molecular Structure Description: The Electrotopological State Academic Press (1999).
Friedman, J. H. Greedy function approximation: a gradient boosting machine . Ann. Stat. 29, 1189–1232 (2001).
Freund, Y. & Schapire, R. E. A decision-theoretic generalization of on-line learning and an application to boosting . J. Comput. Syst. Sci. 55, 119–139 (1997).
Loh, W.-Y. Fifty years of classification and regression trees . Int. Stat. Rev. 82, 329–348 (2014).
Schapire, R. E. The strength of weak learnability . Mach. Learn. 5, 197–227 (1990).
Donohue, J. & Lipscomb, W. N. The crystal structure of hydrazinium dichloride, N2H6Cl2 . J. Chem. Phys. 15, 115–119 (1947).
Dulmage, W. J. & Lipscomb, W. N. The crystal structures of hydrogen cyanide, HCN . Acta Cryst. 4, 330–334 (1951).
Kruszynski, R. & Trzesowska, A. Redetermination of hydrogenhydrazinium dichloride . Acta Crystallogr. Sect. E 63, i179 (2007).
Groom, C. R., Bruno, I. J., Lightfoot, M. P. & Ward, S. C. The Cambridge Structural Database . Acta Crystallogr. Sect. B 72, 171–179 (2016).
Lian, C.-S., Wang, X.-Q. & Wang, J.-T. Hydrogenated K4 carbon: a new stable cubic gauche structure of carbon hydride . J. Chem. Phys. 138, 024702 (2013).
Doll, K., Schön, J. C. & Jansen, M. Structure prediction based on ab initio simulated annealing for boron nitride . Phys. Rev. B 78, 144110 (2008).
Escorcia-Salas, G. E., Sierra-Ortega, J. & Rodríguez Martínez, J. A. Influence of Zr concentration on crystalline structure and its electronic properties in the new compound in wurtzite phase: an ab initio study . Microelectr. J. 39, 579–581 (2008).
Li, Q. et al. A novel low compressible and superhard carbon nitride: body-centered tetragonal CN2 . Phys. Chem. Chem. Phys. 14, 13081–13087 (2012).
Marqués, M., Osorio, J., Ahuja, R., Flórez, M. & Recio, J. M. Pressure effects on the structure and vibrations of β- and γ-C3N4 . Phys. Rev. B 70, 104114 (2004).
Hastie, T., Tibshirani, R. & Friedman, J. H. The Elements of Statistical Learning: Data Mining, Inference, and Prediction Springer-Verlag (2001).
Winter, M. Periodic properties: periodicity. WebElements. Available at https://www.webelements.com/periodicity/ (accessed on 8 August 2016).
Minnich, A. J. Phonon heat conduction in layered anisotropic crystals . Phys. Rev. B 91, 085206 (2015).
Shimahara, H. & Kohmoto, M. Anisotropic superconductivity mediated by phonons in layered compounds with weak screening effects . Phys. Rev. B 65, 174502 (2002).
Jha, S. S. Pairing mechanisms and anisotropic superconductivity in layered crystals . Phase Transit. 19, 3–13 (1989).
Klein, J. et al. Superconductivity in high Debye temperature material . Solid State Commun. 33, 1091–1095 (1980).
Figge, S., Kröncke, H., Hommel, D. & Epelbaum, B. M. Temperature dependence of the thermal expansion of AlN . Appl. Phys. Lett. 94, 101915 (2009).
Degtyareva, O., McMahon, M. I. & Nelmes, R. J. Crystal structure of the high pressure phase of bismuth Bi-III . in European Powder Diffraction EPDIC 7 469–475Trans Tech Publications (2001).
Kocak, B., Ciftci, Y. O., Colakoglu, K. & Deligoz, E. Structural, elastic, electronic, and thermodynamic properties of PrN from first principles calculations . Physica B 405, 4139–4144 (2010).
Zwijnenburg, M. A., Corá, F. & Bell, R. G. Isomorphism of anhydrous tetrahedral halides and silicon chalcogenides: energy landscape of crystalline BeF2, BeCl2, SiO2, and SiS2 . J. Am. Chem. Soc. 130, 11082–11087 (2008).
Morelli, D. T. & Slack, G. A. in High Thermal Conductivity Materials (eds Shindé S. L., Goela J. S. Springer (2006).
Madelung O. (ed.) Semiconductors—Basic Data, 2nd edn Springer (1996).
Zhou, Y. & Xiang, H. Al5BO9: a wide band gap, damage-tolerant, and thermal insulating lightweight material for high-temperature applications . J. Am. Ceram. Soc. 99, 2742–2751 (2016).
Curtarolo, S. et al. AFLOW: An automatic framework for high-throughput materials discovery . Comput. Mater. Sci. 58, 218–226 (2012).
Setyawan, W. & Curtarolo, S. High-throughput electronic band structure calculations: challenges and tools . Comput. Mater. Sci. 49, 299–312 (2010).
Jahnátek, M. et al. Ordered phases in ruthenium binary alloys from high-throughput first-principles calculations . Phys. Rev. B 84, 214110 (2011).
Hart, G. L. W., Curtarolo, S., Massalski, T. B. & Levy, O. Comprehensive search for new phases and compounds in binary alloy systems based on platinum-group metals, using a computational first-principles approach . Phys. Rev. X 3, 041035 (2013).
Levy, O., Hart, G. L. W. & Curtarolo, S. Uncovering compounds by synergy of cluster expansion and high-throughput methods . J. Am. Chem. Soc. 132, 4830–4833 (2010).
Perdew, J. P. Density functional theory and the band gap problem . Int. J. Quantum Chem. 28, 497–523 (1985).
Calderon, C. E. et al. The AFLOW standard for high-throughput materials science calculations . Comput. Mater. Sci. 108, Part A 233–238 (2015).
Yazyev, O. V., Kioupakis, E., Moore, J. E. & Louie, S. G. Quasiparticle effects in the bulk and surface-state bands of Bi2Se3 and Bi2Te3 topological insulators . Phys. Rev. B 85, 161101 (2012).
Zheng, X., Cohen, A. J., Mori-Sánchez, P., Hu, X. & Yang, W. Improving band gap prediction in density functional theory from molecules to solids . Phys. Rev. Lett. 107, 026403 (2011).
Blöchl, P. E. Projector augmented-wave method . Phys. Rev. B 50, 17953–17979 (1994).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method . Phys. Rev. B 59, 1758–1775 (1999).
de Jong, M. et al. Charting the complete elastic properties of inorganic crystalline compounds . Sci. Data 2, 150009 (2015).
Blanco, M. A., Francisco, E. & Luaña, V. GIBBS: isothermal-isobaric thermodynamics of solids from energy curves using a quasi-harmonic Debye model . Comput. Phys. Commun. 158, 57–72 (2004).
A.T. and O.I. acknowledge support from DOD-ONR (N00014-13-1-0028 and N00014-16-1-2311) and ITS Research Computing Center at UNC. Development of the web service was supported by the Russian Scientific Foundation (# 14-43-00024). O.I. acknowledges Extreme Science and Engineering Discovery Environment (XSEDE) award DMR-110088, which is supported by National Science Foundation grant number ACI-1053575. S.C. and C.T. acknowledge support from DOD-ONR (N00014-13-1-0030, N00014-13-1-0635), DOE (DE-AC02-05CH11231, specifically BES Grant # EDCBEE) and the Duke University Center for Materials Genomics. C.O. acknowledges support from the National Science Foundation Graduate Research Fellowship under Grant No. DGF-1106401. AFLOW calculations were performed at the Duke University Center for Materials Genomics.
A.T. and S.C. designed the study. O.I. developed and implemented the method. C.O. and C.T. prepared the data and worked with the AFLOW database. E.G. developed the open-access online application available at aflow.org/aflow-ml leveraging the ML models. O.I and C.O. contributed equally to the work. All authors discussed the results and their implications and contributed to the writing of the article.
Correspondence to Olexandr Isayev or Stefano Curtarolo.

matches: V. V. V. V. V. V. V. V. V.


text:
Fort Lauderdale, Florida Business and Corporate, Taxation, & Estate Planning Attorney |Larry V. Bishins, P.A. - Larry V. Bishins, P.A.
The world of business can be as dangerous as it is exciting. If you are planning to start a business venture or want to expand your current business, it is important to know how to protect yourself and your assets. For the last 41 years, Larry V. Bishins has been dedicated to tax planning for business. When companies in Florida need a business law attorney to help grow their businesses and protect their best interests, they turn to Larry V. Bishins to get results. Larry V. Bishins’ philosphy is that it not how much you make, but how much you keep that counts.
Larry V. Bishins firmly stands beside his clients from the formation of a business to its dissolution. Larry V. Bishins reviews clients’ business plans and advises on all matters involved in properly forming and maintaining a business, addressing topics such as potential liability issues, corporate entity formation and tax structures.
Larry V. Bishins drafts contracts and buy-sell agreements on behalf of his clients and protects their interests in shareholder disputes. You can be confident that Larry V. Bishins is a seasoned negotiator, who negotiates contracts for his clients.
Larry V. Bishins stands ready to protect your business and help turn it into a success. Call the firm’s Fort Lauderdale office at 954-772-7900 or contact us online to schedule your initial consultation.

matches: V. V. V. V. V. V. V. V. V. V.


text:
The Complainants are Credit Industriel et Commercial S.A. and Conf�d�ration Nationale du Credit Mutuel, both of France and both represented by Meyer & Partenaires, Strasbourg, France.
The Respondent is Spiral Matrix, Eldoret, Kenya.
The disputed domain names <cicenterprise.com>, <creditmutuelocean.com> and <filban.com> are registered with Intercosmos Media Group d/b/a directNIC.com.
The Complaint was filed with the WIPO Arbitration and Mediation Center (the “Center”) on March 3, 2006. On March 6, 2006, the Center transmitted by email to Intercosmos Media Group d/b/a directNIC.com a request for registrar verification in connection with the domain names at issue. On the same date, Intercosmos Media Group d/b/a directNIC.com transmitted by email to the Center its verification response confirming that the Respondent is listed as the registrant and providing the contact details for the administrative and technical contact. The Center verified that the Complaint satisfied the formal requirements of the Uniform Domain Name Dispute Resolution Policy (the “Policy”), the Rules for Uniform Domain Name Dispute Resolution Policy (the “Rules”), and the WIPO Supplemental Rules for Uniform Domain Name Dispute Resolution Policy (the “Supplemental Rules”).
In accordance with the Rules, paragraphs 2(a) and 4(a), the Center formally notified the Respondent of the Complaint, and the proceedings commenced on March 14, 2006. In accordance with the Rules, paragraph 5(a), the due date for Response was April�3,�2006. The Respondent did not submit any response. Accordingly, the Center notified the Respondent’s default on April 7, 2006.
The Center appointed Adam Samuel as the sole panelist in this matter on April�19,�2006. The Panel finds that it was properly constituted. The Panel has submitted the Statement of Acceptance and Declaration of Impartiality and Independence, as required by the Center to ensure compliance with the Rules, paragraph 7.
The three domain names in dispute were registered on July 21, 2005, (<cicenterprise.com>), August 6, 2005 (<creditmutuelocean.com>), and August�12,�2005 (<filban.com>).
The first Complainant is a French banking group with branches in a number of different countries around the world. It uses CIC as a common acronym and trade name. It offers Internet facilities through its online banking service, named “Filbanque” which it has offered since 1998, through the domain name, <filbanque.com>.
The first Complainant registered the trademark C.I.C. in 1976, in France, renewing it in 1986 and 1996. It has registered the trademark CIC BANQUES in 1991, renewing it in 2001. It has a trademark registration for FILBANQUE that it last renewed in 2002.
The second Complainant offers banking and insurance services in France. Its domain name <creditmutuel.com> describes those services. Its western France subsidiary, Credit Mutuel Ocean has been registrant of the domain name <cmocean.fr> since 1998.
The second Complainant owns trademark registrations for CREDIT MUTUAL dating back to 1988. One of the second Complainant’s branches, F�d�ration du Credit Mutuel Ocean has owned a trademark registration for CREDIT MUTUEL OCEAN CMO since July�1999.
In 1998, the Complainants merged to form the Group Credit Mutuel-CIC.
On October 26, 2005, the Complainants’ trademark attorneys wrote to the Respondent asking it to transfer the domain names in dispute to the Complainants. They did not receive a reply to this letter and e-mail.
The domain name, <cicenterprise.com> is confusingly similar to the first Complainant’s trademark CIC, which is reproduced in the domain name with the generic word “enterprise” added to it. The mere addition of a generic or descriptive word to a trademark does not prevent a domain name from being confusingly similar to the trademark. The domain name, <cicenterprise.com> provides hyperlinks to banking and financial services providers.
The domain name, <filban.com> is confusingly similar to the first Complainant’s trademark FILBANQUE reproducing as it does the first six letters of the trademark. Again, the domain name provides hyperlinks to banking and financial services providers.
Filban has no independent meaning in English or French. Considering the other domain names that seem to refer to the two other members of the same banking group, one must assume that the domain name, <filban.com> is aimed at diverting traffic from people seeking the “Filbanque” service of the first Complainant.
The second Complainant claims that this domain name is confusingly similar to its trademark CREDIT MUTUEL, which is reproduced at the start of the domain name. The addition of the word “ocean” to the trademark adds to the confusion because it identifies the domain name with the second Complainant’s western France subsidiary’s name.
As regards <creditmutuelocean.com>, the Respondent could not have selected the second Complainant’s subsidiary’s name randomly. It must have been done to attract users seeking the website associated with the subsidiary concerned.
The Complainants argue that the Respondent has no rights or legitimate interests in any of the domain names in question. The Respondent has no connection with the Complainants’ business. The Respondent is not known nor has it ever been known by the names “CIC”, “CICenterprise”, “Credit Mutuel”, “Credit Mutuel Ocean”, “Filbanque” or “Filban”. The Complainants have not licensed the Respondent to use or apply for registration of the domain names in issue.
The domain names have been registered to take advantage of the Complainants’ renown in order to divert Internet users to competing websites through a “pay-per-click” domain parking solution. The Respondent must have known of the Complainants’ trademarks since the company’s names are well known.
The Complainants further contend that many of the hyperlinks on the domain names in question relate to banking products or services indicating clearly that the Respondent was aware of the Complainants’ activities. It can be deduced from this that the names were selected to divert the Complainants’ customers to the Respondent’s website and through it to companies that would pay it by pay-for-clicks. The Respondent has registered the disputed domain names only for the purpose of taking advantage of Internet traffic from Internet users seeking the Complainants’ products and services on the Internet.
The Respondent also has a lengthy history of cybersquatting as indicated by UDRP domain name cases reported in the past.
The Complainants claim that the Respondent uses the domain names in issue to divert Internet users to parking webpages on which several hyperlinks are displayed, some of which lead to hyperlinks in French covering banking and financial services. By clicking on the hyperlinks, users are led to websites of businesses that compete with the Complainants.
The domain name, <cicenterprise.com> is confusingly similar to the first Complainant’s trademark CIC, which is reproduced in the domain name with the generic word “enterprise” added to it. The mere addition of a generic or descriptive word to a trademark name does not prevent a domain name from being confusingly similar to the trademark.
“There can be no serious dispute that the prefix “BARCLAYCARDBUSINESS” is confusingly similar to the word mark “BARCLAYCARD” which the Complainant has registered”.
The domain name, <filban.com> is confusingly similar to the first Complainant’s trademark FILBANQUE reproducing as it does the first six letters of the trademark. The word “filban” has no meaning in French or English. In auditory terms, there is only one consonant sound missing from the domain name for it to reproduce the trademark concerned.
In that case, Volvo was associated with an ocean race. Here, the second Complainant has a subsidiary with the name “ocean” added to its name.
The Respondent is not called any of the domain names in question or by a similar name and does not appear to trade under those names. There is no evidence that the owners of the trademarks in question, the Complainants, have authorized it to use the trademarks concerned. The Respondent has never asserted any rights or legitimate interests in the name. For these reasons and the on the basis of the case file, notably the absence of a Response, the Panel concludes that the Respondent has no rights or legitimate interests in respect of the domain names.
The domain name <cicenterprise.com> currently links to a website promoting mainly holidays but also some financial services products and services. The website gives no indication that the site has no connection with the Complainant.
The domain name <filban.com> currently links to a website promoting a variety of lifestyle products and services including those relating to financial services.
The domain name, <creditmutuelocean.com> currently links to a website promoting mainly financial services and banking-related products and services although it also contains links to lifestyle and travel sites.
Of the three trademarks in question, two are reasonably well known even outside France. “Filban” has no independent meaning and must be aimed at people who mistype “filbanque”. It seems highly likely that the Respondent knew that it was registering the domain names in issue in probable breach of a trademark. At best, it closed its eyes to the possibility.
The Respondent is using the Complainants’ trademarks to attract visitors to a website to profit, presumably through “click-throughs” to a number of other websites. It is irrelevant to the question of bad faith as to whether this is likely to be successful.
The Respondent’s motive in registering and using the sites seems to be either to attract Internet users for commercial gain or simply to disrupt the Complainants’ relationship with its customers or potential customers. Both these motives constitute evidence of registration and use in bad faith: paragraphs 4(b)(iii) and (iv) of the Policy.
As the Complainants point out, the Respondent has an extensive record of being involved in cybersquatting. In the following WIPO UDRP cases, domain name panels have made findings of bad faith against the Respondent: The Sutton Corporation v. Spiral Matrix, WIPO Case No. D2006-0167; PRL USA Holdings, Inc. v. Spiral Matrix, WIPO Case No. D2006-0189; Soci�t� des H�tels Meridien v. Spiral Matrix / Kentech Inc., WIPO Case No.�D2005-1196; Finaxa S.A. v. Spiral Matrix, WIPO Case No. D2005-1044; Dr. Ing. h.c. F. Porsche AG v. Kentech, Inc. a.k.a. Helois Lab a.k.a. Orion Web a.k.a. Titan Net a.k.a. Panda Ventures a.k.a. Spiral Matrix and Domain Purchase, NOLDC, Inc., WIPO Case No. D2005-0890; Staples, Inc., Staples The Office Superstore, LLC, and Staples Contract & Commercial,�Inc. v. Spiral Matrix a/k/a Kenyatech.com a/k/a Kentech.com, WIPO Case No.�D2006-0227; Consitex S.A., Lanificio Ermenegildo Zegna & Figli S.p.A., and Ermenegildo Zegna Corporation v. Spiral Matrix, WIPO Case No. D2006-0169; Deutsche Telekom AG v. Spiral Matrix, WIPO Case No.�D2005-1145.
For all these reasons, the Panel concludes that the domain names in dispute have been registered and are being used in bad faith.
For all the foregoing reasons, in accordance with paragraphs 4(i) of the Policy and 15 of the Rules, the Panel orders that the domain names, <cicenterprise.com>, <creditmutuelocean.com> and <filban.com> be transferred to the Complainant.

matches: v. v. v. v. v. v. v. v.


U.S. Title 50 -- War and National Defense, Chapter 15 -- National Security contains, in one place, a collation of most of the law as it pertains to the Central Intelligence Agency. Most people who write about the CIA and who talk about the CIA -- indeed, many who have served with the CIA -- have never read this law. It is most significant that the legislation that pertains to war and national defense is the same legislation that includes all reference to the CIA. It is almost as if the bomb contained its own live fuse or the gun came with the trigger cocked for action. As we have seen, during the past twenty-five years the CIA has become the active agent that ignites the military establishment whenever that great mass becomes supercritical.
§ 401. Congressional declaration of purpose.
§ 402. National Security Council.
(a) Establishment, presiding officer; functions, composition.
There is established a council to be known as the National Security Council (hereinafter in this section referred to as the "Council").
The President of the United States shall preside over meetings of the Council: Provided, That in his absence he may designate a member of the Council to preside in his place.
The function of the Council shall be to advise the President with respect to the integration of domestic, foreign, and military policies relating to the national security so as to enable the military services and the other departments and agencies of the Government to cooperate more effectively in matters involving the national security.
(7) the Secretaries and Under Secretaries of other executive departments and of the military departments, the Chairman of the Munitions Board, and the Chairman of the Research and Development Board, when appointed by the President by and with the advice and consent of the Senate, to serve at his pleasure.
(2) to consider policies on matters of common interest to the departments and agencies of the Government concerned with the national security, and to make recommendations to the President in connection therewith.
(c) Executive secretary; appointment and compensation; staff employees.
The Council shall have a staff to be headed by a civilian executive secretary who shall be appointed by the President. The executive secretary, subject to the direction of the Council, is authorized, subject to the civil-service laws and chapter 51 and subchapter III of chapter 53 of title %, to appoint and fix the compensation of such personnel as may be necessary to perform such duties as may be prescribed by the Council in connection with the performance of its functions.
1949 -- Subsec. (a) Act Aug. 10, 1949, added the Vice President to the Council, removed the Secretaries of the military departments, to authorize the President to add, with the consent of the Senate, Secretaries and Under Secretaries of other executive departments, to authorize the President to add, with the consent of the Senate, Secretaries and Under Secretaries of other executive departments and of the military department, and the Chairmen of the Munitions Board and the Research and Development Board.
The National Security Council, together with its functions, records, property, personnel, and unexpended balances of appropriations, allocations, and other funds (available or to be made available) were transferred to the Executive Office of the President by 1949 Reorg. Plan No. 4, eff. Aug. 19, 1949, 14 F.R. 5227, 63 Stat. 1067, set out in the Appendix to Title 5, Government Organization and Employees.
§ 403. Central Intelligence Agency.
(a) Establishment; Director and Deputy Director; appointment.
There is established under the National Security Council a Central Intelligence Agency with a Director of Central Intelligence who shall be the head thereof, and with a Deputy Director of Central Intelligence who shall act for, and exercise the powers of, the Director during his absence or disability. The Director and the Deputy Director shall be appointed by the President, by and with the advice and consent of the Senate, from among the commissioned officers of the armed services, whether in an active or retired status, or from among individuals in civilian life: Provided, however, That at no time shall the two positions of the Director and Deputy Director be occupied simultaneously by commissioned officers of the armed services, whether in active or retired status.
(b) Commissioned officer as Director or Deputy Director; powers and limitations, effect on commissioned status.
(B) he shall not possess or exercise any supervision, control, powers, or functions (other than such as he possesses, or is authorized or directed to exercise, as Director, or Deputy Director) with respect to the armed services or any component thereof, the Department of the Army, the Department of the Navy, or the Department of the Air Force, or any branch, bureau, unit, or division thereof, or with respect to any of the personnel (military or civilian) of any of the foregoing.
(2) Except as provided in paragraph (1) of this subsection, the appointment to the office of Director, or Deputy Director, of a commissioned officer of the armed services, and his acceptance of and service in such office, shall in no way affect any status, office, rank, or grade he may occupy or hold in the armed services, or any emolument, perquisite, right, privilege, or benefit incident to or arising out of any such status, office, rank, or grade. Any such commissioned officer shall, while serving in the office of Director, or Deputy Director, continue to hold rank and grade not lower than that in which serving at the time of his appointment and to receive the military pay and allowances (active or retired, as the case may be, including personal money allowance) payable to a commissioned officer of his grade and length of service for which the appropriate department shall be reimbursed from any funds available to defray the expenses of the Central Intelligence Agency. He also shall be paid by the Central Intelligence Agency from such funds an annual compensation at a rate equal to the amount by which the compensation established for such position exceeds the amount of his annual military pay and allowances.
(3) The rank or grade of any such commissioned officer shall, during the period in which such commissioned officer occupies the office of Director of Central Intelligence, or Deputy Director of Central Intelligence, be in addition to the numbers and percentages otherwise authorized and appropriated for the armed service of which he is a member.
(c) Termination of employment of officers and employees; effect on right of subsequent employment.
Notwithstanding the provisions of section 652 of Title 5, or the provisions of any other law, the Director of Central Intelligence may, in his discretion, terminate the employment of any officer or employee of the Agency whenever he shall deem such termination necessary or advisable in the interests of the United States, but such termination shall not affect the right of such officer or employee to seek or accept employment in any other department or agency of the Government if declared eligible for such employment by the United States Civil Service Commission.
(5) to perform such other functions and duties related to intelligence affecting the national security as the National Security Council may from time to time direct.
(e) Inspection of intelligence of other departments.
To the extent recommended by the National Security Council and approved by the President, such intelligence of the departments and agencies of the Government, except as hereinafter provided, relating to the national security shall be open to the inspection of the Director of Central Intelligence, and such intelligence as relates to the national security and is possessed by such departments and other agencies of the Government, except as hereinafter provided, shall be made available to the Director of Central Intelligence for correlation, evaluation, and dissemination: Provided, however, That upon the written request of the Director of Central Intelligence, the Director of the Federal Bureau of Investigation shall make available to the Director of Central Intelligence such information for correlation, evaluation, and dissemination as may be essential to the national security.
(f) Termination of National Intelligence Authority; transfer of personnel, property, records, and unexpended funds.
(2) the personnel, property, and records of the Central Intelligence Group are transferred to the Central Intelligence Agency, and such Group shall cease to exist. Any unexpended balances of appropriations, allocations, or other funds available or authorized to be made available for such Group shall be available and shall be authorized to be made available in like manner for expenditure by the Agency.
(4) report to the President concerning the Boards findings and appraisals, and make appropriate recommendations for actions to achieve increased effectiveness of the Governments foreign intelligence effort in meeting national intelligence needs.
SEC. 2. In order to facilitate performance of the Boards functions, the Director of Central Intelligence and the heads of all other departments and agencies shall make available to the Board all information with respect to foreign intelligence and related matters which the Board may require for the purpose of carrying out its responsibilities to the President in accordance with the terms of this Order. Such information made available to the Board shall be given all necessary security protection in accordance with the terms and provisions of applicable laws and regulations.
SEC. 3. Members of the Board shall be appointed by the President from among persons outside the Government, qualified on the basis of knowledge and experience in matters relating to the national defense and security, or possessing other knowledge and abilities which may be expected to contribute to the effective performance of the Boards duties. The members of the Board shall receive such compensation and allowances, consonant with law, as may be prescribed hereafter.
SEC. 4. The Board shall have a staff headed by an Executive Secretary, who shall be appointed by the President and shall receive such compensation and allowances, consonant with law, as may be prescribed by the Board. The Executive Secretary shall be authorized, subject to the approval of the Board and consonant with law, to appoint and fix the compensation of such personnel as may be necessary for performance of the Boards duties.
SEC. 6. Executive Order No. 10938 of May 4, 1961, is hereby revoked.
Act June 20, 1949, § 10, formerly § 12, 63 Stat. 212, renumbered July 7, 1958, Pub. L. 85-507, § 21(b) (2), 72 Stat. 337, provided that Act June 20, 1949, which is classified to sections 403a-403j of this title, should be popularly known as the "Central Intelligence Agency Act of 1949".
§ 403c. Same; procurement authority.
(a) In the performance of its functions the Central Intelligence Agency is authorized to exercise the authorities contained in sections 151 (c) (1)-(6), (10), (1), (15), (17), 155 and 159 of Title 41.
(b) In the exercise of the authorities granted in subsection (a) of this section, the term "Agency head" shall mean the Director, the Deputy Director, or the Executive of the Agency.
(c) The determinations and decisions provided in subsection (a) of this section to be made by the Agency head may be made with respect to individual purchases and contracts or with respect to classes of purchases or contracts, and shall be final. Except as provided in subsection (d) of this section, the Agency head is authorized to delegate his powers provided in this section, including the making of such determinations and decisions, in his discretion and subject to his direction to any other officer or officers or officials of the Agency.
§ 403f. Same; general authorities of Agency.
The act of June 30, 1932, as amended, referred to in subsec. (c), is the Legislative Branch Appropriation Act, 1933, act June 30, 1932, ch. 314, 47 Stat. 393, and is classified to section 303b of Title 40, Public Buildings, Property, and Works.
Section was not enacted as a part of the National Security Act of 1947 which comprises this chapter.
§ 403g. Same; protection of nature of Agencys functions.
§ 403h. Same; admission of essential aliens; limitation on number.
§ 403j. Central Intelligence Agency: appropriations; expenditures.
(2) supplies, equipment, and personnel and contractual services otherwise authorized by law and regulations, when approved by the Director.
Sections 259 and 267 of Title 40, referred to in text, was repealed by Pub. L. 86-249, § 17 (12), Sept. 9, 1959, 73 Stat. 485. See chapter 12 of Title 40, Public Buildings, Property and Works.
§ 405 Advisory Committees; appointment; compensation of part-time personnel; applicability of other laws.
(a) The Director of the Office of Emergency Preparedness, the Director of Central Intelligence, and the National Security Council, acting through its Executive Secretary, are authorized to appoint such advisory committees and to employ, consistent with other provisions of this Act, such part-time advisory personnel as they may deem necessary in carrying out their respective functions and the functions of agencies under their control. Persons holding other offices or positions under the United States for which they receive compensation, while serving as members of such committees, shall receive no additional compensation for such service. Other members of such committees and other part-time advisory personnel so employed may serve without compensation or may receive compensation at a rate not to exceed $50 for each day of service, as determined by the appointing authority.
§ 407. Study or plan of surrender; use of appropriations.
Section was not enacted as part of the National Security Act of 1947, which comprises this chapter.
§ 409. Definitions of military departments.
(a) The term "Department of the Army" as used in this Act shall be construed to mean the Department of the Army at the seat of the government and all field headquarters, forces, reserve components, installations, activities, and functions under the control or supervision of the Department of the Army.
(b) The term "Department of the Navy" as used in this Act shall be construed to mean the Department of the Navy at the seat of the government; the headquarters, United States Marine Corps; the entire operating forces of the united States Navy, including naval aviation, and of the United States Marine Corps, including the reserve components of such forces; all field activities, headquarters, forces, bases, installations, activities, and functions under the control or supervision of the Department of the Navy; and the United States Coast Guard when operating as a part of the Navy pursuant to law.


text:
Vibration Analysis Based on HJ-Biplots. International Journal of Prognostics and Health Management, v. 9 (2), n. 030, p. 01-10, 2018.
Leveraging Intellectual Capital through Lewin's Force Field Analysis: The case of Software Development Companies. Journal of Innovation & Knowledge (JIK), v. 2, n. 3, p. 125-133, 2017.
Knowledge management as a factor for the formulation and implementation of organization strategy. Journal of Knowledge Management (JKM), v. 21, n. 2, p. 308-329, 2017.
Special Issue: Intellectual capital management as a driver of competitiveness and sustainability. Journal of Intellectual Capital (JIC), v. 18, n. 3, 2017.
Intellectual Capital Management as a Driver of Competitiveness and Sustainability (Guest Editorial). Special Issue: Intellectual capital management as a driver of competitiveness and sustainability. Journal of Intellectual Capital (JIC), v. 18, n. 3, p. 466-469, 2017.
Methods of evaluation of intangible assets and intellectual capital. Special Issue: Intellectual capital management as a driver of competitiveness and sustainability. Journal of Intellectual Capital (JIC), v. 18, n. 3, p. 470-485, 2017.
Increasing smart city competitiveness and sustainability through managing structural capital. Special Issue: Intellectual capital management as a driver of competitiveness and sustainability. Journal of Intellectual Capital (JIC), v. 18, n. 3, p. 693-707, 2017.
Felicidade e Redes nas Empresas: Aplicação no Norte de Portugal e Galiza (Espanha). International Journal of Knowledge Engineering and Management (IJKEM), v. 5, n. 11, p. 28-50, 2016.
Os Processos de Gestão do Conhecimento Influentes na Inovação Organizacional. International Journal of Knowledge Engineering and Management (IJKEM), v. 5, n. 12, p. 1-28, 2016.
Knowledge Management in Brazilian, Portuguese and Polish Organizations: a Comparative Analysis. Online Journal of Applied Knowledge Management v. 4, n. 1, p. 29 – 41, 2016.
A Theoretical Model for the Report of Intellectual Capital. Electronic Journal of Knowledge Management, v. 11, n. 4, p. 339-360, 2013.
Biplot Methodology Applied to an Intellectual Capital Model. Electronic Journal of Knowledge Management, v. 11, n. 1, p. 81-92, 2013.
Linking Knowledge Management and Intellectual Capital: A Study in Portuguese Banks. International Journal of Sociotechnology and Knowledge Development, v. 5, n. 4, p. 1-13.
Why Intellectual Capital Management Accreditation Is a Tool for Organizational Development? Electronic Journal of Knowledge Management, v. 8, n. 2, p. 235-244, 2012.

matches: v. v. v. v. v. v. v. v. v. v. v. v. v. v.


Paul J. TRAFFICANTE et al., Petitioners, v. METROPOLITAN LIFE INSURANCE COMPANY et al.
Two tenants of an apartment complex filed complaints with the Secretary of Housing and Urban Development alleging that their landlord racially discriminated against nonwhites, that the tenants thereby lost the social benefits of living in an integrated community, missed business and professional advantages that would have accrued from living with members of minority groups, and suffered from being 'stigmatized' as residents of a 'white ghetto.' The District Court, not reaching the merits, held that the complaining tenants were not within the class of persons entitled to sue under § 810(a) of the Civil Rights Act of 1968. The Court of Appeals, in affirming, construed § 810(a) to permit complaints only by persons who are the objects of discriminatory housing practices. Held: The definition in § 810(a) of 'person aggrieved,' as 'any person who claims to have been injured by a discriminatory housing practice,' shows a congressional intention to define standing as broadly as is permitted by Article III of the Constitution, and petitioners, being tenants of the apartment complex, have standing to sue under § 810(a). Pp. 208212.
Two tenants of Parkmerced, an apartment complex in San Francisco housing about 8,200 residents, filed separate complaints with the Secretary of Housing and Urban Development (HUD) pursuant to § 810(a) 1 of the Civil Rights Act of 1968, 82 Stat. 85, 42 U.S.C. 3610(a). One tenant is black, one white. Each alleged that the owner 2 of Parkmerced had discriminated against nonwhites on the basis of race in the rental of apartments within the complex in violation of § 804 of the Act.
The Act gives the Secretary of HUD power to receive and investigate complaints regarding discriminatory housing practices. The Secretary, however, must defer to state agencies that can provide relief against the named practice. If the state agency does not act, the Secretary may seek to resolve the controversy by conference, conciliation, or persuasion. If these attempts fail, the complainant may proceed to court pursuant to § 810(d). 7 Moreover, these rights may be enforced 'by civil actions in appropriate United States district courts without regard to the amount in controversy,' if brought within 180 days 'after the alleged discriminatory housing practice occurred.' § 812(a). In addition, § 813 gives the Attorney General authority to bring a civil action in any appropriate United States district court when he has reasonable cause to believe 'that any person or group of persons is engaged in a pattern or practice of resistance to the full enjoyment of any of the rights granted' by the Act.
It is apparent, as the Solicitor General says, that complaints by private persons are the primary method of obtaining compliance with the Act. Hackett v. McGuire Bros., Inc., 445 F.2d 442 (CA 3), which dealt with the phrase that allowed a suit to be started 'by a person claiming to be aggrieved' under the Civil Rights Act of 1964, 42 U.S.C. 2000e5(a), concluded that the words used showed 'a congressional intention to define standing as broadly as is permitted by Article III of the Constitution.' Id., at 446. With respect to suits brought under the 1968 Act, 8 we reach the same conclusion, insofar as tenants of the same housing unit that is charged with discrimination are concerned.
The Assistant Regional Administrator for HUD wrote petitioners' counsel on November 5, 1970, that 'it is the determination of this office that the complainants are aggrieved persons and as such are within the jurisdiction' of the Act. We are told that that is the consistent administrative construction of the Act. Such construction is entitled to great weight. Udall v. Tallman, 380 U.S. 1, 16, 85 S.Ct. 792, 801, 13 L.Ed. 616; Griggs v. Duke Power Co., 401 U.S. 424, 433434, 91 S.Ct. 849, 854 855, 28 L.Ed.2d 158.
Absent the Civil Rights Act of 1968, I would have great difficulty in concluding that petitioners' complaint in this case presented a case or controversy within the jurisdiction of the District Court under Art. III of the Constitution. But with that statute purporting to give all those who are authorized to complain to the agency the right also to sue in court, I would sustain the statute insofar as it extends standing to those in the position of the petitioners in this case. Cf. Katzenbach v. Morgan, 384 U.S. 641, 648649, 86 S.Ct. 1717, 1722, 16 L.Ed.2d 828 (1966); Oregon v. Mitchell, 400 U.S. 112, 240, 248249, 91 S.Ct. 260, 322, 326327, 27 L.Ed.2d 272 (1970). Consequently, I join the Court's opinion and judgment.
We find it unnecessary to reach the question of standing to sue under 42 U.S.C. 1982 which is the basis of the third cause of action alleged in the petition but based on the same allegations as those made under the Civil Rights Act of 1968.
The Dirksen substitute, 114 Cong.Rec. 45704573 retained the present language of § 810(a) which Senator Mondale had previously introduced, id., at 2270, and it was in the bill passed by the Senate, id., at 5992, which the House subsequently passed, id., at 9621.