Abstract
Many butterflies, birds, beetles, and chameleons owe their spectacular colors to the microscopic patterns within their wings, feathers, or skin. When these patterns, or photonic crystals, result in the omnidirectional reflection of commensurate wavelengths of light, it is due to a complete photonic band gap (PBG). The number of natural crystal structures known to have a PBG is relatively small, and those within the even smaller subset of notoriety, including diamond and inverse opal, have proven difficult to synthesize. Here, we report more than 150,000 photonic band calculations for thousands of natural crystal templates from which we predict 351 photonic crystal templates â including nearly 300 previously-unreported structures â that can potentially be realized for a multitude of applications and length scales, including several in the visible range via colloidal self-assembly. With this large variety of 3D photonic crystals, we also revisit and discuss oft-used primary design heuristics for PBG materials.
Similar content being viewed by others
Introduction
When will a 3D crystal possess an omnidirectional photonic band gap (PBG)? By definition, when there exists a range of frequencies which are not transmittable through the crystal due to the interactions of waves moving through media of different permittivity. For a material whose spatially dependent dielectric constant is given by ε(r), the transmittable frequencies are given by the eigenvalues of
as derived from Maxwellâs equations and solved in reciprocal space, where H(r) is the magnetic field, Ï is the frequency, and c is the speed of light. Before the introduction of efficient computational packages such as MPB, MEEP, or COMSOL1,2,3, solving this equation was mathematically burdensome, thus the photonics community developed several heuristics to guide the design of photonic crystals, both from the more straightforward 2D and select 3D photonic crystals.
Yablonovitch4, John5, and Ho, et al.6 pointed to the sphericity of the Brillouin Zone (BZ), which is the reciprocal space analog to the real space primitive unit cell, as an important feature for predicting which photonic crystals have a PBG. They reasoned that cubic face-centered (cF) and body-centered (cI) lattices are the most likely targets for PBGs, as they have the most spherical BZs.
The study of 2D PBGs has also provided useful insight into photonic crystal design. Meade, et al. explain the origins of PBGs in 2D crystals by looking at the modes of the photonic band structure defining the PBG. They showed that PBGs occur when there is a large shift in where the electric energy density is localized, where the electric energy density is given by integrating the product of the electric field E(r) and displacement field D(r)â=âε(r)E(r). In the band below the PBG, the electric energy density primarily resides in the medium with the higher dielectric constant (hereon âdielectricâ); this band is known as the âdielectric bandâ where dielectric constant ε is the square of the index of refraction n. Above the PBG, the electric energy density is in the medium with the lower dielectric constant (hereon "airâ); this band is known as the âair band.â Meade, et al. suggest one must aim to decrease the energy of the dielectric band in order to increase the size of the PBG. For a PBG between the transverse electric (TE) modes, this requires that regions of dielectric be connected7. Likewise, the importance of energy localization suggests that the greater the difference in ε of the two regions, the larger the PBG, as this will decrease the similarity between the dielectric and air band. These 2D principles are used to understand and design 3D photonic crystals, the terminology of electric and air band has become conventional in 3D PBG crystals.8
In this work, we revisit these heuristicsâconnectivity of the dielectric, increasing gap size with increasing ε, and spherical BZâgiven a large dataset of 3D photonic band structure calculations, 151,593 in total, and discuss, despite their usefulness, the inadvertent restrictions that they may have imposed on the design search space of PBG crystals to structures near diamond and inverse opal, both of which are face-centered cubic structures that exhibit PBGs when the dielectric forms a continuous network6,8,9,10,11,12,13,14,15,16,17. From our calculations, we demonstrate that an omnidirectional PBG can be supported in all Bravais lattices, in systems of connected and disconnected high dielectric media, and that in some systems, the largest band gap is possible with an intermediate dielectric constant. Furthermore, we provide a comprehensive dataset of PBG crystal possibilities, opening multiple avenues for further study.
Results
Our data set consists of 2714 crystal structures from multiple sources18,19. Each structure is used as a template, with identical dielectric or air spheres on every lattice site for the "directâ or "inverseâ versions, respectively. We screened each structure for PBGs between the first 20 bands across the structural parameters of volume filling fraction Ï (from 0â1) and dielectric constant ε (from 4â16 in reduced units), resulting in 151,593 band structure calculations via MIT Photonic Bands (MPB)1 and managed with the signac data framework20. PBGs are reported as dimensionless percentages, i.e., the range of frequencies within the PBG (ÎÏ) divided by the mid-gap frequency (Ï*)1. Of the photonic band structures computed, 12,778 contain PBGs of size 0.1% or larger for 0.022ââ¤âÏââ¤â0.711. The lowest ε resulting in a PBG is εâ=â4 (agreeing with literature)6,13, and many structures have PBGs with values of ε as low as 6. In total, we find 474 unique PBGs >0.1% in 351 structures (some structures exhibit PBGs at more than one band location at different filling fractions). A summary of the PBGs is shown in Fig. 1, with the maximum gap size for each location (above band 2 and higher) indicated by the size and color of the circles. PBGs under the first band or between bands 1 and 2 are physically impossible, as the first two bands will approach zero frequency at the center of the BZ. The gap atlas, selected photonic band structures, and isosurface representations of structures mentioned in the main text are provided in the Supplementary Information. Photonic band structures for templates that are of special interest or mentioned later in the text are shown in Fig. 1b. Results were validated where possible against previous literature1,6,10,12,21,22. Validation data sets are provided in the Supplementary Information. The maximum PBG size (for any Ï and with εâ=â16) for the well-known diamond, inverse diamond, inverse opal, and inverse simple cubic are approximately 15, 34, 8, and 12% above bands 2, 2, 8, and 5, respectively.
We find PBGs in all 14 Bravais lattices, including the more asymmetric monoclinic and triclinic lattices. In Fig. 2a, we find that PBGs most often exist in cF and cI lattices, concurring with the earlier postulate, as well as rhombohedral (hR) and cubic primitive (cP) lattices, despite their aspherical BZs. As expected, large PBGs (>25%) tend to occur in structures with diamond or gyroid-like topologies, but less expectedly, they also occur in tI lattices, which can have highly aspherical BZs (Fig. 2b(i)). We can compare the sphericity of the BZs using the isoperimetric quotient, defined as \(\frac{36\pi {V}^{2}}{{A}^{3}}\), where IQâ=â1 for perfect spheres. The IQ of hR, cP, and tI lattices are at most Ï/6â=â0.523, whereas the IQ of cF and cI lattices are 0.7533 and 0.7404, respectively. Thus a nearly spherical BZ is neither necessary nor sufficient to produce a PBG, although it will often correlate with the largest PBGs. Discussions of the point group symmetry of the structures and Wyckoff sites, sphericity of the BZ, space group, and angles between nearest-neighbor vectors are provided in the Supplementary Information.
The space group of the crystal strongly constrains the presence and locations of band crossings. The relationship between the crystal space group and mandatory band crossings is well understood23. However, the relationship between the space group and the absence of band crossings is nontrivial and of interest for band gap design24. Our dataset shows strong correlations between the Bravais lattice and gap locations. We find that a large fraction of PBGs between bands 2 and 3 are in cF and tI structures, and PBGs between bands 5 and 6 are primarily found in cP structures, as shown in Fig. 2a.
Connectivity
Next, we look at the 2D heuristic of connectivity of the dielectric7. For any structure realized with monoatomic spheres on every lattice site, there are two connectivity thresholds defined by the geometry. The first, which we denote ÏN, is min(Ï) such that the dielectric forms a continuous network, defined for both direct and inverse structures. The second, which we denote ÏT, is min(Ï) such that the dielectric spheres on any two lattice sites touch and can only be defined for direct structures. A schematic of these thresholds is provided in the Supplementary Information. For all structures, ÏTââ¤âÏN, and for some direct structures, ÏNâ=âÏT as the lattice sites are evenly distributed. As shown in Fig. 2b, we find that many PBGs, including large PBGs, occur for Ïâ<âÏN and Ïâ<âÏT, especially in PBGs at high frequencies, as these modes can more easily travel between disconnected regions of the dielectric.
Variation with dielectric strength (Epsilon)
Finally, in our dataset, we find 41 of 474 cases where the PBG size is non-monotonic with ε, as seen in Fig. 2d(iii). The most dramatic non-monotonicity occurs for lithium oxide (also known as the fluorite or the C1 structure), whose PBG between bands 17 and 18 increases for 6ââ¤âεââ¤â9 but decreases to half its maximum size for 9ââ¤âεââ¤â16, as seen in Fig. 2d(ii). The other PBG, between bands 8â9, was previously theorized in Maldovan et al.11 and is monotonic with ε. This non-monotonic relationship is potentially due to various factors, e.g., the location of the bandgap; whether the dielectric spheres are commensurate with the frequency of the PBG; occurring within a density regime where Mie scattering may increase the complexity of the mode configurations25,26. Also possible is that as the dielectric constant changes, the optimal topology changes, as suggested in Men, et al.13, and only for intermediate ε is an optimized geometry accessible. Similar non-monotonicity was observed in an experimental system and mentioned but not thoroughly investigated27. We note that all structures studied in this work are studied as perfect crystals, not results from self-assembly.
Field analysis
In simple terms, a photonic band gap will occur when two neighboring electromagnetic modes exhibit modes with sufficiently different frequencies, which is proportional to their energy. The ways to engineer this difference in frequency are highly constrained with only 2D geometries. In 3D, the picture is understandably more complex. Here we parallel the field analysis done for 2D by Meade, et al.7 to understand the 3D analogs to the âdielectricâ and âairâ bands.
We calculated the electric displacement field D above and below the gap for select structures found to have a PBG. From these vector fields, we simplified them into flow diagrams of the electric energy localization using the plotting package Mayavi.28 Several motifs (formally defined as general patterns or distinctive features) emerged from these flow diagrams, as shown in Fig. 3a(iâviii). More details on the reduction from vector fields to flow diagrams are provided in the Supplementary Information. We quantify the fraction of electric energy found in the dielectric with the concentration factor f, given by:
where Vε is the high dielectric region. For each PBG, we highlight four properties of the eigenmodes delineating the PBG: (1) band number, (2) Î representation of the bands, (3) f, and (4) the motif that most closely resembles the displacement field lines. These properties are computed for the filling fraction that maximizes the PBG of a given structure and for the wavevectors k corresponding to the maxima and minima of the bands below and above the PBG, respectively. The dotted line signifies ÏN (for the structures shown, ÏNâ=âÏT).
In 2D, the dielectric band typically has 50% or more of the electric energy density localized in the dielectric compared to the air band7. In contrast, from band 2 to 3 in diamond (Fig. 3b), the concentration factor f decreases by 0.1, i.e., 10% more of the electric energy is in the dielectric in band 2 than band 3. There is no change in concentration factor between bands 8 and 9, and the concentration factor of band 15 is higher than that of band 14. Thus, band 15 is more of a âdielectric bandâ than band 14.
The motifs of electric energy localization in 3D delineating a band gap do not fall neatly into the categories of âdielectricâ and âairâ bands. Instead, they fall into two broad categories: those that travel between adjacent unit cells (Fig. 3a(iâiii)) and those that remain confined to a single region of dielectric (Fig. 3a(ivâviii)). For modes with field conformations similar to Fig. 3a(iâiii), increasing the connectivity of the dielectric reduces the mode frequency. Therefore if this motif is found in the band below a gap (but not above the gap), increasing connectivity will enlarge the PBG, as seen in the PBG between bands 2â3 in Fig. 3b. However, if both modes or the mode above the gap exhibit this type of motif, increasing the connectivity will reduce or close the PBG, such as the 8â9 and 14â15 PBG in Fig. 3b.
We can use a similar analysis to explain the relationship between ε and PBG size in the lithium oxide crystal structure shown in Fig. 2b. The heuristic of monotonicity of gap size with ε assumes a mode configuration where the PBG is enlarged by higher localization in the dielectric, e.g., when there is a distinct air band and dielectric band. However, as the PBG in lithium oxide decreases as the filling fraction approaches the connectivity threshold, this assumption is invalid, and the relationship of ÎÏ/Ï* and ε is complex.
Discussion
Our data-driven exploration of the possible space of photonic band gap crystal structures shows that for many photonic crystals, it is clear that no single design rule applies to all PBGs. Many current heuristics provide general guidance to designing crystals for photonic band gaps, but we find many interesting counter-examples for every heuristic. Consequently, further inquiry is needed to predict PBGs in 3D, especially at higher frequencies, where these design rules most typically fail.
The 291 new PBG structures reported here should make for interesting targets, in particular for colloidal systems where diamond and similar structures have been difficult to synthesize, only recently realized with tetrahedral clusters of patchy particles29. Several of these structures shown to have PBGs are known to assemble on colloidal length scales: lithium oxide (6% in direct form) and simple chiral cubic (3.4% PBG in direct form, 27% PBG in inverse form) can be self-assembled using polyhedral nanoparticles30,31, the Clathrate II structure (33% PBG in inverse form) was recently found using DNA-programmed self-assembly32, and the AB13 structure (13.3% PBG in inverse form) is easily achievable with binary mixtures of hard spheres33,34. Furthermore, there are structures, previously unstudied in the colloidal community, that can be new targets for colloidal self-assembly, such as tI835 (6.0% PBG in direct form, 25% PBG in inverse form) or cI1636 (6.4% PBG in direct form, 18.4% PBG in inverse form, the latter of which was previously reported37,38). However, the choice of structure is largely fabrication- and application-dependent and is often best decided by those trying to realize the structure. We have developed an open-source database of photonic structures, which can be used to browse and download the structures detailed in the study so that further researchers can gain inspiration and insight.
Methods
Data management
The data for this project was managed using signac, and the workflow was managed by signac-flow20 in a multi-level project. The top level of the project consists of statepoints associated with the structural data. Inside each structure statepoint was an additional project managing the statepoints containing radii and dielectric constant.
Structure retrieval and conversion
Crystallographic Information Files (".cifâ format) were downloaded from the Crystallographic Open Database (COD) and the Inorganic Crystallographic Structure Database (ICSD)18,19. Information regarding atomic constituents and structure name was taken directly from the .cif files, regardless of the correctness or conditions of the original authorsâ publication or data. Additionally, we included structures found in unrelated simulations in the Glotzer group, some of which have no natural analog. These position files were then converted to a set of lattice vectors and a fractional basis, with the first lattice vector normalized to unit length.
Input parameters for MPB
All PBG calculations were conducted using MPB, a scheme-based eigenmode solver developed at MIT1. MPB calculates the photonic band structure through iterative planewave eigenmode searching. In MPB, the input parameters required are: (1) lattice vectors and fractional basis, (2) particle radius and dielectric constant, (3) fractional k-point path in reciprocal space, (4) k-point interpolation, (5) resolution, (6) mesh size, and (7) number of bands to calculate.
(1) and (2) define the statepoints for this study. (3) was generated using the procedure below. (4) was set to 16 in the initial screening, and set to 32 in complex structures or during the second round of screening. (5) was set to 1 in the initial screening, and 6 during the second round of screening. (6) was set to 5. (7) was set to 20. Only those statepoints with PBGs during initial screening were included in the second round of screening. Those structures found not to have a PBG with minimal k-points would not have a PBG with additional k-points, which would only reduce or close the calculated PBG.
Radius screening
For each structure, particle radii were initialized to râââ(0.0,â1.5) and Îrâ=â0.01 and calculations run. MPB outputs were then queried, and radii with filling fractions outside of [0.0,â1.0] for a structure were removed. For each structure, radii were added uniformly at smaller intervals where necessary, such that each structure was screened with at least 20 radii for both dielectric spheres and air spheres.
Dielectric constants
For each structure, we initially assumed a dielectric constant εâ=â16, the highest theorized dielectric constant for solid, translucent media39. It has been postulated that any crystal will exhibit a PBG in the limit of infinite dielectric constant8, thus we tested a wide variety of crystals, even if difficult to realize in experiment. After the initial computations at εâ=â16, we performed computations for structures and statepoints found to have PBGs at lower ε, down to εâ=â2 or such that ÎÏâ=â0.
k-point path generation
The first Brillouin zone was computed by finding the Voronoi tessellation of the reciprocal lattice for each structure using the Voronoi functionality of SciPy40. The path through the BZ was taken to be a highly redundant traversal along the edges of each face of the surface of the BZ, including to and from the face center and the Î-point. This path, though redundant, is guaranteed to exhaustively sample the high symmetry points at which band extrema occur. An example of a highly redundant k-point path is provided in the Supplementary Information.
For field analysis, the irreducible BZ was calculated to ensure correct k-point labels in the calculation of Î representation.
Structural analysis
In order to calculate the connectivity thresholds and bond angles, we used the Python package pymatgen41, an open-source materials analysis package.
For every structure, we computed the radius r and filling fraction Ï for up to two different thresholds, (1) corresponding to the first peak in the radial distribution function (RDF), denoted ÏT, and (2) corresponding to the radius at which the spheres on the lattice sites are connected in a continuous network, denoted ÏN. ÏT is ill-defined for inverse structures. (2) was found for direct structures using pymatgen and by analyzing the voxelizations of inverse structures using the skimage package40. For some direct structures, these thresholds are the same.
Space groups and symmetry information were found using the open-source package spglib42 and 3D visualization was done with the open-source package mayavi28.
Field analysis
For MPB, there is an option to output field files for a given computation at every \(\overrightarrow{k}\) and band. We took the electric fields in these files and visualized them using mayavi28. The flow reduction was made using 1â2 unit cells of the structure with a resolution of 24â40 grid points. Motifs were drawn based upon analysis of these flow diagrams. A schematic of this process is provided in the Supplementary Information.
The concentration factor was computed using (2) as described in Meade, et al.7.
Î classification was done using character tables found in the Bilbao Crystal Database43 and Dresselhaus, et al.23.
Data availability
Data for all PBG structures and MPB input parameters is publicly available at https://glotzerlab.engin.umich.edu/photonics/. A summary of the structures analyzed in this study is provided in Supplementary Data 1. Further data generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Johnson, S. & Joannopoulos, J. Block-iterative frequency-domain methods for Maxwellâs equations in a planewave basis. Opt. Express 8, 173 (2001).
Oskooi, A. F. et al. Meep: a flexible free-software package for electromagnetic simulations by the FDTD method. Computer Phys. Commun. 181, 687â702 (2010).
Zimmermann, W. B. J. Introduction to COMSOL Multiphysics. 1â26 (Stockholm, Sweden, 2006).
Yablonovitch, E. Inhibited spontaneous emission in solid-state physics and electronics. Phys. Rev. Lett. 58, 2059â2062 (1987).
John, S. Strong localization of photons in certain disordered dielectric superlattices. Phys. Rev. Lett. 58, 2486â2489 (1987).
Ho, K. M., Chan, C. T. & Soukoulis, C. M. Existence of a photonic gap in periodic dielectric structures. Phys. Rev. Lett. 65, 3152â3155 (1990).
Meade, R. D., Rappe, A. M., Brommer, K. D. & Joannopoulos, J. D. Nature of the photonic band gap: some insights from a field analysis. J. Optical Soc. Am. B 10, 328 (1993).
Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. Photonic Crystals (Princeton University Press, 2011).
Chan, C. T., Ho, K. M. & Soukoulis, C. M. Photonic band gaps in experimentally realizable periodic dielectric structures. Europhys. Lett. 16, 563â568 (1991).
Sözüer, H. S. & Haus, J. W. Photonic bands: simple-cubic lattice. J. Optical Soc. Am. B 10, 296 (1993).
Maldovan, M., Ullal, C. K., Carter, W. C. & Thomas, E. L. Exploring for 3D photonic bandgap structures in the 11 f.c.c. space groups. Nat. Mater. 2, 664â667 (2003).
Sözüer, H. S., Haus, J. W. & Inguva, R. Photonic bands: convergence problems with the plane-wave method. Phys. Rev. B 45, 13962â13972 (1992).
Men, H., Lee, K. Y. K., Freund, R. M., Peraire, J. & Johnson, S. G. Robust topology optimization of three-dimensional photonic-crystal band-gap structures. Opt. Express 22, 22632 (2014).
Cersonsky, R. K., Dshemuchadse, J., Antonaglia, J., van Anders, G. & Glotzer, S. C. Pressure-tunable photonic band gaps in an entropic colloidal crystal. Phys. Rev. Mater. 2, 125201 (2018).
Ducrot, Ã., He, M., Yi, G.-R. & Pine, D. J. Colloidal alloys with preassembled clusters and spheres. Nat. Mater. 16, 652â657 (2017).
Leung, K. M. Diamondlike photonic band-gap crystal with a sizable band gap. Phys. Rev. B 56, 3517â3519 (1997).
Maldovan, M. & Thomas, E. L. Diamond-structured photonic crystals. Nat. Mater. 3, 593â600 (2004).
Graulis, S. et al. Crystallography open databaseâan open-access collection of crystal structures. J. Appl. Crystallogr. 42, 726â729 (2009).
Bergerhoff, G., Hundt, R., Sievers, R. & Brown, I. D. The inorganic crystal structure data base. J. Chem. Inf. Modeling 23, 66â69 (1983).
Adorf, C. S., Dodd, P. M., Ramasubramani, V. & Glotzer, S. C. Simple data and workflow management with the signac framework. Comput. Mater. Sci. 146, 220â229 (2018).
Moroz, A. Three-dimensional complete photonic-band-gap structures in the visible. Phys. Rev. Lett. 83, 5274â5277 (1999).
Chen, H., Zhang, W. & Wang, Z. Comparative studies on photonic band structures of diamond and hexagonal diamond using the multiple scattering method. J. Phys. 16, 741â748 (2004).
Dresselhaus, M. S., Dresselhaus, G. & Jorio, A. A. Group Theory (Springer Berlin Heidelberg, 2008).
Watanabe, H. & Lu, L. Space group theory of photonic bands. Phys. Rev. Lett. 121, 263903 (2018).
Antonoyiannakis, M. I. & Pendry, J. B. Mie resonances and bonding in photonic crystals. Europhys. Lett. 40, 613â618 (1997).
Vandenbem, C. & Vigneron, J. P. Mie resonances of dielectric spheres in face-centered cubic photonic crystals. J. Opt. Soc. Am. A 22, 1042 (2005).
Rosenberg, A. et al. Guided resonances in asymmetrical GaN photonic crystal slabs observed in the visible spectrum. Opt. Express 13, 6564 (2005).
Ramachandran, P. & Varoquaux, G. Mayavi: 3D visualization of scientific data. Comput. Sci. Eng. 13, 40â51 (2011).
He, M. et al. Colloidal diamond. Nature 585, 524â529 (2020).
Cadotte, A. T., Dshemuchadse, J., Damasceno, P. F., Newman, R. S. & Glotzer, S. C. Self-assembly of a space-tessellating structure in the binary system of hard tetrahedra and octahedra. Soft Matter 12, 7073â7078 (2016).
Damasceno, P. F., Karas, A. S., Schultz, B. A., Engel, M. & Glotzer, S. C. Controlling chirality of entropic crystals. Phys. Rev. Lett. 115, 158303 (2015).
Lin, H. et al. Clathrate colloidal crystals. Science 355, 931â935 (2017).
Hachisu, S. & Yoshimura, S. Optical demonstration of crystalline superstructures in binary mixtures of latex globules. Nature 283, 188â189 (1980).
Bartlett, P., Ottewill, R. H. & Pusey, P. N. Superlattice formation in binary mixtures of hard-sphere colloids. Phys. Rev. Lett. 68, 3801â3804 (1992).
Cywinski, R., Kilcoyne, S. H. & Scott, C. A. Magnetic order and moment stability in YMn2. J. Phys. 3, 6473â6488 (1991).
Kasper, J. S. & Richards, S. M. The crystal structures of new forms of silicon and germanium. Acta Crystallogr. 17, 752â755 (1964).
Dyachenko, P. N., Miklyaev, Y. V. & Dmitrienko, V. E. Three-dimensional photonic quasicrystal with a complete band gap. JETP Lett. 86, 240 (2007).
Dyachenko, P. N., Dmitrienko, V. E., Miklyaev, Y. V. & Pavelyev, V. S. Complete photonic band gap in icosahedral quasicrystals with a body-centered six-dimensional lattice. Photonic Crystal Materials and Devices VIII. Vol. 6989, 410â417 (International Society for Optics and Photonics, 2008).
Li, H. H. Refractive index of silicon and germanium and its wavelength and temperature derivatives. J. Phys. Chem. Ref. Data 9, 561â658 (1980).
Virtanen, P. et al. SciPy 1.0âfundamental algorithms for scientific computing in Python. Nature methods 17.3, 261â272 (2020).
Ong, S. P. et al. Python Materials Genomics (pymatgen): a robust, open-source python library for materials analysis. Comput. Mater. Sci. 68, 314â319 (2013).
Togo, A. & Tanaka, I. Spglib : a software library for crystal symmetry search. Preprint at http://arxiv.org/abs/1808.01590 (2018).
Aroyo, M. I., Kirov, A., Capillas, C., Perez-Mato, J. M. & Wondratschek, H. Bilbao crystallographic server. II. Representations of crystallographic point groups and space groups. Acta Crystallogr. Sect. A Found. Crystallogr. 62, 115â128 (2006).
Acknowledgements
The authors thank R. Merlin, J. Dshemuchadse, C.S. Adorf, B. VanSaders, C. Dai, and E.F. Hayman for constructive feedback during the writing process. This work was sponsored by the Department of the Navy, Office of Naval Research under ONR award number N00014-18-1-2497 (large scale evaluation of structures). Research and calculation of the EM fields was supported as part of the Center for Bio-Inspired Energy Science, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award # DE-SC0000989. R.K.C. and J.A. acknowledge support from the University of Michigan Rackham Predoctoral Fellowship Program. B.D.D. acknowledges support from the National Science Foundation Graduate Research Fellowship, Grant DGE 1256260. This research was also supported in part through computational resources and services supported by Advanced Research Computing at the University of Michigan, Ann Arbor.
Author information
Authors and Affiliations
Contributions
R.K.C. designed and conducted computations and analyzed data. J.A. conducted field analysis and contributed to data analysis. B.D.D. contributed to data analysis and designed and constructed the database. S.C.G. supervised the work. All authors contributed to the writing of the manuscript and the interpretation of the results.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cersonsky, R.K., Antonaglia, J., Dice, B.D. et al. The diversity of three-dimensional photonic crystals. Nat Commun 12, 2543 (2021). https://doi.org/10.1038/s41467-021-22809-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-021-22809-6
This article is cited by
-
Spin-wave mode coupling in the presence of the demagnetizing field in cobalt-permalloy magnonic crystals
Scientific Reports (2024)
-
Biosensors based on novel nonlinear delta-function photonic crystals comprising weak nonlinearities
Scientific Reports (2022)
-
Evolution of optical Tamm states in a 1D photonic crystal comprising a nanocomposite layer for optical filtering and reflecting purposes
Optical and Quantum Electronics (2022)