Multilevel design and construction in nanomembrane rolling for three-dimensional angle-sensitive photodetection

Releasing pre-strained two-dimensional nanomembranes to assemble on-chip three-dimensional devices is crucial for upcoming advanced electronic and optoelectronic applications. However, the release process is affected by many unclear factors, hindering the transition from laboratory to industrial applications. Here, we propose a quasistatic multilevel finite element modeling to assemble three-dimensional structures from two-dimensional nanomembranes and offer verification results by various bilayer nanomembranes. Take Si/Cr nanomembrane as an example, we confirm that the three-dimensional structural formation is governed by both the minimum energy state and the geometric constraints imposed by the edges of the sacrificial layer. Large-scale, high-yield fabrication of three-dimensional structures is achieved, and two distinct three-dimensional structures are assembled from the same precursor. Six types of three-dimensional Si/Cr photodetectors are then prepared to resolve the incident angle of light with a deep neural network model, opening up possibilities for the design and manufacturing methods of More-than-Moore-era devices.

Many analytical and numerical methods have been used for morphological prediction and design guidelines for pre-strained nanomembrane structures [27][28][29][30][31][32] .As early as 1925, Timoshenko has developed the curvature radius formula of the bi-metal thermostat, which is applicable to substrates of similar thickness to nanomembranes 33 .
In the later development, Nikishkov 34 and Hsueh 35 have presented the analytical solutions of the curvature formula suitable for self-rolling nanomembranes.In addition, finite element modeling (FEM), as a numerical method, can intuitively simulate and predict the behavior of nanomembranes after release [28][29][30][31] .Huang et al. have demonstrated a quasistatic FEM model with moving boundaries to simulate the selfrolling process of nanomembranes 29 .Generally, the rolling behavior of completely released nanomembranes from long sides is regarded as a preference in minimal elastic energy 6,28,31,36 .Nevertheless, short-side rolling has been observed in many experimental results, which is interpreted as the dependence of rolling trajectory 6,28 .Idealizing the boundary movement of the sacrificial layer for specific structural formation or analysis in a local area lacks broad applicability and an accurate boundary condition model.
It is crucial to adopt a design approach that takes multiple factors into account to ensure efficient processing and optimal application of micro/nanoscale self-assembled devices.Here, we propose a quasistatic multilevel FEM that emphasizes the geometry and boundary condition of nanomembrane for 3D self-assembly structure.By considering the non-uniformity and geometric dependency of the etching process, we reproduce the anisotropy of boundary movement of the sacrificial layer.The multilevel FEM method is successfully applied in a wide range of material systems, nanomembrane thicknesses, pattern types, and pattern sizes, demonstrating excellent generalizability.We take Si/Cr rectangular nanomembrane as an example to study the mechanism and application of the self-assembly process in detail.An idealized model is used to simulate and calculate the elastic energy changes corresponding to different geometric characteristics quantitatively.The gradual shift observed in the rolling direction reflects the transition of minimal elastic energy and explains the change in 3D structural forming behavior with varying aspect ratios.Based on this FEM model, large-scale, high-yield, and high-uniformity 3D configurable structures are successfully achieved.In addition, we develop a range of 3D structure photodetectors for angle detection of incident light at an accuracy of 10°to demonstrate the potential in manufacturing electronic and photoelectronic devices.

