Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

147

148

149

150

151

152

153

154

155

156

157

158

159

160

161

162

163

164

165

166

167

168

169

170

171

172

173

174

175

176

177

178

179

180

181

182

183

184

185

""" 

This module introduces the ForceConstantPotential object which acts as the 

finalized force constant model. 

""" 

 

import pickle 

import numpy as np 

from collections import Counter 

 

from .force_constant_model import ForceConstantModel 

from .core.orbits import Orbit 

from .core.orbits import OrientationFamily 

from .core.tensors import rotation_to_cart_coord, rotate_tensor 

 

 

# TODO: Fix the relation with cluster space 

class ForceConstantPotential: 

""" A finalized force constant model. Can produce force constants for any 

structure compatible with the structure for which the model was set up. 

 

Parameters 

---------- 

cs : ClusterSpace 

The cluster space the model is based upon 

parameters : numpy.ndarray 

The fitted paramteres 

""" 

 

def __init__(self, cs, parameters): 

 

self._parameters = parameters.copy() 

self._prim = cs.primitive_structure.copy() 

self.cluster_list = cs.cluster_list.copy() 

self.atom_list = cs.atom_list.copy() 

self.orbits = [] 

self.spacegroup = cs.spacegroup 

self._config = cs._config 

 

# Extract the eigentensors from the cluster space and use the paramters 

# to construct the finalized force constants 

parameters = cs._map_parameters(parameters) 

p = 0 

for orb in cs.orbits: 

new_orbit = Orbit() 

fc = np.zeros(orb.eigentensors[0].shape) 

for et, a in zip(orb.eigentensors, parameters[p:]): 

fc += et * a 

new_orbit.force_constant = fc 

new_orbit.order = orb.order 

new_orbit.radius = orb.radius 

new_orbit.maximum_distance = orb.maximum_distance 

for of in orb.orientation_families: 

new_of = OrientationFamily() 

new_of.cluster_indices = of.cluster_indices.copy() 

sym_ind = of.symmetry_index 

R = rotation_to_cart_coord(cs.rotation_matrices[sym_ind], 

self.primitive_structure.cell) 

fc = rotate_tensor(new_orbit.force_constant, R.T) 

perm = cs.permutations[of.permutation_indices[0]] 

new_of.force_constant = fc.transpose(perm) 

new_orbit.orientation_families.append(new_of) 

self.orbits.append(new_orbit) 

p += len(orb.eigentensors) 

 

@property 

def parameters(self): 

return self._parameters.copy() 

 

@property 

def symprec(self): 

return self._config['symprec'] 

 

@staticmethod 

def read(f): 

"""Reads a force constant potential from file. 

 

Parameters 

---------- 

f : str or file object 

name of input file (str) or stream to load from (file object) 

 

Returns 

------- 

ForceConstantPotential 

the original object as stored in the file 

""" 

if isinstance(f, str): 

with open(f, 'rb') as fobj: 

return pickle.load(fobj) 

else: 

try: 

return pickle.load(f) 

except Exception: 

raise Exception('Failed loading from file.') 

 

def write(self, f): 

"""Writes a force constant potential to file. 

 

Parameters 

---------- 

f : str or file object 

name of input file (str) or stream to write to (file object) 

""" 

if isinstance(f, str): 

with open(f, 'wb') as fobj: 

pickle.dump(self, fobj) 

else: 

try: 

pickle.dump(self, f) 

except Exception: 

raise Exception('Failed writing to file.') 

 

@property 

def primitive_structure(self): 

""" ase.Atoms : atomic structure """ 

return self._prim.copy() 

 

@property 

def orbit_data(self): 

""" list : list of dictionaries containing detailed information for 

each orbit, e.g. cluster radius and force constant 

""" 

data = [] 

for orbit_index, orbit in enumerate(self.orbits): 

d = {} 

d['index'] = orbit_index 

d['order'] = orbit.order 

d['radius'] = orbit.radius 

d['maximum_distance'] = orbit.maximum_distance 

d['n_clusters'] = len(orbit.orientation_families) 

 

types = [] 

for atom_ind in self.cluster_list[orbit.prototype_index]: 

types.append(self.primitive_structure.numbers[ 

self.atom_list[atom_ind].site]) 

d['prototype_cluster'] = self.cluster_list[orbit.prototype_index] 

d['prototype_atom_types'] = types 

 

d['geometrical_order'] = len(set(d['prototype_cluster'])) 

d['force_constant'] = orbit.force_constant 

d['force_constant_norm'] = np.linalg.norm(orbit.force_constant) 

data.append(d) 

return data 

 

def get_force_constants(self, atoms): 

""" Return the force constants of a compatible structure. 

 

Parameters 

---------- 

atoms : ase.Atoms 

input structure 

 

Returns 

------- 

ForceConstants 

force constants 

""" 

return ForceConstantModel(atoms, self).get_force_constants() 

 

def __str__(self): 

orbits = self.orbit_data 

orbit_counts = Counter([orbit['order'] for orbit in orbits]) 

cluster_counts = Counter() 

for orbit in orbits: 

cluster_counts[orbit['order']] += orbit['n_clusters'] 

 

n = 54 

s = [] 

s.append(' ForceConstantPotential '.center(n, '=')) 

s.append('Spacegroup {}'.format(self.spacegroup)) 

s.append('Cell:\n{}'.format(self.primitive_structure.cell)) 

s.append('Basis:\n{}'.format(self.primitive_structure.basis)) 

s.append('Numbers: {}'.format(self.primitive_structure.numbers)) 

for order in sorted(orbit_counts.keys()): 

s.append('Order {}, #orbits {}, #cluster {}'.format( 

order, orbit_counts[order], cluster_counts[order])) 

s.append('Total number of orbits: {} '.format(len(orbits))) 

s.append('total number of clusters: {} ' 

.format(sum(cluster_counts.values()))) 

s.append(''.center(n, '=')) 

return '\n'.join(s) 

 

def __repr__(self): 

return 'ForceConstantPotential(ClusterSpace({!r}, ...), [...])'.format( 

self.primitive_structure)