六角形のパターンの原子ABABABの配列がある場合、 Mathematica を使用して作成するにはどうすればよいですか。この配列で六角形の格子（無限）を使用するため、各原子AはB原子のみに囲まれ、その逆も同様です。
コメント
- Hola Jose 、Mathematica.SEへようこそ。グラフィックラティス、必然的に有限のプロット、またはラティスの分析的記述を意味しますか？おそらく、それをどうするつもりかについてもっと詳しく説明できるので、簡単に手伝うことができます。
- たとえば、この google.es / … ですが、代わりに2つの原子が1つ（グラフェン）
- 関連： 19165 、 14632 。
- 関連： Wolframデモ
- 固体物理学の知識があればそれを容易にします。
回答
2Dで
unitCell[x_, y_] := { Red , Disk[{x, y}, 0.1] , Blue , Disk[{x, y + 2/3 Sin[120 Degree]}, 0.1] , Gray, , Line[{{x, y}, {x, y + 2/3 Sin[120 Degree]}}] , Line[{{x, y}, {x + Cos[30 Degree]/2, y - Sin[30 Degree]/2}}] , Line[{{x, y}, {x - Cos[30 Degree]/2, y - Sin[30 Degree]/2}}] } これによりユニットセルが作成されます
Graphics[unitCell[0, 0], ImageSize -> 100] 
ラティスに配置します
Graphics[ Block[ { unitVectA = {Cos[120 Degree], Sin[120 Degree]} ,unitVectB = {1, 0} }, Table[ unitCell @@ (unitVectA j + unitVectB k) , {j, 1, 12} , {k, Ceiling[j/2], 20 + Ceiling[j/2]} ] ], ImageSize -> 500 ] 
3Dで
unitCell3D[x_, y_, z_] := { Red , Sphere[{x, y, z}, 0.1] , Blue , Sphere[{x, y + 2/3 Sin[120 Degree], z}, 0.1] , Gray , Cylinder[{{x, y, z}, {x, y +2/3 Sin[120 Degree], z}}, 0.05] , Cylinder[{{x, y, z}, {x + Cos[30 Degree]/2, y - Sin[30 Degree]/2, z}}, 0.05] , Cylinder[{{x, y, z}, {x - Cos[30 Degree]/2, y - Sin[30 Degree]/2, z}}, 0.05] } Graphics3D[ Block[ {unitVectA = {Cos[120 Degree], Sin[120 Degree], 0}, unitVectB = {1, 0, 0} }, Table[unitCell3D @@ (unitVectA j + unitVectB k), {j, 20}, {k, 20}]] , PlotRange -> {{0, 10}, {0, 10}, {-1, 1}} ] 
コメント
- わかりました…：D
- 2Dと3Dの両方を検討するのが好きです！
回答
2Dでは、
Manipulate[( basis = {{s, 0}, {s/2, s Sqrt[3]/2}}; points = Tuples[Range[0, max], 2].basis; Graphics[Point[points], Frame -> True, AspectRatio -> Automatic]) , {s, 0.1, 1} , {max, 2, 10} ] 回答
別の方法は、GeometricTransformationを使用することです。これはレンダリングが速くなる可能性がありますが、。
With[{base = Line[{ {{-(1/2), -(1/(2 Sqrt[3]))}, {0, 0}}, {{0, 0}, {0, 1/Sqrt[3]}}, {{0, 0}, {1/2, -(1/(2 Sqrt[3]))}} }] }, Graphics[{ GeometricTransformation[ base, Flatten@Array[ TranslationTransform[ {1/2, -(1/(2 Sqrt[3]))} + {#1 + If[OddQ[#2], 1/2, 0], #2 Sqrt[3]/2} ] &, {16, 16} ] ] }] ] 
{16, 16}をivに置き換えると、$IterationLimitを増やすと機能しませんid = “b4797a8173″>
。
回答
六角形のグリッドの作成に役立つリソース関数はほとんどありません。 。以下のコードは、 HextileBins の例からのものです。
HextileBins
hexes2 = Keys[ ResourceFunction["HextileBins"][ Flatten[Table[{x, y}, {x, 0, 16}, {y, 0, 12}], 1], 2]]; Graphics[{EdgeForm[Blue], FaceForm[Opacity[0.1]], hexes2}] 
lsBCoords = Union[Flatten[First /@ hexes2, 1]]; Graphics[{EdgeForm[Blue], hexes2 /. Polygon[p_] :> Line[Append[p, First[p]]], Red, PointSize[0.02], Point[lsBCoords]}] 
HexagonalGridGraph
（この関数はWolfram Researchによって送信されていることに注意してください。）
grHex = ResourceFunction["HexagonalGridGraph"][{16, 12}] 
lsVCoords = GraphEmbedding[grHex]; lsVCoords[[1 ;; 12]] (* {{0, 0}, {0, 2}, {Sqrt[3], -1}, {Sqrt[3], 3}, {2 Sqrt[3], 0}, {Sqrt[ 3], 5}, {2 Sqrt[3], 2}, {2 Sqrt[3], 6}, {3 Sqrt[3], -1}, {3 Sqrt[3], 3}, {2 Sqrt[3], 8}, {3 Sqrt[3], 5}} *) grHexPolygons = Map[Polygon@(List @@@ #)[[All, 1]] &, FindCycle[grHex, {6, 6}, All]] /. v_Integer :> lsVCoords[[v]]; Graphics[{EdgeForm[Blue], FaceForm[Opacity[0.2]], grHexPolygons}] 