Release of pre-strained bilayer nanomembrane
We established a standard strain nanomembrane release model consisting of a rectangular nanomembrane with width W, length L, and inplane strain gradient Δε, sacrificial layer, and substrate, as shown in Fig. 1a.A fixed edge is deployed to ensure the on-chip self-assembly of nanomembrane, resulting in the morphological difference between opposite edge rolling and adjacent edge rolling.Once the sacrificial layer is etched, the pre-strained nanomembrane will release from the substrate and then self-assemble due to strain gradient, as shown in the right panel of Fig. 1b.During the release process, etching of the sacrificial layer can be modeled as a gradual movement of boundary, in which the velocity and direction of movement depend on the chemical reaction parameters at the solid-liquid interface, as shown in Fig. 1c.Geometric size of nanomembrane will lead to a diversity of etching velocity and direction at the solid-liquid interface, which can be intuitively represented by v W v L , the velocity ratio of the velocity perpendicular to the length direction and width direction of a coordinate on nanomembrane, which is shown in the inset of Fig. 1d.The method of etching FEM is provided in Methods and Supplementary Fig. 1.In previous studies 20,32,[37][38][39][40][41][42][43][44][45][46][47][48] , the structures of nanomembranes exhibit a transition from ring to arch, taper, and finally tubes with an increase in aspect ratios.When the long side is much longer than the short side (W»L or W«L), the etching from the long side will dominate the etching direction, and short-side etching is negligible due to the short distance between opposite edges.When the W/L is close to 1, etching from both sides needs to be considered because the rolling direction may be competitive in similar geometric size (Fig. 1d).For example, Chalapat et al. 38 discovered the transition between ring and arch structures via fabrication of Ti/Al/Cr nanomembrane with increasing aspect ratios.
The above phenomenon shows the potential of structural design with different aspect ratios.In pre-strained nanomembrane with various aspect ratios, there are differences in elastic energy after release from two directions.We set up a series of FEM models that release from different directions to simulate the change in elastic energy during etching.As shown in Fig. 1e, when the etching direction is fixed to adjacent and opposite edge, the gap in elastic energy from FEM calculation indicates that the distinct release direction may result in two types of 3D structures.E adjacent and E opposite represent elastic energy of adjacent and opposite edge release during etching, and E 0 is the elastic energy of unreleased nanomembrane.As for W/L = 0.1 and W/L = 1, the elastic energy decreases fast in the adjacent side release at the initial state, which will cause the adjacent side rolling in the same etching depth.Compared to the common tubular structure rolling from opposite edge, the obvious decrease from adjacent side indicates the possibility of adjacent side rolling and potential of multi-morphic 3D assembly from one pattern.In the situation of W/L = 2, the elastic energy of the opposite rolling is in a low energy state throughout the entire process, proving the change in the prior rolling direction.As mentioned above, etching velocity, etching direction, and aspect ratio are both important factors for 3D structural assembly, and they need to be considered in FEM modeling.

Quasistatic multilevel finite element method modeling
In order to reproduce the formation process and find the mechanism of boundary condition during the release process, a multi-field coupling FEM model is utilized to simulate the etching process in the sacrificial layer via COMSOL Multiphysics software package [49][50][51] .In the etching FEM model, the bilayer is initially fixed with constraints to simulate the unreleased state.Considering the etchant flowing and boundary movement of the etching system, Fick's law, Navier-Stokes equations, and liquid-solid reaction model are introduced to establish the FEM model: where c is etchant concentration, u is velocity of etchant flow, D is diffusion coefficient, k is reaction rate coefficient, n is normal vector, ρ is density, μ is dynamic viscosity coefficient, v is normal velocity of mesh, and M is molar mass.The relationship between etchant concentration and reaction rate depends on the balance of diffusion and chemical reaction process, and the effect of each parameter is discussed (Supplementary Notes 1, 2, and Supplementary Figs.2-5).As can be seen in Fig. 2a, since the FEM model is made up of a set of small units, the etching process of the sacrificial layer is realized by normal movement of unit meshes with a velocity v.
We perform a series of analytical steps in a time sequence and then realize the whole dynamic release process of bilayer nanomembranes.Geometrical features of boundary variation obtained from previous simulations are applied to the feature division along the thickness direction of the Si/Cr bilayer in coordinates (x i , y i ), i = 1, 2,…, node numbers, as shown in Fig. 2b.After geometric partitioning, boundaries under multiple discrete time points are divided into boundary conditions.In the elastic mechanical analysis step, the boundary conditions will be set as a series of constraints in chronological order to realize the continuous release of the nanomembrane (Fig. 2c).
Due to the asymmetric boundary conditions, complex practical release process, and large geometric deformation of the one-ended fixed bilayer nanomembrane structure, it is inaccurate to apply the large deflection equations of thin elastic plates (such as the Föppl-von Kármán equation 52 ) for structural design.Considering only the release from the opposite edge, the 3D assembly of the nanomembranes will exhibit only tubular structures (Supplementary Note 3 and Supplementary Fig. 6).In contrast, the simulation results of the bilayer nanomembranes are guided by the multilevel design model during the release process, demonstrating better accuracy and successfully reproducing morphologies varying with width (Fig. 2d).Quasistatic FEM results prove that the etching trajectory of the sacrificial layer plays a crucial role in tuning the release of bilayer nanomembranes.In addition to the Si/Cr rectangular patterns, the multilevel FEM model has been designed for various types of patterns such as semicircular, triangular, and parallelogram patterns (Fig. 2d, Supplementary Note 4, Supplementary Figs.7-14, and Supplementary Movies 1, 2), as well as in the high-frequency/low-frequency SiN x (LF/HF SiN x ) nanomembrane, NiTi nanomembrane, and VO 2 /Cr nanomembrane systems, which demonstrate the potential for a wide range of applications (Fig. 2e-h, Supplementary Note 5, and Supplementary Figs.15-17).Meanwhile, the method demonstrates cross-scale compatibility for models from the hundred nanometer scale to the hundred micrometer scale (Supplementary Note 6 and Supplementary Figs.18-22).In this study, we choose the rectangular Si/Cr rectangular nanomembrane for detailed analysis and discussion of the self-assembly behavior due to their moderate parameter complexity and high geometrical symmetry.

