We use a combination of symmetry analysis and high-throughput density functional theory calculations to search for new ferroelectric materials. We use two search strategies to identify candidate materials. In the first strategy, we start with non-polar materials and look for unrecognized energylowering polar distortions. In the second strategy, we consider polar materials and look for related higher symmetry structures. In both cases, if we find new structures with the correct symmetries that are also close in energy to experimentally known structures, then the material is likely to be switchable in an external electric field, making it a candidate ferroelectric. We find sixteen candidate materials, with variety of properties that are rare in typical ferroelectrics, including large polarization, hyperferroelectricity, antiferroelectricity, and multiferroism.Ferroelectrics, which are materials that have a ground state polar phase that can be switched to a symmetryequivalent structure by the application of an external electric field, have been studied and used in applications for many years. Much of the attention on ferroelectrics has focused on prototypical examples from the perovskite oxides, like PbTiO 3 and BiFeO 3 . However, in recent years, there has been a renewed theoretical and experimental interest in discovering and understanding the properties of new ferroelectric materials [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. These new materials have shown a variety of new or rare behaviors, including hybrid improper ferroelectricity, hyperferroelectricity, topological defects/domain walls, etc. In addition, there is interest in and need for ferroelectrics with improved functionality for various applications, including magnetoelectrics, Pb-free piezoelectrics, room temperature multiferroics, solar energy converters, silicon compatible ferroelectrics, ferroelectric catalysts, etc. [15][16][17][18][19][20] High-throughput first principles density functional theory (DFT) calculations have been used increasingly in recent years as a tool for materials discovery and screening [21][22][23][24][25][26]. DFT calculations using semilocal functionals are generally reliable tools for computing ground state properties and the energy differences between closely related phases, which is the primary screening tool needed to identify ferroelectrics. The main obstacle to finding new ferroelectric materials computationally is identifying the relevant polar and non-polar structures to consider. Any insulator known experimentally to have both a polar and a non-polar phase has likely already been recognized as a potential ferroelectric, necessitating a search for new structures of known compounds. However, systematically searching all insulating compounds for possible new phases is too computationally demanding to attempt systematically.In this work, we employ two strategies to identify new ferroelectrics, which we apply to a much larger range of materials than related previous searches [1,10,14,[27][28][29][30][31][32][33][34...