2d:null 2e:["$","title",null,{"children":"QH | Vinícius Mello | Observable"}] 33:T62e,{ const tileWidth = 30; const svg = d3.select(DOM.svg(3 * 8 * tileWidth, 7 * 8 * tileWidth)); const n = ['c', 'd', 'e', 'f', 'g', 'a', 'b']; const data = []; const data2 = []; for (let k = 0; k < 21; ++k) { data2.push({ kj: k % 3, ki: Math.floor(k / 3), text: 'energy: ' + d3.format('.3f')(orbitals[k].energy) }); for (let i = 0; i < 7; ++i) { for (let j = 0; j < 7; ++j) { data.push({ kj: k % 3, ki: Math.floor(k / 3), i: i, j: j, value: orbitals[k].orbital[i][j], text: n[i] + n[j] }); } } } const tile = svg .selectAll("g") .data(data) .enter() .append("g") .attr("transform", function(d, i) { return `translate(${d.kj * 8 * tileWidth + d.j * tileWidth},${d.ki * 8 * tileWidth + (6 - d.i) * tileWidth})`; }); tile .append("rect") .attr("x", 0) .attr("y", 0) .attr("width", tileWidth) .attr("height", tileWidth) .attr("fill", d => d3.interpolateReds(Math.abs(d.value))) .attr("stroke", "white") .attr("stroke-width", 2); tile .append("text") .attr('x', 5) .attr('y', (tileWidth * 2) / 3) .text(d => d.text); svg .selectAll('.orb') .data(data2) .enter() .append('g') .attr('class', 'orb') .attr('transform', function(d, i) { return `translate(${d.kj * 8 * tileWidth},${d.ki * 8 * tileWidth})`; }) .append("text") .attr('x', 10) .attr('y', 8 * tileWidth - 10) .text(d => d.text); return svg.node(); }34:T7ee,{ const tileWidth = 30; const svg = d3.select(DOM.svg(7 * tileWidth, 7 * tileWidth)); const n = ['c', 'd', 'e', 'f', 'g', 'a', 'b']; const data = []; for (let i = 0; i < 7; ++i) { for (let j = 0; j < 7; ++j) { data.push({ i: i, j: j, value: 0, text: n[i] + n[j] }); } } const tile = svg .selectAll("g") .data(data) .enter() .append("g") .attr("transform", function(d, i) { return `translate(${d.j * tileWidth},${(6 - d.i) * tileWidth})`; }); tile .append("rect") .attr("x", 0) .attr("y", 0) .attr("width", tileWidth) .attr("height", tileWidth) .attr("fill", d => d3.interpolateReds(Math.abs(d.value))) .attr("stroke", "white") .attr("stroke-width", 2); tile .append("text") .attr('x', 5) .attr('y', (tileWidth * 2) / 3) .text(d => d.text); const A = orbitals.map(o => o.orbital[0][2]); while (true) { const re = d3.range(7).map(i => d3.range(7).map(j => 0)); const im = d3.range(7).map(i => d3.range(7).map(j => 0)); const de = d3.range(7).map(i => d3.range(7).map(j => 0)); const t = Date.now(); for (let i = 0; i < 7; ++i) { for (let j = 0; j < 7; ++j) { for (let n = 0; n < 49; ++n) { re[i][j] += A[n] * orbitals[n].orbital[i][j] * Math.cos(orbitals[n].energy * t); im[i][j] += A[n] * orbitals[n].orbital[i][j] * Math.cos(orbitals[n].energy * t); } } } let s = 0; for (let i = 0; i < 7; ++i) { for (let j = 0; j < 7; ++j) { de[i][j] = re[i][j] * re[i][j] + im[i][j] * im[i][j]; s += de[i][j]; } } for (let i = 0; i < 7; ++i) { for (let j = 0; j < 7; ++j) { de[i][j] /= s; } } for (let l = 0; l < 49; ++l) { data[l].value = de[data[l].i][data[l].j]; } svg .selectAll("rect") .attr("fill", d => d3.interpolateReds(Math.abs(d.value))); yield Promises.delay(200, svg.node()); } }35:[] 31:{"id":"87de901ec7353d78","slug":"qh","title":"QH","thumbnail":"cca88fb4f77f359a290787c2faa7c65ac11db290b8a4043dae7dd82243791282","access_level":"public","create_time":"$D2021-03-10T10:51:29.912Z","update_time":"$D2021-03-11T04:37:18.224Z","view_time":null,"first_public_time":"$D2021-03-10T18:16:04.562Z","has_unpublished_edits":null,"fork":false,"stars":0,"starred":false,"owner":{"type":"team","id":"7a4ecd430a93cd6f","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":"starter_2024"},"version":153,"publish_version":153,"role":null,"authored":false,"authors":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null}],"members":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"QH","cells":[{"id":0,"value":"md`# QH`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":31,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":133,"value":"$34","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":13,"value":"orbitals = {\n const g = (i, j) => 7 * ((i + 7) % 7) + ((j + 7) % 7);\n const l = k => [Math.floor(k / 7), k % 7];\n const L = d3.range(49).map(i => d3.range(49).map(j => 0));\n const factor = 1.0;\n for (let i = 0; i < 7; ++i) {\n for (let j = 0; j < 7; ++j) {\n L[g(i, j)][g(i, j)] = 4 + factor * V(i, j);\n L[g(i, j)][g(i + 1, j)] = -1;\n L[g(i, j)][g(i - 1, j)] = -1;\n L[g(i, j)][g(i, j + 1)] = -1;\n L[g(i, j)][g(i, j - 1)] = -1;\n }\n }\n const r = matrixEig.eig(L);\n const eig = [];\n for (let k = 0; k < 49; ++k) {\n const e = {};\n e.energy = r.eigenvalues.real[k];\n const v = r.eigenvectors.right.slice(49 * k, 49 * (k + 1));\n e.orbital = d3.range(7).map(i => d3.range(7).map(j => v[g(i, j)]));\n eig.push(e);\n }\n return eig.sort((ea, eb) => (ea.energy < eb.energy ? -1 : 1));\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":7,"value":"V = function(i, j) {\n const n = [0, 2, 4, 5, 7, 9, 11];\n const int = (i, j) => (n[j] - n[i] + 12) % 12;\n const intNameEnergy = [\n ['unison', 3],\n ['semitone', 7],\n ['tone', 5],\n ['minor third', 3],\n ['major third', 2],\n ['fourth', 4],\n ['tritone', 7],\n ['fifth', 1],\n ['minor sixth', 2],\n ['major sixth', 3],\n ['minor seventh', 6],\n ['major seventh', 5]\n ];\n return intNameEnergy[int(i, j)][1];\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4,"value":"import { matrixEig } from '@fil/hello-matrix-eig'","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":9,"value":"import { pt } from '@fil/pt'","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":15,"value":"d3 = require('d3')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":38,"value":"require(\"d3-scale-chromatic@1\", \"d3-interpolate@1\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 32:"$31:body"