Experimental verification
To verify the design model, a series of silicon (Si)/chromium (Cr)/ germanium (Ge) heterostructure with gradient aspect ratio was deposited by electron beam evaporation (Supplementary Note 7 and Supplementary Fig. 23).Si/Cr bilayer nanomembrane system is selected for experimentation, in which Cr can introduce large pre-strain and Si is the most widely used semiconductor materials.The lithography and reaction ion etching then define the one-end fixed rectangular sacrificial layer area by photoresist to produce a window for sacrificial layer etching.As shown in Fig. 3a, a series of one-end fixed Si/Cr bilayer nanomembranes with length L = 40 μm are selected as an example, and the width W is increased from 2 μm to 80 μm with a step of 2 μm.The photoresist was removed via acetone immersion and supersonic cleaning to avoid surface residue of organic solvents.After rinsing the surface in deionized water and drying, the Ge sacrificial layer is etched to realize Si/Cr bilayer release.Ge sacrificial layer was removed by wet etching (30% H 2 O 2 , 75 °C).When Si/Cr bilayers in the etching area are released, they are transferred to acetone and dried in critical point dryer in liquid CO 2 to prevent morphological defects caused by conventional drying.When the length L = 40 μm, there are obvious transition regions (W = 8-10 μm, 16-18 μm, 28-30 μm, and 56-58 μm) and stable regions (W = 2-6 μm, 12-14 μm, 20-26 μm, 32-54 μm, and 60-68 μm) in the distribution map of various morphology, which confirms that 3D assembly is stable and related to the aspect ratio of geometric shapes (Fig. 3a).The morphologies of the released nanomembrane exhibit a transitional change as the width increases: five types of 3D structures of ring, arch, helix, taper, and tube (Fig. 3b).For patterns in transition regions, the unstable structure formation     demonstrates the competition of two rolling directions and a relatively small gap in elastic energy.To systematically analyze the relationship between structural transition and aspect ratio, patterns with different lengths L from 20 μm to 80 μm (step of 5 μm) were prepared.As shown in Fig. 3c, the probability distribution of five types of morphology and pattern sizes is discussed from a statistical point of view.With the increase in length, the transition width of each structure increases accordingly due to the variation in aspect ratio.For patterns with large lengths (L = 80 μm), the structure will lack the transition from taper to tube, but it can be fabricated by a design with larger width.The forming morphology of nanomembranes also follows the order of ring-arch-helix-taper-tube with increasing aspect ratio, corresponding to the structure distribution in Fig. 1d and 7.06 ± 1.93 μm (t Si = 90 nm), respectively.We used the above experiments and the force-moment balance formula to estimate the strain magnitude in the Si/Cr layer, that is where E is Young's modules, t is thickness of nanomembrane, and Δε = ε Cr À ε Si is the initial strain difference between Si and Cr layer.We substitute the thicknesses of nanomembranes and tube radius (Fig. 3d,e) into this formula for fitting, and the prestrain of the Cr layer obtained by fitting is ε Cr ~0.65%,while the pre-strain of the Si layer is determined to be ε Si ∼ À 0:55%.The corresponding stress gradient is approximately Δσ = ε Cr E Cr À ε Si E Si ≈2600 MPa, which is consistent with the preset condition in our simulation (Δσ = 2000 MPa), providing an important guarantee for high-accuracy structural design.In addition, we used the thickness and pre-strain of nanomembranes as parameters to calculate the strain relationship in the Si/Cr nanomembrane.The average strain of Si and Cr nanomembranes are also in good agreement with the experiment, providing an effective method for strain characterization of multilayer nanomembrane systems (Supplementary Fig. 31).

