Water freezes in a wide variety of low-temperature environments, from meteors and atmospheric clouds to soil and biological cells. In nature, ice usually nucleates at or near interfaces, because homogenous nucleation in the bulk can only be observed at deep supercoolings. Although the effect of proximal surfaces on freezing has been extensively studied, major gaps in understanding remain regarding freezing near vapor-liquid interfaces, with earlier experimental studies being mostly inconclusive. The question of how a vapor-liquid interface affects freezing in its vicinity is therefore still a major open question in ice physics. Here, we address this question computationally by using the forward-flux sampling algorithm to compute the nucleation rate in a freestanding nanofilm of supercooled water. We use the TIP4P/ice force field, one of the best existing molecular models of water, and observe that the nucleation rate in the film increases by seven orders of magnitude with respect to bulk at the same temperature. By analyzing the nucleation pathway, we conclude that freezing in the film initiates not at the surface, but within an interior region where the formation of doublediamond cages (DDCs) is favored in comparison with the bulk. This, in turn, facilitates freezing by favoring the formation of nuclei rich in cubic ice, which, as demonstrated by us earlier, are more likely to grow and overcome the nucleation barrier. The films considered here are ultrathin because their interior regions are not truly bulk-like, due to their subtle structural differences with the bulk.ice | nucleation | molecular simulations | surface freezing | statistical mechanics F reezing of water into ice is ubiquitous in nature and can occur in a wide range of environments, from biological cells (1) and soil (2) to meteors (3) and atmospheric clouds (4). Ice formation is particularly important in the atmospheric processes that exert a major influence on our weather and climate. Indeed, the radiative properties (4) and precipitation propensity (5) of a cloud are both determined by the amount of ice that it harbors. Not surprisingly, the liquid fraction of clouds is an extremely important input parameter in meteorological models (6). Freezing is a first-order phase transition and proceeds through a mechanism known as nucleation and growth. In mixed-phase clouds, which are responsible for the bulk of land surface precipitation (7), freezing is usually nucleationlimited, and individual freezing events are rare. Predicting the timing of such rare freezing events is important in modeling cloud microphysics. However, nucleation rates, which are volume-or area-normalized inverse nucleation times, can only be measured experimentally over a narrow range of temperatures and pressures (8), and extrapolations to other temperatures are nontrivial (9). In nature, homogeneous nucleation of ice is only likely at very low temperatures, and the majority of freezing events in mixed-phase clouds occur heterogeneously, at interfaces provided by external insol...