Elastic energy of released bilayer nanomembrane
The assembly of nanomembranes is well simulated.However, the essence of formation and the size dependence of different structures need to be explored.The optimized elastic energy calculation is derive from the Nikishkov's analytical solution model 34 (Supplementary Note 9 and Supplementary Figs.32-34).We approximate the geometric parameters of the model as boundary conditions with symmetrical specific divide angles α to combine the effect of etching direction α = arctan v W v L and aspect ratio W/L during release.Patterns with a length L of 40 μm and a width W of 10-80 μm are set with a series of direction angles for unilateral (Supplementary Note 10 and Supplementary Fig. 35) and bilateral (Fig. 4a) release models.The divided boundaries are parallel, and the divided area is classified as Region I.The constraint of the nanomembrane in the Region I will be removed inward from the boundary to simulate directional release.
When Region I is completely released, the nanomembrane in Region II starts to release.Finally, the elastic energy of the released nanomembrane is collected via FEM.The geometric characteristics of FEM results are provided to idealized analytical elastic energy equations to verify the results of simulations.
It can be observed in Fig. 4b that among the relative strain energy of different widths, we can find the direction angle α min with the minimal relative elastic energy shifts to a larger direction angle with increasing width.The α min of the bilateral-released model shifts from ~70°to ~45°when the width increases from 10 μm to 80 μm.This decreased transition of α min indicates long-side rolling along the opposite edge.When the rolling direction decreases to 0°, the rolling will become the classic model of opposite rolling in Fig. 1d.Then, the α min shifts back from ~70°to ~50°after the width decreases from 10 μm to 4 μm.The shift back of α min is due to the dominant of long-side rolling from length direction, which is both restricted by high v W v L and a small aspect ratio W/L ≤ 0.25.Therefore, the pattern of W = 4 μm with α min = ~50°assembles into a ring rather than rolling from adjacent edges.Generally, the existence of the α min verifies the dependency between aspect ratio and etching direction and the relationship to the minimum energy theorem.For other α with higher relative elastic energy, the bilayer nanomembranes show reduced possibilities for rolling.The elastic energy calculation shows the same increasing and decreasing trend of the α min shift.Fluctuations in relative elastic energy derived from estimation errors in selecting idealized regions.
Moreover, the transition of α contributes to the differences in the 3D structures of nanomembranes (Supplementary Note 11 and Supplementary Figs.36-39).Taking L = 40 μm as an example, the probability of occurrence of various structures exhibits gradient change from 4 μm to 80 μm, which corresponds to the shift of α min (Fig. 4c).For each structure, there are five corresponding areas of α min to assemble: ring (30°-40°), arch (50°-80°), helix (65°-75°), taper (50°-65°), and tube (40°-55°).It is confirmed that the structure with the minimal relative elastic energy possesses a higher probability of rolling up along the α min assembling into the corresponding structure (Supplementary Note 12 and Supplementary Fig. 40).As shown in Fig. 4d, the quasistatic FEM model simulates the change of α with etching time during the assembly process, and generally restores the relationship between α and etching time during the sacrificial layer etching process in the experiment (Fig. 4e, Supplementary Movie 3).When the etching finished, the distribution of α are as follows: ring (25°-35°), arch (70°-80°), helix (50°-65°), taper (60°-75°), and tube (40°-55°).The α angle obtained at the end of etching between model calculations and experimental results is highly consistent with the calculated α min corresponding to each structure in Fig. 4c, providing an observation scheme for tracking intermediate states in the assembly process.We noticed that there is an obvious difference of the α between the multilevel model and the experiments in the early stage of etching, which is mainly due to the non-ideal release of sacrificial layer release and local non-uniform etching caused by insufficient accuracy in micro-nano processing (Supplementary Note 13 and Supplementary Fig. 41).
Regarding the Si/Cr bilayer with different widths, we anticipate that the α min will also shift as the width increases with the instruction of FEM-based multilevel design model.With the guidance of the elastic energy map in Fig. 4c, we employed the multilevel design with a single stable energy state to prepare large-scale arrays with only one structure.As shown in Table 1, all five structure arrays (0.8 cm × 0.8 cm) exhibit high yield (Supplementary Note 14 and Supplementary Fig. 42), demonstrating good accuracy in designing 3D structures from selfrolling.

Configurable structures with controlled etching edges
By introducing sequential etching in designed spatiotemporal order at targeted locations, diverse structural configurations can be achieved for identical patterns.The orientation of the patterns and adjustment of the etchant level position are utilized to control etching at specific locations, resulting in the formation of two distinct 3D structures.During the etching process, the solid-liquid interface rises slowly due to the capillarity between the nanomembrane and the etchant (Supplementary Movie 4).The first method of etching control is bottom etch, in which the fixed end is in contact with the etchant first.Then, the liquid level of the etchant will rise due to the capillary force in the etched pattern (Fig. 5a).In this situation, etching happens at the adjacent edges first.(Fig. 5b).Then, another method for etching control is top etch, and the opposite edge is in contact with the etchant first, as shown in Fig. 5c.For etching starting from the opposite edge (Si/Cr top etch, Fig. 5d, Supplementary Note 15, and Supplementary Movie 5), the etchant interface prefers to contact the free end initially than adjacent edges.A quasistatic multilevel FEM model with selective release is designed to guide the 3D assembly (Fig. 5e).In the bottom etch, the boundary of the sacrificial layer will change from adjacent sides rather than the overall layer, artificially reducing the The formation of tubular structures benefits from the wrinkle in the partial release process by bottom etch.In local wrinkle structures, the bilayer tends to roll from the directions of wrinkling, when the assembly for other structures faces a higher energy barrier (Supplementary Fig. 43, Supplementary Movie 6).Compared to models in the normal etching process, the movement of the opposite edge in the top etch is larger at the beginning.More area will be released correspondingly in mechanical FEM modeling, leading to a larger α and v W v L in the local area (Supplementary Note 16 and Supplementary Fig. 44).Therefore, there is a significant variation in bilayer morphology in bottom etching and top etching.Selective etching introduces controllable regulation of etching direction and etching area, playing an important role in structural assembly.The successful preparation of the multi-morphic structure also proves the correctness of elastic energy analysis.By introducing extra boundary conditions, etching direction, and velocity ratio v W v L can also be simulated and realized.With proper control of surface tension, reconfiguration of different structure can be achieved (Supplementary Movie 7).The local etching method is also applicable to the polymorphic design of other patterns (Supplementary Note 17 and Supplementary Figs.45-57).As shown in Fig. 5f, in parallelogram patterns, there is an additional parameter tilt angle that define the geometric feature of patterns.In such patterns, the bottom etch will result in the free-end of the patterns being released at the end of the etch, when the structure will tend to assemble along the shorter diagonal or the free edge of width.Conversely, when the top etching start at the free-end first, and the localized etching will allow the assembly to develop along a longer diagonal section to satisfy the lowest elastic energy 16 .We have successfully utilized inclined etching to assemble the same parallelogram patterns along different diagonals, thus forming two types of structures with significant differences (Fig. 5f, Supplementary Figs.53-57).It is worth noting that the premise for the preparation of polymorphic microstructures is that there are multi-stable strain energy states when the pattern is assembled, and there is an enough elastic energy barrier between the stable states.Otherwise, the pattern will only assemble into one type of structures (Supplementary Figs.46-52).To summarize, we employed the multilevel design method based on the FEM etching and elastic mechanical model to achieve the 3D assembly of pre-strained nanomembranes, offering fresh opportunities for prospective device design through strain engineering.

Angle-sensitive detection of incident light
The multilevel structural design for nanomembranes is also applicable to guide the fabrication of micro/nanoscale electronic devices.A wide range of incident light angle detection of tubular photodetector has been observed, indicating that anisotropic photodetection can be realized via structural design.Here, incident light angle detection can be achieved by combining the angular response difference of a series of photodetectors.As shown in Fig. 6a, a hemispherical omnidirectional incident light controller is designed to achieve a photodetector with incident light traversing the spherical surface at the selected angle (θ, φ).The controller consists of a transparent PMMA hemispherical shell and an array of optical fiber interfaces, in which the laser can be incident at a specific angle through the interface connected to the corresponding coordinates of the spherical surface.The prepared Si/ Cr photodetector is placed on a platform at the same height as the bottom of the incident light controller.Then the controller is calibrated to ensure that the projected coordinate on YZ plane of (90°,0°) input laser port is aligned with the center of photodetectors (Fig. 6a,b, Supplementary Note 20 and Supplementary Figs.60-61).Then, the collected photodetection data will be imported into neural network to analysis the incident angle (Fig. 6c).The Si/Cr photodetectors are prepared by a similar method to the Si/Cr bilayer while Cr layer is defined as a 5 μm channel in the middle of the pattern to realize the photodetection of the Si layer (Supplementary Figs.58, 59).Si/Cr photodetectors exhibit a maximum responsivity of 60 mA/W, response time of 100 ~700 μs, and external quantum efficiency of 7 ~12%, which can effectively respond to 520 nm incident light to achieve photodetection (Fig. 6d, e, and Supplementary Figs.62-64).After measurements of photoresponse of light incident from different coordinates, the photocurrent of each coordinate will be normalized and classified by structure types.To facilitate visual data comparison, the normalization photocurrent of each structure will be displayed through the projection of light controller on YZ plane.To aid in identifying the position of incident light on the projection sphere, we established spherical coordinates denoted by θ and φ.Compared to the unreleased Si/Cr planar photodetector of Fig. 6f, the ring (Fig. 6g) and tube photodetectors exhibit a wider detection angle of high photocurrent 8,23 , confirming that the 3D microstructure contributes to the anisotropic detection of incident light and enhancing the detection of φ in area of θ = 70°-110°.For the arch photodetector (Fig. 6h), the increasing trend of photocurrent from θ = −0°to θ = 90°, which derives from the out-of-plane bending of the arch structure, exhibit angle detection feature in θ direction, and significantly distinguish the incident light from θ < 90°and θ > 90°.Helix photodetector exhibits an inclined high photocurrent detection angle from φ = −30°to φ = 0°b ecause of the asymmetric rolling behavior via instability during the release process helping the angle identifying whether φ < 0°or φ > 0°(Fig.6i).In the taper photodetector, obvious high photocurrent appears in and φ = 130°-140°, 40°-50°, benefiting from multiple reflections within a specific range (Fig. 6j), focusing on the detection in above area.Compared to ring structure, the tube photodetector exhibit a better stability in wide angle detection, which contributes to the detection of φ in area of θ = 50°-130°(Fig.6k).The results above illustrate the sensitivity of photoelectric coupling in 3D configurations and can be extended and utilized for functionalizing 3D-assembled nanomembrane electronic devices.
Neural networks have made significant advances in handling complex environmental information and ambiguous inference rules [53][54][55] .Thus, based on the difference in incident angle detection between the above structures, the deep neural network (DNN) algorithm can establish the angle-sensitive detection model based on photodetectors with different structures.A data set containing 275 photocurrents for each channel from Fig. 6f-k is collected to establish the incident angle prediction model.After training initial DNN models, two sets of photocurrent data are input for additional validation.As shown in Fig. 6l, training sets from the longitudinal (φ = 0°) and latitudinal (θ = 40°) direction exhibit an accuracy of 95% and 78%, respectively.In Fig. 6m, the relationship between DNN-predicted incident angle and experimental angle demonstrates an accuracy of 83% and 71% for incident angle detection in an angle resolution of 10°f rom the longitudinal and latitudinal directions, respectively.As can be seen in Fig. 6n, 83% of the weight of the corresponding angle (φ = 0°) exhibit >10 times compared to the weight other angles, indicating the effectiveness of incident light detection.Profited from high accuracy prediction, these 3D self-assemble photodetectors show potential in wearable devices, intelligent furniture, and intelligent driving systems.

Discussion
In summary, we have established a quasistatic multilevel FEM to simulate and predict the self-assembly 3D structure of pre-strain nanomembrane, solving the formation design problem under the coupling of complex etching trajectory and geometric pattern of 2D nanomembranes.Quasistatic multilevel FEM reproduces the dynamic release process by collecting coordinates of etching trajectory in chronological order and predict the transition of 3D structure varies with aspect ratio.This FEM method shows good applicability in a variety of material systems, precursor pattern types, pattern sizes, and nanomembrane thicknesses, which can also conduct preliminary geometric feature analysis of the self-assembly process and results.Si/ Cr bilayer rectangular nanomembrane is adopted to verify the gradient change of structures.The analysis of elastic energy proves the correlation between aspect ratio, etching trajectory, and final morphologies.By controlling the spatiotemporal order of the etching process, diverse 3D structures can be realized in the same 2D pattern.The mechanism of multi-morphic construction can be explained by FEM modeling.Furthermore, as a proof-of-concept demonstration, six types of 3D nanomembranes are prepared, realizing incident light angle detection with the assistance of DNN.Our established quasistatic multilevel FEM lays a theoretical foundation for 2D nanomembrane self-assembly behavior, providing a construction paradigm for future optoelectronics, micro-electro-mechanical systems, and microrobots.

Fig. 1 |
Fig. 1 | Releasing of pre-strain rectangle bilayer nanomembranes.a Basic model of pre-strain nanomembrane.b Wet release process of pre-strain bilayer nanomembrane with a fixed edge.c Schematic of sacrificial layer wet etching model: etchant diffusion from liquid phase to solid-liquid interface and chemical reaction between etchant and sacrificial layer in interface.d Self-assembled 3D structuretypes change respect to aspect ratio.Schematic of velocity ratio and aspect ratio (inset).e Elastic energy difference of nanomembranes that is released from adjacent and opposite edge for different aspect ratios.Source data are provided as a Source Data file.

Fig. 2 |
Fig. 2 | Quasistatic multilevel model FEM simulation of bilayer pre-strain nanomembrane release.The modeling process of multilevel quasistatic FEM: a simulating sacrificial layer boundary movement of the release process, b collecting coordinates of boundary in discrete time points, c importing coordinates into the dynamic simulation as a boundary condition.d Modeling and the result of quasistatic multilevel FEM modeling, in which moving boundaries are derived from coordinates in the reaction-diffusion model.Quasistatic multilevel FEM simulation and experimental results of e 60 nm Si/40 nm Cr triangle pattern, f 25 nm LF SiN x /25 nm HF SiN x rectangle pattern, g 15 nm NiTi (0.2 Å s −1 )/ 15 nm NiTi (1.5 Å s −1 ) rectangle pattern, and h 120 nm VO 2 /30 nm Cr rectangle patterns.Scale bars of e-g are 10 μm, and scale bar of h is 100 μm.Pink areas in SEM images are self-assembled structures.Source data are provided as a Source Data file.
Fig. 3 | 3D structures released from bilayer nanomembrane.a SEM image of released Si/Cr bilayer nanomembrane with fixed length L of 40 μm and various widths W from 2 to 80 μm with step of 2 μm.Scale bars, 100 μm.b Five structures: i. ring, ii.arch, iii.helix, iv.taper, and v. tube.Scale bars, 10 μm.c Statistical table of the probability of morphological features in pattern with different sizes of L from 20 to 80 μm with the step of 5 μm, and W from 2 to 80 μm with the step of 2 μm, where the order numbers represent corresponding structures in b. d Experimental

Fig. 4 |c
Fig. 4 | Assembly mechanism and elastic energy of self-assembly nanomembrane.a Diagram of idealized bilateral release models that consist of two regions released chronologically.b Relative elastic energy of idealized bilateral release models with different direction angles and pattern widths via FEM simulation.c Relationship between direction angle, distribution of various morphology, and relative strain energy.d Multilevel FEM simulations and e experimental results of relationship between direction angle α and etching time.Source data are provided as a Source Data file.

Fig. 5 |Fig. 6 |
Fig. 5 | Configurable structure assembly by controlling etching boundary.a Etchant rising and b etchant etching direction in adjacent side preferential etching.c Etchant rising and d etching direction in opposite side preferential etching.e Experimental and quasistatic multilevel FEM simulation results of Si/Cr bilayer with different sizes in unreleased condition, bottom etching condition, and top etching condition.f Mechanism, experimental results, and quasistatic FEM simulation results of multimorphic assembly in Si/Cr bilayer parallelogram patterns.Scale bars, 20 μm.Pink areas in SEM images are self-assembled structures.Source data are provided as a Source Data file.

Table 1 |
Yield of large-scale prepared different 3D structures